Mother-child transmission of epigenetic information by tunable polymorphic imprinting.

Genomic imprinting mediated by DNA methylation restricts gene expression to a single allele determined by parental origin and is not generally considered to be under genetic or environmental influence. Here, we focused on a differentially methylated region (DMR) of approximately 1.9 kb that includes a 101-bp noncoding RNA gene (nc886/VTRNA2-1), which is maternally imprinted in ∼75% of humans. This is unlike other imprinted genes, which demonstrate monoallelic methylation in 100% of individuals. The DMR includes a CTCF binding site on the centromeric side defining the DMR boundary and is flanked by a CTCF binding site on the telomeric side. The centromeric CTCF binding site contains an A/C polymorphism (rs2346018); the C allele is associated with less imprinting. The frequency of imprinting of the nc886 DMR in infants was linked to at least two nongenetic factors, maternal age at delivery and season of conception. In a separate cohort, nc886 imprinting was associated with lower body mass index in children at 5 y of age. Thus, we propose that the imprinting status of the nc886 DMR is "tunable" in that it is associated with maternal haplotype and prenatal environment. This provides a potential mechanism for transmitting information, with phenotypic consequences, from mother to child.

Epigenome-wide association studies (EWASs) have used the Infinium Human Methylation 450 (HM450) array to identify CpG sites that are differentially methylated in the DNA of individuals with differing phenotypes or diseases. Whereas most studies are cross-sectional, some have sought to identify risk-associated epigenetic marks that hold predictive value (1, 2). It has been known for many years that genetic polymorphisms can alter the frequency of DNA methylation, that is, they can alter allele-specific DNA methylation (ASM) (3–5). DNA methylation is also associated with environmental perturbations (6–10): for example, maternal and paternal exposures such as nutrition and smoking are linked to the methylation states of CpG sites in children (11–15).

Genomic imprinting is a type of ASM that is usually associated with extended regions of parental-specific differentially methylated DNA [i.e., differentially methylated regions (DMRs)] and is established in the gametes. Most of the approximately 100 imprinted genes are monoallelically methylated at imprinting control regions, where the gene is expressed from the unmethylated allele in nearly all somatic tissues at a given developmental stage (16). The exceptions include a small subset of genes that are imprinted in a tissue- or isoform-specific manner, in brain and placenta, across most individuals (16, 17). Importantly, unlike nonimprinted ASMs, imprinting is not generally altered by genetic or nongenetic factors, but is instead fixed in the population (18). Tight control of nearly all imprinted genes across the human population is consistent with the known importance of imprinting in regulating embryonic development and neuronal function (19, 20). Disruptions in imprinting can lead to human disorders, including Beckwith–Wiedemann, Silver–Russel, Prader–Willi, and Angelman syndromes, all of which have severe phenotypes impacting human growth and development (21, 22).

Work performed in the P.A.J. and K.G. laboratories uncovered a unique pattern of DNA methylation at the promoter of a noncoding gene (nc886) whereby ∼75% of individuals exhibited ASM and 25% were completely unmethylated (23). Romanelli et al. (24) expanded on this finding to determine that this region of DNA was polymorphically imprinted in a small sample set. Before the conclusion that nc886 was polymorphically imprinted, interindividual variation of imprinting status had been reported only for placental-specific imprinted genes (25, 26).

Genomic imprinting is considered an all-or-none process, in which DMRs on the maternal or paternal autosomes cause monoallelic expression of the corresponding gene (27, 28). Imprinting is commonly, but not always, established by silencing of a single allele through DNA methylation (29). Maternally methylated imprinted genes are often associated with decreased size of the developing fetus (30), whereas paternally expressed imprinted genes are often associated with increased fetal weight (30). Importantly, maternal exposures during pregnancy shift the percentage of DNA methylation or level of expression at imprinted genes, primarily in the placenta (24, 31, 32). The best studied example of genetic and environmental impacts on imprinted genes is that of the H19/IGF2 locus, where quantitatively small but additive effects on levels of DNA methylation are seen (33). However, the possibility that a parental cue from the environment or parental genetics might shift the likelihood that genomic imprinting is established in a human child has yet to be investigated.

Here, we studied the frequency of imprinting of a RNA polymerase III-transcribed, noncoding RNA of 101 bp called nc886 (also referred to as VTRNA2-1), which can be silenced by DNA methylation (23, 34). The complete nc886 sequence is found in higher-order primates, although portions of it are present in other mammals (), which may indicate its importance in primate evolution. ASM at the nc886 locus has previously been studied and described as a human “metastable epiallale,” which is an allele that shows epigenetic heterogeneity in a population and is sensitive to environmental conditions (6, 35). However, an appreciation of the imprinting biology of this locus has not been emphasized in prior studies. In this work, we show that the region surrounding nc886 acts as an imprinted DMR, as opposed to an intermediately methylated region. This locus represents an instance of nonplacental polymorphic imprinting in humans (24), and we show that variation of imprinting in the population is associated with prenatal environment.

nc886 Is Part of a 1.9-kb Polymorphically Imprinted Region

Genomic imprinting occurs as parental-specific DMRs that tightly regulate gene dosage (28). Nonimprinted ASMs, on the contrary, depend on genomic context rather than parental origin, and on local polymorphisms to affect CpG methylation states, which are generally confined to shorter regions of DNA (3, 18, 36). We and others have reported that nc886 expression is tightly regulated by DNA methylation and that the nc886 locus shows polymorphic imprinting in the human population when measured in peripheral blood cells: approximately 75% of individuals worldwide are monoallelically methylated and approximately 25% are biallelically unmethylated at nc886 (23, 24, 34, 37).

We first used HM450 data to examine DNA methylation at known imprinted DMRs (Fig. 1) (38). We calculated the variance of β-values in 26 imprinted DMRs that had a minimum of 5 HM450 probes by using data from peripheral blood of 2,664 European and Indian Asian subjects from the London Life Sciences Prospective Population (LOLIPOP) study (39). We confirm that all imprinted DMRs investigated, except for probes adjacent to nc886, have limited interindividual variance in this cohort. This variance is much higher for nc886 than for other included imprinted DMRs (Fig. 1).

Fig. 1.

The nc886 DMR displays maternal polymorphic imprinting. (A) Populational variance for β-values in known imprinting control regions in peripheral blood. Each dot corresponds to one CpG interrogated by an HM450 probe. Populational SDs (y-axis) are plotted for CpGs located in known imprinting DMRs and nc886 (x-axis). Although all imprinting DMRs were examined, only imprinting control regions with more than five probes are included in the display and sorted by the mean SD in the population. Imprinting control regions were retrieved from an earlier study (38). Probes with multimodal distribution of β-values in the population are colored red. DNA methylation data were retrieved from the LOLIPOP study (39). (B, E, and H) DNA methylation data from the Infinium HM450 BeadChip platform in 66 sperm samples (GEO accession numbers GSE47627 and GSE64096) for a paternally methylated DMR for the H19 gene (B), a maternally methylated DMR for the PEG3 gene (E), and the nc886 DMR (H). Black lines point to the genomic location of each probe relative to the respective DMR. (C, F, and I) Heat map of HM450 β-values from peripheral blood for H19 (C), PEG3 (F), and nc886 (I); data obtained from Wahl et al. (39). (D, G, and J) β-Values of each sample are plotted for Infinium HM450 probes for H19 (D) PEG3 (G), and nc886 (J). HM450 probe IDs are provided in . Only probes with sequence mapped optimally are included.

We chose to survey DNA methylation in more detail for known imprinted loci containing a maternally methylated DMR (PEG3), a paternally methylated DMR (H19), and, to expand our window of investigation of nc886, a polymorphically imprinted locus. We analyzed data from two independent studies of sperm DNA (GSE47627 and GSE64096) and also peripheral blood from the LOLIPOP study (39–41). The H19 DMR is fully methylated and the PEG3 DMR generally lacks DNA methylation in 66 sperm samples (Fig. 1 , respectively). We found that the nc886 locus and its flanking regions were unmethylated in all sperm samples, similar to PEG3, a maternally methylated DMR (Fig. 1, dark blue), with partial methylation at a CTCF site on the telomeric side. As expected of imprinted loci, we saw ∼50% methylation in all peripheral blood samples for H19 and PEG3 (Fig. 1 , green). For nc886 in peripheral blood, the majority of samples showed approximately 50% methylation (green; Fig. 1) and the remainder showed methylation of nearly 0% (blue; Fig. 1). Based on these HM450 data, we defined the nc886 region as a maternally derived DMR in the human population with boundaries marked by two CTCF sites (spaced ∼3 kb apart; Fig. 1 ). Consistent with previous work from Monk and coworkers (24), whole-genome bisulfite sequencing data of various normal tissues from The Cancer Genome Atlas (TCGA) (42) and BLUEPRINT (43, 44) () better resolves the boundaries of the nc886 DMR, which is 1,979 bp in size and includes the nc886 locus. nc886 is overlapped by a CpG island, and the DMR is bounded on the centromeric side by repetitive sequences just outside a variably methylated CTCF site (Fig. 1 ) (24, 45). The nc886 DMR is flanked on the telomeric side by a CTCF site that is unmethylated in peripheral blood (Fig. 1).

We noted that three type II probes (cg04515200, cg13581155, and cg11978884) within the nc886 DMR on the HM450 platform displayed a consistent bias toward hypomethylation relative to the rest of the probes in the DMR (Fig. 1). We believe this is inherent to the Infinium probe design because we could still distinguish a separation in the methylation β-values for individuals who had monoallelic methylation at these probes, and, based on whole-genome bisulfite sequencing data, the three CpGs do not appear biased from an averaged methylation level of 0.5 (Fig. 1 and ). The high correlation of β-values among these three probes and probes across the nc886 DMR also support the idea that these three CpGs bear consistent methylation with the entire DMR and that this DMR is acting as a contiguous unit ().

Imprinting is essentially dichotomous: DMRs are monoallelically methylated or not. As such, methylation data cannot be considered a continuous variable in which the methylation values of all of the assayed DNA molecules are averaged. We reexamined the LOLIPOP data set, taking this biology into account. For the H19 and PEG3 DMR, all data points were near a β-value of 0.5, as expected (Fig. 1 ). However, for nc886, the dichotomous nature of the data became clear, in that the nc886 DMR showed β-values primarily at 50% or near 0% with the exception of the three CpGs previously mentioned (Fig. 1). Thus, by using the data from the large cohort (39), we have confirmed the previous observations that the nc886 DMR is polymorphically imprinted in humans (23, 24).

An A/C SNP Is Associated with Local DNA Methylation Density

It has been unclear whether the polymorphic nature of imprinting of the nc886 DMR is governed by genetic or nongenetic factors, given the lack of evidence in the literature that genomic imprinting can be associated with the local haplotype (5). Similar to Van Baak et al. (35), we reanalyzed DNA methylation data from twins generated by Grundberg et al. (46). We confirmed their findings that monozygotic twins, but not dizygotic twins, are concordant for DNA methylation across the nc886 DMR (), strongly suggesting that the nc886 DMR could be influenced by genetic factors. The nc886 locus is located within a haplotype block ranging in size from ∼5 kb to 35 kb, depending on the population (). The 2 kb spanning the nc886 DMR has been examined for potentially causal SNPs without success or with the conclusion that SNPs did not impact DNA methylation (6, 23, 24, 37, 45). Given the central role of CTCF in imprinting, we focused on an A/C SNP (rs2346018; A allele frequency range, 31–45%; C allele frequency range, 54–68%; from the 1000 Genomes Project) located in the variably methylated centromeric CTCF binding site (36, 47, 48). By using the R package motifbreakR, we found that having a C SNP “breaks” the CTCF binding motif of the centromeric CTCF site and could therefore impact the ability of CTCF to bind its target (49). To determine if SNP status was associated with local DNA methylation density, we analyzed DNA methylation at the centromeric CTCF site and across the nc886 locus in peripheral blood DNA of 31 cancer-free individuals by using bisulfite conversion followed by clonal sequencing (primer locations shown in Fig. 2). Eight examples are shown in Fig. 2, with the data from the remaining 23 individuals shown in . Clearly, the A and C alleles can be densely methylated in homozygotes and heterozygotes, but there was also sporadic methylation in predominantly methylated or unmethylated DNA strands (Fig. 2).

Fig. 2.

An A polymorphism (rs2346018) is associated with higher local DNA methylation density at the centromeric CTCF site, which has a strong positive correlation to DNA methylation at nc886. (A) Diagram of the nc886 DMR. Black arrows indicate primer locations for bisulfite sequencing in Figs. 2 and 3 and . (B) Bisulfite conversion and clonal sequencing of genomic DNA from white blood cells of eight individuals at the centromeric CTCF site and nc886. Asterisk indicates SNP location. (C) Percentage of CpGs methylated for each allele at the centromeric CTCF site in homozygotes and heterozygotes (n = 31; P < 1 × 10−9). (D) DNA methylation of the centromeric CTCF plotted against DNA methylation of nc886 as a discrete variable, with 0 being biallelically unmethylated and 1 being monoallelically methylated (n = 31; P = 4.7 × 10−7). Gray ribbon represents 95% CI.

Fig. 3.

Maternal imprinting of the nc886 region extends to the centromeric CTCF site. Locus-specific bisulfite sequencing was performed on genomic DNA at the centromeric CTCF site. Genomic DNA was isolated from benign lymphoblastoid cell lines derived from parent–offspring trios (mother, father, child; A and B are representative of two different trios) of disease-free individuals, in which SNP status can be determined from the sequence. Two representative parent–offspring trios are shown. Clones are sorted based first on SNP status and then on DNA methylation.

We chose to test the hypothesis that genetic background could “tune” the likelihood that an individual displays imprinting at the nc886 DMR, thereby switching from imprinted to not imprinted. Examining data from heterozygous and homozygous individuals, we determined that the A allele has higher average methylation density than the C allele (Fig. 2; P < 1 × 10−9). Furthermore, the odds ratio that a clone of the centromeric CTCF site is fully methylated in an individual with the A polymorphism (vs. C) was 2.86. These data demonstrate that DNA methylation of the centromeric CTCF site, 900 bp from nc886, correlates with a common SNP, with an A allele more likely to be methylated than a C allele. Given the role of CTCF in imprinting, the presence of the A allele might directly alter the likelihood of establishing DNA methylation at the nc886 DMR through reduced sequence binding affinity of CTCF or other DNA-binding proteins (e.g., DNMTs, ZFP57, Kaiso) (50, 51). Alternatively, the A allele might indirectly impact imprinting through increased DNA methylation density, and therefore binding affinity of CTCF (52) (see Fig. 5).

Fig. 5.

Working model for polymorphic imprinting of the nc886 locus. All individuals we examined have the single paternally derived unmethylated allele at the nc886 DMR. In 75% of the human population, the nc886 DMR is monoallelically methylated on the maternally derived allele (Bottom), with the remaining 25% having no methylation on either allele (Top). Here, we demonstrate that the presence of an A or C polymorphism in the centromeric CTCF binding site is associated with local DNA methylation density in the nc886 DMR. We hypothesize that the presence of the A allele or increased DNA methylation density reduces CTCF binding, allowing DNA methylation to spread into the DMR. We reanalyzed data from Silver et al. (6) and Markunas et al. (54) and found that maternal age and season of conception contribute to the likelihood that a child will have imprinting at the nc886 DMR. Therefore, we propose that a combination of maternal factors shifts the balance between imprinted and not-imprinted status in children, with downstream phenotypic consequences such as BMI.

Linked DNA Methylation of the Centromeric CTCF Binding Site and nc886

Previous work has shown that DNA methylation at the centromeric CTCF site serves as the boundary of ASM for the nc886 region (24, 45). Because the entire 1.9-kb region between the CTCF sites acts as a DMR (Fig. 1 and ), we hypothesized that DNA methylation of the centromeric CTCF site, or of any region within the DMR, might predict DNA methylation at nc886. Although other studies have suggested that the entire region can be imprinted (23, 32), they have not shown a lack of DNA methylation of this region in sperm, as we have (Fig. 1), and it is technically challenging to show that DNA methylation is linked across a single DNA strand.

Fig. 2 shows clonal bisulfite sequencing analysis of genomic DNA from individuals who have monoallelic methylation or no methylation, as indicated by SNP status in the centromeric CTCF site. We classified 50% methylation, representative of monoallelic methylation at nc886, as being “1” and 0% methylation, or biallelic lack of methylation, as being “0,” and modeled the data against the percentage of DNA methylation at the centromeric CTCF site for all 31 individuals (Fig. 2). Lack of DNA methylation at the centromeric CTCF site was associated with a very low probability of having DNA methylation at nc886. With increasing density of DNA methylation at the centromeric CTCF site, the probability of nc886 being classified as “imprinted” increased; thus, there was a strong positive correlation between DNA methylation at the centromeric CTCF site and at nc886 (P < 0.0001; R2 = 0.77; Fig. 2). DNA methylation of the centromeric CTCF site therefore explains much of the variation in DNA methylation at nc886. Although we were not able to sequence individual DNA strands greater than 1 kb in length, these data and the high correlation between β-values across the nc886 DMR () suggest that DNA methylation of this region is indeed present on individual alleles (Fig. 2).

DNA Methylation of the Centromeric CTCF Binding Site Is Maternally Derived

As summarized in , previous studies analyzing SNPs have independently concluded that DNA methylation across the nc886 region is not dependent on genetic context but is maternally derived (23, 24, 37, 45). Additionally, when analyzing HM450 data from sperm, we find that the nc886 DMR lacks DNA methylation, indicating that methylation of this region is likely maternally derived (Fig. 1). To confirm that DNA methylation of the centromeric CTCF site could be included in this region of maternally derived methylation, we used parent–offspring trios informative for the A/C SNP (rs2346018) and performed bisulfite conversion and clonal sequencing. As shown in Fig. 3, DNA methylation of the centromeric CTCF site was maternally derived. Thus, the existing literature and the lack of DNA methylation across the nc886 DMR in sperm (Fig. 1), along with this result, supports that imprinting of the nc886 DMR is not paternally derived.

Maternal Age and Nutrition Are Associated with the nc886 DMR

Recent EWAS studies have identified nc886 as a “metastable epiallele” that is altered in the population and is dependent on maternal nutritional status (6, 35, 53). However, these studies did not consider the data in terms of the percentage of each cohort demonstrating imprinting, which might account for the biology. Therefore, we reanalyzed the results from three independent studies in terms of percentage of each subpopulation that exhibits imprinting to determine if the likelihood that an individual is called “imprinted” shifts in each experimental group (). To address global shifts in β-values introduced by different methods used to preprocess the HM450 data or by residual experimental batch effects, we applied unsupervised hierarchical clustering, using the DNA methylation β-values from the entire region to determine whether an individual had imprinting at the nc886 DMR and whether this call is conclusive.

We used the valuable data from Silver et al. (6), who examined DNA methylation in 114 children in Gambia conceived in the rainy season (low-calorie, nutritionally rich) or dry season (high-calorie, nutritionally poor). Hierarchical clustering identified three groups of individuals, two of which showed the expected dichotomy of being imprinted or not imprinted. The third cluster, consisting of five children (4.4%) who had evidence of DNA methylation in only approximately half of the probes in the nc886 DMR, was classified as inconclusive (). When comparing imprinted individuals and nonimprinted individuals, we found that season of conception was significantly (P = 0.0417) associated with the frequency of imprinting of the nc886 DMR in infants (Fig. 4). In this data set, we performed likelihood ratio tests (LRTs) to adjust for available confounding factors, including exposure to aflatoxin (a common foodborne toxin in Africa) and sex. We found that neither variable was a significant predictor of imprinting at nc886, and, after adjusting for these factors, the significance of the season of conception slightly improves (P = 0.03902). Thus, season of conception and, by association, maternal nutrition, may contribute to “tuning” the likelihood that an individual has imprinting at the nc886 DMR. Furthermore, recent analysis of a Scottish birth cohort identified differential DNA methylation of nc886 in response to folate supplementation during pregnancy (53).

Fig. 4.

Maternal environment in utero impacts DNA methylation of nc886 in infants. HM450 data for the nc886 DMR was analyzed from three independent studies of peripheral blood of infants at birth. Each group is separated by the percentage of the population that is imprinted vs. not imprinted based on hierarchical clustering. Individuals with inconclusive evidence of imprinting have been removed from this analysis. Mat, maternal. *P < 0.05.

The same hierarchical clustering analysis was performed on data from Markunas et al. (54) for DNA methylation of the nc886 DMR in children born to 855 mothers of various ages (). We found that children born to mothers younger than 20 y had significantly less imprinting than those born to older mothers (P = 0.0167; Fig. 4). After adjustment for sex, maternal age is still significantly associated with imprinting of the nc886 DMR (LRT P = 0.0281). Our data suggest that teenage pregnancy is associated with a decreased likelihood that a child will have imprinting at the nc886 DMR.

We replotted data from van Dijk et al. (55) in which HM450 analysis was performed on blood spots obtained from children at birth. This study identified lack of imprinting at nc886 measured at birth as being strongly associated with increased body mass index (BMI) in children assessed at 5 y of age. Hierarchical clustering of the nc886 DMR was used to define imprinted, not-imprinted, and inconclusive groups, and BMI was analyzed as a continuous variable (). We confirmed that imprinting of the nc886 DMR was associated with lower BMI in children when analyzed as discrete quartiles (Fig. 4) and as a continuous variable (χ2 test P = 1.75 × 10−3 and logistic regression P = 1.94 × 10−3, respectively). We conclude that imprinting of the nc886 DMR is associated with a lower childhood BMI. Collectively, these three data sets demonstrate that maternal environment is linked to the establishment of imprinting across the nc886 DMR in children and that this epigenetic mark potentially impacts human phenotypes later in life.

Finally, we calculated the percentage of each population, across all samples, that exhibited imprinting (Fig. 4). When examining all samples from the BMI data set, we found that 76% of individuals were imprinted and 24% were not imprinted at the nc886 DMR, similar to the maternal age study and consistent with our previous findings in which 75% of individuals were imprinted and 25% were not (23). In the smallest data set, which examined maternal nutrition in Gambia, we found 82% of the population with imprinting. Whether this is a result of environmental differences or discrepancies in the genetic makeup of different populations needs to be investigated.

Discussion

We conclude that the nc886 DMR, when imprinted, is maternally methylated, given our evidence that this region behaves as a contiguous DMR that is not methylated in sperm (Fig. 1). The fact that the entire 1.9-kb nc886 DMR is subject to tunable polymorphic imprinting has been largely overlooked, emphasizing the need to scrutinize data obtained from the HM450 platform for possible dichotomization. We find that some of the polymorphic nature of the imprinting can be explained by local genetic makeup, as measured by using an A/C SNP in the centromeric CTCF site. The mechanism by which this polymorphism alters the likelihood of imprinting is not completely clear but is possibly associated with chromatin conformation variation triggered by CTCF binding, given reports that CTCF binding sites are involved in genomic imprinting (56–58). This possibility is supported by the observation that mutations in CTCF binding sites in the XIST promoter alter CTCF binding efficiency and choice of X chromosome inactivation (51). Additionally, the presence of a genetic influence is supported by the fact that monozygotic twins are found to be more concordant in DNA methylation at the nc886 DMR than dizygotic twins (35). However, Van Baak et al. (35) concluded that DNA methylation at nc886 could not be explained by genetics alone and classified it as a region of “epigenetic supersimilarity.” Although our data provide high-resolution analysis of local DNA methylation in 31 individuals, a very recent paper by Zink et al. (37) found no genetic polymorphisms associated with imprinting of nc886 in an Icelandic population, suggesting a need for further genotyping and concurrent DNA methylation analysis in individuals from different populations and environmental contexts.

We find that the frequency of imprinting of the nc886 DMR in children is also associated with the mothers’ age and season of conception, indicating parental environment as potentially another mechanism for tuning the likelihood that the nc886 DMR will be imprinted (6, 54). Furthermore, the frequency of imprinting of the nc886 DMR at birth is maintained and is directly associated with the BMI of children at the age of 5 y (55). Collectively, these results suggest that genetic and environmental factors may affect the establishment of imprinting of a DMR, which is closely associated with human physiology (Fig. 5). Unfortunately, there are no known SNPs in nc886 that would allow us to examine tissues for allele-specific expression. However, we have shown that DNA methylation can silence transcription of nc886, making it very likely that differential methylation of this region is functionally important (23, 34).

A further understanding of causality from all of these observations will require an unraveling of potential biological roles of nc886 or the nc886 DMR. nc886 has been variously described as a regulator of the dsRNA-dependent protein kinase R (PKR) (23, 59–61) and Dicer (62). It was also suggested as a possible tumor suppressor (63) or oncogene (64). Alternatively, some other genetic or epigenetic factors within or around the 1.9-kb DMR sequence may be drivers of phenotypic effects associated with nc886. Thus far, EWAS studies have primarily focused on variations in DNA methylation levels without considering the essential role of imprinting in human development as we have done here.

Methods

DNA Methylation Analysis of Illumina HM450 Data.

Normalized β-values were downloaded from Gene Expression Omnibus (GEO), with accession numbers GSE47627 and GES64096 for sperm and GSE55763 for genomic DNA from peripheral blood of 2,664 individuals assayed on the Infinium HM450 platform (i.e., LOLIPOP data). Probes falling in the region chr5:135413937–135419936 (Genome Build GRCh37) were extracted, and their signals were visualized by using the R package SeSAMe (10). Genomic features were extracted from the University of California, Santa Cruz, Genome Browser. Probes annotated with potential cross-hybridization issues and SNP bias were excluded according to an earlier study (10, 65). This analysis was also performed on matched normal tissues from the TCGA database, as shown in .

Locus-Specific Bisulfite Sequencing Analysis of the Centromeric CTCF Site and nc886.

Genomic DNA isolated from white blood cells of disease-free individuals was obtained from the Van Andel Research Institute Pathology and Biorepository Core. For trio analysis, the following DNA samples were obtained from the National Human Genome Research Institute Sample Repository for Human Genetic Research at the Coriell Institute for Medical Research: HG00731, HG00732, HG00733, HG01938, HG01939, and HG01940. Bisulfite conversion and clean-up of 2 μg of genomic DNA was performed by using the EZ DNA Methylation Kit (Zymo Research) according to the manufacturer’s protocol. Locus-specific PCR was performed by using primers specific for bisulfite-converted DNA and PCR products cloned by using the pGEM-T Easy vector and NEB 5-alpha–competent Escherichia coli (New England Biolabs). Colonies were screened for positive inserts by PCR and sequencing performed using the M13 promoter.

Statistical Analysis.

Determining whether SNP status is associated with DNA methylation in the centromeric CTCF region.

β-Binomial mixed-effects regression, weighted by the number of CpG sites in a region (10 for SNP A, 11 for C), was used to model data obtained from bisulfite sequencing and PCR via the R package (version 3.4.3; https://www.r-project.org/) glmmTMB (66). A random intercept for each individual was used to account for the paired heterozygous and homozygous measures; a second random intercept was included to account for DNA methylation of the nc886 region as determined by the presence of at least 10% of clones being completely methylated. β-Binomial was chosen because it is well-suited for modeling the binary methylated/unmethylated status of each CpG in the region, but the probability of a given site being methylated increases with each additional site methylated in the region. An LRT confirmed that weighted β-binomial was a better fit for these data than a weighted binomial, and Akaike information criterion and Bayesian information criterion showed this regression was better than negative binomial regression.

Determining the correlation between DNA methylation at the centromeric CTCF and nc886 regions.

The percentage of DNA methylation was calculated for each region analyzed by bisulfite conversion and sequencing, and the data were analyzed via logistic regression. Bisulfite sequencing was used to estimate DNA methylation in the centromeric CTCF and nc886 regions. Logistic regression was used to determine whether the probability of nc886 being methylated changed as the percentage of methylation of the CTCF region increased. The reported P value was calculated by using an LRT. The mean percent methylation of the centromeric CTCF site was estimated by using a mixed-effects negative binomial regression via the R package glmmTMB (66). Data were plotted as the mean change in probability per percent increase of methylation in the centromeric CTCF, with a 95% confidence band.

Population Analysis of DNA Methylation at nc886 from Published Datasets.

We downloaded from GEO Infinium HM450 data from genomic DNA of infants in studies of season of conception in Gambia (6) and maternal age at delivery (54) with accession numbers GSE59592 and GSE82273, respectively. For the childhood BMI study (55), the GEO Infinium HM450 data (accession no. GSE103657) was not optimal because of technical effects introduced by batch correction. Instead, the raw array data were renormalized, consistent with the methods of the published study, and BMI measurement data were accessed through the authors of the paper of van Dijk et al. (55).

Based on bisulfite sequencing analysis (), which provides single-base resolution of individual DNA strands for DNA methylation, and previous studies, we concluded that there are only two possibilities for DNA methylation in this region: near 0% DNA methylation (biallelically unmethylated) or near 50% DNA methylation (monoallelically methylated) (6, 23, 24, 67). However, when data were plotted, it was clear that not all individuals exhibited this expected dichotomy of 50% methylated or not methylated (). Upon further investigation, we found that individuals who did not fit this dichotomy comprised a third cluster based on hierarchical clustering via Manhattan distances. Therefore, we identified three clusters of individuals, which we called imprinted, not imprinted, or inconclusive based on partial methylation readouts throughout the region (). To remain consistent with previously observed data, we restricted our analyses to focus on the patients not in this third, inconsistent cluster. Results without this restriction are reported in the .

Treating imprinting as the response variable, we applied logistic regression by using the glm function in R. The season of conception, maternal age (dichotomized by teenage pregnancy), and childhood BMI Z-scores were used as the regressors of the three data sets, respectively. P values for the coefficients were obtained by using a Wald test.