Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP.

RNA transcripts are subject to posttranscriptional gene regulation involving hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) expressed in a cell-type dependent fashion. We developed a cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs. The crosslinked sites are revealed by thymidine to cytidine transitions in the cDNAs prepared from immunopurified RNPs of 4-thiouridine-treated cells. We determined the binding sites and regulatory consequences for several intensely studied RBPs and miRNPs, including PUM2, QKI, IGF2BP1-3, AGO/EIF2C1-4 and TNRC6A-C. Our study revealed that these factors bind thousands of sites containing defined sequence motifs and have distinct preferences for exonic versus intronic or coding versus untranslated transcript regions. The precise mapping of binding sites across the transcriptome will be critical to the interpretation of the rapidly emerging data on genetic variation between individuals and how these variations contribute to complex genetic diseases.

A role for Fkbp6 and the chaperone machinery in piRNA amplification and transposon silencing.

Epigenetic silencing of transposons by Piwi-interacting RNAs (piRNAs) constitutes an RNA-based genome defense mechanism. Piwi endonuclease action amplifies the piRNA pool by generating new piRNAs from target transcripts by a poorly understood mechanism. Here, we identified mouse Fkbp6 as a factor in this biogenesis pathway delivering piRNAs to the Piwi protein Miwi2. Mice lacking Fkbp6 derepress LINE1 (L1) retrotransposon and display reduced DNA methylation due to deficient nuclear accumulation of Miwi2. Like other cochaperones, Fkbp6 associates with the molecular chaperone Hsp90 via its tetratricopeptide repeat (TPR) domain. Inhibition of the ATP-dependent Hsp90 activity in an insect cell culture model results in the accumulation of short antisense RNAs in Piwi complexes. We identify these to be byproducts of piRNA amplification that accumulate only in nuage-localized Piwi proteins. We propose that the chaperone machinery normally ejects these inhibitory RNAs, allowing turnover of Piwi complexes for their continued participation in piRNA amplification.

Miwi catalysis is required for piRNA amplification-independent LINE1 transposon silencing.

Repetitive-element-derived Piwi-interacting RNAs (piRNAs) act together with Piwi proteins Mili (also known as Piwil2) and Miwi2 (also known as Piwil4) in a genome defence mechanism that initiates transposon silencing via DNA methylation in the mouse male embryonic germ line. This silencing depends on the participation of the Piwi proteins in a slicer-dependent piRNA amplification pathway and is essential for male fertility. A third Piwi family member, Miwi (also known as Piwil1), is expressed in specific postnatal germ cells and associates with a unique set of piRNAs of unknown function. Here we show that Miwi is a small RNA-guided RNase (slicer) that requires extensive complementarity for target cleavage in vitro. Disruption of its catalytic activity in mice by a single point mutation causes male infertility, and mutant germ cells show increased accumulation of LINE1 retrotransposon transcripts. We provide evidence for Miwi slicer activity directly cleaving transposon messenger RNAs, offering an explanation for the continued maintenance of repeat-derived piRNAs long after transposon silencing is established in germline stem cells. Furthermore, our study supports a slicer-dependent silencing mechanism that functions without piRNA amplification. Thus, Piwi proteins seem to act in a two-pronged mammalian transposon silencing strategy: one promotes transcriptional repression in the embryo, the other reinforces silencing at the post-transcriptional level after birth.

Cell type-selective disease-association of genes under high regulatory load.

We previously showed that disease-linked metabolic genes are often under combinatorial regulation. Using the genome-wide ChIP-Seq binding profiles for 93 transcription factors in nine different cell lines, we show that genes under high regulatory load are significantly enriched for disease-association across cell types. We find that transcription factor load correlates with the enhancer load of the genes and thereby allows the identification of genes under high regulatory load by epigenomic mapping of active enhancers. Identification of the high enhancer load genes across 139 samples from 96 different cell and tissue types reveals a consistent enrichment for disease-associated genes in a cell type-selective manner. The underlying genes are not limited to super-enhancer genes and show several types of disease-association evidence beyond genetic variation (such as biomarkers). Interestingly, the high regulatory load genes are involved in more KEGG pathways than expected by chance, exhibit increased betweenness centrality in the interaction network of liver disease genes, and carry longer 3' UTRs with more microRNA (miRNA) binding sites than genes on average, suggesting a role as hubs integrating signals within regulatory networks. In summary, epigenetic mapping of active enhancers presents a promising and unbiased approach for identification of novel disease genes in a cell type-selective manner.

Tead2 expression levels control the subcellular distribution of Yap and Taz, zyxin expression and epithelial-mesenchymal transition.

The cellular changes during an epithelial-mesenchymal transition (EMT) largely rely on global changes in gene expression orchestrated by transcription factors. Tead transcription factors and their transcriptional co-activators Yap and Taz have been previously implicated in promoting an EMT; however, their direct transcriptional target genes and their functional role during EMT have remained elusive. We have uncovered a previously unanticipated role of the transcription factor Tead2 during EMT. During EMT in mammary gland epithelial cells and breast cancer cells, levels of Tead2 increase in the nucleus of cells, thereby directing a predominant nuclear localization of its co-factors Yap and Taz via the formation of Tead2-Yap-Taz complexes. Genome-wide chromatin immunoprecipitation and next generation sequencing in combination with gene expression profiling revealed the transcriptional targets of Tead2 during EMT. Among these, zyxin contributes to the migratory and invasive phenotype evoked by Tead2. The results demonstrate that Tead transcription factors are crucial regulators of the cellular distribution of Yap and Taz, and together they control the expression of genes critical for EMT and metastasis.

Integrated analysis of transcript-level regulation of metabolism reveals disease-relevant nodes of the human metabolic network.

Metabolic diseases and comorbidities represent an ever-growing epidemic where multiple cell types impact tissue homeostasis. Here, the link between the metabolic and gene regulatory networks was studied through experimental and computational analysis. Integrating gene regulation data with a human metabolic network prompted the establishment of an open-sourced web portal, IDARE (Integrated Data Nodes of Regulation), for visualizing various gene-related data in context of metabolic pathways. Motivated by increasing availability of deep sequencing studies, we obtained ChIP-seq data from widely studied human umbilical vein endothelial cells. Interestingly, we found that association of metabolic genes with multiple transcription factors (TFs) enriched disease-associated genes. To demonstrate further extensions enabled by examining these networks together, constraint-based modeling was applied to data from human preadipocyte differentiation. In parallel, data on gene expression, genome-wide ChIP-seq profiles for peroxisome proliferator-activated receptor (PPAR) γ, CCAAT/enhancer binding protein (CEBP) α, liver X receptor (LXR) and H3K4me3 and microRNA target identification for miR-27a, miR-29a and miR-222 were collected. Disease-relevant key nodes, including mitochondrial glycerol-3-phosphate acyltransferase (GPAM), were exposed from metabolic pathways predicted to change activity by focusing on association with multiple regulators. In both cell types, our analysis reveals the convergence of microRNAs and TFs within the branched chain amino acid (BCAA) metabolic pathway, possibly providing an explanation for its downregulation in obese and diabetic conditions.

An epigenetic profile of early T-cell development from multipotent progenitors to committed T-cell descendants.

Cellular differentiation of the T-cell branch of the immune system begins with the HSC, which undergoes a series of stages characterized by progressive restriction in multipotency and acquisition of specific lineage identity At the molecular level, the restriction of cell potential, commitment, and differentiation to a specific lineage is achieved through the coordinated control of gene expression and epigenetic mechanisms. Here, we analyzed and compared the gene expression profiles and the genome-wide histone modification marks H3K4me3 (H3 lysine 4 trimethylation) and H3K27me3 (H3 lysine 27 trimethylation) in (i) in vitro propagated HSCs, (ii) in vitro generated and propagated pro-T cells derived from these stem cells, and (iii) double-positive thymocytes derived from these pro-T cells after injection into Rag-deficient mice. The combined analyses of the different datasets in this unique experimental system highlighted the importance of both transcriptional and epigenetic repression in shaping the early phases of T-cell development.

A Telemedicine Center Reduces the Comprehensive Carbon Footprint in Primary Care: A Monocenter, Retrospective Study.

INTRODUCTION: Telemedicine reduces greenhouse gas emissions (CO2eq); however, results of studies vary extremely in dependence of the setting. This is the first study to focus on effects of telemedicine on CO2 imprint of primary care. METHODS: We conducted a comprehensive retrospective study to analyze total CO2eq emissions of kilometers (km) saved by telemedical consultations. We categorized prevented and provoked patient journeys, including pharmacy visits. We calculated CO2eq emission savings through primary care telemedical consultations in comparison to those that would have occurred without telemedicine. We used the comprehensive footprint approach, including all telemedical cases and the CO2eq emissions by the telemedicine center infrastructure. In order to determine the net ratio of CO2eq emissions avoided by the telemedical center, we calculated the emissions associated with the provision of telemedical consultations (including also the total consumption of physicians' workstations) and subtracted them from the total of avoided CO2eq emissions. Furthermore, we also considered patient cases in our calculation that needed to have an in-person visit after the telemedical consultation. We calculated the savings taking into account the source of the consumed energy (renewable or not). RESULTS: 433 890 telemedical consultations overall helped save 1 800 391 km in travel. On average, 1 telemedical consultation saved 4.15 km of individual transport and consumed 0.15 kWh. We detected savings in almost every cluster of patients. After subtracting the CO2eq emissions caused by the telemedical center, the data reveal savings of 247.1 net tons of CO2eq emissions in total and of 0.57 kg CO2eq per telemedical consultation. The comprehensive footprint approach thus indicated a reduced footprint due to telemedicine in primary care. DISCUSSION: Integrating a telemedical center into the health care system reduces the CO2 footprint of primary care medicine; this is true even in a densely populated country with little use of cars like Switzerland. The insight of this study complements previous studies that focused on narrower aspects of telemedical consultations.

miRNA in situ hybridization in formaldehyde and EDC-fixed tissues.

MicroRNAs are small regulatory RNAs with many biological functions and disease associations. We showed that in situ hybridization (ISH) using conventional formaldehyde fixation results in substantial microRNA loss from mouse tissue sections, which can be prevented by fixation with 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide that irreversibly immobilizes the microRNA at its 5' phosphate. We determined optimal hybridization parameters for 130 locked nucleic acid probes by recording nucleic acid melting temperature during ISH.

Exploiting natural variation to uncover candidate genes that control element accumulation in Arabidopsis thaliana.

The plant ionome varies both inter- and intraspecifically despite the highly conserved roles for particular elements across the plant kingdom. Element storage requires transport across the plasma membrane and commonly deposition within the central vacuole. Therefore, tonoplast transport characteristics can be highly influential in controlling the plant ionome. As a result, individual cell types of the same plant, each with unique transcriptomes and vacuolar proteomes, can display very different elemental profiles. Here we address the use of natural variation in Arabidopsis thaliana for identifying genes involved in elemental accumulation. We present a conceptual framework, exploiting publicly available leaf ionomic and transcriptomic data across 31 Arabidopsis accessions, that promises to accelerate conventional forward genetics approaches for candidate gene discovery. Utilizing this framework, we identify numerous genes with documented roles in accumulation of calcium, magnesium and zinc and implicate additional candidate genes. Where appropriate, we discuss their role in cell-specific elemental accumulation. Currently, this framework could represent an alternate approach for identifying genes suitable for element biofortification of plants. Integration of additional cell-specific and whole-plant 'omics' datasets across Arabidopsis accessions under diverse environmental conditions should enable this concept to be developed into a scalable and robust tool for linking genotype and phenotype.

Conserved generation of short products at piRNA loci.

BACKGROUND: The piRNA pathway operates in animal germ lines to ensure genome integrity through retrotransposon silencing. The Piwi protein-associated small RNAs (piRNAs) guide Piwi proteins to retrotransposon transcripts, which are degraded and thereby post-transcriptionally silenced through a ping-pong amplification process. Cleavage of the retrotransposon transcript defines at the same time the 5' end of a secondary piRNA that will in turn guide a Piwi protein to a primary piRNA precursor, thereby amplifying primary piRNAs. Although several studies provided evidence that this mechanism is conserved among metazoa, how the process is initiated and what enzymatic activities are responsible for generating the primary and secondary piRNAs are not entirely clear. RESULTS: Here we analyzed small RNAs from three mammalian species, seeking to gain further insight into the mechanisms responsible for the piRNA amplification loop. We found that in all these species piRNA-directed targeting is accompanied by the generation of short sequences that have a very precisely defined length, 19 nucleotides, and a specific spatial relationship with the guide piRNAs. CONCLUSIONS: This suggests that the processing of the 5' product of piRNA-guided cleavage occurs while the piRNA target is engaged by the Piwi protein. Although they are not stabilized through methylation of their 3' ends, the 19-mers are abundant not only in testes lysates but also in immunoprecipitates of Miwi and Mili proteins. They will enable more accurate identification of piRNA loci in deep sequencing data sets.

MirZ: an integrated microRNA expression atlas and target prediction resource.

MicroRNAs (miRNAs) are short RNAs that act as guides for the degradation and translational repression of protein-coding mRNAs. A large body of work showed that miRNAs are involved in the regulation of a broad range of biological functions, from development to cardiac and immune system function, to metabolism, to cancer. For most of the over 500 miRNAs that are encoded in the human genome the functions still remain to be uncovered. Identifying miRNAs whose expression changes between cell types or between normal and pathological conditions is an important step towards characterizing their function as is the prediction of mRNAs that could be targeted by these miRNAs. To provide the community the possibility of exploring interactively miRNA expression patterns and the candidate targets of miRNAs in an integrated environment, we developed the MirZ web server, which is accessible at www.mirz.unibas.ch. The server provides experimental and computational biologists with statistical analysis and data mining tools operating on up-to-date databases of sequencing-based miRNA expression profiles and of predicted miRNA target sites in species ranging from Caenorhabditis elegans to Homo sapiens.

ChIP-seq profiling of the active chromatin marker H3K4me3 and PPARγ, CEBPα and LXR target genes in human SGBS adipocytes.

Transcription factors (TFs) represent key factors to establish a cellular phenotype. It is known that several TFs could play a role in disease, yet less is known so far how their targets overlap. We focused here on identifying the most highly induced TFs and their putative targets during human adipogenesis. Applying chromatin immunoprecipitation coupled with deep sequencing (ChIP-Seq) in the human SGBS pre-adipocyte cell line, we identified genes with binding sites in their vicinity for the three TFs studied, PPARγ, CEBPα and LXR. Here we describe the experimental design and quality controls in detail for the deep sequencing data and related results published by Galhardo et al. in Nucleic Acids Research 2014 [1] associated with the data uploaded to NCBI Gene Expression Omnibus (GSE41578).

Transcriptomics profiling of human SGBS adipogenesis.

Obesity is an ever-growing epidemic where tissue homeostasis is influenced by the differentiation of adipocytes that function in lipid metabolism, endocrine and inflammatory processes. While this differentiation process has been well-characterized in mice, limited data is available from human cells. Applying microarray expression profiling in the human SGBS pre-adipocyte cell line, we identified genes with differential expression during differentiation in combination with constraint-based modeling of metabolic pathway activity. Here we describe the experimental design and quality controls in detail for the gene expression and related results published by Galhardo et al. in Nucleic Acids Research 2014 associated with the data uploaded to NCBI Gene Expression Omnibus (GSE41352).

MicroRNA activity is suppressed in mouse oocytes.

MicroRNAs (miRNAs) are small endogenous RNAs that typically imperfectly base pair with 3' untranslated regions (3'UTRs) and mediate translational repression and mRNA degradation. Dicer, which generates small RNAs in the miRNA and RNA interference (RNAi) pathways, is essential for meiotic maturation of mouse oocytes. We found that 3'UTRs of transcripts upregulated in Dicer1(-/-) oocytes are not enriched in miRNA binding sites, implicating a weak impact of miRNAs on the maternal transcriptome. Therefore, we tested the ability of endogenous miRNAs to mediate RNA-like cleavage or translational repression of reporter mRNAs. In contrast to somatic cells, endogenous miRNAs in oocytes poorly repressed translation of mRNA reporters, whereas their RNAi-like activity was much less affected. Reporter mRNA carrying let-7-binding sites failed to localize to P body-like structures in oocytes. Our data suggest that miRNA function is downregulated during oocyte development, an idea supported by normal meiotic maturation of oocytes lacking Dgcr8, which is required for the miRNA but not the RNAi pathway (Suh et al. [1], this issue of Current Biology). Suppressing miRNA function during oocyte growth is likely an early event in reprogramming gene expression during the transition of a differentiated oocyte into pluripotent blastomeres of the embryo.

PAR-CliP--a method to identify transcriptome-wide the binding sites of RNA binding proteins.

RNA transcripts are subjected to post-transcriptional gene regulation by interacting with hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) that are often expressed in a cell-type dependently. To understand how the interplay of these RNA-binding factors affects the regulation of individual transcripts, high resolution maps of in vivo protein-RNA interactions are necessary. A combination of genetic, biochemical and computational approaches are typically applied to identify RNA-RBP or RNA-RNP interactions. Microarray profiling of RNAs associated with immunopurified RBPs (RIP-Chip) defines targets at a transcriptome level, but its application is limited to the characterization of kinetically stable interactions and only in rare cases allows to identify the RBP recognition element (RRE) within the long target RNA. More direct RBP target site information is obtained by combining in vivo UV crosslinking with immunoprecipitation followed by the isolation of crosslinked RNA segments and cDNA sequencing (CLIP). CLIP was used to identify targets of a number of RBPs. However, CLIP is limited by the low efficiency of UV 254 nm RNA-protein crosslinking, and the location of the crosslink is not readily identifiable within the sequenced crosslinked fragments, making it difficult to separate UV-crosslinked target RNA segments from background non-crosslinked RNA fragments also present in the sample. We developed a powerful cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs that we term PAR-CliP (Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation) (see Fig. 1A for an outline of the method). The method relies on the incorporation of photoreactive ribonucleoside analogs, such as 4-thiouridine (4-SU) and 6-thioguanosine (6-SG) into nascent RNA transcripts by living cells. Irradiation of the cells by UV light of 365 nm induces efficient crosslinking of photoreactive nucleoside-labeled cellular RNAs to interacting RBPs. Immunoprecipitation of the RBP of interest is followed by isolation of the crosslinked and coimmunoprecipitated RNA. The isolated RNA is converted into a cDNA library and deep sequenced using Solexa technology. One characteristic feature of cDNA libraries prepared by PAR-CliP is that the precise position of crosslinking can be identified by mutations residing in the sequenced cDNA. When using 4-SU, crosslinked sequences thymidine to cytidine transition, whereas using 6-SG results in guanosine to adenosine mutations. The presence of the mutations in crosslinked sequences makes it possible to separate them from the background of sequences derived from abundant cellular RNAs. Application of the method to a number of diverse RNA binding proteins was reported in Hafner et al.

Computational analysis of small RNA cloning data.

Cloning and sequencing is the method of choice for small regulatory RNA identification. Using deep sequencing technologies one can now obtain up to a billion nucleotides--and tens of millions of small RNAs--from a single library. Careful computational analyses of such libraries enabled the discovery of miRNAs, rasiRNAs, piRNAs, and 21U RNAs. Given the large number of sequences that can be obtained from each individual sample, deep sequencing may soon become an alternative to oligonucleotide microarray technology for mRNA expression profiling. In this report we present the methods that we developed for the annotation and expression profiling of small RNAs obtained through large-scale sequencing. These include a fast algorithm for finding nearly perfect matches of small RNAs in sequence databases, a web-accessible software system for the annotation of small RNA libraries, and a Bayesian method for comparing small RNA expression across samples.

MicroRNAs control de novo DNA methylation through regulation of transcriptional repressors in mouse embryonic stem cells.

Loss of microRNA (miRNA) pathway components negatively affects differentiation of embryonic stem (ES) cells, but the underlying molecular mechanisms remain poorly defined. Here we characterize changes in mouse ES cells lacking Dicer (Dicer1). Transcriptome analysis of Dicer-/- cells indicates that the ES-specific miR-290 cluster has an important regulatory function in undifferentiated ES cells. Consistently, many of the defects in Dicer-deficient cells can be reversed by transfection with miR-290 family miRNAs. We demonstrate that Oct4 (also known as Pou5f1) silencing in differentiating Dicer-/- ES cells is accompanied by accumulation of repressive histone marks but not by DNA methylation, which prevents the stable repression of Oct4. The methylation defect correlates with downregulation of de novo DNA methyltransferases (Dnmts). The downregulation is mediated by Rbl2 and possibly other transcriptional repressors, potential direct targets of miR-290 cluster miRNAs. The defective DNA methylation can be rescued by ectopic expression of de novo Dnmts or by transfection of the miR-290 cluster miRNAs, indicating that de novo DNA methylation in ES cells is controlled by miRNAs.