Structural basis of genomic RNA (gRNA) dimerization and packaging determinants of mouse mammary tumor virus (MMTV)
© Aktar et al.; licensee BioMed Central Ltd. 2014
Received: 8 May 2014
Accepted: 23 October 2014
Published: 14 November 2014
One of the hallmarks of retroviral life cycle is the efficient and specific packaging of two copies of retroviral gRNA in the form of a non-covalent RNA dimer by the assembling virions. It is becoming increasingly clear that the process of dimerization is closely linked with gRNA packaging, and in some retroviruses, the latter depends on the former. Earlier mutational analysis of the 5’ end of the MMTV genome indicated that MMTV gRNA packaging determinants comprise sequences both within the 5’ untranslated region (5’ UTR) and the beginning of gag.
The RNA secondary structure of MMTV gRNA packaging sequences was elucidated employing selective 2’hydroxyl acylation analyzed by primer extension (SHAPE). SHAPE analyses revealed the presence of a U5/Gag long-range interaction (U5/Gag LRI), not predicted by minimum free-energy structure predictions that potentially stabilizes the global structure of this region. Structure conservation along with base-pair covariations between different strains of MMTV further supported the SHAPE-validated model. The 5’ region of the MMTV gRNA contains multiple palindromic (pal) sequences that could initiate intermolecular interaction during RNA dimerization. In vitro RNA dimerization, SHAPE analysis, and structure prediction approaches on a series of pal mutants revealed that MMTV RNA utilizes a palindromic point of contact to initiate intermolecular interactions between two gRNAs, leading to dimerization. This contact point resides within pal II (5’ CGGCCG 3’) at the 5’ UTR and contains a canonical “GC” dyad and therefore likely constitutes the MMTV RNA dimerization initiation site (DIS). Further analyses of these pal mutants employing in vivo genetic approaches indicate that pal II, as well as pal sequences located in the primer binding site (PBS) are both required for efficient MMTV gRNA packaging.
Employing structural prediction, biochemical, and genetic approaches, we show that pal II functions as a primary point of contact between two MMTV RNAs, leading to gRNA dimerization and its subsequent encapsidation into the assembling virus particles. The results presented here enhance our understanding of the MMTV gRNA dimerization and packaging processes and the role of structural motifs with respect to RNA-RNA and possibly RNA-protein interactions that might be taking place during MMTV life cycle.
An essential step in retroviral life cycle is the efficient and specific packaging of two copies of plus-strand full-length genomic RNA (gRNA) into the virus particle from a large pool of cellular and other viral RNAs in the cytoplasm (reviewed in -). Retroviral RNA packaging process involves the recognition of cis-acting packaging signal(s) by the zinc finger(s) in the nucleocapsid (NC) domain of the Gag polyprotein (reviewed in -,). Over the years, it has been shown that the packaging signal(s) of different retroviruses reside at the 5’ end of the gRNA and include continuous or discontinuous sequences from R, U5, and 5′ untranslated (UTR) regions and extend into the 5’ end of the gag gene (reviewed in -,-).
The retroviral RNA genome is packaged as a non-covalent dimer and the processes of dimerization and packaging are closely interlinked (reviewed in -,). The sequences responsible for retroviral gRNA packaging and dimerization of a variety of retroviruses map to the same ~100-400 nucleotides (nt) at the 5’ end of the gRNA, and in most cases are genetically indistinguishable (reviewed in -).
Recent retroviral RNA cross-packaging studies have shown that the specificity of retroviral gRNA packaging can be exchanged by substituting the packaging signals of genetically diverse retroviruses ,. This has further been substantiated by the fact that heterodimers involving RNAs from two divergent retroviruses with no sequence homology can also be packaged (reviewed in ,,). In addition, a number of RNA cross- and co-packaging studies among diverse retroviruses suggest that the process of gRNA dimerization and packaging is likely to involve recognition of structural motifs rather than primary sequences ,,-. Consistent with this, the packaging and dimerization sequences of almost all retroviruses have been shown to assume higher order structures comprising of various structural motifs that have been shown to mediate RNA-RNA and RNA-protein interactions during retroviral RNA dimerization and packaging (reviewed in -,.
Retroviral RNA dimerization is usually mediated by a palindromic (pal) sequence known as the dimerization initiation site (DIS) present in the 5’ region of the gRNA. In a number of retroviruses, the DIS has been shown to assume a hairpin structure containing a canonical “GC” dyad which interacts with the DIS loop on the second gRNA copy, resulting in a kissing loop interaction -; (further reviewed in ). Pal sequences involved in RNA dimerization have been identified in human and simian immunodeficiency viruses (HIV-1, HIV-2, and SIV) (reviewed in ,,,-), feline immunodeficiency virus (FIV) ,, Mason-Pfizer monkey virus (MPMV) ,, and avian leukosis virus (ALV) . In addition to gRNA dimerization and packaging, mutations in the DIS have been shown to affect other steps in retroviral life cycle such as reverse transcription ,,,, and recombination -. Since the DIS has been shown to regulate the retroviral life cycle, it is therefore an attractive target for antiretroviral drugs, such as aminoglycosides ,.
Despite having been studied extensively, little is known about the molecular mechanisms of gRNA dimerization and packaging during the mouse mammary tumor virus (MMTV) life cycle ,. An earlier study by Salmons et al.,  suggested that MMTV harbors sequences responsible for gRNA packaging in the 5’ region of its genome. Employing a biologically relevant in vivo packaging and transduction assay for MMTV , we have recently identified a continuous region, spanning the beginning of R to 120 nt of gag that was found to be critical for MMTV gRNA packaging and propagation . Folding algorithms predicted that these sequences fold into a higher order structure comprising of a number of stable structural motifs.
In this study, we employed the technique of selective 2’hydroxyl acylation analyzed by primer extension (SHAPE) - to validate and further refine the predicted higher order features of the MMTV packaging signal. The SHAPE-validated RNA secondary structure revealed four prominent pal sequences (pal I, pal II, primer binding site (PBS) pal, and pal III), and one or more of them could potentially be involved in initiating intermolecular interactions during MMTV gRNA dimerization. A systematic deletion/substitution analysis of different pal mutants using in vitro dimerization, in vivo packaging and transduction assays revealed that RNA dimerization, packaging and propagation processes of pal II mutants was greatly compromised. Results obtained from these complementary approaches indicate that pal II serves as the primary palindromic point of contact during intermolecular interactions, leading to the dimerization and RNA packaging processes during MMTV life cycle.
Predicted secondary structure of the MMTV 5’ gRNA
Validation of the secondary structure by SHAPE
To experimentally validate the predicted structure of the MMTV packaging signal RNA, we employed the SHAPE methodology -. Since ~50 nt at the 3’ end of the RNA molecule is lost due to primer binding for cDNA synthesis; therefore, to ensure that the entire 432 nt region of the 5’ end of the MMTV genome could be properly analyzed by SHAPE, we cloned a larger piece of RNA upto 712 nt in the wild type clone (SA035; Figure 1B) for in vitro transcription. Briefly, RNA transcribed from SA035 was purified, incubated in the dimer buffer, and modified with BzCN, as described earlier . BzCN selectively acylates the 2’ hydroxyl group of ribose within nucleotides in flexible (unpaired) regions. Reactivity at each nucleotide was calculated by subtraction of the reverse transcription product of unmodified RNA from that of BzCN-modified RNA ,. The SHAPE reactivity data obtained for each nucleotide from 3-6 independent experiments were then applied as pseudo-energy constraints in the structure prediction program “RNAstructure” in order to develop a refined RNA structural model of the sequences required for MMTV gRNA dimerization and packaging.
Consistent with the predicted structure, most of the loops and bulges containing unpaired nucleotides in the MMTV packaging signal RNA showed SHAPE reactivity (Figure 1C). We also observed moderate and high reactivities of nucleotides forming G-U base pairs, especially at the end of helices (G120, U320, G193, U294, G324, U356, U357, U358, U415, Figure 1C). Some base paired nucleotides (G21, C25, G84) located on either side of a bulge or loop also showed certain amount of reactivity, as expected. The refined structural model corroborated well with the predicted structure, supporting the existence of major structural motifs (SL1-SL6; Figure 1C). SHAPE also confirmed the existence of the bifurcated SL4 structure containing the 9 nt stretch of ssPurines and pal II, as predicted by Mfold (Figure 1C).
Even though a high conformity was observed between the predicted and the SHAPE-validated structural models, some conspicuous dissimilarity could be noted. In the SHAPE validated structure, SL1 and SL6 were shorter, resulting in the formation of SL7 which was not observed in the predicted structure (compare Figure 1A with Figure 1C). In addition, the sequences 5’ to the shortened SL6 form a new stem loop which is different from the predicted SL5 (Figure 1A) and therefore was labeled as SL5’ (Figure 1C). Interestingly, the single-stranded bulges and loops in and around SL3 did not show significant SHAPE reactivity despite being unpaired in the secondary structure model. This could possibly suggest that these nucleotides are involved in tertiary or inter/intramolecular interactions. Furthermore, the SHAPE-validated structure harbored a long-range interaction (LRI) involving complementary sequences from U5 and Gag (U5/Gag LRI) which was not predicted by Mfold (compare Figure 1A with Figure 1C). This U5/Gag LRI is made up of three short helices separated by a bulge and an internal loop (Figure 1C).
The SHAPE-validated RNA structural model is supported by phylogeny
Sequence alignment of eight strains of MMTV (sequence from +1R-120 nt gag) revealed a high degree of conservation relative to the wild type mtv-1 (AF228550.1) strain, especially of sequences of the major structural motifs such as U5/Gag LRI, ssPurines, pal II, and PBS pal (Additional file 1A) once superimposed on the SHAPE-validated structure (Additional file 1B). Only a small number of nucleotides were observed to be variable (<75% conservation as shown in Additional file 1B). Finally, we generated a consensus structure of this region using RNAalifold (Additional file 2). Consistent with the SHAPE-validated structure (Figure 1C), consensus RNA secondary structure predictions (Additional file 2) revealed that the major structural motifs were consistently maintained. In addition, the consensus structures predicted by RNAalifold also supported the existence of SHAPE-validated short SL6 and SL5’ over the Mfold predicted long SL6 and SL5 motifs (compare Figure 1A, Figure 1C, and Additional file 2). The consensus structures also revealed the existence of LRI which is very similar to the one observed in the SHAPE-validated structure (Figure 1C and Additional file 2).
Identification of the pal sequence mediating MMTV gRNA dimerization and packaging
The SHAPE-validated region found to be crucial for MMTV gRNA packaging revealed four pal sequences (Figure 1C) which could possibly play role(s) in intermolecular RNA interaction(s) to augment dimerization and subsequent packaging of the dimerized genome. Therefore, to ascertain which pal may mediate such interactions, we performed a systematic mutational analysis of these pal sequences and tested their effects on both RNA dimerization and packaging employing a combination of genetic and biochemical approaches.
Mutational analysis of pal I reveals that it does not play a role in MMTV gRNA dimerization and packaging
To further evaluate the functional impact of mutations in pal I, we tested these mutants in a biologically-relevant in vivo packaging and propagation assay (Additional file 3). This was accomplished by cloning these mutations into the MMTV subgenomic transfer vector, DA024, creating SA031P-SA033P pal I mutant transfer vectors (Figure 2A). These transfer vectors were tested by transfecting human embryonic kidney (HEK) 293T cells along with the packaging and vesicular stomatitis virus glycoprotein (VSV-G) envelope expression constructs (Additional file 3). The wild type and mutant viral particles produced were used to: 1) measure the packaged viral RNA using qPCR, and 2) assess the propagation of the packaged RNA in the target cells which was monitored by the appearance of hygromycin resistant colonies. We first checked the transfection efficiency of the cultures containing mutant transfer vectors. The transfection efficiencies from at least three independent experiments were observed to be within 2 folds of each other, suggesting that all cultures were efficiently transfected (Figure 2C). These transfection efficiencies were taken into account to normalize and calculate the RNA packaging efficiency and propagation data for each mutant relative to the wild type.
To determine the RNA packaging efficiency, cDNAs were prepared from cytoplasmic and pelleted virion RNAs. Before making the cDNAs these RNAs were subjected to DNase treatment to ensure that there was no contaminating plasmid DNA, followed by PCR using virus-specific primers. Following amplification, the lack of a positive signal indicated that the DNA contamination in our RNA preparations was below the detection level (Figure 2D, panels i and ii). After having confirmed this, RNA preparations were reverse transcribed and cDNAs were prepared. Since the relative packaging efficiency data is expressed in relation to the efficiently and stably expressed viral RNAs that are exported from the nucleus to the cytoplasm, we monitored the integrity of the cytoplasmic RNA fraction by checking for the absence of unspliced β-actin mRNA by RT-PCR (Figure 2D, panel iii). To ensure that each cytoplasmic RNA fraction contained amplifiable cDNAs, the unspliced β-actin PCRs were conducted in the presence of primers/competimer for 18S ribosomal RNAs as an internal control ,,. In a separate amplification reaction, the same cytoplasmic cDNA samples were amplified using spliced β-actin primers (Figure 2D, panel iv). Figure 2D shows that while 18S ribosomal RNA and spliced β-actin mRNA were amplifiable, unspliced β-actin mRNA was not. These results indicated that the fractionation technique employed was robust and that the nuclear membrane integrity was maintained during the fractionation process.
In order to monitor the effect of pal I mutants on gRNA packaging, the relative expression of transfer vector RNAs in the cytoplasm was tested using a custom-made MMTV real time PCR assay in combination with a commercially available-endogenous β-actin Taqman assay as described earlier . Following normalization of the RNA packaged into the virus particles to the cytoplasmic transfer vector RNA expression and to transfection efficiencies, the packaging efficiency of the transfer vector RNAs into the virus particles was measured (see Methods). The in vivo gRNA packaging data from multiple experiments revealed that none of the pal I mutations significantly decreased gRNA packaging (P >0.01; Figure 2E). This corroborated well with the dimerization data for the same mutants which showed no significant dimerization defect (Figure 2B). Interestingly, when viral particles containing RNAs from these mutants were used to transduce target cells with a hygromycin resistance gene, a drastic reduction of RNA propagation (6-20 fold reduction of the relative CFU/ml, P ≤0.001) was observed for all mutants (Figure 2F). Lack of RNA propagation of the packaged mutant transfer vector RNAs suggests that the introduced mutations impaired post packaging events of viral life cycle such as reverse transcription and/or integration since our readout assay is dependent on successful completion of these steps. Consistent with this hypothesis, mutations in close vicinity of the PBS (as has been the case for these pal I mutants) have been shown to greatly impinge reverse transcription ,,.
Mutational analysis of pal II and the PBS pal suggests that there could be two points of contact between the gRNAs leading to MMTV gRNA dimerization
Role of pal II and PBS pal in regulating MMTV gRNA packaging and propagation
Next, we analyzed the RNA packaging and propagation efficiencies of the pal II and PBS pal mutants by cloning these mutations into the MMTV subgenomic transfer vector, DA024, similar for the pal I mutants, creating SA041P-SA047P and SA051P (Figure 3A and Figure 5A). Data from several experiments showed that pal II and PBS pal mutants’ transfection efficiencies were within two folds and that the transfer vector RNAs were expressed efficiently (data not shown). When similar amounts of virions were used to isolate packaged RNA, all of the mutants tested except for SA051P, were found to be severely impaired for gRNA packaging (27-150 fold reduction, P < 0.001, Figure 3D and Figure 5D). Although packaged, the relative packaging efficiency of SA051P mutant was reduced by 9-fold compared to the wild type (Figure 3D). Interestingly, deletion of either pal II or the PBS pal had a pronounced effect on RNA packaging (even though the effect of pal II deletion was more dramatic on packaging; Figure 3), suggesting that the MMTV Gag precursor initially recognizes an RNA dimer with two points of contact.
Propagation of these mutants, as measured by counting the number of hygromycin resistant colonies (CFU/ml) in the infected cultures correlated well with the significantly-reduced RNA packaging data (Figure 3E and Figure 5E). Lack of propagation of SA051P mutant can be explained despite residual RNA packaging by the deletion in the PBS, since our assay requires the packaged RNA to be successfully reverse transcribed (initiated by an intact PBS) before integration to allow expression of the hygromycin resistant gene and subsequent appearance of hygromycin-resistant colonies.
Altogether, the dimerization and packaging results suggest that the pal II and PBS pal play important roles in both MMTV gRNA dimerization and packaging. However, the case of mutant SA047P containing substitution of pal II with HIV-1 pal along with flanking purines was intriguing as its RNA packaging and propagation were abolished (Figure 5D and Figure 5E), despite displaying wild type RNA dimerization levels (Figure 5B and Figure 5C). These results suggest that RNA dimerization may not be sufficient to ensure RNA packaging and propagation and further indicate that pal II could also be involved in RNA-protein interactions during genome encapsidation.
SHAPE analyses of pal II mutants support their role in MMTV gRNA dimerization
Next we analysed the structure of the SL2 domain, which contains the PBS pal, in the wild type and mutant RNAs. Although the PBS is predicted to be unpaired in monomeric RNA, the overlapping PBS pals were weakly reactive (Figure 7A), especially the second pal (5’ GGCGCC 3’), consistent with their potential role in RNA dimerization by forming loose dimers. The reactivity of this region was not significantly affected by deletion or substitutions in pal II (Figure 7A, mutants SA042 and SA047), suggesting the second pal within the PBS could potentially play a role in RNA dimerization (Figure 3B, 3C and Figure 5B, 5C).
Structural prediction of mutant dimers supports dimerization data
Additional mutational analysis of pal II, PBS pal, and pal III confirms the primary role of pal II in MMTV gRNA dimerization
Finally, structural analysis of the double deletion mutant (SA046; pal II and PBS pal) predicted that this mutant RNA could homodimerize via pal III (Figure 4D). Since dimerization was found to be abrogated for this mutant (Figure 3C), it is unlikely that these dimers are being formed during MMTV replication. To test this experimentally, we created a mutant in which pal III was deleted (SA061; Figure 10A). Test of this mutant revealed that deletion of pal III did not have any effect on dimerization (Figure 10B and Figure 10C), suggesting pal III has no role in MMTV RNA dimerization. Overall, these results suggest that pal II is the primary RNA dimerization determinant in MMTV.
In an attempt to enhance our understanding of the gRNA dimerization and packaging processes during MMTV life cycle, we first validated the predicted structure of MMTV packaging signal RNA employing SHAPE. The SHAPE-validated structure was similar, but not identical to the Mfold-predicted structure. An important feature revealed by SHAPE was a U5/Gag LRI, which was also supported by phylogeny and in silico RNA structure predictions (Figure 1, Additional files 1 and 2). The U5/Gag LRI anchors the overall structure of MMTV packaging signal RNA and is consistent with an earlier report showing that the MMTV packaging sequence extends into gag . Similar LRIs exist in the packaging signal RNAs of a number of retroviruses ,,- and in several cases, disruption of LRIs have been shown to affect RNA dimerization and packaging ,,. In addition to LRI, SHAPE validated the existence of the SL2 and SL4 domains. The SL2 domain contains the PBS and the pal I sequence, and a similar large SL2 domain is found in a number of other retroviral packaging signal RNAs ,,,,. The bifurcated SL4 hairpin contains a pal helix loop (pal II) and an ssPurines loop (Figure 1). Both SL2 and SL4 structures are conserved in all the 8 MMTV strains studied (Additional files 1 and 2).
The presence of a stretch of purines in the packaging sequences on retroviral gRNA has been proposed to facilitate RNA packaging by functioning as a potential NC binding site ,,,. Thus, it is interesting to note that the bifurcated SL4 domain contains a purine-rich loop (5’ GGAGAAGAG 3’) adjacent to the pal II helix loop (Figure 1), reminiscent to the situation found in MPMV ,. In the case of MPMV, genetic and structure-prediction analyses have suggested that either ssPurines or its partial repeat base-paired sequence in an adjacent region plays a crucial role in gRNA packaging, possibly by functioning as a NC binding site . Therefore, it will be interesting to investigate the role of ssPurines in MMTV gRNA packaging.
Retroviral gRNA dimerization has been suggested as a prerequisite for packaging into the budding virions (reviewed in -,). This process is usually initiated by conserved palindromic sequences that initiate kissing loop interactions ,,,,,-,-,,. Nucleotide sequence and structural analyses of MMTV packaging signal RNA initially revealed four pal sequences (pal I, pal II, PBS pal, and pal III, Figure 1). Our mutational analysis clearly indicates that neither pal I (Figure 2) nor pal III (Figure 10) is involved in RNA dimerization and packaging. Nevertheless, pal I mutants had significant propagation defects (Figure 2), suggesting that mutations close to the PBS might affect annealing of the primer tRNA or initiation of reverse transcription itself. On the other hand, initial mutational analysis suggested that both pal II and PBS pal play a role in MMTV RNA dimerization (Figure 3 and Figure 5). Indeed, mutations in pal II (SA042) as well as PBS pal (SA051) reduced RNA dimerization by 2-3-folds, while deleting both pal II and PBS pal (SA046) resulted in nearly a complete abrogation of dimerization (Figure 3). However, additional PBS pal mutants (Figure 9) showed that unlike PBS pal deletion (mutant SA046; Figure 3B, 3C), substitutions in PBS pal had no effect on RNA dimerization, even when the wild type sequence was replaced by a non-palindromic sequence (mutants SA057, SA058, SA059 and SA060; Figure 9). These results indicate that the PBS pal does not mediate intermolecular interactions by direct base pairing and strongly suggest that the effect of the PBS pal deletion on RNA dimerization is indirect. Interestingly, deletion of either pal II or the PBS pal had a pronounced effect on RNA packaging (even though the effect of pal II deletion was more dramatic on packaging; Figure 3), suggesting that the MMTV Gag precursor initially recognizes an RNA dimer with two points of contact.
A classical way of demonstrating the existence of any intermolecular interaction between viral RNAs is by using trans-complementary mutants ,-,. Our initial attempts to restore RNA dimerization using trans-complementary pal II mutants were not successful (Figure 5). This unexpected result could however be explained by misfolding of one of the mutants (SA045; Figure 6B and Figure 8B). These two mutants were predicted to form heterodimer via the PBS pal rather than via their trans complementary substituted sequences (Figure 8C). Previous work with HIV-1 that showed that purines flanking the palindromic sequence can enhance the stability of the RNA kissing-loop dimer  led us to design additional trans-complementary mutants with such purines (Figure 10, mutants SA052 and SA053). While dimerization of these RNA mutants was severely affected when they were incubated alone, they dimerized in trans to wild type levels, demonstrating that pal II mediates MMTV RNA dimerization by direct base pairing between two genomic RNA molecules. Along the same lines, results obtained with mutant SA047 (in which pal II was replaced by the HIV-1 pal with the flanking purines) displayed wild type RNA dimerization levels; however, its RNA packaging and propagation were abolished (Figure 5), suggesting that the pal II hairpin could also be involved in RNA-protein interactions during genome encapsidation, perhaps at the primary sequence level. The overall 3D structure of the packaging domain would thus probably be important for Gag binding. This could explain why deletion of 5’ UTR sequences forming the basal part in the SHAPE-validated structure of SL2 abrogated MMTV gRNA packaging and propagation .
Results presented in this study validate the higher order structural features of the MMTV packaging signal RNA using the SHAPE methodology, which includes LRI involving complementary sequences from U5 and the gag (U5/Gag LRI), ssPurines, and four palindromic sequences (Figure 1). Employing a combination of genetic, biochemical, and structural prediction approaches, we further show that the MMTV genome uses pal II as the primary point of contact, by likely functioning as the DIS, to initiate intermolecular interactions augmenting efficient gRNA dimerization, and leading to its encapsidation in the assembling virus particles (Figure 3 and Figure 4, and Figure 10). SHAPE analyses of pal II mutants further support its role in MMTV gRNA dimerization (Figure 6). These results are in agreement with earlier studies on determinants of MMTV gRNA packaging, and further reveal important functional relationships between gRNA dimerization and packaging processes, suggesting that MMTV prefers a dimeric genome as a packaging substrate.
The secondary structure of the 5’ region of the MMTV gRNA containing either 120 nt of gag (432 nt of the 5’ end from R) or 400 nt of gag (containing 712 nt of the 5’ end from R) were predicted using the Mfold server ,. The SHAPE reactivity data of the first 432 nt region was applied in RNAstructure software (version 5.3)  as constraints to develop the final structural model. The Mfold-predicted and SHAPE-validated structures were redrawn with the XRNA software (http://rna.ucsc.edu/rnacenter/xrna/xrna.html). MMTV nucleotide positions refer to HYBMTV a molecular clone created by Shackleford and Varmus .
Nucleotides 1-432 of the gRNA of different strains of MMTV were aligned using Clustal Omega , and a consensus or “conservation annotation” structure for these strains was predicted using RNAalifold server ,. The accession numbers for different MMTV strains used for sequence alignment are AF228550.1 (mtv-1) , L37517.1 (mtv-6) , AF033807.1 , M15122.1 (BR6) , D16249.1 (JYG) , AF228551.1 (C3H/HeJ) , X00018.1 (GR)  and AF228552.1 (C3H) .
Construction of plasmids
The 5’ end of MMTV gRNA corresponding to nt 1-712 was amplified by polymerase chain reaction (PCR) using the MMTV sub-genomic transfer vector DA024  as the template and the sense (S) OTR 984 (5’ CCCAAGCT T AATACGACTCACTATAGGGCAACAGTCCTAATATTCACG 3’) and the antisense (AS) OTR 985 (5’ AAACCCGGGTTCCCCTGGTCCCATAAG 3’) primers. These primers include the T7 promoter sequence (in bold) at the 5’ end along with flanking restriction sites (HindIII and XmaI/SmaI, shown in italics). The HindIII and XmaI sites were utilized to introduce the MMTV sequence into a pUC-based cloning vector (pIC19R) . Introduction of these sequences (from OTRs 984 and 985) resulted in the addition of GGG and CCC complementary sequences at the ends of the RNA which created base pairing potential; however, folding of the 5’ end of the MMTV genome containing these artificially-introduced complementary sequences in Mfold did not predict base pairing among these sequences nor alter the overall RNA secondary structure of this region (data not shown). The sequence of the resulting clone, SA035, was confirmed by sequencing. SA035 (containing 400 nt of gag) was also used as the wild type clone for in vitro gRNA dimerization assays.
A number of mutations, including deletions, substitutions and substitution containing heterologous trans-complementary sequences, were introduced in pal I (5’ GUCGGCCGAC 3’), pal II (5’ CUGCAG 3’), PBS pal (5’ CAGCUGGCGCC 3’), and pal III (5’ CCUCUCAGAGAGG 3’) through splice overlap extension (SOE) PCR as described previously ,,,. Briefly, in round one, the MMTV sub-genomic vector DA024 was amplified in two separate PCRs using S and AS primers generating products having overlapping complementary sequences. In round two, the products from round one PCRs, which anneal via their overlapping complementary sequences, were amplified using outer primers OTR 984 (S) and OTR 985 (AS), generating a final product containing the desired mutation(s). The sequences of S and AS primers used to introduce mutations are provided in Additional file 4, along with the clone names and the necessary description. Some mutations were also cloned in the MMTV subgenomic transfer vector, DA024, to study the effect of these mutations on RNA packaging and propagation. All clones were confirmed by sequencing. Further details of cloning can be obtained from the authors upon request.
In vitrotranscription assay
After linearization with the SmaI restriction enzyme, the wild type (SA035) or mutant (SA031-SA033, SA041-SA045, SA046, SA047 and SA051) clones were transcribed in vitro to generate the RNAs to be used in dimerization assays, as previously described . Following DNase treatment, RNA extraction and purification, the RNAs were concentrated using Amicon Ultra-4 10 K devices (Millipore) and their concentration was determined using nanodrop (ThermoScientific), as described previously . Finally, the quality of the purified RNAs was checked by denaturing polyacrylamide gel electrophoresis.
SHAPE - was performed on in vitro transcribed RNAs using benzoyl cyanide (BzCN) following the protocol described recently . Briefly, RNAs were modified in the dimer buffer (final concentrations: 50 mM sodium cacodylate, (pH 7.5), 300 mM KCl and 5 mM MgCl2) in the presence of 2 μg total yeast tRNA (Sigma Aldrich). The dimer buffer was used for the SHAPE analysis to help identify points of contact involved in the dimerization process. The modified RNAs were reverse transcribed using two sets of AS primers: OTRs 9, 10, 11 and 12 (5’ AACAGATTTGGCTTCTGCGG 3’; MMTV nt 614-633) and OTRs 13, 14, 15 and 16 (5’ AGTTTCTGCCCTTTTGAGCC 3’; MMTV nt 325-344), labelled with VIC, FAM, NED, or PET, respectively. These two sets of primers enabled the revere transcription of the entire RNA sequence except for a few nucleotides at the 3’ end due to primer binding. The primer extension products were loaded onto an Applied Biosystems 3130xl genetic analyzer and the electropherograms were analysed with the SHAPEfinder ,. The SHAPE reactivities were determined and normalized as described earlier ,. The normalized SHAPE reactivity data from at least 3-6 independent experiments (submitted as Additional file 5) were used as constraints to validate the RNA secondary structure of the MMTV packaging signal with RNAstructure (version 5.3)  and also to determine the effects of pal mutations on the overall higher order structure.
In vitrodimerization assay
In vitro dimerization assays of wild type and mutant MMTV RNAs were performed as described previously ,. Briefly, 300 nM RNAs were heated at 90°C and incubated in ice for 2 min followed by a 37°C incubation for 30 min in either the dimer (final concentrations: 50 mM sodium cacodylate, pH 7.5; 300 mM KCl; 5 mM MgCl2) or monomer (final concentrations: 50 mM sodium cacodylate, pH 7.5; 40 mM KCl; 0.1 mM MgCl2) buffer. The samples were divided equally and analysed by electrophoresis in TBM (50 mM Tris base, 45 mM boric acid, 0.1 mM MgCl2) and TB (50 mM Tris base, 45 mM boric acid) ethidium bromide prestained agarose gels at 4°C or 20°C, respectively. Density of the dimeric and monomeric species was quantitated using Quantity One software, and percentage of dimerization was calculated and expressed as dimerization relative to the wild type for each mutant, as described previously . For the in vitro dimerization assay of trans complementary mutants, equal amounts (150 nM) of in vitro transcribed RNAs from both the mutants (for a total of 300 nM) were incubated together before electrophoresis, as described above.
In vivogenetic complementation assay
The in vivo genetic trans-complementation assay used two expression plasmids, JA10 and MD.G, to provide the structural and enzymatic genes from MMTV Gag/Pol and vesicular stomatitis virus (VSV-G) envelope regions, respectively ,. A third plasmid (MMTV sub-genomic transfer vector-DA024) containing minimum cis-acting sequences required for RNA packaging, reverse transcription, and integration, served as a source of the packageable vector RNA . Design and rationale of the MMTV 3-plasmid genetic complementation assay developed earlier by our group is schematically shown in Additional file 3 . Briefly, the particles produced by JA10 and MD.G allow packaging of the transfer vector RNA by virtue of the presence of the conventional packaging signal (ψ) on the MMTV transfer vector (DA024). The pal I, pal II, and PBS pal mutations introduced into the T7 promoter-based vector were further cloned into the MMTV sub-genomic transfer vector, DA024, and tested in this assay to monitor their effects on gRNA packaging and propagation.
Human embryonic kidney (HEK) 293T producer cells were transfected using the calcium phosphate transfection protocol to produce virus particles containing the wild type or mutated RNAs, as described previously ,. The luciferase expression vector, pGL3 (Promega), was used as an internal control to measure transfection efficiency, as described earlier -,. Supernatants containing viral particles were collected and spun down at 4000 rpm on a table top centrifuge for 10 min. A portion of the supernatant was used to infect the human cervical cancer cell line (HeLaT4) to determine the transduction efficiency by selection with medium containing 200 μg/ml of hygromycin B (Hyclone). Cells transduced with virus particles containing the packageable transfer vector RNA express the marker Hyg r gene, resulting in colonies which are measured as colony forming units (CFU/ml). The CFU/ml values were then normalized to the transfection efficiency. The normalized propagation efficiencies were then divided by the wild type values to represent the propagation of the mutants relative to the wild type (Relative CFU/ml).
Nucleocytoplasmic fractionation, isolation of RNA, and cDNA preparation
Cells from the transfected cultures were fractionated into nuclear and cytoplasmic fractions as described previously ,. Cellular RNAs from the cytoplasmic fractions and packaged viral RNA were isolated from the pelleted viral particles using Trizol and Trizol LS reagents, respectively (Invitrogen), as described earlier ,. The extracted RNAs were amplified for 30 cycles by conventional PCR using vector-specific primers OTR 671(S) or and OTR 672 (AS) (Additional file 4) after DNase (Promega) treatment to ensure that there was no contaminating DNA present in the RNA samples. The PCR conditions were as follows: denaturation for 2 min at 94°C, followed by 30 cycles of 45 sec at 90°C, annealing for 45 sec at 60°C and extension for 45 sec at 72°C followed by incubation at 4°C. The DNase-treated cytoplasmic and viral RNA samples were reverse transcribed into cDNA and PCR amplified to determine the quality of the cDNA samples. To determine the integrity of the nuclear membrane during cytoplasmic RNA extraction, a multiplex PCR was performed on the cDNAs prepared from the cytoplasmic RNA fractions using OTR 582 (S) and OTR 581 (AS) (Additional file 4) designed to amplify unspliced β-actin mRNA ,. Amplifications of the cDNAs using primers/competimer for 18S ribosomal RNA (18S Quantum competimer control, Ambion) was also performed as an ancillary control for the presence of amplifiable cDNA in the multiplex PCRs ,.
Real time quantitative PCR (qPCR) for transfer vector RNA packaging efficiency
To determine the packaging efficiency of the transfer vector RNA, quantitative qPCR was performed on the cDNA prepared from the cytoplasmic and viral RNAs as described earlier . The endogenous control for qPCR was β-actin (Human β-actin Endogenous Control assay, VIC⁄MGB probe, Applied Biosystems). The primers and probes for MMTV transfer vector RNAs in the FAM-based Applied Biosystems custom expression assay anneal at a 68 nt region in the MMTV U5 (nt 1192-1259) common to wild type and mutant RNAs. The relative amplification efficiencies of the MMTV and β-actin expression assays were determined to be 0.0126, thus, validating the assays for quantitative PCR analysis using the ΔΔCt method . Equal amounts of viral and cytoplasmic cDNA were tested in triplicates. The amplification conditions used were as follows: Uracil-N-glycosylase incubation 2 min at 50°C, denaturation for 10 min at 95°C, followed by 40 cycles of denaturation and annealing/extension steps at 95°C for 15 sec and 60°C for 1 min. The relative packaging efficiencies were derived following a modification of our earlier described method . Briefly, the raw Ct (cycle threshold) values for the MMTV transfer vector were normalized by the Ct values of the β-actin endogenous control. The normalized values were then calculated relative to the mock sample (relative quantification, RQ) that was transfected with all the plasmids except the transfer vector plasmid forming the empty virus particles without any packaged RNA. For final relative packaging efficiency calculations, the relative viral RQ values were divided by the relative cytoplasmic RQ values, and the resulting values were normalized to the transfection efficiencies. The final values are presented relative to the wild type.
The standard paired, two-tailed Student’s t-test was performed between the wild type and mutant clones to determine statistically significant differences. For in vitro dimerization assays, and in vivo packaging and propagation studies, a p value of ≤0.05 was considered significant.
This research was funded primarily by grants from Terry Fox Fund for Cancer Research (2009-06) and the United Arab Emirates University (UAEU-NRF 26) and, in part, by a grant from the College of Medicine and Health Sciences (NP NP/10/22) to TAR. RM is supported by the Centre National de la Recherche Scientifique (C.N.R.S.). SJA was supported through a scholarship from the Office of the Deputy Vice Chancellor for Research and Graduate Studies at the United Arab Emirates University (UAEU).
- D’Souza V, Summers MF: How retroviruses select their genomes. Nat Rev Microbiol. 2005, 8: 643-655. 10.1038/nrmicro1210.View ArticleGoogle Scholar
- Johnson SF, Telesnitsky A: Retroviral RNA dimerization and packaging: the what, how, when, where, and why. PLoS Pathog 2010, 6:e1001007.,Google Scholar
- Lever AML: HIV RNA packaging. Adv Pharmacol. 2007, 55: 1-32. 10.1016/S1054-3589(07)55001-5.View ArticlePubMedGoogle Scholar
- Lu K, Heng X, Garyu L, Monti S, Garcia EL, Kharytonchyk S, Dorjsuren B, Kulandaivel G, Jones S, Hiremath A, Divakaruni SS, LaCotti C, Barton S, Tummillo D, Hosic A, Edme K, Albrecht S, Telesnitsky A, Summers MF: NMR detection of structures in the HIV-1 5’-leader RNA that regulate genome packaging. Science. 2011, 334: 242-245. 10.1126/science.1210460.PubMed CentralView ArticlePubMedGoogle Scholar
- Paillart JC, Shehu-Xhilaga M, Marquet R, Mak J: Dimerization of retroviral RNA genomes: An inseparable pair. Nat Rev Microbiol. 2004, 2: 461-472. 10.1038/nrmicro903.View ArticlePubMedGoogle Scholar
- Russell RS, Liang C, Wainberg MA: Is HIV-1 RNA dimerization a prerequisite for packaging? Yes, no, probably? Retrovirology 2004, 1:23.,Google Scholar
- Abd El-Wahab EW, Smyth RP, Mailler E, Bernacchi S, Vivet-Boudou V, Hijnen M, Jossinet F, Mak J, Paillart JC, Marquet R: Specific recognition of the HIV-1 genomic RNA by the Gag precursor. Nat Commun 2014, 5:4304. doi:10.1038/ncomms5304.,Google Scholar
- Jaballah SA, Aktar SJ, Ali J, Phillip PS, Al Dhaheri NS, Jabeen A, Rizvi TA: A G-C rich palindromic structural motif and a stretch of single stranded purines are required for optimal packaging of Mason-Pfizer monkey virus (MPMV) genomic RNA. J Mol Biol. 2010, 401: 996-1014. 10.1016/j.jmb.2010.06.043.View ArticlePubMedGoogle Scholar
- Mustafa F, Ghazawi A, Jayanth P, Phillip PS, Ali J, Rizvi TA: Sequences intervening between the core packaging determinants are dispensable for maintaining the packaging potential and propagation of feline immunodeficiency virus transfer vector RNAs. J Virol. 2005, 79: 13817-13821. 10.1128/JVI.79.21.13817-13821.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Mustafa F, Al Amri D, Al Ali F, Al Sari N, Al Suwaidi S, Jayanth P, Philips PS, Rizvi TA: Sequences within both the 5’ UTR and Gag are required for optimal in vivo packaging and propagation of mouse mammary tumor virus (MMTV) genomic RNA. PLoS One 2012, 7:e47088.,Google Scholar
- Moore MD, Hu WS: HIV-1 RNA dimerization: It takes two to tango. AIDS Rev. 2009, 11: 91-102.PubMed CentralPubMedGoogle Scholar
- Al Dhaheri NS, Phillip PS, Ghazawi A, Ali J, Beebi E, Jaballah SA, Rizvi TA: Cross-packaging of genetically distinct mouse and primate retroviral RNAs. Retrovirology 2009, 6:66.,Google Scholar
- Al Shamsi IR, Al Dhaheri NS, Phillip PS, Mustafa F, Rizvi TA: Reciprocal cross-packaging of primate lentiviral (HIV-1 and SIV) RNAs by heterologous non-lentiviral MPMV proteins. Virus Res. 2011, 155: 352-357. 10.1016/j.virusres.2010.09.018.View ArticlePubMedGoogle Scholar
- Motomura K, Chen J, Hu WS: Genetic recombination between human immunodeficiency virus type 1 (HIV-1) and HIV-2, two distinct human lentiviruses. J Virol. 2008, 82: 1923-1933. 10.1128/JVI.01937-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore MD, Fu W, Nikolaitchik O, Chen J, Ptak RG, Hu WS: Dimer initiation signal of human immunodeficiency virus type 1: its role in partner selection during RNA copackaging and its effects on recombination. J Virol. 2007, 81: 4002-4011. 10.1128/JVI.02589-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Parveen Z, Mukhtar M, Goodrich A, Acheampong E, Dornburg R, Pomerantz RJ: Cross-packaging of human immunodeficiency virus type 1 vector RNA by spleen necrosis virus proteins: construction of a new generation of spleen necrosis virus-derived retroviral vectors. J Virol. 2004, 78: 6480-6488. 10.1128/JVI.78.12.6480-6488.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Rizvi TA, Panganiban AT: Simian immunodeficiency virus RNA is efficiently encapsidated by human immunodeficiency virus type 1 particles. J Virol. 1993, 67: 2681-2688.PubMed CentralPubMedGoogle Scholar
- White SM, Renda M, Nam NY, Klimatcheva E, Zhu Y, Fisk J, Halterman M, Rimel BJ, Federoff H, Pandya S, Rosenblatt JD, Planelles V: Lentivirus vectors using human and simian immunodeficiency virus elements. J Virol. 1999, 73: 2832-2840.PubMed CentralPubMedGoogle Scholar
- Yin PD, Hu WS: RNAs from genetically distinct retroviruses can copackage and exchange genetic information in vivo. J Virol. 1997, 71: 6237-6242.PubMed CentralPubMedGoogle Scholar
- Miyazaki Y, Miyake A, Nomaguchi M, Adachi A: Structural dynamics of retroviral genome and the packaging. Frontiers in Microbiology. 2011, 2: 1-9. 10.3389/fmicb.2011.00264.View ArticleGoogle Scholar
- Aktar SJ, Jabeen A, Ali LM, Vivet-Boudou V, Marquet R, Rizvi TA: SHAPE analysis of the 5’ end of the Mason-Pfizer monkey virus (MPMV) genomic RNA reveals structural elements required for genome dimerization. RNA. 2013, 19: 1648-1658. 10.1261/rna.040931.113.PubMed CentralView ArticlePubMedGoogle Scholar
- Clever JL, Wong ML, Parslow TG: Requirements for kissing-loop-mediated dimerization of human immunodeficiency virus RNA. J Virol. 1996, 70: 5902-5908.PubMed CentralPubMedGoogle Scholar
- Kenyon JC, Tanner S, Legiewicz M, Phillip PS, Rizvi TA, LeGrice S, Lever AM: SHAPE analysis of the FIV leader RNA reveals a structural switch potentially controlling viral packaging and genome dimerization. Nucleic Acids Res. 2011, 39: 6692-6704. 10.1093/nar/gkr252.PubMed CentralView ArticlePubMedGoogle Scholar
- Paillart JC, Marquet R, Skripkin E, Ehresmann B, Ehresmann C: Mutational analysis of the bipartite dimer linkage structure of human immunodeficiency virus type 1 genomic RNA. J Biol Chem. 1994, 269: 27486-27493.PubMedGoogle Scholar
- Paillart JC, Berthoux L, Ottmann M, Darlix JL, Marquet R, Ehresmann B, Ehresmann C: A dual role of the putative RNA dimerization initiation site of human immunodeficiency virus type 1 in genomic RNA packaging and proviral DNA synthesis. J Virol. 1996, 70: 8348-8354.PubMed CentralPubMedGoogle Scholar
- Paillart JC, Westhof E, Ehresmann C, Ehresmann B, Marquet R: Non-canonical interactions in a kissing loop complex: the dimerization initiation site of HIV-1 genomic RNA. J Mol Biol. 1997, 270: 36-49. 10.1006/jmbi.1997.1096.View ArticlePubMedGoogle Scholar
- Hussein ITM, Ni N, Galli A, Chen J, Moore MD, Hu WS: Delineation of the preferences and requirements of the human immunodeficiency virus type 1 dimerization initiation signal by using an in vivo cell-based selection approach. J Virol. 2010, 84: 6866-6875. 10.1128/JVI.01930-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Baig TT, Lanchy JM, Lodmell JS: HIV-2 RNA dimerization is regulated by intramolecular interactions in vitro . RNA. 2007, 8: 1341-1354. 10.1261/rna.483807.View ArticleGoogle Scholar
- Berkhout B, van Wamel JL: Role of the DIS hairpin in replication of human immunodeficiency virus type 1. J Virol. 1996, 70: 6723-6732.PubMed CentralPubMedGoogle Scholar
- Lanchy JM, Ivanovitch JD, Lodmell JS: A structural linkage between the dimerization and encapsidation signals in HIV-2 leader RNA. RNA. 2003, 9: 1007-1018. 10.1261/rna.5590603.PubMed CentralView ArticlePubMedGoogle Scholar
- Lanchy JM, Lodmell JS: An extended stem-loop 1 is necessary for human immunodeficiency virus type 2 replication and affects genomic RNA encapsidation. J Virol. 2007, 7: 3285-3292. 10.1128/JVI.02025-06.View ArticleGoogle Scholar
- Laughrea M, Jetté L, Mak J, Kleiman L, Liang C, Wainberg MA: Mutations in the kissing-loop hairpin of human immunodeficiency virus type 1 reduce viral infectivity as well as genomic RNA packaging and dimerization. J Virol. 1997, 71: 3397-3406.PubMed CentralPubMedGoogle Scholar
- Rizvi TA, Kenyon JC, Ali J, Aktar SJ, Phillip PS, Ghazawi A, Mustafa F, Lever AML: Optimal packaging of FIV genomic RNA depends upon a conserved long-range interaction and a palindromic sequence within gag. J Mol Biol. 2010, 403: 103-119. 10.1016/j.jmb.2010.08.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Polge E, Darlix JL, Paoletti J, Fossé P: Characterization of loose and tight dimer forms of avian leukosis virus RNA. J Mol Biol. 2000, 300: 41-56. 10.1006/jmbi.2000.3832.View ArticlePubMedGoogle Scholar
- Houzet L, Paillart JC, Smagulova F, Maurel S, Morichaud Z, Marquet R, Mougel M: HIV controls the selective packaging of genomic, spliced viral and cellular RNAs into virions through different mechanisms. Nucleic Acids Res. 2007, 35: 2695-2704. 10.1093/nar/gkm153.PubMed CentralView ArticlePubMedGoogle Scholar
- Shen N, Jette L, Liang C, Wainberg MA, Laughrea M: Impact of human immunodeficiency virus type 1 RNA dimerization on viral infectivity and of stem-loop B on RNA dimerization and reverse transcription and dissociation of dimerization from packaging. J Virol. 2000, 74: 5729-5735. 10.1128/JVI.74.12.5729-5735.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersen ES, Jeeninga RE, Damgaard CK, Berkhout B, Kjems J: Dimerization and template switching in the 5’ untranslated region between various subtypes of human immunodeficiency virus type 1. J Virol. 2003, 77: 3020-3030. 10.1128/JVI.77.5.3020-3030.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Balakrishnan M, Fay PJ, Bambara RA: The kissing hairpin sequence promotes recombination within the HIV-I 5’ leader region. J Biol Chem. 2001, 276: 36482-36492. 10.1074/jbc.M102860200.View ArticlePubMedGoogle Scholar
- Chin MP, Rhodes TD, Chen J, Fu W, Hu WS: Identification of a major restriction in HIV-1 intersubtype recombination. Proc Natl Acad Sci. 2005, 102: 9002-9007. 10.1073/pnas.0502522102.PubMed CentralView ArticlePubMedGoogle Scholar
- Ennifar E, Paillart JC, Marquet R, Ehresmann B, Ehresmann C, Dumas P, Walter P: HIV-1 RNA dimerization initiation site is structurally similar to the ribosomal A site and binds aminoglycoside antibiotics. J Biol Chem. 2003, 278: 2723-2730. 10.1074/jbc.M205726200.View ArticlePubMedGoogle Scholar
- Ennifar E, Paillart JC, Bodlenner A, Walter P, Weibel JM, Aubertin AM, Pale P, Dumas P, Marquet R: Targeting the dimerization initiation site of HIV-1 RNA with aminoglycosides: from crystal to cell. Nucleic Acids Res. 2006, 34: 2328-2339. 10.1093/nar/gkl317.PubMed CentralView ArticlePubMedGoogle Scholar
- Cardiff RD, Kenney N: Mouse Mammary Tumor Biology: A Short History. Adv Cancer Res. 2007, 98: 53-116. 10.1016/S0065-230X(06)98003-8.View ArticlePubMedGoogle Scholar
- Ross SR: Mouse mammary tumor virus molecular biology and oncogenesis. Viruses. 2010, 2: 2000-2012. 10.3390/v2092000.PubMed CentralView ArticlePubMedGoogle Scholar
- Salmons B, Moritz-Legrand S, Garcha I, Günzburg WH: Construction and characterization of a packaging cell line for MMTV-based conditional retroviral vectors. Biochem Biophys Res Commun. 1989, 159: 1191-1198. 10.1016/0006-291X(89)92236-5.View ArticlePubMedGoogle Scholar
- Rizvi TA, Ali J, Phillip PS, Ghazawi A, Jayanth P, Mustafa F: Role of a heterologous retroviral transport element in the development of genetic complementation assay for mouse mammary tumor virus (MMTV) replication. Virology. 2009, 385: 464-472. 10.1016/j.virol.2008.12.027.View ArticlePubMedGoogle Scholar
- Merino EJ, Wilkinson KA, Coughlan JL, Weeks KM: RNA structure analysis at single nucleotide resolution by selective 2’-hydroxyl acylation and primer extension (SHAPE). J Am Chem Soc. 2005, 127: 4223-4231. 10.1021/ja043822v.View ArticlePubMedGoogle Scholar
- Mortimer SA, Weeks KM: A fast-acting reagent for accurate analysis of RNA secondary and tertiary structure by SHAPE chemistry. J Am Chem Soc. 2007, 129: 4144-4145. 10.1021/ja0704028.View ArticlePubMedGoogle Scholar
- Mortimer SA, Weeks KM: Time-resolved RNA SHAPE chemistry: quantitative RNA structure analysis in one-second snapshots and at single-nucleotide resolution. Nat Protoc. 2009, 4: 1413-1421. 10.1038/nprot.2009.126.View ArticlePubMedGoogle Scholar
- Vasa SM, Guex N, Wilkinson KA, Weeks KM, Giddings MC: ShapeFinder: a software system for high-throughput quantitative analysis of nucleic acid reactivity information resolved by capillary electrophoresis. RNA. 2008, 14: 1979-1990. 10.1261/rna.1166808.PubMed CentralView ArticlePubMedGoogle Scholar
- Laughrea M, Jetté L: Kissing-loop model of HIV-1 genome dimerization: HIV-1 RNAs can assume alternative dimeric forms, and all sequences upstream or downstream of hairpin 248-271 are dispensable for dimer formation. Biochemistry. 1996, 35: 1589-1598. 10.1021/bi951838f.View ArticlePubMedGoogle Scholar
- Mikkelsen JG, Lund AH, Kristensen KD, Duch M, Sørensen MS, Jørgensen P, Pedersen FS: A preferred region for recombinational patch repair in the 5’ untranslated region of primer binding site-impaired murine leukemia virus vectors. J Virol. 1996, 70: 1439-1447.PubMed CentralPubMedGoogle Scholar
- Oh J, McWilliams MJ, Julias JG, Stephen Hughes H: Mutations in the U5 region adjacent to the primer binding site affect tRNA cleavage by human immunodeficiency virus type 1 reverse transcriptase in vivo. J Virol. 2008, 82: 719-727. 10.1128/JVI.02611-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen J, Nikolaitchik O, Singh J, Wright A, Bencsics CE, Coffin JM, Ni N, Lockett S, Pathak VK, Hu WS: High efficiency of HIV-1 genomic RNA packaging and heterozygote formation revealed by single virion analysis. Proc Natl Acad Sci U S A. 2009, 106: 13535-13540. 10.1073/pnas.0906822106.PubMed CentralView ArticlePubMedGoogle Scholar
- Jossinet F, Lodmell JS, Ehresmann C, Ehresmann B, Marquet R: Identification of the in vitro HIV-2/SIV RNA dimerization site reveals striking differences with HIV-1. J Biol Chem. 2001, 276: 5598-5604. 10.1074/jbc.M008642200.View ArticlePubMedGoogle Scholar
- Paillart JC, Skripkin E, Ehresmann B, Ehresmann C, Marquet R: A loop-loop “kissing” complex is the essential part of the dimer linkage of genomic HIV-1 RNA. Proc Natl Acad Sci U S A. 1996, 93: 5572-5577. 10.1073/pnas.93.11.5572.PubMed CentralView ArticlePubMedGoogle Scholar
- Kenyon JC, Ghazawi A, Cheung WK, Phillip PS, Rizvi TA, Lever AM: The secondary structure of the 5’ end of the FIV genome reveals a long-range interaction between R/U5 and gag sequence, and a large, stable stem loop. RNA 2008, 14:2597.,Google Scholar
- Paillart JC, Skripkin E, Ehresmann B, Ehresmann C, Marquet R: In vitro evidence for a long range pseudoknot in the 5’-untranslated and matrix coding regions of HIV-1 genomic RNA. J Biol Chem. 2002, 277: 5995-6004. 10.1074/jbc.M108972200.View ArticlePubMedGoogle Scholar
- Song R, Kafaie J, Laughrea M: Role of the 5’ TAR stem-loop and the U5-AUG duplex in dimerization of HIV-1 genomic RNA. Biochemistry. 2008, 47: 3283-3293. 10.1021/bi7023173.View ArticlePubMedGoogle Scholar
- Spriggs S, Garyu L, Connor R, Summers MF: Potential intra- and intermolecular interactions involving the unique 5’ region of the HIV-1 5’ UTR. Biochemistry. 2008, 47: 13064-13073. 10.1021/bi8014373.PubMed CentralView ArticlePubMedGoogle Scholar
- Baudin F, Marquet R, Isel C, Darlix JL, Ehresmann B, Ehresmann C: Functional sites in the 5’ region of human immunodeficiency virus type 1 RNA form defined structural domains. J Mol Biol. 1993, 229: 382-397. 10.1006/jmbi.1993.1041.View ArticlePubMedGoogle Scholar
- Lever AML: RNA packaging in lentiviruses. Retrovirology 2009, 6:113.,Google Scholar
- Zeffman A, Hassard S, Varani G, Lever A: The major HIV-1 packaging signal is an extended bulged stem loop whose structure is altered on interaction with the Gag polyprotein. J Mol Biol. 2000, 297: 877-893. 10.1006/jmbi.2000.3611.View ArticlePubMedGoogle Scholar
- Nikolaitchik OA, Dilley KA, Fu W, Gorelick RJ, Tai SH, Soheilian F, Ptak RG, Nagashima K, Pathak VK, Hu WS: Dimeric RNA recognition regulates HIV-1 genome packaging. PLoS Pathog 2013, 9:e1003249.,Google Scholar
- Leitner T, Foley B, Hahn B, Marx P, McCutchan F, Mellors J, Wolinsky S, Korber B (Eds): HIV Sequence Compendium 2005. Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, NM, LA-UR 06-0680; 2005.Google Scholar
- Skripkin E, Paillart JC, Marquet R, Ehresmann B, Ehresmann C: Identification of the primary site of the human immunodeficiency virus type 1 RNA dimerization in vitro. Proc Natl Acad Sci U S A. 1994, 91: 4945-4949. 10.1073/pnas.91.11.4945.PubMed CentralView ArticlePubMedGoogle Scholar
- Gavazzi C, Yver M, Isel C, Smyth RP, Rosa-Calatrava M, Lina B, Moulès V, Marquet R: A functional sequence-specific interaction between influenza A virus genomic RNA segments. Proc Natl Acad SciUSA. 2013, 110: 16604-16609. 10.1073/pnas.1314419110.View ArticleGoogle Scholar
- Mathews DH, Sabina J, Zuker M, Turner DH: Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J Mol Biol. 1999, 288: 911-940. 10.1006/jmbi.1999.2700.View ArticlePubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMed CentralView ArticlePubMedGoogle Scholar
- Reuter JS, Mathews DH: RNA structure: software for RNA secondary structure prediction and analysis. BMC Bioinformatics 2010, 11:129.,Google Scholar
- Shackleford GM, Varmus HE: Construction of a clonable, infectious, and tumorigenic mouse mammary tumor virus provirus and a derivative genetic vector. Proc Natl Acad Sci. 1988, 85: 9655-9659. 10.1073/pnas.85.24.9655.PubMed CentralView ArticlePubMedGoogle Scholar
- Goujon M, McWilliam H, Li W, Valentin F, Squizzato S, Paern J, Lopez R: A new bioinformatics analysis tools framework at EMBL-EBI. Nucleic Acids Res. 2010, 38 (Web Server issue): W695-W699. 10.1093/nar/gkq313. doi:10.1093/nar/gkq313. Epub 2010 May 3PubMed CentralView ArticlePubMedGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 2011, 7:539. doi:10.1038/msb.2011.75.,Google Scholar
- Hofacker IL, Fekete M, Stadler PF: Secondary Structure Prediction for Aligned RNA Sequences. J Mol Biol. 2002, 319: 1059-1066. 10.1016/S0022-2836(02)00308-X.View ArticlePubMedGoogle Scholar
- Bernhart SH, Hofacker IL, Will S, Gruber AR, Stadler PF: RNAalifold: improved consensus structure prediction for RNA alignments. BMC Bioinformatics. 2008, 9: 1471-2105. 10.1186/1471-2105-9-474.View ArticleGoogle Scholar
- Hook LM, Agafonova Y, Ross SR, Turner SJ, Golovkina TV: Genetics of mouse mammary tumor virus-induced mammary tumors: linkage of tumor induction to the gag gene. J Virol. 2000, 74: 8876-8883. 10.1128/JVI.74.19.8876-8883.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Cho K, Ferrick DA, Morris DW: Structure and biological activity of the subgenomic Mtv-6 endogenous provirus. Virology. 1995, 206: 395-402. 10.1016/S0042-6822(95)80055-7.View ArticlePubMedGoogle Scholar
- Petropoulos CJ: Retroviral taxonomy, protein structure, sequences, and genetic maps. Retroviruses. Edited by: Coffin JM. 1997, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, 757-Google Scholar
- Moore R, Dixon M, Smith R, Peters G, Dickson C: Complete nucleotide sequence of a milk-transmitted mouse mammary tumor virus: two frameshift suppression events are required for translation of gag and pol. J Virol. 1987, 61: 480-490.PubMed CentralPubMedGoogle Scholar
- Fasel N, Buetti E, Firzlaff J, Pearson K, Diggelmann H: Nucleotide sequence of the 5’ noncoding region and part of the gag gene of mouse mammary tumor virus; identification of the 5’ splicing site for subgenomic mRNAs. Nucleic Acids Res. 1983, 11: 6943-6955. 10.1093/nar/11.20.6943.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh JL, Erfle M, Wykes EJ: The pIC plasmid and phage vectors with versatile cloning sites for recombinant selection by insertional inactivation. Gene. 1984, 32: 481-485. 10.1016/0378-1119(84)90022-2.View ArticlePubMedGoogle Scholar
- Gibbs JS, Regier DA, Desrosiers RC: Construction and in vitro properties of SIV mac mutants with deletions in “nonessential” genes. AIDS Res Hum Retrovirus. 1994, 10: 607-616. 10.1089/aid.1994.10.607.View ArticleGoogle Scholar
- Marquet R, Baudin F, Gabus C, Darlix JL, Mougel M, Ehresmann C, Ehresmann B: Dimerization of human immunodeficiency virus (type 1) RNA: stimulation by cations and possible mechanism. Nucleic Acids Res. 1991, 9: 2349-2357. 10.1093/nar/19.9.2349.View ArticleGoogle Scholar
- Wilkinson KA, Gorelick RJ, Vasa SM, Guex N, Rein A, Mathews DH, Giddings MC, Weeks KM: High-throughput SHAPE analysis reveals structures in HIV-1 genomic RNA strongly conserved across distinct biological states. PLoS Biol 2008, 6:e96. doi:10.1371/journal.pbio.0060096.,Google Scholar
- Naldini L, Blömer U, Gallay P, Ory D, Mulligan R, Gage FH, Verma IM, Trono D: In vivo gene delivery and stable transduction of nondividing cells by a lentiviral vector. Science. 1996, 272: 263-267. 10.1126/science.272.5259.263.View ArticlePubMedGoogle Scholar
- Ghazawi A, Mustafa F, Phillip PS, Jayanth P, Ali J, Rizvi TA: Both the 5’ and 3’ LTRs of FIV contain minor RNA encapsidation determinants compared to the two core packaging determinants within the 5’ untranslated region and gag. Microbes Infect. 2006, 8: 767-778. 10.1016/j.micinf.2005.09.015.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.