Xiaoyun Xing

Conserved role of intragenic DNA methylation in regulating alternative promoters.

Although it is known that the methylation of DNA in 5′ promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5′ CpG island (CGI) promoters, whereas a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences. Tissue-specific intragenic methylation might reduce, or, paradoxically, enhance transcription elongation efficiency. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes. To investigate the role of intragenic methylation, we generated a map of DNA methylation from the human brain encompassing 24.7 million of the 28 million CpG sites. From the dense, high-resolution coverage of CpG islands, the majority of methylated CpG islands were shown to be in intragenic and intergenic regions, whereas less than 3% of CpG islands in 5′ promoters were methylated. The CpG islands in all three locations overlapped with RNA markers of transcription initiation, and unmethylated CpG islands also overlapped significantly with trimethylation of H3K4, a histone modification enriched at promoters. The general and CpG-island-specific patterns of methylation are conserved in mouse tissues. An in-depth investigation of the human SHANK3 locus and its mouse homologue demonstrated that this tissue-specific DNA methylation regulates intragenic promoter activity in vitro and in vivo. These methylation-regulated, alternative transcripts are expressed in a tissue- and cell type-specific manner, and are expressed differentially within a single cell type from distinct brain regions. These results support a major role for intragenic methylation in regulating cell context-specific alternative promoters in gene bodies.

Principles of regulatory information conservation between mouse and human

To broaden our understanding of the evolution of gene regulation mechanisms, we generated occupancy profiles for 34 orthologous transcription factors (TFs) in human–mouse erythroid progenitor, lymphoblast and embryonic stem-cell lines. By combining the genome-wide transcription factor occupancy repertoires, associated epigenetic signals, and co-association patterns, here we deduce several evolutionary principles of gene regulatory features operating since the mouse and human lineages diverged. The genomic distribution profiles, primary binding motifs, chromatin states, and DNA methylation preferences are well conserved for TF-occupied sequences. However, the extent to which orthologous DNA segments are bound by orthologous TFs varies both among TFs and with genomic location: binding at promoters is more highly conserved than binding at distal elements. Notably, occupancy-conserved TF-occupied sequences tend to be pleiotropic; they function in several tissues and also co-associate with many TFs. Single nucleotide variants at sites with potential regulatory functions are enriched in occupancy-conserved TF-occupied sequences.

DNA hypomethylation within specific transposable element families associates with tissue-specific enhancer landscape

Transposable element (TE)-derived sequences comprise half of the human genome and DNA methylome and are presumed to be densely methylated and inactive. Examination of genome-wide DNA methylation status within 928 TE subfamilies in human embryonic and adult tissues identified unexpected tissue-specific and subfamily-specific hypomethylation signatures. Genes proximal to tissue-specific hypomethylated TE sequences were enriched for functions important for the relevant tissue type, and their expression correlated strongly with hypomethylation within the TEs. When hypomethylated, these TE sequences gained tissue-specific enhancer marks, including monomethylation of histone H3 at lysine 4 (H3K4me1) and occupancy by p300, and a majority exhibited enhancer activity in reporter gene assays. Many such TEs also harbored binding sites for transcription factors that are important for tissue-specific functions and showed evidence of evolutionary selection. These data suggest that sequences derived from TEs may be responsible for wiring tissue type–specific regulatory networks and may have acquired tissue-specific epigenetic regulation.

Exuberant expression of chemokine genes by adult human articular chondrocytes in response to IL-1β

OBJECTIVE To provide a more complete picture of the effect of interleukin-1 beta (IL-1beta) on adult human articular chondrocyte gene expression, in contrast to the candidate gene approach. DESIGN Chondrocytes from human knee cartilage were cultured in medium containing IL-1beta. Changes in gene expression were analyzed by microarray and reverse transcriptase-polymerase chain reaction analysis. The ability of transforming growth factor beta-1 (TGF-beta1), fibroblast growth factor (FGF)-18, and bone morphogenetic protein 2 (BMP-2) to alter the effects of IL-1beta was analyzed. Computational analysis of the promoter regions of differentially expressed genes for transcription factor binding motifs was performed. RESULTS IL-1beta-treated human chondrocytes showed significant increases in the expression of granulocyte colony stimulating factor-3, endothelial leukocyte adhesion molecule 1 and leukemia inhibitory factor as well as for a large group of chemokines that include CXCL1, CXCL2, CXCL3, CXCL5, CXCL6, CXCL8, CCL2, CCL3, CCL4, CCL5, CCL8, CCL20, CCL3L1, CX3CL1 and the cytokine IL-6. As expected, the mRNA for matrix metalloproteinase (MMP)-13 and BMP-2 also increased while mRNA for the matrix genes COL2A1 and aggrecan was down-regulated. A subset of chemokines increased rapidly at very low levels of IL-1beta. The phenotype induced by IL-1beta was partially reversed by TGF-beta1, but not by BMP-2. In the presence of IL-1beta, FGF-18 increased expression of ADAMTS-4, aggrecan, BMP-2, COL2A1, CCL3, CCL4, CCL20, CXCL1, CXCL3, CXCL6, IL-1beta, IL-6, and IL-8 and decreased ADAMTS-5, MMP-13, CCL2, and CCL8. Computational analysis revealed a high likelihood that the most up-regulated chemokines are regulated by the transcription factors myocyte enhancer binding factor-3 (MEF-3), CCAAT/enhancer binding protein (C/EBP) and nuclear factor-kappa B (NF-kappaB). CONCLUSION IL-1beta has a diverse effect on gene expression profile in human chondrocytes affecting matrix genes as well as chemokines and cytokines. TGF-beta1 has the ability to antagonize some of the phenotype induced by IL-1beta.

Widespread contribution of transposable elements to the innovation of gene regulatory networks

Transposable elements (TEs) have been shown to contain functional binding sites for certain transcription factors (TFs). However, the extent to which TEs contribute to the evolution of TF binding sites is not well known. We comprehensively mapped binding sites for 26 pairs of orthologous TFs in two pairs of human and mouse cell lines (representing two cell lineages), along with epigenomic profiles, including DNA methylation and six histone modifications. Overall, we found that 20% of binding sites were embedded within TEs. This number varied across different TFs, ranging from 2% to 40%. We further identified 710 TF-TE relationships in which genomic copies of a TE subfamily contributed a significant number of binding peaks for a TF, and we found that LTR elements dominated these relationships in human. Importantly, TE-derived binding peaks were strongly associated with open and active chromatin signatures, including reduced DNA methylation and increased enhancer-associated histone marks. On average, 66% of TE-derived binding events were cell type-specific with a cell type-specific epigenetic landscape. Most of the binding sites contributed by TEs were species-specific, but we also identified binding sites conserved between human and mouse, the functional relevance of which was supported by a signature of purifying selection on DNA sequences of these TEs. Interestingly, several TFs had significantly expanded binding site landscapes only in one species, which were linked to species-specific gene functions, suggesting that TEs are an important driving force for regulatory innovation. Taken together, our data suggest that TEs have significantly and continuously shaped gene regulatory networks during mammalian evolution.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Resistin induces expression of proinflammatory cytokines and chemokines in human articular chondrocytes via transcription and messenger RNA stabilization

OBJECTIVE To elucidate the effects of resistin on human articular chondrocytes and to generate a picture of their regulation at the transcriptional and posttranscriptional levels. METHODS Human articular chondrocytes were cultured with resistin. Changes in gene expression were analyzed at various doses and times. Cells were also treated with the transcription inhibitor actinomycin D after resistin treatment or with the NF-kappaB inhibitor IKK-NBD before resistin treatment. Gene expression was tested by quantitative real-time polymerase chain reaction. Computational analysis for transcription factor binding motifs was performed on the promoter regions of differentially expressed genes. TC-28 chondrocytes were transfected with CCL3 and CCL4 promoter constructs, pNF-kappaB reporter, and NF-kappaB and CCAAT/enhancer binding protein beta (C/EBPbeta) expression vectors with or without resistin. RESULTS Resistin-treated human articular chondrocytes increased the expression of cytokines and chemokines. Levels of messenger RNA (mRNA) for matrix metalloproteinase 1 (MMP-1), MMP-13, and ADAMTS-4 also increased, while type II collagen alpha1 (COL2A1) and aggrecan were down-regulated. The cytokine and chemokine genes could be categorized into 3 groups according to the pattern of mRNA expression over a 24-hour time course. One pattern suggested rapid regulation by mRNA stability. The second and third patterns were consistent with transcriptional regulation. Computational analysis suggested the transcription factors NF-kappaB and C/EBPbeta were involved in the resistin-induced up-regulation. This prediction was confirmed by the cotransfection of NF-kappaB and C/EBPbeta and the IKK-NBD inhibition. CONCLUSION Resistin has diverse effects on gene expression in human chondrocytes, affecting chemokines, cytokines, and matrix genes. Messenger RNA stabilization and transcriptional up-regulation are involved in resistin-induced gene expression in human chondrocytes.

Combining MeDIP-seq and MRE-seq to investigate genome-wide CpG methylation

DNA CpG methylation is a widespread epigenetic mark in high eukaryotes including mammals. DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and control of gene expression. Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to measure DNA methylation in a genome-wide fashion, making it possible to comprehensively investigate the role of DNA methylation. Several methods have been developed, such as Whole Genome Bisulfite Sequencing (WGBS), Reduced Representation Bisulfite Sequencing (RRBS), and enrichment-based methods including Methylation Dependent ImmunoPrecipitation followed by sequencing (MeDIP-seq), methyl-CpG binding domain (MBD) protein-enriched genome sequencing (MBD-seq), methyltransferase-directed Transfer of Activated Groups followed by sequencing (mTAG), and Methylation-sensitive Restriction Enzyme digestion followed by sequencing (MRE-seq). These methods differ by their genomic CpG coverage, resolution, quantitative accuracy, cost, and software for analyzing the data. Among these, WGBS is considered the gold standard. However, it is still a cost-prohibitive technology for a typical laboratory due to the required sequencing depth. We found that by integrating two enrichment-based methods that are complementary in nature (i.e., MeDIP-seq and MRE-seq), we can significantly increase the efficiency of whole DNA methylome profiling. By using two recently developed computational algorithms (i.e., M&M and methylCRF), the combination of MeDIP-seq and MRE-seq produces genome-wide CpG methylation measurement at high coverage and high resolution, and robust predictions of differentially methylated regions. Thus, the combination of the two enrichment-based methods provides a cost-effective alternative to WGBS. In this article we describe both the experimental protocols for performing MeDIP-seq and MRE-seq, and the computational protocols for running M&M and methylCRF.

Intermediate DNA methylation is a conserved signature of genome regulation

The role of intermediate methylation states in DNA is unclear. Here, to comprehensively identify regions of intermediate methylation and their quantitative relationship with gene activity, we apply integrative and comparative epigenomics to 25 human primary cell and tissue samples. We report 18,452 intermediate methylation regions located near 36% of genes and enriched at enhancers, exons and DNase I hypersensitivity sites. Intermediate methylation regions average 57% methylation, are predominantly allele-independent and are conserved across individuals and between mouse and human, suggesting a conserved function. These regions have an intermediate level of active chromatin marks and their associated genes have intermediate transcriptional activity. Exonic intermediate methylation correlates with exon inclusion at a level between that of fully methylated and unmethylated exons, highlighting gene context-dependent functions. We conclude that intermediate DNA methylation is a conserved signature of gene regulation and exon usage.

Mapping of Variable DNA Methylation across Multiple Cell Types Defines a Dynamic Regulatory Landscape of the Human Genome

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Many studies have mapped DNA methylation changes associated with embryogenesis, cell differentiation, and cancer at a genome-wide scale. Our understanding of genome-wide DNA methylation changes in a developmental or disease-related context has been steadily growing. However, the investigation of which CpGs are variably methylated in different normal cell or tissue types is still limited. Here, we present an in-depth analysis of 54 single-CpG-resolution DNA methylomes of normal human cell types by integrating high-throughput sequencing-based methylation data. We found that the ratio of methylated to unmethylated CpGs is relatively constant regardless of cell type. However, which CpGs made up the unmethylated complement was cell-type specific. We categorized the 26,000,000 human autosomal CpGs based on their methylation levels across multiple cell types to identify variably methylated CpGs and found that 22.6% exhibited variable DNA methylation. These variably methylated CpGs formed 660,000 variably methylated regions (VMRs), encompassing 11% of the genome. By integrating a multitude of genomic data, we found that VMRs enrich for histone modifications indicative of enhancers, suggesting their role as regulatory elements marking cell type specificity. VMRs enriched for transcription factor binding sites in a tissue-dependent manner. Importantly, they enriched for GWAS variants, suggesting that VMRs could potentially be implicated in disease and complex traits. Taken together, our results highlight the link between CpG methylation variation, genetic variation, and disease risk for many human cell types.