Compensatory Base Pair Changes Research Articles

Does rapid sequence divergence preclude RNA structure conservation in vertebrates?

Accelerated evolution of any portion of the genome is of significant interest, potentially signaling positive selection of phenotypic traits and adaptation. Accelerated evolution remains understudied for structured RNAs, despite the fact that an RNA’s structure is often key to its function. RNA structures are typically characterized by compensatory (structure-preserving) basepair changes that are unexpected given the underlying sequence variation, i.e., they have evolved through negative selection on structure. We address the question of how fast the primary sequence of an RNA can change through evolution while conserving its structure. Specifically, we consider predicted and known structures in vertebrate genomes. After careful control of false discovery rates, we obtain 13 de novo structures (and three known Rfam structures) that we predict to have rapidly evolving sequences—defined as structures where the primary sequences of human and mouse have diverged at least twice as fast (1.5 times for Rfam) as nearby neutrally evolving sequences. Two of the three known structures function in translation inhibition related to infection and immune response. We conclude that rapid sequence divergence does not preclude RNA structure conservation in vertebrates, although these events are relatively rare.

Structural Alignment and Covariation Analysis of RNA Sequences.

RNA molecules adopt defined structural conformations that are essential to exert their function. During the course of evolution, the structure of a given RNA can be maintained via compensatory base-pair changes that occur among covarying nucleotides in paired regions. Therefore, for comparative, structural, and evolutionary studies of RNA molecules, numerous computational tools have been developed to incorporate structural information into sequence alignments and a number of tools have been developed to study covariation. The bioinformatic protocol presented here explains how to use some of these tools to generate a secondary-structure-aware multiple alignment of RNA sequences and to annotate the alignment to examine the conservation and covariation of structural elements among the sequences.

RNAscClust: clustering RNA sequences using structure conservation and graph based motifs.

MotivationClustering RNA sequences with common secondary structure is an essential step towards studying RNA function. Whereas structural RNA alignment strategies typically identify common structure for orthologous structured RNAs, clustering seeks to group paralogous RNAs based on structural similarities. However, existing approaches for clustering paralogous RNAs, do not take the compensatory base pair changes obtained from structure conservation in orthologous sequences into account.ResultsHere, we present RNAscClust, the implementation of a new algorithm to cluster a set of structured RNAs taking their respective structural conservation into account. For a set of multiple structural alignments of RNA sequences, each containing a paralog sequence included in a structural alignment of its orthologs, RNAscClust computes minimum free-energy structures for each sequence using conserved base pairs as prior information for the folding. The paralogs are then clustered using a graph kernel-based strategy, which identifies common structural features. We show that the clustering accuracy clearly benefits from an increasing degree of compensatory base pair changes in the alignments.Availability and Implementation RNAscClust is available at http://www.bioinf.uni-freiburg.de/Software/RNAscClust.Supplementary information Supplementary data are available at Bioinformatics online.

Cross talk between ABC transporter mRNAs via a target mRNA-derived sponge of the GcvB small RNA.

There is an expanding list of examples by which one mRNA can posttranscriptionally influence the expression of others. This can involve RNA sponges that sequester regulatory RNAs of mRNAs in the same regulon, but the underlying molecular mechanism of such mRNA cross talk remains little understood. Here, we report sponge-mediated mRNA cross talk in the posttranscriptional network of GcvB, a conserved Hfq-dependent small RNA with one of the largest regulons known in bacteria. We show that mRNA decay from the gltIJKL locus encoding an amino acid ABC transporter generates a stable fragment (SroC) that base-pairs with GcvB. This interaction triggers the degradation of GcvB by RNase E, alleviating the GcvB-mediated mRNA repression of other amino acid-related transport and metabolic genes. Intriguingly, since the gltIJKL mRNA itself is a target of GcvB, the SroC sponge seems to enable both an internal feed-forward loop to activate its parental mRNA in cis and activation of many trans-encoded mRNAs in the same pathway. Disabling this mRNA cross talk affects bacterial growth when peptides are the sole carbon and nitrogen sources.

Genetic variation in the mitochondrial 16S ribosomal RNA gene of Ixodes scapularis (Acari: Ixodidae).

BackgroundIxodes scapularis is a vector of several human pathogens in the United States, and there is geographical variation in the relative number of persons infected with these pathogens. Geographically isolated populations of I. scapularis have established or are in the process of establishing in southern Canada. Knowledge of the genetic variation within and among these populations may provide insight into their geographical origins in the United States and the potential risk of exposure of Canadians to the different pathogens carried by I. scapularis.MethodsPart of the mitochondrial (mt) 16S ribosomal (r) RNA gene was amplified by PCR from 582 ticks collected from southern Canada, and Minnesota and Rhode Island in the United States. Sequence variation was examined in relation to the predicted secondary structure of the gene. Genetic diversity among populations was also determined.ResultsDNA sequence analyses revealed 52 haplotypes. Most mutational alterations in DNA sequence occurred at unpaired sites or represented partial compensatory base pair changes that maintained the stability of the secondary structure. Significant genetic variation was detected within and among populations in different geographical regions. A greater proportion of the haplotypes of I. scapularis from the Canadian Prairie Provinces were found in the Midwest of the United States than in other regions, whereas more of the haplotypes of I. scapularis from the Canadian Central and Atlantic Provinces occurred in the Northeast of the United States. Nonetheless, 58% of I. scapularis were of a haplotype that occurs in the Midwest and Northeast of the United States; thus, their geographical origins could not be determined.ConclusionsThere is considerable genetic variation in the mt 16S rRNA gene of I. scapularis. There is some evidence to support the hypothesis that some lineages of I. scapularis in the Atlantic and Central Provinces of Canada may be derived from colonizing individuals originating in the Northeast of the United States, whereas those in the Prairie Provinces may be derived from individuals originating in the Midwest of the United States. However, additional genetic markers are needed to test hypotheses concerning the geographical origins of I. scapularis in Canada.Electronic supplementary materialThe online version of this article (doi:10.1186/s13071-014-0530-6) contains supplementary material, which is available to authorized users.

Mitochondrial and nuclear phylogenetic analysis with Sanger and next-generation sequencing shows that, in Área de Conservación Guanacaste, northwestern Costa Rica, the skipper butterfly named Urbanus belli (family Hesperiidae) comprises three morphologically cryptic species.

BackgroundSkipper butterflies (Hesperiidae) are a relatively well-studied family of Lepidoptera. However, a combination of DNA barcodes, morphology, and natural history data has revealed several cryptic species complexes within them. Here, we investigate three DNA barcode lineages of what has been identified as Urbanus belli (Hesperiidae, Eudaminae) in Área de Conservación Guanacaste (ACG), northwestern Costa Rica.ResultsAlthough no morphological traits appear to distinguish among the three, congruent nuclear and mitochondrial lineage patterns show that “Urbanus belli” in ACG is a complex of three sympatric species. A single strain of Wolbachia present in two of the three cryptic species indicates that Urbanus segnestami Burns (formerly Urbanus belliDHJ01), Urbanus bernikerni Burns (formerly Urbanus belliDHJ02), and Urbanus ehakernae Burns (formerly Urbanus belliDHJ03) may be biologically separated by Wolbachia, as well as by their genetics. Use of parallel sequencing through 454-pyrosequencing improved the utility of ITS2 as a phylogenetic marker and permitted examination of the intra- and interlineage relationships of ITS2 variants within the species complex. Interlineage, intralineage and intragenomic compensatory base pair changes were discovered in the secondary structure of ITS2.ConclusionThese findings corroborate the existence of three cryptic species. Our confirmation of a novel cryptic species complex, initially suggested by DNA barcode lineages, argues for using a multi-marker approach coupled with next-generation sequencing for exploration of other suspected species complexes.

Comparison of the DNA sequences and secondary structure of the mitochondrial 16S rRNA gene of Ixodes kingi, Ixodes sculptus and Ixodes angustus

The DNA sequences and predicted secondary structure of domains IV and V of the mitochondrial (mt) 16S rRNA gene were compared for three species of Ixodes: Ixodes kingi, Ixodes sculptus and Ixodes angustus. Each species had a unique set of DNA sequences for the 16S gene. Many of the differences in DNA sequence within and among species occurred in a “hypervariable” region of domain V, and either represented partial or full compensatory base pair changes that maintained the helices within the secondary structure, or nucleotide alterations at unpaired positions that had no effect on the secondary structure. The results of the phylogenetic analyses revealed that I. kingi, I. sculptus and I. angustus were placed in a clade with some other species of the subgenera Pholeoixodes and Ixodiopsis. In addition, individuals of I. sculptus from Saskatchewan (Canada) and Colorado (USA) did not form a monophyletic clade, suggesting the possible existence of cryptic species.

Molecular subdivision of the marine diatom Thalassiosira rotulain relation to geographic distribution, genome size, and physiology

BackgroundMarine phytoplankton drift passively with currents, have high dispersal potentials and can be comprised of morphologically cryptic species. To examine molecular subdivision in the marine diatom Thalassiosira rotula, variations in rDNA sequence, genome size, and growth rate were examined among isolates collected from the Atlantic and Pacific Ocean basins. Analyses of rDNA included T. gravida because morphological studies have argued that T. rotula and T. gravida are conspecific.ResultsCulture collection isolates of T. gravida and T. rotula diverged by 7.0 ± 0.3% at the ITS1 and by 0.8 ± 0.03% at the 28S. Within T. rotula, field and culture collection isolates were subdivided into three lineages that diverged by 0.6 ± 0.3% at the ITS1 and 0% at the 28S. The predicted ITS1 secondary structure revealed no compensatory base pair changes among lineages. Differences in genome size were observed among isolates, but were not correlated with ITS1 lineages. Maximum acclimated growth rates of isolates revealed genotype by environment effects, but these were also not correlated with ITS1 lineages. In contrast, intra-individual variation in the multi-copy ITS1 revealed no evidence of recombination amongst lineages, and molecular clock estimates indicated that lineages diverged 0.68 Mya. The three lineages exhibited different geographic distributions and, with one exception, each field sample was dominated by a single lineage.ConclusionsThe degree of inter- and intra-specific divergence between T. gravida and T. rotula suggests they should continue to be treated as separate species. The phylogenetic distinction of the three closely-related T. rotula lineages was unclear. On the one hand, the lineages showed no physiological differences, no consistent genome size differences and no significant changes in the ITS1 secondary structure, suggesting there are no barriers to interbreeding among lineages. In contrast, analysis of intra-individual variation in the multicopy ITS1 as well as molecular clock estimates of divergence suggest these lineages have not interbred for significant periods of time. Given the current data, these lineages should be considered a single species. Furthermore, these T. rotula lineages may be ecologically relevant, given their differential abundance over large spatial scales.

Metapopulation Structure in the Planktonic Diatom Ditylum brightwellii (Bacillariophyceae)

Approximately 200,000 diatom species are thought to exist and yet the underlying processes of speciation in diatoms are unknown. Because genetic subdivision within species can reveal potential speciation mechanisms, we examined genetic differentiation and patterns of gene flow among four populations of the diatom Ditylum brightwellii. Single-cell isolates were examined at two microsatellite markers and two rDNA loci (18S and internal transcribed spacer region I (ITSI)). Among isolates, rDNA sequences varied by 0.08+/-0.04% (18S) and 0.7+/-0.3% (ITSI) and there were no compensatory base pair changes in the predicted ITSI secondary structure, all suggesting that a single species was represented. Two numerically dominant ITSI sequence types were detected and their distribution among isolates from genetically distinct populations was significantly different. Two populations shared ITSI sequence type 1 and two shared ITSI sequence type 2, indicating differences in relatedness among populations. The signature of unequal gene flow among populations suggested that D. brightwellii exhibited a metapopulation structure: the species was subdivided into populations of populations. The identification of metapopulations suggests a possible mechanism of speciation through reduced levels of gene flow, providing newly evolved taxa with a large repository of genetic and physiological diversity and perhaps significant adaptive potential.

Nucleotide alterations in the D3 domain of the large subunit of ribosomal DNA among 21 species of equine strongyle

The expansion segments or divergent (D) domains in the large subunit (LSU) of the ribosomal DNA have been suggested as genetic markers for taxonomic and/or phylogenetic studies of parasites. In the present study, we assessed the degree of sequence variation in the D3 domain and flanking core regions of the LSU for 21 species of equine strongyles (Strongylida: Strongylidae) and determined which positions in the secondary structure of the LSU were associated with the nucleotide alterations. No intraspecific sequence variation was detected in 17 species, for which multiple individual worms were available. Mutations in sequence among species were detected at 19 nucleotide positions, most of which were located in the D3 domain. Fifteen alterations were transitions, three were transversions and one represented a site of multiple mutations. In relation to the secondary structure element of D3, 26% of these mutations were located in unpaired regions (i.e., end of loops, or in bulges of helices) and thus did not appear to alter the pairing arrangement in the helices of the secondary structure. Many of the other mutations represented partial or complete compensatory base pair changes. The magnitude of interspecific nucleotide variation in the D3 domain (0-4%) was considerably less than that recorded for some other nematode groups (enoplids and thelastomatoids), indicating that this region alone is of limited value for taxonomic and phylogenetic studies for strongyles of equids but is interesting in relation to the evolution of ribosomal DNA.

Molecular evolution of the trnTUGU-trnFGAA region in Bryophytes.

Structure, variability, and molecular evolution of the trnT-F region in the Bryophyta (mosses and liverworts) is analyzed based on about 200 sequences of the trnT-L spacer and trnL 5' exon, 1000 sequences of the trnL intron, and 800 sequences of the trnL 3' exon and trnL-F spacer, including comparisons of lengths, GC contents, sequence similarities, and functional elements. Mutations occurring in the trnL 5' and 3' exons, including compensatory base pair changes, and a transition in the trnL anticodon in Takakia lepidozioides, are discussed. All three non-coding regions display a mosaic structure of highly variable elements (V1 - V3 in the trnT-L spacer, V4/V5 corresponding to stem-loop regions P6/P8 in the trnL intron, and V6/V7 in the trnL-F spacer) and more conserved elements. In the trnL intron this structure is a consequence of the defined secondary structure necessary for correct splicing, whereas in both spacers conserved regions are restricted to promoter elements. At least the highly variable regions in the trnT-L spacer and stem-loop region P8 of the trnL intron seem to evolve independently in the major bryophyte lineages and are therefore not suitable for high taxonomic level phylogenetic reconstructions. In mosses, a trend of length reduction towards the more derived lineages is observed in all three non-coding regions. GC contents are mostly linked to sequence variability, with the conserved regions being more GC rich and the more variable AT rich. The lowest GC values (< 10 %) are found in the trnT-L spacer of mosses. In addition to two putative sigma (70)-type promoters in the trnT-L spacer, a third putative promoter is present in the trnL-F spacer, although trnL and trnF are assumed to be co-transcribed. Consensus sequences are provided for the -35 and -10 sequences of the major bryophyte lineages. The third promoter is part of a hairpin secondary structure, whose loop region is highly homoplastic in mosses due to an inversion occurring independently in non-related taxa, even at the intraspecific level.

Molecular circumscription of the hornworts (Anthocerotophyta) based on the chloroplast DNA trnL-trnF region.

In phylogenetic trees generated from partial trnL(UAA) intron sequences, the hornworts (represented by nine species from the genera Anthoceros, Dendroceros, Megaceros, Notothylas and Phaeoceros) are resolved as a monophyletic group and are separated from the clades of mosses, liverworts and tracheophytes. A secondary structure of the trnL(UAA) intron of Anthoceros agrestis is presented, displaying the arrangement of the stem-loop regions P1-P9. Compensatory base-pair changes (coevolutionary sites) are detected in regions P4/5 and P9 within the hornwort sequences. The original homology of the most variable region, P8, cannot be detected anymore due to the extremely fast divergent evolution of this segment in the major land plant groups. Similarly, a high sequence divergence occurs in the trnL-trnF intergenic spacer. Apart from synapomorphic substitutions in the trnL(UAA) intron, the hornworts are characterised by a large P6 region consisting of many repetitive elements. The molecular data therefore support the hornworts as representing an independent land plant lineage (Anthocerotophyta). Although relationships between hornworts and the other land plant groups remain unresolved in the trnL(UAA) intron trees, it is rather unlikely that bryophytes are monophyletic in their traditional circumscription, i.e. comprising hornworts, mosses and liverworts.

First complete large subunit ribosomal RNA sequence and secondary structure for a parasitic nematode: phylogenetic and diagnostic implications

The complete sequence and secondary structure of the large subunit of nuclear ribosomal RNA(LSUrRNA) were determined for the parasitic nematode Labiostrongylus bipapillosus (order Strongylida). Its LSU rRNA sequence was shorter (by 18 bp) than that of the free-living nematode, Caenorhabditis elegans (order Rhabditida), the only other species within the Nematoda for which a complete LSU rRNA sequence has been determined. Interspecific differences in sequence were greater in the 12 D domains compared with the core segments, with the secondary structure being maintained by partial or complete compensatory base pair changes. The magnitude of interspecific sequence difference in each D domain (except for D6 and D12) was similar, suggesting that several domains contain informative genetic markers for phylogenetic studies of the phylum Nematoda at different taxonomic levels. The LSU rRNA may also provide species-specific markers for the identification of some bursate nematodes of veterinary and medical importance.

In vitro evolution used to define a protein recognition site within a large RNA domain

A minimum of 460 nucleotides of 16S ribosomal RNA are needed to fold the target site for E. coli ribosomal protein S4, although a much smaller region within this large domain is protected from chemical reagents by the protein. Starting with a 531-nucleotide rRNA fragment, cycles of mutagenesis, selection with S4, and amplification (‘in vitro evolution’) were used to obtain a pool of 30 RNA sequences selected for S4 recognition but ∼30% different from wild type. Numerous compensatory base pair changes have largely preserved the same secondary structure among these RNAs as found in wild-type sequences. A 20-base deletion and a single nucleotide insertion are among several unusual features found in most of the selected sequences and also prevalent among other prokaryotic rRNAs. Most of the compensatory base changes and selected features are located outside of the region protected by S4 from chemical reagents. It was unexpected that S4 would select for RNA structures throughout such a large domain; the selected features are probably contributing indirectly to S4 recognition by promoting correct tertiary folding of the region actually contacted by S4. The role of S4 may be to stabilize this domain (nearly one-third of the 16S rRNA) in its proper conformation for ribosome function.

Structural analysis of the peptidyl transferase region in ribosomal RNA of the eukaryote Xenopus laevis

Accessible single-strand bases in Xenopus laevis 28 S ribosomal RNA (rRNA) Domain V, the peptidyl transferase region, were determined by chemical modification with dimethylsulfate, 1-cyclohexyl-3-(2-morpholinoethyl-carbodiimide metho- p-toluene sulfonate and kethoxal, followed by primer extension. The relative accessibilities of three rRNA substrates were compared: deproteinized 28 S rRNA under non-denaturing conditions (free 28 S rRNA), 60 S subunits and 80 S ribosomes. Overall, our experimental results support the theoretical secondary structure model of Domain V derived by comparative sequence analysis and compensatory base-pair changes, and support some theoretical tertiary interactions previously suggested by covariation. The 60 S subunits and 80 S ribosomes generally show increasing resistance to chemical modification. Bases which are sensitive in free 28 S rRNA but protected in 60 S subunits may be sites for ribosomal protein binding or induced structural rearrangements. Another class of nucleotides is distinguished by its sensitivity in 60 S subunits but protection in 80 S ribosomes; these nucleotides may be involved in subunit-subunit interactions or located at the interface of the ribosome. We found a third class of bases, which is protected in free 28 S rRNA but sensitive in 60 S subunits and/or 80 S ribosomes, suggesting that structural changes occur in Domain V as a result of subunit assembly and ribosome formation. One such region is uniquely hypersensitive in eukaryotic ribosomes but is absent in Escherichia coli ribosomes. Sites that we determined to be accessible on empty 80 S ribosomes could serve as recognition sites for translation components.

Unusual mRNA pseudoknot structure is recognized by a protein translational repressor

Translation of ribosomal proteins in the α operon of E. coli is repressed by one of the encoded proteins, S4; it specifically recognizes an RNA fragment containing the translational initiation site for the first gene in the operon. RNA structure mapping experiements have suggested a pseudoknot structure for the S4 binding site: the loop of a hairpin is base paired to sequences downstream of the hairpin. Here, we systematically test this proposed structure by measuring S4 binding to an extensive set of site-directed mutations that create compensatory base pair changes in potential helices. The pseudoknot folding is confirmed, and two additional, unexpected interactions within the pseudoknot are also detected. The overall structure is an unusual “double pseudoknot” linking a hairpin up-stream of the ribosome binding site with sequences 2–10 codons downstream of the initiation codon. Stabilization of this structure by S4 could account for translational repression.