The RNase III enzyme DICER generates both microRNAs (miRNAs) and endogenous short interfering RNAs (endo-siRNAs). Both small RNA species silence gene expression post-transcriptionally in association with the ARGONAUTE (AGO) family of proteins. In mammals, there are four AGO proteins (AGO1-4), of which only AGO2 possesses endonucleolytic activity. siRNAs trigger endonucleolytic cleavage of target mRNAs, mediated by AGO2, whereas miRNAs cause translational repression and mRNA decay through association with any of the four AGO proteins. Dicer deletion in mouse oocytes leads to female infertility due to defects during meiosis I. Because mouse oocytes express both miRNAs and endo-siRNAs, this phenotype could be due to the absence of either class of small RNA, or both. However, we and others demonstrated that miRNA function is suppressed in mouse oocytes, which suggested that endo-siRNAs, not miRNAs, are essential for female meiosis. To determine if this was the case we generated mice that express a catalytically inactive knock-in allele of Ago2 (Ago2ADH) exclusively in oocytes and thereby disrupted the function of siRNAs. Oogenesis and hormonal response are normal in Ago2ADH oocytes, but meiotic maturation is impaired, with severe defects in spindle formation and chromosome alignment that lead to meiotic catastrophe. The transcriptome of these oocytes is widely perturbed and shows a highly significant correlation with the transcriptome of Dicer null and Ago2 null oocytes. Expression of the mouse transcript (MT), the most abundant transposable element in mouse oocytes, is increased. This study reveals that endo-siRNAs are essential during meiosis I in mouse females, demonstrating a role for endo-siRNAs in mammals.

TAR DNA-binding protein 43 (TDP-43) is normally a nuclear RNA-binding protein that exhibits a range of functions including regulation of alternative splicing, RNA trafficking, and RNA stability. However, in amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration with TDP-43 inclusions (FTLD-TDP), TDP-43 is abnormally phosphorylated, ubiquitinated, and cleaved, and is mislocalized to the cytoplasm where it forms distinctive aggregates. We previously developed a mouse model expressing human TDP-43 with a mutation in its nuclear localization signal (ΔNLS-hTDP-43) so that the protein preferentially localizes to the cytoplasm. These mice did not exhibit a significant number of cytoplasmic aggregates, but did display dramatic changes in gene expression as measured by microarray, suggesting that cytoplasmic TDP-43 may be associated with a toxic gain-of-function. Here, we analyze new RNA-sequencing data from the ΔNLS-hTDP-43 mouse model, together with published RNA-sequencing data obtained previously from TDP-43 antisense oligonucleotide (ASO) knockdown mice to investigate further the dysregulation of gene expression in the ΔNLS model. This analysis reveals that the transcriptomic effects of the overexpression of the ΔNLS-hTDP-43 transgene are likely due to a gain of cytoplasmic function. Moreover, cytoplasmic TDP-43 expression alters transcripts that regulate chromatin assembly, the nucleolus, lysosomal function, and histone 3' untranslated region (UTR) processing. These transcriptomic alterations correlate with observed histologic abnormalities in heterochromatin structure and nuclear size in transgenic mouse and human brains.

Small non-coding RNAs (sncRNAs) are highly abundant RNAs, typically <100 nucleotides long, that act as key regulators of diverse cellular processes. Although thousands of sncRNA genes are known to exist in the human genome, no single database provides searchable, unified annotation, and expression information for full sncRNA transcripts and mature RNA products derived from these larger RNAs. Here, we present the Database of small human noncoding RNAs (DASHR). DASHR contains the most comprehensive information to date on human sncRNA genes and mature sncRNA products. DASHR provides a simple user interface for researchers to view sequence and secondary structure, compare expression levels, and evidence of specific processing across all sncRNA genes and mature sncRNA products in various human tissues. DASHR annotation and expression data covers all major classes of sncRNAs including microRNAs (miRNAs), Piwi-interacting (piRNAs), small nuclear, nucleolar, cytoplasmic (sn-, sno-, scRNAs, respectively), transfer (tRNAs), and ribosomal RNAs (rRNAs). Currently, DASHR (v1.0) integrates 187 smRNA high-throughput sequencing (smRNA-seq) datasets with over 2.5 billion reads and annotation data from multiple public sources. DASHR contains annotations for ∼ 48,000 human sncRNA genes and mature sncRNA products, 82% of which are expressed in one or more of the curated tissues. DASHR is available at http://lisanwanglab.org/DASHR.

Enhancer elements are essential for tissue-specific gene regulation during mammalian development. Although these regulatory elements are often distant from their target genes, they affect gene expression by recruiting transcription factors to specific promoter regions. Because of this long-range action, the annotation of enhancer element-target promoter pairs remains elusive. Here, we developed a novel analysis methodology that takes advantage of Hi-C data to comprehensively identify these interactions throughout the human genome. To do this, we used a geometric distribution-based model to identify DNA-DNA interaction hotspots that contact gene promoters with high confidence. We observed that these promoter-interacting hotspots significantly overlap with known enhancer-associated histone modifications and DNase I hypersensitive sites. Thus, we defined thousands of candidate enhancer elements by incorporating these features, and found that they have a significant propensity to be bound by p300, an enhancer binding transcription factor. Furthermore, we revealed that their target genes are significantly bound by RNA Polymerase II and demonstrate tissue-specific expression. Finally, we uncovered that these elements are generally found within 1 Mb of their targets, and often regulate multiple genes. In total, our study presents a novel high-throughput workflow for confident, genome-wide discovery of enhancer-target promoter pairs, which will significantly improve our understanding of these regulatory interactions.

ATP-dependent chromatin remodeling is involved in all DNA transactions and is linked to numerous human diseases. We explored functions of chromatin remodelers during cellular aging. Deletion of ISW2, or mutations inactivating the Isw2 enzyme complex, extends yeast replicative lifespan. This extension by ISW2 deletion is epistatic to the longevity effect of calorie restriction (CR), and this mechanism is distinct from suppression of TOR signaling by CR. Transcriptome analysis indicates that isw2Δ partially mimics an upregulated stress response in CR cells. In particular, isw2Δ cells show an increased response to genotoxic stresses, and the DNA repair enzyme Rad51 is important for isw2Δ-mediated longevity. We show that lifespan is also extended in C. elegans by reducing levels of athp-2, a putative ortholog of Itc1/ACF1, a critical subunit of the enzyme complex. Our findings demonstrate that the ISWI class of ATP-dependent chromatin remodeling complexes plays a conserved role during aging and in CR.

MicroRNAs (miRNAs) are essential for regulation of gene expression. Though numerous miRNAs have been identified by high-throughput sequencing, few precursor miRNAs (pre-miRNAs) are experimentally validated. Here we report a strategy for constructing high-throughput sequencing libraries enriched for full-length pre-miRNAs. We find widespread and extensive uridylation of Argonaute (Ago)-bound pre-miRNAs, which is primarily catalyzed by two terminal uridylyltransferases: TUT7 and TUT4. Uridylation by TUT7/4 not only polishes pre-miRNA 3' ends, but also facilitates their degradation by the exosome, preventing clogging of Ago with defective species. We show that the exosome exploits distinct substrate preferences of DIS3 and RRP6, its two catalytic subunits, to distinguish productive from defective pre-miRNAs. Furthermore, we identify a positive feedback loop formed by the exosome and TUT7/4 in triggering uridylation and degradation of Ago-bound pre-miRNAs. Our study reveals a pre-miRNA surveillance system that comprises TUT7, TUT4, and the exosome in quality control of miRNA synthesis.

Two features of eukaryotic RNA molecules that regulate their post-transcriptional fates are RNA secondary structure and RNA-binding protein (RBP) interaction sites. However, a comprehensive global overview of the dynamic nature of these sequence features during erythropoiesis has never been obtained. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approach to reveal the global landscape of RNA secondary structure and RBP-RNA interaction sites and the dynamics of these features during this important developmental process. We identify dynamic patterns of RNA secondary structure and RBP binding throughout the process and determine a set of corresponding protein-bound sequence motifs along with their dynamic structural and RBP-binding contexts. Finally, using these dynamically bound sequences, we identify a number of RBPs that have known and putative key functions in post-transcriptional regulation during mammalian erythropoiesis. In total, this global analysis reveals new post-transcriptional regulators of mammalian blood cell development.

Next-generation RNA-sequencing is an incredibly powerful means of generating a snapshot of the transcriptomic state within a cell, tissue, or whole organism. As the questions addressed by RNA-sequencing (RNA-seq) become both more complex and greater in number, there is a need to simplify RNA-seq processing workflows, make them more efficient and interoperable, and capable of handling both large and small datasets. This is especially important for researchers who need to process hundreds to tens of thousands of RNA-seq datasets. To address these needs, we have developed a scalable, user-friendly, and easily deployable analysis suite called RMTA (Read Mapping, Transcript Assembly). RMTA can easily process thousands of RNA-seq datasets with features that include automated read quality analysis, filters for lowly expressed transcripts, and read counting for differential expression analysis. RMTA is containerized using Docker for easy deployment within any compute environment [cloud, local, or high-performance computing (HPC)] and is available as two apps in CyVerse's Discovery Environment, one for normal use and one specifically designed for introducing undergraduates and high school to RNA-seq analysis. For extremely large datasets (tens of thousands of FASTq files) we developed a high-throughput, scalable, and parallelized version of RMTA optimized for launching on the Open Science Grid (OSG) from within the Discovery Environment. OSG-RMTA allows users to utilize the Discovery Environment for data management, parallelization, and submitting jobs to OSG, and finally, employ the OSG for distributed, high throughput computing. Alternatively, OSG-RMTA can be run directly on the OSG through the command line. RMTA is designed to be useful for data scientists, of any skill level, interested in rapidly and reproducibly analyzing their large RNA-seq data sets.

Type II topoisomerases orchestrate proper DNA topology, and they are the targets of anti-cancer drugs that cause treatment-related leukemias with balanced translocations. Here, we develop a high-throughput sequencing technology to define TOP2 cleavage sites at single-base precision, and use the technology to characterize TOP2A cleavage genome-wide in the human K562 leukemia cell line. We find that TOP2A cleavage has functionally conserved local sequence preferences, occurs in cleavage cluster regions (CCRs), and is enriched in introns and lincRNA loci. TOP2A CCRs are biased toward the distal regions of gene bodies, and TOP2 poisons cause a proximal shift in their distribution. We find high TOP2A cleavage levels in genes involved in translocations in TOP2 poison-related leukemia. In addition, we find that a large proportion of genes involved in oncogenic translocations overall contain TOP2A CCRs. The TOP2A cleavage of coding and lincRNA genes is independently associated with both length and transcript abundance. Comparisons to ENCODE data reveal distinct TOP2A CCR clusters that overlap with marks of transcription, open chromatin, and enhancers. Our findings implicate TOP2A cleavage as a broad DNA damage mechanism in oncogenic translocations as well as a functional role of TOP2A cleavage in regulating transcription elongation and gene activation.

Aging is a major risk factor for many neurodegenerative disorders. A key feature of aging biology that may underlie these diseases is cellular senescence. Senescent cells accumulate in tissues with age, undergo widespread changes in gene expression, and typically demonstrate altered, pro-inflammatory profiles. Astrocyte senescence has been implicated in neurodegenerative disease, and to better understand senescence-associated changes in astrocytes, we investigated changes in their transcriptome using RNA sequencing. Senescence was induced in human fetal astrocytes by transient oxidative stress. Brain-expressed genes, including those involved in neuronal development and differentiation, were downregulated in senescent astrocytes. Remarkably, several genes indicative of astrocytic responses to injury were also downregulated, including glial fibrillary acidic protein and genes involved in the processing and presentation of antigens by major histocompatibility complex class II proteins, while pro-inflammatory genes were upregulated. Overall, our findings suggest that senescence-related changes in the function of astrocytes may impact the pathogenesis of age-related brain disorders.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

The N6-methyladenosine (m6A) modification is the most prevalent internal RNA modification in eukaryotes. The majority of m6A sites are found in the last exon and 3' UTRs. Here we show that the nuclear m6A reader YTHDC1 is essential for embryo viability and germline development in mouse. Specifically, YTHDC1 is required for spermatogonial development in males and for oocyte growth and maturation in females; Ythdc1-deficient oocytes are blocked at the primary follicle stage. Strikingly, loss of YTHDC1 leads to extensive alternative polyadenylation in oocytes, altering 3' UTR length. Furthermore, YTHDC1 deficiency causes massive alternative splicing defects in oocytes. The majority of splicing defects in mutant oocytes are rescued by introducing wild-type, but not m6A-binding-deficient, YTHDC1. YTHDC1 is associated with the pre-mRNA 3' end processing factors CPSF6, SRSF3, and SRSF7. Thus, YTHDC1 plays a critical role in processing of pre-mRNA transcripts in the oocyte nucleus and may have similar non-redundant roles throughout fetal development.

Most regulatory chromatin interactions are mediated by various transcription factors (TFs) and involve physically interacting elements such as enhancers, insulators or promoters. To map these elements and interactions at a fine scale, we developed HIPPIE2 that analyzes raw reads from high-throughput chromosome conformation (Hi-C) experiments to identify precise loci of DNA physically interacting regions (PIRs). Unlike standard genome binning approaches (e.g. 10-kb to 1-Mb bins), HIPPIE2 dynamically infers the physical locations of PIRs using the distribution of restriction sites to increase analysis precision and resolution. We applied HIPPIE2 to in situ Hi-C datasets across six human cell lines (GM12878, IMR90, K562, HMEC, HUVEC, NHEK) with matched ENCODE/Roadmap functional genomic data. HIPPIE2 detected 1042 738 distinct PIRs, with high resolution (average PIR length of 1006 bp) and high reproducibility (92.3% in GM12878). PIRs are enriched for epigenetic marks (H3K27ac, H3K4me1) and open chromatin, suggesting active regulatory roles. HIPPIE2 identified 2.8 million significant PIR-PIR interactions, 27.2% of which were enriched for TF binding sites. 50 608 interactions were enhancer-promoter interactions and were enriched for 33 TFs, including known DNA looping/long-range mediators. These findings demonstrate that the novel dynamic approach of HIPPIE2 (https://bitbucket.com/wanglab-upenn/HIPPIE2) enables the characterization of chromatin and regulatory interactions with high resolution and reproducibility.

Hematopoietic stem cell (HSC) transplantation using umbilical cord blood (UCB) is a potentially life-saving treatment for leukemia and bone marrow failure but is limited by the low number of HSCs in UCB. The loss of HSCs after ex vivo manipulation is also a major obstacle to gene editing for inherited blood disorders. HSCs require a low rate of translation to maintain their capacity for self-renewal, but hematopoietic cytokines used to expand HSCs stimulate protein synthesis and impair long-term self-renewal. We previously described cytokine-free conditions that maintain but do not expand human and mouse HSCs ex vivo. Here we performed a high throughput screen and identified translation inhibitors that allow ex vivo expansion of human HSCs while minimizing cytokine exposure. Transplantation assays show a ~5-fold expansion of long-term HSCs from UCB after one week of culture in low cytokine conditions. Single cell transcriptomic analysis demonstrates maintenance of HSCs expressing mediators of the unfolded protein stress response, further supporting the importance of regulated proteostasis in HSC maintenance and expansion. This expansion method maintains and expands human HSCs after CRISPR/Cas9 editing of the BCL11A+58 enhancer, overcoming a major obstacle to ex vivo gene correction for human hemoglobinopathies.

N6-methyladenosine (m6A) is a dynamic, reversible, covalently modified ribonucleotide that occurs predominantly toward 3' ends of eukaryotic mRNAs and is essential for their proper function and regulation. In Arabidopsis thaliana, many RNAs contain at least one m6A site, yet the transcriptome-wide function of m6A remains mostly unknown. Here, we show that many m6A-modified mRNAs in Arabidopsis have reduced abundance in the absence of this mark. The decrease in abundance is due to transcript destabilization caused by cleavage occurring 4 or 5 nt directly upstream of unmodified m6A sites. Importantly, we also find that, upon agriculturally relevant salt treatment, m6A is dynamically deposited on and stabilizes transcripts encoding proteins required for salt and osmotic stress response. Overall, our findings reveal that m6A generally acts as a stabilizing mark through inhibition of site-specific cleavage in plant transcriptomes, and this mechanism is required for proper regulation of the salt-stress-responsive transcriptome.

Type 2 diabetes mellitus (T2DM) is a complex disease characterized by the inability of the insulin-producing β cells in the endocrine pancreas to overcome insulin resistance in peripheral tissues. To determine if microRNAs are involved in the pathogenesis of human T2DM, we sequenced the small RNAs of human islets from diabetic and nondiabetic organ donors. We identified a cluster of microRNAs in an imprinted locus on human chromosome 14q32 that is highly and specifically expressed in human β cells and dramatically downregulated in islets from T2DM organ donors. The downregulation of this locus strongly correlates with hypermethylation of its promoter. Using HITS-CLIP for the essential RISC-component Argonaute, we identified disease-relevant targets of the chromosome 14q32 microRNAs, such as IAPP and TP53INP1, that cause increased β cell apoptosis upon overexpression in human islets. Our results support a role for microRNAs and their epigenetic control by DNA methylation in the pathogenesis of T2DM.

Persistent infection of the mouse central nervous system (CNS) with mouse hepatitis virus (MHV) induces a demyelinating disease pathologically similar to multiple sclerosis and is therefore used as a model system. There is little information regarding the host factors that correlate with and contribute to MHV-induced demyelination. Here, we detail the genes and pathways associated with MHV-induced demyelinating disease in the spinal cord. High-throughput sequencing of the host transcriptome revealed that demyelination is accompanied by numerous transcriptional changes indicative of immune infiltration as well as changes in the cytokine milieu and lipid metabolism. We found evidence that a Th1-biased cytokine/chemokine response and eicosanoid-derived inflammation accompany persistent MHV infection and that antigen presentation is ongoing. Interestingly, increased expression of genes involved in lipid transport, processing, and catabolism, including some with known roles in neurodegenerative diseases, coincided with demyelination. Lastly, expression of several genes involved in osteoclast or bone-resident macrophage function, most notably TREM2 and DAP12, was upregulated in persistently infected mouse spinal cord. This study highlights the complexity of the host antiviral response, which accompany MHV-induced demyelination, and further supports previous findings that MHV-induced demyelination is immune-mediated. Interestingly, these data suggest a parallel between bone reabsorption by osteoclasts and myelin debris clearance by microglia in the bone and the CNS, respectively. To our knowledge, this is the first report of using an RNA-seq approach to study the host CNS response to persistent viral infection.

Senescence is a stable proliferation arrest, associated with an altered secretory pathway, thought to promote tumor suppression and tissue aging. While chromatin regulation and lamin B1 down-regulation have been implicated as senescence effectors, functional interactions between them are poorly understood. We compared genome-wide Lys4 trimethylation on histone H3 (H3K4me3) and H3K27me3 distributions between proliferating and senescent human cells and found dramatic differences in senescence, including large-scale domains of H3K4me3- and H3K27me3-enriched "mesas" and H3K27me3-depleted "canyons." Mesas form at lamin B1-associated domains (LADs) in replicative senescence and oncogene-induced senescence and overlap DNA hypomethylation regions in cancer, suggesting that pre-malignant senescent chromatin changes foreshadow epigenetic cancer changes. Hutchinson-Gilford progeria syndrome fibroblasts (mutant lamin A) also show evidence of H3K4me3 mesas, suggesting a link between premature chromatin changes and accelerated cell senescence. Canyons mostly form between LADs and are enriched in genes and enhancers. H3K27me3 loss is correlated with up-regulation of key senescence genes, indicating a link between global chromatin changes and local gene expression regulation. Lamin B1 reduction in proliferating cells triggers senescence and formation of mesas and canyons. Our data illustrate profound chromatin reorganization during senescence and suggest that lamin B1 down-regulation in senescence is a key trigger of global and local chromatin changes that impact gene expression, aging, and cancer.

Validation of physiologic miRNA targets has been met with significant challenges. We employed HITS-CLIP to identify which miRNAs participate in liver regeneration, and to identify their target mRNAs.

Wild populations of the model organism Drosophila melanogaster experience highly heterogeneous environments over broad geographical ranges as well as over seasonal and annual timescales. Diapause is a primary adaptation to environmental heterogeneity, and in D. melanogaster the propensity to enter diapause varies predictably with latitude and season. Here we performed global transcriptomic profiling of naturally occurring variation in diapause expression elicited by short day photoperiod and moderately low temperature in two tissue types associated with neuroendocrine and endocrine signaling, heads, and ovaries. We show that diapause in D. melanogaster is an actively regulated phenotype at the transcriptional level, suggesting that diapause is not a simple physiological or reproductive quiescence. Differentially expressed genes and pathways are highly distinct in heads and ovaries, demonstrating that the diapause response is not uniform throughout the soma and suggesting that it may be comprised of functional modules associated with specific tissues. Genes downregulated in heads of diapausing flies are significantly enriched for clinally varying single nucleotide polymorphism (SNPs) and seasonally oscillating SNPs, consistent with the hypothesis that diapause is a driving phenotype of climatic adaptation. We also show that chromosome location-based coregulation of gene expression is present in the transcriptional regulation of diapause. Taken together, these results demonstrate that diapause is a complex phenotype actively regulated in multiple tissues, and support the hypothesis that natural variation in diapause propensity underlies adaptation to spatially and temporally varying selective pressures.