Matt Silver

Maternal nutrition at conception modulates DNA methylation of human metastable epialleles

In experimental animals, maternal diet during the periconceptional period influences the establishment of DNA methylation at metastable epialleles in the offspring, with permanent phenotypic consequences. Pronounced naturally occurring seasonal differences in the diet of rural Gambian women allowed us to test this in humans. We show that significant seasonal variations in methyl-donor nutrient intake of mothers around the time of conception influence 13 relevant plasma biomarkers. The level of several of these maternal biomarkers predicts increased/decreased methylation at metastable epialleles in DNA extracted from lymphocytes and hair follicles in infants postnatally. Our results demonstrate that maternal nutritional status during early pregnancy causes persistent and systemic epigenetic changes at human metastable epialleles.

False positives in neuroimaging genetics using voxel-based morphometry data

Voxel-wise statistical inference is commonly used to identify significant experimental effects or group differences in both functional and structural studies of the living brain. Tests based on the size of spatially extended clusters of contiguous suprathreshold voxels are also widely used due to their typically increased statistical power. In “imaging genetics”, such tests are used to identify regions of the brain that are associated with genetic variation. However, concerns have been raised about the adequate control of rejection rates in studies of this type. A previous study tested the effect of a set of ‘null’ SNPs on brain structure and function, and found that false positive rates were well-controlled. However, no similar analysis of false positive rates in an imaging genetic study using cluster size inference has yet been undertaken. We measured false positive rates in an investigation of the effect of 700 pre-selected null SNPs on grey matter volume using voxel-based morphometry (VBM). As VBM data exhibit spatially-varying smoothness, we used both non-stationary and stationary cluster size tests in our analysis. Image and genotype data on 181 subjects with mild cognitive impairment were obtained from the Alzheimers Disease Neuroimaging Initiative (ADNI). At a nominal significance level of 5%, false positive rates were found to be well-controlled (3.9–5.6%), using a relatively high cluster-forming threshold, αc = 0.001, on images smoothed with a 12 mm Gaussian kernel. Tests were however anticonservative at lower cluster-forming thresholds (αc = 0.01, 0.05), and for images smoothed using a 6 mm Gaussian kernel. Here false positive rates ranged from 9.8 to 67.6%. In a further analysis, false positive rates using simulated data were observed to be well-controlled across a wide range of conditions. While motivated by imaging genetics, our findings apply to any VBM study, and suggest that parametric cluster size inference should only be used with high cluster-forming thresholds and smoothness. We would advocate the use of nonparametric methods in other cases.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment.

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.

Identification of gene pathways implicated in Alzheimer's disease using longitudinal imaging phenotypes with sparse regression.

We present a new method for the detection of gene pathways associated with a multivariate quantitative trait, and use it to identify causal pathways associated with an imaging endophenotype characteristic of longitudinal structural change in the brains of patients with Alzheimers disease (AD). Our method, known as pathways sparse reduced-rank regression (PsRRR), uses group lasso penalised regression to jointly model the effects of genome-wide single nucleotide polymorphisms (SNPs), grouped into functional pathways using prior knowledge of gene–gene interactions. Pathways are ranked in order of importance using a resampling strategy that exploits finite sample variability. Our application study uses whole genome scans and MR images from 99 probable AD patients and 164 healthy elderly controls in the Alzheimers Disease Neuroimaging Initiative (ADNI) database. 66,182 SNPs are mapped to 185 gene pathways from the KEGG pathway database. Voxel-wise imaging signatures characteristic of AD are obtained by analysing 3D patterns of structural change at 6, 12 and 24 months relative to baseline. High-ranking, AD endophenotype-associated pathways in our study include those describing insulin signalling, vascular smooth muscle contraction and focal adhesion. All of these have been previously implicated in AD biology. In a secondary analysis, we investigate SNPs and genes that may be driving pathway selection. High ranking genes include a number previously linked in gene expression studies to β-amyloid plaque formation in the AD brain (PIK3R3, PIK3CG, PRKCA and PRKCB), and to AD related changes in hippocampal gene expression (ADCY2, ACTN1, ACACA, and GNAI1). Other high ranking previously validated AD endophenotype-related genes include CR1, TOMM40 and APOE.

Exposure to aflatoxin B1 in utero is associated with DNA methylation in white blood cells of infants in The Gambia

BACKGROUND Exposure to environmental toxins during embryonic development may lead to epigenetic changes that influence disease risk in later life. Aflatoxin is a contaminant of staple foods in sub-Saharan Africa, is a known human liver carcinogen and has been associated with stunting in infants. METHODS We have measured aflatoxin exposure in 115 pregnant women in The Gambia and examined the DNA methylation status of white blood cells from their infants at 2-8 months old (mean 3.6 ± 0.9). Aflatoxin exposure in women was assessed using an ELISA method to measure aflatoxin albumin (AF-alb) adducts in plasma taken at 1-16 weeks of pregnancy. Genome-wide DNA methylation of infant white blood cells was measured using the Illumina Infinium HumanMethylation450beadchip. RESULTS AF-alb levels ranged from 3.9 to 458.4 pg/mg albumin. We found that aflatoxin exposure in the mothers was associated to DNA methylation in their infants for 71 CpG sites (false discovery rate < 0.05), with an average effect size of 1.7% change in methylation. Aflatoxin-associated differential methylation was observed in growth factor genes such as FGF12 and IGF1, and immune-related genes such as CCL28, TLR2 and TGFBI. Moreover, one aflatoxin-associated methylation region (corresponding to the miR-4520b locus) was identified. CONCLUSIONS This study shows that maternal exposure to aflatoxin during the early stages of pregnancy is associated with differential DNA methylation patterns of infants, including in genes related to growth and immune function. This reinforces the need for interventions to reduce aflatoxin exposure, especially during critical periods of fetal and infant development.

Vitamin D binding protein genotype is associated with plasma 25OHD concentration in West African children

Vitamin D is well known for its role in promoting skeletal health. Vitamin D status is determined conventionally by circulating 25-dihydroxyvitamin D (25OHD) concentration. There is evidence indicating that circulating 25OHD concentration is affected by variation in Gc, the gene encoding the vitamin D binding protein (DBP). The composite genotype of two single nucleotide polymorphisms (rs7041 and rs4588) results in different DBP isotypes (Gc1f, Gc1s and Gc2). The protein configurational differences among DBP isotypes affect DBP substrate binding affinity. The aims of this study were to determine 1) Gc variant frequencies in a population from an isolated rural region of The Gambia, West Africa (n = 3129) with year-round opportunity for cutaneous vitamin D synthesis and 2) the effects of Gc variants on 25OHD concentration (n = 237) in a genetically representative sub-group of children (mean (SD) age: 11.9 (4.8) years). The distribution of Gc variants was Gc1f: 0.86, Gc1s: 0.11 and Gc2: 0.03. The mean (SD) concentration of 25OHD was 59.6 (12.9) nmol/L and was significantly higher in those homozygous for Gc1f compared to other Gc variants (60.7 (13.1) vs. 56.6 (12.1) nmol/L, P = 0.03). Plasma 25OHD and 1,25(OH)2D concentration was significantly associated with parathyroid hormone in Gc1f-1f but not in the other Gc variants combined. This study demonstrates that different Gc variants are associated with different 25OHD concentrations in a rural Gambian population. Gc1f-1f, thought to have the highest affinity for 25OHD, had the highest 25OHD concentration compared with lower affinity Gc variants. The considerable difference in Gc1f frequency observed in Gambians compared with other non-West African populations and associated differences in plasma 25OHD concentration, may have implications for the way in which vitamin D status should be interpreted across different ancestral groups.

Possible relationship between common genetic variation and white matter development in a pilot study of preterm infants.

The consequences of preterm birth are a major public health concern with high rates of ensuing multisystem morbidity, and uncertain biological mechanisms. Common genetic variation may mediate vulnerability to the insult of prematurity and provide opportunities to predict and modify risk.

Epigenetic supersimilarity of monozygotic twin pairs

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.

Evidence for negative selection of gene variants that increase dependence on dietary choline in a Gambian cohort.

Choline is an essential nutrient, and the amount needed in the diet is modulated by several factors. Given geographical differences in dietary choline intake and disparate frequencies of single‐nucleotide polymorphisms (SNPs) in choline metabolism genes between ethnic groups, we tested the hypothesis that 3 SNPs that increase dependence on dietary choline would be under negative selection pressure in settings where choline intake is low: choline dehydrogenase (CHDH) rs12676, methylenetetrahydrofolate reductase 1 (MTHFD1) rs2236225, and phosphatidylethanolamine‐N‐methyltransferase (PEMT) rs12325817. Evidence of negative selection was assessed in 2 populations: one in The Gambia, West Africa, where there is historic evidence of a choline‐poor diet, and the other in the United States, with a comparatively choline‐rich diet. We used 2 independent methods, and confirmation of our hypothesis was sought via a comparison with SNP data from the Maasai, an East African population with a genetic background similar to that of Gambians but with a traditional diet that is higher in choline. Our results show that frequencies of SNPs known to increase dependence on dietary choline are significantly reduced in the low‐choline setting of The Gambia. Our findings suggest that adequate intake levels of choline may have to be reevaluated in different ethnic groups and highlight a possible approach for identifying novel functional SNPs under the influence of dietary selective pressure.—Silver, M. J., Corbin, K. D., Hellenthal, G., da Costa, K.‐A., Dominguez‐Salas, P., Moore, S. E., Owen, J., Prentice, A. M., Hennig, B. J., Zeisel, S. H. Evidence for negative selection of gene variants that increase dependence on dietary choline in a Gambian cohort. FASEB J. 29, 3426‐3435 (2015). www.fasebj.org

Estimation of a significance threshold for epigenome-wide association studies.

Epigenome‐wide association studies (EWAS) are designed to characterise population‐level epigenetic differences across the genome and link them to disease. Most commonly, they assess DNA‐methylation status at cytosine‐guanine dinucleotide (CpG) sites, using platforms such as the Illumina 450k array that profile a subset of CpGs genome wide. An important challenge in the context of EWAS is determining a significance threshold for declaring a CpG site as differentially methylated, taking multiple testing into account. We used a permutation method to estimate a significance threshold specifically for the 450k array and a simulation extrapolation approach to estimate a genome‐wide threshold. These methods were applied to five different EWAS datasets derived from a variety of populations and tissue types. We obtained an estimate of α=2.4×10−7 for the 450k array, and a genome‐wide estimate of α=3.6×10−8 . We further demonstrate the importance of these results by showing that previously recommended sample sizes for EWAS should be adjusted upwards, requiring samples between ∼10% and ∼20% larger in order to maintain type‐1 errors at the desired level.