The HIV-1 accessory protein Vif is essential for viral replication by counteracting the host restriction factor APOBEC3G (A3G), and balanced levels of both proteins are required for efficient viral replication. Noncoding exons 2/2b contain the Vif start codon between their alternatively used splice donors 2 and 2b (D2 and D2b). For vif mRNA, intron 1 must be removed while intron 2 must be retained. Thus, splice acceptor 1 (A1) must be activated by U1 snRNP binding to either D2 or D2b, while splicing at D2 or D2b must be prevented. Here, we unravel the complex interactions between previously known and novel components of the splicing regulatory network regulating HIV-1 exon 2/2b inclusion in viral mRNAs. In particular, using RNA pulldown experiments and mass spectrometry analysis, we found members of the heterogeneous nuclear ribonucleoparticle (hnRNP) A/B family binding to a novel splicing regulatory element (SRE), the exonic splicing silencer ESS2b, and the splicing regulatory proteins Tra2/SRSF10 binding to the nearby exonic splicing enhancer ESE2b. Using a minigene reporter, we performed bioinformatics HEXplorerguided mutational analysis to narrow down SRE motifs affecting splice site selection between D2 and D2b. Eventually, the impacts of these SREs on the viral splicing pattern and protein expression were exhaustively analyzed in viral particle production and replication experiments. Masking of these protein binding sites by use of locked nucleic acids (LNAs) impaired Vif expression and viral replication.IMPORTANCE Based on our results, we propose a model in which a dense network of SREs regulates vif mRNA and protein expression, crucial to maintain viral replication within host cells with varying A3G levels and at different stages of infection. This regulation is maintained by several serine/arginine-rich splicing factors (SRSF) and hnRNPs binding to those elements. Targeting this cluster of SREs with LNAs may lead to the development of novel effective therapeutic strategies.KEYWORDS HEXplorer score, exon recognition, host restriction factor, human immunodeficiency virus, pre-mRNA processing, splicing regulatory elements D uring long terminal repeat (LTR)-driven transcription, over 50 mRNA isoforms emerge by alternative splicing of the HIV-1 precursor mRNA (1, 2). According to their distinct sizes, mRNA isoforms can be divided into three different classes: 2-kb mRNAs (intronless), encoding Tat, Rev, and Nef; intron-containing 4-kb mRNAs, encoding Vif, Vpr, Vpu, and Env; and 9-kb unspliced mRNAs, encoding Gag and Gag-Pol (3). Viral gene expression follows a strict chronological order (4-6). In the early phase, only intronless mRNAs are transported out of the nucleus and translated, whereas intron-

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5’ end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bβ-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3’ exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3’ exon end definition.

Most human pathogenic mutations in 5′ splice sites affect the canonical GT in positions +1 and +2, leading to noncanonical dinucleotides. On the other hand, noncanonical dinucleotides are observed under physiological conditions in ∼1% of all human 5′ss. It is therefore a challenging task to understand the pathogenic mutation mechanisms underlying the conditions under which noncanonical 5′ss are used. In this work, we systematically examined noncanonical 5′ splice site selection, both experimentally using splicing competition reporters and by analyzing a large RNA-seq data set of 54 fibroblast samples from 27 subjects containing a total of 2.4 billion gapped reads covering 269,375 exon junctions. From both approaches, we consistently derived a noncanonical 5′ss usage ranking GC > TT > AT > GA > GG > CT. In our competition splicing reporter assay, noncanonical splicing was strictly dependent on the presence of upstream or downstream splicing regulatory elements (SREs), and changes in SREs could be compensated by variation of U1 snRNA complementarity in the competing 5′ss. In particular, we could confirm splicing at different positions (i.e., −1, +1, +5) of a splice site for all noncanonical dinucleotides “weaker” than GC. In our comprehensive RNA-seq data set analysis, noncanonical 5′ss were preferentially detected in weakly used exon junctions of highly expressed genes. Among high-confidence splice sites, they were 10-fold overrepresented in clusters with a neighboring, more frequently used 5′ss. Conversely, these more frequently used neighbors contained only the dinucleotides GT, GC, and TT, in accordance with the above ranking.

No abstract