RNA pseudoknots downstream of the frameshift sites of retroviruses.

RNA pseudoknot structural motifs could have implications for a wide range of biological processes of RNAs. In this study, the potential RNA pseudoknots just downstream from the known and suspected retroviral frame-shift sites were predicted in the Rous sarcoma virus, primate immunodeficiency viruses (HIV-1, HIV-2, and SIV), equine infectious anemia virus, visna virus, bovine leukemia virus, human T-cell leukemia virus (types I and II), mouse mammary tumor virus, Mason-Pfizer monkey virus, and simian SRV-1 type-D retrovirus. Also, the putative RNA pseudoknots were detected in the gag-pol overlaps of two retrotransposons of Drosophila, 17.6 and gypsy, and the mouse intracisternal A particle. For each sequence, the thermodynamic stability and statistical significance of the secondary structure involved in the predicted tertiary structure were assessed and compared. Our results show that the stem-loop structures in the pseudoknots are both thermodynamically highly stable and statistically significant relative to other such configurations that potentially occur in the gag-pol or gag-pro and pro-pol junction domains of these viruses (300 nucleotides upstream and downstream from the possible frameshift sites are included). Moreover, the structural features of the predicted pseudoknots following the frameshift site of pro-pol overlaps of the HTLV-1 and HTLV-2 retroviruses are structurally well conserved. The occurrence of eight compensatory base changes in the tertiary interaction of the two related sequences allow the conservation of their tertiary structures in spite of the sequence divergence. The results support the possible control mechanism for frameshifting proposed by Brierley et al. and Jacks et al.     

An algorithm for searching RNA motifs in genomic sequences

RNA molecules, which are found in all living cells, fold into characteristic structures that account for their diverse functional activities. Many of these RNA structures consist of a collection of fundamental RNA motifs. The various combinations of RNA basic components form different RNA classes and define their unique structural and functional properties. The availability of many genome sequences makes it possible to search computationally for functional RNAs. Biological experiments indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop regions. The searching for those well-ordered RNA structures and their homologues in genomic sequences is very helpful for the understanding of RNA-based gene regulation. In this paper, we consider the following problem: given an RNA sequence with a known secondary structure, efficiently determine candidate segments in genomic sequences that can potentially form RNA secondary structures similar to the given RNA secondary structure. Our new bottom-up approach searches all potential stem-loops similar to ones of the given RNA secondary structure first, and then based on located stem-loops, detects potential homologous structural RNAs in genomic sequences.     

Stability of RNA stem-loop structure and distribution of non-random structure in the human immunodeficiency virus (HIV-I).

The stability of potential RNA stem-loop structures in human immunodeficiency virus isolates, HTLV-III and ARV, has been calculated, and the relevance to the local significant secondary structures in the sequence has been tested statistically using a Monte Carlo simulation method. Potentially significant structures exist in the 5'non-coding region, the boundary regions between the protein coding frames, and the 3' non-coding region. The locally optimal secondary structure occurring in the 5' terminal region has been assessed using different overlapping segment sizes and the Monte Carlo method. The results show that the most favorable structure for the 5' mRNA leader sequence of HIV has two stem-loops folded at nucleotides 5-104 in the R region (stem-loop I, 5-54 and stem-loop II, 58-104). A large fluctuation of segment score of the local optimal secondary structure also occurs in the boundary between the exterior glycosylated protein or outer membrane protein and transmembrane protein coding region. This finding is surprising since no RNA signals or RNA processing are expected to occur at this site. In addition, regions of the genome predicted to have significantly more open structure at the RNA level correlate closely with hypervariable sites found in these viral genomes. The possible importance of local secondary structure to the biological function of the human immunodeficiency virus genome is discussed.     

Human T-cell leukemia virus type II Rex binding and activity require an intact splice donor site and a specific RNA secondary structure.

The human T-cell leukemia virus type II (HTLV-II) regulatory protein Rex augments cytoplasmic levels of unspliced gag-pol mRNA by acting through a Rex-responsive element (RxRE) in the long terminal repeat. Purified Rex protein binds to long terminal repeat mRNA. Here, using an immunobinding assay to measure the binding of Rex protein to mutated RxRE RNAs, we show that efficient Rex binding requires a stem-bulge-loop RNA secondary structure (nucleotides [nt] 465 to 500) and specific sequences both within the stem-bulge (nt 470 to 476) and within a conserved upstream splice donor site (nt 449 to 455). Rex function in a transient transfection expression system correlates with Rex protein-RxRE RNA binding. The ability of HTLV-II Rex to interact directly with the HTLV-II splice donor site suggests that HTLV-II Rex may increase expression of unspliced gag-pol mRNA, in part, by inhibiting splicing.     

An RNA secondary structure juxtaposes two remote genetic signals for human T-cell leukemia virus type I RNA 3'-end processing.

Sequence analysis of the human T-cell leukemia virus type I (HTLV-I) long terminal repeat (LTR) does not reveal a polyadenylation consensus sequence, AAUAAA, close to the polyadenylation site at the 3' end of the viral RNA. Using site-directed mutagenesis, we demonstrated that two cis-acting signals are required for efficient RNA processing in HTLV-I LTR: (i) a remote AAUAAA hexamer at a distance of 276 nucleotides upstream of the polyadenylation site, and (ii) the 20-nucleotide GU-rich sequence immediately downstream from the poly(A) site. It has been postulated that the folding of RNA into a secondary structure juxtaposes the AAUAAA sequence, in a noncontiguous manner, to within 14 nucleotides of the polyadenylation site. To test this hypothesis, we introduced deletions and point mutations within the U3 and R regions of the LTR. RNA 3'-end processing occurred efficiently at the authentic HTLV-I poly(A) site after deletion of the sequences predicted to form the secondary structure. Thus, the genetic analysis supports the hypothesis that folding of the HTLV-I RNA in the U3 and R regions juxtaposes the AAUAAA sequence and the poly(A) site to the correct functional distance. This unique arrangement of RNA-processing signals is also found in the related retroviruses HTLV-II and bovine leukemia virus.     

稀有的黑檀体基因(ebony)5'UTR自发突变导致黑腹果蝇的黑条体突变

由于果蝇Drosophila群体中有很多自发突变其中包括多种体色突变,因此它是一个研究自发突变的优秀的模式体系。本研究证实我们实验室发现的一个可以引起果蝇体色突变的自发突变(bsr)是一个黑檀体(e)的等位基因,将其命名为ebsr。序列分析显示ebsr的5′端缺失了953个碱基,其中包括外显子1后端的206个碱基及相连的内含子1的747个碱基。逆转录PCR结果显示5′端的缺失导致内含子1不能从mRNA中剪接掉,由此导致该mRNA的翻译起始密码子AUG前端增加了一个3.2kb的序列。该序列导致ebsr的mRNA的5′UTR(5′-untranslated region)区较野生型基因增加近3kb的长度。通过mRNA二级结构分析发现这个增加的3kb的片段可以形成复杂的颈环结构(stem-loop)。免疫印迹结果显示该突变基因没有基因产物产生。本研究进一步证实了由于mRNA的5′UTR序列结构的改变可以影响到蛋白质的翻译。     

Thermodynamic stability and statistical significance of potential stem-loop structures situated at the frameshift sites of retroviruses.

RNA stem-loop structures situated just 3' to the frameshift sites of the retroviral gag-pol or gag-pro and pro-pol regions may make important contributions to frame-shifting in retroviruses. In this study, the thermodynamic stability and statistical significance of such secondary structural features relative to others in the sequence have been assessed using a newly developed method that combines calculations of the lowest free energy of formation of RNA secondary structures and the Monte Carlo simulations. Our results show that stem-loop structures situated just 3' to the frameshift sites are both highly stable and statistically significant relative to others in the gag-pol or gag-pro and pro-pol junction domains (both 300 nucleotides upstream and downstream from the possible frameshift sites are included) of Rous sarcoma virus (RSV), human immunodeficiency virus (HIV-1), bovine leukemia virus (BLV), human T-cell leukemia virus type II (HTLV-II), and mouse mammary tumor virus (MMTV). No other more stable, or significant folding regions are predicted in these domains.     

More than half of yeast U1 snRNA is dispensable for growth.

Yeast U1 snRNA (568 nucleotides) is 3.5-fold larger than its mammalian counterpart (164 nucleotides) and contains apparent sequence homology only at the 5' and 3' ends. We have used deletion analysis to determine whether the yeast-specific U1 sequences play essential roles in vivo. Yeast cells carrying a deletion of more than 60% (355 nucleotides) of the single-copy U1 gene are viable, though slow-growing, while a deletion of 316 nucleotides allows essentially wild-type growth. The boundaries of the viable deletions define a dispensable internal domain which comprises sequences unique to yeast. In contrast, the essential 5' and 3' terminal domains correspond to phylogenetically conserved sequences and/or structures previously implicated in RNA:RNA and RNA:protein interactions. The minimal essential sequences of yeast U1 can be drawn in a secondary structure which resembles metazoan U1 in four of seven structural domains.     

A highly conserved RNA folding region coincident with the Rev response element of primate immunodeficiency viruses.

A series of unusual folding regions (UFR) immediately 3' to the cleavage site of the outer membrane protein (OMP) and transmembrane protein (TMP) were detected in the envelope gene RNA of the human immunodeficiency virus (HIV-1, HIV-2) and simian immunodeficiency virus (SIV) by an extensive Monte Carlo simulation. These RNA secondary structures were predicted to be both highly stable and statistically significant. In the calculation, twenty-five different sequence isolates of HIV-1, three isolates of HIV-2 and eight sequences of SIV were included. Although significant sequence divergence occurs in the env coding regions of these viruses, a distinct UFR of 234-nt is consistently located ten nucleotides 3' to the cleavage site of the OMP/TMP in HIV-1, and a 216-nt UFR occurs forty-six and forty-nine nucleotides downstream from the OMP/TMP cleavage site of HIV-2 and SIV, respectively. Compensatory base changes in the helical stem regions of these conserved RNA secondary structures are identified. These results support the hypothesis that these special RNA folding regions are functionally important and suggest that the role of this sequence as the Rev response element (RRE) is mediated by secondary structure as well as primary RNA sequence.     

Multimer Formation Is Not Essential for Nuclear Export of Human T-Cell Leukemia Virus Type 1 Rex trans-Activator Protein

The Rex trans-regulatory protein of human T-cell leukemia virus type 1 (HTLV-1) is required for the nuclear export of incompletely spliced and unspliced viral mRNAs and is therefore essential for virus replication. Rex is a nuclear phosphoprotein that directly binds to its cis-acting Rex response element RNA target sequence and constantly shuttles between the nucleus and cytoplasm. Moreover, Rex induces nuclear accumulation of unspliced viral RNA. Three protein domains which mediate nuclear import-RNA binding, nuclear export, and Rex oligomerization have been mapped within the 189-amino-acid Rex polypeptide. Here we identified a different region in the carboxy-terminal half of Rex which is also required for biological activity. In inactive mutants with mutations that map within this region, as well as in mutants that are deficient in Rex-specific multimerization, Rex trans activation could be reconstituted by fusion to a heterologous leucine zipper dimerization interface. The intracellular trafficking capabilities of wild-type and mutant Rex proteins reveal that biologically inactive and multimerization-deficient Rex mutants are still efficiently translocated from the nucleus to the cytoplasm. This observation indicates that multimerization is essential for Rex function but is not required for nuclear export. Finally, we are able to provide an improved model of the HTLV-1 Rex domain structure.     

Complete nucleotide sequence of the new simian T-lymphotropic virus, STLV-PH969 from a Hamadryas baboon, and unusual features of its long terminal repeat.

A third type of primate T-lymphotropic virus, PTLV-L, with STLV-PH969 as a prototype, has recently been isolated from an African baboon (Papio hamadryas). Classification of this virus has been based on partial sequence analysis of cDNA from a virus-producing cell line, PH969. We obtained the complete nucleotide sequence of this virus with a proviral genome of 8,916 bp. All major genes, homologous in all human T-cell lymphotropic virus (HTLV)-related viruses, and their corresponding mRNAs, including appropriate splicing, were identified. One additional nonhomologous open reading frame in the proximal pX region is accessible for translation through alternative splicing. Sequence comparison shows that STLV-PH969 is equidistantly related to HTLV type 1 (HTLV-1) and HTLV-2. In all coding regions, the similarity tends to be the lowest between STLV-PH969 and HTLV-1. However, in the long terminal repeat (LTR) region, the lowest similarity was found between STLV-PH969 and HTLV-2. The U3-R and R-U5 boundaries of the STLV-PH969 LTR were experimentally determined at nucleotides 268 and 524, respectively. This 695-bp LTR is 60 and 73 bp shorter than the LTRs of HTLV-1 and HTLV-2, respectively, but its general organization is similar to the one found in the HTLV-bovine leukemia virus genus. In the long region between the polyadenylation signal and the poly(A) site, sequence similarity with the HTLV-1 Rex-responsive element (RexRE) core and secondary structure prediction suggest the presence of a RexRE. The presence of three 21-bp repeats is conserved within the U3 region of HTLV-1, HTLV-2, and BLV. Only two direct repeats with similarity to these Tax-responsive elements were found in the STLV-PH969 LTR, which might suggest differences in the Tax-mediated transactivation of this virus. We conclude that STLV-PH969 has all the genes and genomic regions to suggest a replication cycle comparable to that of HTLV-1 and HTLV-2.     

Thermodynamic and phylogenetic prediction of RNA secondary structures in the coding region of hepatitis C virus

The existence and functional importance of RNA secondary structure in the replication of positive-stranded RNA viruses is increasingly recognized. We applied several computational methods to detect RNA secondary structure in the coding region of hepatitis C virus (HCV), including thermodynamic prediction, calculation of free energy on folding, and a newly developed method to scan sequences for covariant sites and associated secondary structures using a parsimony-based algorithm. Each of the prediction methods provided evidence for complex RNA folding in the core- and NS5B-encoding regions of the genome. The positioning of covariant sites and associated predicted stem-loop structures coincided with thermodynamic predictions of RNA base pairing, and localized precisely in parts of the genome with marked suppression of variability at synonymous sites. Combined, there was evidence for a total of six evolutionarily conserved stem-loop structures in the NS5B-encoding region and two in the core gene. The virus most closely related to HCV, GB virus-B (GBV-B) also showed evidence for similar internal base pairing in its coding region, although predictions of secondary structures were limited by the absence of comparative sequence data for this virus. While the role(s) of stem-loops in the coding region of HCV and GBV-B are currently unknown, the structure predictions in this study could provide the starting point for functional investigations using recently developed self-replicating clones of HCV.     

Specific secondary structures in the capsid-coding region of giardiavirus transcript are required for its translation in Giardia lamblia

Enhanced translation of giardiavirus (GLV)-luciferase chimeric mRNA in Giardia lamblia requires the presence of the initial 264 nucleotides of the viral capsid-coding region. A 13 nt downstream box (DB) sequence within this region, complementary to a 15 nt sequence near the 3' end of G. lamblia 16 S-like ribosomal RNA (rRNA), was found to be essential for the enhanced translation. However, DB is located 64-78 nt downstream of the initiation codon, whereas an exponential increase of translation efficiency depends on a further increment of the coding region from nucleotides 111 to 264. Thus, there could be additional structural requirements for translation enhancement in the region downstream from DB. Four major stem-loop structures, designated I to IV, were identified in the MFOLD-predicted secondary structure of the 264 nt capsid-coding region with an estimated minimum free energy (DeltaG degrees ) of -77.16 kcal x mol(-1). Our chemical probing analysis of the free 264 nt RNA molecule in solution supports the predicted presence of stem-loops I, II and III, but casts doubts on stem-loop IV. It suggests, instead, the presence of a stem-loop IVA at a nearby location in the molecule. Site-directed mutagenesis designed to disrupt stem-loop structures I, II, III or IVA resulted in drastic reduction of translation efficiency, which was restored by compensatory sequence changes to regenerate individual stem-loop structures. Mutations disrupting the originally designated stem-loop IV did not exert any detectable effect on translation. However, alterations of the sequence UCUCC between nucleotides 216 and 220 in the flexible loop region of the revised secondary structure led to a precipitous drop in translation. Another stem-loop predicted by MFOLD that consists of a major portion of the DB sequence was examined by chemical probing but found little experimental support. Changes of the DB sequence without affecting the postulated stem structure led to drastic losses of translation efficiency. Thus, a simple structural basis for the enhanced translation could be that stem-loops I, II, III and IVA and the UCUCC sequence may facilitate the interaction between DB and the anti-DB in 16 S-like rRNA in initiating translation of GLV mRNA in G. lamblia.     

Identification of potential conserved RNA secondary structure throughout influenza A coding regions

Influenza A is a negative sense RNA virus of significant public health concern. While much is understood about the life cycle of the virus, knowledge of RNA secondary structure in influenza A virus is sparse. Predictions of RNA secondary structure can focus experimental efforts. The present study analyzes coding regions of the eight viral genome segments in both the (+) and (-) sense RNA for conserved secondary structure. The predictions are based on identifying regions of unusual thermodynamic stabilities and are correlated with studies of suppression of synonymous codon usage (SSCU). The results indicate that secondary structure is favored in the (+) sense influenza RNA. Twenty regions with putative conserved RNA structure have been identified, including two previously described structured regions. Of these predictions, eight have high thermodynamic stability and SSCU, with five of these corresponding to current annotations (e.g., splice sites), while the remaining 12 are predicted by the thermodynamics alone. Secondary structures with high conservation of base-pairing are proposed within the five regions having known function. A combination of thermodynamics, amino acid and nucleotide sequence comparisons along with SSCU was essential for revealing potential secondary structures.     

Determination of the secondary structure of and cellular protein binding to the 3'-untranslated region of the hepatitis C virus RNA genome.

Hepatitis C virus (HCV) contains a positive-stranded RNA genome of approximately 9.5 kb. Despite the overall sequence diversity among individual HCV isolates, the 3'-end 98 nucleotides (nt) of the HCV RNA, which constitute part of the 3'-untranslated region (3'-UTR), are highly conserved. This conserved region may contain the cis-acting signals for RNA replication involving possibly both viral and cellular proteins. We carried out RNase digestion studies, which revealed that this 98-nt region contains three stem-loops but may also assume alternative structures. We further performed UV cross-linking experiments to detect cellular proteins that bound to this region. A 58-kDa cellular protein (p58) was detected. Its binding site was mapped to the stem-loops 2 and 3, which are the most conserved region of the 3'-UTR. Site-directed mutagenesis studies revealed that both stem structures and specific nucleotide sequence within the two loops are important for p58 binding. Mutations that disrupted stem structures abolished protein binding, while the compensatory mutations restored its binding. This region also contains partial sequence similarity to the reported consensus binding sequence for polypyrimidine tract-binding protein (PTB) (a 57-kDa protein). The UV-cross-linked protein could be immunoprecipitated with the anti-PTB antibody, and the recombinant PTB bound to the HCV 3'-UTR with the same binding specificity as p58, establishing that this protein is PTB. However, the reported PTB-binding sequence was not sufficient, but rather the entire stem-loops 2 and 3 were required, for PTB binding; thus, its binding specificity is significantly different from the reported PTB-binding sequence requirement. This protein was detected in both the nuclei and cytoplasm of most mammalian cell lines tested and human primary hepatocytes. PTB may participate in the regulation of HCV RNA synthesis or translation.     

Structural and functional analysis of the RNA transport element, a member of an extensive family present in the mouse genome

We previously identified an RNA transport element (RTE), present in a subclass of rodent intracisternal A particle retroelements (F. Nappi, R. Schneider, A. Zolotukhin, S. Smulevitch, D. Michalowski, J. Bear, B. Felber, and G. Pavlakis, J. Virol. 75:4558-4569, 2001), that is able to replace Rev-responsive element regulation in human immunodeficiency virus type 1. RTE-directed mRNA export is mediated by a still-unknown cellular factor(s), is independent of the CRM1 nuclear export receptor, and is conserved among vertebrates. Here we show that this RTE folds into an extended RNA secondary structure and thus does not resemble any known RTEs. Computer searches revealed the presence of 105 identical elements and more than 3,000 related elements which share at least 70% sequence identity with the RTE and which are found on all mouse chromosomes. These related elements are predicted to fold into RTE-like structures. Comparison of the sequences and structures revealed that the RTE and related elements can be divided into four groups. Mutagenesis of the RTE revealed that the minimal element contains four internal stem-loops, which are indispensable for function in mammalian cells. In contrast, only part of the element is essential to mediate RNA transport in microinjected Xenopus laevis oocyte nuclei. Importantly, the minimal RTE able to promote RNA transport has key structural features which are preserved in all the RTE-related elements, further supporting their functional importance. Therefore, RTE function depends on a complex secondary structure that is important for the interaction with the cellular export factor(s).     

The Bovine Leukemia Virus Encapsidation Signal Is Composed of RNA Secondary Structures

The encapsidation signal of bovine leukemia virus (BLV) was previously shown by deletion analysis to be discontinuous and to extend into the 5′ end of the gag gene (L. Mansky et al., J. Virol. 69:3282–3289, 1995). The global minimum-energy optimal folding for the entire BLV RNA, including the previously mapped primary and secondary encapsidation signal regions, was analyzed. Two stable stem-loop structures (located just downstream of the gag start codon) were predicted within the primary signal region, and one stable stem-loop structure (in the gag gene) was predicted in the secondary signal region. Based on these predicted structures, we introduced a series of mutations into the primary and secondary encapsidation signals in order to explore the sequence and structural information contained within these regions. The replication efficiency and levels of cytoplasmic and virion RNA were analyzed for these mutants. Mutations that disrupted either or both of the predicted stem-loop structures of the primary signal reduced the replication efficiency by factors of 7 and 40, respectively; similar reductions in RNA encapsidation efficiency were observed. The mutant with both stem-loop structures disrupted had a phenotype similar to that of a mutant containing a deletion of the entire primary signal region. Mutations that disrupted the predicted stem-loop structure of the secondary signal led to similar reductions (factors of 4 to 6) in both the replication and RNA encapsidation efficiencies. The introduction of compensatory mutations into mutants from both the primary and secondary signal regions, which restored the predicted stem-loop structures, led to levels of replication and RNA encapsidation comparable to those of virus containing the wild-type encapsidation signal. Replacement of the BLV RNA region containing the primary and secondary encapsidation signals with a similar region from human T-cell leukemia virus (HTLV) type 1 or type 2 led to virus replication at three-quarters or one-fifth of the level of the parental virus, respectively. The results from both the compensatory mutants and BLV-HTLV chimeras indicate that the encapsidation sequences are recognized largely by their secondary or tertiary structures.