The objective of this study was to understand and measure epigenetic changes associated with the occurrence of CHDs by utilizing the discordant monozygotic twin model. A unique set of monozygotic twins discordant for double-outlet right ventricles (DORVs) was used for this multiomics study. The cardiac and muscle tissue samples from the twins were subjected to whole genome sequencing, whole genome bisulfite sequencing, RNA-sequencing and liquid chromatography-tandem mass spectrometry analysis. Sporadic DORV cases and control fetuses were used for validation. Global hypomethylation status was observed in heart tissue samples from the affected twins. Among 36,228 differentially methylated regions (DMRs), 1097 DMRs involving 1039 genes were located in promoter regions. A total of 419 genes, and lncRNA–mRNA pairs involved 30 genes, and 62 proteins were significantly differentially expressed. Multiple omics integrative analysis revealed that five genes, including BGN, COL1A1, COL3A1, FBLN5, and FLAN, and three pathways, including ECM-receptor interaction, focal adhesion and TGF-β signaling pathway, exhibited differences at all three levels. This study demonstrates a multiomics profile of discordant twins and explores the possible mechanism of DORV development. Global hypomethylation might be associated with the risk of CHDs. Specific genes and specific pathways, particularly those involving ECM–receptor interaction, focal adhesion and TGF–β signaling, might be involved in the occurrence of CHDs.

Phenotypic variation is the result of gene expression based on complex interaction between genetic and environmental factors. It is well known that genetic and environmental factors influence gene expression, but our understanding of their relative importance remains limited. To obtain a hint for the understanding of their contributions, we took advantage of monozygotic twins, as they share genetic and shared environmental factors but differ in nonshared factors, such as environmental differences and stochastic factors. In this study, we performed cap analysis of gene expression on three pairs of twins and clustered each individual based on their expression profiles of annotated genes. The dendrogram of annotated gene transcripts showed a monophyletic clade for each twin pair. We also analyzed the expression of retrotransposons, such as human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs), given their abundance in the genome. Clustering analyses demonstrated that HERV and LINE expression diverged even within monozygotic twin pairs. Thus, HERVs and LINEs are more susceptible to nonshared factors than annotated genes. Motif analysis of differentially expressed annotated genes suggests that specificity protein/Krüppel-like factor family transcription factors are involved in the expression divergence of annotated gene influenced by nonshared factors. Collectively, our findings suggest that expressions of annotated genes and retrotransposons are differently regulated, and that the expression of retrotransposons is more susceptible to nonshared factors than annotated genes.

Schizophrenia (SZ) is a severe, complex, and common mental disorder with high heritability (80%), an adult age of onset, and high discordance (∼50%) in monozygotic twins (MZ). Extensive studies on familial and non-familial cases have implicated a number of segregating mutations and de novo changes in SZ that may include changes to the mitochondrial genome. Yet, no single universally causal variant has been identified, highlighting its extensive genetic heterogeneity. This report specifically focuses on the assessment of changes in the mitochondrial genome in a unique set of monozygotic twins discordant (MZD) for SZ using blood. Genomic DNA from six pairs of MZD twins and two sets of parents (N = 16) was hybridized to the Affymetrix Human SNP Array 6.0 to assess mitochondrial DNA copy number (mtDNA-CN). Whole genome sequencing (WGS) and quantitative polymerase chain reaction (qPCR) was performed for a subset of MZD pairs and their parents and was also used to derive mtDNA-CN estimates. The WGS data were further analyzed to generate heteroplasmy (HP) estimates. Our results show that mtDNA-CN estimates for within-pair and mother-child differences were smaller than comparisons involving unrelated individuals, as expected. MZD twins showed discordance in mtDNA-CN estimates and displayed concordance in directionality of differences for mtDNA-CN across all technologies. Further, qPCR performed better than Affymetrix in estimating mtDNA-CN based on relatedness. No reliable differences in HP were detected between MZD twins. The within-MZD differences in mtDNA-CN observed represent postzygotic somatic changes that may contribute to discordance of MZ twins for diseases, including SZ.

Monozygotic twins (MZT) are 2.5 times more frequent in ART than in natural conceptions. A number of ART-related mechanisms have been probably linked with MZT. Studies that retrospectively analyze the time-lapse (TL) records resulting in MZT suggest that some morphokinetic traits of the inner cell mass and the trophectoderm could be predictors of MZT, but results are controversial. We present the complete TL record of one case of MZT that split itself at the very moment of the division into two cells, with one of the cells coming out through a hole in the zona pellucida (ZP). Both resulting embryos developed normally, and were vitrified. It is suggested that the hole in the ZP may facilitate the extrusion of some cells of the <day 4 embryo and that this cell development is not constrained by being inside the ZP. Despite the lack of the inhibition of the ZP itself or the influence of the other embryo cells, the totipotent cell was then able to develop correctly from the start. Moreover, the embryo inside the ZP compensated for the loss of this cell apparently without problems. Our findings are discussed in the context of previous literature and ethical problems are addressed.

Dermatoglyphs are epidermal ridge configurations on the fingers, palms and soles that are formed during fetal development, and therefore only the intrauterine environment can have any influence on their formation. This study aims at investigating the genetic and environmental contribution in determining quantitative dermatoglyphic traits in 32 monozygotic (MZ) and 35 dizygotic (DZ) same-sex twins from the Albanian population of Kosovo. All genetic analyses were run in the statistical program Mx. After assumptions testing, based on the pattern of MZ–DZ correlations, univariate models were fitted to the data in order to estimate additive genetic (A), common (C) and individual (E) environmental influences for all variables. The exception was the atd-angle for which a model with nonadditive genetic (D) influences was tested, since DZ correlations were less than half of MZ correlations. Goodness of fit of the full ACE or ADE model was compared to the saturated model. The fit of nested models (AE, CE, DE or E) was compared to the full models (ACE or ADE). Our results indicate that additive genetic component strongly contributes to individual differences in finger ridge counts (49−81%), and weakly (0−50%) on the formation of the palmar ridge counts between the palmar triradii a, b, c, and d. The specific pattern found for the atd-angle implies the impact of a nonadditive genetic component, possibly the effect of a major gene. Further, more powered studies are needed to confirm this pattern, especially for resolving the issue of the huge difference in MZ and DZ twin similarity for the atd-angle palmar trait.

Strong associations between neural tube defects (NTDs) and monozygotic (MZ) twinning have long been noted, and it has been suggested that NTD cases who do not present as MZ twins may be the survivors of MZ twinning events. We have recently shown that MZ twins carry a strong, distinctive DNA methylation signature and have developed an algorithm based on genomewide DNA methylation array data that distinguishes MZ twins from dizygotic twins and other relatives at well above chance level. We have applied this algorithm to published methylation data from five fetal tissues (placental chorionic villi, kidney, spinal cord, brain and muscle) collected from spina bifida cases (n = 22), anencephalic cases (n = 15) and controls (n = 19). We see no difference in signature between cases and controls, providing no support for a common etiological role of MZ twinning in NTDs. The strong associations therefore continue to await elucidation.

T1 and T2 relaxation times were examined in four pairs of monozygotic (MZ) twins discordant and concordant for schizophrenia with low and high genetic loading for the illness and five healthy control MZ twin pairs. Patients with schizophrenia (n = 11) showed significant prolongation in T1 relaxation times in the globus pallidus (GP) bilaterally (P < 0.005, Bonferroni corrected) when compared to 14 healthy MZ twins.

The giant panda, native to mountains of south-west China, is one of the world's rarest bear species and is subject to considerable conservation effort. In captivity, the proportion of twins accounts for 54% of the total number of births. To date, little is known about zygosity in panda populations — specifically, the proportion of monozygotic and dizygotic twins. In this study, we used 10 microsatellite markers for reliable zygosity testing, and the probability of monozygotic twins was 99.963% when all 10 markers were concordant. Out of 43 studied twin pairs, no MZ twins were found, indicating that there may be no identical panda twins (or the incidence is very low). We speculate that the fertilized eggs of giant pandas do not have the capability to split into two identical embryos, or that this ability is very poor, which is likely due to delayed implantation that is common in bear species. The results of this study deepen our understanding of giant panda breeding, yield insight into panda twins’ likely mechanism of formation, and reduce the uncertainty of individual identity in wild population surveys.

Familial monozygotic (MZ) twinning reports are rare around the world, and we report a four-generation pedigree with seven recorded pairs of female MZ twins. Whole-genome sequencing of seven family members was performed to explore the featured genetic factors in MZ twins. For variations specific to MZ twins, five novel variants were observed in the X chromosome. These candidates were used to explain the seemingly X-linked dominant inheritance pattern, and only one variant was exonic, located at the 5′UTR region of ZCCHC12 (chrX: 117958597, G > A). Besides, consistent mitochondrial DNA composition in the maternal linage precluded roles of mitochondria for this trait. In this pedigree, autosomes also contain diverse variations specific to MZ twins. Pathway analysis revealed a significant enrichment of genes carrying novel SNVs in the epithelial adherens junction-signaling pathway (p = .011), contributed by FGFR1, TUBB6, and MYH7B. Meanwhile, TBC1D22A, TRIOBP, and TUBB6, also carrying similar SNVs, were involved in the GTPase family-mediated signal pathway. Furthermore, gene-set enrichment analysis for 533 genes covered by copy number variations specific to MZ twins illustrated that the tight junction-signaling pathway was significantly enriched (p < .001). Therefore, the novel changes in the X chromosome and the provided candidate variants across autosomes may be responsible for MZ twinning, giving clues to increase our understanding about the underlying mechanism.

Hundreds of penetrant risk loci have been identified across different neurodevelopmental disorders (NDDs), and these often involve rare (<1% frequency) copy number variations (CNVs), which can involve one or more genes. Monozygotic (MZ) twin pairs are long thought to share 100% of their genomic information. However, genetic differences in the form of postzygotic somatic variants have been reported recently both in typically developing (TD) and in clinically discordant MZ pairs. We sought to investigate the contribution of rare CNVs in 100 twin pairs enriched for NDD phenotypes with a particular focus on postzygotic CNVs in MZ pairs discordant for autism spectrum disorder (ASD) using the Illumina Infinium PsychArray. In our sample, no postzygotic de novo CNVs were found in 55 MZ twin pairs, including the 13 pairs discordant for ASD. We did detect a higher rate of CNVs overlapping genes involved in disorders of the nervous system, such as a rare deletion affecting HNRNPU, in MZ pairs discordant and concordant for ASD in comparison with TD pairs (p = .02). Our results are in concordance with earlier findings that postzygotic de novo CNV events are typically rare in genomic DNA derived from saliva or blood, and suggests that the discordance of NDDs in our sample of twins is not explained by discordant CNVs. Still, studies investigating postzygotic variation in MZ discordant twins using DNA from different tissues and single cells and higher resolution genomics are needed in the future.

CDKN1C and KCNQ1OT1 are imprinted genes that might be potential regulators of placental development. This study investigated placental expressions of CDKN1C and KCNQ1OT1 in monozygotic twins with and without selective intrauterine growth restriction (sIUGR). Seventeen sIUGR and fifteen normal monozygotic(MZ) twin pairs were examined. Placental mRNA expressions of CDKN1C and KCNQ1OT1 were detected by real-time fluorescent quantitative PCR. CDKN1C protein expression was detected by immunohistochemical assay and Western-blotting. In the sIUGR group, smaller fetuses had a smaller share of the placenta, and CDKN1C protein expression was significantly increased while KCNQ1OT1 mRNA expression was significantly decreased. The CDKN1C/KCNQ1OT1 mRNA ratio was lower in the larger fetus than in the smaller fetus (p < .05). In the control group, CDKN1C protein expression showed no difference between larger and smaller fetuses, while KCNQ1OT1 mRNA expression was significantly lower in the larger fetus, and the CDKN1C/KCNQ1OT1 mRNA ratio was higher in the larger fetus than in the smaller fetus (p < .05). Our findings showed that pathogenesis of sIUGR may be related to the co-effect of the up-regulated protein expression of CDKN1C and down-regulated mRNA expression of KCNQ1OT1 in the placenta.

Gestational weight gain (GWG) has in numerous studies been associated with offspring birth weight (BW) and childhood weight. However, these associations might be explained by genetic confounding as offspring inherit their mother's genetic potential to gain weight. Furthermore, little is known about whether particular periods of pregnancy could influence offspring body weight differently. We therefore aimed to explore total and trimester-specific effects of GWG in monozygotic (MZ) twin mother-pairs on their offspring's BW, weight at 1 year and body mass index (BMI) at 5 and 10 years. MZ twin mothers born 1962–1975 were identified in national Swedish registers, and data on exposure and outcome variables was collected from medical records. We analyzed associations within and between twin pairs. We had complete data on the mothers’ GWG and offspring BW for 82 pairs. The results indicated that total, and possibly also second and third trimester GWG were associated with offspring BW within the twin pairs in the fully adjusted model (β = 0.08 z-score units, 95% CI: 0.001, 0.17; β = 1.32 z-score units, 95% CI: -0.29, 2.95; and β = 1.02 z-score units, 95% CI: -0.50, 2.54, respectively). Our findings, although statistically weak, suggested no associations between GWG and offspring weight or BMI during infancy or childhood. Our study suggests that total, and possibly also second and third trimester, GWG are associated with offspring BW when taking shared genetic and environmental factors within twin pairs into account. Larger family-based studies with long follow-up are needed to confirm our findings.

Monozygotic (MZ) twins stem from the same single fertilized egg and therefore share all their inherited genetic variation. This is one of the unequivocal facts on which genetic epidemiology and twin studies are based. To what extent this also implies that MZ twins share genotypes in adult tissues is not precisely established, but a common pragmatic assumption is that MZ twins are 100% genetically identical also in adult tissues. During the past decade, this view has been challenged by several reports, with observations of differences in post-zygotic copy number variations (CNVs) between members of the same MZ pair. In this study, we performed a systematic search for differences of CNVs within 38 adult MZ pairs who had been misclassified as dizygotic (DZ) twins by questionnaire-based assessment. Initial scoring by PennCNV suggested a total of 967 CNV discordances. The within-pair correlation in number of CNVs detected was strongly dependent on confidence score filtering and reached a plateau of r = 0.8 when restricting to CNVs detected with confidence score larger than 50. The top-ranked discordances were subsequently selected for validation by quantitative polymerase chain reaction (qPCR), from which one single ~120kb deletion in NRXN1 on chromosome 2 (bp 51017111–51136802) was validated. Despite involving an exon, no sign of cognitive/mental consequences was apparent in the affected twin pair, potentially reflecting limited or lack of expression of the transcripts containing this exon in nerve/brain.

The large body of literature on the association between blood pressure (BP) and cognitive functioning has yielded mixed results, possibly due to the presence of non-linear effects across age, or because BP affects specific brain areas differently, impacting more on some cognitive skills than on others. If a robust association was detected among BP and specific cognitive tasks, the causal nature of reported associations between BP and cognition could be investigated in twin data, which allow a test of alternative explanations, including genetic pleiotropy. The present study first examines the association between BP and cognition in a sample of 1,140 participants with an age range between 10 and 86 years. Linear and quadratic effects of systolic BP (SBP) and diastolic BP (DBP) on cognitive functioning were examined for 17 tests across five functions. Associations were corrected for effects of sex and linear and quadratic effects of age. Second, to test a causal model, data from 123 monozygotic (MZ) twin pairs were analyzed to test whether cognitive functioning of the twins with the higher BP was different from that of the co-twins with lower BP. Associations between BP and cognitive functioning were absent for the majority of the cognitive tests, with the exception of a lower speed of emotion identification and verbal reasoning in subjects with high diastolic BP. In the MZ twin pair analyses, no effects of BP on cognition were found. We conclude that in the population at large, BP level is not associated with cognitive functioning in a clinically meaningful way.

Monozygotic (MZ) twins provide a natural system for investigating developmental plasticity and the potential epigenetic origins of disease. A major difference in the intrauterine environment between MZ pairs is whether they share a common placenta or have separate placentas. Using DNA methylation measured at >400,000 points in the genome on the Illumina HumanMethylation450 array, we demonstrate that the co-twins of MZ pairs (average age of 14) that shared a common placenta (n = 18 pairs) have more similar DNA methylation levels in blood throughout the genome relative to those with separate placentas (n = 16 pairs). Functional annotation of the genomic regions that show significantly different correlation between monochorionic (MC) and dichorionic (DC) MZ pairs found an over-representation of genes involved in the regulation of transcription, neuronal development, and cellular differentiation. These results support the idea that prenatal environmental exposures may have a lasting effect on an individual's epigenetic landscape, and the potential for these changes to have functional consequences.

Diurnal preference is an individual's preference for daily activities and sleep timing and is strongly correlated with the underlying circadian clock and the sleep-wake cycle validating its use as an indirect circadian measure in humans. Recent research has implicated DNA methylation as a mechanism involved in the regulation of the circadian clock system in humans and other mammals. In order to evaluate the extent of epigenetic differences associated with diurnal preference, we examined genome-wide patterns of DNA methylation in DNA from monozygotic (MZ) twin-pairs discordant for diurnal preference. MZ twins were selected from a longitudinal twin study designed to investigate the interplay of genetic and environmental factors in the development of emotional and behavioral difficulties. Fifteen pairs of MZ twins were identified in which one member scored considerably higher on the Horne–Ostberg Morningness–Eveningness Questionnaire (MEQ) than the other. Genome-wide DNA methylation patterns were assessed in twins’ buccal cell DNA using the Illumina Infinium HumanMethylation450 BeadChips. Quality control and data pre-processing was undertaken using the wateRmelon package. Differentially methylated probes (DMPs) were identified using an analysis strategy taking into account both the significance and the magnitude of DNA methylation differences. Our data indicate that DNA methylation differences are detectable in MZ twins discordant for diurnal preference. Moreover, downstream gene ontology (GO) enrichment analysis on the top-ranked diurnal preference associated DMPs revealed significant enrichment of pathways that have been previously associated with circadian rhythm regulation, including cell adhesion processes and calcium ion binding.

Identification of individuals within pairs of monozygotic (MZ) twins remains unresolved using common forensic DNA typing technology. For some criminal cases involving MZ twins as suspects, the twins had to be released due to inability to identify which of the pair was the perpetrator. In this study, we performed a genome-wide scan on whole blood-derived DNA from four pairs of healthy phenotypically concordant MZ twins using the methylated DNA immunoprecipitation sequencing technology to identify candidate DNA methylation markers with capacity to distinguish MZ twins within a pair. We identified 38 differential methylation regions showing within-pair methylation differences in all four MZ pairs. These are all located in CpG islands, 17 of which are promoter-associated, 17 are intergenic islands, and four are intragenic islands. Genes associated with these markers are related with cell proliferation, differentiation, and growth and development, including zinc finger proteins, PRRX2, RBBP9, or are involved in G-protein signaling, such as the regulator of G-protein signaling 16. Further validation studies on additional MZ twins are now required to evaluate the broader utility of these 38 markers for forensic use.