On the power of epigenome-wide association studies using a disease-discordant twin design

An epigenome-wide association study of waist circumference in Chinese monozygotic twins

OBJECTIVES

Central obesity poses significant health risks because it increases susceptibility to multiple chronic diseases. Epigenetic features such as DNA methylation may be associated with specific obesity traits, which could help us understand how genetic and environmental factors interact to influence the development of obesity. This study aims to identify DNA methylation sites associated with the waist circumference (WC) in Northern Han Chinese population, and to elucidate potential causal relationships.

METHODS

A total of 59 pairs of WC discordant monozygotic twins (ΔWC >0) were selected from the Qingdao Twin Registry in China. Generalized estimated equation model was employed to estimate the methylation levels of CpG sites on WC. Causal relationships between methylation and WC were assessed through the examination of family confounding factors using FAmiliaL CONfounding (ICE FALCON). Additionally, the findings of the epigenome-wide analysis were corroborated in the validation stage.

RESULTS

We identified 26 CpG sites with differential methylation reached false discovery rate (FDR) < 0.05 and 22 differentially methylated regions (slk-corrected p < 0.05) strongly linked to WC. These findings provided annotations for 26 genes, with notable emphasis on MMP17, ITGA11, COL23A1, TFPI, A2ML1-AS1, MRGPRE, C2orf82, and NINJ2. ICE FALCON analysis indicated the DNA methylation of ITGA11 and TFPI had a causal effect on WC and vice versa (p < 0.05). Subsequent validation analysis successfully replicated 10 (p < 0.05) out of the 26 identified sites.

CONCLUSIONS

Our research has ascertained an association between specific epigenetic variations and WC in the Northern Han Chinese population. These DNA methylation features can offer fresh insights into the epigenetic regulation of obesity and WC as well as hints to plausible biological mechanisms.

Mediation role of DNA methylation in association between handgrip strength and cognitive function in monozygotic twins

Handgrip strength is a crucial indicator to monitor the change of cognitive function over time, but its mechanism still needs to be further explored. We sampled 59 monozygotic twin pairs to explore the potential mediating effect of DNA methylation (DNAm) on the association between handgrip strength and cognitive function. The initial step was the implementation of an epigenome-wide association analysis (EWAS) in the study participants, with the aim of identifying DNAm variations that are associated with handgrip strength. Following that, we conducted an assessment of the mediated effect of DNAm by the use of mediation analysis. In order to do an ontology enrichment study for CpGs, the GREAT program was used. There was a significant positive association between handgrip strength and cognitive function (β = 0.194, P < 0.001). The association between handgrip strength and DNAm of 124 CpGs was found to be statistically significant at a significance level of P < 1 × 10-4. Fifteen differentially methylated regions (DMRs) related to handgrip strength were found in genes such as SNTG2, KLB, CDH11, and PANX2. Of the 124 CpGs, 4 within KRBA1, and TRAK1 mediated the association between handgrip strength and cognitive function: each 1 kg increase in handgrip strength was associated with a potential decrease of 0.050 points in cognitive function scores, mediated by modifications in DNAm. The parallel mediating effect of these 4 CpGs was -0.081. The presence of DNAm variation associated with handgrip strength may play a mediated role in the association between handgrip strength and cognitive function.

Epigenome-wide association study identifies DNA methylation loci associated with handgrip strength in Chinese monozygotic twins

Background: The decline in muscle strength and function with aging is well recognized, but remains poorly characterized at the molecular level. Here, we report the epigenetic relationship between genome-wide DNA methylation and handgrip strength (HGS) among Chinese monozygotic (MZ) twins. Methods: DNA methylation (DNAm) profiling was conducted in whole blood samples through Reduced Representation Bisulfite Sequencing method. Generalized estimating equation was applied to regress the DNAm of each CpG with HGS. The Genomic Regions Enrichment of Annotations Tool was used to perform enrichment analysis. Differentially methylated regions (DMRs) were detected using comb-p. Causal inference was performed using Inference about Causation through Examination of Familial Confounding method. Finally, we validated candidate CpGs in community residents. Results: We identified 25 CpGs reaching genome-wide significance level. These CpGs located in 9 genes, especially FBLN1, RXRA, and ABHD14B. Many enriched terms highlighted calcium channels, neuromuscular junctions, and skeletal muscle organ development. We identified 21 DMRs of HGS, with several DMRs within FBLN1, SLC30A8, CST3, and SOCS3. Causal inference indicated that the DNAm of 16 top CpGs within FBLN1, RXRA, ABHD14B, MFSD6, and TYW1B might influence HGS, while HGS influenced DNAm at two CpGs within FBLN1 and RXRA. In validation analysis, methylation levels of six CpGs mapped to FLBN1 and one CpG mapped to ABHD14B were negatively associated with HGS weakness in community population. Conclusion: Our study identified multiple DNAm variants potentially related to HGS, especially CpGs within FBLN1 and ABHD14B. These findings provide new clues to the epigenetic modification underlying muscle strength decline.

Differential methylation patterns from clusters associated with glucose metabolism: evidence from a Shanghai twin study

Aim: To assess the associations between genome-wide DNA methylation (DNAm) and glucose metabolism among a Chinese population, in particular the multisite correlation. Materials & methods: Epigenome-wide associations with fasting plasma glucose (FPG) and hemoglobin A1c (HbA1c) were analyzed among 100 Shanghai monozygotic (MZ) twin pairs using the Infinium HumanMethylationEPIC v2.0 BeadChip. We conducted a Pearson's correlation test, hierarchical cluster and pairwise analysis to examine the differential methylation patterns from clusters. Results: Cg01358804 (TXNIP) was identified as the most significant site associated with FPG and HbA1c. Two clusters with hypermethylated and hypomethylated patterns were observed for both FPG and HbA1c. Conclusion: Differential methylation patterns from clusters may provide new clues for epigenetic changes and biological mechanisms in glucose metabolism.

Co-methylation analyses identify CpGs associated with lipid traits in Chinese discordant monozygotic twins

To control genetic background and early life milieu in genome-wide DNA methylation analysis for blood lipids, we recruited Chinese discordant monozygotic twins to explore the relationships between DNA methylations and total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). 132 monozygotic (MZ) twins were included with discordant lipid levels and completed data. A linear mixed model was conducted in Epigenome-wide association study (EWAS). Generalized estimating equation model was for gene expression analysis. We conducted Weighted correlation network analysis (WGCNA) to build co-methylated interconnected network. Additional Qingdao citizens were recruited for validation. Inference about Causation through Examination of Familial Confounding (ICE FALCON) was used to infer the possible direction of these relationships. A total of 476 top CpGs reached suggestively significant level (P < 10-4), of which, 192 CpGs were significantly associated with TG (FDR < 0.05). They were used to build interconnected network and highlight crucial genes from WGCNA. Finally, four CpGs in GATA4 were validated as risk factors for TC; six CpGs at ITFG2-AS1 were negatively associated with TG; two CpGs in PLXND1 played protective roles in HDL-C. ICE FALCON indicated abnormal TC was regarded as the consequence of DNA methylation in CpGs at GATA4, rather than vice versa. Four CpGs in ITFG2-AS1 were both causes and consequences of modified TG levels. Our results indicated that DNA methylation levels of 12 CpGs in GATA4, ITFG2-AS1, and PLXND1 were relevant to TC, TG, and HDL-C, respectively, which might provide new epigenetic insights into potential clinical treatment of dyslipidemia.

A Mass Spectrometry-Based Proteome Study of Twin Pairs Discordant for Incident Acute Myocardial Infarction within Three Years after Blood Sampling Suggests Novel Biomarkers

Acute myocardial infarction (AMI) is a major cause of mortality and morbidity worldwide, yet biomarkers for AMI in the short- or medium-term are lacking. We apply the discordant twin pair design, reducing genetic and environmental confounding, by linking nationwide registry data on AMI diagnoses to a survey of 12,349 twins, thereby identifying 39 twin pairs (48-79 years) discordant for their first-ever AMI within three years after blood sampling. Mass spectrometry of blood plasma identified 715 proteins. Among 363 proteins with a call rate > 50%, imputation and stratified Cox regression analysis revealed seven significant proteins (FDR < 0.05): FGD6, MCAM, and PIK3CB reflected an increased level in AMI twins relative to their non-AMI co-twins (HR > 1), while LBP, IGHV3-15, C1RL, and APOC4 reflected a decreased level in AMI twins relative to their non-AMI co-twins (HR < 1). Additional 50 proteins were nominally significant (p < 0.05), and bioinformatics analyses of all 57 proteins revealed biology within hemostasis, coagulation cascades, the immune system, and the extracellular matrix. A protein-protein-interaction network revealed Fibronectin 1 as a central hub. Finally, technical validation confirmed MCAM, LBP, C1RL, and APOC3. We put forward novel biomarkers for incident AMI, a part of the proteome field where markers are surprisingly rare and where additional studies are highly needed.

DNA Methylation Mediated the Association of Body Mass Index With Blood Pressure in Chinese Monozygotic Twins

Obesity is an established risk factor for hypertension, but the mechanisms are only partially understood. We examined whether body mass index (BMI)-related DNA methylation (DNAm) variation would mediate the association of BMI with blood pressure (BP). We first conducted a genomewide DNA methylation analysis in monozygotic twin pairs to detect BMI-related DNAm variation and then evaluated the mediating effect of DNAm on the relationship between BMI and BP levels using the causal inference test (CIT) method and mediation analysis. Ontology enrichment analysis was performed for CpGs using the GREAT tool. A total of 60 twin pairs for BMI and systolic blood pressure (SBP) and 58 twin pairs for BMI and diastolic blood pressure (DBP) were included. BMI was positively associated with SBP (β = 1.86, p = .0004). The association between BMI and DNAm of 85 CpGs reached p < 1×10-4 level. Eleven BMI-related differentially methylated regions (DMRs) within LNCPRESS1, OGDHL, RNU1-44P, NPHS1, ECEL1P2, LLGL2, RNY4P15, MOGAT3, PHACTR3, and BAI2 were found. Of the 85 CpGs, 9 mapped to C10orf71-AS1, NDUFB5P1, KRT80, BAI2, ABCA2, PEX11G and FGF4 were significantly associated with SBP levels. Of the 9 CpGs, 2 within ABCA2 negatively mediated the association between BMI and SBP, with a mediating effect of -0.24 (95% CI [-0.65, -0.01]). BMI was also positively associated with DBP (β = 0.60, p = .0495). The association between BMI and DNAm of 193 CpGs reached p < 1×10-4 level. Twenty-five BMI-related DMRs within OGDHL, POU4F2, ECEL1P2, TTC6, SMPD4, EP400, TUBA1C and AGAP2 were found. Of the 193 CpGs, 33 mapped to ABCA2, ADORA2B, CTNNBIP1, KDM4B, NAA60, RSPH6A, SLC25A19 and STIL were significantly associated with DBP levels. Of the 33 CpGs, 12 within ABCA2, SLC25A19, KDM4B, PTPRN2, DNASE1, TFCP2L1, LMNB2 and C10orf71-AS1 negatively mediated the association between BMI and DBP, with a total mediation effect of -0.66 (95% CI [-1.07, -0.30]). Interestingly, BMI might also negatively mediate the association between the DNAm of most CpG mediators mentioned above and BP. The mediating effect of DNAm was also found when stratified by sex. In conclusion, DNAm variation may partially negatively mediate the association of BMI with BP. Our findings may provide new clues to further elucidate the pathogenesis of obesity to hypertension and identify new diagnostic biomarkers and therapeutic targets for hypertension.

Genome-wide DNA methylation analysis of body composition in Chinese monozygotic twins

BACKGROUND

Little is currently known about epigenetic alterations associated with body composition in obesity. Thus, we aimed to explore epigenetic relationships between genome-wide DNA methylation levels and three common traits of body composition as measured by body fat percentage (BF%), fat mass (FM) and lean body mass (LBM) among Chinese monozygotic twins.

METHODS

Generalized estimated equation model was used to regress the methylation level of CpG sites on body composition. Inference about Causation Through Examination Of Familial Confounding was used to explore the evidence of a causal relationship. Gene expression analysis was further performed to validate the results of differentially methylated genes.

RESULTS

We identified 32, 22 and 28 differentially methylated CpG sites (p < 10-5 ) as well as 20, 17 and eight differentially methylated regions (slk-corrected p < 0.05) significantly associated with BF%, FM and LBM which were annotated to 65 genes, showing partially overlapping. Causal inference demonstrated bidirectional causality between DNA methylation and body composition (p < 0.05). Gene expression analysis revealed significant correlations between expression levels of five differentially methylated genes and body composition (p < 0.05).

CONCLUSIONS

These DNA methylation signatures will contribute to increased knowledge about the epigenetic basis of body composition and provide new strategies for early prevention and treatment of obesity and its related diseases.

Identification of key DNA methylation changes on fasting plasma glucose: a genome-wide DNA methylation analysis in Chinese monozygotic twins

BACKGROUND

Elevated fasting plasma glucose (FPG) levels can increase morbidity and mortality even when it is below the diagnostic threshold of type 2 diabetes mellitus (T2DM). We conducted a genome-wide DNA methylation analysis to detect DNA methylation (DNAm) variants potentially related to FPG in Chinese monozygotic twins.

METHODS

Genome-wide DNA methylation profiling in whole blood of twins was performed using Reduced Representation Bisulfite Sequencing (RRBS), yielding 551,447 raw CpGs. Association between DNAm of single CpG and FPG was tested using a generalized estimation equation. Differentially methylated regions (DMRs) were identified using comb-P approach. ICE FALCON method was utilized to perform the causal inference. Candidate CpGs were quantified and validated using Sequenom MassARRAY platform in a community population. Weighted gene co-expression network analysis (WGCNA) was conducted using gene expression data from twins.

RESULTS

The mean age of 52 twin pairs was 52 years (SD: 7). The relationship between DNAm of 142 CpGs and FPG reached the genome-wide significance level. Thirty-two DMRs within 24 genes were identified, including TLCD1, MRPS31P5, CASZ1, and CXADRP3. The causal relationship of top CpGs mapped to TLCD1, MZF1, PTPRN2, SLC6A18, ASTN2, IQCA1, GRIN1, and PDE2A genes with FPG were further identified using ICE FALCON method. Pathways potentially related to FPG were also identified, such as phospholipid-hydroperoxide glutathione peroxidase activity and mitogen-activated protein kinase p38 binding. Three CpGs mapped to SLC6A18 gene were validated in a community population, with a hypermethylated direction in diabetic patients. The expression levels of 18 genes (including SLC6A18 and TLCD1) were positively correlated with FPG levels.

CONCLUSIONS

We detect many DNAm variants that may be associated with FPG in whole blood, particularly the loci within SLC6A18 gene. Our findings provide important reference for the epigenetic regulation of elevated FPG levels and diabetes.

Epigenome-wide association study in Chinese monozygotic twins identifies DNA methylation loci associated with blood pressure

BACKGROUND

Hypertension is a crucial risk factor for developing cardiovascular disease and reducing life expectancy. We aimed to detect DNA methylation (DNAm) variants potentially related to systolic blood pressure (SBP) and diastolic blood pressure (DBP) by conducting epigenome-wide association studies in 60 and 59 Chinese monozygotic twin pairs, respectively.

METHODS

Genome-wide DNA methylation profiling in whole blood of twins was performed using Reduced Representation Bisulfite Sequencing, yielding 551,447 raw CpGs. Association between DNAm of single CpG and blood pressure was tested by applying generalized estimation equation. Differentially methylated regions (DMRs) were identified by comb-P approach. Inference about Causation through Examination of Familial Confounding was utilized to perform the causal inference. Ontology enrichment analysis was performed using Genomic Regions Enrichment of Annotations Tool. Candidate CpGs were quantified using Sequenom MassARRAY platform in a community population. Weighted gene co-expression network analysis (WGCNA) was conducted using gene expression data.

RESULTS

The median age of twins was 52 years (95% range 40, 66). For SBP, 31 top CpGs (p < 1 × 10^-4) and 8 DMRs were identified, with several DMRs within NFATC1, CADM2, IRX1, COL5A1, and LRAT. For DBP, 43 top CpGs (p < 1 × 10^-4) and 12 DMRs were identified, with several DMRs within WNT3A, CNOT10, and DAB2IP. Important pathways, such as Notch signaling pathway, p53 pathway by glucose deprivation, and Wnt signaling pathway, were significantly enriched for SBP and DBP. Causal inference analysis suggested that DNAm at top CpGs within NDE1, MYH11, SRRM1P2, and SMPD4 influenced SBP, while SBP influenced DNAm at CpGs within TNK2. DNAm at top CpGs within WNT3A influenced DBP, while DBP influenced DNAm at CpGs within GNA14. Three CpGs mapped to WNT3A and one CpG mapped to COL5A1 were validated in a community population, with a hypermethylated and hypomethylated direction in hypertension cases, respectively. Gene expression analysis by WGCNA further identified some common genes and enrichment terms.

CONCLUSION

We detect many DNAm variants that may be associated with blood pressure in whole blood, particularly the loci within WNT3A and COL5A1. Our findings provide new clues to the epigenetic modification underlying hypertension pathogenesis.

DNA methylation and waist-to-hip ratio: an epigenome-wide association study in Chinese monozygotic twins

PURPOSE

Epigenetic signatures such as DNA methylation may be associated with specific obesity traits. We performed an epigenome-wide association study (EWAS) by combining with the waist-to-hip ratio (WHR)-discordant monozygotic (MZ) twin design in an attempt to identify genetically independent DNA methylation marks associated with abdominal obesity in Northern Han Chinese and to determine the causation underlying.

METHODS

A total of 60 WHR discordant MZ twin pairs were selected from the Qingdao Twin Registry, China. Generalized estimated equation (GEE) model was used to regress the methylation level of CpG sites on WHR. The Inference about Causation through Examination of FAmiliaL CONfounding (ICE FALCON) was used to assess the temporal relationship between methylation and WHR. Gene expression analysis was conducted to validate the results of differentially methylated analyses.

RESULTS

EWAS identified 92 CpG sites with the level of P < 10^- 4 which were annotated to 32 genes, especially CADPS2, TUSC5, ZCCHC14, CORO7, COL23A1, CACNA1C, CYP26B1, and BCAT1. ICE FALCON showed significant causality between DNA methylation of several genes and WHR (P < 0.05). In region-based analysis, 14 differentially methylated regions (DMRs) located at 15 genes (slk-corrected P < 0.05) were detected. The gene expression analysis identified the significant correlation between expression levels of 5 differentially methylated genes and WHR (P < 0.05).

CONCLUSIONS

Our study identifies the associations between specific epigenetic variations and WHR in Northern Han Chinese. These DNA methylation signatures may have value as diagnostic biomarkers and provide novel insights into the molecular mechanisms of pathogenesis.

Genome-wide DNA methylation analysis of pulmonary function in middle and old-aged Chinese monozygotic twins

BACKGROUND

Previous studies have determined the epigenetic association between DNA methylation and pulmonary function among various ethnics, whereas this association is largely unknown in Chinese adults. Thus, we aimed to explore epigenetic relationships between genome-wide DNA methylation levels and pulmonary function among middle-aged Chinese monozygotic twins.

METHODS

The monozygotic twin sample was drawn from the Qingdao Twin Registry. Pulmonary function was measured by three parameters including forced expiratory volume the first second (FEV1), forced vital capacity (FVC), and FEV1/FVC ratio. Linear mixed effect model was used to regress the methylation level of CpG sites on pulmonary function. After that, we applied Genomic Regions Enrichment of Annotations Tool (GREAT) to predict the genomic regions enrichment, and used comb-p python library to detect differentially methylated regions (DMRs). Gene expression analysis was conducted to validate the results of differentially methylated analyses.

RESULTS

We identified 112 CpG sites with the level of P < 1 × 10^-4 which were annotated to 40 genes. We identified 12 common enriched pathways of three pulmonary function parameters. We detected 39 DMRs located at 23 genes, of which PRDM1 was related to decreased pulmonary function, and MPL, LTB4R2, and EPHB3 were related to increased pulmonary function. The gene expression analyses validated DIP2C, ASB2, SLC6A5, and GAS6 related to decreased pulmonary function.

CONCLUSION

Our DNA methylation sequencing analysis on identical twins provides new references for the epigenetic regulation on pulmonary function. Several CpG sites, genes, biological pathways and DMRs are considered as possible crucial to pulmonary function.

Plant epigenomics for extenuation of abiotic stresses: challenges and future perspectives

Climate change has escalated abiotic stresses, leading to adverse effects on plant growth and development, eventually having deleterious consequences on crop productivity. Environmental stresses induce epigenetic changes, namely cytosine DNA methylation and histone post-translational modifications, thus altering chromatin structure and gene expression. Stable epigenetic changes are inheritable across generations and this enables plants to adapt to environmental changes (epipriming). Hence, epigenomes serve as a good source of additional tier of variability for development of climate-smart crops. Epigenetic resources such as epialleles, epigenetic recombinant inbred lines (epiRILs), epigenetic quantitative trait loci (epiQTLs), and epigenetic hybrids (epihybrids) can be utilized in epibreeding for improving stress tolerance of crops. Epigenome engineering is also gaining momentum for developing sustainable epimarks associated with important agronomic traits. Different epigenome editing tools are available for creating, erasing, and reading such epigenetic codes in plant genomes. However, epigenome editing is still understudied in plants due to its complex nature. Epigenetic interventions such as epi-fingerprinting can be exploited in the near future for health and quality assessment of crops under stress conditions. Keeping in view the challenges and opportunities associated with this important technology, the present review intends to enhance understanding of stress-induced epigenetic changes in plants and its prospects for development of climate-ready crops.

Differential regulation of the DNA methylome in adults born during the Great Chinese Famine in 1959-1961

BACKGROUND

Extensive epidemiological studies have established the association between exposure to early-life adversity and health status and diseases in adults. Epigenetic regulation is considered as a key mediator for this phenomenon but analysis on humans is sparse. The Great Chinese Famine lasting from 1958 to 1961 is a natural string of disasters offering a precious opportunity for elucidating the underlying epigenetic mechanism of the long-term effect of early adversity.

METHODS

Using a high-throughput array platform for DNA methylome profiling, we conducted a case-control epigenome-wide association study on early-life exposure to Chinese famine in 79 adults born during 1959-1961 and compared to 105 unexposed subjects born 1963-1964.

RESULTS

The single CpG site analysis of whole epigenome revealed a predominant pattern of decreased DNA methylation levels associated with fetal exposure to famine. Four CpG sites were detected with p < 1e-06 (linked to EHMT1, CNR1, UBXN7 and ESM1 genes), 16 CpGs detected with 1e-06 < p < 1e-05 and 157 CpGs with 1e-05 < p < 1e-04, with a predominant pattern of hypomethylation. Functional annotation to genes and their enriched biological pathways mainly involved neurodevelopment, neuropsychological disorders and metabolism. Multiple sites analysis detected two top-rank differentially methylated regions harboring RNF39 on chromosome 6 and PTPRN2 on chromosome 7, both showing epigenetic association with stress-related conditions.

CONCLUSION

Early-life exposure to famine could mediate DNA methylation regulations that persist into adulthood with broad impacts in the activities of genes and biological pathways. Results from this study provide new clues to the epigenetic embedding of early-life adversity and its impacts on adult health.

Genome-wide DNA methylation and gene expression analyses in monozygotic twins identify potential biomarkers of depression

Depression is currently the leading cause of disability around the world. We conducted an epigenome-wide association study (EWAS) in a sample of 58 depression score-discordant monozygotic twin pairs, aiming to detect specific epigenetic variants potentially related to depression and further integrate with gene expression profile data. Association between the methylation level of each CpG site and depression score was tested by applying a linear mixed effect model. Weighted gene co-expression network analysis (WGCNA) was performed for gene expression data. The association of DNA methylation levels of 66 CpG sites with depression score reached the level of P < 1 × 10^-4. These top CpG sites were located at 34 genes, especially PTPRN2, HES5, GATA2, PRDM7, and KCNIP1. Many ontology enrichments were highlighted, including Notch signaling pathway, Huntington disease, p53 pathway by glucose deprivation, hedgehog signaling pathway, DNA binding, and nucleic acid metabolic process. We detected 19 differentially methylated regions (DMRs), some of which were located at GRIK2, DGKA, and NIPA2. While integrating with gene expression data, HELZ2, PTPRN2, GATA2, and ZNF624 were differentially expressed. In WGCNA, one specific module was positively correlated with depression score (r = 0.62, P = 0.002). Some common genes (including BMP2, PRDM7, KCNIP1, and GRIK2) and enrichment terms (including complement and coagulation cascades pathway, DNA binding, neuron fate specification, glial cell differentiation, and thyroid gland development) were both identified in methylation analysis and WGCNA. Our study identifies specific epigenetic variations which are significantly involved in regions, functional genes, biological function, and pathways that mediate depression disorder.

Early-life antibiotic use and risk of attention-deficit hyperactivity disorder and autism spectrum disorder: results of a discordant twin study

Background

Development of the gut-brain axis in early life may be disturbed by antibiotic use. It has been hypothesized that this disturbance may contribute to development of neurodevelopmental disorders, including autism spectrum disorder and attention-deficit hyperactivity disorder. We aimed to assess the association between antibiotic use in early life and the risk of developing attention-deficit hyperactivity disorder or autism spectrum disorder, while controlling for shared genetic and environmental factors in a discordant twin design.

Methods

We conducted a cohort study in twins (7–12 years; 25 781 twins) from the Netherlands Twin Register (NTR) and a replication study in the Childhood and Adolescent Twin Study in Sweden (CATSS; 7946 9-year-old twins). Antibiotic use was recorded before age 2 years. Attention-deficit hyperactivity disorder and autism spectrum disorder were parent-reported in the Netherlands Twin Register and register-based in the Childhood and Adolescent Twin Study in Sweden.

Results

Early-life antibiotic use was associated with increased risk of attention-deficit hyperactivity disorder development [pooled odds ratio (OR) 1.10, 95% confidence interval (CI) 1.02-1.17] and autism spectrum disorder (pooled OR 1.15, 95% CI 1.06-1.25) in a case-control design. When restricting to monozygotic twin pairs discordant for the outcome, associations disappeared for both disorders in both cohorts (attention-deficit hyperactivity disorder OR 0.90, 95% CI 0.48-1.69 and OR 0.80, 95% CI 0.37-1.76, and autism spectrum disorder OR 0.66, 95% CI 0.38-1.16 and OR 0.29, 95% CI 0.02-4.50, respectively).

Conclusions

Our findings suggest that the association between early-life antibiotic use and risk of attention-deficit hyperactivity and autism spectrum disorder may be confounded by shared familial environment and genetics.

Epigenome-wide association study of level and change in cognitive abilities from midlife through late life

Background

Epigenetic mechanisms are important in aging and may be involved in late-life changes in cognitive abilities. We conducted an epigenome-wide association study of leukocyte DNA methylation in relation to level and change in cognitive abilities, from midlife through late life in 535 Swedish twins.

Results

Methylation levels were measured with the Infinium Human Methylation 450 K or Infinium MethylationEPIC array, and all sites passing quality control on both arrays were selected for analysis (n = 250,816). Empirical Bayes estimates of individual intercept (age 65), linear, and quadratic change were obtained from latent growth curve models of cognitive traits and used as outcomes in linear regression models. Significant sites (p < 2.4 × 10^–7) were followed up in between-within twin pair models adjusting for familial confounding and full-growth modeling. We identified six significant associations between DNA methylation and level of cognitive abilities at age 65: cg18064256 (PPP1R13L) with processing speed and spatial ability; cg04549090 (NRXN3) with spatial ability; cg09988380 (POGZ), cg25651129 (-), and cg08011941 (ENTPD8) with working memory. The genes are involved in neuroinflammation, neuropsychiatric disorders, and ATP metabolism. Within-pair associations were approximately half that of between-pair associations across all sites. In full-growth curve models, associations between DNA methylation and cognitive level at age 65 were of small effect sizes, and associations between DNA methylation and longitudinal change in cognitive abilities of very small effect sizes.

Conclusions

Leukocyte DNA methylation was associated with level, but not change in cognitive abilities. The associations were substantially attenuated in within-pair analyses, indicating they are influenced in part by genetic factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01075-9.

Genome-wide DNA methylation analysis of cognitive function in middle and old-aged Chinese monozygotic twins

Cognitive ability plays an important role in mental and physical well-beings in the increasingly ageing populations. Here, based on a sample of 30 cognitive function-discordant monozygotic twin pairs, we aimed to detect specific epigenetic variants potentially related to cognitive function by conducting an epigenome-wide association study (EWAS). Association between methylation level of single CpG site with cognitive function score was tested by linear mixed effect model. Functions of cis-regulatory regions and ontology enrichments were predicted by Genomic Regions Enrichment of Annotations Tool (GREAT). Differentially methylated regions (DMRs) were detected by comb-p python library. A list of 28 CpG sites were identified to reach the level of P < 1 × 10^-4, and the strongest association (cor = 0.138, P = 2.549 × 10^-6) was detected for DNA CpG site (Chr17: 40,700,490 bp) located at HSD17B1P1. The identified 14,065 genomic CpG sites (P < 0.05) were mapped to 2646 genes, especially HSD17B1P1, CUL4A, INTS8, GFI1B, ZNF467, CDH15, and PSMA1. GREAT ontology enrichments mainly highlighted nicotine pharmacodynamics pathway, GABA-B receptor II/nicotinic acetylcholine receptor/hedgehog/endothelin/Wnt signaling pathways, Parkinson disease, Huntington disease, glycolysis, neuronal system, and toll-like receptor binding. We detected 15 DMRs located at/near 16 genes, especially LINC01551, LINC02282, and FAM32A. And 32 cognitive function-associated differentially methylated genes could be replicated, such as SHANK2, ABCA2, PRDM16, NCOR2, and INPP5A. Our EWAS in monozygotic twins identify specific epigenetic variations which are significantly involved in functional genes, biological function and pathways that mediate cognitive function. The findings provide clues to further identify new diagnostic biomarkers and therapeutic targets for cognitive dysfunction.

Targeted DNA methylation profiling reveals epigenetic signatures in peanut allergy

DNA methylation (DNAm) has been shown to play a role in mediating food allergy; however, the mechanism by which it does so is poorly understood. In this study, we used targeted next-generation bisulfite sequencing to evaluate DNAm levels in 125 targeted highly informative genomic regions containing 602 CpG sites on 70 immune-related genes to understand whether DNAm can differentiate peanut allergy (PA) versus nonallergy (NA). We found PA-associated DNAm signatures associated with 12 genes (7 potentially novel to food allergy, 3 associated with Th1/Th2, and 2 associated with innate immunity), as well as DNAm signature combinations with superior diagnostic potential compared with serum peanut–specific IgE for PA versus NA. Furthermore, we found that, following peanut protein stimulation, peripheral blood mononuclear cell (PBMCs) from PA participants showed increased production of cognate cytokines compared with NA participants. The varying responses between PA and NA participants may be associated with the interaction between the modification of DNAm and the interference of environment. Using Euclidean distance analysis, we found that the distances of methylation profile comprising 12 DNAm signatures between PA and NA pairs in monozygotic (MZ) twins were smaller than those in randomly paired genetically unrelated individuals, suggesting that PA-related DNAm signatures may be associated with genetic factors.

Systematic integrated analysis of genetic and epigenetic variation in diabetic kidney disease

Poor metabolic control and host genetic predisposition are critical for diabetic kidney disease (DKD) development. The epigenome integrates information from sequence variations and metabolic alterations. Here, we performed a genome-wide methylome association analysis in 500 subjects with DKD from the Chronic Renal Insufficiency Cohort for DKD phenotypes, including glycemic control, albuminuria, kidney function, and kidney function decline. We show distinct methylation patterns associated with each phenotype. We define methylation variations that are associated with underlying nucleotide variations (methylation quantitative trait loci) and show that underlying genetic variations are important drivers of methylation changes. We implemented Bayesian multitrait colocalization analysis (moloc) and summary data-based Mendelian randomization to systematically annotate genomic regions that show association with kidney function, methylation, and gene expression. We prioritized 40 loci, where methylation and gene-expression changes likely mediate the genotype effect on kidney disease development. Functional annotation suggested the role of inflammation, specifically, apoptotic cell clearance and complement activation in kidney disease development. Our study defines methylation changes associated with DKD phenotypes, the key role of underlying genetic variations driving methylation variations, and prioritizes methylome and gene-expression changes that likely mediate the genotype effect on kidney disease pathogenesis.

Differential long noncoding RNA profiling of BMI in twins

Aim: Many efforts have been deployed to identify genetic variants associated with BMI. Alternatively, we explore epigenetic contribution to BMI variation by focusing on long noncoding RNAs (lncRNAs) which represents a key layer of epigenetic control. Materials & methods: We analyzed lncRNA expression in whole blood of 229 monozygotic twin pairs in association with BMI using generalized estimating equations. Results & conclusion: Six lncRNA probes were identified as significant (false discovery rate <0.05), with BMI showing causal effects on the expression of the significant lncRNAs. Functional annotation of differential profiles identified Gene ontology biological processes including kidney development, regulations of lipid biosynthetic process, circadian rhythm, notch signaling, etc. Whole blood lncRNAs are significantly expressed in response to BMI variation.

An epigenome-wide association study of early-onset major depression in monozygotic twins

Major depression (MD) is a debilitating mental health condition with peak prevalence occurring early in life. Genome-wide examination of DNA methylation (DNAm) offers an attractive complement to studies of allelic risk given it can reflect the combined influence of genes and environment. The current study used monozygotic twins to identify differentially and variably methylated regions of the genome that distinguish twins with and without a lifetime history of early-onset MD. The sample included 150 Caucasian monozygotic twins between the ages of 15 and 20 (73% female; Mage = 17.52 SD = 1.28) who were assessed during a developmental stage characterized by relatively distinct neurophysiological changes. All twins were generally healthy and currently free of medications with psychotropic effects. DNAm was measured in peripheral blood cells using the Infinium Human BeadChip 450 K Array. MD associations with early-onset MD were detected at 760 differentially and variably methylated probes/regions that mapped to 428 genes. Genes and genomic regions involved neural circuitry formation, projection, functioning, and plasticity. Gene enrichment analyses implicated genes related to neuron structures and neurodevelopmental processes including cell-cell adhesion genes (e.g., PCDHA genes). Genes previously implicated in mood and psychiatric disorders as well as chronic stress (e.g., NRG3) also were identified. DNAm regions associated with early-onset MD were found to overlap genetic loci identified in the latest Psychiatric Genomics Consortium meta-analysis of depression. Understanding the time course of epigenetic influences during emerging adulthood may clarify developmental phases where changes in the DNA methylome may modulate individual differences in MD risk.

Exploratory analysis of age and sex dependent DNA methylation patterns on the X-chromosome in whole blood samples

BACKGROUND

Large numbers of autosomal sites are found differentially methylated in the aging genome. Due to analytical difficulties in dealing with sex differences in X-chromosome content and X-inactivation (XCI) in females, this has not been explored for the X chromosome.

METHODS

Using data from middle age to elderly individuals (age 55+ years) from two Danish cohorts of monozygotic twins and the Scottish Lothian Birth Cohort 1921, we conducted an X-chromosome-wide analysis of age-associated DNA methylation patterns with consideration of stably inferred XCI status.

RESULTS

Through analysing and comparing sex-specific X-linked DNA methylation changes over age late in life, we identified 123, 293 and 55 CpG sites significant (FDR < 0.05) only in males, only in females and in both sexes of Danish twins. All findings were significantly replicated in the two Danish twin cohorts. CpG sites escaping XCI are predominantly de-methylated with increasing age across cohorts. In contrast, CpGs highly methylated in both sexes are methylated even further with increasing age. Among the replicated sites in Danish samples, 16 (13%), 24 (8.2%) and 3 (5.5%) CpGs were further validated in LBC1921 (FDR < 0.05).

CONCLUSIONS

The X-chromosome of whole blood leukocytes displays age- and sex-dependent DNA methylation patterns in relation to XCI across cohorts.

Chimerism in health and potential implications on behavior: A systematic review

In this review, we focus on the phenomenon of chimerism and especially microchimerism as one of the currently underexplored explanations for differences in health and behavior. Chimerism is an amalgamation of cells from two or more unique zygotes within a single organism, with microchimerism defined by a minor cell population of <1%. This article first presents an overview of the primary techniques employed to detect and quantify the presence of microchimerism and then reviews empirical studies of chimerism in mammals including primates and humans. In women, male microchimerism, a condition suggested to be the result of fetomaternal exchange in utero, is relatively easily detected by polymerase chain reaction molecular techniques targeting Y-chromosomal markers. Consequently, studies of chimerism in human diseases have largely focused on diseases with a predilection for females including autoimmune diseases, and female cancers. We detail studies of chimerism in human diseases and also discuss some potential implications in behavior. Understanding the prevalence of chimerism and the associated health outcomes will provide invaluable knowledge of human biology and guide novel approaches for treating diseases.

Replicating associations between DNA methylation and body mass index in a longitudinal sample of older twins

BACKGROUND

There is an important interplay between epigenetic factors and body weight, and previous work has identified ten sites where DNA methylation is robustly associated with body mass index (BMI) cross-sectionally. However, interpretation of the associations is complicated by the substantial changes in BMI often occurring in late-life, and the fact that methylation is often driven by genetic variation. This study therefore investigated the longitudinal association between these ten sites and BMI from midlife to late-life, and whether associations persist after controlling for genetic factors.

METHODS

We used data from 535 individuals (mean age 68) in the Swedish Adoption/Twin Study of Aging (SATSA) with longitudinal measures of both DNA methylation from blood samples and BMI, spanning up to 20 years. Methylation levels were measured with the Infinium Human Methylation 450K or Infinium MethylationEpic array, with seven of the ten sites passing quality control. Latent growth curve models were applied to investigate longitudinal associations between methylation and BMI, and between-within models to study associations within twin pairs, thus adjusting for genetic factors.

RESULTS

Baseline DNA methylation levels at five of the seven sites were associated with BMI level at age 65 (cg00574958 [CPT1A]; cg11024682 [SREBF1]), and/or change (cg06192883 [MYO5C]; cg06946797 [RMI2]; cg08857797 [VPS25]). For four of the five sites, the associations remained comparable within twin pairs. However, the effects of cg06192883 were substantially attenuated within pairs. No change in DNA methylation was detected for any of the seven evaluated sites.

CONCLUSION

Five of the seven sites investigated were associated with late-life level and/or change in BMI. The effects for four of the sites remained similar when examined within twin pairs, indicating that the associations are mainly environmentally driven. However, the substantial attenuation in the association between cg06192883 and late-life BMI within pairs points to the importance of genetic factors in this association.

The epigenome of twins as a perfect laboratory for studying behavioural traits

The debate upon the relative importance of nature vs nurture in the development of human behaviour can be traced back to ancient times. Traditional epidemiology and genetic epidemiology have confirmed the association of environmental and genetic factors with behavioural traits. Current genomic studies are identifying genetic variants associated with various behavioural traits. However, exploring the relationship of abundant environmental factors with the complex epigenome that mediates human behaviour is just at its beginning. Identical twins can serve as perfect experiments for studying the environmental impact on behavioural epigenetics advantaged by enriched power in association analysis due to controlling of their genetic make-ups. Recent development in causal inference using twin-based models adds more values in twins. This review briefly introduces the various approaches in making use of twins in studying behavioural epigenetics from experiment design to practical applications. Exploring the epigenome of twins using the powerful twin-based study designs and analytical approaches will help identifying causal epigenetic markers mediating environmental exposures and behavioural traits enabling both pharmaceutical intervention and effective prevention.

Twin Registries Moving Forward and Meeting the Future: A Review

Twin registries have developed as a valuable resource for the study of many aspects of disease and society over the years in many different countries. A number of these registries include large numbers of twins with data collected at varying information levels for twin cohorts over the past several decades. More recent expansion of twin datasets has allowed for the collection of genetic data, together with many other levels of 'omic' information along with multiple demographic, physiological, health outcomes and other measures typically used in epidemiologic research. Other twin data sources outside these registries reflect research interests in particular aspects of disease or specific phenotypic assessment. Twin registries have the potential to play a key role in many aspects of the artificial intelligence/machine learning-driven projects of the future and will continue to keep adapting to the changing research landscape.