Retrofitting the genome: L1 extinction follows endogenous retroviral expansion in a group of muroid rodents

Long interspersed nuclear element 1 (LINE-1; L1) retrotransposons are the most common retroelements in mammalian genomes. Unlike individual families of endogenous retroviruses (ERVs), they have remained active throughout the mammalian radiation and are responsible for most of the retroelement movement and much genome rearrangement within mammals. They can be viewed as occupying a substantial niche within mammalian genomes. Our previous demonstration that L1s and B1 short interspersed nuclear elements (SINEs) are inactive in a group of South American rodents led us to ask if other elements have amplified to fill the empty niche. We identified a novel and highly active family of ERVs (mysTR). To determine whether loss of L1 activity was correlated with expansion of mysTR, we examined mysTR activity in four South American rodent species that have lost L1 and B1 activity and four sister species with active L1s. The copy number of recent mysTR insertions was extremely high, with an average of 4,200 copies per genome. High copy numbers exist in both L1-active and L1-extinct species, so the mysTR expansion appears to have preceded the loss of both SINE and L1 activity rather than to have filled an empty niche created by their loss. It may be coincidental that two unusual genomic events--loss of L1 activity and massive expansion of an ERV family--occur in the same group of mammals. Alternatively, it is possible that this large ERV expansion set the stage for L1 extinction.     

Predicting prokaryotic ecological niches using genome sequence analysis

Automated DNA sequencing technology is so rapid that analysis has become the rate-limiting step. Hundreds of prokaryotic genome sequences are publicly available, with new genomes uploaded at the rate of approximately 20 per month. As a result, this growing body of genome sequences will include microorganisms not previously identified, isolated, or observed. We hypothesize that evolutionary pressure exerted by an ecological niche selects for a similar genetic repertoire in those prokaryotes that occupy the same niche, and that this is due to both vertical and horizontal transmission. To test this, we have developed a novel method to classify prokaryotes, by calculating their Pfam protein domain distributions and clustering them with all other sequenced prokaryotic species. Clusters of organisms are visualized in two dimensions as 'mountains' on a topological map. When compared to a phylogenetic map constructed using 16S rRNA, this map more accurately clusters prokaryotes according to functional and environmental attributes. We demonstrate the ability of this map, which we term a "niche map", to cluster according to ecological niche both quantitatively and qualitatively, and propose that this method be used to associate uncharacterized prokaryotes with their ecological niche as a means of predicting their functional role directly from their genome sequence.     

On the Origins of a Vibrio Species

Thirty-two genome sequences of various Vibrionaceae members are compared, with emphasis on what makes V. cholerae unique. As few as 1,000 gene families are conserved across all the Vibrionaceae genomes analysed; this fraction roughly doubles for gene families conserved within the species V. cholerae. Of these, approximately 200 gene families that cluster on various locations of the genome are not found in other sequenced Vibrionaceae; these are possibly unique to the V. cholerae species. By comparing gene family content of the analysed genomes, the relatedness to a particular species is identified for two unspeciated genomes. Conversely, two genomes presumably belonging to the same species have suspiciously dissimilar gene family content. We are able to identify a number of genes that are conserved in, and unique to, V. cholerae. Some of these genes may be crucial to the niche adaptation of this species.     

Evolutionary genomics of a temperate bacteriophage in an obligate intracellular bacteria (Wolbachia)

Genome evolution of bacteria is usually influenced by ecology, such that bacteria with a free-living stage have large genomes and high rates of horizontal gene transfer, while obligate intracellular bacteria have small genomes with typically low amounts of gene exchange. However, recent studies indicate that obligate intracellular species that host-switch frequently harbor agents of horizontal transfer such as mobile elements. For example, the temperate double-stranded DNA bacteriophage WO in Wolbachia persistently transfers between bacterial coinfections in the same host. Here we show that despite the phage's rampant mobility between coinfections, the prophage's genome displays features of constraint related to its intracellular niche. First, there is always at least one intact prophage WO and usually several degenerate, independently-acquired WO prophages in each Wolbachia genome. Second, while the prophage genomes are modular in composition with genes of similar function grouping together, the modules are generally not interchangeable with other unrelated phages and thus do not evolve by the Modular Theory. Third, there is an unusual core genome that strictly consists of head and baseplate genes; other gene modules are frequently deleted. Fourth, the prophage recombinases are diverse and there is no conserved integration sequence. Finally, the molecular evolutionary forces acting on prophage WO are point mutation, intragenic recombination, deletion, and purifying selection. Taken together, these analyses indicate that while lateral transfer of phage WO is pervasive between Wolbachia with occasional new gene uptake, constraints of the intracellular niche obstruct extensive mixture between WO and the global phage population. Although the Modular Theory has long been considered the paradigm of temperate bacteriophage evolution in free-living bacteria, it appears irrelevant in phages of obligate intracellular bacteria.     

The temporal dynamics of slightly deleterious mutations in Escherichia coli and Shigella spp

The Shigella are recently emerged clones of Escherichia coli, which have independently adopted an intracellular pathogenic lifestyle. We examined the molecular evolutionary consequences of this niche specialization by comparing the normalized, directional frequency profiles of unique polymorphisms within 2,098 orthologues representing the intersection of five E. coli and four Shigella genomes. We note a surfeit of AT-enriching changes (GC-->AT), transversions, and nonsynonymous changes in the Shigella genomes. By examining these differences within a temporal framework, we conclude that our results are consistent with relaxed or inefficient selection in Shigella owing to a reduced effective population size. Alternative interpretations, and the interesting exception of Shigella sonnei, are discussed. Finally, this analysis lends support to the view that nucleotide composition typically does not lie at mutational equilibrium but that selection plays a role in maintaining a higher GC content than would result solely from mutation bias. This argument sheds light on the enrichment of adenine and thymine in the genomes of bacterial endosymbionts where purifying selection is very weak.     

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined “species”, whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.     

Heterozygosity in Yellowstone Park elk,Cervus canadensis

Protein products of 24 loci from the genomes of Yellowstone Park elk were analyzed by electrophoresis. Heterozygosity was detected in only one system, making elk much less polymorphic than eastern whitetailed deer. Data for several other large mammals are compared with those for elk and reveal similarly low levels of isozymic variation. The data are consistent with the fine-grained niche theory but difficult to reconcile with bottlenecks and genetic drift.     

Implications of streamlining theory for microbial ecology

Whether a small cell, a small genome or a minimal set of chemical reactions with self-replicating properties, simplicity is beguiling. As Leonardo da Vinci reportedly said, ‘simplicity is the ultimate sophistication''. Two diverging views of simplicity have emerged in accounts of symbiotic and commensal bacteria and cosmopolitan free-living bacteria with small genomes. The small genomes of obligate insect endosymbionts have been attributed to genetic drift caused by small effective population sizes (N_e). In contrast, streamlining theory attributes small cells and genomes to selection for efficient use of nutrients in populations where N_e is large and nutrients limit growth. Regardless of the cause of genome reduction, lost coding potential eventually dictates loss of function. Consequences of reductive evolution in streamlined organisms include atypical patterns of prototrophy and the absence of common regulatory systems, which have been linked to difficulty in culturing these cells. Recent evidence from metagenomics suggests that streamlining is commonplace, may broadly explain the phenomenon of the uncultured microbial majority, and might also explain the highly interdependent (connected) behavior of many microbial ecosystems. Streamlining theory is belied by the observation that many successful bacteria are large cells with complex genomes. To fully appreciate streamlining, we must look to the life histories and adaptive strategies of cells, which impose minimum requirements for complexity that vary with niche.     

A Novel Computational Method Identifies Intra- and Inter-Species Recombination Events in Staphylococcus aureus and Streptococcus pneumoniae

Advances in high-throughput DNA sequencing technologies have determined an explosion in the number of sequenced bacterial genomes. Comparative sequence analysis frequently reveals evidences of homologous recombination occurring with different mechanisms and rates in different species, but the large-scale use of computational methods to identify recombination events is hampered by their high computational costs. Here, we propose a new method to identify recombination events in large datasets of whole genome sequences. Using a filtering procedure of the gene conservation profiles of a test genome against a panel of strains, this algorithm identifies sets of contiguous genes acquired by homologous recombination. The locations of the recombination breakpoints are determined using a statistical test that is able to account for the differences in the natural rate of evolution between different genes. The algorithm was tested on a dataset of 75 genomes of Staphylococcus aureus and 50 genomes comprising different streptococcal species, and was able to detect intra-species recombination events in S. aureus and in Streptococcus pneumoniae. Furthermore, we found evidences of an inter-species exchange of genetic material between S. pneumoniae and Streptococcus mitis, a closely related commensal species that colonizes the same ecological niche. The method has been implemented in an R package, Reco, which is freely available from supplementary material, and provides a rapid screening tool to investigate recombination on a genome-wide scale from sequence data.     

Gene loss through pseudogenization contributes to the ecological diversification of a generalist Roseobacter lineage

Ecologically relevant genes generally show patchy distributions among related bacterial genomes. This is commonly attributed to lateral gene transfer, whereas the opposite mechanism—gene loss—has rarely been explored. Pseudogenization is a major mechanism underlying gene loss, and pseudogenes are best characterized by comparing closely related genomes because of their short life spans. To explore the role of pseudogenization in microbial ecological diversification, we apply rigorous methods to characterize pseudogenes in the 279 newly sequenced Ruegeria isolates of the globally abundant Roseobacter group collected from two typical coastal habitats in Hong Kong, the coral Platygyra acuta and the macroalga Sargassum hemiphyllum. Pseudogenes contribute to ~16% of the accessory genomes of these strains. Ancestral state reconstruction reveals that many pseudogenization events are correlated with ancestral niche shifts. Specifically, genes related to resource scavenging and energy acquisition were often pseudogenized when roseobacters inhabiting carbon-limited and energy-poor coral skeleton switched to other resource-richer niches. For roseobacters inhabiting the macroalgal niches, genes for nitrogen regulation and carbohydrate utilization were important but became dispensable upon shift to coral skeleton where nitrate is abundant but carbohydrates are less available. Whereas low-energy-demanding secondary transporters are more favorable in coral skeleton, ATP-driven primary transporters are preferentially kept in the energy-replete macroalgal niches. Moreover, a large proportion of these families mediate organismal interactions, suggesting their rapid losses by pseudogenization as a potential response to host and niche shift. These findings illustrate an important role of pseudogenization in shaping genome content and driving ecological diversification of marine roseobacters.Subject terms: Ecology, Evolution, Microbiology     

Species-specific protein sequence and fold optimizations

Background

An organism's ability to adapt to its particular environmental niche is of fundamental importance to its survival and proliferation. In the largest study of its kind, we sought to identify and exploit the amino-acid signatures that make species-specific protein adaptation possible across 100 complete genomes.

Results

Environmental niche was determined to be a significant factor in variability from correspondence analysis using the amino acid composition of over 360,000 predicted open reading frames (ORFs) from 17 archae, 76 bacteria and 7 eukaryote complete genomes. Additionally, we found clusters of phylogenetically unrelated archae and bacteria that share similar environments by amino acid composition clustering. Composition analyses of conservative, domain-based homology modeling suggested an enrichment of small hydrophobic residues Ala, Gly, Val and charged residues Asp, Glu, His and Arg across all genomes. However, larger aromatic residues Phe, Trp and Tyr are reduced in folds, and these results were not affected by low complexity biases. We derived two simple log-odds scoring functions from ORFs (C_G) and folds (C_F) for each of the complete genomes. C_F achieved an average cross-validation success rate of 85 ± 8% whereas the C_G detected 73 ± 9% species-specific sequences when competing against all other non-redundant C_G. Continuously updated results are available at http://genome.mshri.on.ca.

Conclusion

Our analysis of amino acid compositions from the complete genomes provides stronger evidence for species-specific and environmental residue preferences in genomic sequences as well as in folds. Scoring functions derived from this work will be useful in future protein engineering experiments and possibly in identifying horizontal transfer events.     

Evolutionary analysis of glycosyl hydrolase family 28 (GH28) suggests lineage-specific expansions in necrotrophic fungal pathogens

Glycosyl hydrolase family 28 (GH28) is a set of structurally related enzymes that hydrolyze glycosidic bonds in pectin, and are important extracellular enzymes for both pathogenic and saprotrophic fungi. Yet, very little is understood about the evolutionary forces driving the diversification of GH28s in fungal genomes. We reconstructed the evolutionary history of family GH28 in fungi by examining the distribution of GH28 copy number across the phylogeny of fungi, and by reconstructing the phylogeny of GH28 genes. We also examined the relationship between lineage-specific GH28 expansions and fungal ecological strategy, testing the hypothesis that GH28 evolution in fungi is driven by ecological strategy (pathogenic vs. non-pathogenic) and pathogenic niche (necrotrophic vs. biotrophic). Our results showed that GH28 phylogeny of Ascomycota and Basidiomycota sequences was structured by specific biochemical function, with endo-polygalacturonases and endo-rhamnogalacturonases forming distinct, apparently ancient clades, while exo-polygalacturonases are more widely distributed. In contrast, Mucoromycotina and Stramenopile sequences formed taxonomically-distinct clades. Large, lineage-specific variation in GH28 copy number indicates that the evolution of this gene family is consistent with the birth-and-death model of gene family evolution, where diversity of GH28 loci within genomes was generated through multiple rounds of gene duplication followed by functional diversification and loss of some gene family members. Although GH28 copy number was correlated with genome size, our findings suggest that ecological strategy also plays an important role in determining the GH28 repertoire of fungi. Both necrotrophic and biotrophic fungi have larger genomes than non-pathogens, yet only necrotrophs possess more GH28 enzymes than non-pathogens. Hence, lineage-specific GH28 expansion is the result of both variation in genome size across fungal species and diversifying selection within the necrotrophic plant pathogen ecological niche. GH28 evolution among necrotrophs has likely been driven by a co-evolutionary arms race with plants, whereas the need to avoid plant immune responses has resulted in purifying selection within biotrophic fungi.     

Transposable elements: powerful facilitators of evolution

Transposable elements (TEs) are powerful facilitators of genome evolution, and hence of phenotypic diversity as they can cause genetic changes of great magnitude and variety. TEs are ubiquitous and extremely ancient, and although harmful to some individuals, they can be very beneficial to lineages. TEs can build, sculpt, and reformat genomes by both active and passive means. Lineages with active TEs or with abundant homogeneous inactive populations of TEs that can act passively by causing ectopic recombination are potentially fecund, adaptable, and taxonate readily. Conversely, taxa deficient in TEs or possessing heterogeneous populations of inactive TEs may be well adapted in their niche, but tend to prolonged stasis and may risk extinction by lacking the capacity to adapt to change, or diversify. Because of recurring intermittent waves of TE infestation, available data indicate a compatibility with punctuated equilibrium, in keeping with widely accepted interpretations of evidence from the fossil record. We propose a general and holistic synthesis on how the presence of TEs within genomes makes them flexible and dynamic, so that genomes themselves are powerful facilitators of their own evolution     

Congruent evolution of different classes of non-coding DNA in prokaryotic genomes

Prokaryotic genomes are considered to be 'wall-to-wall' genomes, which consist largely of genes for proteins and structural RNAs, with only a small fraction of the genomic DNA allotted to intergenic regions, which are thought to typically contain regulatory signals. The majority of bacterial and archaeal genomes contain 6-14% non-coding DNA. Significant positive correlations were detected between the fraction of non-coding DNA and inter- and intra-operonic distances, suggesting that different classes of non-coding DNA evolve congruently. In contrast, no correlation was found between any of these characteristics of non-coding sequences and the number of genes or genome size. Thus, the non-coding regions and the gene sets in prokaryotes seem to evolve in different regimes. The evolution of non-coding regions appears to be determined primarily by the selective pressure to minimize the amount of non-functional DNA, while maintaining essential regulatory signals, because of which the content of non-coding DNA in different genomes is relatively uniform and intra- and inter-operonic non-coding regions evolve congruently. In contrast, the gene set is optimized for the particular environmental niche of the given microbe, which results in the lack of correlation between the gene number and the characteristics of non-coding regions.     

Expanding our view of genomic diversity in Candidatus Accumulibacter clades

Enhanced biological phosphorus removal (EBPR) is an important industrial wastewater treatment process mediated by polyphosphate‐accumulating organisms (PAOs). Members of the genus Candidatus Accumulibacter are one of the most extensively studied PAO as they are commonly enriched in lab‐scale EBPR reactors. Members of different Accumulibacter clades are often enriched through changes in reactor process conditions; however, the two currently sequenced Accumulibacter genomes show extensive metabolic similarity. Here, we expand our understanding of Accumulibacter genomic diversity through recovery of eight population genomes using deep metagenomics, including seven from phylogenetic clades with no previously sequenced representative. Comparative genomic analysis revealed a core of shared genes involved primarily in carbon and phosphorus metabolism; however, each Accumulibacter genome also encoded a substantial number of unique genes (> 700 genes). A major difference between the Accumulibacter clades was the type of nitrate reductase encoded and the capacity to perform subsequent steps in denitrification. The Accumulibacter clade IIF genomes also contained acetaldehyde dehydrogenase that may allow ethanol to be used as carbon source. These differences in metabolism between Accumulibacter genomes provide a molecular basis for niche differentiation observed in lab‐scale reactors and may offer new opportunities for process optimization.     

Intragenomic enzyme complements

Researchers in applied biocatalysis are now reaping the rewards of intensive effort and technological developments in the sequencing of the genomes of microbial and plant species. The genomic resource contains the sequences of millions of new genes with potential application in industrial biotechnology and includes families of enzymes within discrete genomes that potentially catalyze equivalent chemical reactions. One of the key emerging characteristics of these intragenomic complements of enzymes is the impressive breadth of catalytic diversity that is observed within them. This diversity may have been acquired either in order to combat the spectrum of metabolic challenges with which the organism may be presented in its natural environment or as part of the biosynthetic machinery evolved to produce a spectrum of secondary metabolites that will prove to be advantageous in establishing a niche. Attempts have been made to functionally characterize the intragenomic complements of enzyme families catalyzing diverse reactions including carbonyl reduction, ester hydrolysis and Baeyer–Villiger oxidation, in Gram-positive bacteria, yeasts, filamentous fungi and the plant Arabidopsis thaliana. These studies are beginning to describe in detail for the first time the impressive range of catalytic potential within single organisms for attributes such as substrate range, enantioselectivity or thermostability, each of which is of interest from an enzyme discovery perspective.     

Widespread distribution of prophages signaling the potential for adaptability and pathogenicity evolution of Ralstonia solanacearum species complex

Integrated bacteriophages (prophages) can impact host cells, affecting their lifestyle, genomic diversity, and fitness. However, many basic aspects of how these organisms affect the host cell remain poorly understood. Ralstonia solanacearum is a gram-negative plant pathogenic bacterium that encompasses a great diversity of ecotypes regarded as a species complex (R. solanacearum Species Complex - RSSC). RSSC genomes have a mosaic structure containing numerous elements, signaling the potential for its evolution through horizontal gene transfer. Here, we analyzed 120 Ralstonia spp. genomes from the public database to identify prophage sequences. In total, 379 prophage-like elements were found in the chromosome and megaplasmid of Ralstonia spp. These elements encode genes related to host fitness, virulence factors, antibiotic resistance, and niche adaptation, which might contribute to RSSC adaptability. Prophage-like elements are widespread into the complex in different species and geographic origins, suggesting that the RSSC phages are ancestrally acquired. Complete prophages belonging to the families Inoviridae, Myoviridae, and Siphoviridae were found, being the members of Inoviridae the most abundant. Analysis of CRISPR-Cas spacer sequences demonstrated the presence of prophages sequences that indicate successive infection events during bacterial evolution. Besides complete prophages, we also demonstrated 14 novel putative prophages integrated into Ralstonia spp. genomes. Altogether, our results provide insights into the diversity of prophages in RSSC genomes and suggest that these elements may deeply affect the virulence and host adaptation and shaping the genomes among the strains of this important pathogen.     

Predicting Ecological Roles in the Rhizosphere Using Metabolome and Transportome Modeling

The ability to obtain complete genome sequences from bacteria in environmental samples, such as soil samples from the rhizosphere, has highlighted the microbial diversity and complexity of environmental communities. However, new algorithms to analyze genome sequence information in the context of community structure are needed to enhance our understanding of the specific ecological roles of these organisms in soil environments. We present a machine learning approach using sequenced Pseudomonad genomes coupled with outputs of metabolic and transportomic computational models for identifying the most predictive molecular mechanisms indicative of a Pseudomonad’s ecological role in the rhizosphere: a biofilm, biocontrol agent, promoter of plant growth, or plant pathogen. Computational predictions of ecological niche were highly accurate overall with models trained on transportomic model output being the most accurate (Leave One Out Validation F-scores between 0.82 and 0.89). The strongest predictive molecular mechanism features for rhizosphere ecological niche overlap with many previously reported analyses of Pseudomonad interactions in the rhizosphere, suggesting that this approach successfully informs a system-scale level understanding of how Pseudomonads sense and interact with their environments. The observation that an organism’s transportome is highly predictive of its ecological niche is a novel discovery and may have implications in our understanding microbial ecology. The framework developed here can be generalized to the analysis of any bacteria across a wide range of environments and ecological niches making this approach a powerful tool for providing insights into functional predictions from bacterial genomic data.