De novo and inherited CNVs in MZ twin pairs selected for discordance and concordance on Attention Problems

Technical strategy for monozygotic twin discrimination by single-nucleotide variants

Monozygotic (MZ) twins are theoretically genetically identical. Although they are revealed to accumulate mutations after the zygote splits, discriminating between twin genomes remains a formidable challenge in the field of forensic genetics. Single-nucleotide variants (SNVs) are responsible for a substantial portion of genetic variation, thus potentially serving as promising biomarkers for the identification of MZ twins. In this study, we sequenced the whole genome of a pair of female MZ twins when they were 27 and 33 years old to approximately 30 × coverage using peripheral blood on an Illumina NovaSeq 6000 Sequencing System. Potentially discordant SNVs supported by whole-genome sequencing were validated extensively by amplicon-based targeted deep sequencing and Sanger sequencing. In total, we found nine bona fide post-twinning SNVs, all of which were identified in the younger genomes and found in the older genomes. None of the SNVs occurred within coding exons, three of which were observed in introns, supported by whole-exome sequencing results. A double-blind test was employed, and the reliability of MZ twin discrimination by discordant SNVs was endorsed. All SNVs were successfully detected when input DNA amounts decreased to 0.25 ng, and reliable detection was limited to seven SNVs below 0.075 ng input. This comprehensive analysis confirms that SNVs could serve as cost-effective biomarkers for MZ twin discrimination.

Genetic and phenotypic evidence of the predictive validity of preschool parent reports of hyperactivity/impulsivity and inattention

To determine the validity of parent reports (PRs) of ADHD in preschoolers, we assessed hyperactivity/impulsivity (HI) and inattention (IN) in 1114 twins with PRs at 1.5, 2.5, 4, 5, 14, 15, and 17 years, and teacher-reports at 6, 7, 9, 10, and 12. We examined if preschool PRs (1) predict high HI/IN trajectories, and (2) capture genetic contributions to HI/IN into adolescence. Group-based trajectory analyses identified three 6-17 years trajectories for both HI and IN, including small groups with high HI (N = 88, 10.4%, 77% boys) and IN (N = 158, 17.3%, 75% boys). Controlling for sex, each unit of HI PRs starting at 1.5 years and at 4 years for IN, increased more than 2-fold the risk of belonging to the high trajectory, with incremental contributions (Odds Ratios = 2.5-4.5) at subsequent ages. Quantitative genetic analyses showed that genetic contributions underlying preschool PRs accounted for up to a quarter and a third of the heritability of later HI and IN, respectively. Genes underlying 1.5-year HI and 4-year IN contributed to 6 of 8 later HI and IN time-points and largely explained the corresponding phenotypic correlations. Results provide phenotypic and genetic evidence that preschool parent reports of HI and IN are valid means to predict developmental risk of ADHD.

UK IBD Twin Registry: Concordance and Environmental Risk Factors of Twins with IBD

BACKGROUND AND AIMS

Twin studies have long been used to infer heritability. Within the 'omics era, twin cohorts have even greater research potential. This study describes the formation of the UK IBD Twin Registry and analysis of concordance and environmental factors.

METHOD

Twin pairs with IBD were recruited by advertising via IBD charities and social media, re-tracing a dormant IBD database and clinician referral. Details of zygosity, concordance, disease history and environmental factors were assessed. Pair concordance was calculated, and environmental factors were analysed with logistic regression models adjusted for zygosity and concordance.

RESULTS

Ninety-one twin pairs were included in the analysis; forty-two with CD and forty-nine with UC. More MZ twin pairs with CD were concordant compared with DZ pairs, thus inferring heritability (Chi-sq. 15.6. P < 0.001). In UC, MZ concordance was also numerically greater. Cigarette smoking was predictive of CD (OR 2.66, 95% CI 1.16 to 6.07 P = 0.02); there may be an independent association with cannabis smoking (OR 2.59 95% CI 0.89 to 7.55 P = 0.08). Breastfeeding was protective against UC (OR 0.48, 95% CI 0.25-0.93, P = 0.03), but not CD. Self-reports of less occurrences of gastroenteritis than peers were protective against future UC onset (OR 0.33 95% CI 0.15 to 0.74, P = 0.01). Method of delivery, parental attitudes towards hygiene and recall of diet did not impact future IBD concordance.

CONCLUSIONS

This study supports the heritability of IBD. Twin study analysis was able to elucidate environmental factors associated with IBD.

Differences between germline genomes of monozygotic twins

Despite the important role that monozygotic twins have played in genetics research, little is known about their genomic differences. Here we show that monozygotic twins differ on average by 5.2 early developmental mutations and that approximately 15% of monozygotic twins have a substantial number of these early developmental mutations specific to one of them. Using the parents and offspring of twins, we identified pre-twinning mutations. We observed instances where a twin was formed from a single cell lineage in the pre-twinning cell mass and instances where a twin was formed from several cell lineages. CpG>TpG mutations increased in frequency with embryonic development, coinciding with an increase in DNA methylation. Our results indicate that allocations of cells during development shapes genomic differences between monozygotic twins.

Phenotypic differences and similarities of monozygotic twins with maturity-onset diabetes of the young type 5

Here, we report phenotypic differences and similarities of monozygotic twins with maturity-onset diabetes of the young type 5 harboring a partial deletion of chromosome 17q12. The proband and her twin sister manifested complete aplasia and marked hypoplasia, respectively, of the body and tail of the pancreas. Whereas both twins showed marked hypoplasia of the right kidney and multiple cysts in both kidneys, only the proband's sister showed hydronephrosis in the left kidney. The proband had profound defects in insulin and glucagon secretion, as well as mild renal dysfunction, whereas her sister had pronounced renal dysfunction accompanied by mild defects in insulin and glucagon secretion. Both twins manifested hypomagnesemia and hyperuricemia, but no apparent liver dysfunction or intellectual disability. The severity of renal and pancreatic defects differed between monozygotic twins with maturity-onset diabetes of the young type 5, suggesting that the phenotypes of this condition are determined not solely by genetic factors.

De novo single-nucleotide and copy number variation in discordant monozygotic twins reveals disease-related genes

Recent studies have demonstrated genetic differences between monozygotic (MZ) twins. To test the hypothesis that early post-twinning mutational events associate with phenotypic discordance, we investigated a cohort of 13 twin pairs (n = 26) discordant for various clinical phenotypes using whole-exome sequencing and screened for copy number variation (CNV). We identified a de novo variant in PLCB1, a gene involved in the hydrolysis of lipid phosphorus in milk from dairy cows, associated with lactase non-persistence, and a variant in the mitochondrial complex I gene MT-ND5 associated with amyotrophic lateral sclerosis (ALS). We also found somatic variants in multiple genes (TMEM225B, KBTBD3, TUBGCP4, TFIP11) in another MZ twin pair discordant for ALS. Based on the assumption that discordance between twins could be explained by a common variant with variable penetrance or expressivity, we screened the twin samples for known pathogenic variants that are shared and identified a rare deletion overlapping ARHGAP11B, in the twin pair manifesting with either schizotypal personality disorder or schizophrenia. Parent-offspring trio analysis was implemented for two twin pairs to assess potential association of variants of parental origin with susceptibility to disease. We identified a de novo variant in RASD2 shared by 8-year-old male twins with a suspected diagnosis of autism spectrum disorder (ASD) manifesting as different traits. A de novo CNV duplication was also identified in these twins overlapping CD38, a gene previously implicated in ASD. In twins discordant for Tourette's syndrome, a paternally inherited stop loss variant was detected in AADAC, a known candidate gene for the disorder.

Copy Number Variants and Exome Sequencing Analysis in Six Pairs of Chinese Monozygotic Twins Discordant for Congenital Heart Disease

Congenital heart disease (CHD) is one of the most common birth defects. More than 200 susceptibility loci have been identified for CHDs, yet a large part of the genetic risk factors remain unexplained. Monozygotic (MZ) twins are thought to be completely genetically identical; however, discordant phenotypes have been found in MZ twins. Recent studies have demonstrated genetic differences between MZ twins. We aimed to test whether copy number variants (CNVs) and/or genetic mutation differences play a role in the etiology of CHDs by using single nucleotide polymorphism (SNP) genotyping arrays and whole exome sequencing of twin pairs discordant for CHDs. Our goal was to identify mutations present only in the affected twins, which could identify novel candidates for CHD susceptibility loci. We present a comprehensive analysis for the CNVs and genetic mutation results of the selected individuals but detected no consistent differences within the twin pairs. Our study confirms that chromosomal structure or genetic mutation differences do not seem to play a role in the MZ twins discordant for CHD.

Genetic instrumental variable regression: Explaining socioeconomic and health outcomes in nonexperimental data

We propose genetic instrumental variable (GIV) regression—a method that controls for pleiotropic effects of genes on two variables. GIV regression is broadly applicable to study outcomes for which polygenic scores from large-scale genome-wide association studies are available. We explore the performance of GIV regression in the presence of pleiotropy across a range of scenarios and find that it yields more accurate estimates than alternative approaches such as ordinary least-squares regression or Mendelian randomization. When GIV regression is combined with proper controls for purely environmental sources of bias (e.g., using control variables and sibling fixed effects), it improves our understanding of the causal relationships between genetically correlated variables.

Identifying causal effects in nonexperimental data is an enduring challenge. One proposed solution that recently gained popularity is the idea to use genes as instrumental variables [i.e., Mendelian randomization (MR)]. However, this approach is problematic because many variables of interest are genetically correlated, which implies the possibility that many genes could affect both the exposure and the outcome directly or via unobserved confounding factors. Thus, pleiotropic effects of genes are themselves a source of bias in nonexperimental data that would also undermine the ability of MR to correct for endogeneity bias from nongenetic sources. Here, we propose an alternative approach, genetic instrumental variable (GIV) regression, that provides estimates for the effect of an exposure on an outcome in the presence of pleiotropy. As a valuable byproduct, GIV regression also provides accurate estimates of the chip heritability of the outcome variable. GIV regression uses polygenic scores (PGSs) for the outcome of interest which can be constructed from genome-wide association study (GWAS) results. By splitting the GWAS sample for the outcome into nonoverlapping subsamples, we obtain multiple indicators of the outcome PGSs that can be used as instruments for each other and, in combination with other methods such as sibling fixed effects, can address endogeneity bias from both pleiotropy and the environment. In two empirical applications, we demonstrate that our approach produces reasonable estimates of the chip heritability of educational attainment (EA) and show that standard regression and MR provide upwardly biased estimates of the effect of body height on EA.

Copy Number Variation Analysis of 100 Twin Pairs Enriched for Neurodevelopmental Disorders

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.

Post-zygotic genomic changes in glutamate and dopamine pathway genes may explain discordance of monozygotic twins for schizophrenia

Background

Monozygotic twins are valuable in assessing the genetic vs environmental contribution to diseases. In the era of complete genome sequences, they allow identification of mutational mechanisms and specific genes and pathways that offer predisposition to the development of complex diseases including schizophrenia.

Methods

We sequenced the complete genomes of two pairs of monozygotic twins discordant for schizophrenia (MZD), including one representing a family tetrad. The family specific complete sequences have allowed identification of post zygotic mutations between MZD genomes. It allows identification of affected genes including relevant network and pathways that may account for the diseased state in pair specific patient.

Results

We found multiple twin specific sequence differences between co-twins that included small nucleotides [single nucleotide variants (SNV), small indels and block substitutions], copy number variations (CNVs) and structural variations. The genes affected by these changes belonged to a number of canonical pathways, the most prominent ones are implicated in schizophrenia and related disorders. Although these changes were found in both twins, they were more frequent in the affected twin in both pairs. Two specific pathway defects, glutamate receptor signaling and dopamine feedback in cAMP signaling pathways, were uniquely affected in the two patients representing two unrelated families.

Conclusions

We have identified genome-wide post zygotic mutations in two MZD pairs affected with schizophrenia. It has allowed us to use the threshold model and propose the most likely cause of this disease in the two patients studied. The results support the proposition that each schizophrenia patient may be unique and heterogeneous somatic de novo events may contribute to schizophrenia threshold and discordance of the disease in monozygotic twins.

Electronic supplementary material

The online version of this article (10.1186/s40169-017-0174-1) contains supplementary material, which is available to authorized users.

Establishment of the Avera Twin Register in the Midwest USA

The Avera Twin Register (ATR) aims to study environmental and genetic influences on health and disease using a longitudinal repository of biological specimens, survey data, and health information provided by multiples and their family members. The ATR is located in Sioux Falls, South Dakota, which is a rural and frontier area in the Midwestern United States with a density of four people per square kilometer. The target area of the ATR is South Dakota and the four surrounding states: Minnesota, Iowa, North Dakota, and Nebraska. Enrollment of twins and higher-order multiples of all ages and their family members started on May 18, 2016. A description of the first 13 months of enrollment in this longitudinal register will be provided. The ATR will collect longitudinal data on lifestyle, including diet and activity levels, aging, complex traits, and diseases. Upon registration, all participants are genotyped on the Illumina Global Screening Array (GSA) and twins and higher order multiples receive information on their zygosity. The ATR aims to contribute to large international GWAS consortia and collaborates closely with the Netherlands Twin Register, allowing for the comparison of collected data and analyses of results. In addition, the ATR will address twin-specific questions.

The etiology of autistic traits in preschoolers: a population-based twin study

BACKGROUND

Autism Spectrum Disorders (ASD) are highly heritable, but the exact etiological mechanisms underlying the condition are still unclear.

METHODS

Using a multiple rater twin design in a large sample of general population preschool twins, this study aimed to (a) estimate the contribution of genetic and environmental factors to autistic traits, controlling for the possible effects of rater bias, (b) to explore possible sex differences in etiology and (c) to investigate the discordance in autistic traits in monozygotic and same-sex dizygotic twin pairs. The Netherlands Twin Register collected maternal and paternal ratings on autistic traits from a general population of 38,798 three-year-old twins. Autistic traits were assessed with the DSM-oriented Pervasive Developmental Problems scale of the Child Behavior Check List for preschoolers (1½-5 years).

RESULTS

Mother and fathers showed high agreement in their assessment of autistic traits (r = .60-.66). Differences between children in autistic traits were largely accounted for by genetic effects (boys: 78% and girls: 83%). Environmental effects that are unique to a child also played a modest role. Environmental effects shared by children growing up in the same family were negligible, once rater bias was controlled for. While the prevalence for clinical ASD is higher in boys than in girls, this study did not find evidence for striking differences in the etiology of autistic traits across the sexes. Even though the heritability was high, 29% of MZ twin pairs were discordant for high autistic traits (clinical range vs. normal development), suggesting that despite high genetic risk, environmental factors might lead to resilience, unaffected status in the context of genetic risk, in some children.

CONCLUSIONS

It is important to focus future research on risk factors that might interplay with a genetic disposition for ASD, but also on protective factors that make a difference in the lives of children at genetic risk.

A Genome-Wide Association Study to Identify Potential Germline Copy Number Variants for Sporadic Breast Cancer Susceptibility

Breast cancer (BC) predisposition in populations arises from both genetic and nongenetic risk factors. Structural variations such as copy number variations (CNVs) are heritable determinants for disease susceptibility. The primary objectives of this study are (1) to identify CNVs associated with sporadic BC using a genome-wide association study (GWAS) design; (2) to utilize 2 distinct CNV calling algorithms to identify concordant CNVs as a strategy to reduce false positive associations in the hypothesis-generating GWAS discovery phase, and (3) to identify potential candidate CNVs for follow-up replication studies. We used Affymetrix SNP Array 6.0 data profiled on Caucasian subjects (422 cases/348 controls) to call CNVs using algorithms implemented in Nexus Copy Number and Partek Genomics Suite software. Nexus algorithm identified CNVs associated with BC (731 autosomal CNVs with >5% frequency in the total sample and Q < 0.05). Thirteen CNVs were identified when Partek algorithm-called CNVs were overlapped with Nexus-identified CNVs; these CNVs showed concordances for frequency, effect size, and direction. Coding genes present within BC-associated CNVs were known to play a role in disease etiology and prognosis. Long noncoding RNAs identified within CNVs showed tissue-specific expression, indicating potential functional relevance of the findings. The identified candidate CNVs warrant independent replication.

A review of the mechanisms and evidence for typical and atypical twinning

The mechanisms responsible for twinning and disorders of twin gestations have been the subject of considerable interest by physicians and scientists, and cases of atypical twinning have called for a reexamination of the fundamental theories invoked to explain twin gestations. This article presents a review of the literature focusing on twinning and atypical twinning with an emphasis on the phenomena of chimeric twins, phenotypically discordant monozygotic twins, mirror-image twins, polar body twins, complete hydatidiform mole with a coexistent twin, vanishing twins, fetus papyraceus, fetus in fetu, superfetation, and superfecundation. The traditional models attributing monozygotic twinning to a fission event, and more recent models describing monozygotic twinning as a fusion event, are critically reviewed. Ethical restrictions on scientific experimentation with human embryos and the rarity of cases of atypical twinning have limited opportunities to elucidate the exact mechanisms by which these phenomena occur. Refinements in the modeling of early embryonic development in twin pregnancies may have significant clinical implications. The article includes a series of figures to illustrate the phenomena described.

Whole-Exome Sequencing in Nine Monozygotic Discordant Twins

By definition, monozygotic (MZ) twins carry an identical set of genetic information. The observation of early post-twinning mutational events was shown to cause phenotypic discordance among MZ twin pairs. These mutational events comprise genomic alterations at different scales, ranging from single nucleotide changes to larger copy-number variations (CNVs) of varying sizes, as well as epigenetic changes. Here, we performed whole-exome sequencing (WES) in nine discordant MZ twins to identify somatic mutational events in the affected twin that might exert a dominant negative effect. Five of these MZ twin pairs were discordant for congenital heart defects (CHD), two for endocrine disorders, one for omphalocele, and one for congenital diaphragmatic hernia (CDH). Analysis of WES data from all nine MZ twin pairs using the de novo probability tool DeNovoGear detected only one apparent de novo variation in TMPRSS13 in one of the CHD-affected twins. Analysis of WES data from all nine MZ twin pairs by using standard filter criteria without the de novo probability tool DeNovoGear revealed a total of 6,657 variations in which both the twin pairs differed. After filtering for variations only present in the affected twins and absent in in-house controls, 722 variations remained. Visual inspection for read quality decreased this number to 12, present only in the affected twin. However, Sanger sequencing of the overall 13 variations failed to confirm the variation in the affected twin. These results suggest that somatic mutational events in coding regions do not seem to play a major role in the phenotypic expression of MZ discordant twin pairs.

Global patterns of large copy number variations in the human genome reveal complexity in chromosome organization

Global patterns of copy number variations (CNVs) in chromosomes are required to understand the dynamics of genome organization and complexity. For this study, analysis was performed using the Affymetrix Genome-Wide Human SNP Array 6.0 chip and CytoScan High-Density arrays. We identified a total of 44 109 CNVs from 1715 genomes with a mean of 25 CNVs in an individual, which established the first drafts of population-specific CNV maps providing a rationale for prioritizing chromosomal regions. About 19 905 ancient CNVs were identified across all chromosomes and populations at varying frequencies. CNV count, and sometimes CNV size, contributed to the bulk CNV size of the chromosome. Population specific lengthening and shortening of chromosomal length was observed. Sex bias for CNV presence was largely dependent on ethnicity. Lower CNV inheritance rate was observed for India, compared to YRI and CEU. A total of 33 candidate CNV hotspots from 5382 copy number (CN) variable region (CNVR) clusters were identified. Population specific CNV distribution patterns in p and q arms disturbed the assumption that CNV counts in the p arm are less common compared to long arms, and the CNV occurrence and distribution in chromosomes is length independent. This study unraveled the force of independent evolutionary dynamics on genome organization and complexity across chromosomes and populations.

Heritability in inflammatory bowel disease: from the first twin study to genome-wide association studies

Since Tysk et al's pioneering analysis of the Swedish twin registry, twin and family studies continue to support a strong genetic basis of the inflammatory bowel diseases. The coefficient of heritability for siblings of inflammatory bowel disease probands is 25 to 42 for Crohn's disease and 4 to 15 for ulcerative colitis. Heritability estimates for Crohn's disease and ulcerative colitis from pooled twin studies are 0.75 and 0.67, respectively. However, this is at odds with the much lower heritability estimates from Genome-Wide Association Studies (GWAS). This "missing heritability" is likely due to shortfalls in both family studies and GWAS. The coefficient of heritability fails to account for familial shared environment. Heritability calculations from twin data are based on Falconer's method, with premises that are increasingly understood to be flawed. GWAS based heritability estimates may underestimate heritability due to incomplete linkage disequilibrium, and because some single nucleotide polypeptides (SNPs) do not reach a level of significance to allow detection. SNPs missed by GWAS include common SNPs with low penetrance and rare SNPs with high penetrance. All methods of heritability estimation regard genetic and environmental variance as separate entities, although it is now understood that there is a complex multidirectional interplay between genetic are environmental factors mediated by the microbiota, the epigenome, and the innate and acquired immune systems. Due to the limitations of heritability estimates, it is unlikely that a true value for heritability will be reached. Further work aimed at quantifying the variance explained across GWAS, epigenome-wide, and microbiota-wide association studies will help to define factors leading to inflammatory bowel disease.

Genome-wide blood DNA methylation alterations at regulatory elements and heterochromatic regions in monozygotic twins discordant for obesity and liver fat

Background

The current epidemic of obesity and associated diseases calls for swift actions to better understand the mechanisms by which genetics and environmental factors affect metabolic health in humans. Monozygotic (MZ) twin pairs showing discordance for obesity suggest that epigenetic influences represent one such mechanism. We studied genome-wide leukocyte DNA methylation variation in 30 clinically healthy young adult MZ twin pairs discordant for body mass index (BMI; average within-pair BMI difference: 5.4 ± 2.0 kg/m²).

Results

There were no differentially methylated cytosine-guanine (CpG) sites between the co-twins discordant for BMI. However, stratification of the twin pairs based on the level of liver fat accumulation revealed two epigenetically highly different groups. Significant DNA methylation differences (n = 1,236 CpG sites (CpGs)) between the co-twins were only observed if the heavier co-twins had excessive liver fat (n = 13 twin pairs). This unhealthy pattern of obesity was coupled with insulin resistance and low-grade inflammation. The differentially methylated CpGs included 23 genes known to be associated with obesity, liver fat, type 2 diabetes mellitus (T2DM) and metabolic syndrome, and potential novel metabolic genes. Differentially methylated CpG sites were overrepresented at promoters, insulators, and heterochromatic and repressed regions. Based on predictions by overlapping histone marks, repressed and weakly transcribed sites were significantly more often hypomethylated, whereas sites with strong enhancers and active promoters were hypermethylated. Further, significant clustering of differentially methylated genes in vitamin, amino acid, fatty acid, sulfur, and renin-angiotensin metabolism pathways was observed.

Conclusions

The methylome in leukocytes is altered in obesity associated with metabolic disturbances, and our findings indicate several novel candidate genes and pathways in obesity and obesity-related complications.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-015-0073-5) contains supplementary material, which is available to authorized users.

Large Autosomal Copy-Number Differences within Unselected Monozygotic Twin Pairs are Rare

Monozygotic (MZ) twins form an important system for the study of biological plasticity in humans. While MZ twins are generally considered to be genetically identical, a number of studies have emerged that have demonstrated copy-number differences within a twin pair, particularly in those discordant for disease. The rate of autosomal copy-number variation (CNV) discordance within MZ twin pairs was investigated using a population sample of 376 twin pairs genotyped on Illumina Human610-Quad arrays. After CNV calling using both QuantiSNP and PennCNV followed by manual annotation, only a single CNV difference was observed within the MZ twin pairs, being a 130 KB duplication of chromosome 5. Five other potential discordant CNV were called by the software, but excluded based on manual annotation of the regions. It is concluded that large CNV discordance is rare within MZ twin pairs, indicating that any CNV difference found within phenotypically discordant MZ twin pairs has a high probability of containing the causal gene(s) involved.

CNV Concordance in 1,097 MZ Twin Pairs

Monozygotic (MZ) twins are genetically identical at conception, making them informative subjects for studies on somatic mutations. Copy number variants (CNVs) are responsible for a substantial part of genetic variation, have relatively high mutation rates, and are likely to be involved in phenotypic variation. We conducted a genome-wide survey for post-twinning de novo CNVs in 1,097 MZ twin pairs. Comparisons between MZ twins were made by CNVs measured in DNA from blood or buccal epithelium with the Affymetrix 6.0 microarray and two calling algorithms. In addition, CNV concordance rates were compared between the different sources of DNA, and gene-enrichment association analyses were conducted for thought problems (TP) and attention problems (AP) using CNVs concordant within MZ pairs. We found a total of 153 putative post-twinning de novo CNVs >100 kb, of which the majority resided in 15q11.2. Based on the discordance of raw intensity signals a selection was made of 20 de novo CNVs for a qPCR validation experiments. Two out of 20 post-twinning de novo CNVs were validated with qPCR in the same twin pair. The 13-year-old MZ twin pair that showed two discordances in CN in 15q11.2 in their buccal DNA did not show large phenotypic differences. From the remaining 18 putative de novo CNVs, 17 were deletions or duplications that were concordant within MZ twin pairs. Concordance rates within twin pairs of CNV calls with CN ≠ 2 were ~80%. Buccal epithelium-derived DNA showed a slightly but significantly higher concordance rate, and blood-derived DNA showed significantly more concordant CNVs per twin pair. The gene-enrichment analyses on concordant CNVs showed no significant associations between CNVs overlapping with genes involved in neuronal processes and TP or AP after accounting for the source of DNA.

Attention-Deficit/Hyperactivity Disorder: Focus upon Aberrant N-Methyl-D-Aspartate Receptors Systems

Attention-deficit/hyperactivity disorder (ADHD) pathophysiology persists in an obscure manner with complex interactions between symptoms, staging, interventions, genes, and environments. Only on the basis of increasing incidence of the disorder, the need for understanding is greater than ever. The notion of an imbalance between central inhibitory/excitatory neurotransmitters is considered to exert an essential role. In this chapter, we first review how the default mode network functions and dysfunction in individuals diagnosed with ADHD. We also present and briefly review some of the animal models used to examine the neurobiological aspects of ADHD. There is much evidence indicating that compounds/interventions that antagonize/block glutamic acid receptors and/or block the glutamate signal during the "brain growth spurt" or in the adult animal may induce functional and biomarker deficits. Additionally, we present evidence suggesting that animals treated with glutamate blockers at the period of the "brain growth spurt" fail to perform the exploratory activity, observed invariably with control mice, that is associated with introduction to a novel environment (the test cages). Later, when the control animals show less locomotor and rearing activity, i.e., interest in the test cages, the MK-801, ketamine and ethanol treated mice showed successively greater levels of locomotion and rearing (interest), i.e., they fail to "habituate" effectively, implying a cognitive dysfunction. These disturbances of glutamate signaling during a critical period of brain development may contribute to the ADHD pathophysiology. As a final addition, we have briefly identified new research venues in the interaction between ADHD, molecular studies, and personality research.

Copy number variation burden on asthma subgenome in normal cohorts identifies susceptibility markers

PURPOSE

Asthma is a complex disease caused by interplay of genes and environment on the genome of an individual. Copy number variations (CNVs) are more common compared to the other variations that disrupt genome organization. The effect of CNVs on asthma subgenome has been less studied compared to studies on the other variations. We report the assessments of CNV burden in asthma genes of normal cohorts carried out in different geographical areas of the world and discuss the relevance of the observation with respect to asthma pathogenesis.

METHODS

CNV analysis was performed using Affymerix high-resolution arrays, and various bioinformatics tools were used to understand the influence of genes on asthma pathogenesis.

RESULTS

This study identified 61 genes associated with asthma and provided various mechanisms and pathways underlying asthma pathogenesis. CCL3L1, ADAM8, and MUC5B were the most prevalent asthma genes. Among them, CCL3L1 was found across all 12 populations in varying copy number states. This study also identified the inheritance of asthma-CNVs from parents to offspring creating the latent period for manifestation of asthma.

CONCLUSIONS

This study revealed CNV burden with varying copy number states and identified susceptibility towards the disease manifestation. It can be hypothesized that primary CNVs may not be the initiating event in the pathogenesis of asthma and additional preceding mutations or CNVs may be required. The initiator or primary CNVs sensitize normal cohorts leading to an increased probability of accumulating mutations or exposure to allergic stimulating agents that can augment the development of asthma.

Genomic screening of testicular germ cell tumors from monozygotic twins

BACKGROUND

Testicular germ cell tumors (TGCTs) account for 1-2% of all tumors in young and middle aged men. A 75-fold increase in TCGT development has been reported for monozygotic (MZ) twins. Therefore, the occurrence of simultaneous tumors in MZ twins emphasizes the importance of genetic factors that influence the risk of developing these tumors. Genomic screening was performed for one family containing MZ twins with testicular germ cell tumors, in order to define alterations associated with risk of tumor development.

METHODS

Copy number alterations were evaluated using array-CGH (4x44K, Agilent Technologies) in one seminoma and one embryonal carcinoma (EC) from MZ twins. In addition, genomic alterations from the tumors and peripheral blood cells of the twins were compared to the parental genomes via their peripheral blood cells.

RESULTS

Embryonal carcinoma (Twin-1 t) presented a lower frequency of genomic alterations compared to the seminoma (Twin-2 t). One minimal common region of loss was observed in 9p13.1-p12 in the comparison between DNA from blood samples for Twin-1 and Twin-2. In this region is mapped the CNTNAP3 gene which was confirmed as involved in losses by qPCR. Comparative analysis of novel CNVs between the Twin-1 t and Twin-2 t showed five minimal common regions involving gain at chromosomes 12 (12p12.3-p11.1 and 12p13.33-p12.3), while losses were observed at 10p15.3-p15.2, 13q21.1-q21.2 and 15q11.1-q11.2. In addition, one exclusive rare copy number alteration was detected in Twin-1 t and Twin-2 t, and 19 novel alterations were identified in the Twin-2 t.

CONCLUSION

Distinct genomic profiles for MZ twins with phenotypically different TGCT were described. Of particular interest, 12p gains were detected exclusively in tumor samples. In peripheral blood samples, loss of 9p13.1-p12 was the unique novel CNV shared by the twins, confirming the involvement of CNTNAP3 gene in TGCTs development. Although similar CNV profiles were shared by both the peripheral blood and tumor samples of the twins, tumor-specific CNV loci were identified for seminoma and non-seminomatous tumors. These findings suggest the presence of de novo germline structural alterations and TGCT predisposition.

Behavior genetics: past, present, future

The disciplines of developmental psychopathology and behavior genetics are concerned with many of the same questions about the etiology and course of normal and abnormal behavior and about the factors that promote typical development despite the presence of risk. The goal of this paper is to summarize how research in behavior genetics has shed light on questions that are central to developmental psychopathology. We briefly review the origins of behavior genetics, summarize the findings that have been gleaned from several decades of quantitative and molecular genetics research, and describe future directions for research that will delineate gene function as well as pathways from genes to brain to behavior. The importance of environmental contributions, at both genetic and epigenetic levels, will be discussed. We conclude that behavior genetics has made significant contributions to developmental psychopathology by documenting the interplay among risk and protective factors at multiple levels of the organism, by clarifying the causal status of risk exposures, and by identifying factors that account for change and stability in psychopathology. As the tools to identify gene function become increasingly sophisticated, and as behavioral geneticists become increasingly interdisciplinary in their scope, the field is poised to make ever greater contributions to our understanding of typical and atypical development.

Copy number variations burden on miRNA genes reveals layers of complexities involved in the regulation of pathways and phenotypic expression

MicroRNAs are involved in post-transcriptional down-regulation of gene expression. Variations in miRNA genes can severely affect downstream-regulated genes and their pathways. However, population-specific burden of CNVs on miRNA genes and the complexities created towards the phenotype is not known. From a total of 44109 CNVs investigated from 1715 individuals across 12 populations using high-throughput arrays, 4007 miRNA-CNVs (∼ 9%) consisting 6542 (∼ 5%) miRNA genes with a total of 333 (∼ 5%) singleton miRNA genes were identified. We found miRNA-CNVs across the genomes of individuals showing multiple hits in many targets, co-regulated under the same pathway. This study proposes four mechanisms unraveling the many complexities in miRNA genes, targets and co-regulated miRNA genes towards establishment of phenotypic diversity.

Copy number variation distribution in six monozygotic twin pairs discordant for schizophrenia

We have evaluated copy number variants (CNVs) in six monozygotic twin pairs discordant for schizophrenia. The data from Affymetrix® Human SNP 6.0 arrays™ were analyzed using Affymetrix® Genotyping Console™, Partek® Genomics Suite™, PennCNV, and Golden Helix SVS™. This yielded both program-specific and overlapping results. Only CNVs called by Affymetrix Genotyping Console, Partek Genomics Suite, and PennCNV were used in further analysis. This analysis included an assessment of calls in each of the six twin pairs towards identification of unique CNVs in affected and unaffected co-twins. Real time polymerase chain reaction (PCR) experiments confirmed one CNV loss at 7q11.21 that was found in the affected patient but not in the unaffected twin. The results identified CNVs and genes that were previously implicated in mental abnormalities in four of the six twin pairs. It included PYY (twin pairs 1 and 5), EPHA3 (twin pair 3), KIAA1211L (twin pair 4), and GPR139 (twin pair 5). They represent likely candidate genes and CNVs for the discordance of four of the six monozygotic twin pairs for this heterogeneous neurodevelopmental disorder. An explanation for these differences is ontogenetic de novo events that differentiate in the monozygotic twins during development.

Absence of substantial copy number differences in a pair of monozygotic twins discordant for features of autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable disease (~0.9) with a complex genetic etiology. It is initially characterized by altered cognitive ability which commonly includes impaired language and communication skills as well as fundamental deficits in social interaction. Despite the large amount of studies described so far, the high clinical diversity affecting the autism phenotype remains poorly explained. Recent studies suggest that rare genomic variations, in particular copy number variation (CNV), may account for a significant proportion of the genetic basis of ASD. The use of disease-discordant monozygotic twins represents a powerful strategy to identify de novo and inherited CNV in the disorder. Here we present the results of a comparative genome hybridization (CGH) analysis with a pair of monozygotic twins affected of ASD with significant differences in their clinical manifestations that specially affect speech language impairment and communication skills. Array CGH was performed in three different tissues: blood, saliva, and hair follicle, in an attempt to identify germinal and somatic CNV regions that may explain these differences. Our results argue against a role of large CNV rearrangements as a molecular etiology of the observed differences. This forwards future research to explore de novo point mutation and epigenomic alterations as potential explanations of the observed clinical differences.

De novo copy number variations in cloned dogs from the same nuclear donor

Background

Somatic mosaicism of copy number variants (CNVs) in human body organs and de novo CNV event in monozygotic twins suggest that de novo CNVs can occur during mitotic recombination. These de novo CNV events are important for understanding genetic background of evolution and diverse phenotypes. In this study, we explored de novo CNV event in cloned dogs with identical genetic background.

Results

We analyzed CNVs in seven cloned dogs using the nuclear donor genome as reference by array-CGH, and identified five de novo CNVs in two of the seven clones. Genomic qPCR, dye-swap array-CGH analysis and B-allele profile analysis were used for their validation. Two larger de novo CNVs (5.2 Mb and 338 Kb) on chromosomes X and 19 in clone-3 were consistently validated by all three experiments. The other three smaller CNVs (sized from 36.1 to76.4 Kb) on chromosomes 2, 15 and 32 in clone-3 and clone-6 were verified by at least one of the three validations. In addition to the de novo CNVs, we identified a 37 Mb-sized copy neutral de novo loss of heterozygosity event on chromosome 2 in clone-6.

Conclusions

To our knowledge, this is the first report of de novo CNVs in the cloned dogs which were generated by somatic cell nuclear transfer technology. To study de novo genetic events in cloned animals can help understand formation mechanisms of genetic variants and their biological implications.

Theoretical framework of population genetics with somatic mutations taken into account: application to copy number variations in humans

Traditionally, population genetics focuses on the dynamics of frequencies of alleles acquired by mutations on germ-lines, because only such mutations are heritable. Typical genotyping experiments, however, use DNA from some somatic tissues such as blood, which harbors somatic mutations at the current generation in addition to germ-line mutations accumulated since the most recent common ancestor of the sample. This common practice may sometimes cause erroneous interpretations of polymorphism data, unless we properly understand the role of somatic mutations in population genetics. We here introduce a very basic theoretical framework of population genetics with somatic mutations taken into account. It is easy to imagine that somatic mutations at the current generation simply add individual-specific variations, as errors in mutation detection do. Our theory quantifies this increment under various conditions. We find that the major contribution of somatic mutations plus errors is to very rare variants, particularly to singletons. The relative contribution is markedly large when mutations are deleterious. Because negative selection also increases rare variants, it is important to distinguish the roles of these mutually confounding factors when we interpret the data, even after correcting for demography. We apply this theory to human copy number variations (CNVs), for which the composite effect of somatic mutations and errors may not be negligible. Using genome-wide CNV data, we demonstrate how the joint action of the two factors, selection and somatic mutations plus errors, shapes the observed pattern of polymorphism.

Twin Mitochondrial Sequence Analysis

When applying genome-wide sequencing technologies to disease investigation, it is increasingly important to resolve sequence variation in regions of the genome that may have homologous sequences. The human mitochondrial genome challenges interpretation given the potential for heteroplasmy, somatic variation, and homologous nuclear mitochondrial sequences (numts). Identical twins share the same mitochondrial DNA (mtDNA) from early life, but whether the mitochondrial sequence remains similar is unclear. We compared an adult monozygotic twin pair using high-throughput sequencing and evaluated variants with primer extension and mitochondrial preenrichment. Thirty-seven variants were shared between the twin individuals, and the variants were verified on the original genomic DNA. These studies support highly identical genetic sequence in this case. Certain low-level variant calls were of high quality and homology to the mtDNA, and they were further evaluated. When we assessed calls in preenriched mtDNA templates, we found that these may represent numts, which can be differentiated from mtDNA variation. We conclude that twin identity extends to mtDNA, and it is critical to differentiate between numts and mtDNA in genome sequencing, particularly as significant heteroplasmy could influence genome interpretation. Further studies on mtDNA and numts will aid in understanding how variation occurs and persists.

Comparison of Genomic and Epigenomic Expression in Monozygotic Twins Discordant for Rett Syndrome

Monozygotic (identical) twins have been widely used in genetic studies to determine the relative contributions of heredity and the environment in human diseases. Discordance in disease manifestation between affected monozygotic twins has been attributed to either environmental factors or different patterns of X chromosome inactivation (XCI). However, recent studies have identified genetic and epigenetic differences between monozygotic twins, thereby challenging the accepted experimental model for distinguishing the effects of nature and nurture. Here, we report the genomic and epigenomic sequences in skin fibroblasts of a discordant monozygotic twin pair with Rett syndrome, an X-linked neurodevelopmental disorder characterized by autistic features, epileptic seizures, gait ataxia and stereotypical hand movements. The twins shared the same de novo mutation in exon 4 of the MECP2 gene (G269AfsX288), which was paternal in origin and occurred during spermatogenesis. The XCI patterns in the twins did not differ in lymphocytes, skin fibroblasts, and hair cells (which originate from ectoderm as does neuronal tissue). No reproducible differences were detected between the twins in single nucleotide polymorphisms (SNPs), insertion-deletion polymorphisms (indels), or copy number variations. Differences in DNA methylation between the twins were detected in fibroblasts in the upstream regions of genes involved in brain function and skeletal tissues such as Mohawk Homeobox (MKX), Brain-type Creatine Kinase (CKB), and FYN Tyrosine Kinase Protooncogene (FYN). The level of methylation in these upstream regions was inversely correlated with the level of gene expression. Thus, differences in DNA methylation patterns likely underlie the discordance in Rett phenotypes between the twins.

Distribution of disease-associated copy number variants across distinct disorders of cognitive development

OBJECTIVE

The purpose of the present study was to discover the extent to which distinct DSM disorders share large, highly recurrent copy number variants (CNVs) as susceptibility factors. We also sought to identify gene mechanisms common to groups of diagnoses and/or specific to a given diagnosis based on associations with CNVs.

METHOD

Systematic review of 820 PubMed articles on autism spectrum disorder (ASD), intellectual disability (ID), schizophrenia, and epilepsy produced 54 CNVs associated with one or several disorders. Pathway analysis on genes implicated by CNVs in different groupings was conducted.

RESULTS

The majority of CNVs were found in ID with the other disorders somewhat subsumed, yet certain CNVs were associated with isolated or groups of disorders. Based on genes implicated by CNVs, ID encompassed 96.8% of genes in ASD, 92.8% of genes in schizophrenia, and 100.0% of genes in epilepsy. Pathway analysis revealed that synapse processes were enriched in ASD, ID, and schizophrenia. Disease-specific processes were identified in ID (actin cytoskeleton processes), schizophrenia (ubiquitin-related processes), and ASD (synaptic vesicle transport and exocytosis).

CONCLUSIONS

Intellectual disability may arise from the broadest range of genetic pathways, and specific subsets of these pathways appear to be relevant to other disorders or combinations of these disorders. It is clear that statistically significant CNVs across disorders of cognitive development are highly enriched for biological processes related to the synapse. There are also disorder-specific processes that may aid in understanding the distinct presentations and pathophysiology of these disorders.

The Adult Netherlands Twin Register: twenty-five years of survey and biological data collection

Over the past 25 years, the Adult Netherlands Twin Register (ANTR) has collected a wealth of information on physical and mental health, lifestyle, and personality in adolescents and adults. This article provides an overview of the sources of information available, the main research findings, and an outlook for the future. Between 1991 and 2012, longitudinal surveys were completed by twins, their parents, siblings, spouses, and offspring. Data are available for 33,957 participants, with most individuals having completed two or more surveys. Smaller projects provided in-depth phenotyping, including measurements of the autonomic nervous system, neurocognitive function, and brain imaging. For 46% of the ANTR participants, DNA samples are available and whole genome scans have been obtained in more than 11,000 individuals. These data have resulted in numerous studies on heritability, gene x environment interactions, and causality, as well as gene finding studies. In the future, these studies will continue with collection of additional phenotypes, such as metabolomic and telomere length data, and detailed genetic information provided by DNA and RNA sequencing. Record linkage to national registers will allow the study of morbidity and mortality, thus providing insight into the development of health, lifestyle, and behavior across the lifespan.

Genetic and environmental stability in attention problems across the lifespan: evidence from the Netherlands twin register

OBJECTIVE

To review findings on attention-deficit/hyperactivity disorder and attention problems (AP) in children, adolescents, and adults, as established in the database of the Netherlands Twin Register and increase the understanding of stability in AP across the lifespan as a function of genetic and environmental influences.

METHOD

A longitudinal model was fitted on Netherlands Twin Register AP scores from 44,607 child (<12-year-old), adolescent (12- to 18-year-old), and adult (>18-year-old) twins.

RESULTS

Mean AP showed a downward trend with age. Age-to-age correlations ranged from 0.33 (50-≥60 years old) to 0.73 (10-12 years old). Stability in individual differences in AP was due to genetic and environmental factors, and change was due primarily to environmental factors. Nonadditive genetic influences were present from childhood to adulthood. Total genetic variance decreased slightly throughout aging, whereas environmental variance increased substantially with the switch from maternal to self-ratings at 12 years of age. As a result, heritability coefficients decreased from 0.70 to 0.74 in childhood (maternal ratings) to 0.51 to 0.56 in adolescence (self-ratings), and 0.40 to 0.54 in adulthood (self-ratings). In childhood, male subjects scored higher than female subjects. After the rater switch at 12 years of age, female subjects tended to score higher than male subjects.

CONCLUSIONS

Stability of AP is the result of genetic and environmental stability. The decrease in estimated heritability at 12 years of age is due to an increase in occasion-specific environmental variance and likely reflects a methodologic effect. Because environmental influences have lasting effects on AP, their early detection is crucial.

The Young Netherlands Twin Register (YNTR): longitudinal twin and family studies in over 70,000 children

The Netherlands Twin Register (NTR) began in 1987 with data collection in twins and their families, including families with newborn twins and triplets. Twenty-five years later, the NTR has collected at least one survey for 70,784 children, born after 1985. For the majority of twins, longitudinal data collection has been done by age-specific surveys. Shortly after giving birth, mothers receive a first survey with items on pregnancy and birth. At age 2, a survey on growth and achievement of milestones is sent. At ages 3, 7, 9/10, and 12 parents and teachers receive a series of surveys that are targeted at the development of emotional and behavior problems. From age 14 years onward, adolescent twins and their siblings report on their behavior problems, health, and lifestyle. When the twins are 18 years and older, parents are also invited to take part in survey studies. In sub-groups of different ages, in-depth phenotyping was done for IQ, electroencephalography , MRI, growth, hormones, neuropsychological assessments, and cardiovascular measures. DNA and biological samples have also been collected and large numbers of twin pairs and parents have been genotyped for zygosity by either micro-satellites or sets of short nucleotide polymorphisms and repeat polymorphisms in candidate genes. Subject recruitment and data collection is still ongoing and the longitudinal database is growing. Data collection by record linkage in the Netherlands is beginning and we expect these combined longitudinal data to provide increased insights into the genetic etiology of development of mental and physical health in children and adolescents.

Twins, tissue, and time: an assessment of SNPs and CNVs

With the desire to assess genetic variation across the lifespan in large-scale collaborative projects, one question is whether inference of copy number (CN) is sensitive to the source of material for deoxyribonucleic acid (DNA) analysis (e.g., blood and buccal) and another question is whether CN is stable as individual sage. Here, we address these questions by applying Affymetrix 6.0 single nucleotide polymorphism (SNP)micro-arrays to 1,472 DNA samples from 710 individuals from the Netherlands Twin Register, including twin and non-twin individuals (372 with buccal and blood derived DNA and 388 with longitudinal data).Similar concordance for CN and genotype inference between samples from the same individual [or from the monozygotic (MZ) co-twins] was found for blood and buccal tissues. There was a small but statistically significant decrease in across-tissue concordance compared with concordance of samples from the same tissue type. No temporal effect was seen on CN variation from the 388 individuals sampled at two time points ranging from 1 to 12 years apart. The majority of our individuals were sampled at age younger than 20 years. Genotype concordance was very high (~ > 99%) between co-twins from 43 MZ pairs. For75 dizygotic (DZ) pairs, ~was ~65%. CN estimates were highly consistent between co-twins from MZ pairs for both deletions (f?2 ~ 90%) and duplications (~ ~ 86%). For DZ, these were similar for within-individual comparisons, but naturally lower between co-twins (~ ~ 50-60%). These results suggest that DNA from buccal samples perform as well as DNA from blood samples on the current generation of micro-array technologies.

The continuing value of twin studies in the omics era

The classical twin study has been a powerful heuristic in biomedical, psychiatric and behavioural research for decades. Twin registries worldwide have collected biological material and longitudinal phenotypic data on tens of thousands of twins, providing a valuable resource for studying complex phenotypes and their underlying biology. In this Review, we consider the continuing value of twin studies in the current era of molecular genetic studies. We conclude that classical twin methods combined with novel technologies represent a powerful approach towards identifying and understanding the molecular pathways that underlie complex traits.