BackgroundInterindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.ResultsFirst, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.ConclusionsThe non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-015-0660-y) contains supplementary material, which is available to authorized users.

Background DNA methylation is thought to be an important determinant of human phenotypic variation, but its inherent cell type specificity has impeded progress on this question. At exceptional genomic regions, interindividual variation in DNA methylation occurs systemically. Like genetic variants, systemic interindividual epigenetic variants are stable, can influence phenotype, and can be assessed in any easily biopsiable DNA sample. We describe an unbiased screen for human genomic regions at which interindividual variation in DNA methylation is not tissue-specific. Results For each of 10 donors from the NIH Genotype-Tissue Expression (GTEx) program, CpG methylation is measured by deep whole-genome bisulfite sequencing of genomic DNA from tissues representing the three germ layer lineages: thyroid (endoderm), heart (mesoderm), and brain (ectoderm). We develop a computational algorithm to identify genomic regions at which interindividual variation in DNA methylation is consistent across all three lineages. This approach identifies 9926 correlated regions of systemic interindividual variation (CoRSIVs). These regions, comprising just 0.1% of the human genome, are inter-correlated over long genomic distances, associated with transposable elements and subtelomeric regions, conserved across diverse human ethnic groups, sensitive to periconceptional environment, and associated with genes implicated in a broad range of human disorders and phenotypes. CoRSIV methylation in one tissue can predict expression of associated genes in other tissues. Conclusions In addition to charting a previously unexplored molecular level of human individuality, this atlas of human CoRSIVs provides a resource for future population-based investigations into how interindividual epigenetic variation modulates risk of disease. Electronic supplementary material The online version of this article (10.1186/s13059-019-1708-1) contains supplementary material, which is available to authorized users.

BackgroundMonozygotic twins have long been studied to estimate heritability and explore epigenetic influences on phenotypic variation. The phenotypic and epigenetic similarities of monozygotic twins have been assumed to be largely due to their genetic identity.ResultsHere, by analyzing data from a genome-scale study of DNA methylation in monozygotic and dizygotic twins, we identified genomic regions at which the epigenetic similarity of monozygotic twins is substantially greater than can be explained by their genetic identity. This “epigenetic supersimilarity” apparently results from locus-specific establishment of epigenotype prior to embryo cleavage during twinning. Epigenetically supersimilar loci exhibit systemic interindividual epigenetic variation and plasticity to periconceptional environment and are enriched in sub-telomeric regions. In case-control studies nested in a prospective cohort, blood DNA methylation at these loci years before diagnosis is associated with risk of developing several types of cancer.ConclusionsThese results establish a link between early embryonic epigenetic development and adult disease. More broadly, epigenetic supersimilarity is a previously unrecognized phenomenon that may contribute to the phenotypic similarity of monozygotic twins.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1374-0) contains supplementary material, which is available to authorized users.

We have tracked DNA methylation dynamics at metastable epialleles across multiple developmental stages in the early human embryo.

Coupling bisulfite conversion with next-generation sequencing (Bisulfite-seq) enables genome-wide measurement of DNA methylation, but poses unique challenges for mapping. However, despite a proliferation of Bisulfite-seq mapping tools, no systematic comparison of their genomic coverage and quantitative accuracy has been reported. We sequenced bisulfite-converted DNA from two tissues from each of two healthy human adults and systematically compared five widely used Bisulfite-seq mapping algorithms: Bismark, BSMAP, Pash, BatMeth and BS Seeker. We evaluated their computational speed and genomic coverage and verified their percentage methylation estimates. With the exception of BatMeth, all mappers covered >70% of CpG sites genome-wide and yielded highly concordant estimates of percentage methylation (r2 ≥ 0.95). Fourfold variation in mapping time was found between BSMAP (fastest) and Pash (slowest). In each library, 8–12% of genomic regions covered by Bismark and Pash were not covered by BSMAP. An experiment using simulated reads confirmed that Pash has an exceptional ability to uniquely map reads in genomic regions of structural variation. Independent verification by bisulfite pyrosequencing generally confirmed the percentage methylation estimates by the mappers. Of these algorithms, Bismark provides an attractive combination of processing speed, genomic coverage and quantitative accuracy, whereas Pash offers considerably higher genomic coverage.

Iron deficiency anemia (IDA) is the most prevalent nutritional condition worldwide. We studied the contribution of hepcidin-mediated iron blockade to IDA in African children. We measured hepcidin and hemoglobin weekly, and hematological, inflammatory, and iron biomarkers at baseline, 7 weeks, and 12 weeks in 407 anemic (hemoglobin < 11 g/dl), otherwise healthy Gambian children (6 to 27 months). Each child maintained remarkably constant hepcidin levels (P < 0.0001 for between-child variance), with half consistently maintaining levels that indicate physiological blockade of iron absorption. Hepcidin was strongly predicted by nurse-ascribed adverse events with dominant signals from respiratory infections and fevers (all P < 0.0001). Diarrhea and fecal calprotectin were not associated with hepcidin. In multivariate analysis, C-reactive protein was the dominant predictor of hepcidin and contributed to iron blockade even at very low levels. We conclude that even low-grade inflammation, especially associated with respiratory infections, contributes to IDA in African children.

Transposable elements (TEs) are enriched in cytosine methylation, preventing their mobility within the genome. We previously identified a genome-wide repertoire of candidate intracisternal A particle (IAP) TEs in mice that exhibit inter-individual variability in this methylation (VM-IAPs) with implications for genome function. Here we validate these metastable epialleles and discover a novel class that exhibit tissue specificity (tsVM-IAPs) in addition to those with uniform methylation in all tissues (constitutive- or cVM-IAPs); both types have the potential to regulate genes in cis. Screening for variable methylation at other TEs shows that this phenomenon is largely limited to IAPs, which are amongst the youngest and most active endogenous retroviruses. We identify sequences enriched within cVM-IAPs, but determine that these are not sufficient to confer epigenetic variability. CTCF is enriched at VM-IAPs with binding inversely correlated with DNA methylation. We uncover dynamic physical interactions between cVM-IAPs with low methylation ranges and other genomic loci, suggesting that VM-IAPs have the potential for long-range regulation. Our findings indicate that a recently evolved interplay between genetic sequence, CTCF binding, and DNA methylation at young TEs can result in inter-individual variability in transcriptional outcomes with implications for phenotypic variation.

Septicemia is a leading cause of death among neonates in low-income settings, a situation that is deteriorating due to high levels of antimicrobial resistance. Novel interventions are urgently needed. Iron stimulates the growth of most bacteria and hypoferremia induced by the acute phase response is a key element of innate immunity. Cord blood, which has high levels of hemoglobin, iron and transferrin saturation, has hitherto been used as a proxy for the iron status of neonates. We investigated hepcidin-mediated redistribution of iron in the immediate post-natal period and tested the effect of the observed hypoferremia on the growth of pathogens frequently associated with neonatal sepsis. Healthy, vaginally delivered neonates were enrolled in a cohort study at a single center in rural Gambia (N = 120). Cord blood and two further blood samples up to 96 hours of age were analyzed for markers of iron metabolism. Samples pooled by transferrin saturation were used to conduct ex-vivo growth assays with Staphylococcus aureus, Streptococcus agalactiae, Escherichia coli and Klebsiella pneumonia. A profound reduction in transferrin saturation occurred within the first 12 h of life, from high mean levels in cord blood (47.6% (95% CI 43.7–51.5%)) to levels at the lower end of the normal reference range by 24 h of age (24.4% (21.2–27.6%)). These levels remained suppressed to 48 h of age with some recovery by 96 h. Reductions in serum iron were associated with high hepcidin and IL-6 levels. Ex-vivo growth of all sentinel pathogens was strongly associated with serum transferrin saturation. These results suggest the possibility that the hypoferremia could be augmented (e.g. by mini-hepcidins) as a novel therapeutic option that would not be vulnerable to antimicrobial resistance. Trial registration: The original trial in which this study was nested is registered at ISRCTN, number 93854442.