Coupled nucleotide covariations reveal dynamic RNA interaction patterns

Profiling the intragenic toxicity determinants of toxin-antitoxin systems: revisiting hok/Sok regulation

Type I toxin-antitoxin systems (T1TAs) are extremely potent bacterial killing systems difficult to characterize using classical approaches. To assess the killing capability of type I toxins and to identify mutations suppressing the toxin expression or activity, we previously developed the FASTBAC-Seq (Functional AnalysiS of Toxin-Antitoxin Systems in BACteria by Deep Sequencing) method in Helicobacter pylori. This method combines a life and death selection with deep sequencing. Here, we adapted and improved our method to investigate T1TAs in the model organism Escherichia coli. As a proof of concept, we revisited the regulation of the plasmidic hok/Sok T1TA system. We revealed the death-inducing phenotype of the Hok toxin when it is expressed from the chromosome in the absence of the antitoxin and recovered previously described intragenic toxicity determinants of this system. We identified nucleotides that are essential for the transcription, translation or activity of Hok. We also discovered single-nucleotide substitutions leading to structural changes affecting either the translation or the stability of the hok mRNA. Overall, we provide the community with an easy-to-use approach to widely characterize TA systems from diverse types and bacteria.

Variation Profile of the Orthotospovirus Genome

Orthotospoviruses are plant-infecting members of the family Tospoviridae (order Bunyavirales), have a broad host range and are vectored by polyphagous thrips in a circulative-propagative manner. Because diverse hosts and vectors impose heterogeneous selection constraints on viral genomes, the evolutionary arms races between hosts and their pathogens might be manifested as selection for rapid changes in key genes. These observations suggest that orthotospoviruses contain key genetic components that rapidly mutate to mediate host adaptation and vector transmission. Using complete genome sequences, we profiled genomic variation in orthotospoviruses. Results show that the three genomic segments contain hypervariable areas at homologous locations across species. Remarkably, the highest nucleotide variation mapped to the intergenic region of RNA segments S and M, which fold into a hairpin. Secondary structure analyses showed that the hairpin is a dynamic structure with multiple functional shapes formed by stems and loops, contains sites under positive selection and covariable sites. Accumulation and tolerance of mutations in the intergenic region is a general feature of orthotospoviruses and might mediate adaptation to host plants and insect vectors.

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.

Genetic variability and evolutionary dynamics of viruses of the family Closteroviridae

RNA viruses have a great potential for genetic variation, rapid evolution and adaptation. Characterization of the genetic variation of viral populations provides relevant information on the processes involved in virus evolution and epidemiology and it is crucial for designing reliable diagnostic tools and developing efficient and durable disease control strategies. Here we performed an updated analysis of sequences available in Genbank and reviewed present knowledge on the genetic variability and evolutionary processes of viruses of the family Closteroviridae. Several factors have shaped the genetic structure and diversity of closteroviruses. (I) A strong negative selection seems to be responsible for the high genetic stability in space and time for some viruses. (2) Long distance migration, probably by human transport of infected propagative plant material, have caused that genetically similar virus isolates are found in distant geographical regions. (3) Recombination between divergent sequence variants have generated new genotypes and plays an important role for the evolution of some viruses of the family Closteroviridae. (4) Interaction between virus strains or between different viruses in mixed infections may alter accumulation of certain strains. (5) Host change or virus transmission by insect vectors induced changes in the viral population structure due to positive selection of sequence variants with higher fitness for host-virus or vector-virus interaction (adaptation) or by genetic drift due to random selection of sequence variants during the population bottleneck associated to the transmission process.

Structural basis for the biological relevance of the invariant apical stem in IRES-mediated translation

RNA structure plays a fundamental role in internal initiation of translation. Picornavirus internal ribosome entry site (IRES) are long, efficient cis-acting elements that recruit the ribosome to internal mRNA sites. However, little is known about long-range constraints determining the IRES RNA structure. Here, we sought to investigate the functional and structural relevance of the invariant apical stem of a picornavirus IRES. Mutation of this apical stem revealed better performance of G:C compared with C:G base pairs, demonstrating that the secondary structure solely is not sufficient for IRES function. In turn, mutations designed to disrupt the stem abolished IRES activity. Lack of tolerance to accept genetic variability in the apical stem was supported by the presence of coupled covariations within the adjacent stem-loops. SHAPE structural analysis, gel mobility-shift and microarrays-based RNA accessibility revealed that the apical stem contributes to maintain IRES RNA structure through the generation of distant interactions between two adjacent stem-loops. Our results demonstrate that a highly interactive structure constrained by distant interactions involving invariant G:C base pairs plays a key role in maintaining the RNA conformation necessary for IRES-mediated translation.

RNA antitoxins

Recent genomic analyses revealed a surprisingly large number of toxin-antitoxin loci in free-living prokaryotes. The antitoxins are proteins or antisense RNAs that counteract the toxins. Two antisense RNA-regulated toxin-antitoxin gene families, hok/sok and ldr, are unrelated sequence-wise but have strikingly similar properties at the level of gene and RNA organization. Recently, two SOS-induced toxins were found to be regulated by RNA antitoxins. One such toxin, SymE, exhibits similarity with MazE antitoxin and, surprisingly, inhibits translation. Thus, it is possible that an ancestral antitoxin gene evolved into the present toxin gene (symE) whose translation is repressed by an RNA antitoxin (SymR).

Competitive inhibition of natural antisense Sok-RNA interactions activates Hok-mediated cell killing in Escherichia coli

Short regulatory RNAs are widespread in bacteria, and many function through antisense recognition of mRNA. Among the best studied antisense transcripts are RNA antitoxins that repress toxin mRNA translation. The hok/sok locus of plasmid R1 from Escherichia coli is an established model for RNA antitoxin action. Base-pairing between hok mRNA and Sok-antisense-RNA increases plasmid maintenance through post-segregational-killing of plasmid-free progeny cells. To test the model and the idea that sequestration of Sok-RNA activity could provide a novel antimicrobial strategy, we designed anti Sok peptide nucleic acid (PNA) oligomers that, according to the model, would act as competitive inhibitors of hok mRNA::Sok-RNA interactions. In hok/sok-carrying cells, anti Sok PNAs were more bactericidal than rifampicin. Also, anti Sok PNAs induced ghost cell morphology and an accumulation of mature hok mRNA, consistent with cell killing through synthesis of Hok protein. The results support the sense/antisense model for hok mRNA repression by Sok-RNA and demonstrate that antisense agents can be used to out-compete RNA::RNA interactions in bacteria. Finally, BLAST analyses of approximately 200 prokaryotic genomes revealed that many enteric bacteria have multiple hok/sok homologous and analogous RNA-regulated toxin-antitoxin loci. Therefore, it is possible to activate suicide in bacteria by targeting antitoxins.

Free energy landscapes of RNA/RNA complexes: with applications to snRNA complexes in spliceosomes

We develop a statistical mechanical model for RNA/RNA complexes with both intramolecular and intermolecular interactions. As an application of the model, we compute the free energy landscapes, which give the full distribution for all the possible conformations, for U4/U6 and U2/U6 in major spliceosome and U4atac/U6atac and U12/U6atac in minor spliceosome. Different snRNA experiments found contrasting structures, our free energy landscape theory shows why these structures emerge and how they compete with each other. For yeast U2/U6, the model predicts that the two distinct experimental structures, the four-helix junction structure and the helix Ib-containing structure, can actually coexist and specifically compete with each other. In addition, the energy landscapes suggest possible mechanisms for the conformational switches in splicing. For instance, our calculation shows that coaxial stacking is essential for stabilizing the four-helix junction in yeast U2/U6. Therefore, inhibition of the coaxial stacking possibly by protein-binding may activate the conformational switch from the four-helix junction to the helix Ib-containing structure. Moreover, the change of the energy landscape shape gives information about the conformational changes. We find multiple (native-like and misfolded) intermediates formed through base-pairing rearrangements in snRNA complexes. For example, the unfolding of the U2/U6 undergoes a transition to a misfolded state which is functional, while in the unfolding of U12/U6atac, the functional helix Ib is found to be the last one to unfold and is thus the most stable structural component. Furthermore, the energy landscape gives the stabilities of all the possible (functional) intermediates and such information is directly related to splicing efficiency.

Evidence for a base triple in the free HIV-1 TAR RNA

We propose the existence of a novel base triple in the HIV-1 TAR hairpin. This triple is supported by covariation of loop residue 31 with residue 22, which is part of an unusual base pair with U40 below the 3-nucleotide bulge. A set of mutants was constructed to test the involvement of bases A22, U31, and U40 in a triple interaction. RNA structure probing, trans-activation assays, and structure modeling are consistent with the existence of this base triple in a bent conformation of the free TAR element. However, disruption of the base triple does not affect binding of a Tat-derived peptide. We therefore compared the structure of free and Tat-bound TAR RNA by footprinting and site-specific cross-linking analyses. These studies indicate that the Tat arginine-rich motif, in addition to its known binding site at the bulge, is in close contact with U31 in the TAR loop. Because binding of Tat to TAR is known to coincide with the formation of a base triple with residues U23, A27, and U38, we hypothesize that Tat binding and the associated straightening of TAR triggers the disruption of the (A22-U40)U31 triple.

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.