RESUMEN

The beetle suborder Adephaga is traditionally divided into two sections on the basis of habitat, terrestrial Geadephaga and aquatic Hydradephaga. Monophyly of both groups is uncertain, and the relationship of the two groups has implications for inferring habitat transitions within Adephaga. Here we examine phylogenetic relationships of these groups using evidence provided by DNA sequences from all four suborders of beetles, including 60 species of Adephaga, four Archostemata, three Myxophaga, and ten Polyphaga. We studied 18S ribosomal DNA and 28S ribosomal DNA, aligned with consideration of secondary structure, as well as the nuclear protein-coding gene wingless. Independent and combined Bayesian, likelihood, and parsimony analyses of all three genes supported placement of Trachypachidae in a monophyletic Geadephaga, although for analyses of 28S rDNA and some parsimony analyses only if Coleoptera is constrained to be monophyletic. Most analyses showed limited support for the monophyly of Hydradephaga. Outside of Adephaga, there is support from the ribosomal genes for a sister group relationship between Adephaga and Polyphaga. Within the small number of sampled Polyphaga, analyses of 18S rDNA, wingless, and the combined matrix supports monophyly of Polyphaga exclusive of Scirtoidea. Unconstrained analyses of the evolution of habitat suggest that Adephaga was ancestrally aquatic with one transition to terrestrial. However, in analyses constrained to disallow changes from aquatic to terrestrial habitat, the phylogenies imply two origins of aquatic habit within Adephaga.

RESUMEN

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

Asunto(s)

RESUMEN

We report the entire sequence (2864 nts) and secondary structure of the nuclear small subunit ribosomal RNA (SSU rRNA) gene (18S) from the twisted-wing parasite Caenocholax fenyesi texensis Kathirithamby & Johnston (Strepsiptera: Myrmecolacidae). The majority of the base pairings in this structural model map on to the SSU rRNA secondary and tertiary helices that were previously predicted with comparative analysis. These regions of the core rRNA were unambiguously aligned across all Arthropoda. In contrast, many of the variable regions, as previously characterized in other insect taxa, had very large insertions in C. f. texensis. The helical base pairs in these regions were predicted with a comparative analysis of a multiple sequence alignment (that contains C. f. texensis and 174 published arthropod 18S rRNA sequences, including eleven strepsipterans) and thermodynamic-based algorithms. Analysis of our structural alignment revealed four unusual insertions in the core rRNA structure that are unique to animal 18S rRNA and in general agreement with previously proposed insertion sites for strepsipterans. One curious result is the presence of a large insertion within a hairpin loop of a highly conserved pseudoknot helix in variable region 4. Despite the extraordinary variability in sequence length and composition, this insertion contains the conserved sequences 5'-AUUGGCUUAAA-3' and 5'-GAC-3' that immediately flank a putative helix at the 5'- and 3'-ends, respectively. The longer sequence has the potential to form a nine base pair helix with a sequence in the variable region 2, consistent with a recent study proposing this tertiary interaction. Our analysis of a larger set of arthropod 18S rRNA sequences has revealed possible errors in some of the previously published strepsipteran 18S rRNA sequences. Thus we find no support for the previously recovered heterogeneity in the 18S molecules of strepsipterans. Our findings lend insight to the evolution of RNA structure and function and the impact large insertions pose on genome size. We also provide a novel alignment template that will improve the phylogenetic placement of the Strepsiptera among other insect taxa.

Asunto(s)

RESUMEN

RNA molecules fold into characteristic secondary and tertiary structures that account for their diverse functional activities. Many of these RNA structures are assembled from a collection of RNA structural motifs. These basic building blocks are used repeatedly, and in various combinations, to form different RNA types and define their unique structural and functional properties. Identification of recurring RNA structural motifs will therefore enhance our understanding of RNA structure and help associate elements of RNA structure with functional and regulatory elements. Our goal was to develop a computer program that can describe an RNA structural element of any complexity and then search any nucleotide sequence database, including the complete prokaryotic and eukaryotic genomes, for these structural elements. Here we describe in detail a new computational motif search algorithm, RNAMotif, and demonstrate its utility with some motif search examples. RNAMotif differs from other motif search tools in two important aspects: first, the structure definition language is more flexible and can specify any type of base-base interaction; second, RNAMotif provides a user controlled scoring section that can be used to add capabilities that patterns alone cannot provide.

Asunto(s)

RESUMEN

How group I introns originate in nuclear ribosomal (r)RNA genes is an important question in evolutionary biology. Central to this issue is the multitude of group I introns present in evolutionarily distantly related plant, fungal, and protist lineages, together with an understanding of their origin and lateral transfer from one exon to another, between cell organelles, and between cells. These introns vary considerably in primary and secondary structure; and their provenance from a few or perhaps many mobile elements that have spread in rRNAs is unknown. Here we show that a novel lineage of group IC1 introns inserted at position 516 (Escherichia coli gene numbering) in the small subunit rRNA in bangiophyte red algae and a brown alga (Aureoumbra lagunensis) are specifically related, although their host cells are not. These bangiophyte and Aureoumbra introns are the only known cases that have a helical insertion in the P5b helix. The highly conserved primary and secondary structure of the extra P5b helix suggests that it is important, although its specific function is unknown. Our study attempts to understand the origin and movement of these IC1 introns.

Asunto(s)

RESUMEN

Our previous study of the North American biogeography of Bangia revealed the presence of two introns inserted at positions 516 and 1506 in the nuclear-encoded SSU rRNA gene. We subsequently sequenced nuclear SSU rRNA in additional representatives of this genus and the sister genus Porphyra in order to examine the distribution, phylogeny, and structural characteristics of these group I introns. The lengths of these introns varied considerably, ranging from 467 to 997 nt for intron 516 and from 509 to 1,082 nt for intron 1506. The larger introns contained large insertions in the P2 domain of intron 516 and the P1 domain of intron 1506 that correspond to open reading frames (ORFs) with His-Cys box homing endonuclease motifs. These ORFs were found on the complementary strand of the 1506 intron in Porphyra fucicola and P. umbilicalis (HG), unlike the 516 intron in P. abbottae, P. kanakaensis, P. tenera (SK), Bangia fuscopurpurea (Helgoland), and B. fuscopurpurea (MA). Frameshifts were noted in the ORFs of the 516 introns in P. kanakaensis and B. fuscopurpurea (HL), and all ORFs terminated prematurely relative to the amino acid sequence for the homing endonuclease I-Ppo I. This raises the possibility that these sequences are pseudogenes. Phylogenies generated using sequences of both introns and the 18S rRNA gene were congruent, which indicated long-term immobility and vertical inheritance of the introns followed by subsequent loss in more derived lineages. The introns within the florideophyte species Hildenbrandia rubra (position 1506) were included to determine relationships with those in the Bangiales. The two sequences of intron 1506 analyzed in Hildenbrandia were positioned on a well-supported branch associated with members of the Bangiales, indicating possible common ancestry. Structural analysis of the intron sequences revealed a signature structural feature in the P5b domain of intron 516 that is unique to all Bangialean introns in this position and not seen in intron 1506 or other group IC1 introns.

Asunto(s)

RESUMEN

Polytoma obtusum and Polytoma uvella are members of a clade of nonphotosynthetic chlorophyte algae closely related to Chlamydomonas humicola and other photosynthetic members of the Chlamydomonadaceae. Descended from a nonphotosynthetic mutant, these obligate heterotrophs retain a plastid (leucoplast) with a functional protein synthetic system, and a plastid genome (lpDNA) with functional genes encoding proteins required for transcription and translation. Comparative studies of the evolution of genes in chloroplasts and leucoplasts can identify modes of selection acting on the plastid genome. Two plastid genes--rrn16, encoding the plastid small-subunit rRNA, and tufA, encoding elongation factor Tu--retain their functions in protein synthesis after the loss of photosynthesis in two nonphotosynthetic Polytoma clades but show a substantially accelerated rate of base substitution in the P. uvella clade. The accelerated evolution of tufA is due, at least partly, to relaxed codon bias favoring codons that can be read without wobble, mainly in three amino acids. Selection for these codons may be relaxed because leucoplasts are required to synthesize fewer protein molecules per unit time than are chloroplasts (reduced protein synthetic load) and thus require a lower rate of synthesis of elongation factor Tu. Relaxed selection due to a lower protein synthetic load is also a plausible explanation for the accelerated rate of evolution of rrn16, but the available data are insufficient to test the hypothesis for this gene. The tufA and rrn16 genes in Polytoma oviforme, the sole member of a second nonphotosynthetic clade, are also functional but show no sign of relaxed selection.

Asunto(s)

RESUMEN

This study reveals that AA and AG oppositions occur frequently at the ends of helices in RNA crystal and NMR structures in the PDB database and in the 16 S and 23 S rRNA comparative structure models, with the G usually 3' to the helix for the AG oppositions. In addition, these oppositions are frequently base-paired and usually in the sheared conformation, although other conformations are present in NMR and crystal structures. These A:A and A:G base-pairs are present in a variety of structural environments, including GNRA tetraloops, E and E-like loops, interfaced between two helices that are coaxially stacked, tandem G:A base-pairs, U-turns, and adenosine platforms. Finally, given structural studies that reveal conformational rearrangements occurring in regions of the RNA with AA and AG oppositions at the ends of helices, we suggest that these conformationally unique helix extensions might be associated with functionally important structural rearrangements.

Asunto(s)

RESUMEN

Biomass of the protistan parasite QPX (quahaug parasite X) of hard-shell clam Mercenaria mercenaria was enriched from in vitro culture. The nuclear gene encoding the 18S RNA of the small-subunit ribosomal (ssu-rDNA) was recovered using the polymerase chain reaction (PCR) and sequenced. Phylogenetic analysis clearly showed that QPX is a member of phylum Labyrinthulomycota, within which it appears as a specific relative of Thraustochytrium pachydermum. These results confirm the provisional assignment of QPX to the Labyrinthulomycota made previously on the basis of morphological and ultrastructural characters found in some, but not all, geographic isolates.

Asunto(s)

RESUMEN

In 1985 an analysis of the Escherichia coli 16 S rRNA covariation-based structure model revealed a strong bias for unpaired adenosines. The same analysis revealed that the majority of the G, C, and U bases were paired. These biases are (now) consistent with the high percentage of unpaired adenosine nucleotides in several structure motifs. An analysis of a larger set of bacterial comparative 16 S and 23 S rRNA structure models has substantiated this initial finding and revealed new biases in the distribution of adenosine nucleotides in loop regions. The majority of the adenosine nucleotides are unpaired, while the majority of the G, C, and U bases are paired in the covariation-based structure model. The unpaired adenosine nucleotides predominate in the middle and at the 3' end of loops, and are the second most frequent nucleotide type at the 5' end of loops (G is the most common nucleotide). There are additional biases for unpaired adenosine nucleotides at the 3' end of loops and adjacent to a G at the 5' end of the helix. The most prevalent consecutive nucleotides are GG, GA, AG, and AA. A total of 70 % of the GG sequences are within helices, while more than 70 % of the AA sequences are unpaired. Nearly 50 % of the GA sequences are unpaired, and approximately one-third of the AG sequences are within helices while another third are at the 3' loop.5' helix junction. Unpaired positions with an adenosine nucleotide in more than 50 % of the sequences at the 3' end of 16 S and 23 S rRNA loops were identified and arranged into the A-motif categories XAZ, AAZ, XAG, AAG, and AAG:U, where G or Z is paired, G:U is a base-pair, and X is not an A and Z is not a G in more than 50 % of the sequences. These sequence motifs were associated with several structural motifs, such as adenosine platforms, E and E-like loops, A:A and A:G pairings at the end of helices, G:A tandem base-pairs, GNRA tetraloop hairpins, and U-turns.

Asunto(s)

RESUMEN

We used an Escherichia coli genetic assay based on the phage T4 td intron to test the ability of the Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (CYT-18 protein) to suppress mutations that cause structural defects around its binding site in the P4-P6 domain of the group I intron catalytic core. We analyzed all possible combinations of nucleotides at either P4 bp-1 or P6 bp-1, which together form the junction of the P4-P6 stacked helices, and looked for synergistic effects in double mutants. Most mutations at either position inhibit self-splicing, but can be suppressed by CYT-18. CYT-18 can compensate efficiently for mutations that disrupt base-pairing at either P4 bp-1 or P6 bp-1, for mutations at P6 bp-1 that disrupt the base-triple interaction with J3/4-3, and for nucleotide substitutions at either position that are predicted to be suboptimal for base stacking, based on the analysis of DNA four-way junctions. However, CYT-18 has difficulty suppressing combinations of mutations at P4 bp-1 and P6 bp-1 that simultaneously disrupt base-pairing and base stacking. Thermal denaturation and Fe(II)-EDTA analysis showed that mutations at the junction of the P4-P6 stacked helices lead to grossly impaired tertiary-structure formation centered in the P4-P6 domain. CYT-18-suppressible mutants bind the protein with K(d) values up to 79-fold higher than that for the wild-type intron, but in all cases tested, the k(off) value for the complex remains within twofold of the wild-type value, suggesting that the binding site can be formed properly and that the increased K(d) value reflects primarily an increased k(on) value for the binding of CYT-18 to the misfolded intron. Our results indicate that the P4/P6 junction is a linchpin region, where even small nucleotide substitutions grossly disrupt the catalytically-active group I intron tertiary structure, and that CYT-18 binding induces the formation of the correct structure in this region, leading to folding of the group I intron catalytic core.

Asunto(s)

RESUMEN

A previously described unusual form of the protistan parasite Ichthyophonus, differing in morphological and developmental features from I. hoferi sensu Plehn & Mulsow, was recovered from yellowtail flounder Limanda ferruginea Storer from the Brown's Bank area of the Nova Scotia shelf. The nuclear gene encoding the rRNA of the small ribosomal subunit was amplified from this unusual form of Ichthyophonus using the polymerase chain reaction, sequenced and aligned with other eukaryote small subunit (ssu)-rDNAs. Inferred phylogenetic trees clearly show that its ssu-rDNA is distinct from those of 2 isolates of I. hoferi sensu Plehn & Mulsow from different hosts and geographical locations (herring in the North Sea, and yellowtail flounder from the Nova Scotia shelf). We consider the unusual form to be a separate species, I. irregularis. The occurrence of a second, distinct type of Ichthyophonus within a single host species raises the possibility that ichthyophoniasis could be produced by different (although related) pathogens, and in some cases, by concurrent infections of the two.

Asunto(s)

RESUMEN

The U-turn is a well-known RNA motif characterized by a sharp reversal of the RNA backbone following a single-stranded uridine base. In experimentally determined U-turn motifs, the nucleotides 3' to the turn are frequently involved in tertiary interactions, rendering this motif particularly attractive in RNA modeling and functional studies. The U-turn signature is composed of an UNR sequence pattern flanked by a Y:Y, Y:A (Y=pyrimidine) or G:A base juxtaposition. We have identified 33 potential UNR-type U-turns and 25 related GNRA-type U-turns in a large set of aligned 16 S and 23 S rRNA sequences. U-turn candidates occur in hairpin loops (34 times) as well as in internal and multi-stem loops (24 times). These are classified into ten families based on loop type, sequence pattern (UNR or GNRA) and the nature of the closing base juxtaposition. In 13 cases, the bases on the 3' side of the turn, or on the immediate 5' side, are involved in tertiary covariations, making these sites strong candidates for tertiary interactions.

Asunto(s)

RESUMEN

Mollusks are an extraordinarily diverse group of animals with an estimated 200,000 species, second only to the phylum Arthropoda. We conducted a comparative analysis of complete mitochondrial ribosomal large subunit sequences (LSU) of a chiton, two bivalves, six gastropods, and a cephalopod. In addition, we determined secondary structure models for each of them. Comparative analyses of nucleotide variation revealed substantial length variation among the taxa, with stylommatophoran gastropods possessing the shortest lengths. Phylogenetic analyses of the nucleotide sequence data supported the monophyly of Albinaria, Euhadra herklotsi + Cepaea nemoralis, Stylommatophora, Cerithioidea, and when only transversions are included, the Bivalvia. The phylogenetic limits of the mitochondrial LSU rRNA gene within mollusks appear to be up to 400 million years, although this estimate will have to be tested further with additional taxa. Our most novel finding was the discovery of phylogenetic signal in the secondary structure of rRNA of mollusks. The absence of entire stem/loop structures in Domains II, III, and V can be viewed as three shared derived characters uniting the stylommatophoran gastropods. The absence of the aforementioned stem/loop structure explains much of the observed length variation of the mitochondrial LSU rRNA found within mollusks. The distribution of these unique secondary structure characters within mollusks should be examined.

Asunto(s)

RESUMEN

Comparative sequence analysis complements experimental methods for the determination of RNA three-dimensional structure. This approach is based on the concept that different sequences within the same gene family form similar higher-order structures. The large number of rRNA sequences with sufficient variation, along with improved covariation algorithms, are providing us with the opportunity to identify new base triples in 16S rRNA. The three-dimensional conformations for one of our strongest candidates involving U121 (C124:G237) and/or U121 (U125:A236) (Escherichia coli sequence and numbering) are analyzed here with different molecular modeling tools. Molecular modeling shows that U121 interacts with C124 in the U121 (C124:G237) base triple. This arrangement maintains isomorphic structures for the three most frequent sequence motifs (approximately 93% of known bacterial and archaeal sequences), is consistent with chemical reactivity of U121 in E. coli ribosomes, and is geometrically favorable. Further, the restricted set of observed canonical (GU, AU, GC) base-pair types at positions 124:237 and 125:236 is consistent with the fact that the canonical base-pair sets (for both base pairs) that are not observed in nature prevent the formation of the 121 (124:237) base triple. The analysis described here serves as a general scheme for the prediction of specific secondary and tertiary structure base pairing where there is a network of correlated base changes.

Asunto(s)

RESUMEN

Comparative sequence analysis reveals a coordinated set of nucleotide exchanges between the base pair 1092/1099 and the unpaired position 1072 [(1092/1099)1072] in the L11 binding domain of 23S ribosomal RNA. This set of exchanges has occurred at least 4 times during evolution, suggesting that these positions form a base triple. The analysis further suggests an important role for positions (1065/1073), adjacent to 1072. The covariation at positions (1092/1099)1072 is studied here by analysis of RNA variants using UV melting and binding of ribosomal protein L11 and thiostrepton to assay for tertiary folding of this domain. The tertiary structure of the RNA is eliminated by alteration of the unpaired nucleotide (C1072 to U mutation), and binding of L11 and thiostrepton are reduced 10-fold compared to the wild type. In contrast, substitution of the base pair (CG1092/1099 to UA mutation) allows formation of the tertiary structure but dramatically alters the pH dependence of tertiary folding. The fully compensated set of mutations, (CG)C to (UA)U, restores the tertiary structure of the RNA to a state almost identical to the wild type. The nature of this base triple and its implications for the folding of the RNA and ligand interactions are discussed.

Asunto(s)

RESUMEN

Adenosines are present at a disproportionately high frequency within several RNA structural motifs. To explore the importance of individual adenosine functional groups for group I intron activity, we performed Nucleotide Analog Interference Mapping (NAIM) with a collection of adenosine analogues. This paper reports the synthesis, transcriptional incorporation, and the observed interference pattern throughout the Tetrahymena group I intron for eight adenosine derivatives tagged with an alpha-phosphorothioate linkage for use in NAIM. All of the analogues were accurately incorporated into the transcript as an A. The sites that interfere with the 3'-exon ligation reaction of the Tetrahymena intron are coincident with the sites of phylogenetic conservation, yet the interference patterns for each analogue are different. These interference data provide several biochemical constraints that improve our understanding of the Tetrahymena ribozyme structure. For example, the data support an essential A-platform within the J6/6a region, major groove packing of the P3 and P7 helices, minor groove packing of the P3 and J4/5 helices, and an axial model for binding of the guanosine cofactor. The data also identify several essential functional groups within a highly conserved single-stranded region in the core of the intron (J8/7). At four sites in the intron, interference was observed with 2'-fluoro A, but not with 2'-deoxy A. Based upon comparison with the P4-P6 crystal structure, this may provide a biochemical signature for nucleotide positions where the ribose sugar adopts an essential C2'-endo conformation. In other cases where there is interference with 2'-deoxy A, the presence or absence of 2'-fluoro A interference helps to establish whether the 2'-OH acts as a hydrogen bond donor or acceptor. Mapping of the Tetrahymena intron establishes a basis set of information that will allow these reagents to be used with confidence in systems that are less well understood.

Asunto(s)

RESUMEN

Tertiary interactions are important in the higher-order folding of catalytic RNAs. Recently, a base triple, joining the two major domains of the catalytic core, was determined in group I introns from the cyanobacterium Anabaena PCC7120 and the eukaryote Tetrahymena thermophila. This base triple involves the fifth base pair of P4 and the fifth base of the single-stranded region J8/7. We made base pair and single-nucleotide substitutions in the fifth base pair of P4, a G-C in the wild-type Anabaena intron, and tested them for self-splicing activity. The results suggest a hydrogen bonding model in which only the C of the base pair interacts directly with the fifth base of J8/7. Comparative sequence analysis was used to determine the different combinations of base triples that occur in approximately 450 natural group I introns identified to date. About 94% of the base triples analyzed are compatible with the proposed hydrogen bonding model. Disrupting this base triple in the Tetrahymena intron resulted in the disappearance of splicing intermediates (intron 3' exon and 5' exon), even though the first step of splicing was not affected. Restoration of the base triple by a compensatory mutation reverted the intermediates to wild-type levels. These results suggest that disruption of the base triple increases the rate of the second step of splicing or of a conformational change preceding the second step. Repositioning of the base triple to form a new set of interactions may be required for the second step of splicing.

Asunto(s)

RESUMEN

The human malaria parasite Plasmodium vivax has been shown to regulate the transcription of two distinct 18 RNAs during development. Here we show a third and distinctive type of ribosome that is present shortly after zygote formation, a transcriptional pattern of ribosome types that relates closely to the developmental state of the parasite and a phenomenon that separates ribosomal types at a critical phase of maturation. The A-type ribosome is predominantly found in infected erythrocytes of the vertebrate and the mosquito blood meal. Transcripts from the A gene are replaced by transcripts from another locus, the O gene, shortly after fertilization and increase in number as the parasite develops on the mosquito midgut. Transcripts from another locus, the S gene, begins as the oocyst form of the parasite matures. RNA transcripts from the S gene are preferentially included in sporozoites that bud off from the oocyst and migrate to the salivary gland while the O gene transcripts are left within the oocyst. Although all three genes are typically eukaryotic in structure, the O gene transcript, described here, varies from the other two in core regions of the rRNA that are involved in mRNA decoding and translational termination. We now can correlate developmental progression of the parasite with changes in regions of rRNA sequence that are broadly conserved, where sequence alterations have been related to function in other systems and whose effects can be studied outside of Plasmodium. This should allow assessment of the role of translational control in parasite development.

Asunto(s)