Evolution of the HIV-1 Rev Response Element during Natural Infection Reveals Nucleotide Changes That Correlate with Altered Structure and Increased Activity over Time.

The HIV-1 Rev response element (RRE) is a cis-acting RNA element characterized by multiple stem-loops. Binding and multimerization of the HIV Rev protein on the RRE promote the nucleocytoplasmic export of incompletely spliced mRNAs, an essential step in HIV replication. Most of our understanding of the Rev-RRE regulatory axis comes from studies of lab-adapted HIV clones. However, in human infection, HIV evolves rapidly, and mechanistic studies of naturally occurring Rev and RRE sequences are essential to understanding this system. We previously described the functional activity of two RREs found in circulating viruses in a patient followed during the course of HIV infection. The early RRE was less functionally active than the late RRE, despite differing in sequence by only 4 nucleotides. In this study, we describe the sequence, function, and structural evolution of circulating RREs in this patient using plasma samples collected over 6 years of untreated infection. RRE sequence diversity varied over the course of infection, with evidence of selection pressure that led to sequence convergence as disease progressed being found. An increase in RRE functional activity was observed over time, and a key mutation was identified that correlates with a major conformational change in the RRE and increased functional activity. Additional mutations were found that may have contributed to increased activity as a result of greater Shannon entropy in RRE stem-loop II, which is key to primary Rev binding.IMPORTANCE HIV-1 replication requires interaction of the viral Rev protein with a cis-acting regulatory RNA, the Rev response element (RRE), whose sequence changes over time during infection within a single host. In this study, we show that the RRE is subject to selection pressure and that RREs from later time points in infection tend to have higher functional activity. Differences in RRE functional activity are attributable to specific changes in RNA structure. Our results suggest that RRE evolution during infection may be important for HIV pathogenesis and that efforts to develop therapies acting on this viral pathway should take this into account.

HIV-1 Rev interacts with HERV-K RcREs present in the human genome and promotes export of unspliced HERV-K proviral RNA.

BACKGROUND: The HERV-K (HML-2) viruses are the youngest of the human endogenous retroviruses. They are present as several almost complete proviral copies and numerous fragments in the human genome. Many HERV-K proviruses express a regulatory protein Rec, which binds to an element present in HERV-K mRNAs called the RcRE. This interaction is necessary for the nucleo-cytoplasmic export and expression of HERV-K mRNAs that retain introns and plays a role analogous to that of Rev and the RRE in HIV replication. There are over 900 HERV-K RcREs distributed throughout the human genome. Thus, it was of interest to determine if Rev could functionally interact with selected RcRE elements that map either to HERV-K proviruses or human gene regions. This interaction would have the potential to alter the expression of both HERV-K mRNAs and cellular mRNAs during HIV-1 infection. RESULTS: In this study we employed a combination of RNAseq, bioinformatics and cell-based functional assays. Potential RcREs were identified through a number of bioinformatic approaches. They were then tested for their ability to promote export and translation of a reporter mRNA with a retained intron in conjunction with Rev or Rec. Some of the selected elements functioned well with either Rev, Rec or both, whereas some showed little or no function. Rev function on individual RcREs varied and was also dependent on the Rev sequence. We also performed RNAseq on total and cytoplasmic RNA isolated from SupT1 cells expressing HIV Rev, with or without Tat, or HERV-K Rec. Proviral mRNA from three HERV-K loci (4p16.1b, 22q11.23 and most significantly 3q12.3) accumulated in the cytoplasm in the presence of Rev or Tat and Rev, but not Rec. Consistent with this, the 3' RcRE from 3q12.3 functioned well with HIV-Rev in our reporter assay. In contrast, this RcRE showed little or no function with Rec. CONCLUSIONS: The HIV Rev protein can functionally interact with many RcREs present in the human genome, depending on the RcRE sequence, as well as the Rev sequence. This leads to export of some of the HERV-K proviral mRNAs and also has the potential to change the expression of non-viral genes.

Characterization of APOBEC3 variation in a population of HIV-1 infected individuals in northern South Africa.

BACKGROUND: The apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3 (APOBEC3) genes A3D, A3F, A3G and A3H have all been implicated in the restriction of human immunodeficiency virus type 1 (HIV-1) replication. Polymorphisms in these genes are likely to impact viral replication and fitness, contributing to viral diversity. Currently, only a few studies indicate that polymorphisms in the A3 genes may be correlated with infection risk and disease progression. METHODS: To characterize polymorphisms in the coding regions of these APOBEC3 genes in an HIV-1 infected population from the Limpopo Province of South Africa, APOBEC3 gene fragments were amplified from genomic DNA of 192 HIV-1 infected subjects and sequenced on an Illumina MiSeq platform. SNPs were confirmed and compared to SNPs in other populations reported in the 1000 Genome Phase III and HapMap databases, as well as in the ExAC exome database. Hardy-Weinberg Equilibrium was calculated and haplotypes were inferred using the LDlink 3.0 web tool. Linkage Disequilibrium (LD) for these SNPS were calculated in the total 1000 genome and AFR populations using the same tool. RESULTS: Known variants compared to the GRCh37 consensus genome sequence were detected at relatively high frequencies (> 5%) in all of the APOBEC3 genes. A3H showed the most variation, with several of the variants present in both alleles in almost all of the patients. Several minor allele variants (< 5%) were also detected in A3D, A3F and A3G. In addition, novel R6K, L221R and T238I variants in A3D and I117I in A3F were observed. Four, five, four, and three haplotypes were identified for A3D, A3F, A3G, and A3H respectively. CONCLUSIONS: The study showed significant polymorphisms in the APOBEC3D, 3F, 3G and 3H genes in our South African HIV1-infected cohort. In the case of all of these genes, the polymorphisms were generally present at higher frequencies than reported in other 1000 genome populations and in the ExAC exome consortium database .

Evolutionary conservation of a molecular machinery for export and expression of mRNAs with retained introns.

Intron retention is one of the least studied forms of alternative splicing. Through the use of retrovirus and other model systems, it was established many years ago that mRNAs with retained introns are subject to restriction both at the level of nucleocytoplasmic export and cytoplasmic expression. It was also demonstrated that specific cis-acting elements in the mRNA could serve to bypass these restrictions. Here we show that one of these elements, the constitutive transport element (CTE), first identified in the retrovirus MPMV and subsequently in the human NXF1 gene, is a highly conserved element. Using GERP analysis, CTEs with strong primary sequence homology, predicted to display identical secondary structure, were identified in NXF genes from >30 mammalian species. CTEs were also identified in the predicted NXF1 genes of zebrafish and coelacanths. The CTE from the zebrafish NXF1 was shown to function efficiently to achieve expression of mRNA with a retained intron in human cells in conjunction with zebrafish Nxf1 and cofactor Nxt proteins. This demonstrates that all essential functional components for expression of mRNA with retained introns have been conserved from fish to man.

High level of HIV-1 drug resistance mutations in patients with unsuppressed viral loads in rural northern South Africa.

BACKGROUND: Combination antiretroviral therapy (cART) has significantly reduced HIV morbidity and mortality in both developed and developing countries. However, the sustainability of cART may be compromised by the emergence of viral drug resistance mutations (DRM) and the cellular persistence of proviruses carrying these DRM. This is potentially a more serious problem in resource limited settings. METHODS: DRM were evaluated in individuals with unsuppressed viral loads after first or multiple lines of cART at two sites in rural Limpopo, South Africa. Seventy-two patients with viral loads of >1000 copies/ml were recruited between March 2014 and December 2015. Complete protease (PR) and partial Reverse Transcriptase (RT) sequences were amplified from both plasma RNA and paired proviral DNA from 35 of these subjects. Amplicons were directly sequenced to determine subtype and DRM using the Stanford HIV Drug Resistance Interpretation algorithm. RESULTS: Among the 72 samples, 69 could be PCR amplified from RNA and 35 from both RNA and DNA. Sixty-five (94.2%) viruses were subtype C, while one was subtype B (1.4%), one recombinant K/C, one recombinant C/B and one unclassified. Fifty-eight (84%) sequences carried at least one DRM, while 11 (15.9%) displayed no DRM. DRM prevalence according to drug class was: NRTI 60.8% NNRTI 65.2%, and PI 5.8%. The most common DRMs were; M184V (51.7%), K103N (50%), V106M (20.6%), D67N (13.3%), K65R (12%). The frequency of the DRM tracked well with the frequency of use of medications to which the mutations were predicted to confer resistance. Interestingly, a significant number of subjects showed predicted resistance to the newer NNRTIs, etravirine (33%) and rilpivirine (42%), both of which are not yet available in this setting. The proportion of DRM in RNA and DNA were mostly similar with the exception of the thymidine analogue mutations (TAMs) D67N, K70R, K219QE; and K103N which were slightly more prevalent in DNA than RNA. Subjects who had received cART for at least 5 years were more likely to harbour >2 DRM (p < 0.05) compared to those treated for a shorter period. DRM were more prevalent in this rural setting compared to a neighbouring urban setting. CONCLUSION: We found a very high prevalence of NRTI and NNRTI DRM in patients from rural Limpopo settings with different durations of treatment. The prevalence was significantly higher than those reported in urban settings in South Africa. The dominance of NNRTI based mutations late in treatment supports the use of PI based regimens for second line treatment in this setting. The slight dominance of TAMs in DNA from infected PBMCs compared to plasma virus requires further studies that should include cART subjects with suppressed virus. Such studies will improve our understanding of the pattern of drug resistance and dynamics of viral persistence in these rural settings.

The HIV-1 Rev response element (RRE) adopts alternative conformations that promote different rates of virus replication.

The HIV Rev protein forms a complex with a 351 nucleotide sequence present in unspliced and incompletely spliced human immunodeficiency virus (HIV) mRNAs, the Rev response element (RRE), to recruit the cellular nuclear export receptor Crm1 and Ran-GTP. This complex facilitates nucleo-cytoplasmic export of these mRNAs. The precise secondary structure of the HIV-1 RRE has been controversial, since studies have reported alternative structures comprising either four or five stem-loops. The published structures differ only in regions that lie outside of the primary Rev binding site. Using in-gel SHAPE, we have now determined that the wt NL4-3 RRE exists as a mixture of both structures. To assess functional differences between these RRE 'conformers', we created conformationally locked mutants by site-directed mutagenesis. Using subgenomic reporters, as well as HIV replication assays, we demonstrate that the five stem-loop form of the RRE promotes greater functional Rev/RRE activity compared to the four stem-loop counterpart.

Characterization and in vitro activity of a branched peptide boronic acid that interacts with HIV-1 RRE RNA.

A branched peptide containing multiple boronic acids was found to bind RRE IIB selectively and inhibit HIV-1 p24 capsid production in a dose-dependent manner. Structure-activity relationship studies revealed that branching in the peptide is crucial for the low micromolar binding towards RRE IIB, and the peptide demonstrates selectivity towards RRE IIB in the presence of tRNA. Footprinting studies suggest a binding site on the upper stem and internal loop regions of the RNA, which induces enzymatic cleavage of the internal loops of RRE IIB upon binding.

Limited nucleotide changes in the Rev response element (RRE) during HIV-1 infection alter overall Rev-RRE activity and Rev multimerization.

HIV-1 Rev and the Rev response element (RRE) enable a critical step in the viral replication cycle by facilitating the nuclear export of intron-containing mRNAs, yet their activities have rarely been analyzed in natural infections. This study characterized their genetic and functional variation in a small cohort of HIV-infected individuals. Multiple Rev and RRE sequences were obtained using single-genome sequencing (SGS) of plasma samples collected within 6 months after seroconversion and at a later time. This allowed the identification of cognate sequences that were linked in vivo in the same viral genome and acted together as a functional unit. Phylogenetic analyses of these sequences indicated that 4/5 infections were founded by a single transmission event. Rev and RRE variants from each time point were subjected to functional analysis as both cognate pairs and as individual components. While a range of Rev-RRE activities were seen, the activity of cognate pairs from a single time point clustered to a discrete level, which was termed the set point. In 3/5 patients, this set point changed significantly over the time period studied. In all patients, RRE activity was more sensitive to sequence variation than Rev activity and acted as the primary driver of the cognate set point. Selected patient RREs were also shown to have differences in Rev multimerization using gel shift binding assays. Thus, rather than acting as a simple on-off switch or maintaining a constant level of activity throughout infection, the Rev-RRE system can fluctuate, presumably to control replication.

Measurement of HIV Rev-Rev Response Element Functional Activity.

Retroviruses must overcome cellular restrictions to the nucleocytoplasmic export of viral mRNAs that retain introns in order to complete their replication cycle. HIV accomplishes this using a system comprised of a trans-acting viral protein, Rev, and a cis-acting RNA secondary structure in the viral genome, the Rev-Response Element (RRE). HIV primary isolates differ with respect to the sequence and functional activity of the Rev-RRE system. Here, we describe a high throughput assay system for analyzing Rev-RRE functional activity using packageable viral vectors.

Xenotropic MLV envelope proteins induce tumor cells to secrete factors that promote the formation of immature blood vessels.

BACKGROUND: Xenotropic Murine leukemia virus-Related Virus (XMRV) is a γ-retrovirus initially reported to be present within familial human prostate tumors and the blood of patients with chronic fatigue syndrome. Subsequent studies however were unable to replicate these findings, and there is now compelling evidence that the virus evolved through rare retroviral recombination events in human tumor cell lines established through murine xenograft experiments. There is also no direct evidence that XMRV infection has any functional effects that contribute to tumor pathogenesis. RESULTS: Herein we describe an additional xenotropic MLV, "B4rv", found in a cell line derived from xenograft experiments with the human prostate cancer LNCaP cell line. When injected subcutaneously in nude mice, LNCaP cells infected with XMRV or B4rv formed larger tumors that were highly hemorrhagic and displayed poor pericyte/smooth muscle cell (SMC) investment, markers of increased metastatic potential. Conditioned media derived from XMRV- or B4rv-infected LNCaPs, but not an amphotropic MLV control virus infected LNCaPs, profoundly decreased expression of marker genes in cultured SMC, consistent with inhibition of SMC differentiation/maturation. Similar effects were seen with a chimeric virus of the amphotropic MLV control virus containing the XMRV env gene, but not with an XMRV chimeric virus containing the amphotropic MLV env gene. UV-inactivated XMRV and pseudovirions that were pseudotyped with XMRV envelope protein also produce conditioned media that down-regulated SMC marker gene expression in vitro. CONCLUSIONS: Together these results indicate that xenotropic MLV envelope proteins are sufficient to induce the production of factors by tumor cells that suppress vascular SMC differentiation, providing evidence for a novel mechanism by which xenotropic MLVs might alter tumor pathogenesis by disrupting tumor vascular maturation. Although it is highly unlikely that either XMRV or B4Rv themselves infect humans and are pathogenic, the results suggest that xenograft approaches commonly used in the study of human cancer promote the evolution of novel retroviruses with pathogenic properties.

The Tpr protein regulates export of mRNAs with retained introns that traffic through the Nxf1 pathway.

Post-transcriptional regulation of mRNA includes restriction mechanisms to prevent export and expression of mRNAs that are incompletely spliced. Here we present evidence that the mammalian protein Tpr is involved in this restriction. To study the role of Tpr in export of mRNA with retained introns, we used reporters in which the mRNA was exported either via the Nxf1/Nxt1 pathway using a CTE or via the Crm1 pathway using Rev/RRE. Our data show that even modest knockdown of Tpr using RNAi leads to a significant increase in export and translation from the mRNA containing the CTE. In contrast, Tpr perturbation has no effect on export of mRNA containing the RRE, either in the absence or presence of Rev. Also, no effects were observed on export of a completely spliced mRNA. Taken together, our results indicate that Tpr plays an important role in quality control of mRNA trafficked on the Nxf1 pathway.

Rev-Rev Response Element Activity Selection Bias at the HIV Transmission Bottleneck.

HIV is not efficiently transmitted between hosts, and selection of viral variants occurs during the process of sexual transmission. The factors that confer selective advantage at the transmission bottleneck remain incompletely understood. We explored whether differences in the Rev-Rev Response Element (RRE) regulatory axis of HIV affect transmission fitness, since functional variation in the Rev-RRE axis in different viral isolates has been shown to affect replication kinetics and relative expression of many HIV proteins. Single genome HIV sequences were identified from nine linked subject pairs near the time of female-to-male transmission. Using a rapid flow-cytometric assay, we found that the functional Rev-RRE activity varied significantly between isolates. Moreover, it was generally lower in recipients' viruses compared to the corresponding donor viruses. In six of nine transmission events, recipient virus Rev-RRE activity clustered at the extreme low end of the range of donor virus activity. Rev-RRE pair activity was an unpredictable product of component Rev and RRE activity variation. These data indicate selection pressure on the Rev-RRE axis during female-to-male sexual transmission. Variation in the activity of the Rev-RRE axis may permit viral adaptation to different fitness landscapes and could play an important role in HIV pathogenesis.

Rev-Rev Response Element Activity Selection Bias at the Human Immunodeficiency Virus Transmission Bottleneck.

Background: Sexual transmission of human immunodeficiency virus (HIV) is inefficient and results in selection of viral variants based on incompletely understood factors. Functional variation in the Rev-Rev response element (RRE) regulatory axis of HIV affect replication kinetics and relative expression of viral proteins. We explored whether differences in this axis among viral isolates affect transmission fitness. Methods: HIV sequences were identified from nine female-to-male transmission pairs. Using a rapid flow cytometric assay, we analyzed Rev-RRE functional activity of primary isolates. Results: Rev-RRE activity was significantly lower in recipient viruses compared with corresponding donor viruses. In most transmission events, recipient virus Rev-RRE activity clustered at the extreme low end of the range of donor virus activity. Conclusions: These data indicate selection pressure on the Rev-RRE axis during female-to-male sexual transmission. Variation in Rev-RRE activity may permit viral adaptation to different fitness landscapes and could play an important role in HIV pathogenesis.

Single-nucleotide changes in the HIV Rev-response element mediate resistance to compounds that inhibit Rev function.

Previously we described the identification of two compounds (3-amino-5-ethyl-4,6-dimethylthieno[2,3-b]pyridine-2-carboxamide [103833] and 4-amino-6-methoxy-2-(trifluoromethyl)-3-quinolinecarbonitrile [104366]) that interfered with HIV replication through the inhibition of Rev function. We now describe resistant viral variants that arose after drug selection, using virus derived from two different HIV proviral clones, NL4-3 and R7/3. With HIV(NL4-3), each compound selected a different single point mutation in the Rev response element (RRE) at the bottom of stem-loop IIC. Either mutation led to the lengthening of the stem-loop IIC stem by an additional base pair, creating an RRE that was more responsive to lower concentrations of Rev than the wild type. Surprisingly, wild-type HIV(R7/3) was also found to be inhibited when tested with these compounds, in spite of the fact this virus already has an RNA stem-loop IIC similar to the one in the resistant NL4-3 variant. When drug resistance was selected in HIV(R7/3), a virus arose with two nucleotide changes that mapped to the envelope region outside the RRE. One of these nucleotide changes was synonymous with respect to env, and one was not. The combination of both nucleotide changes appeared to be necessary for the resistance phenotype as the individual point mutations by themselves did not convey resistance. Thus, although drug-resistant variants can be generated with both viral strains, the underlying mechanism is clearly different. These results highlight that minor nucleotide changes in HIV RNA, outside the primary Rev binding site, can significantly alter the efficiency of the Rev/RRE pathway.

Molecular epidemiology of HIV in two highly endemic areas of northeastern South Africa.

There is paucity of data on the genetic landscape of HIV-1 viruses circulating in the Limpopo Province of northeastern South Africa. Here, we examine the genetic diversity of viruses from Bela-Bela and Musina, two towns with high HIV prevalence. Between June 2007 and March 2008, blood samples were collected from antiretroviral-drug-naïve individuals. Viruses were analyzed for genetic subtypes and drug resistance mutations. All of the viruses in these samples were shown by phylogenetic analysis based on gag p17, gag p24, reverse transcriptase, protease and envelope C2-C3 gene regions to belong to HIV-1 subtype C. Two of 44 reverse transcriptase sequences (4.5%) contained N rather than the consensus K at position 103. The K103N mutation is normally associated with resistance to NNRTIs. No major mutations were observed in the protease gene. However, several polymorphisms and amino acid changes normally considered to be minor drug resistance mutations were observed in the protease sequences. These results suggest that HIV-1 subtype C remains the predominant variant responsible for the epidemic in northeastern South Africa and that the prevalence of drug-resistant viruses among the naïve population is low.

An intron with a constitutive transport element is retained in a Tap messenger RNA.

Alternative splicing is a key factor contributing to genetic diversity and evolution. Intron retention, one form of alternative splicing, is common in plants but rare in higher eukaryotes, because messenger RNAs with retained introns are subject to cellular restriction at the level of cytoplasmic export and expression. Often, retention of internal introns restricts the export of these mRNAs and makes them the targets for degradation by the cellular nonsense-mediated decay machinery if they contain premature stop codons. In fact, many of the database entries for complementary DNAs with retained introns represent them as artefacts that would not affect the proteome. Retroviruses are important model systems in studies of regulation of RNAs with retained introns, because their genomic and mRNAs contain one or more unspliced introns. For example, Mason-Pfizer monkey virus overcomes cellular restrictions by using a cis-acting RNA element known as the constitutive transport element (CTE). The CTE interacts directly with the Tap protein (also known as nuclear RNA export factor 1, encoded by NXF1), which is thought to be a principal export receptor for cellular mRNA, leading to the hypothesis that cellular mRNAs with retained introns use cellular CTE equivalents to overcome restrictions to their expression. Here we show that the Tap gene contains a functional CTE in its alternatively spliced intron 10. Tap mRNA containing this intron is exported to the cytoplasm and is present in polyribosomes. A small Tap protein is encoded by this mRNA and can be detected in human and monkey cells. Our results indicate that Tap regulates expression of its own intron-containing RNA through a CTE-mediated mechanism. Thus, CTEs are likely to be important elements that facilitate efficient expression of mammalian mRNAs with retained introns.

HIV-1 Rev-RRE functional activity in primary isolates is highly dependent on minimal context-dependent changes in Rev.

During HIV infection, intron-containing viral mRNAs are exported from the cell nucleus to the cytoplasm to complete the replication cycle. Cellular restrictions on the export of incompletely spliced transcripts are overcome by a viral protein, Rev, and an RNA structure found in all unspliced and incompletely spliced viral mRNAs, the Rev Response Element (RRE). Primary HIV isolates display substantial variation in the sequence and functional activity of Rev proteins. We analyzed Rev from two primary isolates with disparate activity that resulted in differences in in vitro fitness of replication-competent viral constructs. The results showed that amino acid differences within the oligomerization domain, but not the arginine-rich motif or the nuclear export signal, determined the level of Rev activity. Two specific amino acid substitutions were sufficient to alter the low-activity Rev to a high-activity phenotype. Other mutations in Rev sequences had unpredictable effects on activity that differed between the two Rev backbones. The sensitivity of Rev function level to small sequence changes likely permits modulation of Rev-RRE activity during HIV infection, which may play a role in pathogenesis. The functional consequences of Rev mutations differed between primary isolates, highlighting the challenge of generalizing studies of Rev conducted using laboratory HIV strains.

Probing the HIV-1 genomic RNA trafficking pathway and dimerization by genetic recombination and single virion analyses.

Once transcribed, the nascent full-length RNA of HIV-1 must travel to the appropriate host cell sites to be translated or to find a partner RNA for copackaging to form newly generated viruses. In this report, we sought to delineate the location where HIV-1 RNA initiates dimerization and the influence of the RNA transport pathway used by the virus on downstream events essential to viral replication. Using a cell-fusion-dependent recombination assay, we demonstrate that the two RNAs destined for copackaging into the same virion select each other mostly within the cytoplasm. Moreover, by manipulating the RNA export element in the viral genome, we show that the export pathway taken is important for the ability of RNA molecules derived from two viruses to interact and be copackaged. These results further illustrate that at the point of dimerization the two main cellular export pathways are partially distinct. Lastly, by providing Gag in trans, we have demonstrated that Gag is able to package RNA from either export pathway, irrespective of the transport pathway used by the gag mRNA. These findings provide unique insights into the process of RNA export in general, and more specifically, of HIV-1 genomic RNA trafficking.

Resistance to RevM10 inhibition reflects a conformational switch in the HIV-1 Rev response element.

Nuclear export of certain HIV-1 mRNAs requires an interaction between the viral Rev protein and the Rev response element (RRE), a structured element located in the Env region of its RNA genome. This interaction is an attractive target for both drug design and gene therapy, exemplified by RevM10, a transdominant negative protein that, when introduced into host cells, disrupts viral mRNA export. However, two silent G->A mutations in the RRE (RRE61) confer RevM10 resistance, which prompted us to examine RRE structure using a novel chemical probing strategy. Variations in region III/IV/V of mutant RNAs suggest a stepwise rearrangement to RevM10 resistance. Mass spectrometry was used to directly assess Rev "loading" onto RRE and its variants, indicating that this is unaffected by RNA structural changes. Similarity in chemical footprints with mutant protein implicates additional host factors in RevM10 resistance.

Intron retention and its impact on gene expression and protein diversity: A review and a practical guide.

Intron retention (IR) occurs when a complete and unspliced intron remains in mature mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions, including several that impact human health and disease. Although experimental technologies used to study other forms of mRNA splicing can also be used to investigate IR, a specialized downstream computational analysis is optimal for IR discovery and analysis. Here we provide a review of IR and its biological implications, as well as a practical guide for how to detect and analyze it. Several methods, including long read third generation direct RNA sequencing, are described. We have developed an R package, FakIR, to facilitate the execution of the bioinformatic tasks recommended in this review and a tutorial on how to fit them to users aims. Additionally, we provide guidelines and experimental protocols to validate IR discovery and to evaluate the potential impact of IR on gene expression and protein output. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Regulation/Alternative Splicing RNA Methods > RNA Analyses in vitro and In Silico.