Comparative study of the evolution of nuclear ribosomal spacers incorporating secondary structure analyzes within Dodonaeoideae,Hippocastanoideae and Xanthoceroideae (Sapindaceae)

Ribosomal DNA internal transcribed spacers (ITS) and partial external transcribed spacers (ETSf) are popularly used to infer evolutionary hypotheses. However, there is generally little consideration given to the secondary structures of these small RNA molecules and their potential effects on sequence alignment and phylogenetic analyzes. Intergeneric relationships amongst three of the four major lineages in the Sapindaceae, the Dodonaeoideae, Hippcastanoideae and Xanthoceroideae were assessed by firstly, generating secondary structure predictions for ITS and partial ETSf sequences, and then these predictions were used to assist alignment of the sequences. Secondly, the alignment was analyzed using RNA specific models of sequence evolution that account for the variation in nucleotide evolution in the independent loops and covariating stems regions of the ribosomal spacers. These models and phylogeny drawn from these analyzes were compared with that from analyzes using ‘traditional’ 4-state models and previous plastid analyzes. These analyzes identified that paired-site models developed to deal specifically with stem structures in RNA encoding sequences more appropriately account for the evolutionary history of the sequences than traditional 4-state substitution models.     

Phylogenetic relationships among Myrceugenia, Blepharocalyx, and Luma (Myrtaceae) based on paired-sites models and the secondary structures of ITS and ETS sequences

To establish relationships among the genera Blepharocalyx, Luma, and Myrceugenia, the phylogeny of tribe Myrteae was reconstructed based on secondary structures of the sequences of ITS (ITS1-5.8S-ITS2) and ETS regions from 93 taxa belonging to 29 genera. The DNA sequences were aligned according to their secondary structure and analysed with Bayesian inference (BI) using base substitution models for RNA, which can differentiate between the regions consisting of base pairs (helices) and unpaired bases (loops). The best-fit models were RNA7C for ITS and RNA7B for ETS, which differ in how they predict double substitutions and symmetry of the frequency of base pairs. The phylogenetic hypothesis was compared with results obtained using classical 4-state models, the results of maximum parsimony, and using sequence alignment without considering the secondary structure. Analyses show low support along the backbone of the consensus trees, those using substitution models of paired sites showing more highly supported derived clades than substitution models of independent sites. All tests were consistent and differed in the positioning of certain groups, but they also showed that Blepharocalyx, Luma, and Myrceugenia arise from three independent lineages that are not closely related. The substitution model for RNA shows that Luma is a monophyletic group and sister to the other genera of Myrteae, except for Myrtus communis. Myrceugenia appears as a monophyletic group, except for M. fernandeziana, which is related to Blepharocalyx, a genus that appears to be biphyletic.     

Accurate and efficient reconstruction of deep phylogenies from structured RNAs

Ribosomal RNA (rRNA) genes are probably the most frequently used data source in phylogenetic reconstruction. Individual columns of rRNA alignments are not independent as a consequence of their highly conserved secondary structures. Unless explicitly taken into account, these correlation can distort the phylogenetic signal and/or lead to gross overestimates of tree stability. Maximum likelihood and Bayesian approaches are of course amenable to using RNA-specific substitution models that treat conserved base pairs appropriately, but require accurate secondary structure models as input. So far, however, no accurate and easy-to-use tool has been available for computing structure-aware alignments and consensus structures that can deal with the large rRNAs. The RNAsalsa approach is designed to fill this gap. Capitalizing on the improved accuracy of pairwise consensus structures and informed by a priori knowledge of group-specific structural constraints, the tool provides both alignments and consensus structures that are of sufficient accuracy for routine phylogenetic analysis based on RNA-specific substitution models. The power of the approach is demonstrated using two rRNA data sets: a mitochondrial rRNA set of 26 Mammalia, and a collection of 28S nuclear rRNAs representative of the five major echinoderm groups.     

ITS secondary structure derived from comparative analysis: implications for sequence alignment and phylogeny of the Asteraceae

An RNA secondary structure model is presented for the nuclear ribosomal internal transcribed spacers (ITS) based on comparative analysis of 340 sequences from the angiosperm family Asteraceae. The model based on covariation analysis agrees with structural features proposed in previous studies using mainly thermodynamic criteria and provides evidence for additional structural motifs within ITS1 and ITS2. The minimum structure model suggests that at least 20% of ITS1 and 38% of ITS2 nucleotide positions are involved in base pairing to form helices. The sequence alignment enabled by conserved structural features provides a framework for broadscale molecular evolutionary studies and the first family-level phylogeny of the Asteraceae based on nuclear DNA data. The phylogeny based on ITS sequence data is very well resolved and shows considerable congruence with relationships among major lineages of the family suggested by chloroplast DNA studies, including a monophyletic subfamily Asteroideae and a paraphyletic subfamily Cichorioideae. Combined analyses of ndhF and ITS sequences provide additional resolution and support for relationships in the family.     

Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices.

We have recently shown that isoalloxazine derivatives are able to photocleave RNA specifically at G.U base pairs embedded within a helical stack. The reaction involves the selective molecular recognition of G.U base pairs by the isoalloxazine ring and the removal of one nucleoside downstream of the uracil residue. Divalent metal ions are absolutely required for cleavage. Here we extend our studies to complex natural RNA molecules with known secondary and tertiary structures, such as tRNAs and a group I intron (td). G.U pairs were cleaved in accordance with the phylogenetically and experimentally derived secondary and tertiary structures. Tandem G.U pairs or certain G.U pairs located at a helix extremity were not affected. These new cleavage data, together with the RNA crystal structure, allowed us to perform molecular dynamics simulations to provide a structural basis for the observed specificity. We present a stable structural model for the ternary complex of the G. U-containing helical stack, the isoalloxazine molecule and a metal ion. This model provides significant new insight into several aspects of the cleavage phenomenon, mechanism and specificity for G. U pairs. Our study shows that in large natural RNAs a secondary structure motif made of an unusual base pair can be recognized and cleaved with high specificity by a low molecular weight molecule. This photocleavage reaction thus opens up the possibility of probing the accessibility of G.U base pairs, which are endowed with specific structural and functional roles in numerous structured and catalytic RNAs and interactions of RNA with proteins, in folded RNAs.     

Evolutionarily conserved structural features in the ITS2 of mammalian pre-rRNAs and potential interactions with the snoRNA U8 detected by comparative analysis of new mouse sequences.

Mechanisms of ITS2 excision from pre-rRNA remain largely elusive. In mammals, at least two endonucleolytic cleavages are involved, which result in the transient accumulation of precursors to 5.8S rRNA termed 8S and 12S RNAs. We have sequenced ITS2 in four new species of the Mus genus and investigated its secondary structure using thermodynamic prediction and comparative approach. Phylogenetic evidence supports an ITS2 folding organized in four domains of secondary structure extending from a preserved structural core. This folding is also largely conserved for the previously available mammalian ITS2 sequences, rat and human, despite their extensive sequence divergence relative to the Mus species. Conserved structural features include the structural core, containing the 3' end of 8S pre-rRNA within a single-stranded sequence, and a stem containing the 3' end of the 12S pre-rRNA species. A putative, phylogenetically preserved pseudoknot has been detected 1 nt downstream from the 12S 3' end. Two long complementarities have also been identified, in sequences conserved among vertebrates, between the pre-rRNA 32S and the snoRNA (small nucleolar RNA) U8 which is required for the excision of Xenopus ITS2. The first complementarity involves the 5.8S-ITS2 junction and 13 nt at the 5' end of U8, whereas the other one occurs between a mature 28S rRNA segment known to be required for ITS2 excision and positions 15-25 of snoRNA U8. These two potential interactions, in combination with ITS2 folding, could organize a functional pocket containing three cleavage sites and key elements for pre-rRNA processing, suggesting a chaperone role for the snoRNA U8.     

Evolution of rDNA ITS1 and ITS2 sequences and RNA secondary structures within members of the fungal genera Grosmannia and Leptographium

The two internal transcribed spacers (ITS) of the nuclear ribosomal (r) DNA tandem repeat were examined in ophiostomatoid fungi belonging to the genera Grosmannia and Leptographium and closely-related taxa. Although the DNA sequence of the ITS region evolves rapidly, core features of the RNA secondary structure of the ITS1 and ITS2 segments are conserved. The results demonstrate that structural conservation of GC-rich helical regions is facilitated primarily through compensatory base changes (CBCs), hemi-CBCs, and compensating insertions/deletions (indels), although slippage of the RNA strand is potentially an additional mechanism for maintaining basepairing interactions. The major conclusion of the structural analysis of both ITS segments is that two factors appear to be involved in limiting the type of changes observed: a high GC bias for both ITS1 and ITS2 and structural constraints at the RNA level.     

Demonstration of the GC-rich common arm in yeast ribosomal 5.8S RNA via 500-MHz proton nuclear magnetic resonance and Overhauser enhancements

In this paper we report the first 1H NMR study of the base-paired secondary structure of yeast 5.8S RNA. On the basis of a combination of homonuclear Overhauser enhancements and temperature dependence of the proton 500-MHz NMR spectrum, we are able to identify and assign eight of the nine base pairs in the most thermally stable helical arm: G116.C137-C117.G136-C118.G135- C119.G134-C120.G133-U121.G132- U122.A131-G123.C130. This arm contains an unusually temperature-stable (to 71 degrees C) segment of four consecutive G.C base pairs. This work constitutes the most direct evidence to date for the existence and base-pair sequence of the GC-rich helix, which is common to most currently popular secondary structural models for eukaryotic 5.8S ribosomal RNA.     

Molecular authentication of three Italian melon accessions by ARMS-PCR and ITS1 (internal transcribed spacer 1) secondary structure prediction

Genetic assessment was carried out on three Italian melon accessions by sequence and structural analysis of the internal transcribed spacer 1 (ITS1) from three populations belonging to two Cucumis melo L. varieties (madras and tendral). Alignment of the 18S-5.8S-26S sequences from three melon accessions showed that there were three single-nucleotide polymorphisms (SNPs) and one short insertion-deletion (indel) at the 5'end ITS1. An amplification refractory mutation system (ARMS)-PCR-based analysis was successfully applied to the SNP markers of the ITS1 sequences for the fingerprinting analysis of three melon populations. Secondary structure models for each ITS1 were derived. The prediction of ITS1 RNA secondary structure from each accession was improved by detecting key functional elements shared by all sequences in the alignments. Our results demonstrated that the ITS1secondary structure models can be used to improve the preliminary genetic assessment of the three melon accessions, suggesting a new tool in plant fingerprinting analysis.     

Identification and assignment of base pairs in three helical stems of Torulopsis utilis ribosomal 5S RNA and its RNase T1 and RNase T2 cleaved fragments via 500-MHz proton homonuclear overhauser enhancements

Imino proton resonances in the downfield region (10-14 ppm) of the 500-MHz 1H NMR spectrum of Torulopsis utilis 5S RNA are identified (A X U, G X C, or G X U) and assigned to base pairs in helices I, IV, and V via analysis of homonuclear Overhauser enhancements (NOE) from intact T. utilis 5S RNA, its RNase T1 and RNase T2 digested fragments, and a second yeast (Saccharomyces cerevisiae) 5S RNA whose nucleotide sequence differs at only six residues from that of T. utilis 5S RNA. The near-identical chemical shifts and NOE behavior of most of the common peaks from these four RNAs strongly suggest that helices I, IV, and V retain the same conformation after RNase digestion and that both T. utilis and S. cerevisiae 5S RNAs share a common secondary and tertiary structure. Of the four G X U base pairs identified in the intact 5S RNA, two are assigned to the terminal stem (helix I) and the other two to helices IV and V. Seven of the nine base pairs of the terminal stem have been assigned. Our experimental demonstration of a G X U base pair in helix V supports the 5S RNA secondary structural model of Luehrsen and Fox [Luehrsen, K. R., & Fox, G.E. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 2150-2154]. Finally, the base-pair proton peak assigned to the terminal G X U in helix V of the RNase T2 cleaved fragment is shifted downfield from that in the intact 5S RNA, suggesting that helices I and V may be coaxial in intact T. utilis 5S RNA.     

Secondary structure of MS2 phage RNA and bias in code word usage.

Based on the secondary structural model of MS2 RNA, it is shown that, in base-pairing regions of the RNA, there is a bias in the use of synonymous codons which favours C and/or G over U and/or A in the third codon positions, and that in non-pairing regions, there is an opposite bias which favours U and/or A over C and/or G. This nature is interpreted as a result of selective constraint which stabilises the secondary structure of the single-stranded RNA genome of the MS2 phage.     

长苞铁杉nrDNA内转录间隔区及5.8S的二级结构分析研究

核糖体RNA及相邻区域的二级结构研究,作为一个重要的工具,已在一些分类等级上被应用于系统发育的分析。以长苞铁杉为实验材料,通过克隆、测序,利用最小自由能原理预测nrDNA内转录间隔区及5.8S转录本的二级结构,分析它们的结构特点,探讨假基因化拷贝与功能拷贝结构上的差异。分析结果表明:(1)ITS1区的二级结构主要由几个延展的发夹结构组成,配对的亚重复单位在松科植物特有的保守序列处有部分重叠,未配对的亚重复单位通常能自身折叠,保守序列的部分碱基出现在发夹结构的环中;(2)假基因化拷贝二级结构的自由能比正常拷贝高;(3)与正常拷贝的二级结构相比,假基因化拷贝在进化速率很低的5.8S功能区发生较大的变异,且在5.8 S末端没有和26 S连接配对。     

Evaluation of the internal transcribed spacer 2 (ITS2) as a molecular marker for phylogenetic inference using sequence and secondary structure information in blow flies (Diptera: Calliphoridae)

The internal transcribed spacer 2 (ITS2) is a small non-coding region located inside the nuclear ribosomal DNA cluster. ITS2 sequence variability is thought to be appropriate to differentiate species and for phylogenetic reconstructions analyses, which can be further improved if structural information is considered. We evaluated the potential of ITS2 as a molecular marker for phylogenetic inference in Calliphoridae (Diptera: Brachycera) using a broad range of inference methods and different substitution models, accounting or not for structural information. Sequence analyses revealed a hierarchically organized pattern of sequence variation and a small level of nucleotide substitution saturation. Intragenomic variation due to small sequence repeats was found mainly in the most variable domain (IV), but it has no significant impact on the phylogenetic signal at the species level. Inferred secondary structures revealed that GC pairs are more frequently found flanking bulges and loops regions in more conserved domains, thus ensuring structure stability. In the phylogenetic analyses, the use of substitution models accounting for structural information significantly improves phylogenetic inference in both neighbour-joining and Bayesian analyses, although the former provides limited resolution for dealing with highly divergent sequences. For Bayesian analyses, a significant improvement in likelihood was observed when considering structure information, although with small changes in topology and overall support, probably reflecting better evolutionary rates estimates. Based on these findings, ITS2 is a suitable molecular marker for phylogenetic analyses in Calliphoridae, at both species and generic level.