Inferring explicit weighted consensus networks to represent alternative evolutionary histories

Background: The advent of molecular biology techniques and constant increase in availability of genetic material have triggered the development of many phylogenetic tree inference methods. However, several reticulate evolution processes, such as horizontal gene transfer and hybridization, have been shown to blur the species\ud evolutionary history by causing discordance among phylogenies inferred from different genes.\ud Methods: To tackle this problem, we hereby describe a new method for inferring and representing alternative(reticulate) evolutionary histories of species as an explicit weighted consensus network which can be constructed from a collection of gene trees with or without prior knowledge of the species phylogeny.\ud Results: We provide a way of building a weighted phylogenetic network for each of the following reticulation\ud mechanisms: diploid hybridization, intragenic recombination and complete or partial horizontal gene transfer. We successfully tested our method on some synthetic and real datasets to infer the above-mentioned evolutionary events which may have influenced the evolution of many species.\ud Conclusions: Our weighted consensus network inference method allows one to infer, visualize and validate statistically major conflicting signals induced by the mechanisms of reticulate evolution. The results provided by the new method can be used to represent the inferred conflicting signals by means of explicit and easy-to-interpret phylogenetic networks

Comparative Analysis of RNA Families Reveals Distinct Repertoires for Each Domain of Life

The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes and genetically-encoded protein catalysts, has been central to models for the early evolution of life on Earth. A key part of such models is continuity between the earliest stages in the evolution of life and the RNA repertoires of extant lineages. Some assessments seem consistent with a diverse RNA world, yet direct continuity between modern RNAs and an RNA world has not been demonstrated for the majority of RNA families, and, anecdotally, many RNA functions appear restricted in their distribution. Despite much discussion of the possible antiquity of RNA families, no systematic analyses of RNA family distribution have been performed. To chart the broad evolutionary history of known RNA families, we performed comparative genomic analysis of over 3 million RNA annotations spanning 1446 families from the Rfam 10 database. We report that 99% of known RNA families are restricted to a single domain of life, revealing discrete repertoires for each domain. For the 1% of RNA families/clans present in more than one domain, over half show evidence of horizontal gene transfer, and the rest show a vertical trace, indicating the presence of a complex protein synthesis machinery in the Last Universal Common Ancestor (LUCA) and consistent with the evolutionary history of the most ancient protein-coding genes. However, with limited interdomain transfer and few RNA families exhibiting demonstrable antiquity as predicted under RNA world continuity, our results indicate that the majority of modern cellular RNA repertoires have primarily evolved in a domain-specific manner.Comment: 47 pages, 4 main figures, 3 supplementary figures, 4 supplementary tables. Submitted to PLOS Computational Biolog

ZOMBIES IN BACTERIAL GENOMES: IDENTIFICATION AND ANALYSIS OF PREVIOUSLY VIRULENT PHAGE

Bacteriophage (or ‘phage’) are viruses that infect and reproduce within their bacterial hosts. They have a major global impact on bacterial evolution and ecology, and might influence the pathogenicity of their host bacterium by providing virulence factors. Phage can either be described as “virulent” or “temperate”; the distinguishing feature between the two is their method of replication. This study sought to identify phage sequences within bacterial host genomes and determine the life cycle of the phage, exploring whether there is a connection between defective phage and previously virulent phage. It would normally be expected that any phage sequences identified within a bacterial host would have a temperate life cycle, since only temperate phage enter the lysogenic cycle and insert their DNA into the host as a ‘prophage,’ while virulent phage replicate via the lytic cycle, in which phage DNA replicates separately from that of the host’s and infected cells are lysed. Defective phage–‘zombies’ in bacterial genomes–are dormant phage that have become inactive through mutational decay or some other process. It is possible that some of these defective phage are in fact previously virulent phage that have become accidentally inserted within the host genome. This study detected phage within bacterial genomes using the prophage identification tools PHAge Search Tool (PHAST) and Prophage Finder. Identified sequences were categorized as ‘intact,’ ‘questionable,’ or ‘incomplete’; questionable and incomplete phage were classified as defective. The lifestyles of the uncovered phage sequences were then determined using PHACTS; six phage were identified as possibly virulent. The life cycles of the phage were further analyzed by assessing the genomic signature distances (GSD) and codon adaptation indexes (CAI) for each phage. Three phage were shown to have a GSD consistent with a virulent life cycle, and the CAI values of four phage corresponded with that of virulent phage. Although previous studies have indicated that some virulent phage may have a temperate lineage, identifying prophage as previously virulent is a novel finding. This has implications for our understanding of phage life cycles and the infection process, as it challenges the idea that only temperate phage insert their DNA into the host genome

A reservoir of 'historical' antibiotic resistance genes in remote pristine Antarctic soils

Background: Soil bacteria naturally produce antibiotics as a competitive mechanism, with a concomitant evolution, and exchange by horizontal gene transfer, of a range of antibiotic resistance mechanisms. Surveys of bacterial resistance elements in edaphic systems have originated primarily from human-impacted environments, with relatively little information from remote and pristine environments, where the resistome may comprise the ancestral gene diversity. Methods: We used shotgun metagenomics to assess antibiotic resistance gene (ARG) distribution in 17 pristine and remote Antarctic surface soils within the undisturbed Mackay Glacier region. We also interrogated the phylogenetic placement of ARGs compared to environmental ARG sequences and tested for the presence of horizontal gene transfer elements flanking ARGs. Results: In total, 177 naturally occurring ARGs were identified, most of which encoded single or multi-drug efflux pumps. Resistance mechanisms for the inactivation of aminoglycosides, chloramphenicol and beta-lactam antibiotics were also common. Gram-negative bacteria harboured most ARGs (71%), with fewer genes from Gram-positive Actinobacteria and Bacilli (Firmicutes) (9%), reflecting the taxonomic composition of the soils. Strikingly, the abundance of ARGs per sample had a strong, negative correlation with species richness (r=-0.49, P < 0.05). This result, coupled with a lack of mobile genetic elements flanking ARGs, suggests that these genes are ancient acquisitions of horizontal transfer events. Conclusions: ARGs in these remote and uncontaminated soils most likely represent functional efficient historical genes that have since been vertically inherited over generations. The historical ARGs in these pristine environments carry a strong phylogenetic signal and form a monophyletic group relative to ARGs from other similar environments

A Network Approach to Analyzing Highly Recombinant Malaria Parasite Genes

The var genes of the human malaria parasite Plasmodium falciparum present a challenge to population geneticists due to their extreme diversity, which is generated by high rates of recombination. These genes encode a primary antigen protein called PfEMP1, which is expressed on the surface of infected red blood cells and elicits protective immune responses. Var gene sequences are characterized by pronounced mosaicism, precluding the use of traditional phylogenetic tools that require bifurcating tree-like evolutionary relationships. We present a new method that identifies highly variable regions (HVRs), and then maps each HVR to a complex network in which each sequence is a node and two nodes are linked if they share an exact match of significant length. Here, networks of var genes that recombine freely are expected to have a uniformly random structure, but constraints on recombination will produce network communities that we identify using a stochastic block model. We validate this method on synthetic data, showing that it correctly recovers populations of constrained recombination, before applying it to the Duffy Binding Like-α (DBLα) domain of var genes. We find nine HVRs whose network communities map in distinctive ways to known DBLα classifications and clinical phenotypes. We show that the recombinational constraints of some HVRs are correlated, while others are independent. These findings suggest that this micromodular structuring facilitates independent evolutionary trajectories of neighboring mosaic regions, allowing the parasite to retain protein function while generating enormous sequence diversity. Our approach therefore offers a rigorous method for analyzing evolutionary constraints in var genes, and is also flexible enough to be easily applied more generally to any highly recombinant sequences

Comparative Genomics of Ape Plasmodium Parasites Reveals Key Evolutionary Events Leading to Human Malaria

African great apes are infected with at least six species of P. falciparum-like parasites, including the ancestor of P. falciparum. Comparative studies of these parasites and P. falciparum (collectively termed the Laverania subgenus) will provide insight into the evolutionary origins of P. falciparum and identify genetic features that influence host tropism. Here we show that ape Laverania parasites do not serve as a recurrent source of human malaria and use novel enrichment techniques to derive near full-length genomes of close and distant relatives of P. falciparum. Using a combination of single template amplification and deep sequencing, we observe no evidence of ape Laverania infections in forest dwelling humans in Cameroon. This result supports previous findings that ape Laverania parasites are host specific and have successfully colonized humans only once. To understand the determinants of host specificity and identify genetic characteristics unique to P. falciparum, we develop a novel method for selective enrichment of Plasmodium DNA from sub-microscopically infected whole blood samples. We use this technique to enrich for Laverania genomic DNA from chimpanzee blood samples and assemble near full length genomes for both close (P. reichenowi) and distant (P. gaboni) relatives of P. falciparum. Comparative analyses of these genomes to P. falciparum identify features that are conserved across the Laverania subgenus, including the expansion of the FIKK kinases and the presence of var-like multigene families in all Laverania species. Our analyses also identify genetic features that are unique to P. falciparum, such as a very low within-species diversity and a complex evolutionary history of the essential invasion genes RH5 and CyRPA. This dissertation lays the groundwork for future comparative analyses of the Laverania subgenus including population genomic analyses of ape parasites and comparisons of P. falciparum to its ancestor, P. praefalciparum

At the nexus of three kingdoms: the genome of the mycorrhizal fungus Gigaspora margarita provides insights into plant, endobacterial and fungal interactions.

As members of the plant microbiota, arbuscular mycorrhizal fungi (AMF, Glomeromycotina) symbiotically colonize plant roots. AMF also possess their own microbiota, hosting some uncultivable endobacteria. Ongoing research has revealed the genetics underlying plant responses to colonization by AMF, but the fungal side of the relationship remains in the dark. Here, we sequenced the genome of Gigaspora margarita, a member of the Gigasporaceae in an early diverging group of the Glomeromycotina. In contrast to other AMF, G. margarita may host distinct endobacterial populations and possesses the largest fungal genome so far annotated (773.104 Mbp), with more than 64% transposable elements. Other unique traits of the G. margarita genome include the expansion of genes for inorganic phosphate metabolism, the presence of genes for production of secondary metabolites and a considerable number of potential horizontal gene transfer events. The sequencing of G. margarita genome reveals the importance of its immune system, shedding light on the evolutionary pathways that allowed early diverging fungi to interact with both plants and bacteria

Evolutionary patchwork of an insecticidal toxin shared between plant-associated pseudomonads and the insect pathogens Photorhabdus and Xenorhabdus

Background: Root-colonizing fluorescent pseudomonads are known for their excellent abilities to protect plants against soil-borne fungal pathogens. Some of these bacteria produce an insecticidal toxin (Fit) suggesting that they may exploit insect hosts as a secondary niche. However, the ecological relevance of insect toxicity and the mechanisms driving the evolution of toxin production remain puzzling. Results: Screening a large collection of plant-associated pseudomonads for insecticidal activity and presence of the Fit toxin revealed that Fit is highly indicative of insecticidal activity and predicts that Pseudomonas protegens and P. chlororaphis are exclusive Fit producers. A comparative evolutionary analysis of Fit toxin-producing Pseudomonas including the insect-pathogenic bacteria Photorhabdus and Xenorhadus, which produce the Fit related Mcf toxin, showed that fit genes are part of a dynamic genomic region with substantial presence/absence polymorphism and local variation in GC base composition. The patchy distribution and phylogenetic incongruence of fit genes indicate that the Fit cluster evolved via horizontal transfer, followed by functional integration of vertically transmitted genes, generating a unique Pseudomonas-specific insect toxin cluster. Conclusions: Our findings suggest that multiple independent evolutionary events led to formation of at least three versions of the Mcf/Fit toxin highlighting the dynamic nature of insect toxin evolution

Insights into bacterial genome composition through variable target GC content profiling

This study presents a new computational method for guanine (G) and cytosine (C), or GC, content profiling based on the idea of multiple resolution sampling (MRS). The benefit of our new approach over existing techniques follows from its ability to locate significant regions without prior knowledge of the sequence, nor the features being sought. The use of MRS has provided novel insights into bacterial genome composition. Key findings include those that are related to the core composition of bacterial genomes, to the identification of large genomic islands (in Enterobacterial genomes), and to the identification of surface protein determinants in human pathogenic organisms (e.g., Staphylococcus genomes). We observed that bacterial surface binding proteins maintain abnormal GC content, potentially pointing to a viral origin. This study has demonstrated that GC content holds a high informational worth and hints at many underlying evolutionary processes.&nbsp;For online Supplementary Material, see www.liebertonline.com