Longitudinal, genome-scale analysis of DNA methylation in twins from birth to 18 months of age reveals rapid epigenetic change in early life and pair-specific effects of discordance
- David Martino1,
- Yuk Jin Loke2, 3,
- Lavinia Gordon4,
- Miina Ollikainen5,
- Mark N Cruickshank2, 6,
- Richard Saffery†1, 3 and
- Jeffrey M Craig†2, 3Email author
© Martino et al.; licensee BioMed Central Ltd. 2013
Received: 20 December 2012
Accepted: 22 May 2013
Published: 22 May 2013
The extent to which development- and age-associated epigenetic changes are influenced by genetic, environmental and stochastic factors remains to be discovered. Twins provide an ideal model with which to investigate these influences but previous cross-sectional twin studies provide contradictory evidence of within-pair epigenetic drift over time. Longitudinal twin studies can potentially address this discrepancy.
In a pilot, genome-scale study of DNA from buccal epithelium, a relatively homogeneous tissue, we show that one-third of the CpGs assayed show dynamic methylation between birth and 18 months. Although all classes of annotated genomic regions assessed show an increase in DNA methylation over time, probes located in intragenic regions, enhancers and low-density CpG promoters are significantly over-represented, while CpG islands and high-CpG density promoters are depleted among the most dynamic probes. Comparison of co-twins demonstrated that within-pair drift in DNA methylation in our cohort is specific to a subset of pairs, who show more differences at 18 months. The rest of the pairs show either minimal change in methylation discordance, or more similar, converging methylation profiles at 18 months. As with age-associated regions, sites that change in their level of within-pair discordance between birth and 18 months are enriched in genes involved in development, but the average magnitude of change is smaller than for longitudinal change.
Our findings suggest that DNA methylation in buccal epithelium is influenced by non-shared stochastic and environmental factors that could reflect a degree of epigenetic plasticity within an otherwise constrained developmental program.
Epigenetic modifications such as DNA methylation play an important role in development, ageing and disease [1–3]. However, the factors that influence epigenetic dynamics are poorly understood. Twin studies have the potential to estimate genetic components of epigenetic state [4, 5] and have demonstrated that gene expression and DNA methylation profiles can both be influenced by allelic, stochastic and environmental factors [6–10]. Non-shared environmental and stochastic factors together have been estimated to be the largest influence on promoter methylation in utero .
Studies of epigenetic change over time have predominantly used cross-sectional approaches and have focused on adults [11–17] or on intrauterine development [18, 19]. A small number of such studies have assessed age-associated DNA methylation across wider time-spans, encompassing childhood, adolescence and adulthood [20–23]. Consistently, age-associated changes in DNA methylation are more likely to involve (1) increases in methylation; (2) genes associated with development, signaling and regulation of transcription; and (3) regions involved in epigenetic reprogramming during embryonic stem cell differentiation [12, 14]. Since most of these studies have focused on CpG islands and promoter regions, age-associated epigenetic changes are incompletely characterized in relation to genomic coverage and life-course.
Longitudinal studies investigating aging and longevity have distinct advantages over cross-sectional designs, particularly in relation to controlling for genetic variation. For example, longitudinal studies directly query temporal sequences and pathways and individuals are studied rather than group averages. Longitudinal, array-based studies in blood from children in the first 1 to 5 years of postnatal life have shown similar results to adult studies with respect to gene function, genomic location and direction of age-related changes in DNA methylation, with the majority of age-related changes being observed in regions flanking CpG islands [24–26]. In adults, a high-resolution array-based study found a mixture of age-stable and age-dynamic variability throughout the methylome in adults . Other studies of global DNA methylation in adults also showed a genetic influence on increase and decrease in DNA methylation  and a decrease in interspersed repeat DNA methylation over time .
A small number of epigenetic studies of ageing have focused specifically on twins. Cross-sectional studies have found that older monozygotic (MZ) twins differ more with respect to global and repetitive DNA methylation , a phenomenon referred to as 'epigenetic drift'. In contrast, using DNA from saliva, no evidence for such drift was found within MZ twins aged 21 to 55 years using a promoter/CpG island array . A cross-sectional comparison of DNA methylation at the imprinted IGF2/H19 locus in adolescent and middle-aged MZ twins also found no evidence of epigenetic drift within pairs . To our knowledge, only one study has examined DNA methylation in twins longitudinally, measuring three loci in buccal DNA from 46 MZ pairs and 45 dizygotic (DZ) pairs at 5 and 10 years of age . This study revealed (1) locus-specific variability in DNA methylation; (2) change over time in individuals; (3) within-pair correlation and (4) epigenetic drift. Taken together, these studies show that epigenetic change over time is likely to be regulated by many factors, potentially in a tissue-specific and genome context-dependent manner. Longitudinal epigenetic studies in twins offer tremendous potential to further our understanding of the relationship between genetics and other factors that specify inter-individual temporal change in DNA methylation profile in humans.
We have used the Infinium HumanMethylation450 BeadChip (HM450) platform, which interrogates >485,000 CpG dinucleotides and contains probes from CpG islands, shores (2 kb regions flanking CpG islands), shelves (2 kb regions flanking shores), sites from 1,500 bp upstream of transcription start sites through to gene bodies and 3' UTRs, in addition to intergenic regions, regions involved in epigenetic reprogramming during embryonic stem cell differentiation and enhancers [32, 33]. Although repeats are not covered by these arrays and intergenic regions are not covered to the same depth as genic regions, the platform represents a significant step towards genome-scale coverage. Using the Infinium HM450 platform, we have performed a longitudinal study of DNA methylation at birth and age 18 months in DNA from buccal swabs from 10 MZ and 5 DZ twin pairs from the Peri/postnatal Epigenetic Twins Study (PETS) cohort . We report a large degree of epigenetic change during the first 18 months of postnatal life, with strong regional genomic biases for rate of change over time. We also present evidence for pair-specific levels of epigenetic change, suggesting a complex interplay between environment, non-shared environment and stochastic factors in molding the early postnatal epigenome.
Twin pair characteristics
Twin pair ID number
Twin 1 sex
Twin 2 sex
Birth weight discordance (%)c
Samples removed after QCd
T1_18, T2 _18
Determination of technical versus biological variation
Determining relationships between samples
Unsupervised hierarchical clustering of the entire dataset (Figure S1A in Additional file 1) revealed that most samples cluster according to age. The majority of co-twins also cluster together: 7/9 (78%) MZ co-twins cluster at birth and 6/6 (100%) at 18 months, while with DZ co-twins, 4/5 (80%) cluster at birth and 2/4 (50%) cluster at 18 months. To explore the variation in this dataset attributable to the effect of sequence variation on methylation values via cis genetic effects or probe hybridization, we performed hierarchical clustering selectively for probes overlapping known SNPs as defined by the HM450 SNP manifest (version 3, 103,148 probes). The results compared well to the full dataset: 7/9 MZ co-twins cluster at birth and 6/6 cluster at 18 months, while 4/5 DZ co-twins cluster at birth and 3/4 cluster at 18 months (Figure S1B in Additional file 1). Restricting this analysis to probes with reported SNPs at the CpG site assayed by the probe (2,527 probes in this data set) resulted in 8/9 MZ co-twins clustering at birth and 6/6 at 18 months (Figure S1C in Additional file 1). Interestingly, on average for this set of probes, DZ twins did not cluster with their co-twin; rather, DZ twins at birth clustered with their matched samples at 18 months. Thus, data for such probes are likely to reflect the genotype of the individual rather than representing purely methylation levels. A random sampling of the same number of SNP-associated probes did not reproduce this clustering (data not shown), indicating this effect did not represent a sampling bias. These results suggest that SNP-containing probes account for little variation in the overall data set, with the exception of probes with SNPs at the CpG site assayed.
Identification of age-associated differentially methylated probes
Ontology enrichment analysis of age-associated differentially methylated probes
Cell morphogenesis involved in differentiation
Regulation of small GTPase-mediated signal transduction
Regulation of Ras protein signal transduction
Neuron projection development
Cell morphogenesis involved in neuron development
Neuron projection morphogenesis
Identification of age-associated differentially methylated regions
Ontology enrichment analysis for age-associated differentially methylated regions
-log10 binomial P-value
Cellular development process
Response to organic substance
Regulation of transcription from RNAs pol III promoter
Regulatory region DNA binding
Transcription regulatory region DNA binding
Positive regulation of transcription from RNA pol III promoter
Negative regulation of transcription from RNA pol III promoter
Cardiac chamber development
Embryonic limb morphogenesis
Digestive tract morphogenesis
Epigenetic discordance within twin pairs at birth and 18 months
Within-pair epigenetic discordance resulting from non-shared environmental factors has been postulated to underscore variation in phenotypic traits [40, 41]. We examined discordance in DNA methylation profile within twin pairs at birth and at 18 months of age. We calculated twin discordance as the absolute difference in β methylation values within pairs in birth samples and separately for 18-month samples. We ranked all probes according to average within-pair discordance at each age and performed 'ranked-list' ontology, which differs from 'gene-list' ontology in that there is no requirement for a predefined cutoff. All probes on the array were ranked by their scores for average within-pair discordance at each age (most discordant to least discordant), and the ranked list of probes was analyzed by the GOrilla bioinformatics tool  to identify ontology terms over-represented at the top of the list, compared with the bottom. We found that the most discordant genes at birth were consistently enriched for ontology terms associated with RNA metabolism, including spliceosome components and transcription factors (Table S5 in Additional file 2). At 18 months of age, the most discordant genes were associated with a similar set of gene ontologies as seen at birth (Table S6 in Additional file 2). The genes with discordant probes at both time points include a wide array of spliceosome components (for example, WDR83 and CWC22), zinc finger proteins (for example, ZNF267, ZBTB1, ZNF10), ribosomal proteins (for example, RPS26, RPL15, RPL12) and transcription factors (for example, MAML1, HOXB13).
We next investigated the distribution of DNA methylation discordance across genomic regions to determine whether discordance is more likely to occur at specific genomic locations. We have shown previously, using HM27 arrays, that median within-pair methylation discordance increased with increasing distance from CpG islands in three tissues (cord blood mononuclear cells, human umbilical vein endothelial cells and placenta) in both MZ and DZ twins at birth . As the HM27 array focuses primarily on gene promoters and CpG islands, we repeated this analysis taking advantage of the diversity of genomic locations contained within the HM450 arrays. We calculated absolute within-pair discordance as before, and plotted probe discordance across genomic location at birth and at 18 months. The distribution of discordance values was consistent across all genomic annotations targeted on the array, with no evidence of regional enrichment (Figure S3A in Additional file 1). Similar results were observed selecting the top 10,000 most variable probes, or alternatively when the analysis was performed separately at birth and 18 months separately for both MZ and DZ twins (data not shown). We then filtered the dataset to include only probes present on the HM27 arrays and found evidence of higher levels of discordance around shores and shelves of CpG islands (Figure S3B in Additional file 1), which is consistent with our previously published observation with this platform .
Level of change in epigenetic discordance (drift versus convergence) over the first 18 months is a pair-specific phenomenon
Ontology enrichment analysis for drifting and converging pairs
Anatomical structure development
Cellular developmental process
Anatomical structure morphogenesis
Pattern specification process
Regulation of developmental process
Cell surface receptor signaling pathway
Regulation of multicellular organismal development
Regulation of signaling
Anatomical structure development
Regulation of developmental process
Regulation of cell differentiation
Cellular developmental process
Regulation of signaling
Regulation of multicellular organismal development
This study, examining DNA methylation profiles in buccal cells of young twins, has further confirmed the highly dynamic nature of the human epigenome postnatally, in agreement with previous studies in other tissues [20, 21, 24–26, 43]. Buccal cells represent a key bio-resource for age- and disease-associated epigenetic association studies [8, 12, 44, 45]. From a previous study  and our own unpublished data, this sample type comprises >90% squamous epithelial cells with <10% blood cells. Moreover, we minimized cell heterogeneity due to immune reactions by not collecting from infants with mouth infections. Almost a third of all HM450 probes in our final dataset showed significantly changed DNA methylation levels at FDR <0.05. Furthermore, the average absolute methylation change in these aDMPs was 4.7% (3.2% per year over 18 months) and almost 3% of aDMPs exhibited an absolute methylation change of >20%. These changes are similar in magnitude to those seen in blood from birth to one year of age using HM27 arrays and FDR <0.05 (3.9% of probes with changes >20%; average change of 9.2% per year)  and in T cells from birth to one year of age using HM450 arrays and FDR <0.01 (2.7% of probes with changes >20%; average change of 14% per year) . In addition, our findings are of similar magnitude to a cross-sectional study of DNA methylation in the prefrontal cortex from human cadavers using HM27 arrays . In combination with other cross-sectional studies [20, 22], these cumulative data support the idea that rate of change of DNA methylation in the genome in any one tissue is highest in utero, possibly reflecting extensive cellular differentiation during organogenesis, and then declines in childhood, with a further drop in adulthood [20, 21, 43].
In the current study, most (approximately 90%) age-associated changes involved an increase in DNA methylation over time. This agrees with previous longitudinal studies of early childhood [24–26] and cross-sectional studies of placenta throughout gestation , peripheral blood in children , peripheral blood in adults [12, 13, 16], buccal cells in adults  and in a comparison of differences in various tissues between fetal and adult tissues . A recent study comparing the entire methylomes of a newborn and a centenarian using bisulfite sequencing  observed an age-associated decrease in methylation in most genomic regions, including interspersed repetitive DNA, intergenic and intragenic regions. Although additional studies are needed to reconcile these discrepancies, it is also likely that age-associated methylation changes are dependent on genomic and tissue context, on the methylation analysis platform used and on sample size.
We found that aDMRs are more likely to be in intergenic and intragenic regions (Figure 3). Such regions were also enriched in aDMRs identified in multiple tissues in rats  and mice . The intergenic regions identified in our study are single copy and overlap with enhancers (Figure 3). Such regions undergo the most dynamic changes accompanying differentiation of pluripotent stem cells . Almost a third of aDMPs identified, and almost a half of aDMRs, lay within 5 kb of transcriptional start sites (Figure S2 in Additional file 1), implicating such regions in the regulation of gene expression. The higher proportion of aDMRs (44%) compared to aDMPs (29%) around gene promoters most likely reflects the higher CpG density and co-methylation (locally correlated methylation) within these regions . However, this could also be due at least in part to the relatively wider HM450 probe spacing in intergenic regions.
Ontology and pathways analysis showed that approximately three quarters of the genes associated with stringent aDMPs and all aDMRs are implicated in development and morphogenesis (Tables 2 and 3; Table S2 in Additional file 2). Signaling pathways, including those based on GTPase signal transduction, pathways intrinsic to development, were also enriched, as has been observed in other studies of age-associated changes in methylation [50, 53]. A bias towards genes involved in development has been seen in cross-sectional [13, 20, 54, 55] and longitudinal [25, 27] studies of human aging and in a study of embryonic stem cell differentiation .
It is interesting to note that approximately half of the top ten aDMP ontologies related to neural development (Table 2), a finding shared with previous studies of methylation in saliva  and blood cell fractions and buccal epithelium from adults  and in a large meta-analysis of multiple HM27 datasets from human brain and blood . Further studies are needed to ascertain whether these findings relate to biases related to the large proportion of the genome expressed in the brain.
We found that despite a trend towards increased methylation with age in all regions of the genome, CpG-dense promoter regions were particularly depleted in aDMRs and CpG-poor promoters were moderately increased (Figure 3), contrary to previous cross-sectional [13, 15, 16, 56] and longitudinal  studies of ageing-associated methylation change. This disagreement is most likely because HM27 arrays are enriched in CpG-dense CpG islands. However, our data agree with findings from studies using methods that include intergenic regions, that low CpG density promoters are enriched in mouse aDMRs , during differentiation of human embryonic stem cells  and between birth and very old age .
Epigenetic discordance at birth and 18 months of age
Very few genome-wide studies of methylation or expression have been performed on buccal cells. One such study, of buccal cells collected from 20 twins aged 13 to 14 years using a low resolution CpG island array, found no significant methylation differences within pairs . However, a study of smoking-induced differential gene expression in buccal cells identified a differentially expressed network of genes with, at the hub, transcription factors REL and CREB , which are among the top 10% most discordant genes at birth and 18 months in our data (Tables S5 and S6 in Additional file 2). Despite the extensive longitudinal changes in DNA methylation described above, we found that, in general, probes located within genes associated with RNA metabolism (for example, spliceosome components) and control of gene expression (for example, transcription factors) were consistently more discordant within twin pairs at both birth and 18 months of age. Of interest, this class of genes has previously been shown to have altered levels of transcription in buccal cancer .
Epigenetic drift and convergence
In the current study, we found that a summed value (Euclidean distance) of epigenetic discordance across hundreds of thousands of loci can vary between and within pairs and can increase or decrease over time. In accordance with our genome-scale findings, a longitudinal study of DNA methylation at seven imprinted gene loci in buccal cells between birth and one year of age showed that inter-individual variation similarly increased, decreased or remained similar in singletons and that the direction of change could differ between individuals . A longitudinal study of DNA methylation at three genes in buccal cells in 46 MZ and 45 DZ twin pairs found that methylation drifted in some pairs and converged in others over time . Similar results were found for MZ and DZ twins and a role for genetic, shared and non-shared environmental factors, dependent on genomic location, in these longitudinal changes was postulated . For MZ pairs, changes in within-pair discordance must be influenced solely by stochastic and non-shared environmental factors. Evidence for the latter comes from our previous studies of methylation in newborn twins [6, 7, 60] and from a cross-sectional study of DNA methylation in seven genes in whole blood from >200 MZ twin pairs aged 18 to 89 years . Data from a longitudinal, genome-scale study of DNA methylation (using HM450 arrays) in whole blood from an independent cohort of young adults (aged 22 to 32 years) also provides evidence of genome-scale methylation drift and convergence defined by changes in Euclidean distance over time (Figure S5 in Additional file 1).
Epigenetic drift has been postulated to arise from the cumulative effects of (non-shared) environment and stochastic events [30, 62, 63], the latter influenced by epigenetic events such as promoter occupancy by transcription factors  and by errors made during the maintenance of DNA methylation profile following DNA replication [30, 63]. Recent studies suggest that epigenetic drift may also reflect differing rates of change of methylation among the population . Furthermore, others have argued that epigenetic variability (or noise) is itself genetically programmed and has evolved to mediate some degree of plasticity (via canalization) . In contrast, we suggest that 'convergence' may involve sites of methylation equalization between co-twins, possibly reflecting regression to the mean as a contributing factor. Regression to the mean is a phenomenon in which it is a statistical certainty that individual phenotypes, such as growth patterns , shift to the population mean over time . This explains why twins with birth weight discordance become more similar over time  and can be understood in terms of twin-specific uterine-specific restrictions being replaced postnatally by a greater degree of shared environment [69–71]. Indeed, the twins in the current study had a median weight discordance [(Weight of the heavier twin - Weight of the lighter twin)/Weight of the heavier twin] of 13.3% at birth and 2.8% at 18 months. Although caution is needed with interpretations from a small sample size, we note that 'converging' pairs were more likely to start with a higher within-pair discordance (mean Euclidean distance = 375) than the drifting pairs (mean Euclidean distance = 295) (Figure 5b), although this difference did not reach significance (P = 0.11). Clearly, larger longitudinal twin-based studies are needed to further investigate factors contributing to epigenetic drifting and convergence over time.
We have conducted the first longitudinal study of epigenetic change in buccal cells in twins from birth, using a validated, genome-scale methylation array. We have shown evidence that the epigenetic profile of both MZ and DZ twin pairs can exhibit epigenetic drift or convergence early in postnatal development. As genes involved in development exhibited the largest absolute changes in methylation over time and the largest, smaller-scale changes within twin pairs, we conclude that the epigenetically driven developmental program is influenced to some extent by stochastic and/or non-shared environmental factors. Thus, canalization may be influenced by such factors, in addition to genetic factors as suggested by Waddington [72, 73].
Materials and methods
Subjects, tissues and DNA extraction
Sample collection from twins at the time of delivery was carried out with appropriate human ethics approval from the Royal Women's Hospital (project number 06/21), Mercy Hospital for Women (project number R06/30), and Monash Medical Centre (project number 06117C), Melbourne and the study was conducted according to the Declaration of Helsinki principles. The twin pairs chosen for methylation array analysis are shown in Table 1. The 10 MZ pairs and 5 DZ pairs shared a similar sex ratio, gestational age and birth weight to the full group of 250 pairs. Buccal cells were collected with Catch-all Sample Collection Swabs (EPICENTRE Biotechnologies, Madison, WI, USA) and were stored at -20°C until DNA extraction, which was performed as previously described .
Infinium HumanMethylation450 BeadChip data acquisition and processing
DNA samples (1 μg) were bisulfite converted using the Methyl EasyXceed bisulfite modification kit (Human Genetic Signatures, North Ryde, Australia), according to the manufacturer's instructions. Conversion efficiency was assessed by bisulfite-specific PCR. DNA samples were hybridized to Illumina Infinium Human Methylation450 (HM450) BeadChip arrays according to the manufacturer's instructions. Raw intensity data (IDAT) files were imported into the R environment (version 2.14.1)  and processed using the minfi package . All analyses were performed in R using packages available from the Bioconductor project . Data quality was assessed in minfi using plots derived from various control probes on the array. Poor performing probes defined as those with an average detection P-value >0.001 in one or more samples were removed from the analysis (n = 132,113). Data from five samples with an average detection P-value >0.05 and with evidence of poor bisulfite conversion efficiency were removed completely. Probes on the × and Y chromosomes were also discarded from all samples. The resulting data were pre-processed using the Illumina method within minfi and subset-quantile within-array normalization was performed  for combined normalization of Infinium type I and type II probes. The log2 ratio of methylated probe intensity to unmethylated probe intensity were calculated in minfi and the resulting M-values [77, 78] were quantile normalized between arrays using the limma package . Sample quality was further assessed using hierarchical clustering plots available in minfi and lumi  packages. Following this, three additional samples were removed as outliers constituting a final data set of 330,168 probes and 53 samples.
Exploratory analysis of sample relationships was performed using unsupervised hierarchical clustering analysis with the Euclidean distance and complete linkage algorithm, and dendrogram was created using gplots . Differential methylation analysis was performed on M-values using the limma package using a cutoff of FDR-corrected P-values <0.05  and delta beta values >0.2. To study discordance among co-twins at the probe-level, a linear model was fitted to the M-values with twin-pair as a predictive factor to model the twin relationship. The level of discordance among co-twins was interpreted as the residual measurement for each CpG from the model-fit. For enrichment analysis, gene sets were populated with probe IDs using the annotated regions provided in the Illumina HM450 manifest file (version 1.1). Annotations used were classified as gene-related (TSS1500 and TSS200, regions from -1500 to -200 and -200 to the transcriptional start site respectively, 5' UTRs, first exons, gene bodies, 3' UTR and intergenic (no gene annotation)); CpG island-related (islands (also split into intragenic and intergenic)), shores (0 to 2 kb flanking islands), shelves (2 to 4 kb flanking islands) and open sea (>4 kb from islands) ); DMRs (associated with cancer (CDMRs) and induced pluripotent stem cell reprogramming (RDMRs);  and regulatory regions (promoters, enhancers and DNAse hypersensitivity sites, likely to be a mixture of promoters and enhancers [84, 85]). Boxplots were produced to graph each category by discordance score. The 'bump-hunting' methods described by Jaffe and colleagues  were implemented using the charm package available in Bioconductor . We used the 'dmrFinder' algorithm without covariate adjustment, using the default SPAN settings and specifying a minimum four probes, and a t-statistic cutoff to identify probes as being in a DMR at 0.995. For gene ontologies the GOrilla bioinformatics tool  was used to perform ranked-list ontology using the entire array content ranked by scores for discordance. Gene-list ontology enrichment was performed on significant gene lists (FDR <0.05) using the DAVID bioinformatics tool under the default settings . Pathway analysis data were analyzed through the use of Ingenuity Pathway Analysis (Ingenuity Systems, Redwood City, CA, USA). The analysis tool GREAT (Genomic Regions of Annotations Tool)  was used to analyze the functional significance of aDMRs using the single nearest gene association rule within a 100 kb window.
Sequenom MassArray target validation
Target validation was performed using the Sequenom MassArray EpiTYPER (Sequenom, San Diego, CA, USA) performed as previously described [18, 60]. Amplicons were designed using Sequenom EpiDesigner software. Primers are listed in Table S7 in Additional file 2. In brief, amplification was performed after bisulfite conversion of genomic DNA with the MethylEasy Xceed bisulphite conversion kit (Human Genetic Signatures, North Ryde, Australia). All PCR amplifications and downstream processing were carried out at least in duplicate and the mean methylation level at specific CpG sites determined. Raw data obtained from MassArray EpiTYPING were cleaned systematically using an R-script to remove samples that failed to generate data for more than 70% of CpG sites tested . Also, technical replicates showing ≥10% absolute difference from the median value of the technical replicates were removed and only samples with at least two successful technical replicates were analyzed. Samples were compared across each analyzable CpG site in the amplicon, as well as the mean across the whole amplicon.
Array data described in this manuscript have been submitted to the Gene Expression Omnibus public repository and are freely available under the accession number GSE42700.
age-associated differentially methylated probe/region
differentially methylated region
false discovery rate
Infinium HumanMethylation 450 BeadChip array platform
reprogramming differentially methylated probe/region
single nucleotide polymorphism
We wish to thank Ruth Morley, John Carlin, Clinical Epidemiology and Biostatistics Unit, MCRI, Mark Umstad, Royal Women's Hospital, Melbourne, Euan Wallace, Monash Medical Centre, Melbourne and Mark Permezel, Mercy Hospital For Women, Melbourne for their contributions to establishing the PETS cohort; Sarah Healy, Tina Vaiano, Nicole Brooks, Jennifer Foord, Sheila Holland, Anne Krastev, Siva Illancheran and Joanne Mockler for recruitment and sample collection; Research Assistant Xin Li, Technical officer Anna Czajko, Study Coordinator Geraldine McIlroy, and all mothers and twins that participated in this study. We also thank Jaakko Kaprio, Khadeeja Ismail and Kirsi Pietiläinen for sharing unpublished data and Ben Ong for help with Sequenom analysis. This work was supported by grants from the from the Australian National Health and Medical Research Council (grant numbers 437015 and 607358 to JMC and RS), the Financial Markets Foundation for Children (grant number 032-2007) and by the Victorian Government's Operational Infrastructure Support Program. RS is supported by a NHMRC Senior Research Fellowship. JMC and DM would also like to acknowledge financial support from the Murdoch Childrens Research Institute.
- Feinberg AP: Phenotypic plasticity and the epigenetics of human disease. Nature. 2007, 447: 433-440. 10.1038/nature05919.PubMedView ArticleGoogle Scholar
- Foley DL, Craig JM, Morley R, Olsson CA, Dwyer T, Smith K, Saffery R: Prospects for epigenetic epidemiology. Am J Epidemiol. 2009, 169: 389-400.PubMedPubMed CentralView ArticleGoogle Scholar
- Petronis A: Epigenetics as a unifying principle in the aetiology of complex traits and diseases. Nature. 2010, 465: 721-727. 10.1038/nature09230.PubMedView ArticleGoogle Scholar
- Bell JT, Spector TD: A twin approach to unraveling epigenetics. Trends Genet. 2011, 27: 116-125. 10.1016/j.tig.2010.12.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Bell JT, Saffery R: The value of twins in epigenetic epidemiology. Int J Epidemiol. 2012, 41: 140-150. 10.1093/ije/dyr179.PubMedView ArticleGoogle Scholar
- Gordon L, Joo EJ, Andronikos R, Ollikainen M, Wallace EM, Umstad MP, Permezel M, Oshlack A, Morley R, Carlin JB, Saffery R, Smyth GK, Craig JM: Expression discordance of monozygotic twins at birth: effect of intrauterine environment and a possible mechanism for fetal programming. Epigenetics. 2011, 6: 579-592. 10.4161/epi.6.5.15072.PubMedView ArticleGoogle Scholar
- Gordon L, Joo EJ, Powell JE, Ollikainen M, Novakovic B, Li X, Andronikos R, Cruickshank MN, Conneely KN, Smith AK, Alisch RS, Morley R, Visscher PM, Craig JM, Saffery R: Neonatal DNA methylation profile in human twins is specified by a complex interplay between intrauterine environmental/genetic factors subject to tissue-specific influence. Genome Res. 2012, 22: 1395-1406. 10.1101/gr.136598.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Wong CC, Caspi A, Williams B, Craig IW, Houts R, Ambler A, Moffitt TE, Mill J: A longitudinal study of epigenetic variation in twins. Epigenetics. 2010, 5: 516-526. 10.4161/epi.5.6.12226.PubMedPubMed CentralView ArticleGoogle Scholar
- Dempster EL, Pidsley R, Schalkwyk LC, Owens S, Georgiades A, Kane F, Kalidindi S, Picchioni M, Kravariti E, Toulopoulou T, Murray RM, Mill J: Disease-associated epigenetic changes in monozygotic twins discordant for schizophrenia and bipolar disorder. Hum Mol Genet. 2011, 20: 4786-4796. 10.1093/hmg/ddr416.PubMedPubMed CentralView ArticleGoogle Scholar
- Gervin K, Hammero M, Akselsen HE, Moe R, Nygard H, Brandt I, Gjessing HK, Harris JR, Undlien DE, Lyle R: Extensive variation and low heritability of DNA methylation identified in a twin study. Genome Res. 2011, 21: 1813-1821. 10.1101/gr.119685.110.PubMedPubMed CentralView ArticleGoogle Scholar
- Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, Sommer IE, Kahn RS, Ophoff RA: The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS One. 2009, 4: e6767-10.1371/journal.pone.0006767.PubMedPubMed CentralView ArticleGoogle Scholar
- Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, Beyan H, Whittaker P, McCann OT, Finer S, Valdes AM, Leslie RD, Deloukas P, Spector TD: Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010, 20: 434-439. 10.1101/gr.103101.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Bell JT, Tsai PC, Yang TP, Pidsley R, Nisbet J, Glass D, Mangino M, Zhai G, Zhang F, Valdes A, Shin SY, Dempster EL, Murray RM, Grundberg E, Hedman AK, Nica A, Small KS, Dermitzakis ET, McCarthy MI, Mill J, Spector TD, Deloukas P: Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population. PLoS Genet. 2012, 8: e1002629-10.1371/journal.pgen.1002629.PubMedPubMed CentralView ArticleGoogle Scholar
- Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, Shen H, Campan M, Noushmehr H, Bell CG, Maxwell AP, Savage DA, Mueller-Holzner E, Marth C, Kocjan G, Gayther SA, Jones A, Beck S, Wagner W, Laird PW, Jacobs IJ, Widschwendter M: Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010, 20: 440-446. 10.1101/gr.103606.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Bocklandt S, Lin W, Sehl ME, Sanchez FJ, Sinsheimer JS, Horvath S, Vilain E: Epigenetic predictor of age. PLoS One. 2011, 6: e14821-10.1371/journal.pone.0014821.PubMedPubMed CentralView ArticleGoogle Scholar
- Christensen BC, Houseman EA, Marsit CJ, Zheng S, Wrensch MR, Wiemels JL, Nelson HH, Karagas MR, Padbury JF, Bueno R, Sugarbaker DJ, Yeh RF, Wiencke JK, Kelsey KT: Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context. PLoS Genet. 2009, 5: e1000602-10.1371/journal.pgen.1000602.PubMedPubMed CentralView ArticleGoogle Scholar
- Hernandez MI, Mericq V: Metabolic syndrome in children born small-for-gestational age. Arquivos Brasileiros Endocrinologia Metabologia. 2011, 55: 583-589. 10.1590/S0004-27302011000800012.View ArticleGoogle Scholar
- Novakovic B, Gordon L, Wong NC, Moffett A, Manuelpillai U, Craig JM, Sharkey A, Saffery R: Wide ranging DNA methylation differences of primary trophoblast cell populations and derived-cell lines: implications and opportunities for understanding trophoblast function. Mol Hum Reprod. 2011, 17: 344-353. 10.1093/molehr/gar005.PubMedPubMed CentralView ArticleGoogle Scholar
- Schroeder JW, Conneely KN, Cubells JC, Kilaru V, Newport DJ, Knight BT, Stowe ZN, Brennan PA, Krushkal J, Tylavsky FA, Taylor RN, Adkins RM, Smith AK: Neonatal DNA methylation patterns associate with gestational age. Epigenetics. 2011, 6: 1498-1504. 10.4161/epi.6.12.18296.PubMedPubMed CentralView ArticleGoogle Scholar
- Alisch RS, Barwick BG, Chopra P, Myrick LK, Satten GA, Conneely KN, Warren ST: Age-associated DNA methylation in pediatric populations. Genome Res. 2012, 22: 623-632. 10.1101/gr.125187.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Siegmund KD, Connor CM, Campan M, Long TI, Weisenberger DJ, Biniszkiewicz D, Jaenisch R, Laird PW, Akbarian S: DNA methylation in the human cerebral cortex is dynamically regulated throughout the life span and involves differentiated neurons. PLoS One. 2007, 2: e895-10.1371/journal.pone.0000895.PubMedPubMed CentralView ArticleGoogle Scholar
- Yuen RK, Neumann SM, Fok AK, Penaherrera MS, McFadden DE, Robinson WP, Kobor MS: Extensive epigenetic reprogramming in human somatic tissues between fetus and adult. Epigenet Chromatin. 2011, 4: 7-10.1186/1756-8935-4-7.View ArticleGoogle Scholar
- Heyn H, Li N, Ferreira HJ, Moran S, Pisano DG, Gomez A, Diez J, Sanchez-Mut JV, Setien F, Carmona FJ, Puca AA, Sayols S, Pujana MA, Serra-Musach J, Iglesias-Platas I, Formiga F, Fernandez AF, Fraga MF, Heath SC, Valencia A, Gut IG, Wang J, Esteller M: Distinct DNA methylomes of newborns and centenarians. Proc Natl Acad Sci USA. 2012, 109: 10522-10527. 10.1073/pnas.1120658109.PubMedPubMed CentralView ArticleGoogle Scholar
- Martino DJ, Tulic MK, Gordon L, Hodder M, Richman T, Metcalfe J, Prescott SL, Saffery R: Evidence for age-related and individual-specific changes in DNA methylation profile of mononuclear cells during early immune development in humans. Epigenetics. 2011, 6: 1085-1094. 10.4161/epi.6.9.16401.PubMedView ArticleGoogle Scholar
- Martino D, Maksimovic J, Joo JH, Prescott SL, Saffery R: Genome-scale profiling reveals a subset of genes regulated by DNA methylation that program somatic T-cell phenotypes in humans. Genes Immun. 2012, 13: 388-398. 10.1038/gene.2012.7.PubMedView ArticleGoogle Scholar
- Wang D, Liu X, Zhou Y, Xie H, Hong X, Tsai HJ, Wang G, Liu R, Wang X: Individual variation and longitudinal pattern of genome-wide DNA methylation from birth to the first two years of life. Epigenetics. 2012, 7: 594-605. 10.4161/epi.20117.PubMedPubMed CentralView ArticleGoogle Scholar
- Feinberg AP, Irizarry RA, Fradin D, Aryee MJ, Murakami P, Aspelund T, Eiriksdottir G, Harris TB, Launer L, Gudnason V, Fallin MD: Personalized epigenomic signatures that are stable over time and covary with body mass index. Sci Transl Med. 2010, 2: 49ra67-10.1126/scitranslmed.3001262.PubMedPubMed CentralGoogle Scholar
- Bjornsson HT, Sigurdsson MI, Fallin MD, Irizarry RA, Aspelund T, Cui H, Yu W, Rongione MA, Ekstrom TJ, Harris TB, Launer LJ, Eiriksdottir G, Leppert MF, Sapienza C, Gudnason V, Feinberg AP: Intra-individual change over time in DNA methylation with familial clustering. JAMA. 2008, 299: 2877-2883. 10.1001/jama.299.24.2877.PubMedPubMed CentralView ArticleGoogle Scholar
- Bollati V, Schwartz J, Wright R, Litonjua A, Tarantini L, Suh H, Sparrow D, Vokonas P, Baccarelli A: Decline in genomic DNA methylation through aging in a cohort of elderly subjects. Mech Ageing Dev. 2009, 130: 234-239. 10.1016/j.mad.2008.12.003.PubMedPubMed CentralView ArticleGoogle Scholar
- Fraga MF, Ballestar E, Paz MF, Ropero S, Setien F, Ballestar ML, Heine-Suner D, Cigudosa JC, Urioste M, Benitez J, Boix-Chornet M, Sanchez-Aguilera A, Ling C, Carlsson E, Poulsen P, Vaag A, Stephan Z, Spector TD, Wu YZ, Plass C, Esteller M: Epigenetic differences arise during the lifetime of monozygotic twins. Proc Natl Acad Sci USA. 2005, 102: 10604-10609. 10.1073/pnas.0500398102.PubMedPubMed CentralView ArticleGoogle Scholar
- Heijmans BT, Kremer D, Tobi EW, Boomsma DI, Slagboom PE: Heritable rather than age-related environmental and stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus. Hum Mol Genet. 2007, 16: 547-554. 10.1093/hmg/ddm010.PubMedView ArticleGoogle Scholar
- Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, Fan JB, Shen R: High density DNA methylation array with single CpG site resolution. Genomics. 2011, 98: 288-295. 10.1016/j.ygeno.2011.07.007.PubMedView ArticleGoogle Scholar
- Sandoval J, Heyn H, Moran S, Serra-Musach J, Pujana MA, Bibikova M, Esteller M: Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011, 6: 692-702. 10.4161/epi.6.6.16196.PubMedView ArticleGoogle Scholar
- Saffery R, Morley R, Carlin JB, Joo JH, Ollikainen M, Novakovic B, Andronikos R, Li X, Loke YJ, Carson N, Wallace EM, Umstad MP, Permezel M, Galati JC, Craig JM: Cohort Profile: The Peri/post-natal Epigenetic Twins Study. Int J Epidemiol. 2012, 41: 55-61. 10.1093/ije/dyr140.PubMedView ArticleGoogle Scholar
- Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F: Evaluation of the Infinium Methylation 450K technology. Epigenomics. 2011, 3: 771-784. 10.2217/epi.11.105.PubMedView ArticleGoogle Scholar
- Maksimovic J, Gordon L, Oshlack A: SWAN: Subset quantile Within-Array Normalization for Illumina Infinium HumanMethylation450 BeadChips. Genome Biol. 2012, 13: R44-10.1186/gb-2012-13-6-r44.PubMedPubMed CentralView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- Bibikova M, Le J, Barnes B, Saedinia-Melnyk S, Zhou L, Shen R, Gunderson KL: Genome-wide methylation profiling using Infinium assay. Epigenomics. 2009, 1: 177-200. 10.2217/epi.09.14.PubMedView ArticleGoogle Scholar
- Jaffe AE, Murakami P, Lee H, Leek JT, Fallin MD, Feinberg AP, Irizarry RA: Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int J Epidemiol. 2012, 41: 200-209. 10.1093/ije/dyr238.PubMedPubMed CentralView ArticleGoogle Scholar
- Salsberry PJ, Reagan PB: Effects of heritability, shared environment, and nonshared intrauterine conditions on child and adolescent BMI. Obesity (Silver Spring). 2010, 18: 1775-1780. 10.1038/oby.2009.485.View ArticleGoogle Scholar
- Plomin R: Commentary: Why are children in the same family so different? Non-shared environment three decades later. Int J Epidemiol. 2011, 40: 582-592. 10.1093/ije/dyq144.PubMedPubMed CentralView ArticleGoogle Scholar
- Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009, 10: 48-10.1186/1471-2105-10-48.PubMedPubMed CentralView ArticleGoogle Scholar
- Schneider E, Pliushch G, El Hajj N, Galetzka D, Puhl A, Schorsch M, Frauenknecht K, Riepert T, Tresch A, Muller AM, Coerdt W, Zechner U, Haaf T: Spatial, temporal and interindividual epigenetic variation of functionally important DNA methylation patterns. Nucleic Acids Res. 2010, 38: 3880-3890. 10.1093/nar/gkq126.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaminsky ZA, Tang T, Wang SC, Ptak C, Oh GH, Wong AH, Feldcamp LA, Virtanen C, Halfvarson J, Tysk C, McRae AF, Visscher PM, Montgomery GW, Martin NG, Petronis A: DNA methylation profiles in monozygotic and dizygotic twins. Nat Genet. 2009, 41: 240-245. 10.1038/ng.286.PubMedView ArticleGoogle Scholar
- Murphy SK, Huang Z, Hoyo C: Differentially methylated regions of imprinted genes in prenatal, perinatal and postnatal human tissues. PloS One. 2012, 7: e40924-10.1371/journal.pone.0040924.PubMedPubMed CentralView ArticleGoogle Scholar
- Thiede C, Prange-Krex G, Freiberg-Richter J, Bornhauser M, Ehninger G: Buccal swabs but not mouthwash samples can be used to obtain pretransplant DNA fingerprints from recipients of allogeneic bone marrow transplants. Bone Marrow Transplant. 2000, 25: 575-577.PubMedView ArticleGoogle Scholar
- Numata S, Ye T, Hyde TM, Guitart-Navarro X, Tao R, Wininger M, Colantuoni C, Weinberger DR, Kleinman JE, Lipska BK: DNA methylation signatures in development and aging of the human prefrontal cortex. Am J Hum Genet. 2012, 90: 260-272. 10.1016/j.ajhg.2011.12.020.PubMedPubMed CentralView ArticleGoogle Scholar
- Novakovic B, Yuen RK, Gordon L, Penaherrera MS, Sharkey A, Moffett A, Craig JM, Robinson WP, Saffery R: Evidence for widespread changes in promoter methylation profile in human placenta in response to increasing gestational age and environmental/stochastic factors. BMC Genomics. 2011, 12: 529-10.1186/1471-2164-12-529.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson RF, Atzmon G, Gheorghe C, Liang HQ, Lowes C, Greally JM, Barzilai N: Tissue-specific dysregulation of DNA methylation in aging. Aging Cell. 2010, 9: 506-518. 10.1111/j.1474-9726.2010.00577.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Liang P, Song F, Ghosh S, Morien E, Qin M, Mahmood S, Fujiwara K, Igarashi J, Nagase H, Held WA: Genome-wide survey reveals dynamic widespread tissue-specific changes in DNA methylation during development. BMC Genomics. 2011, 12: 231-10.1186/1471-2164-12-231.PubMedPubMed CentralView ArticleGoogle Scholar
- Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, Zhang X, Bernstein BE, Nusbaum C, Jaffe DB, Gnirke A, Jaenisch R, Lander ES: Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008, 454: 766-770.PubMedPubMed CentralGoogle Scholar
- Li Y, Zhu J, Tian G, Li N, Li Q, Ye M, Zheng H, Yu J, Wu H, Sun J, Zhang H, Chen Q, Luo R, Chen M, He Y, Jin X, Zhang Q, Yu C, Zhou G, Huang Y, Cao H, Zhou X, Guo S, Hu X, Li X, Kristiansen K, Bolund L, Xu J, Wang W, Yang H, et al: The DNA methylome of human peripheral blood mononuclear cells. PLoS Biol. 2010, 8: e1000533-10.1371/journal.pbio.1000533.PubMedPubMed CentralView ArticleGoogle Scholar
- Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, Low HM, Kin Sung KW, Rigoutsos I, Loring J, Wei CL: Dynamic changes in the human methylome during differentiation. Genome Res. 2010, 20: 320-331. 10.1101/gr.101907.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Gentilini D, Mari D, Castaldi D, Remondini D, Ogliari G, Ostan R, Bucci L, Sirchia SM, Tabano S, Cavagnini F, Monti D, Franceschi C, Di Blasio AM, Vitale G: Role of epigenetics in human aging and longevity: genome-wide DNA methylation profile in centenarians and centenarians' offspring. Age (Dordr). 2012Google Scholar
- Horvath S, Zhang Y, Langfelder P, Kahn RS, Boks MP, van Eijk K, van den Berg LH, Ophoff RA: Aging effects on DNA methylation modules in human brain and blood tissue. Genome Biol. 2012, 13: R97-10.1186/gb-2012-13-10-r97.PubMedPubMed CentralView ArticleGoogle Scholar
- Hernandez DG, Nalls MA, Gibbs JR, Arepalli S, van der Brug M, Chong S, Moore M, Longo DL, Cookson MR, Traynor BJ, Singleton AB: Distinct DNA methylation changes highly correlated with chronological age in the human brain. Hum Mol Genet. 2011, 20: 1164-1172. 10.1093/hmg/ddq561.PubMedPubMed CentralView ArticleGoogle Scholar
- Brunner AL, Johnson DS, Kim SW, Valouev A, Reddy TE, Neff NF, Anton E, Medina C, Nguyen L, Chiao E, Oyolu CB, Schroth GP, Absher DM, Baker JC, Myers RM: Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver. Genome Res. 2009, 19: 1044-1056. 10.1101/gr.088773.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Kupfer DM, White VL, Jenkins MC, Burian D: Examining smoking-induced differential gene expression changes in buccal mucosa. BMC Med Genomics. 2010, 3: 24-10.1186/1755-8794-3-24.PubMedPubMed CentralView ArticleGoogle Scholar
- Sajnani MR, Patel AK, Bhatt VD, Tripathi AK, Ahir VB, Shankar V, Shah S, Shah TM, Koringa PG, Jakhesara SJ, Joshi CG: Identification of novel transcripts deregulated in buccal cancer by RNA-seq. Gene. 2012, 507: 152-158. 10.1016/j.gene.2012.07.036.PubMedView ArticleGoogle Scholar
- Ollikainen M, Smith KR, Joo EJ, Ng HK, Andronikos R, Novakovic B, Abdul Aziz NK, Carlin JB, Morley R, Saffery R, Craig JM: DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome. Hum Mol Genet. 2010, 19: 4176-4188. 10.1093/hmg/ddq336.PubMedView ArticleGoogle Scholar
- Talens RP, Christensen K, Putter H, Willemsen G, Christiansen L, Kremer D, Suchiman HE, Slagboom PE, Boomsma DI, Heijmans BT: Epigenetic variation during the adult lifespan: cross-sectional and longitudinal data on monozygotic twin pairs. Aging Cell. 2012, 11: 694-703. 10.1111/j.1474-9726.2012.00835.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin GM: Epigenetic drift in aging identical twins. Proc Natl Acad Sci USA. 2005, 102: 10413-10414. 10.1073/pnas.0504743102.PubMedPubMed CentralView ArticleGoogle Scholar
- Petronis A: Epigenetics and twins: three variations on the theme. Trends Genet. 2006, 22: 347-350. 10.1016/j.tig.2006.04.010.PubMedView ArticleGoogle Scholar
- Whitelaw NC, Chong S, Whitelaw E: Tuning in to noise: epigenetics and intangible variation. Dev Cell. 2010, 19: 649-650. 10.1016/j.devcel.2010.11.001.PubMedView ArticleGoogle Scholar
- Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, Klotzle B, Bibikova M, Fan JB, Gao Y, Deconde R, Chen M, Rajapakse I, Friend S, Ideker T, Zhang K: Genome-wide Methylation Profiles Reveal Quantitative Views of Human Aging Rates. Mol Cell. 2013, 49: 359-367. 10.1016/j.molcel.2012.10.016.PubMedPubMed CentralView ArticleGoogle Scholar
- Pujadas E, Feinberg AP: Regulated noise in the epigenetic landscape of development and disease. Cell. 2012, 148: 1123-1131. 10.1016/j.cell.2012.02.045.PubMedPubMed CentralView ArticleGoogle Scholar
- Touwslager RN, Gerver WJ, Mulder AL, Jansen AJ, de Bruin R: Longitudinal growth during the first years of life: what is normal?. Hormone Res. 2008, 70: 273-277. 10.1159/000157873.PubMedView ArticleGoogle Scholar
- Barnett AG, van der Pols JC, Dobson AJ: Regression to the mean: what it is and how to deal with it. Int J Epidemiol. 2005, 34: 215-220.PubMedView ArticleGoogle Scholar
- Wilson RS: Twin growth: initial deficit, recovery, and trends in concordance from birth to nine years. Ann Hum Biol. 1979, 6: 205-220. 10.1080/03014467900007212.PubMedView ArticleGoogle Scholar
- Keith L, Machin G: Zygosity testing. Current status and evolving issues. J Reprod Med. 1997, 42: 699-707.PubMedGoogle Scholar
- Martin N, Boomsma D, Machin G: A twin-pronged attack on complex traits. Nat Genet. 1997, 17: 387-392. 10.1038/ng1297-387.PubMedView ArticleGoogle Scholar
- Waddington CH: Canalization of development and genetic assimilation of acquired characters. Nature. 1959, 183: 1654-1655. 10.1038/1831654a0.PubMedView ArticleGoogle Scholar
- Waddington CH: The Strategy of the Genes. 1957, London: Geo Allen and UnwinGoogle Scholar
- R_Development_Core_Team: R: A Language and Environment for Statistical Computing. 2009, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Hansen DK, Aryee MJ: Analyze Illumina's 450k methylation arrays. R package version 1.6.0. [http://bioconductor.org/packages/release/bioc/html/minfi.html]
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView ArticleGoogle Scholar
- Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224.PubMedView ArticleGoogle Scholar
- Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, Lin SM: Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2011, 11: 587-View ArticleGoogle Scholar
- Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconducto. Edited by: Gentleman R, Carey V, Dudoit S, Huber W, Irizarry R. 2005, New York: Springer, 397-420.View ArticleGoogle Scholar
- Docherty SJ, Davis OS, Haworth CM, Plomin R, Mill J: DNA methylation profiling using bisulfite-based epityping of pooled genomic DNA. Methods. 2010, 52: 255-258. 10.1016/j.ymeth.2010.06.017.PubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995, 57: 289-300.Google Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.PubMedPubMed CentralView ArticleGoogle Scholar
- Doi A, Park IH, Wen B, Murakami P, Aryee MJ, Irizarry R, Herb B, Ladd-Acosta C, Rho J, Loewer S, Miller J, Schlaeger T, Daley GQ, Feinberg AP: Differential methylation of tissue- and cancer-specific CpG island shores distinguishes human induced pluripotent stem cells, embryonic stem cells and fibroblasts. Nat Genet. 2009, 41: 1350-1353. 10.1038/ng.471.PubMedPubMed CentralView ArticleGoogle Scholar
- Lian H, Thompson WA, Thurman R, Stamatoyannopoulos JA, Noble WS, Lawrence CE: Automated mapping of large-scale chromatin structure in ENCODE. Bioinformatics. 2008, 24: 1911-1916. 10.1093/bioinformatics/btn335.PubMedPubMed CentralView ArticleGoogle Scholar
- Xi H, Shulha HP, Lin JM, Vales TR, Fu Y, Bodine DM, McKay RD, Chenoweth JG, Tesar PJ, Furey TS, Ren B, Weng Z, Crawford GE: Identification and characterization of cell type-specific and ubiquitous chromatin regulatory structures in the human genome. PLoS Genet. 2007, 3: e136-10.1371/journal.pgen.0030136.PubMedPubMed CentralView ArticleGoogle Scholar
- Aryee MJ, Wu Z, Ladd-Acosta C, Herb B, Feinberg AP, Yegnasubramanian S, Irizarry RA: Accurate genome-scale percentage DNA methylation estimates from microarray data. Biostatistics. 2011, 12: 197-210. 10.1093/biostatistics/kxq055.PubMedPubMed CentralView ArticleGoogle Scholar
- Barton G, Abbott J, Chiba N, Huang DW, Huang Y, Krznaric M, Mack-Smith J, Saleem A, Sherman BT, Tiwari B, Tomlinson C, Aitman T, Darlington J, Game L, Sternberg MJ, Butcher SA: EMAAS: an extensible grid-based rich internet application for microarray data analysis and management. BMC Bioinformatics. 2008, 9: 493-10.1186/1471-2105-9-493.PubMedPubMed CentralView ArticleGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G: GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010, 28: 495-501. 10.1038/nbt.1630.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.