Twin methodology in epigenetic studies

Birth weight and anthropometric measurements of twins

BACKGROUND

Genetic and environmental influences on anthropometric measures can be investigated by comparing dizygotic (DZ) versus monozygotic (MZ) twins. Investigating cohorts living in different geographical areas across the globe can identify the variation in heritability versus environment.

AIMS

(1) To investigate the association between birth weight and anthropometric measurements during adulthood; (2) to study the genetic and environmental influences on body measures including birth weight, weight and height among twins; and (3) to assess the variation in heritability versus environment among two cohorts of twins who lived in different geographical areas.

SUBJECTS AND METHODS

Twins were collected from two twin registers. Data on birth weight, adult weight and height in 430 MZ and 170 DZ twins living in two geographically distinct parts of the world were collected. A genetic analysis was performed using MX software.

RESULTS

Birth weight was associated with weight, height and BMI. Both MZ and DZ twins with low birth weight had shorter height during their adult life (p = 0.001), but only MZ twins with lower birth weight were lighter at adulthood (p = 0.001). Intra-pair differences in birth weight were not associated with differences in adult height (p = 0.366) or weight (p = 0.796). Additive genetic effects accounted for 53% of the variance in weight, 43% in height and 55% in birth weight. The remaining variance was attributed to unique environmental effects (15% for weight, 13% for height and 45% for birth weight and only 16% for BMI). Variability was found to be different in the two cohorts. The best fitting model for birth weight and BMI was additive genetic and non-shared environment and for weight and height was additive genetic, non-shared environment (plus common Environment).

CONCLUSIONS

Data suggests that the association between weight at birth and anthropometric measures in later life is influenced by both genetic and environmental factors. Living in different environments can potentially relate to variation found in the environment.

Genetics of Obesity Traits: A Bivariate Genome-Wide Association Analysis

Previous genome-wide association studies on anthropometric measurements have identified more than 100 related loci, but only a small portion of heritability in obesity was explained. Here we present a bivariate twin study to look for the genetic variants associated with body mass index and waist-hip ratio, and to explore the obesity-related pathways in Northern Han Chinese. Cholesky decomposition model for 242 monozygotic and 140 dizygotic twin pairs indicated a moderate genetic correlation (r = 0.53, 95%CI: 0.42-0.64) between body mass index and waist-hip ratio. Bivariate genome-wide association analysis in 139 dizygotic twin pairs identified 26 associated SNPs with p < 10-5. Further gene-based analysis found 291 nominally associated genes (P < 0.05), including F12, HCRTR1, PHOSPHO1, DOCK2, DOCK6, DGKB, GLP1R, TRHR, MMP1, GPR55, CCK, and OR2AK2, as well as 6 enriched gene-sets with FDR < 0.05. Expression quantitative trait loci analysis identified rs2242044 as a significant cis-eQTL in both the normal adipose-subcutaneous (P = 1.7 × 10-9) and adipose-visceral (P = 4.4 × 10-15) tissue. These findings may provide an important entry point to unravel genetic pleiotropy in obesity traits.

Circulating microRNAs disclose biology of normal cognitive function in healthy elderly people - a discovery twin study

Neurobiology is regulated by miRNA. Here circulating plasma miRNAs were assayed on a 754 miRNA OpenArray platform using 90 monozygotic elderly twins (73-95 year of age) and associated with mini mental state examination (MMSE) and a five-component cognitive score (CCS) in an explorative study. Both ordinary individual and twin-pair analyses were performed with level of cognitive scores. Candidate miRNAs were further associated with cognitive decline over 10 years using up to six repeated assessments. A total of 278 miRNAs were expressed in plasma from at least ten participants and 23 miRNAs were nominally associated (i.e., at an uncorrected p < 0.05) with CCS or MMSE in the paired analyses. Generally, elderly individuals with poor cognitive function had increase miRNA expression compared with equivalent individuals who performed better on the cognitive scale. Three miRNAs, miR-151a-3p, miR-212-3p and miR-1274b were associated with CCS both in the paired and the individual analysis. Four miRNAs found to be associated with CCS in cross-sectional analysis were also found to show an association in longitudinal analysis such that increase miRNA expression was associated with steeper cognitive decline. We propose a shared biological path underlies dementia and normative cognitive aging.

The role of DNA methylation in the association between childhood adversity and cardiometabolic disease

Growing evidence suggests that adverse environmental stimuli, especially during sensitive periods in early life, may lead to cardiometabolic disease in later life. However, the underlying biological mechanisms remain a mystery. Recent studies inferred that epigenetic modifications are likely involved. We review recent studies, primarily focused on the findings from human studies, to indicate the role of DNA methylation in the associations between childhood adversity and cardiometabolic disease in adulthood. In particular, we focused on DNA methylation modifications in genes regulating the hypothalamus-pituitary-adrenal axis as well as the immune system.

Weighted gene co-expression network analysis of expression data of monozygotic twins identifies specific modules and hub genes related to BMI

BACKGROUND

The therapeutic management of obesity is challenging, hence further elucidating the underlying mechanisms of obesity development and identifying new diagnostic biomarkers and therapeutic targets are urgent and necessary. Here, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) to identify significant genes and specific modules related to BMI based on gene expression profile data of 7 discordant monozygotic twins.

RESULTS

In the differential gene expression analysis, it appeared that 32 differentially expressed genes (DEGs) were with a trend of up-regulation in twins with higher BMI when compared to their siblings. Categories of positive regulation of nitric-oxide synthase biosynthetic process, positive regulation of NF-kappa B import into nucleus, and peroxidase activity were significantly enriched within GO database and NF-kappa B signaling pathway within KEGG database. DEGs of NAMPT, TLR9, PTGS2, HBD, and PCSK1N might be associated with obesity. In the WGCNA, among the total 20 distinct co-expression modules identified, coral1 module (68 genes) had the strongest positive correlation with BMI (r = 0.56, P = 0.04) and disease status (r = 0.56, P = 0.04). Categories of positive regulation of phospholipase activity, high-density lipoprotein particle clearance, chylomicron remnant clearance, reverse cholesterol transport, intermediate-density lipoprotein particle, chylomicron, low-density lipoprotein particle, very-low-density lipoprotein particle, voltage-gated potassium channel complex, cholesterol transporter activity, and neuropeptide hormone activity were significantly enriched within GO database for this module. And alcoholism and cell adhesion molecules pathways were significantly enriched within KEGG database. Several hub genes, such as GAL, ASB9, NPPB, TBX2, IL17C, APOE, ABCG4, and APOC2 were also identified. The module eigengene of saddlebrown module (212 genes) was also significantly correlated with BMI (r = 0.56, P = 0.04), and hub genes of KCNN1 and AQP10 were differentially expressed.

CONCLUSION

We identified significant genes and specific modules potentially related to BMI based on the gene expression profile data of monozygotic twins. The findings may help further elucidate the underlying mechanisms of obesity development and provide novel insights to research potential gene biomarkers and signaling pathways for obesity treatment. Further analysis and validation of the findings reported here are important and necessary when more sample size is acquired.

The epigenetic mechanism for discordance of autoimmunity in monozygotic twins

Monozygotic twins share an identical DNA sequence but are not truly "identical". In fact, when it comes to health and disease, they may often display some level of phenotypic discordance. The cause of this discordance is often unknown. Epigenetic modifications such as DNA methylation, histone modification, and microRNAs-mediated regulation regulate gene expression and are sensitive to external stimuli. These modifications may be seen to bridge the gap between genetics and the environment. Over the years, the importance of epigenetics as a primary mechanism for the role that the environment plays in defining phenotype has been increasingly appreciated. Mechanisms of epigenetics include DNA methylation, histone modifications and microRNAs. Discordance rates in monozygotic twins vary depending on the specific condition, from 11% in SLE to 64% in psoriasis and 77% in PBC. Other autoimmune diseases in which discordance is found among monozygotic twins has also been studied include type 1 diabetes, multiple sclerosis, rheumatoid arthritis, dermatomyositis and systemic sclerosis. In some cases, the differences in various epigenetic modifications is slight, even though the concordance rate is low, suggesting that epigenetics is not the only factor that needs to be considered. Nonetheless, the study of phenotypic discordance in monozygotic twins may shed light on the pathogenesis of autoimmune diseases and contribute to the development of new methodologies for the diagnosis and treatment of these diseases.

Genome-wide analysis of sperm DNA methylation from monozygotic twin bulls

Monozygotic (MZ) twins are of great interest to elucidate the contributions of pre- and postnatal environmental factors on epigenetics in the expression of complex traits and diseases. Progeny testing recently revealed that MZ twin bulls do not necessarily lead to identical genetic merit estimates (i.e. breeding values). Therefore, to explain differences in offspring productivity of MZ twin bulls despite their identical genetic backgrounds, we hypothesised that paternal sperm epigenomes vary between MZ twin bulls. In the present study, semen characteristics and global sperm DNA methylome were profiled for four pairs of MZ twin bulls. Some MZ twin pairs had divergent semen quality (sperm morphology, motility and viability). Comparative genome-wide DNA methylome surveys were performed using methyl-sensitive enrichment and microarray identification. Between 2% and 10% of all probes (400 000) were differentially methylated between MZ twin pairs. In addition, there were 580 loci differentially methylated across all pairs of MZ twins. Furthermore, enrichment analysis indicated a significant enrichment for fertility associated quantitative trait loci (P = 0.033). In conclusion, differences in the sperm epigenome may contribute to incongruous diverging performances of daughters sired by bulls that are MZ twins.

Disease-Concordant Twins Empower Genetic Association Studies

Genome-wide association studies with moderate sample sizes are underpowered, especially when testing SNP alleles with low allele counts, a situation that may lead to high frequency of false-positive results and lack of replication in independent studies. Related individuals, such as twin pairs concordant for a disease, should confer increased power in genetic association analysis because of their genetic relatedness. We conducted a computer simulation study to explore the power advantage of the disease-concordant twin design, which uses singletons from disease-concordant twin pairs as cases and ordinary healthy samples as controls. We examined the power gain of the twin-based design for various scenarios (i.e., cases from monozygotic and dizygotic twin pairs concordant for a disease) and compared the power with the ordinary case-control design with cases collected from the unrelated patient population. Simulation was done by assigning various allele frequencies and allelic relative risks for different mode of genetic inheritance. In general, for achieving a power estimate of 80%, the sample sizes needed for dizygotic and monozygotic twin cases were one half and one fourth of the sample size of an ordinary case-control design, with variations depending on genetic mode. Importantly, the enriched power for dizygotic twins also applies to disease-concordant sibling pairs, which largely extends the application of the concordant twin design. Overall, our simulation revealed a high value of disease-concordant twins in genetic association studies and encourages the use of genetically related individuals for highly efficiently identifying both common and rare genetic variants underlying human complex diseases without increasing laboratory cost.

Common epigenetic variation in a European population of mentally healthy young adults

DNA methylation represents an important link between structural genetic variation and complex phenotypes. The study of genome-wide CpG methylation and its relation to traits relevant to psychiatry has become increasingly important. Here, we analyzed quality metrics of 394,043 CpG sites in two samples of 568 and 319 mentally healthy young adults. For 25% of all CpGs we observed medium to large common epigenetic variation. These CpGs were overrepresented in open sea and shore regions, as well as in intergenic regions. They also showed a strong enrichment of significant hits in association analyses. Furthermore, a significant proportion of common DNA methylation is at least partially genetically driven and thus may be observed similarly across tissues. These findings could be of particular relevance for studies of complex neuropsychiatric traits, which often rely on proxy tissues.

Differentially Methylated DNA Regions in Monozygotic Twin Pairs Discordant for Rheumatoid Arthritis: An Epigenome-Wide Study

Objectives

In an explorative epigenome-wide association study (EWAS) to search for gene independent, differentially methylated DNA positions and regions (DMRs) associated with rheumatoid arthritis (RA) by studying monozygotic (MZ) twin pairs discordant for RA.

Methods

Genomic DNA was isolated from whole blood samples from 28 MZ twin pairs discordant for RA. DNA methylation was measured using the HumanMethylation450 BeadChips. Smoking, anti-cyclic citrullinated peptide antibodies, and immunosuppressive treatment were included as covariates. Pathway analysis was performed using GREAT.

Results

Smoking was significantly associated with hypomethylation of a DMR overlapping the promoter region of the RNF5 and the AGPAT1, which are implicated in inflammation and autoimmunity, whereas DMARD treatment induced hypermethylation of the same region. Additionally, the promotor region of both S100A6 and EFCAB4B were hypomethylated, and both genes have previously been associated with RA. We replicated several candidate genes identified in a previous EWAS in treatment-naïve RA singletons. Gene-set analysis indicated the involvement of immunologic signatures and cancer-related pathways in RA.

Conclusion

We identified several differentially methylated regions associated with RA, which may represent environmental effects or consequences of the disease and plausible biological pathways pertinent to the pathogenesis of RA.

Heritable DNA Methylation in CD4+ Cells among Complex Families Displays Genetic and Non-Genetic Effects

DNA methylation at CpG sites is both heritable and influenced by environment, but the relative contributions of each to DNA methylation levels are unclear. We conducted a heritability analysis of CpG methylation in human CD4+ cells across 975 individuals from 163 families in the Genetics of Lipid-lowering Drugs and Diet Network (GOLDN). Based on a broad-sense heritability (H2) value threshold of 0.4, we identified 20,575 highly heritable CpGs among the 174,445 most variable autosomal CpGs (SD > 0.02). Tests for associations of heritable CpGs with genotype at 2,145,360 SNPs using 717 of 975 individuals showed that ~74% were cis-meQTLs (< 1 Mb away from the CpG), 6% of CpGs exhibited trans-meQTL associations (>1 Mb away from the CpG or located on a different chromosome), and 20% of CpGs showed no strong significant associations with genotype (based on a p-value threshold of 1e-7). Genes proximal to the genotype independent heritable CpGs were enriched for functional terms related to regulation of T cell activation. These CpGs were also among those that distinguished T cells from other blood cell lineages. Compared to genes proximal to meQTL-associated heritable CpGs, genotype independent heritable CpGs were moderately enriched in the same genomic regions that escape erasure during primordial germ cell development and could carry potential for generational transmission.

Epigenetic drift in the aging genome: a ten-year follow-up in an elderly twin cohort

BACKGROUND

Current epigenetic studies on aging are dominated by the cross-sectional design that correlates subjects' ages or age groups with their measured epigenetic profiles. Such studies have been more aimed at age prediction or building up the epigenetic clock of age rather than focusing on the dynamic patterns in epigenetic changes during the aging process.

METHODS

We performed an epigenome-wide association study of intra-individual longitudinal changes in DNA methylation at CpG (cytosine-phosphate-guanine) sites measured in whole-blood samples of a cohort of 43 elderly twin pairs followed for 10 years (age at intake 73-82 years). Biological pathway analysis and survival analysis were also conducted on CpGs showing longitudinal change in their DNA-methylation levels. Classical twin models were fitted to each CpG site to estimate the genetic and environmental effects on DNA-methylation.

RESULTS

Our analysis identified 2284 CpG sites whose DNA-methylation levels changed longitudinally over the follow-up. Twin modelling revealed that the longitudinal change for 90% of these CpG sites was explained solely by individual unique environmental factors and only for 10% of these sites was it influenced by familial factors (genetic or shared environment). Over 60% of the identified CpG sites were replicated (same direction and replication P < 0.05) in an independent cross-sectional sample of 300 twins aged from 30 to 74 years. The replication rate went up to 91% for the top 53 CpGs with P < 1 × 10-07. Pathway analysis of genes linked to these CpGs identified biologically meaningful gene-sets involved in cellular-signalling events and in transmission across chemical synapses, which are important molecular underpinnings of aging-related degenerative disorders.

CONCLUSION

Our epigenome-wide association studies on a cohort of old twins followed up for 10 years identified highly replicable epigenetic biomarkers predominantly implicated in signalling pathways of degenerative disorders and survival in the elderly.

Epigenetic Modifications of Major Depressive Disorder

Major depressive disorder (MDD) is a chronic disease whose neurological basis and pathophysiology remain poorly understood. Initially, it was proposed that genetic variations were responsible for the development of this disease. Nevertheless, several studies within the last decade have provided evidence suggesting that environmental factors play an important role in MDD pathophysiology. Alterations in epigenetics mechanism, such as DNA methylation, histone modification and microRNA expression could favor MDD advance in response to stressful experiences and environmental factors. The aim of this review is to describe genetic alterations, and particularly altered epigenetic mechanisms, that could be determinants for MDD progress, and how these alterations may arise as useful screening, diagnosis and treatment monitoring biomarkers of depressive disorders.

Contribution of Heritability and Epigenetic Factors to Skeletal Muscle Mass Variation in United Kingdom Twins

CONTEXT

Skeletal muscle mass (SMM) is one of the major components of human body composition, with deviations from normal values often leading to sarcopenia.

OBJECTIVE

Our major aim was to conduct a genome-wide DNA methylation study in an attempt to identify potential genomic regions associated with SMM.

DESIGN

This was a mixed cross-sectional and longitudinal study.

SETTING

Community-based study.

PARTICIPANTS

A total of 1550 middle-aged United Kingdom twins (monozygotic [MZ] and dizygotic [DZ]), 297 of which were repeatedly measured participated in the study.

MAIN OUTCOME MEASURE

Appendicular lean mass assessed using dual-energy X-ray absorptiometry technology, and methylated DNA immunoprecipitation sequencing DNA methylation profiling genome-wide were obtained from each individual.

RESULTS

Heritability estimate of SMM, with simultaneous adjustment for covariates obtained using variance decomposition analysis, was h(2) = 0.809 ± 0.050. After quality control and analysis of longitudinal stability, the DNA methylation data comprised of 723 029 genomic sites, with positive correlations between repeated measurements (Rrepeated = 0.114-0.905). Correlations between MZ and DZ twins were 0.51 and 0.38 at a genome-wide average, respectively, and clearly increased with Rrepeated. Testing for DNA methylation association with SMM in 50 discordant MZ twins revealed 36 081 nominally significant results, of which the top-ranked 134 signals (P < .01 and Rrepeated > 0.40) were subjected to replication in the sample of 1196 individuals. Seven SMM methylation association signals replicated at a false discovery rate less than 0.1, and these were located in or near genes DNAH12, CAND1, CYP4F29P, and ZFP64, which have previously been highlighted in muscle-related studies. Adjusting for age, smoking, and blood cell heterogeneity did not alter significance of these associations.

CONCLUSION

This epigenome-wide study, testing longitudinally stable methylation sites, discovered and replicated a number of associations between DNA methylation at CpG loci and SMM. Four replicated signals were related to genes with potential muscle functions, suggesting that the methylome of whole blood may be informative of SMM variation.

Environmental Epigenetics and a Unified Theory of the Molecular Aspects of Evolution: A Neo-Lamarckian Concept that Facilitates Neo-Darwinian Evolution

Environment has a critical role in the natural selection process for Darwinian evolution. The primary molecular component currently considered for neo-Darwinian evolution involves genetic alterations and random mutations that generate the phenotypic variation required for natural selection to act. The vast majority of environmental factors cannot directly alter DNA sequence. Epigenetic mechanisms directly regulate genetic processes and can be dramatically altered by environmental factors. Therefore, environmental epigenetics provides a molecular mechanism to directly alter phenotypic variation generationally. Lamarck proposed in 1802 the concept that environment can directly alter phenotype in a heritable manner. Environmental epigenetics and epigenetic transgenerational inheritance provide molecular mechanisms for this process. Therefore, environment can on a molecular level influence the phenotypic variation directly. The ability of environmental epigenetics to alter phenotypic and genotypic variation directly can significantly impact natural selection. Neo-Lamarckian concept can facilitate neo-Darwinian evolution. A unified theory of evolution is presented to describe the integration of environmental epigenetic and genetic aspects of evolution.