Compositional properties and thermal adaptation of 18S rRNA in vertebrates

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.

The "Genomic Code": DNA Pervasively Moulds Chromatin Structures Leaving no Room for "Junk"

The chromatin of the human genome was analyzed at three DNA size levels. At the first, compartment level, two "gene spaces" were found many years ago: A GC-rich, gene-rich "genome core" and a GC-poor, gene-poor "genome desert", the former corresponding to open chromatin centrally located in the interphase nucleus, the latter to closed chromatin located peripherally. This bimodality was later confirmed and extended by the discoveries (1) of LADs, the Lamina-Associated Domains, and InterLADs; (2) of two "spatial compartments", A and B, identified on the basis of chromatin interactions; and (3) of "forests and prairies" characterized by high and low CpG islands densities. Chromatin compartments were shown to be associated with the compositionally different, flat and single- or multi-peak DNA structures of the two, GC-poor and GC-rich, "super-families" of isochores. At the second, sub-compartment, level, chromatin corresponds to flat isochores and to isochore loops (due to compositional DNA gradients) that are susceptible to extrusion. Finally, at the short-sequence level, two sets of sequences, GC-poor and GC-rich, define two different nucleosome spacings, a short one and a long one. In conclusion, chromatin structures are moulded according to a "genomic code" by DNA sequences that pervade the genome and leave no room for "junk".

Comparative cytogenetic of six species of Amazonian Peacock bass (Cichla, Cichlinae): intrachromosomal variations and genetic introgression among sympatric species

Cytogenetic data for the genus Cichla Bloch et Schneider, 1801 are still very limited, with only four karyotype descriptions to date. The sum of the available cytogenetic information for Cichla species, points to a maintenance of the diploid number of 48 acrocentric chromosomes, considered a typical ancestral feature in cichlids. In the current study, we performed molecular and classical cytogenetic analyses of the karyotype organization of six species of Cichla, the earliest-diverging genus of Neotropical cichlids. We cytogenetically analysed Cichla kelberi Kullander et Ferreira, 2006, Cichla monoculus Agassiz, 1831, Cichla piquiti Kullander et Ferreira, 2006, Cichla temensis Humboldt, 1821, Cichla vazzoleri Kullander et Ferreira, 2006 and Cichla pinima Kullander et Ferreira, 2006, including three individuals that showed mixed morphological characteristics, likely from different species, suggesting they were hybrid individuals. All individuals analysed showed 2n = 48 acrocentric chromosomes, with centromeric heterochromatic blocks on all chromosomes and a terminal heterochromatic region on the q arm of the 2^nd pair. Mapping 18S rDNA gave hybridization signals, correlated with the nucleolus organizer regions, on the 2^nd pair for all analyzed individuals. However, we found distinct patterns for 5S rDNA: interstitially at the proximal position on 6^th pair of four species (C. kelberi, C. pinima, C. piquiti and C. vazzoleri), and on the distal of the 4^th pair in two (C. monoculus and C. temensis). Accordingly, we present here new data for the genus and discuss the evolutionary trends in the karyotype of this group of fish. In addition, we provide data that supports the occurrence of hybrid individuals in the Uatumã River region, mainly based on 5S rDNA mapping.

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.

Genomic Organization Under Different Environmental Conditions: Hoplosternum Littorale as a Model

The Amazon has abundant rivers, streams, and floodplains in both polluted and nonpolluted environments, which show great adaptability. Thus, the goal of this study was to map repetitive DNA sequences in both mitotic chromosomes and erythrocyte micronuclei of tamoatás from polluted and nonpolluted environments and to assess the possible genotoxic effects of these environments. Individuals were collected in Manaus, Amazonas (AM), and submitted to classical and molecular cytogenetic techniques, as well as to a blood micronucleus test. Diploid number equal to 60 chromosomes are present in all individuals, with 18S ribosomal DNA sites present in one chromosome pair and no interstitial telomeric sites on chromosomes. The micronucleus test showed no significant differences in pairwise comparisons between environments or collection sites, but the Rex3 retroelement was dispersed on the chromosomes of individuals from unpolluted environments and compartmentalized in individuals from polluted environments. Divergent numbers of 5S rDNA sites are present in individuals from unpolluted and polluted environments. The mapping of repetitive sequences revealed that micronuclei have different compositions both intra- and interindividually that suggests different regions are lost in the formation of micronuclei, and no single fragile region undergoes breaks, although repetitive DNA elements are involved in this process.

Highly divergent 18S rRNA gene paralogs in a Cryptosporidium genotype from eastern chipmunks (Tamias striatus)

Cryptosporidium is an apicomplexan parasite that causes the disease cryptosporidiosis in humans, livestock, and other vertebrates. Much of the knowledge on Cryptosporidium diversity is derived from 18S rRNA gene (18S rDNA) phylogenies. Eukaryote genomes generally have multiple 18S rDNA copies that evolve in concert, which is necessary for the accurate inference of phylogenetic relationships. However, 18S rDNA copies in some genomes evolve by a birth-and-death process that can result in sequence divergence among copies. Most notably, divergent 18S rDNA paralogs in the apicomplexan Plasmodium share only 89-95% sequence similarity, encode structurally distinct rRNA molecules, and are expressed at different life cycle stages. In the present study, Cryptosporidium 18S rDNA was amplified from 28/72 (38.9%) eastern chipmunks (Tamias striatus). Phylogenetic analyses showed the co-occurrence of two 18S rDNA types, Type A and Type B, in 26 chipmunks, and Type B clustered with a sequence previously identified as Cryptosporidium chipmunk genotype II. Types A and B had a sister group relationship but shared less than 93% sequence similarity. In contrast, actin and heat shock protein 70 gene sequences were homogeneous in samples with both Types A and B present. It was therefore concluded that Types A and B are divergent 18S rDNA paralogs in Cryptosporidium chipmunk genotype II. Substitution patterns in Types A and B were consistent with functionally constrained evolution; however, Type B evolved more rapidly than Type A and had a higher G+C content (46.3% versus 41.0%). Oocysts of Cryptosporidium chipmunk genotype II measured 4.17 μm (3.73-5.04 μm) × 3.94 μm (3.50-4.98 μm) with a length-to-width ratio of 1.06 ± 0.06 μm, and infection occurred naturally in the jejunum, cecum, and colon of eastern chipmunks. The findings of this study have implications for the use of 18S rDNA sequences to infer phylogenetic 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.

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.

Variability of 18S rDNA locus among Symphysodon fishes: chromosomal rearrangements

Three species of cichlids belonging to the genus Symphysodon have demonstrated interspecific and intraspecific variation in nucleolus organizer regions (NOR) detected with silver nitrate. In order to understand the evolution of this marker in the genus, the structural variability of these sequences in mitotic chromosomes from Symphysodon aequifasciatus, Symphysodon discus and Symphysodon haraldi was investigated using both silver nitrate impregnation and hybridization of the 18S rRNA gene probe. For the three species, the two markers were intraspecifically and interspecifically variable both in the number and in the size of the sites. This polymorphism may stem from duplications and translocations, which suggests that structural chromosome rearrangements effectively act in the karyoevolution of wild Symphysodon species and may have favoured the adaptability of these fishes to diverse aquatic environments in the Amazon.

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.

The evolution of isochore patterns in vertebrate genomes

Background

Previous work from our laboratory showed that (i) vertebrate genomes are mosaics of isochores, typically megabase-size DNA segments that are fairly homogeneous in base composition; (ii) isochores belong to a small number of families (five in the human genome) characterized by different GC levels; (iii) isochore family patterns are different in fishes/amphibians and mammals/birds, the latter showing GC-rich isochore families that are absent or very scarce in the former; (iv) there are two modes of genome evolution, a conservative one in which isochore patterns basically do not change (e.g., among mammalian orders), and a transitional one, in which they do change (e.g., between amphibians and mammals); and (v) isochores are tightly linked to a number of basic biological properties, such as gene density, gene expression, replication timing and recombination.

Results

The present availability of a number of fully sequenced genomes ranging from fishes to mammals allowed us to carry out investigations that (i) more precisely quantified our previous conclusions; (ii) showed that the different isochore families of vertebrate genomes are largely conserved in GC levels and dinucleotide frequencies, as well as in isochore size; and (iii) isochore family patterns can be either conserved or change within both warm- and cold-blooded vertebrates.

Conclusion

On the basis of the results presented, we propose that (i) the large conservation of GC levels and dinucleotide frequencies may reflect the conservation of chromatin structures; (ii) the conservation of isochore size may be linked to the role played by isochores in chromosome structure and replication; (iii) the formation, the maintainance and the changes of isochore patterns are due to natural selection.