Correction: An activator of G protein-coupled receptor and MEK1/2-ERK1/2 signaling inhibits HIV-1 replication by altering viral RNA processing.

[This corrects the article DOI: 10.1371/journal.ppat.1008307.].

The global Protein-RNA interaction map of ESRP1 defines a post-transcriptional program that is essential for epithelial cell function.

The epithelial splicing regulatory proteins, ESRP1 and ESRP2, are essential for mammalian development through the regulation of a global program of alternative splicing of genes involved in the maintenance of epithelial cell function. To further inform our understanding of the molecular functions of ESRP1, we performed enhanced crosslinking immunoprecipitation coupled with high-throughput sequencing (eCLIP) in epithelial cells of mouse epidermis. The genome-wide binding sites of ESRP1 were integrated with RNA-Seq analysis of alterations in splicing and total gene expression that result from epidermal ablation of Esrp1 and Esrp2. These studies demonstrated that ESRP1 functions in splicing regulation occur primarily through direct binding in a position-dependent manner to promote either exon inclusion or skipping. In addition, we also identified widespread binding of ESRP1 in 3' and 5' untranslated regions (UTRs) of genes involved in epithelial cell function, suggesting that its post-transcriptional functions extend beyond splicing regulation.

The mAB 13A4 monoclonal antibody to the mouse PROM1 protein recognizes a structural epitope.

PROM1 (CD133, AC133) is a protein that is required for the maintenance of primary cilia. Mutation in the Prom1 gene in humans and animal models are associated with several forms of retinal degeneration. mAB 13A4 is the main reagent used to detect the mouse PROM1 protein. We endeavored to map the epitope of the rat monoclonal antibody mAB 13A4 to the mouse PROM1 protein. Deletion mutagenesis demonstrated that mAB 13A4 recognizes a structural epitope that is stabilized by two of the extracellular domains of PROM1. Furthermore, the affinity of mAB 13A4 to the major PROM1 isoform in photoreceptor cells is significantly reduced due to the inclusion of a photoreceptor-specific alternative exon in the third extracellular domain. Interestingly, a deletion in the photoreceptor specific isoform of six amino acids adjacent to the alternative exon restored the affinity of mAB 13A4 to PROM1. The results of the mutagenesis are consistent with the computationally predicted helical bundle structure of PROM1 and point to the utility of mAB 13A4 for evaluating the effect of mutations on the PROM1 structure. Our results show that the PROM1 isoform composition needs to be considered when interpreting tissue and developmental expression data produced by mAB 13A4.

The Musashi proteins direct post-transcriptional control of protein expression and alternate exon splicing in vertebrate photoreceptors.

The Musashi proteins, MSI1 and MSI2, are conserved RNA binding proteins with a role in the maintenance and renewal of stem cells. Contrasting with this role, terminally differentiated photoreceptor cells express high levels of MSI1 and MSI2, pointing to a role for the two proteins in vision. Combined knockout of Msi1 and Msi2 in mature photoreceptor cells abrogated the retinal response to light and caused photoreceptor cell death. In photoreceptor cells the Musashi proteins perform distinct nuclear and cytoplasmic functions. In the nucleus, the Musashi proteins promote splicing of photoreceptor-specific alternative exons. Surprisingly, conserved photoreceptor-specific alternative exons in genes critical for vision proved to be dispensable, raising questions about the selective pressures that lead to their conservation. In the cytoplasm MSI1 and MSI2 activate protein expression. Loss of Msi1 and Msi2 lead to reduction in the levels of multiple proteins including proteins required for vision and photoreceptor survival. The requirement for MSI1 and MSI2 in terminally differentiated photoreceptors alongside their role in stem cells shows that, depending on cellular context, these two proteins can control processes ranging from cell proliferation to sensory perception.

SARS-CoV-2 Delta variant induces enhanced pathology and inflammatory responses in K18-hACE2 mice.

The COVID-19 pandemic has been fueled by SARS-CoV-2 novel variants of concern (VOC) that have increased transmissibility, receptor binding affinity, and other properties that enhance disease. The goal of this study is to characterize unique pathogenesis of the Delta VOC strain in the K18-hACE2-mouse challenge model. Challenge studies suggested that the lethal dose of Delta was higher than Alpha or Beta strains. To characterize the differences in the Delta strain's pathogenesis, a time-course experiment was performed to evaluate the overall host response to Alpha or Delta variant challenge. qRT-PCR analysis of Alpha- or Delta-challenged mice revealed no significant difference between viral RNA burden in the lung, nasal wash or brain. However, histopathological analysis revealed high lung tissue inflammation and cell infiltration following Delta- but not Alpha-challenge at day 6. Additionally, pro-inflammatory cytokines were highest at day 6 in Delta-challenged mice suggesting enhanced pneumonia. Total RNA-sequencing analysis of lungs comparing challenged to no challenge mice revealed that Alpha-challenged mice have more total genes differentially activated. Conversely, Delta-challenged mice have a higher magnitude of differential gene expression. Delta-challenged mice have increased interferon-dependent gene expression and IFN-γ production compared to Alpha. Analysis of TCR clonotypes suggested that Delta challenged mice have increased T-cell infiltration compared to Alpha challenged. Our data suggest that Delta has evolved to engage interferon responses in a manner that may enhance pathogenesis. The in vivo and in silico observations of this study underscore the need to conduct experiments with VOC strains to best model COVID-19 when evaluating therapeutics and vaccines.

Transcriptional profiling reveals roles of intercellular Fgf9 signaling in astrocyte maturation and synaptic refinement during brainstem development.

Neural tissue maturation is a coordinated process under tight transcriptional control. We previously analyzed the kinetics of gene expression in the medial nucleus of the trapezoid body (MNTB) in the brainstem during the critical postnatal phase of its development. While this work revealed timed execution of transcriptional programs, it was blind to the specific cells where gene expression changes occurred. Here, we utilized single-cell RNA-Seq to determine transcriptional profiles of each major MNTB cell type. We discerned directional signaling patterns between neuronal, glial, and vascular-associated cells for VEGF, TGFß, and Delta-Notch pathways during a robust period of vascular remodeling in the MNTB. Furthermore, we describe functional outcomes of the disruption of neuron-astrocyte fibroblast growth factor 9 (Fgf9) signaling. We used a conditional KO (cKO) approach to genetically delete Fgf9 from principal neurons in the MNTB, which led to an early onset of glial fibrillary acidic protein (Gfap) expression in astrocytes. In turn, Fgf9 cKO mice show increased levels of astrocyte-enriched brevican (Bcan), a component of the perineuronal net matrix that ensheaths principal neurons in the MNTB and the large calyx of Held terminal, while levels of the neuron-enriched hyaluronan and proteoglycan link protein 1 (Hapln1) were unchanged. Finally, volumetric analysis of vesicular glutamate transporters 1 and 2 (Vglut1/2), which serves as a proxy for terminal size, revealed an increase in calyx of Held volume in the Fgf9 cKO. Overall, we demonstrate a coordinated neuron-astrocyte Fgf9 signaling network that functions to regulate astrocyte maturation, perineuronal net structure, and synaptic refinement.

Rat Spinal Cord Injury Associated with Spasticity Leads to Widespread Changes in the Regulation of Retained Introns.

To determine molecular changes that correlate with long-term physiological changes after spinal cord injury associated with spasticity, we used a complete transection model with an injury at sacral spinal level S2, wherein tail spasms develop in rats weeks to months post-injury. Using Illumina and nanopore sequencing, we found that from 12,266 expressed genes roughly 11% (1,342) change expression levels in the rats with spasticity. The transcription factor PU.1 (Spi-1 proto-oncogene) and several of its known regulated genes were upregulated during injury, possibly reflecting changes in cellular composition. In contrast to widespread changes in gene expression, only a few changes in alternative exon usage could be detected because of injury. There were more than 1,000 changes in retained intron usage, however. Unexpectedly, most of these retained introns have not been described yet but could be validated using direct RNA nanopore sequencing. In addition to changes from injury, our model allowed regional analysis of gene expression. Comparing the segments rostral and caudal to the injury site in naïve animals showed 525 differentially regulated genes and differential regional use of retained introns. We did not detect changes in the serotonin receptor 2C editing that were implicated previously in this spinal cord injury model. Our data suggest that regulation of intron retention of polyadenylated pre-mRNA is an important regulatory mechanism in the spinal cord under both physiological and pathophysiological conditions.

Cwc27, associated with retinal degeneration, functions as a splicing factor in vivo.

Previous in vitro studies indicate that CWC27 functions as a splicing factor in the Bact spliceosome complex, interacting with CWC22 to form a landing platform for eIF4A3, a core component of the exon junction complex. However, the function of CWC27 as a splicing factor has not been validated in any in vivo systems. CWC27 variants have been shown to cause autosomal recessive retinal degeneration, in both syndromic and non-syndromic forms. The Cwc27K338fs/K338fs mouse model was shown to have significant retinal dysfunction and degeneration by 6 months of age. In this report, we have taken advantage of the Cwc27K338fs/K338fs mouse model to show that Cwc27 is involved in splicing in vivo in the context of the retina. Bulk RNA and single cell RNA-sequencing of the mouse retina showed that there were gene expression and splicing pattern changes, including alternative splice site usage and intron retention. Positive staining for CHOP suggests that ER stress may be activated in response to the splicing pattern changes and is a likely contributor to the disease mechanism. Our results provide the first evidence that CWC27 functions as a splicing factor in an in vivo context. The splicing defects and gene expression changes observed in the Cwc27K338fs/K338fs mouse retina provide insight to the potential disease mechanisms, paving the way for targeted therapeutic development.

Differential reprogramming of breast cancer subtypes in 3D cultures and implications for sensitivity to targeted therapy.

Screening for effective candidate drugs for breast cancer has shifted from two-dimensional (2D) to three-dimensional (3D) cultures. Here we systematically compared the transcriptomes of these different culture conditions by RNAseq of 14 BC cell lines cultured in both 2D and 3D conditions. All 3D BC cell cultures demonstrated increased mitochondrial metabolism and downregulated cell cycle programs. Luminal BC cells in 3D demonstrated overall limited reprogramming. 3D basal B BC cells showed increased expression of extracellular matrix (ECM) interaction genes, which coincides with an invasive phenotype not observed in other BC cells. Genes downregulated in 3D were associated with metastatic disease progression in BC patients, including cyclin dependent kinases and aurora kinases. Furthermore, the overall correlation of the cell line transcriptome to the BC patient transcriptome was increased in 3D cultures for all TNBC cell lines. To define the most optimal culture conditions to study the oncogenic pathway of interest, an open source bioinformatics strategy was established.

The Musashi proteins MSI1 and MSI2 are required for photoreceptor morphogenesis and vision in mice.

The Musashi family of RNA-binding proteins is known for its role in stem-cell renewal and is a negative regulator of cell differentiation. Interestingly, in the retina, the Musashi proteins MSI1 and MSI2 are differentially expressed throughout the cycle of retinal development, with MSI2 protein displaying robust expression in the adult retinal tissue. In this study, we investigated the importance of Musashi proteins in the development and function of photoreceptor neurons in the retina. We generated a pan-retinal and rod photoreceptor neuron-specific conditional KO mouse lacking MSI1 and MSI2. Independent of the sex, photoreceptor neurons with simultaneous deletion of Msi1 and Msi2 were unable to respond to light and displayed severely disrupted photoreceptor outer segment morphology and ciliary defects. Mice lacking MSI1 and MSI2 in the retina exhibited neuronal degeneration, with complete loss of photoreceptors within 6 months. In concordance with our earlier studies that proposed a role for Musashi proteins in regulating alternative splicing, the loss of MSI1 and MSI2 prevented the use of photoreceptor-specific exons in transcripts critical for outer segment morphogenesis, ciliogenesis, and synaptic transmission. Overall, we demonstrate a critical role for Musashi proteins in the morphogenesis of terminally differentiated photoreceptor neurons. This role is in stark contrast with the canonical function of these two proteins in the maintenance and renewal of stem cells.

Developmental Attenuation of Neuronal Apoptosis by Neural-Specific Splicing of Bak1 Microexon.

Continuous neuronal survival is vital for mammals because mammalian brains have limited regeneration capability. After neurogenesis, suppression of apoptosis is needed to ensure a neuron's long-term survival. Here we describe a robust genetic program that intrinsically attenuates apoptosis competence in neurons. Developmental downregulation of the splicing regulator PTBP1 in immature neurons allows neural-specific splicing of the evolutionarily conserved Bak1 microexon 5. Exon 5 inclusion triggers nonsense-mediated mRNA decay (NMD) and unproductive translation of Bak1 transcripts (N-Bak mRNA), leading to suppression of pro-apoptotic BAK1 proteins and allowing neurons to reduce apoptosis. Germline heterozygous ablation of exon 5 increases BAK1 proteins exclusively in the brain, inflates neuronal apoptosis, and leads to early postnatal mortality. Therefore, neural-specific exon 5 splicing and depletion of BAK1 proteins uniquely repress neuronal apoptosis. Although apoptosis is important for development, attenuation of apoptosis competence through neural-specific splicing of the Bak1 microexon is essential for neuronal and animal survival.

An activator of G protein-coupled receptor and MEK1/2-ERK1/2 signaling inhibits HIV-1 replication by altering viral RNA processing.

The ability of HIV-1 to evolve resistance to combined antiretroviral therapies (cARTs) has stimulated research into alternative means of controlling this infection. We assayed >60 modulators of RNA alternative splicing (AS) to identify new inhibitors of HIV-1 RNA processing-a segment of the viral lifecycle not targeted by current drugs-and discovered compound N-[4-chloro-3-(trifluoromethyl)phenyl]-7-nitro-2,1,3-benzoxadiazol-4-amine (5342191) as a potent inhibitor of both wild-type (Ba-L, NL4-3, LAI, IIIB, and N54) and drug-resistant strains of HIV-1 (IC50: ~700 nM) with no significant effect on cell viability at doses tested. 5342191 blocks expression of four essential HIV-1 structural and regulatory proteins (Gag, Env, Tat, and Rev) without affecting total protein synthesis of the cell. This response is associated with altered unspliced (US) and singly-spliced (SS) HIV-1 RNA accumulation (~60% reduction) and transport to the cytoplasm (loss of Rev) whereas parallel analysis of cellular RNAs revealed less than a 0.7% of host alternative splicing (AS) events (0.25-0.67% by ≥ 10-20%), gene expression (0.01-0.46% by ≥ 2-5 fold), and protein abundance (0.02-0.34% by ≥ 1.5-2 fold) being affected. Decreased expression of Tat, but not Gag/Env, upon 5342191 treatment was reversed by a proteasome inhibitor, suggesting that this compound alters the synthesis/degradation of this key viral factor. Consistent with an affect on HIV-1 RNA processing, 5342191 treatment of cells altered the abundance and phosphorylation of serine/arginine-rich splicing factor (SRSF) 1, 3, and 4. Despite the activation of several intracellular signaling pathways by 5342191 (Ras, MEK1/2-ERK1/2, and JNK1/2/3), inhibition of HIV-1 gene expression by this compound could be reversed by pre-treatment with either a G-protein α-subunit inhibitor or two different MEK1/2 inhibitors. These observations demonstrate enhanced sensitivity of HIV-1 gene expression to small changes in host RNA processing and highlights the potential of modulating host intracellular signaling as an alternative approach for controlling HIV-1 infection.

Uncovering the signaling landscape controlling breast cancer cell migration identifies novel metastasis driver genes.

Ttriple-negative breast cancer (TNBC) is an aggressive and highly metastatic breast cancer subtype. Enhanced TNBC cell motility is a prerequisite of TNBC cell dissemination. Here, we apply an imaging-based RNAi phenotypic cell migration screen using two highly motile TNBC cell lines (Hs578T and MDA-MB-231) to provide a repository of signaling determinants that functionally drive TNBC cell motility. We have screened ~4,200 target genes individually and discovered 133 and 113 migratory modulators of Hs578T and MDA-MB-231, respectively, which are linked to signaling networks predictive for breast cancer progression. The splicing factors PRPF4B and BUD31 and the transcription factor BPTF are essential for cancer cell migration, amplified in human primary breast tumors and associated with metastasis-free survival. Depletion of PRPF4B, BUD31 and BPTF causes primarily down regulation of genes involved in focal adhesion and ECM-interaction pathways. PRPF4B is essential for TNBC metastasis formation in vivo, making PRPF4B a candidate for further drug development.

ADP-Ribosylation Factor-Like 2 (ARL2) regulates cilia stability and development of outer segments in rod photoreceptor neurons.

Photoreceptor cells are specialized neurons with a sensory cilium carrying an elaborate membrane structure, the outer segment (OS). Inherited mutations in genes involved in ciliogenesis frequently result in OS malformation and blindness. ADP-ribosylation factor-like 2 (ARL2) has recently been implicated in OS formation through its association with Binder of ARL2 (BART or ARL2BP), a protein linked to inherited blinding disease. To test the role of ARL2 in vision we created a transgenic mouse model expressing a tagged-dominant active form of human ARL2 (ARL2-Q70L) under a rod-specific promoter. Transgenic ARL2-Q70L animals exhibit reduced photoreceptor cell function as early as post-natal day 16 and progressive rod degeneration. We attribute loss of photoreceptor function to the defective OS morphogenesis in the ARL2-Q70L transgenic model. ARL2-Q70L expression results in shortened inner and outer segments, shortened and mislocalized axonemes and cytoplasmic accumulation of rhodopsin. In conclusion, we show that ARL2-Q70L is crucial for photoreceptor neuron sensory cilium development. Future research will expand upon our hypothesis that ARL2-Q70L mutant interferes with microtubule maintenance and tubulin regulation resulting in impaired growth of the axoneme and elaboration of the photoreceptor outer segment.

Myeloid Disease Mutations of Splicing Factor SRSF2 Cause G2-M Arrest and Skewed Differentiation of Human Hematopoietic Stem and Progenitor Cells.

Myeloid malignancies, including myelodysplastic syndromes, chronic myelomonocytic leukemia, and acute myeloid leukemia, are characterized by abnormal proliferation and differentiation of hematopoietic stem and progenitor cells (HSPCs). Reports on analysis of bone marrow samples from patients have revealed a high incidence of mutations in splicing factors in early stem and progenitor cell clones, but the mechanisms underlying transformation of HSPCs harboring these mutations remain unknown. Using ex vivo cultures of primary human CD34+ cells as a model, we find that mutations in splicing factors SRSF2 and U2AF1 exert distinct effects on proliferation and differentiation of HSPCs. SRSF2 mutations cause a dramatic inhibition of proliferation via a G2-M phase arrest and induction of apoptosis. U2AF1 mutations, conversely, do not significantly affect proliferation. Mutations in both SRSF2 and U2AF1 cause abnormal differentiation by skewing granulo-monocytic differentiation toward monocytes but elicit diverse effects on megakaryo-erythroid differentiation. The SRSF2 mutations skew differentiation toward megakaryocytes whereas U2AF1 mutations cause an increase in the erythroid cell populations. These distinct functional consequences indicate that SRSF2 and U2AF1 mutations have cell context-specific effects and that the generation of myeloid disease phenotype by mutations in the genes coding these two proteins likely involves different intracellular mechanisms. Stem Cells 2018;36:1663-1675.

SEASTAR: systematic evaluation of alternative transcription start sites in RNA.

Alternative first exons diversify the transcriptomes of eukaryotes by producing variants of the 5' Untranslated Regions (5'UTRs) and N-terminal coding sequences. Accurate transcriptome-wide detection of alternative first exons typically requires specialized experimental approaches that are designed to identify the 5' ends of transcripts. We developed a computational pipeline SEASTAR that identifies first exons from RNA-seq data alone then quantifies and compares alternative first exon usage across multiple biological conditions. The exons inferred by SEASTAR coincide with transcription start sites identified directly by CAGE experiments and bear epigenetic hallmarks of active promoters. To determine if differential usage of alternative first exons can yield insights into the mechanism controlling gene expression, we applied SEASTAR to an RNA-seq dataset that tracked the reprogramming of mouse fibroblasts into induced pluripotent stem cells. We observed dynamic temporal changes in the usage of alternative first exons, along with correlated changes in transcription factor expression. Using a combined sequence motif and gene set enrichment analysis we identified N-Myc as a regulator of alternative first exon usage in the pluripotent state. Our results demonstrate that SEASTAR can leverage the available RNA-seq data to gain insights into the control of gene expression and alternative transcript variation in eukaryotic transcriptomes.

Bardet-Biedl syndrome-8 (BBS8) protein is crucial for the development of outer segments in photoreceptor neurons.

Bardet-Biedl syndrome (BBS) is an autosomal recessive ciliopathy characterized by developmental abnormalities and vision loss. To date, mutations in 21 genes have been linked to BBS. The products of eight of these BBS genes form a stable octameric complex termed the BBSome. Mutations in BBS8, a component of the BBSome, cause early vision loss, but the role of BBS8 in supporting vision is not known. To understand the mechanisms by which BBS8 supports rod and cone photoreceptor function, we generated animal models lacking BBS8. The loss of BBS8 protein led to concomitant decrease in the levels of BBSome subunits, BBS2 and BBS5 and increase in the levels of the BBS1 and BBS4 subunits. BBS8 ablation was associated with severe reduction of rod and cone photoreceptor function and progressive degeneration of each photoreceptor subtype. We observed disorganized and shortened photoreceptor outer segments (OS) at post-natal day 10 as the OS elaborates. Interestingly, loss of BBS8 led to changes in the distribution of photoreceptor axonemal proteins and hyper-acetylation of ciliary microtubules. In contrast to properly localized phototransduction machinery, we observed OS accumulation of syntaxin3, a protein normally found in the cytoplasm and the synaptic termini. In conclusion, our studies demonstrate the requirement for BBS8 in early development and elaboration of ciliated photoreceptor OS, explaining the need for BBS8 in normal vision. The findings from our study also imply that early targeting of both rods and cones in BBS8 patients is crucial for successful restoration of vision.

Mutations in the Spliceosome Component CWC27 Cause Retinal Degeneration with or without Additional Developmental Anomalies.

Pre-mRNA splicing factors play a fundamental role in regulating transcript diversity both temporally and spatially. Genetic defects in several spliceosome components have been linked to a set of non-overlapping spliceosomopathy phenotypes in humans, among which skeletal developmental defects and non-syndromic retinitis pigmentosa (RP) are frequent findings. Here we report that defects in spliceosome-associated protein CWC27 are associated with a spectrum of disease phenotypes ranging from isolated RP to severe syndromic forms. By whole-exome sequencing, recessive protein-truncating mutations in CWC27 were found in seven unrelated families that show a range of clinical phenotypes, including retinal degeneration, brachydactyly, craniofacial abnormalities, short stature, and neurological defects. Remarkably, variable expressivity of the human phenotype can be recapitulated in Cwc27 mutant mouse models, with significant embryonic lethality and severe phenotypes in the complete knockout mice while mice with a partial loss-of-function allele mimic the isolated retinal degeneration phenotype. Our study describes a retinal dystrophy-related phenotype spectrum as well as its genetic etiology and highlights the complexity of the spliceosomal gene network.