DNA-Methylation Patterns in Trisomy 21 Using Cells from Monozygotic Twins

DNA methylation is essential in mammalian development. We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in Down syndrome (DS). In order to determine the methylation differences in T21 without interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. Reduced Representation Bisulfite Sequencing (RRBS) revealed 35 differentially methylated promoter regions (DMRs) (Absolute methylation differences = 25%, FDR < 0.001) in MZ twins discordant for T21 that have also been observed in comparison between unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis (FDR = 1.60 e -03) and include genes of the HOXB and HOXD clusters. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) observed in the iPSC of the T21 versus normal twin. Altogether, the results of this study highlight the epigenetic effects of the extra chromosome 21 in T21 on loci outside of this chromosome that are relevant to DS associated phenotypes.     

Spatial,temporal and interindividual epigenetic variation of functionally important DNA methylation patterns

DNA methylation is an epigenetic modification that plays an important role in gene regulation. It can be influenced by stochastic events, environmental factors and developmental programs. However, little is known about the natural variation of gene-specific methylation patterns. In this study, we performed quantitative methylation analyses of six differentially methylated imprinted genes (H19, MEG3, LIT1, NESP55, PEG3 and SNRPN), one hypermethylated pluripotency gene (OCT4) and one hypomethylated tumor suppressor gene (APC) in chorionic villus, fetal and adult cortex, and adult blood samples. Both average methylation level and range of methylation variation depended on the gene locus, tissue type and/or developmental stage. We found considerable variability of functionally important methylation patterns among unrelated healthy individuals and a trend toward more similar methylation levels in monozygotic twins than in dizygotic twins. Imprinted genes showed relatively little methylation changes associated with aging in individuals who are >25 years. The relative differences in methylation among neighboring CpGs in the generally hypomethylated APC promoter may not only reflect stochastic fluctuations but also depend on the tissue type. Our results are consistent with the view that most methylation variation may arise after fertilization, leading to epigenetic mosaicism.     

A longitudinal twin study of skewed X chromosome-inactivation

X-chromosome inactivation (XCI) is a pivotal epigenetic mechanism involved in the dosage compensation of X-linked genes between males and females. In any given cell, the process of XCI in early female development is thought to be random across alleles and clonally maintained once established. Recent studies, however, suggest that XCI might not always be random and that skewed inactivation may become more prevalent with age. The factors influencing such XCI skewing and its changes over time are largely unknown. To elucidate the influence of stochastic, heritable and environmental factors in longitudinal changes in XCI, we examined X inactivation profiles in a sample of monozygotic (MZ) (n = 23) and dizygotic (DZ) (n = 22) female twin-pairs at ages 5 and 10 years. Compared to MZ twins who were highly concordant for allelic XCI ratios, DZ twins showed much lower levels of concordance. Whilst XCI patterns were moderately stable between ages 5 and 10 years, there was some drift over time with an increased prevalence of more extreme XCI skewing at age 10. To our knowledge, this study represents the earliest longitudinal assessment of skewed XCI patterns, and suggests that skewed XCI may already be established in early childhood. Our data also suggest a link between MZ twinning and the establishment of allelic XCI ratios, and demonstrate that acquired skewing in XCI after establishment is primarily mediated by stochastic mechanisms. These data have implications for our understanding about sex differences in complex disease, and the potential causes of phenotypic discordance between MZ female twins.     

Application of Wavelet Packet Transform to detect genetic polymorphisms by the analysis of inter-Alu PCR patterns

Background  

The analysis of Inter-Alu PCR patterns obtained from human genomic DNA samples is a promising technique for a simultaneous analysis of many genomic loci flanked by Alu repetitive sequences in order to detect the presence of genetic polymorphisms. Inter-Alu PCR products may be separated and analyzed by capillary electrophoresis using an automatic sequencer that generates a complex pattern of peaks. We propose an algorithmic method based on the Haar-Walsh Wavelet Packet Transformation (WPT) for an efficient detection of fingerprint-type patterns generated by PCR-based methodologies. We have tested our algorithmic approach on inter-Alu patterns obtained from the genomic DNA of three couples of monozygotic twins, expecting that the inter-Alu patterns of each twins couple will show differences due to unavoidable experimental variability. On the contrary the differences among samples of different twins are supposed to originate from genetic variability. Our goal is to automatically detect regions in the inter-Alu pattern likely associated to the presence of genetic polymorphisms.     

Differentiating between monozygotic twins through DNA methylation-specific high-resolution melt curve analysis

Although short tandem repeat profiling is extremely powerful in identifying individuals from crime scene stains, it is unable to differentiate between monozygotic (MZ) twins. Efforts to address this include mutation analysis through whole genome sequencing and through DNA methylation studies. Methylation of DNA is affected by environmental factors; thus, as MZ twins age, their DNA methylation patterns change. This can be characterized by bisulfite treatment followed by pyrosequencing. However, this can be time-consuming and expensive; thus, it is unlikely to be widely used by investigators. If the sequences are different, then in theory the melting temperature should be different. Thus, the aim of this study was to assess whether high-resolution melt curve analysis can be used to differentiate between MZ twins. Five sets of MZ twins provided buccal swabs that underwent extraction, quantification, bisulfite treatment, polymerase chain reaction amplification and high-resolution melting curve analysis targeting two markers, Alu-E2F3 and Alu-SP. Significant differences were observed between all MZ twins targeting Alu-E2F3 and in four of five MZ twins targeting Alu-SP (P < 0.05). Thus, it has been demonstrated that bisulfite treatment followed by high-resolution melting curve analysis could be used to differentiate between MZ twins.     

Intra-Monozygotic Twin Pair Discordance and Longitudinal Variation of Whole-Genome Scale DNA Methylation in Adults

Monozygotic twins share identical genomic DNA and are indistinguishable using conventional genetic markers. Increasing evidence indicates that monozygotic twins are epigenetically distinct, suggesting that a comparison between DNA methylation patterns might be useful to approach this forensic problem. However, the extent of epigenetic discordance between healthy adult monozygotic twins and the stability of CpG loci within the same individual over a short time span at the whole-genome scale are not well understood. Here, we used Infinium HumanMethylation450 Beadchips to compare DNA methylation profiles using blood collected from 10 pairs of monozygotic twins and 8 individuals sampled at 0, 3, 6, and 9 months. Using an effective and unbiased method for calling differentially methylated (DM) CpG sites, we showed that 0.087%–1.530% of the CpG sites exhibit differential methylation in monozygotic twin pairs. We further demonstrated that, on whole-genome level, there has been no significant epigenetic drift within the same individuals for up to 9 months, including one monozygotic twin pair. However, we did identify a subset of CpG sites that vary in DNA methylation over the 9-month period. The magnitude of the intra-pair or longitudinal methylation discordance of the CpG sites inside the CpG islands is greater than those outside the CpG islands. The CpG sites located on shores appear to be more suitable for distinguishing between MZ twins.     

Additive,Epistatic, and Environmental Effects Through the Lens of Expression Variability QTL in a Twin Cohort

The expression of a gene can vary across individuals in the general population, as well as between monozygotic twins. This variable expression is assumed to be due to the influence of both genetic and nongenetic factors. Yet little evidence supporting this assumption has been obtained from empirical data. In this study, we used expression data from a large twin cohort to investigate the influences of genetic and nongenetic factors on variable gene expression. We focused on a set of expression variability QTL (evQTL)—i.e., genetic loci associated with the variance, as opposed to the mean, of gene expression. We identified evQTL for 99, 56, and 79 genes in lymphoblastoid cell lines, skin, and fat, respectively. The differences in gene expression, measured by the relative mean difference (RMD), tended to be larger between pairs of dizygotic (DZ) twins than between pairs of monozygotic (MZ) twins, showing that genetic background influenced the expression variability. Furthermore, a more profound RMD was observed between pairs of MZ twins whose genotypes were associated with greater expression variability than the RMD found between pairs of MZ twins whose genotypes were associated with smaller expression variability. This suggests that nongenetic (e.g., environmental) factors contribute to the variable expression. Lastly, we demonstrated that the formation of evQTL is likely due to partial linkages between eQTL SNPs that are additively associated with the mean of gene expression; in most cases, no epistatic effect is involved. Our findings have implications for understanding divergent sources of gene expression variability.     

Genome-Wide DNA Methylation and Gene Expression Analyses of Monozygotic Twins Discordant for Intelligence Levels

Human intelligence, as measured by intelligence quotient (IQ) tests, demonstrates one of the highest heritabilities among human quantitative traits. Nevertheless, studies to identify quantitative trait loci responsible for intelligence face challenges because of the small effect sizes of individual genes. Phenotypically discordant monozygotic (MZ) twins provide a feasible way to minimize the effects of irrelevant genetic and environmental factors, and should yield more interpretable results by finding epigenetic or gene expression differences between twins. Here we conducted array-based genome-wide DNA methylation and gene expression analyses using 17 pairs of healthy MZ twins discordant intelligently. ARHGAP18, related to Rho GTPase, was identified in pair-wise methylation status analysis and validated via direct bisulfite sequencing and quantitative RT-PCR. To perform expression profile analysis, gene set enrichment analysis (GSEA) between the groups of twins with higher IQ and their co-twins revealed up-regulated expression of several ribosome-related genes and DNA replication-related genes in the group with higher IQ. To focus more on individual pairs, we conducted pair-wise GSEA and leading edge analysis, which indicated up-regulated expression of several ion channel-related genes in twins with lower IQ. Our findings implied that these groups of genes may be related to IQ and should shed light on the mechanism underlying human intelligence.     

Novel Epigenetic Changes Unveiled by Monozygotic Twins Discordant for Smoking Habits

Exposure to cigarette smoking affects the epigenome and could increase the risk of developing diseases such as cancer and cardiovascular disorders. Changes in DNA methylation associated with smoking may help to identify molecular pathways that contribute to disease etiology. Previous studies are not completely concordant in the identification of differentially methylated regions in the DNA of smokers. We performed an epigenome-wide DNA methylation study in a group of monozygotic (MZ) twins discordant for smoking habits to determine the effect of smoking on DNA methylation. As MZ twins are considered genetically identical, this model allowed us to identify smoking-related DNA methylation changes independent from genetic components. We investigated the whole blood genome-wide DNA methylation profiles in 20 MZ twin pairs discordant for smoking habits by using the Illumina HumanMethylation450 BeadChip. We identified 22 CpG sites that were differentially methylated between smoker and non-smoker MZ twins by intra-pair analysis. We confirmed eight loci already described by other groups, located in AHRR, F2RL3, MYOG1 genes, at 2q37.1 and 6p21.33 regions, and also identified several new loci. Moreover, pathway analysis showed an enrichment of genes involved in GTPase regulatory activity. Our study confirmed the evidence of smoking-related DNA methylation changes, emphasizing that well-designed MZ twin models can aid the discovery of novel DNA methylation signals, even in a limited sample population.     

Methylomic analysis of monozygotic twins discordant for childhood psychotic symptoms

Childhood psychotic symptoms are associated with increased rates of schizophrenia, other psychiatric disorders, and suicide attempts in adulthood; thus, elucidating early risk indicators is crucial to target prevention efforts. There is considerable discordance for psychotic symptoms between monozygotic twins, indicating that child-specific non-genetic factors must be involved. Epigenetic processes may constitute one of these factors and have not yet been investigated in relation to childhood psychotic symptoms. Therefore, this study explored whether differences in DNA methylation at age 10 were associated with monozygotic twin discordance for psychotic symptoms at age 12. The Environmental Risk (E-Risk) Longitudinal Twin Study cohort of 2,232 children (1,116 twin pairs) was assessed for age-12 psychotic symptoms and 24 monozygotic twin pairs discordant for symptoms were identified for methylomic comparison. Children provided buccal samples at ages 5 and 10. DNA was bisulfite modified and DNA methylation was quantified using the Infinium HumanMethylation450 array. Differentially methylated positions (DMPs) associated with psychotic symptoms were subsequently tested in post-mortem prefrontal cortex tissue from adult schizophrenia patients and age-matched controls. Site-specific DNA methylation differences were observed at age 10 between monozygotic twins discordant for age-12 psychotic symptoms. Similar DMPs were not found at age 5. The top-ranked psychosis-associated DMP (cg23933044), located in the promoter of the C5ORF42 gene, was also hypomethylated in post-mortem prefrontal cortex brain tissue from schizophrenia patients compared to unaffected controls. These data tentatively suggest that epigenetic variation in peripheral tissue is associated with childhood psychotic symptoms and may indicate susceptibility to schizophrenia and other mental health problems.     

A monozygotic twin pair with Rett syndrome

Summary A five-year-old, monozygotic, Turkish female twin pair with Rett syndrome is described. The twins are almost completely concordant in all clinical signs. This observation suggests a genetic cause of Rett syndrome.     

Array-based profiling of reference-independent methylation status (aPRIMES) identifies frequent promoter methylation and consecutive downregulation of ZIC2 in pediatric medulloblastoma

Existing microarray-based approaches for screening of DNA methylation are hampered by a number of shortcomings, such as the introduction of bias by DNA copy-number imbalances in the test genome and negligence of tissue-specific methylation patterns. We developed a method designated array-based profiling of reference-independent methylation status (aPRIMES) that allows the detection of direct methylation status rather than relative methylation. Array-PRIMES is based on the differential restriction and competitive hybridization of methylated and unmethylated DNA by methylation-specific and methylation-sensitive restriction enzymes, respectively. We demonstrate the accuracy of aPRIMES in detecting the methylation status of CpG islands for different states of methylation. Application of aPRIMES to the DNA from desmoplastic medulloblastomas of monozygotic twins showed strikingly similar methylation profiles. Additional analysis of 18 sporadic medulloblastomas revealed an overall correlation between highly methylated tumors and poor clinical outcome and identified ZIC2 as a frequently methylated gene in pediatric medulloblastoma.     

Utilizing Twins as Controls for Non-Twin Case-Materials in Genome Wide Association Studies

Twin registries around the globe have collected DNA samples from large numbers of monozygotic and dizygotic twins. The twin sample collections are frequently used as controls in disease-specific studies together with non-twins. This approach is unbiased under the hypothesis that twins and singletons are comparable in terms of allele frequencies; i.e. there are no genetic variants associated with being a twin per se. To test this hypothesis we performed a genome-wide association study comparing the allele frequency of 572,352 single nucleotide polymorphisms (SNPs) in 1,413 monozygotic (MZ) and 5,451 dizygotic (DZ) twins with 3,720 healthy singletons. Twins and singletons have been genotyped using the same platform. SNPs showing association with being a twin at P-value < 1 × 10^-5 were selected for replication analysis in 1,492 twins (463 MZ and 1,029 DZ) and 1,880 singletons from Finland. No SNPs reached genome-wide significance (P-value < 5 × 10^-8) in the main analysis combining MZ and DZ twins. In a secondary analysis including only DZ twins two SNPs (rs2033541 close to ADAMTSL1 and rs4149283 close to ABCA1) were genome-wide significant after meta-analysis with the Finnish population. The estimated proportion of variance on the liability scale explained by all SNPs was 0.08 (P-value=0.003) when MZ and DZ were considered together and smaller for MZ (0.06, P-value=0.10) compared to DZ (0.09, P-value=0.003) when analyzed separately. In conclusion, twins and singletons can be used in genetic studies together with general population samples without introducing large bias. Further research is needed to explore genetic variances associated with DZ twinning.     

Adult monozygotic twins discordant for intra-uterine growth have indistinguishable genome-wide DNA methylation profiles

Background

Low birth weight is associated with an increased adult metabolic disease risk. It is widely discussed that poor intra-uterine conditions could induce long-lasting epigenetic modifications, leading to systemic changes in regulation of metabolic genes. To address this, we acquire genome-wide DNA methylation profiles from saliva DNA in a unique cohort of 17 monozygotic monochorionic female twins very discordant for birth weight. We examine if adverse prenatal growth conditions experienced by the smaller co-twins lead to long-lasting DNA methylation changes.

Results

Overall, co-twins show very similar genome-wide DNA methylation profiles. Since observed differences are almost exclusively caused by variable cellular composition, an original marker-based adjustment strategy was developed to eliminate such variation at affected CpGs. Among adjusted and unchanged CpGs 3,153 are differentially methylated between the heavy and light co-twins at nominal significance, of which 45 show sensible absolute mean β-value differences. Deep bisulfite sequencing of eight such loci reveals that differences remain in the range of technical variation, arguing against a reproducible biological effect. Analysis of methylation in repetitive elements using methylation-dependent primer extension assays also indicates no significant intra-pair differences.

Conclusions

Severe intra-uterine growth differences observed within these monozygotic twins are not associated with long-lasting DNA methylation differences in cells composing saliva, detectable with up-to-date technologies. Additionally, our results indicate that uneven cell type composition can lead to spurious results and should be addressed in epigenomic studies.     

X-chromosome inactivation in monkey embryos and pluripotent stem cells

Inactivation of one X chromosome in female mammals (XX) compensates for the reduced dosage of X-linked gene expression in males (XY). However, the inner cell mass (ICM) of mouse preimplantation blastocysts and their in vitro counterparts, pluripotent embryonic stem cells (ESCs), initially maintain two active X chromosomes (XaXa). Random X chromosome inactivation (XCI) takes place in the ICM lineage after implantation or upon differentiation of ESCs, resulting in mosaic tissues composed of two cell types carrying either maternal or paternal active X chromosomes. While the status of XCI in human embryos and ICMs remains unknown, majority of human female ESCs show non-random XCI. We demonstrate here that rhesus monkey ESCs also display monoallelic expression and methylation of X-linked genes in agreement with non-random XCI. However, XIST and other X-linked genes were expressed from both chromosomes in isolated female monkey ICMs indicating that ex vivo pluripotent cells retain XaXa. Intriguingly, the trophectoderm (TE) in preimplantation monkey blastocysts also expressed X-linked genes from both alleles suggesting that, unlike the mouse, primate TE lineage does not support imprinted paternal XCI. Our results provide insights into the species-specific nature of XCI in the primate system and reveal fundamental epigenetic differences between in vitro and ex vivo primate pluripotent cells.