Search

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.

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.

RNA molecules can fold into complex structures and interact with trans-acting factors to control their biology. Recent methods have been focused on developing novel tools to measure RNA structure transcriptome-wide, but their utility to study and predict RNA-protein interactions or RNA processing has been limited thus far. Here, we extend these studies with the first transcriptome-wide mapping method for cataloging RNA solvent accessibility, icLASER. By combining solvent accessibility (icLASER) with RNA flexibility (icSHAPE) data, we efficiently predict RNA-protein interactions transcriptome-wide and catalog RNA polyadenylation sites by RNA structure alone. These studies showcase the power of designing novel chemical approaches to studying RNA biology. Further, our study exemplifies merging complementary methods to measure RNA structure inside cells and its utility for predicting transcriptome-wide interactions that are critical for control of and regulation by RNA structure. We envision such approaches can be applied to studying different cell types or cells under varying conditions, using RNA structure and footprinting to characterize cellular interactions and processing involving RNA.

For decades the field of radiation oncology has sought to improve the therapeutic ratio through innovations in physics, chemistry, and biology. To date, technological advancements in image guided beam delivery techniques have provided clinicians with their best options for improving this critical tool in cancer care. Medical physics has focused on the preferential targeting of tumors while minimizing the collateral dose to the surrounding normal tissues, yielding only incremental progress. However, recent developments involving ultra-high dose rate irradiation termed FLASH radiotherapy (FLASH-RT), that were initiated nearly 50 years ago, have stimulated a renaissance in the field of radiotherapy, long awaiting a breakthrough modality able to enhance therapeutic responses and limit normal tissue injury. Compared to conventional dose rates used clinically (0.1-0.2 Gy/s), FLASH can implement dose rates of electrons or X-rays in excess of 100 Gy/s. The implications of this ultra-fast delivery of dose are significant and need to be re-evaluated to appreciate the fundamental aspects underlying this seemingly unique radiobiology. The capability of FLASH to significantly spare normal tissue complications in multiple animal models, when compared to conventional rates of dose-delivery, while maintaining persistent growth inhibition of select tumor models has generated considerable excitement, as well as skepticism. Based on fundamental principles of radiation physics, radio-chemistry, and tumor vs. normal cell redox metabolism, this article presents a series of testable, biologically relevant hypotheses, which may help rationalize the differential effects of FLASH irradiation observed between normal tissue and tumors.

Chronic graft-versus-host disease (cGVHD) is a leading cause of morbidity and mortality following allotransplant. Activated donor effector T cells can differentiate into pathogenic T helper (Th)-17 cells and germinal center (GC)-promoting T follicular helper (Tfh) cells, resulting in cGVHD. Phosphoinositide-3-kinase-δ (PI3Kδ), a lipid kinase, is critical for activated T cell survival, proliferation, differentiation, and metabolism. We demonstrate PI3Kδ activity in donor T cells that become Tfh cells is required for cGVHD in a nonsclerodermatous multiorgan system disease model that includes bronchiolitis obliterans (BO), dependent upon GC B cells, Tfhs, and counterbalanced by T follicular regulatory cells, each requiring PI3Kδ signaling for function and survival. Although B cells rely on PI3Kδ pathway signaling and GC formation is disrupted resulting in a substantial decrease in Ig production, PI3Kδ kinase-dead mutant donor bone marrow-derived GC B cells still supported BO cGVHD generation. A PI3Kδ-specific inhibitor, compound GS-649443, that has superior potency to idelalisib while maintaining selectivity, reduced cGVHD in mice with active disease. In a Th1-dependent and Th17-associated scleroderma model, GS-649443 effectively treated mice with active cGVHD. These data provide a foundation for clinical trials of US Food and Drug Administration (FDA)-approved PI3Kδ inhibitors for cGVHD therapy in patients.

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.

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.

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 DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-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.

RNA has the intrinsic property to base pair, forming complex structures fundamental to its diverse functions. Here, we develop PARIS, a method based on reversible psoralen crosslinking for global mapping of RNA duplexes with near base-pair resolution in living cells. PARIS analysis in three human and mouse cell types reveals frequent long-range structures, higher-order architectures, and RNA-RNA interactions in trans across the transcriptome. PARIS determines base-pairing interactions on an individual-molecule level, revealing pervasive alternative conformations. We used PARIS-determined helices to guide phylogenetic analysis of RNA structures and discovered conserved long-range and alternative structures. XIST, a long noncoding RNA (lncRNA) essential for X chromosome inactivation, folds into evolutionarily conserved RNA structural domains that span many kilobases. XIST A-repeat forms complex inter-repeat duplexes that nucleate higher-order assembly of the key epigenetic silencing protein SPEN. PARIS is a generally applicable and versatile method that provides novel insights into the RNA structurome and interactome. VIDEO ABSTRACT.

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 DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-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.

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.

While much is known about genes that promote aging, little is known about genes that protect against or prevent aging, particularly in human skin. The main objective of this study was to perform an unbiased, whole transcriptome search for genes that associate with intrinsic skin youthfulness. To accomplish this, healthy women (n = 122) of European descent, ages 18-89 years with Fitzpatrick skin type I/II were examined for facial skin aging parameters and clinical covariates, including smoking and ultraviolet exposure. Skin youthfulness was defined as the top 10% of individuals whose assessed skin aging features were most discrepant with their chronological ages. Skin biopsies from sun-protected inner arm were subjected to 3'-end sequencing for expression quantification, with results verified by quantitative reverse transcriptase-polymerase chain reaction. Unbiased clustering revealed gene expression signatures characteristic of older women with skin youthfulness (n = 12) compared to older women without skin youthfulness (n = 33), after accounting for gene expression changes associated with chronological age alone. Gene set analysis was performed using Genomica open-access software. This study identified a novel set of candidate skin youthfulness genes demonstrating differences between SY and non-SY group, including pleckstrin homology like domain family A member 1 (PHLDA1) (p = 2.4x10-5), a follicle stem cell marker, and hyaluronan synthase 2-anti-sense 1 (HAS2-AS1) (p = 0.00105), a non-coding RNA that is part of the hyaluronan synthesis pathway. We show that immunologic gene sets are the most significantly altered in skin youthfulness (with the most significant gene set p = 2.4x10-5), suggesting the immune system plays an important role in skin youthfulness, a finding that has not previously been recognized. These results are a valuable resource from which multiple future studies may be undertaken to better understand the mechanisms that promote skin youthfulness in humans.

While much is known about genes that promote aging, little is known about genes that protect against or prevent aging, particularly in human skin. The main objective of this study was to perform an unbiased, whole transcriptome search for genes that associate with intrinsic skin youthfulness. To accomplish this, healthy women (n = 122) of European descent, ages 18-89 years with Fitzpatrick skin type I/II were examined for facial skin aging parameters and clinical covariates, including smoking and ultraviolet exposure. Skin youthfulness was defined as the top 10% of individuals whose assessed skin aging features were most discrepant with their chronological ages. Skin biopsies from sun-protected inner arm were subjected to 3'-end sequencing for expression quantification, with results verified by quantitative reverse transcriptase-polymerase chain reaction. Unbiased clustering revealed gene expression signatures characteristic of older women with skin youthfulness (n = 12) compared to older women without skin youthfulness (n = 33), after accounting for gene expression changes associated with chronological age alone. Gene set analysis was performed using Genomica open-access software. This study identified a novel set of candidate skin youthfulness genes demonstrating differences between SY and non-SY group, including pleckstrin homology like domain family A member 1 (PHLDA1) (p = 2.4x10-5), a follicle stem cell marker, and hyaluronan synthase 2-anti-sense 1 (HAS2-AS1) (p = 0.00105), a non-coding RNA that is part of the hyaluronan synthesis pathway. We show that immunologic gene sets are the most significantly altered in skin youthfulness (with the most significant gene set p = 2.4x10-5), suggesting the immune system plays an important role in skin youthfulness, a finding that has not previously been recognized. These results are a valuable resource from which multiple future studies may be undertaken to better understand the mechanisms that promote skin youthfulness in humans.

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.

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.

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. However, 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 (RSSs) in a human family trio (mother, father and their child). This provides a comprehensive RSS map of human coding and non-coding RNAs. We identify unique RSS signatures that demarcate open reading frames and 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 1,900 transcribed single nucleotide variants (approximately 15% of all transcribed single nucleotide variants) alter local RNA structure. We discover simple sequence and spacing rules that determine the ability of point mutations to impact RSSs. Selective depletion of 'riboSNitches' versus structurally synonymous variants at precise locations suggests selection for specific RNA shapes at thousands of sites, including 3' untranslated regions, binding sites of microRNAs and RNA-binding proteins genome-wide. These results highlight the potentially broad contribution of RNA structure and its variation to gene regulation.

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. However, 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 (RSSs) in a human family trio (mother, father and their child). This provides a comprehensive RSS map of human coding and non-coding RNAs. We identify unique RSS signatures that demarcate open reading frames and 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 1,900 transcribed single nucleotide variants (approximately 15% of all transcribed single nucleotide variants) alter local RNA structure. We discover simple sequence and spacing rules that determine the ability of point mutations to impact RSSs. Selective depletion of 'riboSNitches' versus structurally synonymous variants at precise locations suggests selection for specific RNA shapes at thousands of sites, including 3' untranslated regions, binding sites of microRNAs and RNA-binding proteins genome-wide. These results highlight the potentially broad contribution of RNA structure and its variation to gene regulation.