Lasse Sinkkonen

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.

Regulation of the human cyclin C gene via multiple vitamin D3-responsive regions in its promoter

The candidate human tumor suppressor gene cyclin C is a primary target of the anti-proliferative hormone 1α,25-dihydroxyvitamin D3 [1α,25(OH)2D3], but binding sites for the 1α,25(OH)2D3 receptor (VDR), so-called 1α,25(OH)2D3 response elements (VDREs), have not yet been identified in the promoter of this gene. We screened various cancer cell lines by quantitative PCR and found that the 1α,25(OH)2D3 inducibility of cyclin C mRNA expression, in relationship with the 24-hydroxylase (CYP24) gene, was best in MCF-7 human breast cancer cells. To characterize the molecular mechanisms, we analyzed 8.4 kb of the cyclin C promoter by using chromatin immunoprecipitation assays (ChIP) with antibodies against acetylated histone 4, VDR and its partner receptor, retinoid X receptor (RXR). The histone 4 acetylation status of all 23 investigated regions of the cyclin C promoter did not change significantly in response to 1α,25(OH)2D3, but four independent promoter regions showed a consistent, 1α,25(OH)2D3-dependent association with VDR and RXR over a time period of 240 min. Combined in silico/in vitro screening identified in each of these promoter regions a VDRE and reporter gene assays confirmed their functionality. Moreover, re-ChIP assays monitored simultaneous association of VDR with RXR, coactivator, mediator and RNA polymerase II proteins on these regions. Since cyclin C protein is associated with those mediator complexes that display transcriptional repressive properties, this study contributes to the understanding of the downregulation of a number of secondary 1α,25(OH)2D3-responding genes.

Neurodegeneration by activation of the microglial complement-phagosome pathway

Systemic inflammatory reactions have been postulated to exacerbate neurodegenerative diseases via microglial activation. We now demonstrate in vivo that repeated systemic challenge of mice over four consecutive days with bacterial LPS maintained an elevated microglial inflammatory phenotype and induced loss of dopaminergic neurons in the substantia nigra. The same total cumulative LPS dose given within a single application did not induce neurodegeneration. Whole-genome transcriptome analysis of the brain demonstrated that repeated systemic LPS application induced an activation pattern involving the classical complement system and its associated phagosome pathway. Loss of dopaminergic neurons induced by repeated systemic LPS application was rescued in complement C3-deficient mice, confirming the involvement of the complement system in neurodegeneration. Our data demonstrate that a phagosomal inflammatory response of microglia is leading to complement-mediated loss of dopaminergic neurons.

Dicer is associated with ribosomal DNA chromatin in mammalian cells.

Background RNA silencing is a common term for pathways utilizing small RNAs as sequence-specific guides to repress gene expression. Components of the RNA silencing machinery are involved in different aspects of chromatin function in numerous organisms. However, association of RNA silencing with chromatin in mammalian cells remains unclear. Methodology/Principal Findings Immunostaining of mitotic chromosomes with antibodies visualizing either endogenous or ectopically expressed Dicer in mammalian cells revealed association of the protein with ribosomal DNA (rDNA) repeats. Chromatin immunoprecipitations and bisulfite sequencing experiments indicated that Dicer is associated with transcribed regions of both active and silenced genes in rDNA arrays of interphase chromosomes. Metabolic labeling of the mouse embryonic stem (ES) cells lacking Dicer did not reveal apparent defect in rRNA biogenesis though pre-rRNA synthesis in these cells was decreased, likely as a consequence of their slower growth caused by the loss of miRNAs. We analyzed in detail chromatin structure of rDNA but did not find any epigenetic changes at rDNA loci in Dicer−/− ES cells. Instead, we found that rDNA methylation is rather low in primary tissues, contrasting with rDNA methylation patterns in transformed cell lines. Conclusion/Significance We found that Dicer, a key component of RNA silencing pathways, can be detected in association with rDNA chromatin in mammalian cells. The role of this particular localization of Dicer is not readily apparent since the enzyme is associated with rDNA genes regardless of their transcriptional activity. However, localization of Dicer to the transcribed region suggests that transcription may contribute to the Dicer deposition at rDNA chromatin. We hypothesize that Dicer functions in maintaining integrity of rDNA arrays.

Gene-pair expression signatures reveal lineage control

The distinct cell types of multicellular organisms arise due to constraints imposed by gene regulatory networks on the collective change of gene expression across the genome, creating self-stabilizing expression states, or attractors. We compiled a resource of curated human expression data comprising 166 cell types and 2,602 transcription regulating genes and developed a data driven method built around the concept of expression reversal defined at the level of gene pairs, such as those participating in toggle switch circuits. This approach allows us to organize the cell types into their ontogenetic lineage-relationships and to reflect regulatory relationships among genes that explain their ability to function as determinants of cell fate. We show that this method identifies genes belonging to regulatory circuits that control neuronal fate, pluripotency and blood cell differentiation, thus offering a novel large-scale perspective on lineage specification.The distinct cell types of multicellular organisms arise owing to constraints imposed by gene regulatory networks on the collective change of gene expression across the genome, creating self-stabilizing expression states, or attractors. We curated human expression data comprising 166 cell types and 2,602 transcription-regulating genes and developed a data-driven method for identifying putative determinants of cell fate built around the concept of expression reversal of gene pairs, such as those participating in toggle-switch circuits. This approach allows us to organize the cell types into their ontogenic lineage relationships. Our method identifies genes in regulatory circuits that control neuronal fate, pluripotency and blood cell differentiation, and it may be useful for prioritizing candidate factors for direct conversion of cell fate.

Dataset integration identifies transcriptional regulation of microRNA genes by PPARγ in differentiating mouse 3T3-L1 adipocytes

Peroxisome proliferator-activated receptor γ (PPARγ) is a key transcription factor in mammalian adipogenesis. Genome-wide approaches have identified thousands of PPARγ binding sites in mouse adipocytes and PPARγ upregulates hundreds of protein-coding genes during adipogenesis. However, no microRNA (miRNA) genes have been identified as primary PPARγ-targets. By integration of four separate datasets of genome-wide PPARγ binding sites in 3T3-L1 adipocytes we identified 98 miRNA clusters with PPARγ binding within 50 kb from miRNA transcription start sites. Nineteen mature miRNAs were upregulated ≥2-fold during adipogenesis and for six of these miRNA loci the PPARγ binding sites were confirmed by at least three datasets. The upregulation of five miRNA genes miR-103-1 (host gene Pank3), miR-148b (Copz1), miR-182/96/183, miR-205 and miR-378 (Ppargc1b) followed that of Pparg. The PPARγ-dependence of four of these miRNA loci was demonstrated by PPARγ knock-down and the loci of miR-103-1 (Pank3), miR-205 and miR-378 (Ppargc1b) were also responsive to the PPARγ ligand rosiglitazone. Finally, chromatin immunoprecipitation analysis validated in silico predicted PPARγ binding sites at all three loci and H3K27 acetylation was analyzed to confirm the activity of these enhancers. In conclusion, we identified 22 putative PPARγ target miRNA genes, showed the PPARγ dependence of four of these genes and demonstrated three as direct PPARγ target genes in mouse adipogenesis.

Controlling the chromatin organization of vitamin D target genes by multiple vitamin D receptor binding sites

An essential prerequisite for the direct modulation of transcription by 1alpha,25-dihydroxy vitamin D(3) (1alpha,25(OH)(2)D(3)) is the location of at least one activated vitamin D receptor (VDR) protein close to the transcription start site of the respective primary 1alpha,25(OH)(2)D(3) target gene. This is achieved through the specific binding of VDR to a 1alpha,25(OH)(2)D(3) response element (VDRE). Although these elements are well characterized in vitro, the function of VDREs in living cells in the context of chromatin is still largely unknown. To resolve this issue, approximately 8kB of the promoter regions of the primary 1alpha,25(OH)(2)D(3) target genes CYP24, cyclin C and p21((Waf1/Cip1)) were screened by chromatin immunoprecipitation (ChIP) assays for VDR binding sites using antibodies against VDR and its partner proteins. This approach identified three to four functional VDREs per gene promoter. In parallel, in silico screening of the extended gene areas (i.e. 10kB of promoter, introns, exons and 10kB of the downstream region) of all six members of the insulin-like growth factor binding protein (IGFBP) gene family was performed. Gel shift, reporter gene and ChIP assays identified in total 10 functional VDREs in the genes IGFBP1, IGFBP3 and IGFBP5. Taken together, both screening approaches suggest that a reasonable proportion of all VDR target genes, if not all, are under the control of multiple VDREs.

Time-resolved expression profiling of the nuclear receptor superfamily in human adipogenesis.

Background The differentiation of fibroblast-like pre-adipocytes to lipid-loaded adipocytes is regulated by a network of transcription factors, the most prominent one being the nuclear receptor peroxisome proliferator-activated receptor (PPAR) γ. However, many of the other 47 members of the nuclear receptor superfamily have an impact on adipogenesis, which in human cells has not been investigated in detail. Methodology/Principal Findings We analyzed by quantitative PCR all human nuclear receptors at multiple time points during differentiation of SGBS pre-adipocytes. The earliest effect was the down-regulation of the genes RARG, PPARD, REV-ERBA, REV-ERBB, VDR and GR followed by the up-regulation of PPARG, LXRA and AR. These observations are supported with data from 3T3-L1 mouse pre-adipocytes and primary human adipocytes. Investigation of the effects of the individual differentiation mix components in short-term treatments and of their omission from the full mix showed that the expression levels of the early-regulated nuclear receptor genes were most affected by the glucocorticoid receptor (GR) ligand cortisol and the phosophodiesterase inhibitor IBMX. Interestingly, the effects of both compounds converged to repress the genes PPARD, REV-ERBA, REV-ERBB, VDR and GR, whereas cortisol and IBMX showed antagonistic interaction for PPARG, LXRA and AR causing a time lag in their up-regulation. We hypothesize that the well-known auto-repression of GR fine-tunes the detected early responses. Consistently, chromatin immunoprecipitation experiments showed that GR association increased on the transcription start sites of the genes RARG, REV-ERBB, VDR and GR. Conclusions/Significance Adipocyte differentiation is a process, in which many members of the nuclear receptor superfamily change their mRNA expression. The actions of cortisol and IBMX converged to repress several nuclear receptors early in differentiation, while up-regulation of other nuclear receptor genes showed a time lag due to antagonisms of the signals. Our results place GR and its ligand cortisol as central regulatory factors controlling early regulatory events in human adipogenesis that precedes the regulation of the later events by PPARG.

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.

Combinatorial regulation of lipoprotein lipase by microRNAs during mouse adipogenesis

MicroRNAs (miRNAs) regulate gene expression directly through base pairing to their targets or indirectly through participating in multi-scale regulatory networks. Often miRNAs take part in feed-forward motifs where a miRNA and a transcription factor act on shared targets to achieve accurate regulation of processes such as cell differentiation. Here we show that the expression levels of miR-27a and miR-29a inversely correlate with the mRNA levels of lipoprotein lipase (Lpl), their predicted combinatorial target, and its key transcriptional regulator peroxisome proliferator-activated receptor gamma (Pparg) during 3T3-L1 adipocyte differentiation. More importantly, we show that Lpl, a key lipogenic enzyme, can be negatively regulated by the two miRNA families in a combinatorial fashion on the mRNA and functional level in maturing adipocytes. This regulation is mediated through the Lpl 3′UTR as confirmed by reporter gene assays. In addition, a small mathematical model captures the dynamics of this feed-forward motif and predicts the changes in Lpl mRNA levels upon network perturbations. The obtained results might offer an explanation to the dysregulation of LPL in diabetic conditions and could be extended to quantitative modeling of regulation of other metabolic genes under similar regulatory network motifs.