Transient RNA structures underlie highly pathogenic avian influenza virus genesis

Highly pathogenic avian influenza viruses (HPAIVs) cause severe disease and high fatality in poultry1. They emerge exclusively from H5 and H7 low pathogenic avian influenza viruses (LPAIVs)2. Although insertion of a furin-cleavable multibasic cleavage site (MBCS) in the hemagglutinin gene was identified decades ago as the genetic basis for LPAIV-to-HPAIV transition3,4, the exact mechanisms underlying said insertion have remained unknown. Here we used an innovative combination of bioinformatic models to predict RNA structures forming around the influenza virus RNA polymerase during replication, and circular sequencing5 to reliably detect nucleotide insertions. We show that transient H5 hemagglutinin RNA structures predicted to trap the polymerase on purine-rich sequences drive nucleotide insertions characteristic of MBCSs, providing the first strong empirical evidence of RNA structure involvement in MBCS acquisition. Insertion frequencies at the H5 cleavage site were strongly affected by substitutions in flanking genomic regions altering predicted transient RNA structures. Introduction of H5-like cleavage site sequences and structures into an H6 hemagglutinin resulted in MBCS-yielding insertions never observed before in H6 viruses. Our results demonstrate that nucleotide insertions that underlie H5 HPAIV emergence result from a previously unknown RNA-structure-driven diversity-generating mechanism, which could be shared with other RNA viruses.

Conserved structured domains in plant non-coding RNA enod40, their evolution and recruitment of sequences from transposable elements

Plant long noncoding RNA enod40 is involved in the regulation of symbiotic associations with bacteria, in particular, in nitrogen-fixing root nodules of legumes, and with fungi in phosphate-acquiring arbuscular mycorrhizae formed by various plants. The presence of enod40 genes in plants that do not form such symbioses indicates its other roles in cell physiology. The molecular mechanisms of enod40 RNA function are poorly understood. Enod40 RNAs form several structured domains, conserved to different extents. Due to relatively low sequence similarity, identification of enod40 sequences in plant genomes is not straightforward, and many enod40 genes remain unannotated even in complete genomes. Here, we used comparative structure analysis and sequence similarity searches in order to locate enod40 genes and determine enod40 RNA structures in nitrogen-fixing clade plants and in grasses. The structures combine conserved features with considerable diversity of structural elements, including insertions of structured domain modules originating from transposable elements. Remarkably, these insertions contain sequences similar to tandem repeats and several stem-loops are homologous to microRNA precursors.

Hemagglutinin Subtype Specificity and Mechanisms of Highly Pathogenic Avian Influenza Virus Genesis

Highly Pathogenic Avian Influenza Viruses (HPAIVs) arise from low pathogenic precursors following spillover from wild waterfowl into poultry populations. The main virulence determinant of HPAIVs is the presence of a multi-basic cleavage site (MBCS) in the hemagglutinin (HA) glycoprotein. The MBCS allows for HA cleavage and, consequently, activation by ubiquitous proteases, which results in systemic dissemination in terrestrial poultry. Since 1959, 51 independent MBCS acquisition events have been documented, virtually all in HA from the H5 and H7 subtypes. In the present article, data from natural LPAIV to HPAIV conversions and experimental in vitro and in vivo studies were reviewed in order to compile recent advances in understanding HA cleavage efficiency, protease usage, and MBCS acquisition mechanisms. Finally, recent hypotheses that might explain the unique predisposition of the H5 and H7 HA sequences to obtain an MBCS in nature are discussed.

In Silico Analyses of the Role of Codon Usage at the Hemagglutinin Cleavage Site in Highly Pathogenic Avian Influenza Genesis

A vast diversity of 16 influenza hemagglutinin (HA) subtypes are found in birds. Interestingly, viruses from only two subtypes, H5 and H7, have so far evolved into highly pathogenic avian influenza viruses (HPAIVs) following insertions or substitutions at the HA cleavage site by the viral polymerase. The mechanisms underlying this striking subtype specificity are still unknown. Here, we compiled a comprehensive dataset of 20,488 avian influenza virus HA sequences to investigate differences in nucleotide and amino acid usage at the HA cleavage site between subtypes and how these might impact the genesis of HPAIVs by polymerase stuttering and realignment. We found that sequences of the H5 and H7 subtypes stand out by their high purine content at the HA cleavage site. In addition, fewer substitutions were necessary in H5 and H7 HAs than in HAs from other subtypes to acquire an insertion-prone HA cleavage site sequence, as defined based on in vitro and in vivo data from the literature. Codon usage was more favorable for HPAIV genesis in sequences of viruses isolated from species or geographical regions in which HPAIV genesis is more frequently observed in nature. The results of the present analyses suggest that the subtype restriction of HPAIV genesis to H5 and H7 influenza viruses might be due to the particular codon usage at the HA cleavage site in these subtypes.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions

Motivation

The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3′ untranslated regions (3′UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies.

Results

We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3′UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3′UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses.

Availability and implementation

The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Xrn1-resistant RNA structures are well-conserved within the genus flavivirus

Subgenomic RNAs are produced by several RNA viruses through incomplete degradation of their genomic RNA by the exoribonuclease Xrn1, and have been shown to be essential for viral growth and pathogenicity. Within the flavivirus genus of the Flaviviridae family, two distinct classes of Xrn1-resistant RNA motifs have been proposed; one for mosquito-borne and insect-specific flaviviruses, and one for tick-borne flaviviruses and no-known-vector flaviviruses. We investigated tick-borne and no-known-vector flavivirus Xrn1-resistant RNA motifs through systematic in vitro mutational analysis and showed that both classes actually possess very similar structural configurations, including a double pseudoknot and a base-triple at identical, conserved locations. For the no-known-vector flavivirus Modoc virus, we show that in vivo generation of subgenomic flaviviral RNA was affected by mutations targeted at nucleotides involved in the structural features of flaviviral Xrn1-resistant RNA motifs that were defined in this work. Our results suggest that throughout the genus flavivirus Xrn1-resistant RNA motifs adopt the same topologically conserved structure.

All genera of Flaviviridae host a conserved Xrn1-resistant RNA motif

After infection by flaviviruses like Zika and West Nile virus, eukaryotic hosts employ the well-conserved endoribonuclease Xrn1 to degrade the viral genomic RNA. Within the 3' untranslated regions, this enzyme encounters intricate Xrn1-resistant structures. This results in the accumulation of subgenomic flaviviral RNAs, an event that improves viral growth and aggravates viral pathogenicity. Xrn1-resistant RNAs have been established throughout the flaviviral genus, but not yet throughout the entire Flaviviridae family. In this work, we use previously determined characteristics of these structures to identify homologous sequences in many members of the genera pegivirus, hepacivirus and pestivirus. We used structural alignment and mutational analyses to establish that these sequences indeed represent Xrn1-resistant RNA and that they employ the general features of the flaviviral xrRNAs, consisting of a double pseudoknot formed by five base-paired regions stitched together by a crucial triple base interaction. Furthermore, we demonstrate that the pestivirus Bungowannah virus produces subgenomic RNA in vivo. Altogether, these results indicate that viruses make use of a universal Xrn1-resistant RNA throughout the Flaviviridae family.

Insertions of codons encoding basic amino acids in H7 hemagglutinins of influenza A viruses occur by recombination with RNA at hotspots near snoRNA binding sites

The presence of multiple basic amino acids in the protease cleavage site of the hemagglutinin (HA) protein is the main molecular determinant of virulence of highly pathogenic avian influenza (HPAI) viruses. Recombination of HA RNA with other RNA molecules of host or virus origin is a dominant mechanism of multibasic cleavage site (MBCS) acquisition for H7 subtype HA. Using alignments of HA RNA sequences from documented cases of MBCS insertion due to recombination, we show that such recombination with host RNAs is most likely to occur at particular hotspots in ribosomal RNAs (rRNAs), transfer RNAs (tRNAs), and viral RNAs. The locations of these hotspots in highly abundant RNAs indicate that RNA recombination is facilitated by the binding of small nucleolar RNA (snoRNA) near the recombination points.

Fungal metabarcoding data integration framework for the MycoDiversity DataBase (MDDB)

Fungi have crucial roles in ecosystems, and are important associates for many organisms. They are adapted to a wide variety of habitats, however their global distribution and diversity remains poorly documented. The exponential growth of DNA barcode information retrieved from the environment is assisting considerably the traditional ways for unraveling fungal diversity and detection. The raw DNA data in association to environmental descriptors of metabarcoding studies are made available in public sequence read archives. While this is potentially a valuable source of information for the investigation of Fungi across diverse environmental conditions, the annotation used to describe environment is heterogenous. Moreover, a uniform processing pipeline still needs to be applied to the available raw DNA data. Hence, a comprehensive framework to analyses these data in a large context is still lacking. We introduce the MycoDiversity DataBase, a database which includes public fungal metabarcoding data of environmental samples for the study of biodiversity patterns of Fungi. The framework we propose will contribute to our understanding of fungal biodiversity and aims to become a valuable source for large-scale analyses of patterns in space and time, in addition to assisting evolutionary and ecological research on Fungi.

Conserved structural RNA domains in regions coding for cleavage site motifs in hemagglutinin genes of influenza viruses

The acquisition of a multibasic cleavage site (MBCS) in the hemagglutinin (HA) glycoprotein is the main determinant of the conversion of low pathogenic avian influenza viruses into highly pathogenic strains, facilitating HA cleavage and virus replication in a broader range of host cells. In nature, substitutions or insertions in HA RNA genomic segments that code for multiple basic amino acids have been observed only in the HA genes of two out of sixteen subtypes circulating in birds, H5 and H7. Given the compatibility of MBCS motifs with HA proteins of numerous subtypes, this selectivity was hypothesized to be determined by the existence of specific motifs in HA RNA, in particular structured domains. In H5 and H7 HA RNAs, predictions of such domains have yielded alternative conserved stem-loop structures with the cleavage site codons in the hairpin loops. Here, potential RNA secondary structures were analyzed in the cleavage site regions of HA segments of influenza viruses of different types and subtypes. H5- and H7-like stem-loop structures were found in all known influenza A virus subtypes and in influenza B and C viruses with homology modeling. Nucleotide covariations supported this conservation to be determined by RNA structural constraints that are stronger in the domain-closing bottom stems as compared to apical parts. The structured character of this region in (sub-)types other than H5 and H7 indicates its functional importance beyond the ability to evolve toward an MBCS responsible for a highly pathogenic phenotype.

Structural features of an Xrn1-resistant plant virus RNA

Xrn1 is a major 5ʹ-3ʹ exoribonuclease involved in the RNA metabolism of many eukaryotic species. RNA viruses have evolved ways to thwart Xrn1 in order to produce subgenomic non-coding RNA that affects the hosts RNA metabolism. The 3ʹ untranslated region of several beny- and cucumovirus RNAs harbors a so-called ‘coremin’ motif that is required for Xrn1 stalling. The structural features of this motif have not been studied in detail yet. Here, by using in vitro Xrn1 degradation assays, we tested over 50 different RNA constructs based on the Beet necrotic yellow vein virus sequence to deduce putative structural features responsible for Xrn1 stalling. We demonstrated that the minimal benyvirus stalling site consists of two hairpins of 3 and 4 base pairs respectively. The 5ʹ proximal hairpin requires a YGAD (Y = U/C, D = G/A/U) consensus loop sequence, whereas the 3ʹ proximal hairpin loop sequence is variable. The sequence of the 10-nucleotide spacer that separates the hairpins is highly conserved and potentially involved in tertiary interactions. Similar coremin motifs were identified in plant virus isolates from other families including Betaflexiviridae, Virgaviridae, Potyviridae and Secoviridae (order of the Picornavirales). We conclude that Xrn1-stalling motifs are more widespread among RNA viruses than previously realized.

A compensatory mutagenesis study of a conserved hairpin in the M gene segment of influenza A virus shows its role in virus replication

RNA structures are increasingly recognized to be of importance during influenza A virus replication. Here, we investigated a predicted conserved hairpin in the M gene segment (nt 967-994) within the region of the vRNA 5′ packaging signal. The existence of this RNA structure and its possible role in virus replication was investigated using a compensatory mutagenesis approach. Mutations were introduced in the hairpin stem, based on natural variation. Virus replication properties were studied for the mutant viruses with disrupted and restored RNA structures. Viruses with structure-disrupting mutations had lower virus titers and a significantly reduced median plaque size when compared with the wild-type (WT) virus, while viruses with structure restoring-mutations replicated comparable to WT. Moreover, virus replication was also reduced when mutations were introduced in the hairpin loop, suggesting its involvement in RNA interactions. Northern blot and FACS experiments were performed to study differences in RNA levels as well as production of M1 and M2 proteins, expressed via alternative splicing. Stem-disruptive mutants caused lower vRNA and M2 mRNA levels and reduced M2 protein production at early time-points. When the RNA structure was restored, vRNA, M2 mRNA and M2 protein levels were increased, demonstrating a compensatory effect. Thus, this study provides evidence for functional importance of the predicted M RNA structure and suggests its role in splicing regulation.

Optimisations and Challenges Involved in the Creation of Various Bioluminescent and Fluorescent Influenza A Virus Strains for In Vitro and In Vivo Applications

Bioluminescent and fluorescent influenza A viruses offer new opportunities to study influenza virus replication, tropism and pathogenesis. To date, several influenza A reporter viruses have been described. These strategies typically focused on a single reporter gene (either bioluminescent or fluorescent) in a single virus backbone. However, whilst bioluminescence is suited to in vivo imaging, fluorescent viruses are more appropriate for microscopy. Therefore, the idea l reporter virus varies depending on the experiment in question, and it is important that any reporter virus strategy can be adapted accordingly. Herein, a strategy was developed to create five different reporter viruses in a single virus backbone. Specifically, enhanced green fluorescent protein (eGFP), far-red fluorescent protein (fRFP), near-infrared fluorescent protein (iRFP), Gaussia luciferase (gLUC) and firefly luciferase (fLUC) were inserted into the PA gene segment of A/PR/8/34 (H1N1). This study provides a comprehensive characterisation of the effects of different reporter genes on influenza virus replication and reporter activity. In vivo reporter gene expression, in lung tissues, was only detected for eGFP, fRFP and gLUC expressing viruses. In vitro, the eGFP-expressing virus displayed the best reporter stability and could be used for correlative light electron microscopy (CLEM). This strategy was then used to create eGFP-expressing viruses consisting entirely of pandemic H1N1, highly pathogenic avian influenza (HPAI) H5N1 and H7N9. The HPAI H5N1 eGFP-expressing virus infected mice and reporter gene expression was detected, in lung tissues, in vivo. Thus, this study provides new tools and insights for the creation of bioluminescent and fluorescent influenza A reporter viruses.

RNA structural constraints in the evolution of the influenza A virus genome NP segment

Conserved RNA secondary structures were predicted in the nucleoprotein (NP) segment of the influenza A virus genome using comparative sequence and structure analysis. A number of structural elements exhibiting nucleotide covariations were identified over the whole segment length, including protein-coding regions. Calculations of mutual information values at the paired nucleotide positions demonstrate that these structures impose considerable constraints on the virus genome evolution. Functional importance of a pseudoknot structure, predicted in the NP packaging signal region, was confirmed by plaque assays of the mutant viruses with disrupted structure and those with restored folding using compensatory substitutions. Possible functions of the conserved RNA folding patterns in the influenza A virus genome are discussed.

Tospovirus ambisense genomic RNA segments use almost complete repertoire of stable tetraloops in the intergenic region

The intergenic regions of the ambisense RNA segments of viruses from the Tospovirus genus form large extended RNA structures that regulate virus replication. Using comparative structure analysis, we show the presence of conserved alternative conformations at the apical parts of these structures. In one conformation, a branched Y-shape, the 5'-proximal hairpin arms are mostly capped by exceptionally stable tetraloop motifs. The tetraloop hairpins are folded in both virus and virus-complementary sense RNAs, and different tetraloops can functionally replace each other. Folding simulations show that the branched Y-shape structures can undergo a conformational transition to alternative extended rod-like conformations. Functional importance of both alternatives is supported by nucleotide covariations. The balanced equilibrium between alternative structures is evidenced by native gel electrophoresis of mutant RNA transcripts with shifted equilibria. The tetraloops play a role in the stability and dynamics of structures but may also be recognized by proteins involved in translation and/or replication.

Upstream start codon in segment 4 of North American H2 avian influenza A viruses

H2N2 influenza A virus was the cause of the 1957 pandemic. Due to its constant presence in birds, the H2 subtype remains a topic of interest. In this work, comparison of H2 leader sequences of influenza A segment 4 revealed the presence of an upstream in-frame start codon in a majority of North American avian strains. This AUG is located seven codons upstream of the conventional start codon and is in a good Kozak context. In vivo experiments, using a luciferase reporter gene fused to leader sequences derived from North American avian H2 strains, support the efficient use of the upstream start codon. These results were corroborated by in vitro translation data using full-length segment 4 mRNA. Phylogenic analyses indicate that the upstream AUG, first detected in 1976, is stably nested in the North American avian lineage of H2 strains nowadays. The possible consequences of the upstream AUG are discussed.

Influenza virus RNA structure: unique and common features

The influenza A virus genome consists of eight negative-sense RNA segments. Here we review the currently available data on structure-function relationships in influenza virus RNAs. Various ideas and hypotheses about the roles of influenza virus RNA folding in the virus replication are also discussed in relation to other viruses.

A family of non-classical pseudoknots in influenza A and B viruses

A very conserved pseudoknot structure has been shown to fold in influenza virus RNA. The pseudoknot encompasses the 3' splice site of segment 8 RNA in both influenza A and B viruses. By sequence comparison of influenza virus strains, we derive a consensus motif that defines a novel RNA pseudoknot family. The orientation of the coaxially stacked stems in the influenza pseudoknot differs from that in classical H-pseudoknots. Apart from the size of the loops, the topology of the influenza pseudoknot resembles that of some long-range pseudoknotted conformations. A seed alignment of the influenza pseudoknot family, containing representative strain sequences together with a consensus structure description, has been submitted to the RNA families (Rfam) database.

PseudoBase++: an extension of PseudoBase for easy searching, formatting and visualization of pseudoknots

Pseudoknots have been recognized to be an important type of RNA secondary structures responsible for many biological functions. PseudoBase, a widely used database of pseudoknot secondary structures developed at Leiden University, contains over 250 records of pseudoknots obtained in the past 25 years through crystallography, NMR, mutational experiments and sequence comparisons. To promptly address the growing analysis requests of the researchers on RNA structures and bring together information from multiple sources across the Internet to a single platform, we designed and implemented PseudoBase++, an extension of PseudoBase for easy searching, formatting and visualization of pseudoknots. PseudoBase++ (http://pseudobaseplusplus.utep.edu) maps the PseudoBase dataset into a searchable relational database including additional functionalities such as pseudoknot type. PseudoBase++ links each pseudoknot in PseudoBase to the GenBank record of the corresponding nucleotide sequence and allows scientists to automatically visualize RNA secondary structures with PseudoViewer. It also includes the capabilities of fine-grained reference searching and collecting new pseudoknot information.

Identification of conserved secondary structures and expansion segments in enod40 RNAs reveals new enod40 homologues in plants

enod40 is a plant gene that participates in the regulation of symbiotic interaction between leguminous plants and bacteria or fungi. Furthermore, it has been suggested to play a general role in non-symbiotic plant development. Although enod40 seems to have multiple functions, being present in many land plants, the molecular mechanisms of its activity are unclear; they may be determined though, by short peptides and/or RNA structures encoded in the enod40 genes. We utilized conserved RNA structures in enod40 sequences to search nucleotide sequence databases and identified a number of new enod40 homologues in plant species that belong to known, but also, to yet unknown enod40-containing plant families. RNA secondary structure predictions and comparative sequence analysis of enod40 RNAs allowed us to determine the most conserved structural features, present in all known enod40 genes. Remarkably, the topology and evolution of one of the conserved structural domains are similar to those of the expansion segments found in structural RNAs such as rRNAs, RNase P and SRP RNAs. Surprisingly, the enod40 RNA structural elements are much more stronger conserved than the encoded peptides. This finding suggests that some general functions of enod40 gene could be determined by the encoded RNA structure, whereas short peptides may be responsible for more diverse functions found only in certain plant families.

An RNA conformational shift in recent H5N1 influenza A viruses

Recent outbreaks of avian influenza are being caused by unusually virulent H5N1 strains. It is unknown what makes these recent H5N1 strains more aggressive than previously circulating strains. Here, we have compared more than 3000 RNA sequences of segment 8 of type A influenza viruses and found a unique single nucleotide substitution typically associated with recent H5N1 strains. By phylogenetic analysis, biochemical and biophysical experiments, we demonstrate that this substitution dramatically affects the equilibrium between a hairpin and a pseudoknot conformation near the 3' splice-site of the NS gene. This conformational shift may have consequences for splicing regulation of segment 8 mRNA. Our data suggest that besides changes at the protein level, changes in RNA secondary structure should be seriously considered when attempting to explain influenza virus evolution.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computer simulation of chaperone effects of Archaeal C/D box sRNA binding on rRNA folding

Archaeal C/D box small RNAs (sRNAs) are homologues of eukaryotic C/D box small nucleolar RNAs (snoRNAs). Their main function is guiding 2'-O-ribose methylation of nucleotides in rRNAs. The methylation requires the pairing of an sRNA antisense element to an rRNA target site with formation of an RNA-RNA duplex. The temporary formation of such a duplex during rRNA maturation is expected to influence rRNA folding in a chaperone-like way, in particular in thermophilic Archaea, where multiple sRNAs with two binding sites are found. Here we investigate possible mechanisms of chaperone function of Archaeoglobus fulgidus and Pyrococcus abyssi C/D box sRNAs using computer simulations of rRNA secondary structure formation by genetic algorithm. The effects of sRNA binding on rRNA structure are introduced as temporary structural constraints during co-transcriptional folding. Comparisons of the final predictions with simulations without sRNA binding and with phylogenetic structures show that sRNAs with two antisense elements may significantly facilitate the correct formation of long-range interactions in rRNAs, in particular at elevated temperatures. The simulations suggest that the main mechanism of this effect is a transient restriction of folding in rRNA domains where the termini are brought together by binding to double-guide sRNAs.

Nemesia ring necrosis virus: a new tymovirus with a genomic RNA having a histidylatable tobamovirus-like 3′ end

The complete nucleotide sequence of the genomic RNA of the new virus Nemesia ring necrosis virus (NeRNV), which is widespread in various ornamental plant species belonging to the Scrophulariaceae and Verbenaceae, has been determined. Based on its gene content, the folding properties of its 5′-untranslated region and in vitro translation experiments, NeRNV RNA is a typical tymovirus RNA. Its 3′ end, however, differs greatly from those of the valine-specific tymoviral RNAs that have been analysed previously. It can be folded into an upstream pseudoknot domain and a histidine-specific tRNA-like structure, a combination that, so far, has been found only in tobamoviral RNAs. The identity elements found in NeRNV RNA for recognition by yeast histidyl-tRNA synthetase are more similar to those of yeast tRNAHis than the ones found in tobacco mosaic virus RNA. As a result NeRNV RNA can be charged with histidine even more efficiently than tobacco mosaic virus RNA.

Discontinuous subgenomic RNA synthesis in arteriviruses is guided by an RNA hairpin structure located in the genomic leader region

Nidoviruses produce an extensive 3'-coterminal nested set of subgenomic (sg) mRNAs, which are used to express structural proteins and sometimes accessory proteins. In arteriviruses and coronaviruses, these mRNAs contain a common 5' leader sequence, derived from the genomic 5' end. The joining of the leader sequence to different segments derived from the 3'-proximal part of the genome (mRNA bodies) presumably involves a unique mechanism of discontinuous minus-strand RNA synthesis in which base pairing between sense and antisense transcription-regulating sequences (TRSs) plays an essential role. The leader TRS is present in the loop of a hairpin structure that functions in sg mRNA synthesis. In this study, the minimal sequences in the 5'-proximal region of the Equine arteritis virus genome that are required for sg RNA synthesis were delimited through mutagenesis. A full-length cDNA clone was engineered in which this domain was duplicated, allowing us to make mutations and monitor their effects on sg RNA synthesis without seriously affecting genome replication and translation. The leader TRS present in the duplicated sequence was used and yielded novel sg mRNAs with significantly extended leaders. Our combined findings suggest that the leader TRS hairpin (LTH) and its immediate flanking sequences are essential for efficient sg RNA synthesis and form an independent functional entity that could be moved 300 nucleotides downstream of its original position in the genome. We hypothesize that a conformational switch in the LTH region regulates the role of the 5'-proximal region of the arterivirus genome in subgenomic RNA synthesis.

Conservation of RNA structures enables TNV and BYDV 5' and 3' elements to cooperate synergistically in cap-independent translation

The subgenomic RNA 2 of tobacco necrosis virus A (TNV sgRNA2) encodes the viral coat protein, is unpolyadenylated and presumably uncapped. Here, we show that TNV sgRNA2 is translated cap independently. This cap-independent translation requires the leader and a 140 nt element of the trailer both in wheat germ extract and in tobacco protoplasts. Similar to barley yellow dwarf virus (BYDV), the TNV 5' and 3' elements stimulate translation synergistically. Computer-aided phylogenetic analysis of the secondary structure of the TNV trailer revealed that the 3' translation element is part of a major conserved stem-loop that contains similarities to structures in the BYDV 3' translation element. These data suggest that the translation mechanisms of TNV sgRNA2 and BYDV RNA are related. To further characterize this relationship, we tested whether cooperativity exists between TNV sgRNA2 and BYDV 5' and 3' elements. We found that the TNV sgRNA2 5' element stimulates translation synergistically with the BYDV 3' element in vitro. This finding is the first evidence for conservation of structures that enable a 5'-3' interaction stimulating cap-independent translation.

Secondary structure and function of the 5'-proximal region of the equine arteritis virus RNA genome

Nidoviruses produce an extensive 3'-coterminal nested set of subgenomic mRNAs, which are used to express their structural proteins. In addition, arterivirus and coronavirus mRNAs contain a common 5' leader sequence, derived from the genomic 5' end. The joining of this leader sequence to different segments (mRNA bodies) from the genomic 3'-proximal region presumably involves a unique mechanism of discontinuous minus-strand RNA synthesis. Key elements in this process are the so-called transcription-regulating sequences (TRSs), which determine a base-pairing interaction between sense and antisense viral RNA that is essential for leader-to-body joining. To identify RNA structures in the 5'-proximal region of the equine arteritis virus genome that may be involved in subgenomic mRNA synthesis, a detailed secondary RNA structure model was established using bioinformatics, phylogenetic analysis, and RNA structure probing. According to this structure model, the leader TRS is located in the loop of a prominent hairpin (leader TRS hairpin; LTH). The importance of the LTH was supported by the results of a mutagenesis study using an EAV molecular clone. Besides evidence for a direct role of the LTH in subgenomic RNA synthesis, indications for a role of the LTH region in genome replication and/or translation were obtained. Similar LTH structures could be predicted for the 5'-proximal region of all arterivirus genomes and, interestingly, also for most coronaviruses. Thus, we postulate that the LTH is a key structural element in the discontinuous subgenomic RNA synthesis and is likely critical for leader TRS function.

Structural motifs in the RNA encoded by the early nodulation gene enod40 of soybean

The plant gene enod40 is highly conserved among legumes and also present in various non-legume species. It is presumed to play a central regulatory role in the Rhizobium-legume interaction, being expressed well before the initiation of cortical cell divisions resulting in nodule formation. Two small peptides encoded by enod40 mRNA as well as its secondary structure have been shown to be key elements in the signalling processes underlying nodule organogenesis. Here results concerning the secondary structure of mRNA of enod40 in soybean are presented. This study combined a theoretical approach, involving structure prediction and comparison, as well as structure probing. Our study indicates five conserved domains in enod40 mRNA among numerous leguminous species. Structure comparison suggests that some domains are also conserved in non-leguminous species and that an additional domain exists that was found only in leguminous species developing indeterminate nodules. Enzymatic and chemical probing data support the structure for three of the domains, and partially for the remaining two. The rest of the molecule appears to be less structured. Some of the domains include motifs, such as U-containing internal loops and bulges, which seem to be conserved. Therefore, they might be involved in the regulatory role of enod40 RNA.

Structural variation and functional importance of a D-loop-T-loop interaction in valine-accepting tRNA-like structures of plant viral RNAs

Valine-accepting tRNA-like structures (TLSs) are found at the 3' ends of the genomic RNAs of most plant viruses belonging to the genera Tymovirus, Furovirus, Pomovirus and Pecluvirus, and of one Tobamovirus species. Sequence alignment of these TLSs suggests the existence of a tertiary D-loop-T-loop interaction consisting of 2 bp, analogous to those in the elbow region of canonical tRNAs. The conserved G(18).Psi(55) pair of regular tRNAs is found to covary in these TLSs between G.U (possibly also modified to G.Psi) and A.G. We have mutated the relevant bases in turnip yellow mosaic virus (TYMV) and examined the mutants for symptom development on Chinese cabbage plants and for accumulation of genetic reversions. Development of symptoms is shown to rely on the presence of either A.G or G.U in the original mutants or in revertants. This finding supports the existence and functional importance of this tertiary interaction. The fact that only G.U and A.G are accepted at this position appears to result from steric and energetic limitations related to the highly compact nature of the elbow region. We discuss the implications of these findings for the various possible functions of the valine-accepting TLS.

A pH-jump approach for investigating secondary structure refolding kinetics in RNA

It has been shown that premature translation of the plasmid-mediated toxin in hok/sok of plasmid R1 and pnd/pndB of plasmid R483 is prevented during transcription of the hok and pnd mRNAs by the formation of metastable hairpins at the 5'-end of the mRNA. Here, an experimental approach is presented, which allows the accurate measurement of the refolding kinetics of the 5'-end RNA fragments in vitro without chemically modifying the RNA. The method is based on acid denaturation followed by a pH-jump to neutral pH as a novel way to trap kinetically favoured RNA secondary structures, allowing the measurement of a wide range of biologically relevant refolding rates, with or without the use of standard stopped-flow equipment. The refolding rates from the metastable to the stable conformation in both the hok74 and pnd58 5'-end RNA fragments were determined by using UV absorbance changes corresponding to the structural rearrangements. The measured energy barriers showed that the refolding path does not need complete unfolding of the metastable structures before the formation of the final structures. Two alternative models of such a pathway are discussed.

Functional evidence for D- and T-loop interactions in tmRNA

During bacterial protein synthesis, stalled ribosomes can be rescued by tmRNA, a molecule with both tRNA and mRNA features. The tRNA region of tmRNA has sequence similarity with tRNA(Ala) and also has a clover-leaf structure folded similarly as in canonical tRNAs. Here we propose the L-shape of tmRNA to be stabilized by two tertiary interactions between its D- and T-loop on the basis of phylogenetic and experimental evidence. Mutational analysis clearly demonstrates a tertiary interaction between G(13) and U(342). Strikingly, this in evolution conserved interaction is not primarily important for tmRNA alanylation and for binding to elongation factor Tu, but especially for a proper functioning of SmpB.

The 5' and 3' extremities of the satellite tobacco necrosis virus translational enhancer domain contribute differentially to stimulation of translation

The translational enhancer domain (TED) of satellite tobacco necrosis virus (STNV) RNA stimulates translation of uncapped RNAs autonomously. Here we set out to identify the 5' and 3' extremities of TED and features of these sequences with respect to translation. We found that both in wheat germ extract and in tobacco protoplasts, the 5' border is confined to 3 nt. Mutational analysis revealed that the autonomous function of TED is sensitive to 5' flanking sequences. At the 3' end of TED, 23 nt have a cumulative, quantitative effect on translation in wheat germ extract, whereas in tobacco protoplasts, the most 3' 14 nt of these 23 nt do not enhance translation. The 5' and 3' sequence requirements triggered the development of a new secondary structure model. In this model, TED folds into a phylogenetically conserved stem-loop structure in which the essential 5' nucleotides base-pair with the 3' nucleotides that stimulate translation both in vitro and in vivo. Importantly, the 14 3' nucleotides in TED that stimulate translation in the wheat germ extract only do not require the predicted base-pairing in order to function. The discrepancy between in vitro and in vivo sequence requirements thus correlates with potential base-pairing requirements, opening the possibility that TED contains two functional domains.

Selective pressures on RNA hairpins in vivo and in vitro

Comparison of the most stable potential hairpins in the sequences of natural ribozymes with those in the randomized sequences has revealed that the hairpin loop energies are lower than expected by chance. Although these hairpins are not necessarily parts of functional structures, there is a selective pressure to diminish the destabilizing free energies of the hairpin loops. In contrast, no significant bias is observed in the stacking values of the most stable stems. In the ribozymes isolated in vitro the loops of potential hairpins are closer to random values, which can result in less efficient folding rates. Furthermore, the effects of kinetic traps seem to be more significant in the folding pathways of the in vitro isolates due to a potential to form stable stacks incompatible with the functional folds. Similarly to natural ribozyme sequences, the untranslated regions of viral RNAs also form hairpins with relatively low loop free energies. These evolutionary trends suggest ways for efficient engineering of improved RNA constructs on the basis of analysis of in vitro isolates and approaches for the search of regions coding for functional RNA structures in large genome sequences.

PseudoBase: structural information on RNA pseudoknots

PseudoBase is a database containing structural, functional and sequence data related to RNA pseudo-knots. It can be reached at http://wwwbio.LeidenUniv.nl/ approximately Batenburg/PKB.html. For each pseudoknot, thirteen items are stored, for example the relevant sequence, the stem positions of the pseudoknot, the EMBL accession number of the sequence and the support that can be given regarding the reliability of the pseudo-knot. Since the last publication, information on sizes of the stems and the loops in the pseudoknots has been added. Also added are alternative entries that produce surveys of where the pseudoknots are, sorted according to stem size or loop size.

Genetic manipulation of arterivirus alternative mRNA leader-body junction sites reveals tight regulation of structural protein expression

To express its structural proteins, the arterivirus Equine arteritis virus (EAV) produces a nested set of six subgenomic (sg) RNA species. These RNA molecules are generated by a mechanism of discontinuous transcription, during which a common leader sequence, representing the 5' end of the genomic RNA, is attached to the bodies of the sg RNAs. The connection between the leader and body parts of an mRNA is formed by a short, conserved sequence element termed the transcription-regulating sequence (TRS), which is present at the 3' end of the leader as well as upstream of each of the structural protein genes. With the exception of RNA3, only one body TRS was previously assumed to be used to join the leader and body of each EAV sg RNA. Here we show that for the synthesis of two other sg RNAs, RNA4 and RNA5, alternative leader-body junction sites that differ substantially in transcriptional activity are used. By site-directed mutagenesis of an EAV infectious cDNA clone, the alternative TRSs used to generate RNA3, -4, and -5 were inactivated, which strongly influenced the corresponding RNA levels and the production of infectious progeny virus. The relative amounts of RNA produced from alternative TRSs differed significantly and corresponded to the relative infectivities of the virus mutants. This strongly suggested that the structural proteins that are expressed from these RNAs are limiting factors during the viral life cycle and that the discontinuous step in sg RNA synthesis is crucial for the regulation of their expression. On the basis of a theoretical analysis of the predicted RNA structure of the 3' end of the EAV genome, we propose that the local secondary RNA structure of the body TRS regions is an important factor in the regulation of the discontinuous step in EAV sg mRNA synthesis.

Coupled nucleotide covariations reveal dynamic RNA interaction patterns

Evolutionarily conserved structures in related RNA molecules contain coordinated variations (covariations) of paired nucleotides. Analysis of covariations is a very powerful approach to deduce phylogenetically conserved (i.e., functional) conformations, including tertiary interactions. Here we discuss conserved RNA folding pathways that are revealed by covariation patterns. In such pathways, structural requirements for alternative pairings cause some nucleotides to covary with two different partners. Such "coupled" covariations between three or more nucleotides were found in various types of RNAs. The analysis of coupled covariations can unravel important features of RNA folding dynamics and improve phylogeny reconstruction in some cases. Importantly, it is necessary to distinguish between multiple covariations determined by mutually exclusive structures and those determined by tertiary contacts.

PseudoBase: a database with RNA pseudoknots

PseudoBase is a database containing structural, functional and sequence data related to RNA pseudo-knots. It can be reached at http://wwwbio. Leiden Univ.nl/ approximately Batenburg/PKB.html. This page will direct the user to a retrieval page from where a particular pseudoknot can be chosen, or to a submission page which enables the user to add pseudoknot information to the database or to an informative page that elaborates on the various aspects of the database. For each pseudoknot, 12 items are stored, e.g. the nucleotides of the region that contains the pseudoknot, the stem positions of the pseudoknot, the EMBL accession number of the sequence that contains this pseudoknot and the support that can be given regarding the reliability of the pseudoknot. Access is via a small number of steps, using 16 different categories. The development process was done by applying the evolutionary methodology for software development rather than by applying the methodology of the classical waterfall model or the more modern spiral model.

Metastable structures and refolding kinetics in hok mRNA of plasmid R1

Programmed cell death by hok/sok of plasmid R1 and pnd/pndB of R483 mediates plasmid maintenance by killing of plasmid-free cells. It has been previously suggested that premature translation of the plasmid-mediated toxin is prevented during transcription of the hok and pnd mRNAs by the formation of metastable hairpins in the mRNA at the 5' end. Here, experimental evidence is presented for the existence of metastable structures in the 5' leader of the hok and pnd mRNAs in vitro. The kinetics of refolding from the metastable to the stable structure in the isolated fragments of the 5' ends of both the hok and pnd mRNAs could be estimated, in agreement with the structural rearrangement in this region, as predicted to occur during transcription and mRNA activation. The refolding rates of hok and pnd structures are slow enough to allow for the formation of downstream hairpin structures during elongation of the mRNAs, which thereby helps to stabilize the metastable structures. Thus, the kinetic refolding parameters of the hok and pnd mRNAs are consistent with the proposal that the metastable structures prevent premature translation and/or antisense RNA binding during transcription.

Arterivirus discontinuous mRNA transcription is guided by base pairing between sense and antisense transcription-regulating sequences

To generate an extensive set of subgenomic (sg) mRNAs, nidoviruses (arteriviruses and coronaviruses) use a mechanism of discontinuous transcription. During this process, mRNAs are generated that represent the genomic 5' sequence, the so-called leader RNA, fused at specific positions to different 3' regions of the genome. The fusion of the leader to the mRNA bodies occurs at a short, conserved sequence element, the transcription-regulating sequence (TRS), which precedes every transcription unit in the genome and is also present at the 3' end of the leader sequence. Here, we have used site-directed mutagenesis of the infectious cDNA clone of the arterivirus equine arteritis virus to show that sg mRNA synthesis requires a base-pairing interaction between the leader TRS and the complement of a body TRS in the viral negative strand. Mutagenesis of the body TRS of equine arteritis virus RNA7 reduced sg RNA7 transcription severely or abolished it completely. Mutations in the leader TRS dramatically influenced the synthesis of all sg mRNAs. The construction of double mutants in which a mutant leader TRS was combined with the corresponding mutant RNA7 body TRS resulted in the specific restoration of mRNA7 synthesis. The analysis of the mRNA leader-body junctions of a number of mutants with partial transcriptional activity provided support for a mechanism of discontinuous minus-strand transcription that resembles similarity-assisted, copy-choice RNA recombination.

An approximation of loop free energy values of RNA H-pseudoknots

A set of free energy values is suggested for RNA H-pseudoknot loops. The parameters are adjusted to be consistent with the theory of polymer thermodynamics and known data on pseudoknots. The values can be used for estimates of pseudoknot stabilities and computer predictions of RNA structures.

The mechanism of actinomycin D-mediated inhibition of HIV-1 reverse transcription

The mechanism of reverse transcription was analyzed in vitro with RNA templates and the reverse transcriptase (RT) enzyme of human immunodeficiency virus type 1 (HIV-1). In particular, we analyzed the mechanism of actinomycin D (ActD) mediated inhibition of the strand transfer step, in which the newly synthesized cDNA, termed the (-) strand strong stop or (-)ssDNA, is transferred from the donor RNA onto the acceptor RNA. This strand transfer reaction is a rather inefficient process in vitro. We found that this is in part due to the presence of an excess donor RNA, and highly efficient strand transfer was achieved by reducing the amount of donor RNA. We suggest that annealing of the (-)ssDNA to the excess donor RNA is preferred over productive binding to the acceptor RNA because of a higher basepair complementarity. ActD remains a potent inhibitor of strand transfer in this optimized assay system. We measured no effect of ActD on the elongation of reverse transcription or the RNase H action of the RT enzyme. Instead, we provide evidence that ActD acts through direct interaction with the (-)ssDNA, thereby blocking the basepairing capacity of this molecule. The possible use of single-stranded DNA binding molecules as antiretroviral agents is discussed.

The nucleotide sequence of satellite tobacco necrosis virus strain C and helper-assisted replication of wild-type and mutant clones of the virus

The complete nucleotide sequence of satellite tobacco necrosis virus strain C (STNV-C) was determined. The genome has a similar overall organization to two STNV isolates studied previously but differs significantly from them in the secondary structure of the translated and untranslated regions (UTRs). STNV-C RNA is naturally uncapped and contains 1221 nt: 101 nt in the 5' UTR, 606 nt in the capsid protein (CP) coding region and 514 nt in the 3' UTR. Using the known sequences of STNV-C and tobacco necrosis virus strain D (TNV-D) RNAs, full-length cDNA clones of both RNAs were constructed. Synthetic transcripts derived from STNV-C cDNA clones only replicated in plants and protoplasts when co-inoculated with TNV-D transcripts. A number of mutant clones in both the 3' and the 5' STNV-C RNA UTRs were constructed which disrupted putative cis-acting elements recognized by helper virus polymerase. Deletion analysis revealed an essential requirement of all 3' and 5' proximal sequences in the STNV-C UTRs for replication. However, an internal region in the 3' UTR could be deleted without loss of infectivity. Likewise, the entire STNV-C CP-encoding region could be deleted and replaced with a marker gene of a similar size without loss of transcript accumulation in plants.

Dynamic competition between alternative structures in viroid RNAs simulated by an RNA folding algorithm

The folding pathways of viroid RNAs were studied using computer simulations by the genetic algorithm for RNA folding. The folding simulations were performed for PSTVd RNAs of both polarities, using the wild-type sequence and some previously known mutants with suggested changes in the stable or metastable structures. It is shown that metastable multihairpin foldings in the minus strand replicative intermediates are established due to the specific folding pathway that ensures the absence of the most stable rod-like structure. Simulations of the PSTVd minus strand folding during transcription reveal a metastable hairpin, formed in the left terminal domain region of the PSTVd. Despite high sequence variability, this hairpin is conserved in all known large viroids of both subgroups of PSTVd type, and is presumably necessary to guide the folding of the HPII hairpin which is functional in the minus strand. The folding simulations are able to demonstrate the changes in the balance between metastable and stable structures in mutant PSTVd RNAs. The stable rod-like structure of the circular viroid (+) RNA is also folded via a dynamic folding pathway. Furthermore, the simulations show that intermediate steps in the forced evolution of a shortened PSTVd replicon may be reconstructed by a mechanistic model of different folding pathway requirements in plus- and minus-strand RNAs. Thus the formation of viroid RNA structure strongly depends on dynamics of competition between alternative RNA structures. This also suggests that the replication efficiency of viroid sequences may be estimated by a simulation of the folding process.

Antisense RNA-regulated programmed cell death

Eubacterial plasmids and chromosomes encode multiple killer genes belonging to the hok gene family. The plasmid-encoded killer genes mediate plasmid stabilization by killing plasmid-free cells. This review describes the genetics, molecular biology, and evolution of the hok gene family. The complicated antisense RNA-regulated control-loop that regulates posttranscriptional and postsegregational activation of killer mRNA translation in plasmid-free cells is described in detail. Nucleotide covariations in the mRNAs reveal metastable stem-loop structures that are formed at the mRNA 5' ends in the nascent transcripts. The metastable structures prevent translation and antisense RNA binding during transcription. Coupled nucleotide covariations provide evidence for a phylogenetically conserved mRNA folding pathway that involves sequential dynamic RNA rearrangements. Our analyses have elucidated an intricate mechanism by which translation of an antisense RNA-regulated mRNA can be conditionally activated. The complex phylogenetic relationships of the plasmid- and chromosome-encoded systems are also presented and discussed.

Programmed cell death by hok/sok of plasmid R1: processing at the hok mRNA 3'-end triggers structural rearrangements that allow translation and antisense RNA binding

The hok/sok locus of plasmid R1 mediates plasmid stabilization by killing of plasmid-free cells. The locus specifies two RNAs, hok mRNA and Sok antisense RNA. The post-segregational killing mediated by hok/sok is governed by a complicated control mechanism that involves both post-transcriptional inhibition of translation by Sok-RNA and activation of hok translation by mRNA 3' processing. Sok-RNA inhibits translation of a reading frame (mok) that overlaps with hok, and translation of hok is coupled to translation of mok. In the inactive full-length hok mRNA, the translational activator element at the mRNA 5'-end (tac) is sequestered by the fold-back-inhibitory element located at the mRNA 3'-end (fbi). The 5' to 3' pairing locks the RNA in an inert configuration in which the SDmok and Sok-RNA target regions are sequestered. Here we show that the 3' processing leads to major structural rearrangements in the mRNA 5'-end. The structure of the refolded RNA explains activation of translation and antisense RNA binding. The refolded RNA contains an antisense RNA target stem-loop that presents the target nucleotides in a single-stranded conformation. The stem of the target hairpin contains SDmok and AUGmok in a paired configuration. Using toeprinting analysis, we show that this pairing keeps SDmok in an accessible configuration. Furthermore, a mutational analysis shows that an internal loop in the target stem is prerequisite for efficient translation and antisense RNA binding.

Programmed cell death by hok/sok of plasmid R1: coupled nucleotide covariations reveal a phylogenetically conserved folding pathway in the hok family of mRNAs

The hok/sok system of plasmid R1 mediates plasmid maintenance by killing of plasmid-free cells. Translation of the stable toxin-encoding hok mRNA is repressed by the unstable Sok antisense RNA. Using genetic algorithm simulations and phylogenetic comparisons, we analyse five plasmid-encoded and two chromosome-encoded hok-homologous mRNAs. A similar folding pathway was found for all mRNAs. Metastable hairpins at the very 5'-ends of the mRNAs were predicted to prevent the formation of structures required for translation and antisense RNA binding. Thus the folding of the mRNA 5'-ends appears to explain the apparent inactivity of the nascent transcripts. In the full-length mRNAs, long-range 5' to 3' interactions were predicted in all cases. The 5' to 3' interactions lock the mRNAs in inactive configurations. Translation of the mRNAs is activated by 3' exonucleolytic processing. Simulation of the 3' processing predicted that it triggers rearrangements of the mRNA 5'-ends with the formation of translational activator and antisense RNA target hairpins. Alignment of the mRNA sequences revealed a large number of nucleotide covariations that support the existence of the proposed secondary structures. Furthermore, coupled covariations support the folding pathway and provide evidence that the mRNA 5'-ends pair with three different partners during the proposed series of dynamic RNA rearrangements.

Secondary structure model of the coat protein gene of turnip yellow mosaic virus RNA: long, C-rich, single-stranded regions

The RNA of all tymoviruses, a group of ssRNA plant viruses, has a base composition that is different from that of most other viruses. The excess of cytosines (35-42%) and the low number of guanosines (15-17%) must impel an RNA structure with a relatively low amount of base pairing and a high incidence of unpaired cytosines. These unpaired cytosines probably function in RNA-protein interactions. To gain Insight into the way the RNA is positioned inside the virion, the secondary structure has been determined of a part of TYMV RNA, including the so-called tymobox, the coat protein gene, and the 3' untranslated region, by structure probing, sequence comparison, and computer predictions. Conservation of secondary structure elements in tymoviruses is not high and does not parallel the conservation of the primary structure. A combination of structure prediction and probing experiments, however, results in a model consisting of structured domains of 100-200 nucleotides interspersed by long unpaired cytosine-rich regions. The latter may interact with the coat protein inside the virion. The structure of some functionally interesting regions of the 3' part of TYMV RNA is also discussed.

Phylogenetic relationships inferred from the sequence and secondary structure of ITS1 rRNA in Albinaria and putative Isabellaria species (Gastropoda, Pulmonata, Clausiliidae)

An analysis of the ITS1 sequence variation among five species of terrestrial pulmonate snails was performed to decide between two conflicting hypotheses concerning the phylogeny of these anatomically similar gastropods. It turned out that the so-called genus Isabellaria is a polyphyletic entity; the diagnostic, apomorphic structure of its clausilial apparatus, enabling a nearly complete obstruction of the shell aperture with the animal at rest, apparently evolved more than once from ancestors currently classified with the speciose genus Albinaria. The classification based on general shell shape and sculpture, and distributional patterns, turns out to be the natural one. This study also provides the first data on ITS1 sequences in gastropods. The recently published ITS1 sequence of another molluscan species, a bivalve, is quite different but similar in length to that of the snails.

The influence of a metastable structure in plasmid primer RNA on antisense RNA binding kinetics

Replication of the ColE1 group plasmids is kinetically regulated by the interaction between plasmid-encoded primer RNA II and antisense RNA I. The binding is dependent on alternative RNA II conformations, formed during the transcription, and effectively inhibits the primer function within some time interval. In this paper, the folding pathways for the wild type and copy number mutants of ColE1 RNA II are studied using simulations by a genetic algorithm. The simulated pathways reveal a transient formation of a metastable structure, which is stabilized by copy number mutations. The folding kinetics of the proposed conformational transitions is calculated using a model of a multistep refolding process with elementary steps of double-helical stem formation or disruption. The approximation shows that the lifetime of the metastable structure is relatively long and is considerably increased in the mutants, resulting in a delay of the formation of the stable RNA II structure, which is the most sensitive to the inhibition by the antisense RNA I. Thus the effect of copy number mutations can be interpreted as a compression of the time window of effective inhibition due to an increased time spent by the RNA II in the metastable state. The implications of metastable foldings in RNA functioning are discussed.

The computer simulation of RNA folding pathways using a genetic algorithm

A procedure for simulating the RNA folding process using the principles of genetic algorithm is proposed. The method allows one to simulate a folding pathway of RNA, including such processes as disruption of temporarily formed structures, the folding of a molecule during its synthesis and pseudoknot formation. The simulations are able to predict functional metastable foldings and kinetically driven transitions to more stable structures. The analysis of free energies for intermediate foldings allows estimation of the ranges of kinetic refolding barriers and suggests that in some RNAs the selective evolutionary pressure suppresses the possibilities for alternative structures that could form in the course of transcription. It is shown that the folding pathway simulation can result in structure predictions that are more consistent with phylogenetically proven structures than minimum energy solutions. This suggest that RNA folding kinetics is very important for the formation of functional RNA structures. Therefore, apart form its value for predictions of RNA structures, the proposed computer simulations can be a powerful tool in the studies of RNA folding features.

An APL-programmed genetic algorithm for the prediction of RNA secondary structure

The possibilities of using a genetic algorithm for the prediction of RNA secondary structure were investigated. The algorithm, using the procedure of stepwise selection of the most fit structures (similarly to natural evolution), allows different models of fitness or driving forces determining RNA structure to be easily introduced. This can be used for simulation of the RNA folding process and for the investigation of possible folding pathways. Such an algorithm needs several modifications before it can predict RNA secondary structures. After modification, a fair number of correct stems are predicted, even when using computationally quick, but very crude, fitness criteria such as stem length and stacking energy, including elements of tertiary structure (pseudoknots). The fact that genetic algorithm simulation includes both stem formations and stem disruption allows one to observe intermediate structures that may be used in combination with phylogenetic or experimental research.