Epigenetic memory via concordant DNA methylation is inversely correlated to developmental potential of mammalian cells.

In storing and transmitting epigenetic information, organisms must balance the need to maintain information about past conditions with the capacity to respond to information in their current and future environments. Some of this information is encoded by DNA methylation, which can be transmitted with variable fidelity from parent to daughter strand. High fidelity confers strong pattern matching between the strands of individual DNA molecules and thus pattern stability over rounds of DNA replication; lower fidelity confers reduced pattern matching, and thus greater flexibility. Here, we present a new conceptual framework, Ratio of Concordance Preference (RCP), that uses double-stranded methylation data to quantify the flexibility and stability of the system that gave rise to a given set of patterns. We find that differentiated mammalian cells operate with high DNA methylation stability, consistent with earlier reports. Stem cells in culture and in embryos, in contrast, operate with reduced, albeit significant, methylation stability. We conclude that preference for concordant DNA methylation is a consistent mode of information transfer, and thus provides epigenetic stability across cell divisions, even in stem cells and those undergoing developmental transitions. Broader application of our RCP framework will permit comparison of epigenetic-information systems across cells, developmental stages, and organisms whose methylation machineries differ substantially or are not yet well understood.

Errors in the bisulfite conversion of DNA: modulating inappropriate- and failed-conversion frequencies.

Bisulfite treatment can be used to ascertain the methylation states of individual cytosines in DNA. Ideally, bisulfite treatment deaminates unmethylated cytosines to uracils, and leaves 5-methylcytosines unchanged. Two types of bisulfite-conversion error occur: inappropriate conversion of 5-methylcytosine to thymine, and failure to convert unmethylated cytosine to uracil. Conventional bisulfite treatment requires hours of exposure to low-molarity, low-temperature bisulfite ('LowMT') and, sometimes, thermal denaturation. An alternate, high-molarity, high-temperature ('HighMT') protocol has been reported to accelerate conversion and to reduce inappropriate conversion. We used molecular encoding to obtain validated, individual-molecule data on failed- and inappropriate-conversion frequencies for LowMT and HighMT treatments of both single-stranded and hairpin-linked oligonucleotides. After accounting for bisulfite-independent error, we found that: (i) inappropriate-conversion events accrue predominantly on molecules exposed to bisulfite after they have attained complete or near-complete conversion; (ii) the HighMT treatment is preferable because it yields greater homogeneity among sites and among molecules in conversion rates, and thus yields more reliable data; (iii) different durations of bisulfite treatment will yield data appropriate to address different experimental questions; and (iv) conversion errors can be used to assess the validity of methylation data collected without the benefit of molecular encoding.

Molecular barcodes detect redundancy and contamination in hairpin-bisulfite PCR.

PCR amplification of limited amounts of DNA template carries an increased risk of product redundancy and contamination. We use molecular barcoding to label each genomic DNA template with an individual sequence tag prior to PCR amplification. In addition, we include molecular 'batch-stamps' that effectively label each genomic template with a sample ID and analysis date. This highly sensitive method identifies redundant and contaminant sequences and serves as a reliable method for positive identification of desired sequences; we can therefore capture accurately the genomic template diversity in the sample analyzed. Although our application described here involves the use of hairpin-bisulfite PCR for amplification of double-stranded DNA, the method can readily be adapted to single-strand PCR. Useful applications will include analyses of limited template DNA for biomedical, ancient DNA and forensic purposes.

Why do fragile X carrier frequencies differ between Asian and non-Asian populations?

Asian and non-Asian populations have been reported to differ substantially in the distribution of fragile X alleles into the normal (< 55 CGG repeats), premutation (55-199 CGG repeats), and full-mutation (> 199 CGG repeats) size classes. Our statistical analyses of data from published general-population studies confirm that Asian populations have markedly lower frequencies of premutation alleles, reminiscent of earlier findings for expanded alleles at the Huntington's Disease locus. To examine historical and contemporary factors that may have shaped and now sustain allele-frequency differences at the fragile X locus, we develop a population-genetic/epigenetic model, and apply it to these published data. We find that founder-haplotype effects likely contribute to observed frequency differences via substantially lower mutation rates in Asian populations. By contrast, any premutation frequency differences present in founder populations would have disappeared in the several millennia since initial establishment of these groups. Differences in the reproductive fitness of female premutation carriers arising from fragile X primary ovarian insufficiency (FXPOI) and from differences in mean maternal age may also contribute to global variation in carrier frequencies.

Statistical inference of in vivo properties of human DNA methyltransferases from double-stranded methylation patterns.

DNA methyltransferases establish methylation patterns in cells and transmit these patterns over cell generations, thereby influencing each cell's epigenetic states. Three primary DNA methyltransferases have been identified in mammals: DNMT1, DNMT3A and DNMT3B. Extensive in vitro studies have investigated key properties of these enzymes, namely their substrate specificity and processivity. Here we study these properties in vivo, by applying novel statistical analysis methods to double-stranded DNA methylation patterns collected using hairpin-bisulfite PCR. Our analysis fits a novel Hidden Markov Model (HMM) to the observed data, allowing for potential bisulfite conversion errors, and yields statistical estimates of parameters that quantify enzyme processivity and substrate specificity. We apply this model to methylation patterns established in vivo at three loci in humans: two densely methylated inactive X (Xi)-linked loci (FMR1 and G6PD), and an autosomal locus (LEP), where methylation densities are tissue-specific but moderate. We find strong evidence for a high level of processivity of DNMT1 at FMR1 and G6PD, with the mean association tract length being a few hundred base pairs. Regardless of tissue types, methylation patterns at LEP are dominated by DNMT1 maintenance events, similar to the two Xi-linked loci, but are insufficiently informative regarding processivity to draw any conclusions about processivity at that locus. At all three loci we find that DNMT1 shows a strong preference for adding methyl groups to hemi-methylated CpG sites over unmethylated sites. The data at all three loci also suggest low (possibly 0) association of the de novo methyltransferases, the DNMT3s, and are consequently uninformative about processivity or preference of these enzymes. We also extend our HMM to reanalyze published data on mouse DNMT1 activities in vitro. The results suggest shorter association tracts (and hence weaker processivity), and much longer non-association tracts than human DNMT1 in vivo.

Testing the FMR1 promoter for mosaicism in DNA methylation among CpG sites, strands, and cells in FMR1-expressing males with fragile X syndrome.

Variability among individuals in the severity of fragile X syndrome (FXS) is influenced by epigenetic methylation mosaicism, which may also be common in other complex disorders. The epigenetic signal of dense promoter DNA methylation is usually associated with gene silencing, as was initially reported for FMR1 alleles in individuals with FXS. A paradox arose when significant levels of FMR1 mRNA were reported for some males with FXS who had been reported to have predominately methylated alleles. We have used hairpin-bisufite PCR, validated with molecular batch-stamps and barcodes, to collect and assess double-stranded DNA methylation patterns from these previously studied males. These patterns enable us to distinguish among three possible forms of methylation mosaicism, any one of which could explain FMR1 expression in these males. Our data indicate that cryptic inter-cell mosaicism in DNA methylation can account for the presence of FMR1 mRNA in some individuals with FXS.

STATISTICAL INFERENCE OF TRANSMISSION FIDELITY OF DNA METHYLATION PATTERNS OVER SOMATIC CELL DIVISIONS IN MAMMALS.

We develop Bayesian inference methods for a recently-emerging type of epigenetic data to study the transmission fidelity of DNA methylation patterns over cell divisions. The data consist of parent-daughter double-stranded DNA methylation patterns with each pattern coming from a single cell and represented as an unordered pair of binary strings. The data are technically difficult and time-consuming to collect, putting a premium on an efficient inference method. Our aim is to estimate rates for the maintenance and de novo methylation events that gave rise to the observed patterns, while accounting for measurement error. We model data at multiple sites jointly, thus using whole-strand information, and considerably reduce confounding between parameters. We also adopt a hierarchical structure that allows for variation in rates across sites without an explosion in the effective number of parameters. Our context-specific priors capture the expected stationarity, or near-stationarity, of the stochastic process that generated the data analyzed here. This expected stationarity is shown to greatly increase the precision of the estimation. Applying our model to a data set collected at the human FMR1 locus, we find that measurement errors, generally ignored in similar studies, occur at a non-trivial rate (inappropriate bisulfite conversion error: 1.6% with 80% CI: 0.9-2.3%). Accounting for these errors has a substantial impact on estimates of key biological parameters. The estimated average failure of maintenance rate and daughter de novo rate decline from 0.04 to 0.024 and from 0.14 to 0.07, respectively, when errors are accounted for. Our results also provide evidence that de novo events may occur on both parent and daughter strands: the median parent and daughter de novo rates are 0.08 (80% CI: 0.04-0.13) and 0.07 (80% CI: 0.04-0.11), respectively.

Encoding PCR products with batch-stamps and barcodes.

Polymerase chain reaction (PCR) has become the mainstay of DNA sequence analysis. Yet there is always uncertainty concerning the source of the template DNA that gave rise to a particular PCR product. The risks of contamination, biased amplification, and product redundancy are especially high when limited amounts of template DNA are used. We have developed and applied molecular encoding principles to solve this source-uncertainty problem for DNA sequences generated by standard PCR. Batch-stamps specify the date and sample identity, and barcodes detect template redundancy. Our approach thus enables classification of each PCR-derived sequence as valid, contaminant, or redundant, and provides a measure of sequence diversity. We recommend that batch-stamps and barcodes be used when amplifying irreplaceable DNAs and cDNAs available for forensic, clinical, single cell, and ancient DNA analyses.

A population-epigenetic model to infer site-specific methylation rates from double-stranded DNA methylation patterns.

Cytosine methylation is an epigenetic mechanism in eukaryotes that is often associated with stable transcriptional silencing, such as in X-chromosome inactivation and genomic imprinting. Aberrant methylation patterns occur in several inherited human diseases and in many cancers. To understand how methylated and unmethylated states of cytosine residues are transmitted during DNA replication, we develop a population-epigenetic model of DNA methylation dynamics. The model is informed by our observation that de novo methylation can occur on the daughter strand while leaving the opposing cytosine unmethylated, as revealed by the patterns of methylation on the two complementary strands of individual DNA molecules. Under our model, we can infer site-specific rates of both maintenance and de novo methylation, values that determine the fidelity of methylation inheritance, from double-stranded methylation data. This approach can be used for populations of cells obtained from individuals without the need for cell culture. We use our method to infer cytosine methylation rates at several sites within the promoter of the human gene FMR1.

Hemimethylation and non-CpG methylation levels in a promoter region of human LINE-1 (L1) repeated elements.

DNA methylation within the promoter region of human LINE1 (L1) transposable elements is important for maintaining transcriptional inactivation and for inhibiting L1 transposition. Determining methylation patterns on the complementary strands of repeated sequences is difficult using standard bisulfite methylation analysis. Evolutionary changes in each repeat and the variations between cells or alleles of the same repeat lead to a heterogeneous population of sequences. Potential sequence biases can arise during analyses that are different for the converted sense and antisense strands. These problems can be avoided with hairpin-bisulfite PCR, a double-stranded PCR method in which complementary strands of individual molecules are attached by a hairpin linker ligated to genomic DNA. Using human L1 elements to study methylation of repeated sequences, (i) we distinguish valid L1 sequences from redundant and contaminant sequences by applying the powerful new method of molecular barcodes, (ii) we resolve a controversy on the level of hemimethylation of L1 sequences in fetal fibroblasts in favor of relatively little hemimethylation, (iii) we report that human L1 sequences in different cell types also have primarily concordant CpG methylation patterns on complementary strands, and (iv) we provide evidence that non-CpG cytosines within the regions analyzed are rarely methylated.

Hairpin-bisulfite PCR: assessing epigenetic methylation patterns on complementary strands of individual DNA molecules.

Epigenetic inheritance, the transmission of gene expression states from parent to daughter cells, often involves methylation of DNA. In eukaryotes, cytosine methylation is a frequent component of epigenetic mechanisms. Failure to transmit faithfully a methylated or an unmethylated state of cytosine can lead to altered phenotypes in plants and animals. A central unresolved question in epigenetics concerns the mechanisms by which a locus maintains, or changes, its state of cytosine methylation. We developed "hairpin-bisulfite PCR" to analyze these mechanisms. This method reveals the extent of methylation symmetry between the complementary strands of individual DNA molecules. Using hairpin-bisulfite PCR, we determined the fidelity of methylation transmission in the CpG island of the FMR1 gene in human lymphocytes. For the hypermethylated CpG island of this gene, characteristic of inactive-X alleles, we estimate a maintenance methylation efficiency of approximately 0.96 per site per cell division. For de novo methylation efficiency (E(d)), remarkably different estimates were obtained for the hypermethylated CpG island (E(d) = 0.17), compared with the hypomethylated island on the active-X chromosome (E(d) < 0.01). These results clarify the mechanisms by which the alternative hypomethylated and hypermethylated states of CpG islands are stably maintained through many cell divisions. We also analyzed a region of human L1 transposable elements. These L1 data provide accurate methylation patterns for the complementary strand of each repeat sequence analyzed. Hairpin-bisulfite PCR will be a powerful tool in studying other processes for which genetic or epigenetic information differs on the two complementary strands of DNA.