Genomic Landscape of Human Allele-Specific DNA Methylation
Total Page:16
File Type:pdf, Size:1020Kb
Genomic landscape of human allele-specific DNA methylation Fang Fanga, Emily Hodgesb, Antoine Molarob, Matthew Deana, Gregory J. Hannonb, and Andrew D. Smitha,1 aMolecular and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, California; and bHoward Hughes Medical Institute, Cold Spring Harbor Laboratory, 1 Bungtown Road, Cold Spring Harbor, NY 11724 Edited by Wing Hung Wong, Stanford University, Stanford, CA, and approved March 8, 2012 (received for review January 24, 2012) DNA methylation mediates imprinted gene expression by passing stages preceding the context in which they become active. Such an epigenomic state across generations and differentially marking methods have been successfully applied to identify unique im- specific regulatory regions on maternal and paternal alleles. Im- printed genes (14–17). printing has been tied to the evolution of the placenta in mammals Advances in DNA sequencing technology have been leveraged and defects of imprinting have been associated with human dis- for high-throughput identification of imprinted genes. The “BS- eases. Although recent advances in genome sequencing have revo- seq” technology couples bisulfite treatment with high-throughput lutionized the study of DNA methylation, existing methylome data short-read sequencing, and has enabled genome-wide profiling of remain largely untapped in the study of imprinting. We present a DNA methylation in mammalian genomes at single-CpG (cyto- statistical model to describe allele-specific methylation (ASM) in sine guanine dinucleotide) resolution (18). Li et al. (19) produced data from high-throughput short-read bisulfite sequencing. Simula- a methylome from peripheral blood of a single individual and tion results indicate technical specifications of existing methylome recognized the potential of using such data to profile ASM. They data, such as read length and coverage, are sufficient for full- employed a method based on associating heterozygous SNPs genome ASM profiling based on our model. We used our model to with differential methylation, and identified hundreds of ASM analyze methylomes for a diverse set of human cell types, including regions. Methods such as this, however, must be applied to data cultured and uncultured differentiated cells, embryonic stem cells from a single individual and for which matching genotypic data and induced pluripotent stem cells. Regions of ASM identified most are available. There are two shortcomings of approaches that BIOPHYSICS AND consistently across methylomes are tightly connected with known depend on genotype. First, they can be confounded by ASM that imprinted genes and precisely delineate the boundaries of several is associated with genotype, but which may not have any regula- COMPUTATIONAL BIOLOGY known imprinting control regions. Predicted regions of ASM com- tory effect. The amount of ASM typically associated with geno- mon to multiple cell types frequently mark noncoding RNA promo- type is not well understood, but recent reports suggest it is ters and represent promising starting points for targeted validation. significant (20). More importantly, because imprinted methyla- More generally, our model provides the analytical complement to tion is not necessarily associated with genotypic variation, these cutting-edge experimental technologies for surveying ASM in spe- methods will be inherently blind to some portion of ASM. cific cell types and across species. We present a probabilistic model to describe ASM based on data from BS-seq experiments. Our model is independent of enomic imprinting refers to genes that are preferentially ex- genotype, and therefore has broad applicability to identify ASM Gpressed from either the maternal or paternal allele without in the context of imprinting. In essence, our model describes the genotype dependence (1). In mammals, such parent-of-origin gene degree to which methylation states in reads appear to reflect expression is believed to have evolved along with the placenta, ser- two distinct patterns, each pattern representing roughly half the ving to mediate resource distribution between a mother and her data. We validated our method using semisimulated data in which offspring (2, 3), though other theories have been proposed (4–6). methylation states were simulated within actual reads from BS- The connection between imprinting and DNA methylation was seq experiments. Our results indicate that technical characteris- uncovered shortly after the first identification of imprinted genes tics of existing public methylomes (i.e., read length and coverage) in mammals (7). Imprinted gene expression, in all known cases, is are sufficient to accurately identify AMRs. By applying our model regulated by allele-specific methylation (ASM) of some cis-acting to 22 human methylomes, emphasizing those from uncultured regulatory regions. We use the term allelically methylated region cells, we identified a set of candidate AMRs involved in im- (AMR) in reference to any genomic interval of ASM, whether or printed gene regulation. Candidates consistently identified across not it is associated with imprinted regulation. Typically, an entire methylomes display remarkable concordance with known im- imprinted locus is organized as a cluster and regulated by an printed genes and allow boundaries of known AMRs to be pre- imprinting control region (ICR) and several other AMRs. The cisely defined. Many candidates not associated with known allelic methylation patterns of ICRs are set during gametogenesis imprinted genes mark the promoters of long noncoding RNAs and stably maintained throughout somatic development in the (lncRNAs) and are also supported by similar analyses at ortho- offspring (8), irrespective of gene expression levels. The remain- logous regions in chimp; these provide a starting point for iden- ing AMRs may be established after fertilization (9), possibly un- tifying additional imprinted genes, ICRs, and possibly imprinted der the control of nearby ICRs or other epigenetic signals. clusters. Our model, therefore, is an essential analytical comple- The identification of imprinted genes and a detailed under- ment to recently emerged experimental methods for understand- standing of their regulation has become increasingly important, ing the role of DNA methylation in genomic imprinting. along with the realization that aberrant genomic imprinting contributes to several complex diseases (10). Much effort has Author contributions: F.F., E.H., A.M., M.D., G.J.H., and A.D.S. designed research; F.F., E.H., been directed toward locating imprinted genes using expression A.M., and A.D.S. performed research; F.F. and A.D.S. analyzed data; and F.F. and A.D.S. screen-based approaches (11, 12). One limitation of such ap- wrote the paper. proaches is that many imprinted genes may only show allele- The authors declare no conflict of interest. specific expressions in particular tissues at appropriate develop- This article is a PNAS Direct Submission. mental stages (13). ASM screen-based approaches might over- 1To whom correspondence should be addressed. E-mail: [email protected]. come the effect of temporal and spatial expression patterns This article contains supporting information online at www.pnas.org/lookup/suppl/ because the ICRs are expected to exist through developmental doi:10.1073/pnas.1201310109/-/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.1201310109 PNAS Early Edition ∣ 1of6 Downloaded by guest on September 27, 2021 Yn Y2 Modeling Allele-Specific Methylation in BS-Seq Data R mðγj;iÞ L2 Θ R; γ j j 0 5jRj θ 1 − θ uðγj;iÞ; [3] We begin this section with a verbal description of our question ð j Þ¼ γ . ij ð ijÞ j 1j i 1 j 1 and the main issues that are addressed by our model. We assume ¼ ¼ any read has been sequenced after bisulfite treatment and where the m and u are as defined for Eq. 1. Because the allele of mapped uniquely to the reference genome. Because we are inter- origin for each read is missing data, we fit the two-allele model ested in mammalian methylation, we restrict our attention to using expectation maximization (21), obtaining expectations on CpG sites both in the genome and in the reads. Reads not map- membership in γ1 and γ2. Details are provided in SI Text. ping over a CpG are ignored. Our goal is to identify intervals of the genome where it appears that the two alleles have different Identifying Intervals of Allele-Specific Methylation. We use Bayesian methylation patterns—typically, in such a case, one allele will be information criterion (BIC) (22) as a model selection criterion highly methylated and the other not. There are two kinds of im- in determining whether a fixed interval is best described using portant information our model must capture: (i) The set of reads a single-allele [Eq. 1] or two-allele [Eq. 2] model. A single-allele mapping into the interval should appear to represent two distinct model has one parameter for each of the n CpGs, and the number methylation patterns, and (ii) the subsets of reads corresponding of observations is equal to jRj: to those two patterns should be in roughly equal proportions be- n R − 2 L Θ R : [4] cause the alleles themselves are present in equal proportions. BICðsingleÞ¼ ln j j ln 1ð j Þ One can consider a methylation pattern as analogous to a hap- For the two-allele model, there are two parameters for each CpG: lotype, but with a strong stochastic component. Therefore, reads that contain only a single CpG will provide us with relatively little