Systematic analysis of alternative exon-dependent interactome remodeling reveals multitasking functions of gene regulatory factors.

Alternative splicing significantly expands biological complexity, particularly in the vertebrate nervous system. Increasing evidence indicates that developmental and tissue-dependent alternative exons often control protein-protein interactions; yet, only a minor fraction of these events have been characterized. Using affinity purification-mass spectrometry (AP-MS), we show that approximately 60% of analyzed neural-differential exons in proteins previously implicated in transcriptional regulation result in the gain or loss of interaction partners, which in some cases form unexpected links with coupled processes. Notably, a neural exon in Chtop regulates its interaction with the Prmt1 methyltransferase and DExD-Box helicases Ddx39b/a, affecting its methylation and activity in promoting RNA export. Additionally, a neural exon in Sap30bp affects interactions with RNA processing factors, modulating a critical function of Sap30bp in promoting the splicing of <100 nt "mini-introns" that control nuclear RNA levels. AP-MS is thus a powerful approach for elucidating the multifaceted functions of proteins imparted by context-dependent alternative exons.

Survival-based CRISPR genetic screens across a panel of permissive cell lines identify common and cell-specific SARS-CoV-2 host factors.

SARS-CoV-2 depends on host cell components for infection and replication. Identification of virus-host dependencies offers an effective way to elucidate mechanisms involved in viral infection and replication. If druggable, host factor dependencies may present an attractive strategy for anti-viral therapy. In this study, we performed genome wide CRISPR knockout screens in Vero E6 cells and four human cell lines including Calu-3, UM-UC-4, HEK-293 and HuH-7 to identify genetic regulators of SARS-CoV-2 infection. Our findings identified only ACE2, the cognate SARS-CoV-2 entry receptor, as a common host dependency factor across all cell lines, while other host genes identified were largely cell line specific, including known factors TMPRSS2 and CTSL. Several of the discovered host-dependency factors converged on pathways involved in cell signalling, immune-related pathways, and chromatin modification. Notably, the chromatin modifier gene KMT2C in Calu-3 cells had the strongest impact in preventing SARS-CoV-2 infection when perturbed.

Regulation of alternative polyadenylation by the C2H2-zinc-finger protein Sp1.

Alternative polyadenylation (APA) enhances gene regulatory potential by increasing the diversity of mRNA transcripts. 3' UTR shortening through APA correlates with enhanced cellular proliferation and is a widespread phenomenon in tumor cells. Here, we show that the ubiquitously expressed transcription factor Sp1 binds RNA in vivo and is a common repressor of distal poly(A) site usage. RNA sequencing identified 2,344 genes (36% of the total mapped mRNA transcripts) with lengthened 3' UTRs upon Sp1 depletion. Sp1 preferentially binds the 3' UTRs of such lengthened transcripts and inhibits cleavage at distal sites by interacting with the subunits of the core cleavage and polyadenylation (CPA) machinery. The 3' UTR lengths of Sp1 target genes in breast cancer patient RNA-seq data correlate with Sp1 expression levels, implicating Sp1-mediated APA regulation in modulating tumorigenic properties. Taken together, our findings provide insights into the mechanism for dynamic APA regulation by unraveling a previously unknown function of the DNA-binding transcription factor Sp1.

Systematic exploration of dynamic splicing networks reveals conserved multistage regulators of neurogenesis.

Alternative splicing (AS) is a critical regulatory layer; yet, factors controlling functionally coordinated splicing programs during developmental transitions are poorly understood. Here, we employ a screening strategy to identify factors controlling dynamic splicing events important for mammalian neurogenesis. Among previously unknown regulators, Rbm38 acts widely to negatively control neural AS, in part through interactions mediated by the established repressor of splicing, Ptbp1. Puf60, a ubiquitous factor, is surprisingly found to promote neural splicing patterns. This activity requires a conserved, neural-differential exon that remodels Puf60 co-factor interactions. Ablation of this exon rewires distinct AS networks in embryonic stem cells and at different stages of mouse neurogenesis. Single-cell transcriptome analyses further reveal distinct roles for Rbm38 and Puf60 isoforms in establishing neuronal identity. Our results describe important roles for previously unknown regulators of neurogenesis and establish how an alternative exon in a widely expressed splicing factor orchestrates temporal control over cell differentiation.

Definition of germ layer cell lineage alternative splicing programs reveals a critical role for Quaking in specifying cardiac cell fate.

Alternative splicing is critical for development; however, its role in the specification of the three embryonic germ layers is poorly understood. By performing RNA-Seq on human embryonic stem cells (hESCs) and derived definitive endoderm, cardiac mesoderm, and ectoderm cell lineages, we detect distinct alternative splicing programs associated with each lineage. The most prominent splicing program differences are observed between definitive endoderm and cardiac mesoderm. Integrative multi-omics analyses link each program with lineage-enriched RNA binding protein regulators, and further suggest a widespread role for Quaking (QKI) in the specification of cardiac mesoderm. Remarkably, knockout of QKI disrupts the cardiac mesoderm-associated alternative splicing program and formation of myocytes. These changes arise in part through reduced expression of BIN1 splice variants linked to cardiac development. Mechanistically, we find that QKI represses inclusion of exon 7 in BIN1 pre-mRNA via an exonic ACUAA motif, and this is concomitant with intron removal and cleavage from chromatin. Collectively, our results uncover alternative splicing programs associated with the three germ lineages and demonstrate an important role for QKI in the formation of cardiac mesoderm.

Systematic mapping of nuclear domain-associated transcripts reveals speckles and lamina as hubs of functionally distinct retained introns.

The nucleus is highly compartmentalized through the formation of distinct classes of membraneless domains. However, the composition and function of many of these structures are not well understood. Using APEX2-mediated proximity labeling and RNA sequencing, we surveyed human transcripts associated with nuclear speckles, several additional domains, and the lamina. Remarkably, speckles and lamina are associated with distinct classes of retained introns enriched in genes that function in RNA processing, translation, and the cell cycle, among other processes. In contrast to the lamina-proximal introns, retained introns associated with speckles are relatively short, GC-rich, and enriched for functional sites of RNA-binding proteins that are concentrated in these domains. They are also highly differentially regulated across diverse cellular contexts, including the cell cycle. Thus, our study provides a resource of nuclear domain-associated transcripts and further reveals speckles and lamina as hubs of distinct populations of retained introns linked to gene regulation and cell cycle progression.

SARS-CoV-2 nucleocapsid protein binds host mRNAs and attenuates stress granules to impair host stress response.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nucleocapsid (N) protein is essential for viral replication, making it a promising target for antiviral drug and vaccine development. SARS-CoV-2 infected patients exhibit an uncoordinated immune response; however, the underlying mechanistic details of this imbalance remain obscure. Here, starting from a functional proteomics workflow, we cataloged the protein-protein interactions of SARS-CoV-2 proteins, including an evolutionarily conserved specific interaction of N with the stress granule resident proteins G3BP1 and G3BP2. N localizes to stress granules and sequesters G3BPs away from their typical interaction partners, thus attenuating stress granule formation. We found that N binds directly to host mRNAs in cells, with a preference for 3' UTRs, and modulates target mRNA stability. We show that the N protein rewires the G3BP1 mRNA-binding profile and suppresses the physiological stress response of host cells, which may explain the imbalanced immune response observed in SARS-CoV-2 infected patients.

Nuclear compartmentalization of TERT mRNA and TUG1 lncRNA is driven by intron retention.

The spatial partitioning of the transcriptome in the cell is an important form of gene-expression regulation. Here, we address how intron retention influences the spatio-temporal dynamics of transcripts from two clinically relevant genes: TERT (Telomerase Reverse Transcriptase) pre-mRNA and TUG1 (Taurine-Upregulated Gene 1) lncRNA. Single molecule RNA FISH reveals that nuclear TERT transcripts uniformly and robustly retain specific introns. Our data suggest that the splicing of TERT retained introns occurs during mitosis. In contrast, TUG1 has a bimodal distribution of fully spliced cytoplasmic and intron-retained nuclear transcripts. We further test the functionality of intron-retention events using RNA-targeting thiomorpholino antisense oligonucleotides to block intron excision. We show that intron retention is the driving force for the nuclear compartmentalization of these RNAs. For both RNAs, altering this splicing-driven subcellular distribution has significant effects on cell viability. Together, these findings show that stable retention of specific introns can orchestrate spatial compartmentalization of these RNAs within the cell. This process reveals that modulating RNA localization via targeted intron retention can be utilized for RNA-based therapies.

A multiplexed, next generation sequencing platform for high-throughput detection of SARS-CoV-2.

Population scale sweeps of viral pathogens, such as SARS-CoV-2, require high intensity testing for effective management. Here, we describe "Systematic Parallel Analysis of RNA coupled to Sequencing for Covid-19 screening" (C19-SPAR-Seq), a multiplexed, scalable, readily automated platform for SARS-CoV-2 detection that is capable of analyzing tens of thousands of patient samples in a single run. To address strict requirements for control of assay parameters and output demanded by clinical diagnostics, we employ a control-based Precision-Recall and Receiver Operator Characteristics (coPR) analysis to assign run-specific quality control metrics. C19-SPAR-Seq coupled to coPR on a trial cohort of several hundred patients performs with a specificity of 100% and sensitivity of 91% on samples with low viral loads, and a sensitivity of >95% on high viral loads associated with disease onset and peak transmissibility. This study establishes the feasibility of employing C19-SPAR-Seq for the large-scale monitoring of SARS-CoV-2 and other pathogens.

Actionable Cytopathogenic Host Responses of Human Alveolar Type 2 Cells to SARS-CoV-2.

Human transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causative pathogen of the COVID-19 pandemic, exerts a massive health and socioeconomic crisis. The virus infects alveolar epithelial type 2 cells (AT2s), leading to lung injury and impaired gas exchange, but the mechanisms driving infection and pathology are unclear. We performed a quantitative phosphoproteomic survey of induced pluripotent stem cell-derived AT2s (iAT2s) infected with SARS-CoV-2 at air-liquid interface (ALI). Time course analysis revealed rapid remodeling of diverse host systems, including signaling, RNA processing, translation, metabolism, nuclear integrity, protein trafficking, and cytoskeletal-microtubule organization, leading to cell cycle arrest, genotoxic stress, and innate immunity. Comparison to analogous data from transformed cell lines revealed respiratory-specific processes hijacked by SARS-CoV-2, highlighting potential novel therapeutic avenues that were validated by a high hit rate in a targeted small molecule screen in our iAT2 ALI system.

Genetic interaction mapping and exon-resolution functional genomics with a hybrid Cas9-Cas12a platform.

Systematic mapping of genetic interactions (GIs) and interrogation of the functions of sizable genomic segments in mammalian cells represent important goals of biomedical research. To advance these goals, we present a CRISPR (clustered regularly interspaced short palindromic repeats)-based screening system for combinatorial genetic manipulation that employs coexpression of CRISPR-associated nucleases 9 and 12a (Cas9 and Cas12a) and machine-learning-optimized libraries of hybrid Cas9-Cas12a guide RNAs. This system, named Cas Hybrid for Multiplexed Editing and screening Applications (CHyMErA), outperforms genetic screens using Cas9 or Cas12a editing alone. Application of CHyMErA to the ablation of mammalian paralog gene pairs reveals extensive GIs and uncovers phenotypes normally masked by functional redundancy. Application of CHyMErA in a chemogenetic interaction screen identifies genes that impact cell growth in response to mTOR pathway inhibition. Moreover, by systematically targeting thousands of alternative splicing events, CHyMErA identifies exons underlying human cell line fitness. CHyMErA thus represents an effective screening approach for GI mapping and the functional analysis of sizable genomic regions, such as alternative exons.

Autism-Misregulated eIF4G Microexons Control Synaptic Translation and Higher Order Cognitive Functions.

Microexons represent the most highly conserved class of alternative splicing, yet their functions are poorly understood. Here, we focus on closely related neuronal microexons overlapping prion-like domains in the translation initiation factors, eIF4G1 and eIF4G3, the splicing of which is activity dependent and frequently disrupted in autism. CRISPR-Cas9 deletion of these microexons selectively upregulates synaptic proteins that control neuronal activity and plasticity and further triggers a gene expression program mirroring that of activated neurons. Mice lacking the Eif4g1 microexon display social behavior, learning, and memory deficits, accompanied by altered hippocampal synaptic plasticity. We provide evidence that the eIF4G microexons function as a translational brake by causing ribosome stalling, through their propensity to promote the coalescence of cytoplasmic granule components associated with translation repression, including the fragile X mental retardation protein FMRP. The results thus reveal an autism-disrupted mechanism by which alternative splicing specializes neuronal translation to control higher order cognitive functioning.

A slow transcription rate causes embryonic lethality and perturbs kinetic coupling of neuronal genes.

The rate of RNA polymerase II (RNAPII) elongation has an important role in the control of alternative splicing (AS); however, the in vivo consequences of an altered elongation rate are unknown. Here, we generated mouse embryonic stem cells (ESCs) knocked in for a slow elongating form of RNAPII We show that a reduced transcriptional elongation rate results in early embryonic lethality in mice. Focusing on neuronal differentiation as a model, we observed that slow elongation impairs development of the neural lineage from ESCs, which is accompanied by changes in AS and in gene expression along this pathway. In particular, we found a crucial role for RNAPII elongation rate in transcription and splicing of long neuronal genes involved in synapse signaling. The impact of the kinetic coupling of RNAPII elongation rate with AS is greater in ESC-differentiated neurons than in pluripotent cells. Our results demonstrate the requirement for an appropriate transcriptional elongation rate to ensure proper gene expression and to regulate AS during development.

Genome-wide CRISPR-Cas9 Interrogation of Splicing Networks Reveals a Mechanism for Recognition of Autism-Misregulated Neuronal Microexons.

Alternative splicing is crucial for diverse cellular, developmental, and pathological processes. However, the full networks of factors that control individual splicing events are not known. Here, we describe a CRISPR-based strategy for the genome-wide elucidation of pathways that control splicing and apply it to microexons with important functions in nervous system development and that are commonly misregulated in autism. Approximately 200 genes associated with functionally diverse regulatory layers and enriched in genetic links to autism control neuronal microexons. Remarkably, the widely expressed RNA binding proteins Srsf11 and Rnps1 directly, preferentially, and frequently co-activate these microexons. These factors form critical interactions with the neuronal splicing regulator Srrm4 and a bi-partite intronic splicing enhancer element to promote spliceosome formation. Our study thus presents a versatile system for the identification of entire splicing regulatory pathways and further reveals a common mechanism for the definition of neuronal microexons that is disrupted in autism.

An atlas of alternative splicing profiles and functional associations reveals new regulatory programs and genes that simultaneously express multiple major isoforms.

Alternative splicing (AS) generates remarkable regulatory and proteomic complexity in metazoans. However, the functions of most AS events are not known, and programs of regulated splicing remain to be identified. To address these challenges, we describe the Vertebrate Alternative Splicing and Transcription Database (VastDB), the largest resource of genome-wide, quantitative profiles of AS events assembled to date. VastDB provides readily accessible quantitative information on the inclusion levels and functional associations of AS events detected in RNA-seq data from diverse vertebrate cell and tissue types, as well as developmental stages. The VastDB profiles reveal extensive new intergenic and intragenic regulatory relationships among different classes of AS and previously unknown and conserved landscapes of tissue-regulated exons. Contrary to recent reports concluding that nearly all human genes express a single major isoform, VastDB provides evidence that at least 48% of multiexonic protein-coding genes express multiple splice variants that are highly regulated in a cell/tissue-specific manner, and that >18% of genes simultaneously express multiple major isoforms across diverse cell and tissue types. Isoforms encoded by the latter set of genes are generally coexpressed in the same cells and are often engaged by translating ribosomes. Moreover, they are encoded by genes that are significantly enriched in functions associated with transcriptional control, implying they may have an important and wide-ranging role in controlling cellular activities. VastDB thus provides an unprecedented resource for investigations of AS function and regulation.

Multilayered Control of Alternative Splicing Regulatory Networks by Transcription Factors.

Networks of coordinated alternative splicing (AS) events play critical roles in development and disease. However, a comprehensive knowledge of the factors that regulate these networks is lacking. We describe a high-throughput system for systematically linking trans-acting factors to endogenous RNA regulatory events. Using this system, we identify hundreds of factors associated with diverse regulatory layers that positively or negatively control AS events linked to cell fate. Remarkably, more than one-third of the regulators are transcription factors. Further analyses of the zinc finger protein Zfp871 and BTB/POZ domain transcription factor Nacc1, which regulate neural and stem cell AS programs, respectively, reveal roles in controlling the expression of specific splicing regulators. Surprisingly, these proteins also appear to regulate target AS programs via binding RNA. Our results thus uncover a large "missing cache" of splicing regulators among annotated transcription factors, some of which dually regulate AS through direct and indirect mechanisms.

Control of embryonic stem cell self-renewal and differentiation via coordinated alternative splicing and translation of YY2.

Translational control of gene expression plays a key role during the early phases of embryonic development. Here we describe a transcriptional regulator of mouse embryonic stem cells (mESCs), Yin-yang 2 (YY2), that is controlled by the translation inhibitors, Eukaryotic initiation factor 4E-binding proteins (4E-BPs). YY2 plays a critical role in regulating mESC functions through control of key pluripotency factors, including Octamer-binding protein 4 (Oct4) and Estrogen-related receptor-ß (Esrrb). Importantly, overexpression of YY2 directs the differentiation of mESCs into cardiovascular lineages. We show that the splicing regulator Polypyrimidine tract-binding protein 1 (PTBP1) promotes the retention of an intron in the 5'-UTR of Yy2 mRNA that confers sensitivity to 4E-BP-mediated translational suppression. Thus, we conclude that YY2 is a major regulator of mESC self-renewal and lineage commitment and document a multilayer regulatory mechanism that controls its expression.

SMN and symmetric arginine dimethylation of RNA polymerase II C-terminal domain control termination.

The carboxy-terminal domain (CTD) of the RNA polymerase II (RNAP II) subunit POLR2A is a platform for modifications specifying the recruitment of factors that regulate transcription, mRNA processing, and chromatin remodelling. Here we show that a CTD arginine residue (R1810 in human) that is conserved across vertebrates is symmetrically dimethylated (me2s). This R1810me2s modification requires protein arginine methyltransferase 5 (PRMT5) and recruits the Tudor domain of the survival of motor neuron (SMN, also known as GEMIN1) protein, which is mutated in spinal muscular atrophy. SMN interacts with senataxin, which is sometimes mutated in ataxia oculomotor apraxia type 2 and amyotrophic lateral sclerosis. Because POLR2A R1810me2s and SMN, like senataxin, are required for resolving RNA-DNA hybrids created by RNA polymerase II that form R-loops in transcription termination regions, we propose that R1810me2s, SMN, and senataxin are components of an R-loop resolution pathway. Defects in this pathway can influence transcription termination and may contribute to neurodegenerative disorders.