The lineage-specific base-pair contents in the stem regions of ribosomal RNAs and their influence on the estimation of evolutionary distances

The base-pair changes in the stem regions of ribosomal RNAs provide a useful measure for resolving the phylogeny of organisms. In the present study, how the biased base-pair content influences the estimation of evolutionary distances is theoretically investigated. By regarding the biased base-pair content as a result of the difference in selective strength between A:U and G:C base pairs, the evolutionary distance empirically obtained by enumerating base-pair changes is theoretically expressed in terms of selective strength, base-pair change rate, and divergence time. Its application to nuclear-coded large subunit ribosomal RNAs (LSU rRNAs) reveals the followings. LSU rRNAs from most organisms have moderate base-pair contents and the empirical evolutionary distances obtained by the comparison of these LSU rRNAs are approximately proportional to their divergence times. In the comparison of these moderate LSU rRNAs with the GC-rich LSU rRNAs such as those from Mycoplasma, Crenarchaeota, and Giardia, however, the empirically calculated distances are considerably smaller than the true evolutionary distances, while the comparison with AU-rich LSU rRNAs from Microsporidia overestimates their distances. With this result in mind, the relative base-pair change probabilities among three kingdoms are carefully estimated from the statistical distribution of base-pair change ratios enumerated for LSU rRNAs showing almost the same base-pair contents, leading to the result that prokaryotes and eukaryotes first diverged and that archaebacteria and eubacteria diverged on the line of prokaryotes slightly later, by about 0.3 billion years.     

Rabbit liver tRNAVal1: II. Unusual secondary structure of T ψ C stem and loop due to a U54:A60 base pair

In contrast to all other known tRNAs, mammalian tRNA^Val₁ contains two adenosines A₅₉ and A₆₀, opposite to U₅₄ and ψ₅₅ in the UψCG sequence of the TψC loop, which could form unusual A:U (or A:ψ) pairs in addition to the five “normal” G:C pairs. In order to measure the number of G:C and A:U (A:ψ) pairs in the TψC stem, we prepared the 30 nucleotide long 3′-terminal fragment of this tRNA by “m⁷G-cleavage”. From differentiated melting curves and temperature jump experiments it was concluded that the TψC stem in this fragment is in fact extended by an additional A₆₀:U₅₄ pair. A dimer of this fragment with 14 base pairs was characterized by gel electrophoresis and by the same physical methods. An additional A:U pair in the tRNA^Val₁ fragment does not necessarily mean that this is also true for intact tRNA. However, we showed that U₅₄ is far less available for enzymatic methylation in mammalian tRNA^Val₁ compared to tRNA from T⁻E. coli. This clear difference in U₅₄ reactivity, together with the identification of an extra A₆₀:U₅₄ pair in the UψCG containing fragment suggests the presence of a 6 base pair TψC stem and a 5 nucleotide TψC loop in this tRNA.     

The G x U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems

The G x U wobble base pair is a fundamental unit of RNA secondary structure that is present in nearly every class of RNA from organisms of all three phylogenetic domains. It has comparable thermodynamic stability to Watson-Crick base pairs and is nearly isomorphic to them. Therefore, it often substitutes for G x C or A x U base pairs. The G x U wobble base pair also has unique chemical, structural, dynamic and ligand-binding properties, which can only be partially mimicked by Watson-Crick base pairs or other mispairs. These features mark sites containing G x U pairs for recognition by proteins and other RNAs and allow the wobble pair to play essential functional roles in a remarkably wide range of biological processes.     

Molecular evolution of phosphoenolpyruvate carboxylase for C4 photosynthesis in maize: comparison of its cDNA sequence with a newly isolated cDNA encoding an isozyme involved in the anaplerotic function.

Molecular events associated with the evolution of an enzyme for C4 photosynthesis were investigated. In maize, at least three isozymes of phosphoenolpyruvate carboxylase [EC 4.1.1.31] are known: the C4-form, the C2-form and the root-form, being named according to their physiological roles and pattern of tissue distribution [Ting, I.P. & Osmond, C.B. (1973) Plant Physiol. 51, 448-453]. A cDNA clone which presumably encodes the C3-form isozyme was newly isolated and analyzed. Comparison of the sequences of the C3-form and C4-form isozymes revealed that (i) the homologies in the nucleotide and deduced amino acid sequences were 71 and 77%, respectively, and (ii) the gene for the C4-form isozyme evolved under strong G/C pressure. The genes for these isozymes were found to be located apart on different chromosomes. A phylogenetic tree was constructed using 8 amino acid sequences of phosphoenolpyruvate carboxylases from various sources. The topology of the tree indicated that, at least in monocots such as maize and sorghum, the genes for the C4-form and C3-form isozymes diverged from their common ancestral gene earlier than the monocot-dicot divergence (about 2 x 10(8) yr ago), though the divergence of maize (C4 plant) from wheat (C3 plant) is supposed to have occurred much later (6 x 10(7) yr ago).     

Effects of flanking G . C base pairs on internal Watson-Crick, G . U, and nonbonded base pairs within a short ribonucleic acid duplex

A series of pentaribonucleotides, ApGpXpGpU (where X identical to A, G, C, or U), was synthesized to investigate the effects of flanking G . C pairs on internal Watson-Crick, G . U, and nonbonded base pairs. Sequences ApGpApCpU (Tm = 26 degrees C) and ApGpCpCpU (Tm = 25 degrees C) were each found to form a duplex with non-base-paired internal residues that stacked with the rest of the sequence but were not looped out. ApGpGpCpU also forms a duplex (Tm = 30 degrees C) but with dangling terminal nonbonded adenosines rather than internal nonbonded guanosines. ApGpUpCpU prefers a stacked single-strand conformation. In addition, contribution to duplex stability from an internal A . U or G . C base pair is enhanced by 6 degrees C when flanked by G . C base pairs as compared to A . U base pairs. G . C base pairs flanking an internal G . U base pair were found to be more tolerant to the altered conformation of a G . U pair and result in an increase to stability comparable with that found for an internal A . U base pair.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Effects of internal nonbonded bases and a G.U base pair on the stability of a short ribonucleic acid helix.

Proton nuclear magnetic resonance has been used to examine the effect of both noncomplementary and G.U oppositions in the duplexes formed by the synthetic pentaribonucleotides CpApApUpG, CpApUpUpG, CpApGpUpG, and CpApCpUpG. The lack of any sigmoidal behavior in the chemical shift vs. temperature plots of the base protons in the individual pentaribonucleotides indicates that duplexes with noncomplementary base oppositions of the type: formula: (see text), (where X = A, U, G, or C) do not form. Variable temperature spectra of the mixture of CpApGpUpG and CpApUpUpG were recorded over the range of 70--10 degrees C. The chemical shift vs. temperature plot of the purine aromatic protons displayed sigmoidal curves. This demonstrated both duplex formation and the presence of a G.U. base pair. The average Tm of the duplex was found to be 23.4 +/- 2.0 degrees C. This is similar to that of the duplex formed by CpApUpG (24.0 +/- 1.0 degrees C) but less than the Tm of the following duplexes: CpApApUpG:CpApUpUpG (Tm = 28.5 +/- 2.1 degrees C), CpApGpUpG:CpApCpUpG (Tm = 38.4 +/- 0.6 degrees C) and CpApUpApUpG (Tm = 41.5 +/- 1.1 degrees C). The G.U base pair has a Tm (20.0 degrees C) significantly lower than the rest of the duplex (24 +/- 1 degree C) and is a region of local instability within the double helix. This 1H NMR study is the first to investigate both the formation and relative stability of an internal G.U. base pair neighboring regular Watson--Crick base pairs.     

Synthesis and hybridization of 2′-O-methyl-RNAs incorporating 2′-O-carbamoyluridine and unique participation of the carbamoyl group in U–G base pair

2′-O-Carbamoyluridine (U_cm) was synthesized and incorporated into DNAs and 2′-O-Me-RNAs. The oligonucleotides incorporating U_cm formed less stable duplexes with their complementary and U_cm–U, U_cm–C single-base mismatched DNAs and RNAs in comparison with those without the carbamoyl group. On the contrary, the T_m analyses revealed that the duplexes with a mismatched U_cm–G base pair showed almost the same thermostability as the corresponding unmodified duplexes. Molecular dynamics (MD) simulations of the U_cm-modified 2′-O-Me-RNA/RNA duplexes with U_cm–G mismatched base pair suggested that the carbamoyl group could participate in the U_cm–G base pair by an additional intermolecular hydrogen bond between the carbamoyl oxygen and the H2 of the guanine base.     

Experimental and theoretical rationalization for the base pairing abilities of inosine,guanosine, adenosine,and their corresponding 8-oxo-7,8-dihydropurine,and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (T_m) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.     

The influence of base identity and base pairing on the function of the alpha-sarcin loop of 23S rRNA.

The alpha-sarcin loop of large subunit rRNAs is one of the sites of interaction of elongation factors with the ribosome, and the target of the cytotoxins alpha-sarcin and ricin. Using a genetic selection for increased frameshifting in a reporter gene, we have isolated a C --> U mutation at position 2666 in the alpha-sarcin loop. In the NMR-derived structure of the loop, bases equivalent to 2666 and 2654 are paired via a non-canonical base pairing interaction. Each of the three base substitutions at C2666 and A2654 was constructed by site-directed mutagenesis of a plasmid borne copy of the rrnB operon of Escherichia coli. Only the C2666 --> U and A2654 --> G mutations that resulted in the formation of canonical A-U and C-G base pairs respectively, increased the levels of stop codon readthrough and frameshifting. The effects of different base pair combinations at positions 2666 and 2654 on ribosome function were then tested by constructing and analyzing all possible base combinations at these sites. All A --> G base substitution mutations at position 2654 and C --> U substitutions at position 2666 increased the levels of translational errors. However, these effects were greatest when G2654 and U2666 had the potential to engage in standard Watson-Crick base pairing interactions. These data indicate that base identity as well as base pairing interactions are important for the function of this essential component of the large subunit rRNA.     

Evolution of the FAD2-1 fatty acid desaturase 5' UTR intron and the molecular systematics of Gossypium (Malvaceae)

The FAD2-1 microsomal omega-6 desaturase gene contains a large intron ( approximately 1133 bp [base pairs]) in the 5' untranslated region that may participate in gene regulation and, in GOSSYPIUM:, is evolving at an evolutionary rate useful for elucidating recently diverged lineages. FAD2-1 is single copy in diploid GOSSYPIUM: species, and two orthologs are present in the allotetraploid species. Among the diploid species, the D-genome FAD2-1 introns have accumulated substitutions 1.4-1.8 times faster than the A-genome introns. In the tetraploids, the difference between the D-subgenome introns and their A-subgenome orthologs is even greater. The substitution rate of the intron in the D-genome diploid G. gossypioides more closely approximates that of the A genome than other D genome species, highlighting its unique evolutionary history. However, phylogenetic analyses support G. raimondii as the closest living relative of the D-subgenome donor. The Australian K-genome species diverged 8-16 million years ago into two clades. One clade comprises the sporadically distributed, erect to suberect coastal species; a second clade comprises the more widely spread, prostrate, inland species. A comparison of published gene trees to the FAD2-1 intron topology suggests that G. bickii arose from an early divergence, but that it carries a G. australe-like rDNA captured via a previously undetected hybridization event.     

Crystal structure of human selenocysteine tRNA

Selenocysteine (Sec) is the 21st amino acid in translation. Sec tRNA (tRNA^Sec) has an anticodon complementary to the UGA codon. We solved the crystal structure of human tRNA^Sec. tRNA^Sec has a 9-bp acceptor stem and a 4-bp T stem, in contrast with the 7-bp acceptor stem and the 5-bp T stem in the canonical tRNAs. The acceptor stem is kinked between the U6:U67 and G7:C66 base pairs, leading to a bent acceptor-T stem helix. tRNA^Sec has a 6-bp D stem and a 4-nt D loop. The long D stem includes unique A14:U21 and G15:C20a pairs. The D-loop:T-loop interactions include the base pairs G18:U55 and U16:U59, and a unique base triple, U20:G19:C56. The extra arm comprises of a 6-bp stem and a 4-nt loop. Remarkably, the D stem and the extra arm do not form tertiary interactions in tRNA^Sec. Instead, tRNA^Sec has an open cavity, in place of the tertiary core of a canonical tRNA. The linker residues, A8 and U9, connecting the acceptor and D stems, are not involved in tertiary base pairing. Instead, U9 is stacked on the first base pair of the extra arm. These features might allow tRNA^Sec to be the target of the Sec synthesis/incorporation machineries.     

Conformational features of the four successive non-Watson-Crick base pairs in RNA duplex.

A tridecaribonucleotide, r(UGAGCUUCGGCUC) doesn't form hairpin or interior loop and forms a double helix of 12 base pairs including the four successive nonstandard base pairs, U.G-U.C-C.U-G.U, in the crystal. Non-Watson-Crick base pairs, G.U and U.C are nicely incorporated in RNA duplex maintaining the regular A-form backbone. There exist the good overlapping between base pairings, U.G and U.C, so as to stabilize the nonstandard base pair track. Hydrogen bond networks involving water molecules in the major and minor grooves to stabilize this mismatch base pairing array, are observed and its conformational features are described.     

Rabbit liver tRNA1Val:II. unusual secondary structure of T psi C stem and loop due to a U54:A60 base pair.

In contrast to all other known tRNAs, mammalian tRNA1Val contains two adenosines A59 and A60, opposite to U54 and psi 55 in the U psi CG sequence of the T psi C loop, which could form unusual A:U (or A: psi pairs in addition to the five "normal" G:C pairs. In order to measure the number of G:C and A:U (A: psi) pairs in the T psi C stem, we prepared the 30 nucleotide long 3'-terminal fragment of this tRNA by "m7G-cleavage". From differentiated melting curves and temperature jump experiments it was concluded that the T psi C stem in this fragment is in fact extended by an additional A60:U54 pair. A dimer of this fragment with 14 base pairs was characterized by gel electrophoresis and by the same physical methods. An additional A:U pair in the tRNA1Val fragment does not necessarily mean that this is also true for intact tRNA. However, we showed that U54 is far less available for enzymatic methylation in mammalian tRNA1Val compared to tRNA from T-E. coli. This clear difference in U54 reactivity, together with the identification of an extra A60:U54 pair in the U psi CG containing fragment suggests the presence of a 6 base pair T psi C stem and a 5 nucleotide T psi C loop in this tRNA.     

Closing of the Tethys Sea and the phylogeny of Eurasian killifishes (Cyprinodontiformes: Cyprinodontidae)

To test vicariant speciation hypotheses derived from geological evidence of the closing of the Tethys Sea, we reconstruct phylogenetic relationships of the predominantly fresh-water killifish genus Aphanius using 3263 aligned base pairs of mitochondrial DNA from samples representing 49 populations of 13 species. We use additional 11 cyprinodontid species as outgroup taxa. Genes analysed include those encoding the partial 12S and 16S ribosomal RNAs; transfer RNAs for valine, leucine, isoleucine, glutamine, methionine, tryptophan, alanine, asparagine, cysteine and tyrosine; and complete nicotinamide adenine dinucleotide dehydrogenase subunit I and II. Molecular substitution rate for this DNA region is estimated at of 8.6 +/- 0.1 x 10(-9) substitutions base pair(-1) year(-1), and is derived from a well dated transgression of the Red Sea into the Wadi Sirhan of Jordan 13 million years ago; an alternate substitution rate of 1.1 +/- 0.2 x 10(-8) substitutions base pair(-1) year(-1) is estimated from fossil evidence. Aphanius forms two major clades which correspond to the former eastern and western Tethys Sea. Within the eastern clade Oligocene divergence into a fresh-water clade inhabiting the Arabian Peninsula and an euhaline clade inhabiting coastal area from Pakistan to Somalia is observed. Within the western Tethys Sea clade we observe a middle Oligocene divergence into Iberian Peninsula and Atlas Mountains, and Turkey and Iran sections. Within Turkey we observe a large amount of genetic differentiation correlated with late Miocene orogenic events. Based on concordance of patterns of phylogenetic relationships and area relationships derived from geological and fossil data, as well as temporal congruence of these patterns, we support a predominantly vicariant-based speciation hypothesis for the genus Aphanius. An exception to this pattern forms the main clade of A. fasciatus, an euhaline circum-Mediterranean species, which shows little genetic differentiation or population structuring, thus providing no support for the hypothesis of vicariant differentiation associated with the Messinian Salinity Crisis. The two phylogenetically deepest events were also likely driven by ecological changes associated with the closing of the Tethys Sea.     

An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA

A previous NMR structure of the duplex 5'GGU GGA GGCU/PCCG AAG CCG5' revealed an unusually stable RNA internal loop with three consecutive sheared GA pairs. Here, we report NMR studies of two duplexes, 5'GGU GGA GGCU/PCCA AAG CCG5' (replacing the UG pair with a UA closing pair) and 5'GGU GAA GGCU/PCCG AAG CCG5' (replacing the middle GA pair with an AA pair). An unusually stable loop with three consecutive sheared GA pairs forms in the duplex 5'GGU GGA GGCU/PCCA AAG CCG5'. The structure contrasts with that reported for this loop in the crystal structure of the large ribosomal subunit of Deinococcus radiodurans [Harms, J., Schluenzen, F., Zarivach, R., Bashan, A., Gat, S., Agmon, I., Bartels, H., Franceschi, F., and Yonath, A. (2001) Cell 107, 679-688]. The middle AA pair in the duplex 5'GGU GAA GGCU/PCCG AAG CCG5' rapidly exchanges orientations, resulting in alternative base stacking and pseudosymmetry with exclusively sheared pairs. The U GAA G/G AAG C internal loop is 2.1 kcal/mol less stable than the U GGA G/G AAG C internal loop at 37 degrees C. Structural, energetic, and dynamic consequences upon functional group substitutions within related 3 x 3 and 3 x 6 internal loops are also reported.     

Pleistocene climatic fluctuations explain the disjunct distribution and complex phylogeographic structure of the Southern Red-backed Salamander, <Emphasis Type="Italic">Plethodon serratus</Emphasis>

The southeastern United States (U.S.) has experienced dynamic climatic changes over the past several million years that have impacted species distributions. In many cases, contiguous ranges were fragmented and a lack of gene flow between allopatric populations led to genetic divergence and speciation. The Southern Red-backed Salamander, Plethodon serratus, inhabits four widely disjunct regions of the southeastern U.S.: the southern Appalachian Mountains, the Ozark Plateau, the Ouachita Mountains, and the Southern Tertiary Uplands of central Louisiana. We integrated phylogenetic analysis of mitochondrial DNA sequences (1399 base pairs) with ecological niche modeling to test the hypothesis that climate fluctuations during the Pleistocene drove the isolation and divergence of disjunct populations of P. serratus. Appalachian, Ozark, and Louisiana populations each formed well-supported clades in our phylogeny. Ouachita Mountain populations sorted into two geographically distinct clades; one Ouachita clade was sister to the Louisiana clade whereas the other Ouachita clade grouped with the Appalachian and Ozark clades but relationships were unresolved. Plethodon serratus diverged from its sister taxon, P. sherando, ~5.4 million years ago (Ma), and lineage diversification within P. serratus occurred ~1.9–0.6 Ma (Pleistocene). Ecological niche models showed that the four geographic isolates of P. serratus are currently separated by unsuitable habitat, but the species was likely more continuously distributed during the colder climates of the Pleistocene. Our results support the hypothesis that climate-induced environmental changes during the Pleistocene played a dominant role in driving isolation and divergence of disjunct populations of P. serratus.     

Origin and evolution of influenza virus hemagglutinin genes

Influenza A, B, and C viruses are the etiological agents of influenza. Hemagglutinin (HA) is the major envelope glycoprotein of influenza A and B viruses, and hemagglutinin-esterase (HE) in influenza C viruses is a protein homologous to HA. Because influenza A virus pandemics in humans appear to occur when new subtypes of HA genes are introduced from aquatic birds that are known to be the natural reservoir of the viruses, an understanding of the origin and evolution of HA genes is of particular importance. We therefore conducted a phylogenetic analysis of HA and HE genes and showed that the influenza A and B virus HA genes diverged much earlier than the divergence between different subtypes of influenza A virus HA genes. The rate of amino acid substitution for A virus HAs from duck, a natural reservoir, was estimated to be 3.19 x 10(-4) per site per year, which was slower than that for human and swine A virus HAs but similar to that for influenza B and C virus HAs (HEs). Using this substitution rate from the duck, we estimated that the divergences between different subtypes of A virus HA genes occurred from several thousand to several hundred years ago. In particular, the earliest divergence time was estimated to be about 2,000 years ago. Also, the A virus HA gene diverged from the B virus HA gene about 4,000 years ago and from the C virus HE gene about 8,000 years ago. These time estimates are much earlier than the previous ones.