Jiayu Wen

A high-resolution map of human evolutionary constraint using 29 mammals

The comparison of related genomes has emerged as a powerful lens for genome interpretation. Here we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and locate constrained elements covering ∼4.2% of the genome. We use evolutionary signatures and comparisons with experimental data sets to suggest candidate functions for ∼60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements and more than 1,000 primate- and human-accelerated elements. Overlap with disease-associated variants indicates that our findings will be relevant for studies of human biology, health and disease.

Diversity and dynamics of the Drosophila transcriptome

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

microRNA-101 is a potent inhibitor of autophagy

Autophagy is an evolutionarily conserved mechanism of cellular self‐digestion in which proteins and organelles are degraded through delivery to lysosomes. Defects in this process are implicated in numerous human diseases including cancer. To further elucidate regulatory mechanisms of autophagy, we performed a functional screen in search of microRNAs (miRNAs), which regulate the autophagic flux in breast cancer cells. In this study, we identified the tumour suppressive miRNA, miR‐101, as a potent inhibitor of basal, etoposide‐ and rapamycin‐induced autophagy. Through transcriptome profiling, we identified three novel miR‐101 targets, STMN1, RAB5A and ATG4D. siRNA‐mediated depletion of these genes phenocopied the effect of miR‐101 overexpression, demonstrating their importance in autophagy regulation. Importantly, overexpression of STMN1 could partially rescue cells from miR‐101‐mediated inhibition of autophagy, indicating a functional importance for this target. Finally, we show that miR‐101‐mediated inhibition of autophagy can sensitize breast cancer cells to 4‐hydroxytamoxifen (4‐OHT)‐mediated cell death. Collectively, these data establish a novel link between two highly important and rapidly growing research fields and present a new role for miR‐101 as a key regulator of autophagy.

Comparative analysis of the transcriptome across distant species.

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a ‘universal model’ based on a single set of organism-independent parameters.

MicroRNA-145 Targets YES and STAT1 in Colon Cancer Cells

Background MicroRNAs (miRNAs) have emerged as important gene regulators and are recognized as key players in tumorigenesis. miR-145 is reported to be down-regulated in several cancers, but knowledge of its targets in colon cancer remains limited. Methodology/Principal Findings To investigate the role of miR-145 in colon cancer, we have employed a microarray based approach to identify miR-145 targets. Based on seed site enrichment analyses and unbiased word analyses, we found a significant enrichment of miRNA binding sites in the 3′-untranslated regions (UTRs) of transcripts down-regulated upon miRNA overexpression. Gene Ontology analysis showed an overrepresentation of genes involved in cell death, cellular growth and proliferation, cell cycle, gene expression and cancer. A number of the identified miRNA targets have previously been implicated in cancer, including YES, FSCN1, ADAM17, BIRC2, VANGL1 as well as the transcription factor STAT1. Both YES and STAT1 were verified as direct miR-145 targets. Conclusions/Significance The study identifies and validates new cancer-relevant direct targets of miR-145 in colon cancer cells and hereby adds important mechanistic understanding of the tumor-suppressive functions of miR-145.

MicroRNA transfection and AGO-bound CLIP-seq data sets reveal distinct determinants of miRNA action.

Microarray expression analyses following miRNA transfection/inhibition and, more recently, Argonaute cross-linked immunoprecipitation (CLIP)-seq assays have been used to detect miRNA target sites. CLIP and expression approaches measure differing stages of miRNA functioning-initial binding of the miRNP complex and subsequent message repression. We use nonparametric predictive models to characterize a large number of known target and flanking features, utilizing miRNA transfection, HITS-CLIP, and PAR-CLIP data. In particular, we utilize the precise spatial information provided by CLIP-seq to analyze the predictive effect of target flanking features. We observe distinct target determinants between expression-based and CLIP-based data. Target flanking features such as flanking region conservation are an important AGO-binding determinant-we hypothesize that CLIP experiments have a preference for strongly bound miRNP-target interactions involving adjacent RNA-binding proteins that increase the strength of cross-linking. In contrast, seed-related features are major determinants in expression-based studies, but less so for CLIP-seq studies, and increased miRNA concentrations typical of transfection studies contribute to this difference. While there is a good overlap between miRNA targets detected by miRNA transfection and CLIP-seq, the detection of CLIP-seq targets is largely independent of the level of subsequent mRNA degradation. Also, models built using CLIP-seq data show strong predictive power between independent CLIP-seq data sets, but are not strongly predictive for expression change. Similarly, models built from expression data are not strongly predictive for CLIP-seq data sets, supporting the finding that the determinants of miRNA binding and mRNA degradation differ. Predictive models and results are available at http://servers.binf.ku.dk/antar/.

Diversity of miRNAs, siRNAs, and piRNAs across 25 Drosophila cell lines

We expanded the knowledge base for Drosophila cell line transcriptomes by deeply sequencing their small RNAs. In total, we analyzed more than 1 billion raw reads from 53 libraries across 25 cell lines. We verify reproducibility of biological replicate data sets, determine common and distinct aspects of miRNA expression across cell lines, and infer the global impact of miRNAs on cell line transcriptomes. We next characterize their commonalities and differences in endo-siRNA populations. Interestingly, most cell lines exhibit enhanced TE-siRNA production relative to tissues, suggesting this as a common aspect of cell immortalization. We also broadly extend annotations of cis-NAT-siRNA loci, identifying ones with common expression across diverse cells and tissues, as well as cell-restricted loci. Finally, we characterize small RNAs in a set of ovary-derived cell lines, including somatic cells (OSS and OSC) and a mixed germline/somatic cell population (fGS/OSS) that exhibits ping-pong piRNA signatures. Collectively, the ovary data reveal new genic piRNA loci, including unusual configurations of piRNA-generating regions. Together with the companion analysis of mRNAs described in a previous study, these small RNA data provide comprehensive information on the transcriptional landscape of diverse Drosophila cell lines. These data should encourage broader usage of fly cell lines, beyond the few that are presently in common usage.

Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data.

BackgroundInitiation and elongation of RNA polymerase II (RNAPII) transcription is regulated by both DNA sequence and chromatin signals. Recent breakthroughs make it possible to measure the chromatin state and activity of core promoters genome-wide, but dedicated computational strategies are needed to progress from descriptive annotation of data to quantitative, predictive models.ResultsHere, we describe a computational framework which with high accuracy can predict the locations of core promoters, the amount of recruited RNAPII at the promoter, the amount of elongating RNAPII in the gene body, the mRNA production originating from the promoter and finally also the stalling characteristics of RNAPII by considering both quantitative and spatial features of histone modifications around the transcription start site (TSS).As the model framework can also pinpoint the signals that are the most influential for prediction, it can be used to infer underlying regulatory biology. For example, we show that the H3K4 di- and tri- methylation signals are strongly predictive for promoter location while the acetylation marks H3K9 and H3K27 are highly important in estimating the promoter usage. All of these four marks are found to be necessary for recruitment of RNAPII but not sufficient for the elongation. We also show that the spatial distributions of histone marks are almost as predictive as the signal strength and that a set of histone marks immediately downstream of the TSS is highly predictive of RNAPII stalling.ConclusionsIn this study we introduce a general framework to accurately predict the level of RNAPII recruitment, elongation, stalling and mRNA expression from chromatin signals. The versatility of the method also makes it ideally suited to investigate other genomic data.

Loss of miR-10a activates lpo and collaborates with activated Wnt signaling in inducing intestinal neoplasia in female mice.

miRNAs are small regulatory RNAs that, due to their considerable potential to target a wide range of mRNAs, are implicated in essentially all biological process, including cancer. miR-10a is particularly interesting considering its conserved location in the Hox cluster of developmental regulators. A role for this microRNA has been described in developmental regulation as well as for various cancers. However, previous miR-10a studies are exclusively based on transient knockdowns of this miRNA and to extensively study miR-10a loss we have generated a miR-10a knock out mouse. Here we show that, in the Apcmin mouse model of intestinal neoplasia, female miR-10a deficient mice develop significantly more adenomas than miR-10+/+ and male controls. We further found that Lpo is extensively upregulated in the intestinal epithelium of mice deprived of miR-10a. Using in vitro assays, we demonstrate that the primary miR-10a target KLF4 can upregulate transcription of Lpo, whereas siRNA knockdown of KLF4 reduces LPO levels in HCT-116 cells. Furthermore, Klf4 is upregulated in the intestines of miR-10a knockout mice. Lpo has previously been shown to have the capacity to oxidize estrogens into potent depurinating mutagens, creating an instable genomic environment that can cause initiation of cancer. Therefore, we postulate that Lpo upregulation in the intestinal epithelium of miR-10a deficient mice together with the predominant abundance of estrogens in female animals mainly accounts for the sex-related cancer phenotype we observed. This suggests that miR-10a could be used as a potent diagnostic marker for discovering groups of women that are at high risk of developing colorectal carcinoma, which today is one of the leading causes of cancer-related deaths.

A non-conserved miRNA regulates lysosomal function and impacts on a human lysosomal storage disorder

Sulfatases are key enzymatic regulators of sulfate homeostasis with several biological functions including degradation of glycosaminoglycans (GAGs) and other macromolecules in lysosomes. In a severe lysosomal storage disorder, multiple sulfatase deficiency (MSD), global sulfatase activity is deficient due to mutations in the sulfatase-modifying factor 1 (SUMF1) gene, encoding the essential activator of all sulfatases. We identify a novel regulatory layer of sulfate metabolism mediated by a microRNA. miR-95 depletes SUMF1 protein levels and suppresses sulfatase activity, causing the disruption of proteoglycan catabolism and lysosomal function. This blocks autophagy-mediated degradation, causing cytoplasmic accumulation of autophagosomes and autophagic substrates. By targeting miR-95 in cells from MSD patients, we can effectively increase residual SUMF1 expression, allowing for reactivation of sulfatase activity and increased clearance of sulfated GAGs. The identification of this regulatory mechanism opens the opportunity for a unique therapeutic approach in MSD patients where the need for exogenous enzyme replacement is circumvented.