Solution structure and functional importance of a conserved RNA hairpin of eel LINE UnaL2

Acylation probing of "generic" RNA libraries reveals critical influence of loop constraints on reactivity

The reactivity of RNA 2'-OH acylation is broadly useful both in probing structure and in preparing conjugates. To date, this reactivity has been analyzed in limited sets of biological RNA sequences, leaving open questions of how reactivity varies inherently without regard to sequence in structured contexts. We constructed and probed "generic" structured RNA libraries using homogeneous loop sequences, employing deep sequencing to carry out a systematic survey of reactivity. We find a wide range of RNA reactivities among single-stranded sequences, with nearest neighbors playing substantial roles. Remarkably, certain small loops are found to be far more reactive on average (up to 4,000-fold) than single-stranded RNAs, due to conformational constraints that enhance reactivity. Among loops, we observe large variations in reactivity based on size, type, and position. The results lend insights into RNA designs for achieving high-efficiency local conjugation and provide new opportunities to refine structure analysis.

Discovery of alternatively spliced isoforms and long non-coding RNA in full length brain transcriptomes of anadromous Hilsa shad, Tenualosa ilisha (Hamilton, 1822)

BACKGROUND

Full length transcriptomes, achieved through long-read sequencing, along with the isoform analysis can reveal complexities in the gene expression profiles, as well as annotate the transcriptomes of non-model organisms.

METHODS AND RESULT

Full length transcripts of brain transcriptome of Tenualosa ilisha, Hilsa shad, were generated through PacBio single molecule real-time sequencing and were characterized. A total of 8.30 Gb clean reads were generated, with PacBio RSII, which resulted in 57,651 high quality consensus transcripts. After removing redundant reads, a total of 19,220 high-quality non-redundant transcripts and 17,341 full length ORF transcripts were classified to 7522 putative ortholog groups. Genes involved in various neural pathways were identified. In addition, isoform clusters and lncRNAs were discovered, along with Hilsa specific transcripts with coding frames and 29,147 SSRs in 944 transcripts (1141 annotated).

CONCLUSION

The present study provided, for the first time, a comprehensive view of the alternative isoforms of genes and transcriptome complexity in Hilsa shad brain and forms a rich resource for functional studies in brain of this anadromous fish.

Recognition of 3'-end L1, Alu, processed pseudogenes, and mRNA stem-loops in the human genome using sequence-based and structure-based machine-learning models

The role of 3′-end stem-loops in retrotransposition was experimentally demonstrated for transposons of various species, where LINE-SINE retrotransposons share the same 3′-end sequences, containing a stem-loop. We have discovered that 62–68% of processed pseduogenes and mRNAs also have 3′-end stem-loops. We investigated the properties of 3′-end stem-loops of human L1s, Alus, processed pseudogenes and mRNAs that do not share the same sequences, but all have 3′-end stem-loops. We have built sequence-based and structure-based machine-learning models that are able to recognize 3′-end L1, Alu, processed pseudogene and mRNA stem-loops with high performance. The sequence-based models use only sequence information and capture compositional bias in 3′-ends. The structure-based models consider physical, chemical and geometrical properties of dinucleotides composing a stem and position-specific nucleotide content of a loop and a bulge. The most important parameters include shift, tilt, rise, and hydrophilicity. The obtained results clearly point to the existence of structural constrains for 3′-end stem-loops of L1 and Alu, which are probably important for transposition, and reveal the potential of mRNAs to be recognized by the L1 machinery. The proposed approach is applicable to a broader task of recognizing RNA (DNA) secondary structures. The constructed models are freely available at github (https://github.com/AlexShein/transposons/).

Distinct RNA recognition mechanisms in closely related LINEs from zebrafish

Long interspersed nuclear element (LINE) is known to be transposed by the reverse transcription using its RNA transcript. Recognition of the 3' stem-loop of LINE RNA by its reverse transcriptase (RT) is an important step of the retrotransposition. We focused on the RNA recognition by RT from two related LINEs, ZfL2-1 and ZfL2-2, from zebrafish. Previous study showed that RT from ZfL2-2 recognizes a single residue in the specific position of the RNA loop. In the present study, it was found that RT from ZfL2-1 recognizes the inserted stem-loop of ZfL2-1 RNA. Thus, these related RTs recognize the same region of LINE RNAs but discriminate them by different mechanism.

Solution structure of a reverse transcriptase recognition site of a LINE RNA from zebrafish

Long interspersed nuclear element (LINE) is known to be transposed by reverse transcription using its RNA transcript. Recognition of the 3' stem-loop of LINE RNA by its reverse transcriptase (RT) is an important step of the retrotransposition. Our previous study revealed that the second G residue (G8) in the GGAUA loop of a 17mer LINE RNA from eel, UnaL2-17, is recognized by its RT and the U residue (U10) in the same loop is required to maintain the loop structure (Baba S, Kajikawa M, Okada N, Kawai G. Solution structure of an RNA stem-loop derived from the 3' conserved region of eel LINE UnaL2. RNA 2004;10:1380-1387). ZfL2-2, a LINE from zebrafish, has the same 3' stem-loop with UnaL2 and ZfL2-1 has similar but distinct 3' stem-loop with an insertion which can form an additional stem-loop. Here, we determined the solution structure of the 34mer RT recognition site of the LINE RNA (ZfL2-1-34). It was found that ZfL2-1-34 forms a hairpin with an internal loop, the tertiary structure of which is superimposed with that of ZfL2-2. It is noted that A10 and the inserted stem-loop, starting with A12, in ZfL2-1-34 located at the positions corresponding to those of G8 and U10, respectively, in UnaL2-17. These results strongly suggest that the two LINEs share the similar recognition mechanism and the A10 in ZfL2-1-34 is the determinant recognized by its RT.

Conserved 3' UTR stem-loop structure in L1 and Alu transposons in human genome: possible role in retrotransposition

BACKGROUND

In the process of retrotransposition LINEs use their own machinery for copying and inserting themselves into new genomic locations, while SINEs are parasitic and require the machinery of LINEs. The exact mechanism of how a LINE-encoded reverse transcriptase (RT) recognizes its own and SINE RNA remains unclear. However it was shown for the stringent-type LINEs that recognition of a stem-loop at the 3'UTR by RT is essential for retrotransposition. For the relaxed-type LINEs it is believed that the poly-A tail is a common recognition element between LINE and SINE RNA. However polyadenylation is a property of any messenger RNA, and how the LINE RT recognizes transposon and non-transposon RNAs remains an open question. It is likely that RNA secondary structures play an important role in RNA recognition by LINE encoded proteins.

RESULTS

Here we selected a set of L1 and Alu elements from the human genome and investigated their sequences for the presence of position-specific stem-loop structures. We found highly conserved stem-loop positions at the 3'UTR. Comparative structural analyses of a human L1 3'UTR stem-loop showed a similarity to 3'UTR stem-loops of the stringent-type LINEs, which were experimentally shown to be recognized by LINE RT. The consensus stem-loop structure consists of 5-7 bp loop, 8-10 bp stem with a bulge at a distance of 4-6 bp from the loop. The results show that a stem loop with a bulge exists at the 3'-end of Alu. We also found conserved stem-loop positions at 5'UTR and at the end of ORF2 and discuss their possible role.

CONCLUSIONS

Here we presented an evidence for the presence of a highly conserved 3'UTR stem-loop structure in L1 and Alu retrotransposons in the human genome. Both stem-loops show structural similarity to the stem-loops of the stringent-type LINEs experimentally confirmed as essential for retrotransposition. Here we hypothesize that both L1 and Alu RNA are recognized by L1 RT via the 3'-end RNA stem-loop structure. Other conserved stem-loop positions in L1 suggest their possible functions in protein-RNA interactions but to date no experimental evidence has been reported.

Mechanism by which a LINE protein recognizes its 3' tail RNA

LINEs mobilize their own copies via retrotransposition. LINEs can be divided into two types. One is a stringent type, which constitutes a majority of LINEs. The other is a relaxed type. To elucidate the molecular mechanism of retrotransposition, we used here two different zebrafish LINEs belonging to the stringent type. By using retrotransposition assays, we demonstrated that proteins (ORF2) encoded by an individual LINE recognize the cognate 3′ tail sequence of the LINE RNA strictly. By conducting in vitro binding assays with a variety of ORF2 proteins, we demonstrated that the region between the endonuclease and reverse transcriptase domains in ORF2 is the site at which the proteins bind the stem-loop structure of the 3′ tail RNA, showing that the strict recognition of the stem-loop structure by the cognate ORF2 protein is an important step in retrotransposition. This recognition can be bipartite, involving the general recognition of the stem by cTBR (conserved tail-binding region) of ORF2 and the specific recognition of the loop by vTBR (variable tail-binding region). This is the first report that clearly characterized the RNA-binding region in ORF2, providing the generality for the recognition mechanism of the RNA tail by the ORF2 protein encoded by LINEs.

Solution structure of a DNA mimicking motif of an RNA aptamer against transcription factor AML1 Runt domain

AML1/RUNX1 is an essential transcription factor involved in the differentiation of hematopoietic cells. AML1 binds to the Runt-binding double-stranded DNA element (RDE) of target genes through its N-terminal Runt domain. In a previous study, we obtained RNA aptamers against the AML1 Runt domain by systematic evolution of ligands by exponential enrichment and revealed that RNA aptamers exhibit higher affinity for the Runt domain than that for RDE and possess the 5'-GCGMGNN-3' and 5'-N'N'CCAC-3' conserved motif (M: A or C; N and N' form Watson-Crick base pairs) that is important for Runt domain binding. In this study, to understand the structural basis of recognition of the Runt domain by the aptamer motif, the solution structure of a 22-mer RNA was determined using nuclear magnetic resonance. The motif contains the AH(+)-C mismatch and base triple and adopts an unusual backbone structure. Structural analysis of the aptamer motif indicated that the aptamer binds to the Runt domain by mimicking the RDE sequence and structure. Our data should enhance the understanding of the structural basis of DNA mimicry by RNA molecules.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

The RNA-binding region of human TRBP interacts with microRNA precursors through two independent domains

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression through RNA interference. Human miRNAs are generated through a series of enzymatic processing steps. The precursor miRNA (pre-miRNA) is recognized and cleaved by a complex containing Dicer and several non-catalytic accessory proteins. HIV TAR element binding protein (TRBP) is a constituent of the Dicer complex, which augments complex stability and potentially functions in substrate recognition and product transfer to the RNA-induced silencing complex. Here we have analysed the interaction between the RNA-binding region of TRBP and an oncogenic human miRNA, miR-155, at different stages in the biogenesis pathway. We show that the region of TRBP that binds immature miRNAs comprises two independent double-stranded RNA-binding domains connected by a 60-residue flexible linker. No evidence of contact between the two double-stranded RNA-binding domains was observed either in the apo- or RNA-bound state. We establish that the RNA-binding region of TRBP interacts with both pre-miR-155 and the miR-155/miR-155* duplex through the same binding surfaces and with similar affinities, and that two protein molecules can simultaneously interact with each immature miRNA. These data suggest that TRBP could play a role before and after processing of pre-miRNAs by Dicer.

Hydrotalcite Intercalated siRNA: Computational Characterization of the Interlayer Environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.

Parallel relaxation of stringent RNA recognition in plant and mammalian L1 retrotransposons

L1 elements are mammalian non-long terminal repeat retrotransposons, or long interspersed elements (LINEs), that significantly influence the dynamics and fluidity of the genome. A series of observations suggest that plant L1-clade LINEs, just as mammalian L1s, mobilize both short interspersed elements (SINEs) and certain messenger RNA by recognizing the 3'-poly(A) tail of RNA. However, one L1 lineage in monocots was shown to possess a conserved 3'-end sequence with a solid RNA structure also observed in maize and sorghum SINEs. This strongly suggests that plant LINEs require a particular 3'-end sequence during initiation of reverse transcription. As one L1-clade LINE was also found to share the 3'-end sequence with a SINE in a green algal genome, I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized the specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Self-interaction, nucleic acid binding, and nucleic acid chaperone activities are unexpectedly retained in the unique ORF1p of zebrafish LINE

Long interspersed elements (LINEs) are mobile elements that comprise a large proportion of many eukaryotic genomes. Although some LINE-encoded open reading frame 1 proteins (ORF1ps) were suggested to be required for LINE mobilization through binding to their RNA, their general role is not known. The ZfL2-1 ORF1p, which belongs to the esterase-type ORF1p, is especially interesting because it has no known RNA-binding domain. Here we demonstrate that ZfL2-1 ORF1p has all the canonical activities associated with known ORF1ps, including self-interaction, nucleic acid binding, and nucleic acid chaperone activities. In particular, we showed that its chaperone activity is reversible, suggesting that the chaperone activities of many other ORF1ps are also reversible. From this discovery, we propose that LINE ORF1ps play a general role in LINE integration by forming a complex with LINE RNA and rearranging its conformation.

Genetically tagged TRE5-A retrotransposons reveal high amplification rates and authentic target site preference in the Dictyostelium discoideum genome

Retrotransposons contribute significantly to the evolution of eukaryotic genomes. They replicate by producing DNA copies of their own RNA, which are integrated at new locations in the host cell genome. In the gene-dense genome of the social amoeba Dictyostelium discoideum, retrotransposon TRE5-A avoids insertional mutagenesis by targeting the transcription factor (TF) IIIC/IIIB complex and integrating ∼ 50 bp upstream of tRNA genes. We generated synthetic TRE5-A retrotransposons (TRE5-A(bsr)) that were tagged with a selection marker that conferred resistance to blasticidin after a complete retrotransposition cycle. We found that the TRE5-A(bsr) elements were efficiently mobilized in trans by proteins expressed from the endogenous TRE5-A population found in D. discoideum cells. ORF1 protein translated from TRE5-A(bsr) elements significantly enhanced retrotransposition. We observed that the 5' untranslated region of TRE5-A could be replaced by an unrelated promoter, whereas the 3' untranslated region of TRE5-A was essential for retrotransposition. A predicted secondary structure in the RNA of the 3' untranslated region of TRE5-A may be involved in the retrotransposition process. The TRE5-A(bsr) elements were capable of identifying authentic integration targets in vivo, including formerly unnoticed, putative binding sites for TFIIIC on the extrachromosomal DNA element that carries the ribosomal RNA genes.

Physics-based de novo prediction of RNA 3D structures

Current experiments on structural determination cannot keep up the pace with the steadily emerging RNA sequences and new functions. This underscores the request for an accurate model for RNA three-dimensional (3D) structural prediction. Although considerable progress has been made in mechanistic studies, accurate prediction for RNA tertiary folding from sequence remains an unsolved problem. The first and most important requirement for the prediction of RNA structure from physical principles is an accurate free energy model. A recently developed three-vector virtual bond-based RNA folding model ("Vfold") has allowed us to compute the chain entropy and predict folding free energies and structures for RNA secondary structures and simple pseudoknots. Here we develop a free energy-based method to predict larger more complex RNA tertiary folds. The approach is based on a multiscaling strategy: from the nucleotide sequence, we predict the two-dimensional (2D) structures (defined by the base pairs and tertiary contacts); based on the 2D structure, we construct a 3D scaffold; with the 3D scaffold as the initial state, we combine AMBER energy minimization and PDB-based fragment search to predict the all-atom structure. A key advantage of the approach is the statistical mechanical calculation for the conformational entropy of RNA structures, including those with cross-linked loops. Benchmark tests show that the model leads to significant improvements in RNA 3D structure prediction.

Flexibility of Short-Strand RNA in Aqueous Solution as Revealed by Molecular Dynamics Simulation: Are A-RNA and A´-RNA Distinct Conformational Structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A′-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA – namely major groove width, inclination and the number of base pairs in a helical twist – over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A′-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors).

Analysis and classification of RNA tertiary structures

There is a fast growing interest in noncoding RNA transcripts. These transcripts are not translated into proteins, but play essential roles in many cellular and pathological processes. Recent efforts toward comprehension of their function has led to a substantial increase in both the number and the size of solved RNA structures. With the aim of addressing questions relating to RNA structural diversity, we examined RNA conservation at three structural levels: primary, secondary, and tertiary structure. Additionally, we developed an automated method for classifying RNA structures based on spatial (three-dimensional [3D]) similarity. Applying the method to all solved RNA structures resulted in a classified database of RNA tertiary structures (DARTS). DARTS embodies 1333 solved RNA structures classified into 94 clusters. The classification is hierarchical, reflecting the structural relationship between and within clusters. We also developed an application for searching DARTS with a new structure. The search is fast and its performance was successfully tested on all solved RNA structures since the creation of DARTS. A user-friendly interface for both the database and the search application is available online. We show intracluster and intercluster similarities in DARTS and demonstrate the usefulness of the search application. The analysis reveals the current structural repertoire of RNA and exposes common global folds and local tertiary motifs. Further study of these conserved substructures may suggest possible RNA domains and building blocks. This should be beneficial for structure prediction and for gaining insights into structure-function relationships.