Thermal Adaptation of the Small Subunit Ribosomal RNA Gene: A Comparative Study

Differential Selection for Translation Efficiency Shapes Translation Machineries in Bacterial Species

Different bacterial species have dramatically different generation times, from 20-30 min in Escherichia coli to about two weeks in Mycobacterium leprae. The translation machinery in a cell needs to synthesize all proteins for a new cell in each generation. The three subprocesses of translation, i.e., initiation, elongation, and termination, are expected to be under stronger selection pressure to optimize in short-generation bacteria (SGB) such as Vibrio natriegens than in the long-generation Mycobacterium leprae. The initiation efficiency depends on the start codon decoded by the initiation tRNA, the optimal Shine-Dalgarno (SD) decoded by the anti-SD (aSD) sequence on small subunit rRNA, and the secondary structure that may embed the initiation signals and prevent them from being decoded. The elongation efficiency depends on the tRNA pool and codon usage. The termination efficiency in bacteria depends mainly on the nature of the stop codon and the nucleotide immediately downstream of the stop codon. By contrasting SGB with long-generation bacteria (LGB), we predict (1) SGB to have more ribosome RNA operons to produce ribosomes, and more tRNA genes for carrying amino acids to ribosomes, (2) SGB to have a higher percentage of genes using AUG as the start codon and UAA as the stop codon than LGB, (3) SGB to exhibit better codon and anticodon adaptation than LGB, and (4) SGB to have a weaker secondary structure near the translation initiation signals than LGB. These differences between SGB and LGB should be more pronounced in highly expressed genes than the rest of the genes. We present empirical evidence in support of these predictions.

A Conserved Ribosomal Protein Has Entirely Dissimilar Structures in Different Organisms

Ribosomes from different species can markedly differ in their composition by including dozens of ribosomal proteins that are unique to specific lineages but absent in others. However, it remains unknown how ribosomes acquire new proteins throughout evolution. Here, to help answer this question, we describe the evolution of the ribosomal protein msL1/msL2 that was recently found in ribosomes from the parasitic microorganism clade, microsporidia. We show that this protein has a conserved location in the ribosome but entirely dissimilar structures in different organisms: in each of the analyzed species, msL1/msL2 exhibits an altered secondary structure, an inverted orientation of the N-termini and C-termini on the ribosomal binding surface, and a completely transformed 3D fold. We then show that this fold switching is likely caused by changes in the ribosomal msL1/msL2-binding site, specifically, by variations in rRNA. These observations allow us to infer an evolutionary scenario in which a small, positively charged, de novo-born unfolded protein was first captured by rRNA to become part of the ribosome and subsequently underwent complete fold switching to optimize its binding to its evolving ribosomal binding site. Overall, our work provides a striking example of how a protein can switch its fold in the context of a complex biological assembly, while retaining its specificity for its molecular partner. This finding will help us better understand the origin and evolution of new protein components of complex molecular assemblies-thereby enhancing our ability to engineer biological molecules, identify protein homologs, and peer into the history of life on Earth.

Genomic profiling and characteristics of a C1 degrading heterotrophic fresh-water bacterium Paracoccus sp. strain DMF

Paracoccus species are metabolically versatile gram-negative, aerobic facultative methylotrophic bacteria showing enormous promise for environmental and bioremediation studies. Here we report, the complete genome analysis of Paracoccus sp. strain DMF (P. DMF) that was isolated from a domestic wastewater treatment plant in Kanpur, India (26.4287 °N, 80.3891 °E) based on its ability to degrade a recalcitrant organic solvent N, N-dimethylformamide (DMF). The results reveal a genome size of 4,202,269 base pairs (bp) with a G + C content of 67.9%. The assembled genome comprises 4141 coding sequences (CDS), 46 RNA sequences, and 2 CRISPRs. Interestingly, catabolic operons related to the conventional marine-based methylated amines (MAs) degradation pathway were functionally annotated within the genome of an obligated aerobic heterotroph that is P. DMF. The genomic data-based characterization presented here for the novel heterotroph P. DMF aims to improve the understanding of the phenotypic gene products, enzymes, and pathways involved with greater emphasis on facultative methylotrophic motility-based latent pathogenicity.

Ribosomal proteins can hold a more accurate record of bacterial thermal adaptation compared to rRNA

Ribosomal genes are widely used as 'molecular clocks' to infer evolutionary relationships between species. However, their utility as 'molecular thermometers' for estimating optimal growth temperature of microorganisms remains uncertain. Previously, some estimations were made using the nucleotide composition of ribosomal RNA (rRNA), but the universal application of this approach was hindered by numerous outliers. In this study, we aimed to address this problem by identifying additional indicators of thermal adaptation within the sequences of ribosomal proteins. By comparing sequences from 2021 bacteria with known optimal growth temperature, we identified novel indicators among the metal-binding residues of ribosomal proteins. We found that these residues serve as conserved adaptive features for bacteria thriving above 40°C, but not at lower temperatures. Furthermore, the presence of these metal-binding residues exhibited a stronger correlation with the optimal growth temperature of bacteria compared to the commonly used correlation with the 16S rRNA GC content. And an even more accurate correlation was observed between the optimal growth temperature and the YVIWREL amino acid content within ribosomal proteins. Overall, our work suggests that ribosomal proteins contain a more accurate record of bacterial thermal adaptation compared to rRNA. This finding may simplify the analysis of unculturable and extinct species.

Microbial genomic trait evolution is dominated by frequent and rare pulsed evolution

On the macroevolutionary time scale, does trait evolution proceed gradually or by rapid bursts (pulses) separated by prolonged periods of stasis or slow evolution? Although studies have shown that pulsed evolution is prevalent in animals, our knowledge about the tempo and mode of evolution across the tree of life is very limited. This long-standing debate calls for a test in bacteria and archaea, the most ancient and diverse forms of life with unique population genetic properties. Using a likelihood-based framework, we show that pulsed evolution is not only present but also prevalent and predominant in microbial genomic trait evolution. We detected two distinct types of pulsed evolution (small frequent and large rare jumps) that are predicted by the punctuated equilibrium and quantum evolution theories. Our findings suggest that major bacterial lineages could have originated in quick bursts and that pulsed evolution is a common theme across the tree of life.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

Thermal adaptation of mRNA secondary structure: stability versus lability

Macromolecular function commonly involves rapidly reversible alterations in three-dimensional structure (conformation). To allow these essential conformational changes, macromolecules must possess higher order structures that are appropriately balanced between rigidity and flexibility. Because of the low stabilization free energies (marginal stabilities) of macromolecule conformations, temperature changes have strong effects on conformation and, thereby, on function. As is well known for proteins, during evolution, temperature-adaptive changes in sequence foster retention of optimal marginal stability at a species' normal physiological temperatures. Here, we extend this type of analysis to messenger RNAs (mRNAs), a class of macromolecules for which the stability-lability balance has not been elucidated. We employ in silico methods to determine secondary structures and estimate changes in free energy of folding (ΔG_fold) for 25 orthologous mRNAs that encode the enzyme cytosolic malate dehydrogenase in marine mollusks with adaptation temperatures spanning an almost 60 °C range. The change in free energy that occurs during formation of the ensemble of mRNA secondary structures is significantly correlated with adaptation temperature: ΔG_fold values are all negative and their absolute values increase with adaptation temperature. A principal mechanism underlying these adaptations is a significant increase in synonymous guanine + cytosine substitutions with increasing temperature. These findings open up an avenue of exploration in molecular evolution and raise interesting questions about the interaction between temperature-adaptive changes in mRNA sequence and in the proteins they encode.

When bacteria meet mitochondria: The strange case of the tick symbiont Midichloria mitochondrii†

Mitochondria are key eukaryotic organelles that perform several essential functions. Not surprisingly, many intracellular bacteria directly or indirectly target mitochondria, interfering with innate immunity, energy production or apoptosis, to make the host cell a more hospitable niche for bacterial replication. The alphaproteobacterium Midichloria mitochondrii has taken mitochondrial targeting to another level by physically colonising mitochondria, as shown by transmission electron micrographs of bacteria residing in the mitochondrial intermembrane space. This unique localization provokes a number of questions around the mechanisms allowing, and reasons driving intramitochondrial tropism. We suggest possible scenarios that could lead to this peculiar localization and hypothesize potential costs and benefits of mitochondrial colonisation for the bacterium and its host.

Improving Phylogenetic Signals of Mitochondrial Genes Using a New Method of Codon Degeneration

Recovering deep phylogeny is challenging with animal mitochondrial genes because of their rapid evolution. Codon degeneration decreases the phylogenetic noise and bias by aiming to achieve two objectives: (1) alleviate the bias associated with nucleotide composition, which may lead to homoplasy and long-branch attraction, and (2) reduce differences in the phylogenetic results between nucleotide-based and amino acid (AA)-based analyses. The discrepancy between nucleotide-based analysis and AA-based analysis is partially caused by some synonymous codons that differ more from each other at the nucleotide level than from some nonsynonymous codons, e.g., Leu codon TTR in the standard genetic code is more similar to Phe codon TTY than to synonymous CTN codons. Thus, nucleotide similarity conflicts with AA similarity. There are many such examples involving other codon families in various mitochondrial genetic codes. Proper codon degeneration will make synonymous codons more similar to each other at the nucleotide level than they are to nonsynonymous codons. Here, I illustrate a "principled" codon degeneration method that achieves these objectives. The method was applied to resolving the mammalian basal lineage and phylogenetic position of rheas among ratites. The codon degeneration method was implemented in the user-friendly and freely available DAMBE software for all known genetic codes (genetic codes 1 to 33).

Hyperdiverse archaea near life limits at the polyextreme geothermal Dallol area

Microbial life has adapted to various individual extreme conditions; yet, organisms simultaneously adapted to very low pH, high salt and high temperature are unknown. We combined environmental 16S/18S rRNA-gene metabarcoding, cultural approaches, fluorescence-activated cell sorting, scanning electron microscopy and chemical analyses to study samples along such unique polyextreme gradients in the Dallol-Danakil area (Ethiopia). We identify two physicochemical barriers to life in the presence of surface liquid water defined by: i) high chaotropicity-low water activity in Mg²⁺/Ca²⁺-dominated brines and ii) hyperacidity-salt combinations (pH~0/NaCl-dominated salt-saturation). When detected, life was dominated by highly diverse ultrasmall archaea widely distributed across phyla with and without previously known halophilic members. We hypothesize that high cytoplasmic K⁺-level was an original archaeal adaptation to hyperthermophily, subsequently exapted during multiple transitions to extreme halophily. We detect active silica encrustment/fossilization of cells but also abiotic biomorphs of varied chemistry. Our work helps circumscribing habitability and calls for cautionary interpretations of morphological biosignatures on Earth and beyond.

Occurrence of randomly recombined functional 16S rRNA genes in Thermus thermophilus suggests genetic interoperability and promiscuity of bacterial 16S rRNAs

Based on the structural complexity of ribosomes, 16S rRNA genes are considered species-specific and hence used for bacterial phylogenetic analysis. However, a growing number of reports suggest the occurrence of horizontal gene transfer, raising genealogical questions. Here we show the genetic interoperability and promiscuity of 16S rRNA in the ribosomes of an extremely thermophilic bacterium, Thermus thermophilus. The gene in this thermophile was systematically replaced with a diverse array of heterologous genes, resulting in the discovery of various genes that supported growth, some of which were from different phyla. Moreover, numerous functional chimeras were spontaneously generated. Remarkably, cold-adapted mutants were obtained carrying chimeric or full-length heterologous genes, indicating that horizontal gene transfer promoted adaptive evolution. The ribosome may well be understood as a patchworked supramolecule comprising patchworked components. We here propose the “random patch model” for ribosomal evolution.

Integrative rDNAomics-Importance of the Oldest Repetitive Fraction of the Eukaryote Genome

Nuclear ribosomal RNA (rRNA) genes represent the oldest repetitive fraction universal to all eukaryotic genomes. Their deeply anchored universality and omnipresence during eukaryotic evolution reflects in multiple roles and functions reaching far beyond ribosomal synthesis. Merely the copy number of non-transcribed rRNA genes is involved in mechanisms governing e.g., maintenance of genome integrity and control of cellular aging. Their copy number can vary in response to environmental cues, in cellular stress sensing, in development of cancer and other diseases. While reaching hundreds of copies in humans, there are records of up to 20,000 copies in fish and frogs and even 400,000 copies in ciliates forming thus a literal subgenome or an rDNAome within the genome. From the compositional and evolutionary dynamics viewpoint, the precursor 45S rDNA represents universally GC-enriched, highly recombining and homogenized regions. Hence, it is not accidental that both rDNA sequence and the corresponding rRNA secondary structure belong to established phylogenetic markers broadly used to infer phylogeny on multiple taxonomical levels including species delimitation. However, these multiple roles of rDNAs have been treated and discussed as being separate and independent from each other. Here, I aim to address nuclear rDNAs in an integrative approach to better assess the complexity of rDNA importance in the evolutionary context.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Structural signatures of thermal adaptation of bacterial ribosomal RNA, transfer RNA, and messenger RNA

Temperature adaptation of bacterial RNAs is a subject of both fundamental and practical interest because it will allow a better understanding of molecular mechanism of RNA folding with potential industrial application of functional thermophilic or psychrophilic RNAs. Here, we performed a comprehensive study of rRNA, tRNA, and mRNA of more than 200 bacterial species with optimal growth temperatures (OGT) ranging from 4°C to 95°C. We investigated temperature adaptation at primary, secondary and tertiary structure levels. We showed that unlike mRNA, tRNA and rRNA were optimized for their structures at compositional levels with significant tertiary structural features even for their corresponding randomly permutated sequences. tRNA and rRNA are more exposed to solvent but remain structured for hyperthermophiles with nearly OGT-independent fluctuation of solvent accessible surface area within a single RNA chain. mRNA in hyperthermophiles is essentially the same as random sequences without tertiary structures although many mRNA in mesophiles and psychrophiles have well-defined tertiary structures based on their low overall solvent exposure with clear separation of deeply buried from partly exposed bases as in tRNA and rRNA. These results provide new insight into temperature adaptation of different RNAs.

The survival mechanisms of thermophiles at high temperatures: an angle of omics

Thermophiles are referred to as microorganisms with optimal growth temperatures of >60 °C. Over the past few years, a number of studies have been conducted regarding thermophiles, especially using the omics strategies. This review provides a systematic view of the survival physiology of thermophiles from an "omics" perspective, which suggests that the adaptive ability of thermophiles is based on a cooperative mode with multi-dimensional regulations integrating genomics, transcriptomics, and proteomics.

Microbial diversity in a submarine carbonate edifice from the serpentinizing hydrothermal system of the Prony Bay (New Caledonia) over a 6-year period

Active carbonate chimneys from the shallow marine serpentinizing Prony Hydrothermal Field were sampled 3 times over a 6 years period at site ST09. Archaeal and bacterial communities composition was investigated using PCR-based methods (clone libraries, Denaturating Gel Gradient Electrophoresis, quantitative PCR) targeting 16S rRNA genes, methyl coenzyme M reductase A and dissimilatory sulfite reductase subunit B genes. Methanosarcinales (Euryarchaeota) and Thaumarchaea were the main archaeal members. The Methanosarcinales, also observed by epifluorescent microscopy and FISH, consisted of two phylotypes that were previously solely detected in two other serpentinitzing ecosystems (The Cedars and Lost City Hydrothermal Field). Surprisingly, members of the hyperthermophilic order Thermococcales were also found which may indicate the presence of a hot subsurface biosphere. The bacterial community mainly consisted of Firmicutes, Chloroflexi, Alpha-, Gamma-, Beta-, and Delta-proteobacteria and of the candidate division NPL-UPA2. Members of these taxa were consistently found each year and may therefore represent a stable core of the indigenous bacterial community of the PHF chimneys. Firmicutes isolates representing new bacterial taxa were obtained by cultivation under anaerobic conditions. Our study revealed diverse microbial communities in PHF ST09 related to methane and sulfur compounds that share common populations with other terrestrial or submarine serpentinizing ecosystems.

The ancient heart of the ribosomal large subunit: a response to Caetano-Anolles

Our recent Accretion Model of ribosomal evolution uses insertion fingerprints and a "trunk-branch" formalism to recapitulate the building up of common core rRNA of the Large Ribosomal Subunit. The Accretion Model is a conservative and natural extension of a method developed by Bokov and Steinberg (Nature 457:977-80, 2009), which confirms the correctness of lower resolution models by Fox and others. In each of these models, the LSU originates with the peptidyl transferase center (PTC), consistent with expectations that the ribosome is the source of defined-sequence functional proteins. In an adjacent note, Caetano-Anolles (J Mol Evol 80:162-165, 2015) disparages the Accretion Model, because it controverts the 'Growth Inferred by Genothermal Ordering' (GIGO) model. GIGO analyzes secondary structures, assigns the origin of the ribosome to a region outside of the PTC, and assumes or deduces that (i) large protein enzymes of defined amino acid sequence predate ribosomal synthesis of proteins, (ii) proteins directly replicate by non-ribosomal mechanisms, (iii) rRNA unfailingly increases in thermodynamic stability over time, and (iv) the Woese and Fox canonical tree of life is mis-rooted. Much of the specific GIGO critique of the Accretion Model is based on confusion about the three-dimensional nature of RNA and trunk-branch polymorphism; the Accretion Model incorporates several types of trunk-branch relationships.

The GC content of LSU rRNA evolves across topological and functional regions of the ribosome in all three domains of life

Large-subunit rRNA is the ribozyme that catalyzes protein synthesis by translation, and many of its features vary along a deep-to-superficial gradient. By measuring the G+C proportions in this rRNA in all three domains of life (60 bacteria, 379 eukaryote, and 23 archaean sequences), we tested whether the proportion of GC nucleotides varies along this in-out gradient. The rRNA regions used were several zones identified by Bokov and Steinberg (2009) as being arranged from deep to superficial within the LSU. To the Bokov-Steinberg zones, we added the most superficial zone of all, the divergent domains (expansion segments), which are greatly enlarged in eukaryotes. Regression lines constructed from the hundreds of species of organisms revealed the expected in-out gradient, showing that species with high %GC (or high %AT) in their rRNA distribute more of these abundant nucleotides into the peripheral zones. This could be explained by the evolutionary rates of replacement of all nucleotides (A, C, G, T), because these latter rates are fastest at the periphery and slowest near the conserved core. As an overall explanation, we propose that when extrinsic factors (whole-genome nucleotide composition, or environmental temperature) demand the percentage of GC in the rRNA of a species be high or low, then the deep-lying zones are buffered against GC variation because they are the slowest to evolve. The deep, conserved zones are also the most involved in translation, hinting that stabilizing selection there prevents a high GC variability that would diminish LSU rRNA's core functions. We found only a few domain-specific trends in rRNA-GC distribution, which relate to many Archaea living at high temperatures or to the highly complex genes and adaptations of Eukaryota. Use of rRNA sequences in molecular phylogenetic studies, for reconstructing the relationships of organisms across the tree of life, requires accurate models of how rRNA evolves. The demonstration that GC distributes in regular patterns across rRNA regions can improve these tree-reconstruction models in the future and should yield phylogenies of greater accuracy.

A comparative study of ribosomal proteins: linkage between amino acid distribution and ribosomal assembly

BACKGROUND

Assembly of the ribosome from its protein and RNA constituents must occur quickly and efficiently in order to synthesize the proteins necessary for all cellular activity. Since the early 1960's, certain characteristics of possible assembly pathways have been elucidated, yet the mechanisms that govern the precise recognition events remain unclear.We utilize a comparative analysis to investigate the amino acid composition of ribosomal proteins (r-proteins) with respect to their role in the assembly process. We compared small subunit (30S) r-protein sequences to those of other housekeeping proteins from 560 bacterial species and searched for correlations between r-protein amino acid content and factors such as assembly binding order, environmental growth temperature, protein size, and contact with ribosomal RNA (rRNA) in the 30S complex.

RESULTS

We find r-proteins have a significantly high percent of positive residues, which are highly represented at rRNA contact sites. An inverse correlation between the percent of positive residues and r-protein size was identified and is mainly due to the content of Lysine residues, rather than Arginine. Nearly all r-proteins carry a net positive charge, but no statistical correlation between the net charge and the binding order was detected. Thermophilic (high-temperature) r-proteins contain increased Arginine, Isoleucine, and Tyrosine, and decreased Serine and Threonine compared to mesophilic (lower-temperature), reflecting a known distinction between thermophiles and mesophiles, possibly to account for protein thermostability. However, this difference in amino acid content does not extend to rRNA contact sites, as the proportions of thermophilic and mesophilic contact residues are not significantly different.

CONCLUSIONS

Given the significantly higher level of positively charged residues in r-proteins and at contact sites, we conclude that ribosome assembly relies heavily on an electrostatic component of interaction. However, the binding order of r-proteins in assembly does not appear to depend on these electrostatics interactions. Additionally, because thermophiles and mesophiles exhibit significantly different amino acid compositions in their sequences but not in the identities of contact sites, we conclude that this electrostatic component of interaction is insensitive to temperature and is not the determining factor differentiating the temperature sensitivity of ribosome assembly.

Growth temperatures of archaeal communities can be estimated from the guanine-plus-cytosine contents of 16S rRNA gene fragments

Prokaryote growth temperatures in environmental samples are difficult to measure because it is hard to culture viable prokaryotes in natural environments. We comprehensively surveyed growth temperatures and 16S rRNA sequences of prokaryotes to estimate their growth temperatures based on guanine-plus-cytosine contents (P(GC)) of their 16S rRNA sequences. We focused on archaea because of the wide range of growth temperatures within this group. Their minimum (Tmin), optimum (Topt) and maximum (Tmax) growth temperatures correlated strongly with PGC of their 16S rRNA genes. Linear regression equations were established to approximate Tmin, Topt and Tmax from P(GC). We also established a linear regression equation for calculating P(GC) of 16S rRNA genes based on the melting temperatures (Tm) of PCR fragments, without using a clone library or sequencing. Environmental samples were obtained from a wide variety of microbial natural habitats. Tm of archaeal 16S rRNA genes amplified by real-time PCR were determined by melting curve analysis. Based on those values, P(GC) of 16S rRNA genes and mean Tmin, Topt and Tmax were calculated using the linear regression equations. These temperatures correlated strongly with the in situ temperatures. Tmax agreed particularly well with these temperatures, suggesting many archaea live at their maximum growth temperatures.

DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution

Since its first release in 2001 as mainly a software package for phylogenetic analysis, data analysis for molecular biology and evolution (DAMBE) has gained many new functions that may be classified into six categories: 1) sequence retrieval, editing, manipulation, and conversion among more than 20 standard sequence formats including MEGA, NEXUS, PHYLIP, GenBank, and the new NeXML format for interoperability, 2) motif characterization and discovery functions such as position weight matrix and Gibbs sampler, 3) descriptive genomic analysis tools with improved versions of codon adaptation index, effective number of codons, protein isoelectric point profiling, RNA and protein secondary structure prediction and calculation of minimum folding energy, and genomic skew plots with optimized window size, 4) molecular phylogenetics including sequence alignment, testing substitution saturation, distance-based, maximum parsimony, and maximum-likelihood methods for tree reconstructions, testing the molecular clock hypothesis with either a phylogeny or with relative-rate tests, dating gene duplication and speciation events, choosing the best-fit substitution models, and estimating rate heterogeneity over sites, 5) phylogeny-based comparative methods for continuous and discrete variables, and 6) graphic functions including secondary structure display, optimized skew plot, hydrophobicity plot, and many other plots of amino acid properties along a protein sequence, tree display and drawing by dragging nodes to each other, and visual searching of the maximum parsimony tree. DAMBE features a graphic, user-friendly, and intuitive interface and is freely available from http://dambe.bio.uottawa.ca (last accessed April 16, 2013).

Evolutionary force in confamiliar marine vertebrates of different temperature realms: adaptive trends in zoarcid fish transcriptomes

BACKGROUND

Studies of temperature-induced adaptation on the basis of genomic sequence data were mainly done in extremophiles. Although the general hypothesis of an increased molecular flexibility in the cold is widely accepted, the results of thermal adaptation are still difficult to detect at proteomic down to the genomic sequence level. Approaches towards a more detailed picture emerge with the advent of new sequencing technologies. Only small changes in primary protein structure have been shown to modify kinetic and thermal properties of enzymes, but likewise for interspecies comparisons a high genetic identity is still essential to specify common principles. The present study uses comprehensive transcriptomic sequence information to uncover general patterns of thermal adaptation on the RNA as well as protein primary structure.

RESULTS

By comparing orthologous sequences of two closely related zoarcid fish inhabiting different latitudinal zones (Antarctica: Pachycara brachycephalum, temperate zone: Zoarces viviparus) we were able to detect significant differences in the codon usage. In the cold-adapted species a lower GC content in the wobble position prevailed for preserved amino acids. We were able to estimate 40-60% coverage of the functions represented within the two compared zoarcid cDNA-libraries on the basis of a reference genome of the phylogenetically closely related fish Gasterosteus aculeatus. A distinct pattern of amino acid substitutions could be identified for the non-synonymous codon exchanges, with a remarkable surplus of serine and reduction of glutamic acid and asparagine for the Antarctic species.

CONCLUSION

Based on the differences between orthologous sequences from confamiliar species, distinguished mainly by the temperature regimes of their habitats, we hypothesize that temperature leaves a signature on the composition of biological macromolecules (RNA, proteins) with implications for the transcription and translation level. As the observed pattern of amino acid substitutions only partly support the flexibility hypothesis further evolutionary forces may be effective at the global transcriptome level.

Amino acid composition in endothermic vertebrates is biased in the same direction as in thermophilic prokaryotes

Background

Among bacteria and archaea, amino acid usage is correlated with habitat temperatures. In particular, protein surfaces in species thriving at higher temperatures appear to be enriched in amino acids that stabilize protein structure and depleted in amino acids that decrease thermostability. Does this observation reflect a causal relationship, or could the apparent trend be caused by phylogenetic relatedness among sampled organisms living at different temperatures? And do proteins from endothermic and exothermic vertebrates show similar differences?

Results

We find that the observed correlations between the frequencies of individual amino acids and prokaryotic habitat temperature are strongly influenced by evolutionary relatedness between the species analysed; however, a proteome-wide bias towards increased thermostability remains after controlling for phylogeny. Do eukaryotes show similar effects of thermal adaptation? A small shift of amino acid usage in the expected direction is observed in endothermic ('warm-blooded') mammals and chicken compared to ectothermic ('cold-blooded') vertebrates with lower body temperatures; this shift is not simply explained by nucleotide usage biases.

Conclusion

Protein homologs operating at different temperatures have different amino acid composition, both in prokaryotes and in vertebrates. Thus, during the transition from ectothermic to endothermic life styles, the ancestors of mammals and of birds may have experienced weak genome-wide positive selection to increase the thermostability of their proteins.

Compositional properties and thermal adaptation of SRP-RNA in bacteria and archaea

Previous studies have reported a positive correlation between the GC content of the double-stranded regions of structural RNAs and the optimal growth temperature (OGT) in prokaryotes. These observations led to the hypothesis that natural selection favors an increase in GC content to ensure the correct folding and the structural stability of the molecule at high temperature. To date these studies have focused mainly on ribosomal and transfer RNAs. Therefore, we addressed the question of the relationship between GC content and OGT in a different and universally conserved structural RNA, the RNA component of the signal recognition particle (SRP). To this end we generated the secondary structures of SRP-RNAs for mesophilic, thermophilic, and hyperthermophilic bacterial and archaeal species. The analysis of the GC content in the stems and loops of the SRP-RNA of these organisms failed to detect a relationship between the GC contents in the stems of this structural RNA and the growth temperature of bacteria. By contrast, we found that in archaea the GC content in the stem regions of SRP-RNA is highest in hyperthermophiles, intermediate in thermophiles, and lower in mesophiles. In these organisms, we demonstrated a clear positive correlation between the GC content of the stem regions of their SRP-RNAs and their OGT. This correlation was confirmed by a phylogenetic nonindependence analysis. Thus we conclude that in archaea the increase in GC content in the stem regions of SRP-RNA is an adaptation response to environmental temperature.

Environmental shaping of ribosomal RNA nucleotide composition

Ribosomal RNA (rRNA) is one of the most important macromolecules in the cell. It is well established that high-temperature environmental conditions destabilize rRNA, leading to a selection for G+C-rich stabilizing structures. Our knowledge about the nucleotide composition effect of other environmental conditions, however, is limited. In the present work, I addressed this by correlating the rRNA nucleotide composition to known environmental habitats for bacteria. The bacterial phyla Firmicutes, Actinobacteria, and Proteobacteria were chosen for in-depth analyses due to the abundance of information available in the databases. Major differences in nucleotide composition were identified between these phyla. In addition to the G+C-->A+T gradients, a main gradient of G+A-->C+T was identified for Firmicutes, while a G+T-->A+C gradient was identified for Actinobacteria. With respect to correlation to environmental conditions, the Firmicutes showed a main structure of high G+C being correlated to thermophilic conditions, high A+T to anaerobic conditions, and high C+T to halophilic conditions. The main patterns detected for Firmicutes can be explained by structural stability for high G+C, chemical instability of G under aerobic conditions, and structural stability by purine/pyrimidine skew for halophilic conditions. On the contrary, the correlations for Actinobacteria cannot easily be explained by chemical and/or structural stability. This may indicate interference with factors not included in my work. Finally, I found a main correlation between high A+T and endosymbiosis for Proteobacteria. High A+T probably reflects adaptation to cell internal growth. Further support for environmentally driven nucleotide composition shaping was found and that polyphyletic bacteria were associated with the same environment/nucleotide correlations. My conclusion is that environmental conditions and habitats have a major effect on rRNA nucleotide composition but that the effects may differ between the bacterial phyla.

Compositional properties and thermal adaptation of 18S rRNA in vertebrates

In order to investigate the influence of temperature on the GC level of the paired sequences of ribosomal 18S RNAs in vertebrates, we have studied their base composition in cold- and warm-blooded vertebrates using a stem-by-stem comparison. We observed that a number of stems of 18S ribosomal RNAs (rRNAs) are variable among species and that the majority of such stems are GC richer in warm-blooded than in cold-blooded vertebrates. We also constructed the secondary structures of the 18S rRNAs of a polar fish, a marsupial, and a monotreme to compare them with those of temperate/tropical fishes and of eutherians, respectively. In these cases, differences similar to those already mentioned were found. We conclude that there is a correlation between stem stability and body temperature even within the relatively limited temperature range of vertebrates.

Hyperthermophilic enzymes − stability, activity and implementation strategies for high temperature applications

Current theories agree that there appears to be no unique feature responsible for the remarkable heat stability properties of hyperthermostable proteins. A concerted action of structural, dynamic and other physicochemical attributes are utilized to ensure the delicate balance between stability and functionality of proteins at high temperatures. We have thoroughly screened the literature for hyperthermostable enzymes with optimal temperatures exceeding 100 degrees C that can potentially be employed in multiple biotechnological and industrial applications and to substitute traditionally used, high-cost engineered mesophilic/thermophilic enzymes that operate at lower temperatures. Furthermore, we discuss general methods of enzyme immobilization and suggest specific strategies to improve thermal stability, activity and durability of hyperthermophilic enzymes.

The neoselectionist theory of genome evolution

The vertebrate genome is a mosaic of GC-poor and GC-rich isochores, megabase-sized DNA regions of fairly homogeneous base composition that differ in relative amount, gene density, gene expression, replication timing, and recombination frequency. At the emergence of warm-blooded vertebrates, the gene-rich, moderately GC-rich isochores of the cold-blooded ancestors underwent a GC increase. This increase was similar in mammals and birds and was maintained during the evolution of mammalian and avian orders. Neither the GC increase nor its conservation can be accounted for by the random fixation of neutral or nearly neutral single-nucleotide changes (i.e., the vast majority of nucleotide substitutions) or by a biased gene conversion process occurring at random genome locations. Both phenomena can be explained, however, by the neoselectionist theory of genome evolution that is presented here. This theory fully accepts Ohta's nearly neutral view of point mutations but proposes in addition (i) that the AT-biased mutational input present in vertebrates pushes some DNA regions below a certain GC threshold; (ii) that these lower GC levels cause regional changes in chromatin structure that lead to deleterious effects on replication and transcription; and (iii) that the carriers of these changes undergo negative (purifying) selection, the final result being a compositional conservation of the original isochore pattern in the surviving population. Negative selection may also largely explain the GC increase accompanying the emergence of warm-blooded vertebrates. In conclusion, the neoselectionist theory not only provides a solution to the neutralist/selectionist debate but also introduces an epigenomic component in genome evolution.

Bushes in the tree of life

A discussion of how homoplasy (the frequency of independently evolved characters) and the spacing of cladogenetic events limit our ability to reconstruct the tree of life using existing phylogenetic methods.