Galectin-3 does not interact with RNA directly

Galectin-3, well characterized as a glycan binding protein, has been identified as a putative RNA binding protein, possibly through participation in pre-mRNA maturation through interactions with splicosomes. Given recent developments with cell surface RNA biology, the putative dual-function nature of galectin-3 evokes a possible non-classical connection between glycobiology and RNA biology. However, with limited functional evidence of a direct RNA interaction, many molecular-level observations rely on affinity reagents and lack appropriate genetic controls. Thus, evidence of a direct interaction remains elusive. We demonstrate that antibodies raised to endogenous human galectin-3 can isolate RNA-protein crosslinks, but this activity remains insensitive to LGALS3 knock-out. Proteomic characterization of anti-galectin-3 IPs revealed enrichment of galectin-3, but high abundance of hnRNPA2B1, an abundant, well-characterized RNA-binding protein with weak homology to the N-terminal domain of galectin-3, in the isolate. Genetic ablation of HNRNPA2B1, but not LGALS3, eliminates the ability of the anti-galectin-3 antibodies to isolate RNA-protein crosslinks, implying either an indirect interaction or cross-reactivity. To address this, we introduced an epitope tag to the endogenous C-terminal locus of LGALS3. Isolation of the tagged galectin-3 failed to reveal any RNA-protein crosslinks. This result suggests that the galectin-3 does not directly interact with RNA and may be misidentified as an RNA-binding protein, at least in HeLa where the putative RNA associations were first identified. We encourage further investigation of this phenomenon employ gene deletions and, when possible, endogenous epitope tags to achieve the specificity required to evaluate potential interactions.

Development and application of GlycanDIA workflow for glycomic analysis

Glycans modify protein, lipid, and even RNA molecules to form the regulatory outer coat on cells called the glycocalyx. The changes in glycosylation have been linked to the initiation and progression of many diseases. Thus, while the significance of glycosylation is well established, a lack of accessible methods to characterize glycans has hindered the ability to understand their biological functions. Mass spectrometry (MS)-based methods have generally been at the core of most glycan profiling efforts; however, modern data-independent acquisition (DIA), which could increase sensitivity and simplify workflows, has not been benchmarked for analyzing glycans. Herein, we developed a DIA-based glycomic workflow, termed GlycanDIA, to identify and quantify glycans with high sensitivity and accuracy. The GlycanDIA workflow combined higher energy collisional dissociation (HCD)-MS/MS and staggered windows for glycomic analysis, which facilitates the sensitivity in identification and the accuracy in quantification compared to conventional data-dependent acquisition (DDA)-based glycomics. To facilitate its use, we also developed a generic search engine, GlycanDIA Finder, incorporating an iterative decoy searching for confident glycan identification and quantification from DIA data. The results showed that GlycanDIA can distinguish glycan composition and isomers from N-glycans, O-glycans, and human milk oligosaccharides (HMOs), while it also reveals information on low-abundant modified glycans. With the improved sensitivity, we performed experiments to profile N-glycans from RNA samples, which have been underrepresented due to their low abundance. Using this integrative workflow to unravel the N-glycan profile in cellular and tissue glycoRNA samples, we found that RNA-glycans have specific forms as compared to protein-glycans and are also tissue-specific differences, suggesting distinct functions in biological processes. Overall, GlycanDIA can provide comprehensive information for glycan identification and quantification, enabling researchers to obtain in-depth and refined details on the biological roles of glycosylation.

Nucleolus activity-dependent recruitment and biomolecular condensation by pH sensing

DEAD-box ATPases are major regulators of biomolecular condensates and orchestrate diverse biochemical processes that are critical for the functioning of cells. How DEAD-box proteins are selectively recruited to their respective biomolecular condensates is unknown. We explored this in the context of the nucleolus and DEAD-box protein DDX21. We find that the pH of the nucleolus is intricately linked to the transcriptional activity of the organelle and facilitates the recruitment and condensation of DDX21. We identify an evolutionarily conserved feature of the C terminus of DDX21 responsible for nucleolar localization. This domain is essential for zebrafish development, and its intrinsically disordered and isoelectric properties are necessary and sufficient for the ability of DDX21 to respond to changes in pH and form condensates. Molecularly, the enzymatic activities of poly(ADP-ribose) polymerases contribute to maintaining the nucleolar pH and, consequently, DDX21 recruitment and nucleolar partitioning. These observations reveal an activity-dependent physicochemical mechanism for the selective recruitment of biochemical activities to biomolecular condensates.

RNA Crossing Membranes: Systems and Mechanisms Contextualizing Extracellular RNA and Cell Surface GlycoRNAs

The subcellular localization of a biopolymer often informs its function. RNA is traditionally confined to the cytosolic and nuclear spaces, where it plays critical and conserved roles across nearly all biochemical processes. Our recent observation of cell surface glycoRNAs may further explain the extracellular role of RNA. While cellular membranes are efficient gatekeepers of charged polymers such as RNAs, a large body of research has demonstrated the accumulation of specific RNA species outside of the cell, termed extracellular RNAs (exRNAs). Across various species and forms of life, protein pores have evolved to transport RNA across membranes, thus providing a mechanistic path for exRNAs to achieve their extracellular topology. Here, we review types of exRNAs and the pores capable of RNA transport to provide a logical and testable path toward understanding the biogenesis and regulation of cell surface glycoRNAs.

Co-transcriptional genome surveillance by HUSH is coupled to termination machinery

The HUSH complex recognizes and silences foreign DNA such as viruses, transposons, and transgenes without prior exposure to its targets. Here, we show that endogenous targets of the HUSH complex fall into two distinct classes based on the presence or absence of H3K9me3. These classes are further distinguished by their transposon content and differential response to the loss of HUSH. A de novo genomic rearrangement at the Sox2 locus induces a switch from H3K9me3-independent to H3K9me3-associated HUSH targeting, resulting in silencing. We further demonstrate that HUSH interacts with the termination factor WDR82 and-via its component MPP8-with nascent RNA. HUSH accumulates at sites of high RNAPII occupancy including long exons and transcription termination sites in a manner dependent on WDR82 and CPSF. Together, our results uncover the functional diversity of HUSH targets and show that this vertebrate-specific complex exploits evolutionarily ancient transcription termination machinery for co-transcriptional chromatin targeting and genome surveillance.

LINE-associated cryptic splicing induces dsRNA-mediated interferon response and tumor immunity

RNA splicing plays a critical role in post-transcriptional gene regulation. Exponential expansion of intron length poses a challenge for accurate splicing. Little is known about how cells prevent inadvertent and often deleterious expression of intronic elements due to cryptic splicing. In this study, we identify hnRNPM as an essential RNA binding protein that suppresses cryptic splicing through binding to deep introns, preserving transcriptome integrity. Long interspersed nuclear elements (LINEs) harbor large amounts of pseudo splice sites in introns. hnRNPM preferentially binds at intronic LINEs and represses LINE-containing pseudo splice site usage for cryptic splicing. Remarkably, a subgroup of the cryptic exons can form long dsRNAs through base-pairing of inverted Alu transposable elements scattered in between LINEs and trigger interferon immune response, a well-known antiviral defense mechanism. Notably, these interferon-associated pathways are found to be upregulated in hnRNPM-deficient tumors, which also exhibit elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity. Targeting hnRNPM in tumors may be used to trigger an inflammatory immune response thereby boosting cancer surveillance.

Glucose dissociates DDX21 dimers to regulate mRNA splicing and tissue differentiation

Glucose is a universal bioenergy source; however, its role in controlling protein interactions is unappreciated, as are its actions during differentiation-associated intracellular glucose elevation. Azido-glucose click chemistry identified glucose binding to a variety of RNA binding proteins (RBPs), including the DDX21 RNA helicase, which was found to be essential for epidermal differentiation. Glucose bound the ATP-binding domain of DDX21, altering protein conformation, inhibiting helicase activity, and dissociating DDX21 dimers. Glucose elevation during differentiation was associated with DDX21 re-localization from the nucleolus to the nucleoplasm where DDX21 assembled into larger protein complexes containing RNA splicing factors. DDX21 localized to specific SCUGSDGC motif in mRNA introns in a glucose-dependent manner and promoted the splicing of key pro-differentiation genes, including GRHL3, KLF4, OVOL1, and RBPJ. These findings uncover a biochemical mechanism of action for glucose in modulating the dimerization and function of an RNA helicase essential for tissue differentiation.

Nuclear accumulation of host transcripts during Zika Virus Infection

Zika virus (ZIKV) infects fetal neural progenitor cells (NPCs) causing severe neurodevelopmental disorders in utero. Multiple pathways involved in normal brain development are dysfunctional in infected NPCs but how ZIKV centrally reprograms these pathways remains unknown. Here we show that ZIKV infection disrupts subcellular partitioning of host transcripts critical for neurodevelopment in NPCs and functionally link this process to the up-frameshift protein 1 (UPF1). UPF1 is an RNA-binding protein known to regulate decay of cellular and viral RNAs and is less expressed in ZIKV-infected cells. Using infrared crosslinking immunoprecipitation and RNA sequencing (irCLIP-Seq), we show that a subset of mRNAs loses UPF1 binding in ZIKV-infected NPCs, consistent with UPF1's diminished expression. UPF1 target transcripts, however, are not altered in abundance but in subcellular localization, with mRNAs accumulating in the nucleus of infected or UPF1 knockdown cells. This leads to diminished protein expression of FREM2, a protein required for maintenance of NPC identity. Our results newly link UPF1 to the regulation of mRNA transport in NPCs, a process perturbed during ZIKV infection.

The human disease gene LYSET is essential for lysosomal enzyme transport and viral infection

Lysosomes are key degradative compartments of the cell. Transport to lysosomes relies on GlcNAc-1-phosphotransferase-mediated tagging of soluble enzymes with mannose 6-phosphate (M6P). GlcNAc-1-phosphotransferase deficiency leads to the severe lysosomal storage disorder mucolipidosis II (MLII). Several viruses require lysosomal cathepsins to cleave structural proteins and thus depend on functional GlcNAc-1-phosphotransferase. Here, we used genome-scale CRISPR screens to identify Lysosomal Enzyme Trafficking factor (LYSET) as essential for infection by cathepsin-dependent viruses including SARS-CoV-2. LYSET deficiency resulted in global loss of M6P tagging and mislocalization of GlcNAc-1-phosphotransferase from the Golgi complex to lysosomes. Lyset knockout mice exhibited MLII-like phenotypes and human pathogenic LYSET alleles failed to restore lysosomal sorting defects. Thus, LYSET is required for correct functioning of the M6P trafficking machinery, and mutations in LYSET can explain the phenotype of the associated disorder.

An atlas of posttranslational modifications on RNA binding proteins

RNA structure and function are intimately tied to RNA binding protein recognition and regulation. Posttranslational modifications are chemical modifications which can control protein biology. The role of PTMs in the regulation RBPs is not well understood, in part due to a lacking analysis of PTM deposition on RBPs. Herein, we present an analysis of posttranslational modifications (PTMs) on RNA binding proteins (RBPs; a PTM RBP Atlas). We curate published datasets and primary literature to understand the landscape of PTMs and use protein-protein interaction data to understand and potentially provide a framework for understanding which enzymes are controlling PTM deposition and removal on the RBP landscape. Intersection of our data with The Cancer Genome Atlas also provides researchers understanding of mutations that would alter PTM deposition. Additional characterization of the RNA-protein interface provided from in-cell UV crosslinking experiments provides a framework for hypotheses about which PTMs could be regulating RNA binding and thus RBP function. Finally, we provide an online database for our data that is easy to use for the community. It is our hope our efforts will provide researchers will an invaluable tool to test the function of PTMs controlling RBP function and thus RNA biology.

The CD22-IGF2R interaction is a therapeutic target for microglial lysosome dysfunction in Niemann-Pick type C

Lysosome dysfunction is a shared feature of rare lysosomal storage diseases and common age-related neurodegenerative diseases. Microglia, the brain-resident macrophages, are particularly vulnerable to lysosome dysfunction because of the phagocytic stress of clearing dying neurons, myelin, and debris. CD22 is a negative regulator of microglial homeostasis in the aging mouse brain, and soluble CD22 (sCD22) is increased in the cerebrospinal fluid of patients with Niemann-Pick type C disease (NPC). However, the role of CD22 in the human brain remains unknown. In contrast to previous findings in mice, here, we show that CD22 is expressed by oligodendrocytes in the human brain and binds to sialic acid–dependent ligands on microglia. Using unbiased genetic and proteomic screens, we identify insulin-like growth factor 2 receptor (IGF2R) as the binding partner of sCD22 on human myeloid cells. Targeted truncation of IGF2R revealed that sCD22 docks near critical mannose 6-phosphate–binding domains, where it disrupts lysosomal protein trafficking. Interfering with the sCD22-IGF2R interaction using CD22 blocking antibodies ameliorated lysosome dysfunction in human NPC1 mutant induced pluripotent stem cell–derived microglia-like cells without harming oligodendrocytes in vitro. These findings reinforce the differences between mouse and human microglia and provide a candidate microglia-directed immunotherapeutic to treat NPC.

fSHAPE, fSHAPE-eCLIP, and SHAPE-eCLIP probe transcript regions that interact with specific proteins

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) structure probing techniques characterize the secondary structure of RNA molecules, which influence their functions and interactions. A variation of SHAPE, footprinting SHAPE (fSHAPE), probes RNA in the presence and absence of protein to identify RNA bases that hydrogen-bond with protein. SHAPE or fSHAPE coupled with enhanced crosslinking and immunoprecipitation (SHAPE-eCLIP or fSHAPE-eCLIP) pulls down RNAs bound by any protein of interest and returns their structure or protein interaction information, respectively. Here, we describe detailed protocols for SHAPE-eCLIP and fSHAPE-eCLIP and an analysis protocol for fSHAPE. For complete details on the use and execution of these protocols, please refer to Corley et al. (2020).

Small RNAs are modified with N-glycans and displayed on the surface of living cells

Glycans modify lipids and proteins to mediate inter- and intramolecular interactions across all domains of life. RNA is not thought to be a major target of glycosylation. Here, we challenge this view with evidence that mammals use RNA as a third scaffold for glycosylation. Using a battery of chemical and biochemical approaches, we found that conserved small noncoding RNAs bear sialylated glycans. These "glycoRNAs" were present in multiple cell types and mammalian species, in cultured cells, and in vivo. GlycoRNA assembly depends on canonical N-glycan biosynthetic machinery and results in structures enriched in sialic acid and fucose. Analysis of living cells revealed that the majority of glycoRNAs were present on the cell surface and can interact with anti-dsRNA antibodies and members of the Siglec receptor family. Collectively, these findings suggest the existence of a direct interface between RNA biology and glycobiology, and an expanded role for RNA in extracellular biology.

Discovery and functional interrogation of SARS-CoV-2 RNA-host protein interactions

SARS-CoV-2 is the cause of a pandemic with growing global mortality. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with ChIRP-MS data from three other RNA viruses defined viral specificity of RNA-host protein interactions. Targeted CRISPR screens revealed that the majority of functional RNA-binding proteins protect the host from virus-induced cell death, and comparative CRISPR screens across seven RNA viruses revealed shared and SARS-specific antiviral factors. Finally, by combining the RNA-centric approach and functional CRISPR screens, we demonstrated a physical and functional connection between SARS-CoV-2 and mitochondria, highlighting this organelle as a general platform for antiviral activity. Altogether, these data provide a comprehensive catalog of functional SARS-CoV-2 RNA-host protein interactions, which may inform studies to understand the host-virus interface and nominate host pathways that could be targeted for therapeutic benefit.

Footprinting SHAPE-eCLIP Reveals Transcriptome-wide Hydrogen Bonds at RNA-Protein Interfaces

Discovering the interaction mechanism and location of RNA-binding proteins (RBPs) on RNA is critical for understanding gene expression regulation. Here, we apply selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) on in vivo transcripts compared to protein-absent transcripts in four human cell lines to identify transcriptome-wide footprints (fSHAPE) on RNA. Structural analyses indicate that fSHAPE precisely detects nucleobases that hydrogen bond with protein. We demonstrate that fSHAPE patterns predict binding sites of known RBPs, such as iron response elements in both known loci and previously unknown loci in CDC34, SLC2A4RG, COASY, and H19. Furthermore, by integrating SHAPE and fSHAPE with crosslinking and immunoprecipitation (eCLIP) of desired RBPs, we interrogate specific RNA-protein complexes, such as histone stem-loop elements and their nucleotides that hydrogen bond with stem-loop-binding proteins. Together, these technologies greatly expand our ability to study and understand specific cellular RNA interactions in RNA-protein complexes.

Progenitor identification and SARS-CoV-2 infection in human distal lung organoids

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange. Three-dimensional in vitro human distal lung culture systems would strongly facilitate the investigation of pathologies such as interstitial lung disease, cancer and coronavirus disease 2019 (COVID-19) pneumonia caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here we describe the development of a long-term feeder-free, chemically defined culture system for distal lung progenitors as organoids derived from single adult human alveolar epithelial type II (AT2) or KRT5+ basal cells. AT2 organoids were able to differentiate into AT1 cells, and basal cell organoids developed lumens lined with differentiated club and ciliated cells. Single-cell analysis of KRT5+ cells in basal organoids revealed a distinct population of ITGA6+ITGB4+ mitotic cells, whose offspring further segregated into a TNFRSF12Ahi subfraction that comprised about ten per cent of KRT5+ basal cells. This subpopulation formed clusters within terminal bronchioles and exhibited enriched clonogenic organoid growth activity. We created distal lung organoids with apical-out polarity to present ACE2 on the exposed external surface, facilitating infection of AT2 and basal cultures with SARS-CoV-2 and identifying club cells as a target population. This long-term, feeder-free culture of human distal lung organoids, coupled with single-cell analysis, identifies functional heterogeneity among basal cells and establishes a facile in vitro organoid model of human distal lung infections, including COVID-19-associated pneumonia.

Proteomic analysis of young and old mouse hematopoietic stem cells and their progenitors reveals post-transcriptional regulation in stem cells

The balance of hematopoietic stem cell (HSC) self-renewal and differentiation is critical for a healthy blood supply; imbalances underlie hematological diseases. The importance of HSCs and their progenitors have led to their extensive characterization at genomic and transcriptomic levels. However, the proteomics of hematopoiesis remains incompletely understood. Here we report a proteomics resource from mass spectrometry of mouse young adult and old adult mouse HSCs, multipotent progenitors and oligopotent progenitors; 12 cell types in total. We validated differential protein levels, including confirmation that Dnmt3a protein levels are undetected in young adult mouse HSCs until forced into cycle. Additionally, through integrating proteomics and RNA-sequencing datasets, we identified a subset of genes with apparent post-transcriptional repression in young adult mouse HSCs. In summary, we report proteomic coverage of young and old mouse HSCs and progenitors, with broader implications for understanding mechanisms for stem cell maintenance, niche interactions and fate determination.

Systematic discovery and functional interrogation of SARS-CoV-2 viral RNA-host protein interactions during infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of a pandemic with growing global mortality. There is an urgent need to understand the molecular pathways required for host infection and anti-viral immunity. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with viral ChIRP-MS data from three other positive-sense RNA viruses defined pan-viral and SARS-CoV-2-specific host interactions. Functional interrogation of these factors with a genome-wide CRISPR screen revealed that the vast majority of viral RNA-binding proteins protect the host from virus-induced cell death, and we identified known and novel anti-viral proteins that regulate SARS-CoV-2 pathogenicity. Finally, our RNA-centric approach demonstrated a physical connection between SARS-CoV-2 RNA and host mitochondria, which we validated with functional and electron microscopy data, providing new insights into a more general virus-specific protein logic for mitochondrial interactions. Altogether, these data provide a comprehensive catalogue of SARS-CoV-2 RNA-host protein interactions, which may inform future studies to understand the mechanisms of viral pathogenesis, as well as nominate host pathways that could be targeted for therapeutic benefit.

HIGHLIGHTS

· ChIRP-MS of SARS-CoV-2 RNA identifies a comprehensive viral RNA-host protein interaction network during infection across two species· Comparison to RNA-protein interaction networks with Zika virus, dengue virus, and rhinovirus identify SARS-CoV-2-specific and pan-viral RNA protein complexes and highlights distinct intracellular trafficking pathways· Intersection of ChIRP-MS and genome-wide CRISPR screens identify novel SARS-CoV-2-binding proteins with pro- and anti-viral function· Viral RNA-RNA and RNA-protein interactions reveal specific SARS-CoV-2-mediated mitochondrial dysfunction during infection.

A memory of eS25 loss drives resistance phenotypes

In order to maintain cellular protein homeostasis, ribosomes are safeguarded against dysregulation by myriad processes. Remarkably, many cell types can withstand genetic lesions of certain ribosomal protein genes, some of which are linked to diverse cellular phenotypes and human disease. Yet the direct and indirect consequences from these lesions are poorly understood. To address this knowledge gap, we studied in vitro and cellular consequences that follow genetic knockout of the ribosomal proteins RPS25 or RACK1 in a human cell line, as both proteins are implicated in direct translational control. Prompted by the unexpected detection of an off-target ribosome alteration in the RPS25 knockout, we closely interrogated cellular phenotypes. We found that multiple RPS25 knockout clones display viral- and toxin-resistance phenotypes that cannot be rescued by functional cDNA expression, suggesting that RPS25 loss elicits a cell state transition. We characterized this state and found that it underlies pleiotropic phenotypes and has a common rewiring of gene expression. Rescuing RPS25 expression by genomic locus repair failed to correct for the phenotypic and expression hysteresis. Our findings illustrate how the elasticity of cells to a ribosome perturbation can drive specific phenotypic outcomes that are indirectly linked to translation and suggests caution in the interpretation of ribosomal protein gene mutation data.

Progenitor identification and SARS-CoV-2 infection in long-term human distal lung organoid cultures

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange and is affected by disorders including interstitial lung disease, cancer, and SARS-CoV-2-associated COVID-19 pneumonia. Investigations of these localized pathologies have been hindered by a lack of 3D in vitro human distal lung culture systems. Further, human distal lung stem cell identification has been impaired by quiescence, anatomic divergence from mouse and lack of lineage tracing and clonogenic culture. Here, we developed robust feeder-free, chemically-defined culture of distal human lung progenitors as organoids derived clonally from single adult human alveolar epithelial type II (AT2) or KRT5 + basal cells. AT2 organoids exhibited AT1 transdifferentiation potential, while basal cell organoids progressively developed lumens lined by differentiated club and ciliated cells. Organoids consisting solely of club cells were not observed. Upon single cell RNA-sequencing (scRNA-seq), alveolar organoids were composed of proliferative AT2 cells; however, basal organoid KRT5 + cells contained a distinct ITGA6 + ITGB4 + mitotic population whose proliferation segregated to a TNFRSF12A hi subfraction. Clonogenic organoid growth was markedly enriched within the TNFRSF12A hi subset of FACS-purified ITGA6 + ITGB4 + basal cells from human lung or derivative organoids. In vivo, TNFRSF12A + cells comprised ~10% of KRT5 + basal cells and resided in clusters within terminal bronchioles. To model COVID-19 distal lung disease, we everted the polarity of basal and alveolar organoids to rapidly relocate differentiated club and ciliated cells from the organoid lumen to the exterior surface, thus displaying the SARS-CoV-2 receptor ACE2 on the outwardly-facing apical aspect. Accordingly, basal and AT2 apical-out organoids were infected by SARS-CoV-2, identifying club cells as a novel target population. This long-term, feeder-free organoid culture of human distal lung alveolar and basal stem cells, coupled with single cell analysis, identifies unsuspected basal cell functional heterogeneity and exemplifies progenitor identification within a slowly proliferating human tissue. Further, our studies establish a facile in vitro organoid model for human distal lung infectious diseases including COVID-19-associated pneumonia.

RNA helicase DDX21 mediates nucleotide stress responses in neural crest and melanoma cells

The availability of nucleotides has a direct impact on transcription. The inhibition of dihydroorotate dehydrogenase (DHODH) with leflunomide impacts nucleotide pools by reducing pyrimidine levels. Leflunomide abrogates the effective transcription elongation of genes required for neural crest development and melanoma growth in vivo¹. To define the mechanism of action, we undertook an in vivo chemical suppressor screen for restoration of neural crest after leflunomide treatment. Surprisingly, we found that alterations in progesterone and progesterone receptor (Pgr) signalling strongly suppressed leflunomide-mediated neural crest effects in zebrafish. In addition, progesterone bypasses the transcriptional elongation block resulting from Paf complex deficiency, rescuing neural crest defects in ctr9 morphant and paf1(aln^z24) mutant embryos. Using proteomics, we found that Pgr binds the RNA helicase protein Ddx21. ddx21-deficient zebrafish show resistance to leflunomide-induced stress. At a molecular level, nucleotide depletion reduced the chromatin occupancy of DDX21 in human A375 melanoma cells. Nucleotide supplementation reversed the gene expression signature and DDX21 occupancy changes prompted by leflunomide. Together, our results show that DDX21 acts as a sensor and mediator of transcription during nucleotide stress.

DNA-PKcs has KU-dependent function in rRNA processing and haematopoiesis

The DNA-dependent protein kinase (DNA-PK), which comprises the KU heterodimer and a catalytic subunit (DNA-PKcs), is a classical non-homologous end-joining (cNHEJ) factor1. KU binds to DNA ends, initiates cNHEJ, and recruits and activates DNA-PKcs. KU also binds to RNA, but the relevance of this interaction in mammals is unclear. Here we use mouse models to show that DNA-PK has an unexpected role in the biogenesis of ribosomal RNA (rRNA) and in haematopoiesis. The expression of kinase-dead DNA-PKcs abrogates cNHEJ2. However, most mice that both expressed kinase-dead DNA-PKcs and lacked the tumour suppressor TP53 developed myeloid disease, whereas all other previously characterized mice deficient in both cNHEJ and TP53 expression succumbed to pro-B cell lymphoma3. DNA-PK autophosphorylates DNA-PKcs, which is its best characterized substrate. Blocking the phosphorylation of DNA-PKcs at the T2609 cluster, but not the S2056 cluster, led to KU-dependent defects in 18S rRNA processing, compromised global protein synthesis in haematopoietic cells and caused bone marrow failure in mice. KU drives the assembly of DNA-PKcs on a wide range of cellular RNAs, including the U3 small nucleolar RNA, which is essential for processing of 18S rRNA4. U3 activates purified DNA-PK and triggers phosphorylation of DNA-PKcs at T2609. DNA-PK, but not other cNHEJ factors, resides in nucleoli in an rRNA-dependent manner and is co-purified with the small subunit processome. Together our data show that DNA-PK has RNA-dependent, cNHEJ-independent functions during ribosome biogenesis that require the kinase activity of DNA-PKcs and its phosphorylation at the T2609 cluster.

Functional significance of U2AF1 S34F mutations in lung adenocarcinomas

The functional role of U2AF1 mutations in lung adenocarcinomas (LUADs) remains incompletely understood. Here, we report a significant co-occurrence of U2AF1 S34F mutations with ROS1 translocations in LUADs. To characterize this interaction, we profiled effects of S34F on the transcriptome-wide distribution of RNA binding and alternative splicing in cells harboring the ROS1 translocation. Compared to its wild-type counterpart, U2AF1 S34F preferentially binds and modulates splicing of introns containing CAG trinucleotides at their 3' splice junctions. The presence of S34F caused a shift in cross-linking at 3' splice sites, which was significantly associated with alternative splicing of skipped exons. U2AF1 S34F induced expression of genes involved in the epithelial-mesenchymal transition (EMT) and increased tumor cell invasion. Finally, S34F increased splicing of the long over the short SLC34A2-ROS1 isoform, which was also associated with enhanced invasiveness. Taken together, our results suggest a mechanistic interaction between mutant U2AF1 and ROS1 in LUAD.

HiChIRP reveals RNA-associated chromosome conformation

Modular domains of long non-coding RNAs can serve as scaffolds to bring distant regions of the linear genome into spatial proximity. Here, we present HiChIRP, a method leveraging bio-orthogonal chemistry and optimized chromosome conformation capture conditions, which enables interrogation of chromatin architecture focused around a specific RNA of interest down to approximately ten copies per cell. HiChIRP of three nuclear RNAs reveals insights into promoter interactions (7SK), telomere biology (telomerase RNA component) and inflammatory gene regulation (lincRNA-EPS).

The novel lncRNA lnc-NR2F1 is pro-neurogenic and mutated in human neurodevelopmental disorders

Long noncoding RNAs (lncRNAs) have been shown to act as important cell biological regulators including cell fate decisions but are often ignored in human genetics. Combining differential lncRNA expression during neuronal lineage induction with copy number variation morbidity maps of a cohort of children with autism spectrum disorder/intellectual disability versus healthy controls revealed focal genomic mutations affecting several lncRNA candidate loci. Here we find that a t(5:12) chromosomal translocation in a family manifesting neurodevelopmental symptoms disrupts specifically lnc-NR2F1. We further show that lnc-NR2F1 is an evolutionarily conserved lncRNA functionally enhances induced neuronal cell maturation and directly occupies and regulates transcription of neuronal genes including autism-associated genes. Thus, integrating human genetics and functional testing in neuronal lineage induction is a promising approach for discovering candidate lncRNAs involved in neurodevelopmental diseases.

Medulloblastoma-associated DDX3 variant selectively alters the translational response to stress

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3′s interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5′UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3′s role in this process. Arsenite-induced stress shifts DDX3 binding from the 5′UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3^R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3^R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.

Inhibition of NGLY1 Inactivates the Transcription Factor Nrf1 and Potentiates Proteasome Inhibitor Cytotoxicity

Proteasome inhibitors are used to treat blood cancers such as multiple myeloma (MM) and mantle cell lymphoma. The efficacy of these drugs is frequently undermined by acquired resistance. One mechanism of proteasome inhibitor resistance may involve the transcription factor Nuclear Factor, Erythroid 2 Like 1 (NFE2L1, also referred to as Nrf1), which responds to proteasome insufficiency or pharmacological inhibition by upregulating proteasome subunit gene expression. This "bounce-back" response is achieved through a unique mechanism. Nrf1 is constitutively translocated into the ER lumen, N-glycosylated, and then targeted for proteasomal degradation via the ER-associated degradation (ERAD) pathway. Proteasome inhibition leads to accumulation of cytosolic Nrf1, which is then processed to form the active transcription factor. Here we show that the cytosolic enzyme N-glycanase 1 (NGLY1, the human PNGase) is essential for Nrf1 activation in response to proteasome inhibition. Chemical or genetic disruption of NGLY1 activity results in the accumulation of misprocessed Nrf1 that is largely excluded from the nucleus. Under these conditions, Nrf1 is inactive in regulating proteasome subunit gene expression in response to proteasome inhibition. Through a small molecule screen, we identified a cell-active NGLY1 inhibitor that disrupts the processing and function of Nrf1. The compound potentiates the cytotoxicity of carfilzomib, a clinically used proteasome inhibitor, against MM and T cell-derived acute lymphoblastic leukemia (T-ALL) cell lines. Thus, NGLY1 inhibition prevents Nrf1 activation and represents a new therapeutic approach for cancers that depend on proteasome homeostasis.

Enhancer connectome in primary human cells identifies target genes of disease-associated DNA elements

The challenge of linking intergenic mutations to target genes has limited molecular understanding of human diseases. Here we show that H3K27ac HiChIP generates high-resolution contact maps of active enhancers and target genes in rare primary human T cell subtypes and coronary artery smooth muscle cells. Differentiation of naive T cells into T helper 17 cells or regulatory T cells creates subtype-specific enhancer-promoter interactions, specifically at regions of shared DNA accessibility. These data provide a principled means of assigning molecular functions to autoimmune and cardiovascular disease risk variants, linking hundreds of noncoding variants to putative gene targets. Target genes identified with HiChIP are further supported by CRISPR interference and activation at linked enhancers, by the presence of expression quantitative trait loci, and by allele-specific enhancer loops in patient-derived primary cells. The majority of disease-associated enhancers contact genes beyond the nearest gene in the linear genome, leading to a fourfold increase in the number of potential target genes for autoimmune and cardiovascular diseases.

The Mammalian Ribo-interactome Reveals Ribosome Functional Diversity and Heterogeneity

During eukaryotic evolution, ribosomes have considerably increased in size, forming a surface-exposed ribosomal RNA (rRNA) shell of unknown function, which may create an interface for yet uncharacterized interacting proteins. To investigate such protein interactions, we establish a ribosome affinity purification method that unexpectedly identifies hundreds of ribosome-associated proteins (RAPs) from categories including metabolism and cell cycle, as well as RNA- and protein-modifying enzymes that functionally diversify mammalian ribosomes. By further characterizing RAPs, we discover the presence of ufmylation, a metazoan-specific post-translational modification (PTM), on ribosomes and define its direct substrates. Moreover, we show that the metabolic enzyme, pyruvate kinase muscle (PKM), interacts with sub-pools of endoplasmic reticulum (ER)-associated ribosomes, exerting a non-canonical function as an RNA-binding protein in the translation of ER-destined mRNAs. Therefore, RAPs interconnect one of life's most ancient molecular machines with diverse cellular processes, providing an additional layer of regulatory potential to protein expression.

Comparative Analysis Reveals Furoyl in Vivo Selective Hydroxyl Acylation Analyzed by Primer Extension Reagents Form Stable Ribosyl Ester Adducts

RNA molecules depend on structural elements that are critical for cellular function. Chemical methods for probing RNA structure have emerged as a necessary component of characterizing RNA function. As such, understanding the limitations and idiosyncrasies of these methods is essential for their utility. Selective hydroxyl acylation has emerged as a common method for analyzing RNA structure. Ester products as a result of 2'-hydroxyl acylation can then be identified through reverse transcription or mutational enzyme profiling. The central aspect of selective hydroxyl acylation analyzed by primer extension (SHAPE) experiments is the fact that stable ester adducts are formed on the 2'-hydroxyl. Despite its importance, there has not been a direct comparison of SHAPE electrophiles for their ability to make stable RNA adducts. Herein, we conduct a systematic analysis of hydrolysis stability experiments to demonstrate that furoyl imidazole SHAPE reagents form stable ester adducts even at elevated temperatures. We also demonstrate that the acylation reaction with the furoyl acylimidaole SHAPE reagent can be controlled with dithiothreitol quenching, even in live cells. These results are important for our understanding of the biochemical details of the SHAPE experiment.

Comparison of SHAPE reagents for mapping RNA structures inside living cells

Recent advances in SHAPE technology have converted the classic primer extension method to next-generation sequencing platforms, allowing transcriptome-level analysis of RNA secondary structure. In particular, icSHAPE and SHAPE-MaP, using NAI-N₃ and 1M7 reagents, respectively, are methods that claim to measure in vivo structure with high-throughput sequencing. However, these compounds have not been compared on an unbiased, raw-signal level. Here, we directly compare several in vivo SHAPE acylation reagents using the simple primer extension assay. We conclude that while multiple SHAPE technologies are effective at measuring purified RNAs in vitro, acylimidazole reagents NAI and NAI-N₃ give markedly greater signals with lower background than 1M7 for in vivo measurement of the RNA structurome.

HiChIP: efficient and sensitive analysis of protein-directed genome architecture

Genome conformation is central to gene control but challenging to interrogate. Here we present HiChIP, a protein-centric chromatin conformation method. HiChIP improves the yield of conformation-informative reads by over 10-fold and lowers the input requirement over 100-fold relative to that of ChIA-PET. HiChIP of cohesin reveals multiscale genome architecture with greater signal-to-background ratios than those of in situ Hi-C.

An inducible long noncoding RNA amplifies DNA damage signaling

Long noncoding RNAs (lncRNAs) are prevalent genes with frequently precise regulation but mostly unknown functions. Here we demonstrate that lncRNAs guide the organismal DNA damage response. DNA damage activated transcription of the DINO (Damage Induced Noncoding) lncRNA via p53. DINO was required for p53-dependent gene expression, cell cycle arrest and apoptosis in response to DNA damage, and DINO expression was sufficient to activate damage signaling and cell cycle arrest in the absence of DNA damage. DINO bound to p53 protein and promoted its stabilization, mediating a p53 auto-amplification loop. Dino knockout or promoter inactivation in mice dampened p53 signaling and ameliorated acute radiation syndrome in vivo. Thus, inducible lncRNA can create a feedback loop with its cognate transcription factor to amplify cellular signaling networks.

Stress from Nucleotide Depletion Activates the Transcriptional Regulator HEXIM1 to Suppress Melanoma

Studying cancer metabolism gives insight into tumorigenic survival mechanisms and susceptibilities. In melanoma, we identify HEXIM1, a transcription elongation regulator, as a melanoma tumor suppressor that responds to nucleotide stress. HEXIM1 expression is low in melanoma. Its overexpression in a zebrafish melanoma model suppresses cancer formation, while its inactivation accelerates tumor onset in vivo. Knockdown of HEXIM1 rescues zebrafish neural crest defects and human melanoma proliferation defects that arise from nucleotide depletion. Under nucleotide stress, HEXIM1 is induced to form an inhibitory complex with P-TEFb, the kinase that initiates transcription elongation, to inhibit elongation at tumorigenic genes. The resulting alteration in gene expression also causes anti-tumorigenic RNAs to bind to and be stabilized by HEXIM1. HEXIM1 plays an important role in inhibiting cancer cell-specific gene transcription while also facilitating anti-cancer gene expression. Our study reveals an important role for HEXIM1 in coupling nucleotide metabolism with transcriptional regulation in melanoma.

7SK-BAF axis controls pervasive transcription at enhancers

RNA functions at enhancers remain mysterious. Here we show that the 7SK small nuclear RNA (snRNA) inhibits enhancer transcription by modulating nucleosome position. 7SK occupies enhancers and super enhancers genome-wide in mouse and human cells, and 7SK is required to limit eRNA initiation and synthesis in a manner distinct from promoter pausing. Clustered elements at super enhancers uniquely require 7SK to prevent convergent transcription and DNA damage signaling. 7SK physically interacts with the BAF chromatin remodeling complex, recruit BAF to enhancers, and inhibits enhancer transcription by modulating chromatin structure. In turn, 7SK occupancy at enhancers coincides with Brd4 and is exquisitely sensitive to the bromodomain inhibitor JQ1. Thus, 7SK employs distinct mechanisms to counteract diverse consequences of pervasive transcription that distinguish super enhancers, enhancers, and promoters.

Leukemia-Associated Cohesin Mutants Dominantly Enforce Stem Cell Programs and Impair Human Hematopoietic Progenitor Differentiation

Recurrent mutations in cohesin complex proteins have been identified in pre-leukemic hematopoietic stem cells and during the early development of acute myeloid leukemia and other myeloid malignancies. Although cohesins are involved in chromosome separation and DNA damage repair, cohesin complex functions during hematopoiesis and leukemic development are unclear. Here, we show that mutant cohesin proteins block differentiation of human hematopoietic stem and progenitor cells (HSPCs) in vitro and in vivo and enforce stem cell programs. These effects are restricted to immature HSPC populations, where cohesin mutants show increased chromatin accessibility and likelihood of transcription factor binding site occupancy by HSPC regulators including ERG, GATA2, and RUNX1, as measured by ATAC-seq and ChIP-seq. Epistasis experiments show that silencing these transcription factors rescues the differentiation block caused by cohesin mutants. Together, these results show that mutant cohesins impair HSPC differentiation by controlling chromatin accessibility and transcription factor activity, possibly contributing to leukemic disease.

Long noncoding RNAs in cell-fate programming and reprogramming

In recent years, long noncoding RNAs (lncRNAs) have emerged as an important class of regulators of gene expression. lncRNAs exhibit several distinctive features that confer unique regulatory functions, including exquisite cell- and tissue-specific expression and the capacity to transduce higher-order spatial information. Here we review evidence showing that lncRNAs exert critical functions in adult tissue stem cells, including skin, brain, and muscle, as well as in developmental patterning and pluripotency. We highlight new approaches for ascribing lncRNA functions and discuss mammalian dosage compensation as a classic example of an lncRNA network coupled to stem cell differentiation.

Structural imprints in vivo decode RNA regulatory mechanisms

Visualizing the physical basis for molecular behavior inside living cells is a grand challenge in biology. RNAs are central to biological regulation, and RNA’s ability to adopt specific structures intimately controls every step of the gene expression program¹. However, our understanding of physiological RNA structures is limited; current in vivo RNA structure profiles view only two of four nucleotides that make up RNA^2,3. Here we present a novel biochemical approach, In Vivo Click SHAPE (icSHAPE), that enables the first global view of RNA secondary structures of all four bases in living cells. icSHAPE of mouse embryonic stem cell transcriptome versus purified RNA folded in vitro shows that the structural dynamics of RNA in the cellular environment distinguishes different classes of RNAs and regulatory elements. Structural signatures at translational start sites and ribosome pause sites are conserved from in vitro, suggesting that these RNA elements are programmed by sequence. In contrast, focal structural rearrangements in vivo reveal precise interfaces of RNA with RNA binding proteins or RNA modification sites that are consistent with atomic-resolution structural data. Such dynamic structural footprints enable accurate prediction of RNA-protein interactions and N6-methyladenosine (m⁶A) modification genome-wide. These results open the door for structural genomics of RNA in living cells and reveal key physiological structures controlling gene expression.

Dissecting noncoding and pathogen RNA-protein interactomes

RNA-protein interactions are central to biological regulation. Cross-linking immunoprecipitation (CLIP)-seq is a powerful tool for genome-wide interrogation of RNA-protein interactomes, but current CLIP methods are limited by challenging biochemical steps and fail to detect many classes of noncoding and nonhuman RNAs. Here we present FAST-iCLIP, an integrated pipeline with improved CLIP biochemistry and an automated informatic pipeline for comprehensive analysis across protein coding, noncoding, repetitive, retroviral, and nonhuman transcriptomes. FAST-iCLIP of Poly-C binding protein 2 (PCBP2) showed that PCBP2-bound CU-rich motifs in different topologies to recognize mRNAs and noncoding RNAs with distinct biological functions. FAST-iCLIP of PCBP2 in hepatitis C virus-infected cells enabled a joint analysis of the PCBP2 interactome with host and viral RNAs and their interplay. These results show that FAST-iCLIP can be used to rapidly discover and decipher mechanisms of RNA-protein recognition across the diversity of human and pathogen RNAs.

RNA helicase DDX21 coordinates transcription and ribosomal RNA processing

DEAD-box RNA helicases are vital for the regulation of various aspects of the RNA life cycle, but the molecular underpinnings of their involvement, particularly in mammalian cells, remain poorly understood. Here we show that the DEAD-box RNA helicase DDX21 can sense the transcriptional status of both RNA polymerase (Pol) I and II to control multiple steps of ribosome biogenesis in human cells. We demonstrate that DDX21 widely associates with Pol I- and Pol II-transcribed genes and with diverse species of RNA, most prominently with non-coding RNAs involved in the formation of ribonucleoprotein complexes, including ribosomal RNA, small nucleolar RNAs (snoRNAs) and 7SK RNA. Although broad, these molecular interactions, both at the chromatin and RNA level, exhibit remarkable specificity for the regulation of ribosomal genes. In the nucleolus, DDX21 occupies the transcribed rDNA locus, directly contacts both rRNA and snoRNAs, and promotes rRNA transcription, processing and modification. In the nucleoplasm, DDX21 binds 7SK RNA and, as a component of the 7SK small nuclear ribonucleoprotein (snRNP) complex, is recruited to the promoters of Pol II-transcribed genes encoding ribosomal proteins and snoRNAs. Promoter-bound DDX21 facilitates the release of the positive transcription elongation factor b (P-TEFb) from the 7SK snRNP in a manner that is dependent on its helicase activity, thereby promoting transcription of its target genes. Our results uncover the multifaceted role of DDX21 in multiple steps of ribosome biogenesis, and provide evidence implicating a mammalian RNA helicase in RNA modification and Pol II elongation control.

m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells

N6-methyl-adenosine (m(6)A) is the most abundant modification on messenger RNAs and is linked to human diseases, but its functions in mammalian development are poorly understood. Here we reveal the evolutionary conservation and function of m(6)A by mapping the m(6)A methylome in mouse and human embryonic stem cells. Thousands of messenger and long noncoding RNAs show conserved m(6)A modification, including transcripts encoding core pluripotency transcription factors. m(6)A is enriched over 3' untranslated regions at defined sequence motifs and marks unstable transcripts, including transcripts turned over upon differentiation. Genetic inactivation or depletion of mouse and human Mettl3, one of the m(6)A methylases, led to m(6)A erasure on select target genes, prolonged Nanog expression upon differentiation, and impaired ESC exit from self-renewal toward differentiation into several lineages in vitro and in vivo. Thus, m(6)A is a mark of transcriptome flexibility required for stem cells to differentiate to specific lineages.

RNA structural analysis by evolving SHAPE chemistry

RNA is central to the flow of biological information. From transcription to splicing, RNA localization, translation, and decay, RNA is intimately involved in regulating every step of the gene expression program, and is thus essential for health and understanding disease. RNA has the unique ability to base-pair with itself and other nucleic acids to form complex structures. Hence the information content in RNA is not simply its linear sequence of bases, but is also encoded in complex folding of RNA molecules. A general chemical functionality that all RNAs have is a 2'-hydroxyl group in the ribose ring, and the reactivity of the 2'-hydroxyl in RNA is gated by local nucleotide flexibility. In other words, the 2'-hydroxyl is reactive at single-stranded and conformationally flexible positions but is unreactive at nucleotides constrained by base-pairing. Recent efforts have been focused on developing reagents that modify RNA as a function of RNA 2' hydroxyl group reactivity. Such RNA structure probing techniques can be read out by primer extension in experiments termed RNA SHAPE (selective 2'- hydroxyl acylation and primer extension). Herein, we describe the efforts devoted to the design and utilization of SHAPE probes for characterizing RNA structure. We also describe current technological advances that are being applied to utilize SHAPE chemistry with deep sequencing to probe many RNAs in parallel. The merging of chemistry with genomics is sure to open the door to genome-wide exploration of RNA structure and function.

Landscape and variation of RNA secondary structure across the human transcriptome

In parallel to the genetic code for protein synthesis, a second layer of information is embedded in all RNA transcripts in the form of RNA structure. RNA structure influences practically every step in the gene expression program¹. Yet the nature of most RNA structures or effects of sequence variation on structure are not known. Here we report the initial landscape and variation of RNA secondary structures (RSS) in a human family Trio, providing a comprehensive RSS map of human coding and noncoding RNAs. We identify unique RSS signatures that demarcate open reading frames, splicing junctions, and define authentic microRNA binding sites. Comparison of native deproteinized RNA isolated from cells versus refolded purified RNA suggests that the majority of the RSS information is encoded within RNA sequence. Over 1900 transcribed single nucleotide variants (~15% of all transcribed SNVs) alter local RNA structure. We discover simple sequence and spacing rules that determine the ability of point mutations to impact RSS. Selective depletion of RiboSNitches versus structurally synonymous variants at precise locations suggests selection for specific RNA shapes at thousands of sites, including 3’UTRs, binding sites of miRNAs and RNA binding proteins genome-wide. These results highlight the potentially broad contribution of RNA structure and its variation to gene regulation.