Deqiang Sun

Dnmt3a is essential for hematopoietic stem cell differentiation

Loss of the de novo DNA methyltransferases Dnmt3a and Dnmt3b in embryonic stem cells obstructs differentiation; however, the role of these enzymes in somatic stem cells is largely unknown. Using conditional ablation, we show that Dnmt3a loss progressively impairs hematopoietic stem cell (HSC) differentiation over serial transplantation, while simultaneously expanding HSC numbers in the bone marrow. Dnmt3a-null HSCs show both increased and decreased methylation at distinct loci, including substantial CpG island hypermethylation. Dnmt3a-null HSCs upregulate HSC multipotency genes and downregulate differentiation factors, and their progeny exhibit global hypomethylation and incomplete repression of HSC-specific genes. These data establish Dnmt3a as a critical participant in the epigenetic silencing of HSC regulatory genes, thereby enabling efficient differentiation.

Epigenomic profiling of young and aged HSCs reveals concerted changes during aging that reinforce self-renewal

To investigate the cell-intrinsic aging mechanisms that erode the function of somatic stem cells during aging, we have conducted a comprehensive integrated genomic analysis of young and aged cells. We profiled the transcriptome, DNA methylome, and histone modifications of young and old murine hematopoietic stem cells (HSCs). Transcriptome analysis indicated reduced TGF-β signaling and perturbation of genes involved in HSC proliferation and differentiation. Aged HSCs exhibited broader H3K4me3 peaks across HSC identity and self-renewal genes and showed increased DNA methylation at transcription factor binding sites associated with differentiation-promoting genes combined with a reduction at genes associated with HSC maintenance. Altogether, these changes reinforce HSC self-renewal and diminish differentiation, paralleling phenotypic HSC aging behavior. Ribosomal biogenesis emerged as a particular target of aging with increased transcription of ribosomal protein and RNA genes and hypomethylation of rRNA genes. This data set will serve as a reference for future epigenomic analysis of stem cell aging.

Large conserved domains of low DNA methylation maintained by Dnmt3a

Gains and losses in DNA methylation are prominent features of mammalian cell types. To gain insight into the mechanisms that promote shifts in DNA methylation and contribute to changes in cell fate, including malignant transformation, we performed genome-wide mapping of 5-methylcytosine and 5-hydroxymethylcytosine in purified mouse hematopoietic stem cells. We discovered extended regions of low methylation (canyons) that span conserved domains frequently containing transcription factors and are distinct from CpG islands and shores. About half of the genes in these methylation canyons are coated with repressive histone marks, whereas the remainder are covered by activating histone marks and are highly expressed in hematopoietic stem cells (HSCs). Canyon borders are demarked by 5-hydroxymethylcytosine and become eroded in the absence of DNA methyltransferase 3a (Dnmt3a). Genes dysregulated in human leukemias are enriched for canyon-associated genes. The new epigenetic landscape we describe may provide a mechanism for the regulation of hematopoiesis and may contribute to leukemia development.

MOABS: model based analysis of bisulfite sequencing data

Bisulfite sequencing (BS-seq) is the gold standard for studying genome-wide DNA methylation. We developed MOABS to increase the speed, accuracy, statistical power and biological relevance of BS-seq data analysis. MOABS detects differential methylation with 10-fold coverage at single-CpG resolution based on a Beta-Binomial hierarchical model and is capable of processing two billion reads in 24 CPU hours. Here, using simulated and real BS-seq data, we demonstrate that MOABS outperforms other leading algorithms, such as Fisher’s exact test and BSmooth. Furthermore, MOABS analysis can be easily extended to differential 5hmC analysis using RRBS and oxBS-seq. MOABS is available at http://code.google.com/p/moabs/.

RRBSMAP: a Fast, Accurate and User-friendly Alignment Tool for Reduced Representation Bisulfite Sequencing

SUMMARY Reduced representation bisulfite sequencing (RRBS) is a powerful yet cost-efficient method for studying DNA methylation on a genomic scale. RRBS involves restriction-enzyme digestion, bisulfite conversion and size selection, resulting in DNA sequencing data that require special bioinformatic handling. Here, we describe RRBSMAP, a short-read alignment tool that is designed for handling RRBS data in a user-friendly and scalable way. RRBSMAP uses wildcard alignment, and avoids the need for any preprocessing or post-processing steps. We benchmarked RRBSMAP against a well-validated MAQ-based pipeline for RRBS read alignment and observed similar accuracy but much improved runtime performance, easier handling and better scaling to large sample sets. In summary, RRBSMAP removes bioinformatic hurdles and reduces the computational burden of large-scale epigenome association studies performed with RRBS. AVAILABILITY http://rrbsmap.computational-epigenetics.org/ http://code.google.com/p/bsmap/ CONTACT wl1@bcm.tmc.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Long Non-Coding RNAs Control Hematopoietic Stem Cell Function

Hematopoietic stem cells (HSCs) possess unique gene expression programs that enforce their identity and regulate lineage commitment. Long non-coding RNAs (lncRNAs) have emerged as important regulators of gene expression and cell fate decisions, although their functions in HSCs are unclear. Here we profiled the transcriptome of purified HSCs by deep sequencing and identified 323 unannotated lncRNAs. Comparing their expression in differentiated lineages revealed 159 lncRNAs enriched in HSCs, some of which are likely HSC specific (LncHSCs). These lncRNA genes share epigenetic features with protein-coding genes, including regulated expression via DNA methylation, and knocking down two LncHSCs revealed distinct effects on HSC self-renewal and lineage commitment. We mapped the genomic binding sites of one of these candidates and found enrichment for key hematopoietic transcription factor binding sites, especially E2A. Together, these results demonstrate that lncRNAs play important roles in regulating HSCs, providing an additional layer to the genetic circuitry controlling HSC function.

BSeQC: quality control of bisulfite sequencing experiments.

MOTIVATION Bisulfite sequencing (BS-seq) has emerged as the gold standard to study genome-wide DNA methylation at single-nucleotide resolution. Quality control (QC) is a critical step in the analysis pipeline to ensure that BS-seq data are of high quality and suitable for subsequent analysis. Although several QC tools are available for next-generation sequencing data, most of them were not designed to handle QC issues specific to BS-seq protocols. Therefore, there is a strong need for a dedicated QC tool to evaluate and remove potential technical biases in BS-seq experiments. RESULTS We developed a package named BSeQC to comprehensively evaluate the quality of BS-seq experiments and automatically trim nucleotides with potential technical biases that may result in inaccurate methylation estimation. BSeQC takes standard SAM/BAM files as input and generates bias-free SAM/BAM files for downstream analysis. Evaluation based on real BS-seq data indicates that the use of the bias-free SAM/BAM file substantially improves the quantification of methylation level. AVAILABILITY AND IMPLEMENTATION BSeQC is freely available at: http://code.google.com/p/bseqc/.

Sexually Dimorphic Genome-Wide Binding of Retinoid X Receptor alpha (RXRα) Determines Male-Female Differences in the Expression of Hepatic Lipid Processing Genes in Mice

Many hepatic functions including lipid metabolism, drug metabolism, and inflammatory responses are regulated in a sex-specific manner due to distinct patterns of hepatic gene expression between males and females. Regulation for the majority of these genes is under control of Nuclear Receptors (NRs). Retinoid X Receptor alpha (RXRα) is an obligate partner for multiple NRs and considered a master regulator of hepatic gene expression, yet the full extent of RXRα chromatin binding in male and female livers is unclear. ChIP-Seq analysis of RXRα and RNA Polymerase2 (Pol2) binding was performed livers of both genders and combined with microarray analysis. Mice were gavage-fed with the RXR ligand LG268 for 5 days (30 mg/kg/day) and RXRα-binding and RNA levels were determined by ChIP-qPCR and qPCR, respectively. ChIP-Seq revealed 47,845 (male) and 46,877 (female) RXRα binding sites (BS), associated with ∼12,700 unique genes in livers of both genders, with 91% shared between sexes. RXRα-binding showed significant enrichment for 2227 and 1498 unique genes in male and female livers, respectively. Correlating RXRα binding strength with Pol2-binding revealed 44 genes being male-dominant and 43 female-dominant, many previously unknown to be sexually-dimorphic. Surprisingly, genes fundamental to lipid metabolism, including Scd1, Fasn, Elovl6, and Pnpla3-implicated in Fatty Liver Disease pathogenesis, were predominant in females. RXRα activation using LG268 confirmed RXRα-binding was 2–3 fold increased in female livers at multiple newly identified RXRα BS including for Pnpla3 and Elovl6, with corresponding ∼10-fold and ∼2-fold increases in Pnpla3 and Elovl6 RNA respectively in LG268-treated female livers, supporting a role for RXRα regulation of sexually-dimorphic responses for these genes. RXRα appears to be one of the most widely distributed transcriptional regulators in mouse liver and is engaged in determining sexually-dimorphic expression of key lipid-processing genes, suggesting novel gender- and gene-specific responses to NR-based treatments for lipid-related liver diseases.

Research Resource: The Estrogen Receptor α Cistrome Defined by DamIP

Gene expression is tightly regulated by transcription factors and cofactors that function by directly or indirectly interacting with DNA of the genome. Understanding how and where these proteins bind provides essential information to uncover genetic regulatory mechanisms. We have developed a new method to study DNA-protein interaction in vivo called DNA adenine methyltransferase (Dam)IP, which is based on fusing a protein of interest to a mutant form of Dam from Escherichia coli. We showed previously that DamIP can efficiently identify in vivo binding sites of Dam-tethered human estrogen receptor (hER)α. In current study, we present the cistrome of hERα determined by DamIP and high throughput sequencing (DamIP-seq). The DamIP-seq-defined hERα cistrome identifies many new binding regions and overlaps with those determined by chromatin immunoprecipitation (ChIP)-chip or ChIP-seq. Elements uniquely identified by DamIP-seq include a unique class of elements that show low, but persistent, hERα binding when reexamined by conventional ChIP. In contrast, DamIP-seq fails to detect some elements with very transient hERα binding. The methyl-adenine modifications introduced by Dam are stable and do not decrease over 12 d. In summary, the current study provides both an alternate view of the hERα cistrome to further understand the mechanism of hERα-mediated transcription and a new tool to explore other transcriptional factors and cofactors that is very different from conventional ChIP.