A UL26-PIAS1 complex antagonizes anti-viral gene expression during Human Cytomegalovirus infection.

Viral disruption of innate immune signaling is a critical determinant of productive infection. The Human Cytomegalovirus (HCMV) UL26 protein prevents anti-viral gene expression during infection, yet the mechanisms involved are unclear. We used TurboID-driven proximity proteomics to identify putative UL26 interacting proteins during infection to address this issue. We find that UL26 forms a complex with several immuno-regulatory proteins, including several STAT family members and various PIAS proteins, a family of E3 SUMO ligases. Our results indicate that UL26 prevents STAT phosphorylation during infection and antagonizes transcriptional activation induced by either interferon α (IFNA) or tumor necrosis factor α (TNFα). Additionally, we find that the inactivation of PIAS1 sensitizes cells to inflammatory stimulation, resulting in an anti-viral transcriptional environment similar to ΔUL26 infection. Further, PIAS1 is important for HCMV cell-to-cell spread, which depends on the presence of UL26, suggesting that the UL26-PIAS1 interaction is vital for modulating intrinsic anti-viral defense.

Regulation of the Drosophila transcriptome by Pumilio and the CCR4-NOT deadenylase complex.

The sequence-specific RNA-binding protein Pumilio (Pum) controls Drosophila development; however, the network of mRNAs that it regulates remains incompletely characterized. In this study, we use knockdown and knockout approaches coupled with RNA-seq to measure the impact of Pum on the transcriptome of Drosophila cells in culture. We also use an improved RNA coimmunoprecipitation method to identify Pum-bound mRNAs in Drosophila embryos. Integration of these data sets with the locations of Pum-binding motifs across the transcriptome reveals novel direct Pum target genes involved in neural, muscle, wing, and germ cell development and in cellular proliferation. These genes include components of Wnt, TGF-ß, MAPK/ERK, and Notch signaling pathways, DNA replication, and lipid metabolism. We identify the mRNAs regulated by the CCR4-NOT deadenylase complex, a key factor in Pum-mediated repression, and observe concordant regulation of Pum:CCR4-NOT target mRNAs. Computational modeling reveals that Pum binding, binding site number, clustering, and sequence context are important determinants of regulation. In contrast, we show that the responses of direct mRNA targets to Pum-mediated repression are not influenced by the content of optimal synonymous codons. Moreover, contrary to a prevailing model, we do not detect a role for CCR4-NOT in the degradation of mRNAs with low codon optimality. Together, the results of this work provide new insights into the Pum regulatory network and mechanisms and the parameters that influence the efficacy of Pum-mediated regulation.

PolyAMiner-Bulk is a deep learning-based algorithm that decodes alternative polyadenylation dynamics from bulk RNA-seq data.

Alternative polyadenylation (APA) is a key post-transcriptional regulatory mechanism; yet, its regulation and impact on human diseases remain understudied. Existing bulk RNA sequencing (RNA-seq)-based APA methods predominantly rely on predefined annotations, severely impacting their ability to decode novel tissue- and disease-specific APA changes. Furthermore, they only account for the most proximal and distal cleavage and polyadenylation sites (C/PASs). Deconvoluting overlapping C/PASs and the inherent noisy 3' UTR coverage in bulk RNA-seq data pose additional challenges. To overcome these limitations, we introduce PolyAMiner-Bulk, an attention-based deep learning algorithm that accurately recapitulates C/PAS sequence grammar, resolves overlapping C/PASs, captures non-proximal-to-distal APA changes, and generates visualizations to illustrate APA dynamics. Evaluation on multiple datasets strongly evinces the performance merit of PolyAMiner-Bulk, accurately identifying more APA changes compared with other methods. With the growing importance of APA and the abundance of bulk RNA-seq data, PolyAMiner-Bulk establishes a robust paradigm of APA analysis.

Coevolution combined with molecular dynamics simulations provides structural and mechanistic insights into the interactions between the integrator complex subunits.

Finding the 3D structure of large, multi-subunit complexes is difficult, despite recent advances in cryo-EM technology, due to remaining challenges to expressing and purifying subunits. Computational approaches that predict protein-protein interactions, including Direct Coupling Analysis (DCA), represent an attractive alternative for dissecting interactions within protein complexes. However, they are readily applicable only to small proteins due to high computational complexity and a high number of false positives. To solve this problem, we proposed a modified DCA approach, a powerful tool to predict the most likely interfaces of protein complexes. Since our modified approach cannot provide structural and mechanistic details of interacting peptides, we combine it with Molecular Dynamics (MD) simulations. To illustrate this novel approach, we predict interacting domains and structural details of interactions of two Integrator complex subunits, INTS9 and INTS11. Our predictions of interacting residues of INTS9/INTS11 are highly consistent with crystallographic structure. We then expand our procedure to two complexes whose structures are not well-studied: 1) The heterodimer formed by the Cleavage and Polyadenylation Specificity Factor 100-kD (CPSF100) and 73-kD (CPSF73); 2) The heterotrimer formed by INTS4/INTS9/INTS11. Experimental data supports our predictions of interactions within these two complexes, demonstrating that combining DCA and MD simulations is a powerful approach to revealing structural insights of large protein complexes.

Regulation of the Drosophila transcriptome by Pumilio and CCR4-NOT deadenylase.

The sequence-specific RNA-binding protein Pumilio controls development of Drosophila; however, the network of mRNAs that it regulates remains incompletely characterized. In this study, we utilize knockdown and knockout approaches coupled with RNA-Seq to measure the impact of Pumilio on the transcriptome of Drosophila cells. We also used an improved RNA co-immunoprecipitation method to identify Pumilio bound mRNAs in Drosophila embryos. Integration of these datasets with the content of Pumilio binding motifs across the transcriptome revealed novel direct Pumilio target genes involved in neural, muscle, wing, and germ cell development, and cellular proliferation. These genes include components of Wnt, TGF-beta, MAPK/ERK, and Notch signaling pathways, DNA replication, and lipid metabolism. Additionally, we identified the mRNAs regulated by the CCR4-NOT deadenylase complex, a key factor in Pumilio-mediated repression, and observed concordant regulation of Pumilio:CCR4-NOT target mRNAs. Computational modeling revealed that Pumilio binding, binding site number, density, and sequence context are important determinants of regulation. Moreover, the content of optimal synonymous codons in target mRNAs exhibits a striking functional relationship to Pumilio and CCR4-NOT regulation, indicating that the inherent translation efficiency and stability of the mRNA modulates their response to these trans-acting regulatory factors. Together, the results of this work provide new insights into the Pumilio regulatory network and mechanisms, and the parameters that influence the efficacy of Pumilio-mediated regulation.

NUDT21 alters glioma migration through differential alternative polyadenylation of LAMC1.

PURPOSE: Gliomas and their surrounding microenvironment constantly interact to promote tumorigenicity, yet the underlying posttranscriptional regulatory mechanisms that govern this interplay are poorly understood. METHODS: Utilizing our established PAC-seq approach and PolyAMiner bioinformatic analysis pipeline, we deciphered the NUDT21-mediated differential APA dynamics in glioma cells. RESULTS: We identified LAMC1 as a critical NUDT21 alternative polyadenylation (APA) target, common in several core glioma-driving signaling pathways. qRT-PCR analysis confirmed that NUDT21-knockdown in glioma cells results in the preferred usage of the proximal polyA signal (PAS) of LAMC1. Functional studies revealed that NUDT21-knockdown-induced 3'UTR shortening of LAMC1 is sufficient to cause translational gain, as LAMC1 protein is upregulated in these cells compared to their respective controls. We demonstrate that 3'UTR shortening of LAMC1 after NUDT21 knockdown removes binding sites for miR-124/506, thereby relieving potent miRNA-based repression of LAMC1 expression. Remarkably, we report that the knockdown of NUDT21 significantly promoted glioma cell migration and that co-depletion of LAMC1 with NUDT21 abolished this effect. Lastly, we observed that LAMC1 3'UTR shortening predicts poor prognosis of low-grade glioma patients from The Cancer Genome Atlas. CONCLUSION: This study identifies NUDT21 as a core alternative polyadenylation factor that regulates the tumor microenvironment through differential APA and loss of miR-124/506 inhibition of LAMC1. Knockdown of NUDT21 in GBM cells mediates 3'UTR shortening of LAMC1, contributing to an increase in LAMC1, increased glioma cell migration/invasion, and a poor prognosis.

Bi-allelic variants in INTS11 are associated with a complex neurological disorder.

The Integrator complex is a multi-subunit protein complex that regulates the processing of nascent RNAs transcribed by RNA polymerase II (RNAPII), including small nuclear RNAs, enhancer RNAs, telomeric RNAs, viral RNAs, and protein-coding mRNAs. Integrator subunit 11 (INTS11) is the catalytic subunit that cleaves nascent RNAs, but, to date, mutations in this subunit have not been linked to human disease. Here, we describe 15 individuals from 10 unrelated families with bi-allelic variants in INTS11 who present with global developmental and language delay, intellectual disability, impaired motor development, and brain atrophy. Consistent with human observations, we find that the fly ortholog of INTS11, dIntS11, is essential and expressed in the central nervous systems in a subset of neurons and most glia in larval and adult stages. Using Drosophila as a model, we investigated the effect of seven variants. We found that two (p.Arg17Leu and p.His414Tyr) fail to rescue the lethality of null mutants, indicating that they are strong loss-of-function variants. Furthermore, we found that five variants (p.Gly55Ser, p.Leu138Phe, p.Lys396Glu, p.Val517Met, and p.Ile553Glu) rescue lethality but cause a shortened lifespan and bang sensitivity and affect locomotor activity, indicating that they are partial loss-of-function variants. Altogether, our results provide compelling evidence that integrity of the Integrator RNA endonuclease is critical for brain development.

PolyAMiner-Bulk: A Machine Learning Based Bioinformatics Algorithm to Infer and Decode Alternative Polyadenylation Dynamics from bulk RNA-seq data.

More than half of human genes exercise alternative polyadenylation (APA) and generate mRNA transcripts with varying 3' untranslated regions (UTR). However, current computational approaches for identifying cleavage and polyadenylation sites (C/PASs) and quantifying 3'UTR length changes from bulk RNA-seq data fail to unravel tissue- and disease-specific APA dynamics. Here, we developed a next-generation bioinformatics algorithm and application, PolyAMiner-Bulk, that utilizes an attention-based machine learning architecture and an improved vector projection-based engine to infer differential APA dynamics accurately. When applied to earlier studies, PolyAMiner-Bulk accurately identified more than twice the number of APA changes in an RBM17 knockdown bulk RNA-seq dataset compared to current generation tools. Moreover, on a separate dataset, PolyAMiner-Bulk revealed novel APA dynamics and pathways in scleroderma pathology and identified differential APA in a gene that was identified as being involved in scleroderma pathogenesis in an independent study. Lastly, we used PolyAMiner-Bulk to analyze the RNA-seq data of post-mortem prefrontal cortexes from the ROSMAP data consortium and unraveled novel APA dynamics in Alzheimer's Disease. Our method, PolyAMiner-Bulk, creates a paradigm for future alternative polyadenylation analysis from bulk RNA-seq data.

Alternative polyadenylation transcriptome-wide association study identifies APA-linked susceptibility genes in brain disorders.

Alternative polyadenylation (APA) plays an essential role in brain development; however, current transcriptome-wide association studies (TWAS) largely overlook APA in nominating susceptibility genes. Here, we performed a 3' untranslated region (3'UTR) APA TWAS (3'aTWAS) for 11 brain disorders by combining their genome-wide association studies data with 17,300 RNA-seq samples across 2,937 individuals. We identified 354 3'aTWAS-significant genes, including known APA-linked risk genes, such as SNCA in Parkinson's disease. Among these 354 genes, ~57% are not significant in traditional expression- and splicing-TWAS studies, since APA may regulate the translation, localization and protein-protein interaction of the target genes independent of mRNA level expression or splicing. Furthermore, we discovered ATXN3 as a 3'aTWAS-significant gene for amyotrophic lateral sclerosis, and its modulation substantially impacted pathological hallmarks of amyotrophic lateral sclerosis in vitro. Together, 3'aTWAS is a powerful strategy to nominate important APA-linked brain disorder susceptibility genes, most of which are largely overlooked by conventional expression and splicing analyses.

Alternative polyadenylation alters protein dosage by switching between intronic and 3'UTR sites.

Alternative polyadenylation (APA) creates distinct transcripts from the same gene by cleaving the pre-mRNA at poly(A) sites that can lie within the 3' untranslated region (3'UTR), introns, or exons. Most studies focus on APA within the 3'UTR; however, here, we show that CPSF6 insufficiency alters protein levels and causes a developmental syndrome by deregulating APA throughout the transcript. In neonatal humans and zebrafish larvae, CPSF6 insufficiency shifts poly(A) site usage between the 3'UTR and internal sites in a pathway-specific manner. Genes associated with neuronal function undergo mostly intronic APA, reducing their expression, while genes associated with heart and skeletal function mostly undergo 3'UTR APA and are up-regulated. This suggests that, under healthy conditions, cells toggle between internal and 3'UTR APA to modulate protein expression.

Integrator is a global promoter-proximal termination complex.

Integrator is a metazoan-specific protein complex capable of inducing termination at all RNAPII-transcribed loci. Integrator recognizes paused, promoter-proximal RNAPII and drives premature termination using dual enzymatic activities: an endonuclease that cleaves nascent RNA and a protein phosphatase that removes stimulatory phosphorylation associated with RNAPII pause release and productive elongation. Recent breakthroughs in structural biology have revealed the overall architecture of Integrator and provided insights into how multiple Integrator modules are coordinated to elicit termination effectively. Furthermore, functional genomics and biochemical studies have unraveled how Integrator-mediated termination impacts protein-coding and noncoding loci. Here, we review the current knowledge about the assembly and activity of Integrator and describe the role of Integrator in gene regulation, highlighting the importance of this complex for human health.

The Integrator complex desensitizes cellular response to TGF-ß/BMP signaling.

Maintenance of stem cells requires the concerted actions of niche-derived signals and stem cell-intrinsic factors. Although Decapentaplegic (Dpp), a Drosophila bone morphogenetic protein (BMP) molecule, can act as a long-range morphogen, its function is spatially limited to the germline stem cell niche in the germarium. We show here that Integrator, a complex known to be involved in RNA polymerase II (RNAPII)-mediated transcriptional regulation in the nucleus, promotes germline differentiation by restricting niche-derived Dpp/BMP activity in the cytoplasm. Further results show that Integrator works in various developmental contexts to desensitize the cellular response to Dpp/BMP signaling during Drosophila development. Mechanistically, our results show that Integrator forms a multi-subunit complex with the type I receptor Thickveins (Tkv) and other Dpp/BMP signaling components and acts in a negative feedback loop to promote Tkv turnover independent of its transcriptional activity. Similarly, human Integrator subunits bind transforming growth factor ß (TGF-ß)/BMP signaling components and antagonize their activity, suggesting a conserved role of Integrator across metazoans.

INTS13 variants causing a recessive developmental ciliopathy disrupt assembly of the Integrator complex.

Oral-facial-digital (OFD) syndromes are a heterogeneous group of congenital disorders characterized by malformations of the face and oral cavity, and digit anomalies. Mutations within 12 cilia-related genes have been identified that cause several types of OFD, suggesting that OFDs constitute a subgroup of developmental ciliopathies. Through homozygosity mapping and exome sequencing of two families with variable OFD type 2, we identified distinct germline variants in INTS13, a subunit of the Integrator complex. This multiprotein complex associates with RNA Polymerase II and cleaves nascent RNA to modulate gene expression. We determined that INTS13 utilizes its C-terminus to bind the Integrator cleavage module, which is disrupted by the identified germline variants p.S652L and p.K668Nfs*9. Depletion of INTS13 disrupts ciliogenesis in human cultured cells and causes dysregulation of a broad collection of ciliary genes. Accordingly, its knockdown in Xenopus embryos leads to motile cilia anomalies. Altogether, we show that mutations in INTS13 cause an autosomal recessive ciliopathy, which reveals key interactions between components of the Integrator complex.

Integrator endonuclease drives promoter-proximal termination at all RNA polymerase II-transcribed loci.

RNA polymerase II (RNAPII) pausing in early elongation is critical for gene regulation. Paused RNAPII can be released into productive elongation by the kinase P-TEFb or targeted for premature termination by the Integrator complex. Integrator comprises endonuclease and phosphatase activities, driving termination by cleavage of nascent RNA and removal of stimulatory phosphorylation. We generated a degron system for rapid Integrator endonuclease (INTS11) depletion to probe the direct consequences of Integrator-mediated RNA cleavage. Degradation of INTS11 elicits nearly universal increases in active early elongation complexes. However, these RNAPII complexes fail to achieve optimal elongation rates and exhibit persistent Integrator phosphatase activity. Thus, only short transcripts are significantly upregulated following INTS11 loss, including transcription factors, signaling regulators, and non-coding RNAs. We propose a uniform molecular function for INTS11 across all RNAPII-transcribed loci, with differential effects on particular genes, pathways, or RNA biotypes reflective of transcript lengths rather than specificity of Integrator activity.

Inositol hexakisphosphate is required for Integrator function.

Integrator is a multi-subunit protein complex associated with RNA polymerase II (Pol II), with critical roles in noncoding RNA 3'-end processing and transcription attenuation of a broad collection of mRNAs. IntS11 is the endonuclease for RNA cleavage, as a part of the IntS4-IntS9-IntS11 Integrator cleavage module (ICM). Here we report a cryo-EM structure of the Drosophila ICM, at 2.74 Å resolution, revealing stable association of an inositol hexakisphosphate (IP6) molecule. The IP6 binding site is located in a highly electropositive pocket at an interface among all three subunits of ICM, 55 Å away from the IntS11 active site and generally conserved in other ICMs. We also confirmed IP6 association with the same site in human ICM. IP6 binding is not detected in ICM samples harboring mutations in this binding site. Such mutations or disruption of IP6 biosynthesis significantly reduced Integrator function in snRNA 3'-end processing and mRNA transcription attenuation. Our structural and functional studies reveal that IP6 is required for Integrator function in Drosophila, humans, and likely other organisms.

3'UTR shortening of HAS2 promotes hyaluronan hyper-synthesis and bioenergetic dysfunction in pulmonary hypertension.

Pulmonary hypertension (PH) comprises a diverse group of disorders that share a common pathway of pulmonary vascular remodeling leading to right ventricular failure. Development of anti-remodeling strategies is an emerging frontier in PH therapeutics that requires a greater understanding of the interactions between vascular wall cells and their extracellular matrices. The ubiquitous matrix glycan, hyaluronan (HA), is markedly elevated in lungs from patients and experimental models with PH. Herein, we identified HA synthase-2 (HAS2) in the pulmonary artery smooth muscle cell (PASMC) layer as a predominant locus of HA dysregulation. HA upregulation involves depletion of NUDT21, a master regulator of alternative polyadenylation, resulting in 3'UTR shortening and hyper-expression of HAS2. The ensuing increase of HAS2 and hyper-synthesis of HA promoted bioenergetic dysfunction of PASMC characterized by impaired mitochondrial oxidative capacity and a glycolytic shift. The resulting HA accumulation stimulated pro-remodeling phenotypes such as cell proliferation, migration, apoptosis-resistance, and stimulated pulmonary artery contractility. Transgenic mice, mimicking HAS2 hyper-synthesis in smooth muscle cells, developed spontaneous PH, whereas targeted deletion of HAS2 prevented experimental PH. Pharmacological blockade of HAS2 restored normal bioenergetics in PASMC, ameliorated cell remodeling phenotypes, and reversed experimental PH in vivo. In summary, our results uncover a novel mechanism of HA hyper-synthesis and downstream effects on pulmonary vascular cell metabolism and remodeling.

Mapping information-rich genotype-phenotype landscapes with genome-scale Perturb-seq.

A central goal of genetics is to define the relationships between genotypes and phenotypes. High-content phenotypic screens such as Perturb-seq (CRISPR-based screens with single-cell RNA-sequencing readouts) enable massively parallel functional genomic mapping but, to date, have been used at limited scales. Here, we perform genome-scale Perturb-seq targeting all expressed genes with CRISPR interference (CRISPRi) across >2.5 million human cells. We use transcriptional phenotypes to predict the function of poorly characterized genes, uncovering new regulators of ribosome biogenesis (including CCDC86, ZNF236, and SPATA5L1), transcription (C7orf26), and mitochondrial respiration (TMEM242). In addition to assigning gene function, single-cell transcriptional phenotypes allow for in-depth dissection of complex cellular phenomena-from RNA processing to differentiation. We leverage this ability to systematically identify genetic drivers and consequences of aneuploidy and to discover an unanticipated layer of stress-specific regulation of the mitochondrial genome. Our information-rich genotype-phenotype map reveals a multidimensional portrait of gene and cellular function.

Sine oculis homeobox homolog 1 plays a critical role in pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a fatal disease with limited treatment options. The role of the developmental transcription factor Sine oculis homeobox homolog 1 (SIX1) in the pathophysiology of lung fibrosis is not known. IPF lung tissue samples and IPF-derived alveolar type II cells (AT2) showed a significant increase in SIX1 mRNA and protein levels, and the SIX1 transcriptional coactivators EYA1 and EYA2 were elevated. Six1 was also upregulated in bleomycin-treated (BLM-treated) mice and in a model of spontaneous lung fibrosis driven by deletion of Telomeric Repeat Binding Factor 1 (Trf1) in AT2 cells. Conditional deletion of Six1 in AT2 cells prevented or halted BLM-induced lung fibrosis, as measured by a significant reduction in histological burden of fibrosis, reduced fibrotic mediator expression, and improved lung function. These effects were associated with increased macrophage migration inhibitory factor (MIF) in lung epithelial cells in vivo following SIX1 overexpression in BLM-induced fibrosis. A MIF promoter-driven luciferase assay demonstrated direct binding of Six1 to the 5'-TCAGG-3' consensus sequence of the MIF promoter, identifying a likely mechanism of SIX1-driven MIF expression in the pathogenesis of lung fibrosis and providing a potentially novel pathway for targeting in IPF therapy.

Suppression of premature transcription termination leads to reduced mRNA isoform diversity and neurodegeneration.

Tight regulation of mRNA isoform expression is essential for neuronal development, maintenance, and function; however, the repertoire of proteins that govern isoform composition and abundance remains incomplete. Here, we show that the RNA kinase CLP1 regulates mRNA isoform expression through suppression of proximal cleavage and polyadenylation. We found that human stem-cell-derived motor neurons without CLP1 or with the disease-associated CLP1 p.R140H variant had distinct patterns of RNA-polymerase-II-associated cleavage and polyadenylation complex proteins that correlated with polyadenylation site usage. These changes resulted in imbalanced mRNA isoform expression of long genes important for neuronal function that were recapitulated in vivo. Strikingly, we observed the same pattern of reduced mRNA isoform diversity in 3' end sequencing data from brain tissues of patients with neurodegenerative disease. Together, our results identify a previously uncharacterized role for CLP1 in mRNA 3' end formation and reveal an mRNA misprocessing signature in neurodegeneration that may suggest a common mechanism of disease.