Trisomy 21 Alters DNA Methylation in Parent-of-Origin-Dependent and -Independent Manners

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

Genomic map of candidate human imprint control regions: the imprintome

Imprinted genes – critical for growth, metabolism, and neuronal function – are expressed from one parental allele. Parent-of-origin-dependent CpG methylation regulates this expression at imprint control regions (ICRs). Since ICRs are established before tissue specification, these methylation marks are similar across cell types. Thus, they are attractive for investigating the developmental origins of adult diseases using accessible tissues, but remain unknown. We determined genome-wide candidate ICRs in humans by performing whole-genome bisulphite sequencing (WGBS) of DNA derived from the three germ layers and from gametes. We identified 1,488 hemi-methylated candidate ICRs, including 19 of 25 previously characterized ICRs (https://humanicr.org/). Gamete methylation approached 0% or 100% in 332 ICRs (178 paternally and 154 maternally methylated), supporting parent-of-origin-specific methylation, and 65% were in well-described CTCF-binding or DNaseI hypersensitive regions. This draft of the human imprintome will allow for the systematic determination of the role of early-acquired imprinting dysregulation in the pathogenesis of human diseases and developmental and behavioural disorders.

Aberrant hypomethylation at imprinted differentially methylated regions is involved in biparental placental mesenchymal dysplasia

BACKGROUND

Placental mesenchymal dysplasia (PMD) is a morphological abnormality resembling partial hydatidiform moles. It is often associated with androgenetic/biparental mosaicism (ABM) and complicated by Beckwith-Wiedemann syndrome (BWS), an imprinting disorder. These phenomena suggest an association between PMD and aberrant genomic imprinting, particularly of CDKN1C and IGF2. The existence of another type of PMD containing the biparental genome has been reported. However, the frequency and etiology of biparental PMD are not yet fully understood.

RESULTS

We examined 44 placental specimens from 26 patients with PMD: 19 of these were macroscopically normal and 25 exhibited macroscopic PMD. Genotyping by DNA microarray or short tandem repeat analysis revealed that approximately 35% of the macroscopic PMD specimens could be classified as biparental, while the remainder were ABM. We performed a DNA methylation analysis using bisulfite pyrosequencing of 15 placenta-specific imprinted differentially methylated regions (DMRs) and 36 ubiquitous imprinted DMRs. As expected, most DMRs in the macroscopic PMD specimens with ABM exhibited the paternal epigenotype. Importantly, the biparental macroscopic PMD specimens exhibited frequent aberrant hypomethylation at seven of the placenta-specific DMRs. Allelic expression analysis using single-nucleotide polymorphisms revealed that five imprinted genes associated with these aberrantly hypomethylated DMRs were biallelically expressed. Frequent aberrant hypomethylation was observed at five ubiquitous DMRs, including GRB10 but not ICR2 or ICR1, which regulate the expression of CDKN1C and IGF2, respectively. Whole-exome sequencing performed on four biparental macroscopic PMD specimens did not reveal any pathological genetic abnormalities. Clinical and molecular analyses of babies born from pregnancies with PMD revealed four cases with BWS, each exhibiting different molecular characteristics, and those between BWS and PMD specimens were not always the same.

CONCLUSION

These data clarify the prevalence of biparental PMD and ABM-PMD and strongly implicate hypomethylation of DMRs in the pathogenesis of biparental PMD, particularly placenta-specific DMRs and the ubiquitous GRB10, but not ICR2 or ICR1. Aberrant hypomethylation of DMRs was partial, indicating that it occurs after fertilization. PMD is an imprinting disorder, and it may be a missing link between imprinting disorders and placental disorders incompatible with life, such as complete hydatidiform moles and partial hydatidiform moles.

Pervasive Inter-Individual Variation in Allele-Specific Expression in Monozygotic Twins

Despite being developed from one zygote, heterokaryotypic monozygotic (MZ) co-twins exhibit discordant karyotypes. Epigenomic studies in biological samples from heterokaryotypic MZ co-twins are of the most significant value for assessing the effects on gene- and allele-specific expression of an extranumerary chromosomal copy or structural chromosomal disparities in otherwise nearly identical germline genetic contributions. Here, we use RNA-Seq data from existing repositories to establish within-pair correlations for the breadth and magnitude of allele-specific expression (ASE) in heterokaryotypic MZ co-twins discordant for trisomy 21 and maternal 21q inheritance, as well as homokaryotypic co-twins. We show that there is a genome-wide disparity at ASE sites between the heterokaryotypic MZ co-twins. Although most of the disparity corresponds to changes in the magnitude of biallelic imbalance, ASE sites switching from either strictly monoallelic to biallelic imbalance or the reverse occur in few genes that are known or predicted to be imprinted, subject to X-chromosome inactivation or A-to-I(G) RNA edited. We also uncovered comparable ASE differences between homokaryotypic MZ twins. The extent of ASE discordance in MZ twins (2.7%) was about 10-fold lower than the expected between pairs of unrelated, non-twin males or females. The results indicate that the observed within-pair dissimilarities in breadth and magnitude of ASE sites in the heterokaryotypic MZ co-twins could not solely be attributable to the aneuploidy and the missing allelic heritability at 21q.

Genetic and epigenetic pathways in Down syndrome: Insights to the brain and immune system from humans and mouse models

The presence of an extra copy of human chromosome 21 (Hsa21) leads to a constellation of phenotypic manifestations in Down syndrome (DS), including prominent effects on the brain and immune system. Intensive efforts to unravel the molecular mechanisms underlying these phenotypes may help developing effective therapies, both in DS and in the general population. Here we review recent progress in genetic and epigenetic analysis of trisomy 21 (Ts21). New mouse models of DS based on syntenic conservation of segments of the mouse and human chromosomes are starting to clarify the contributions of chromosomal subregions and orthologous genes to specific phenotypes in DS. The expression of genes on Hsa21 is regulated by epigenetic mechanisms, and with recent findings of highly recurrent gene-specific changes in DNA methylation patterns in brain and immune system cells with Ts21, the epigenomics of DS has become an active research area. Here we highlight the value of combining human studies with mouse models for defining DS critical genes and understanding the trans-acting effects of a simple chromosomal aneuploidy on genome-wide epigenetic patterning. These genetic and epigenetic studies are starting to uncover fundamental biological mechanisms, leading to insights that may soon become therapeutically relevant.

Genome-wide allele-specific methylation is enriched at gene regulatory regions in a multi-generation pedigree from the Norfolk Island isolate

Background

Allele-specific methylation (ASM) occurs when DNA methylation patterns exhibit asymmetry among alleles. ASM occurs at imprinted loci, but its presence elsewhere across the human genome is indicative of wider importance in terms of gene regulation and disease risk. Here, we studied ASM by focusing on blood-based DNA collected from 24 subjects comprising a 3-generation pedigree from the Norfolk Island genetic isolate. We applied a genome-wide bisulphite sequencing approach with a genotype-independent ASM calling method to map ASM across the genome. Regions of ASM were then tested for enrichment at gene regulatory regions using Genomic Association Test (GAT) tool.

Results

In total, we identified 1.12 M CpGs of which 147,170 (13%) exhibited ASM (P ≤ 0.05). When including contiguous ASM signal spanning ≥ 2 CpGs, this condensed to 12,761 ASM regions (AMRs). These AMRs tagged 79% of known imprinting regions and most (98.1%) co-localised with known single nucleotide variants. Notably, miRNA and lncRNA showed a 3.3- and 1.8-fold enrichment of AMRs, respectively (P < 0.005). Also, the 5′ UTR and start codons each showed a 3.5-fold enrichment of AMRs (P < 0.005). There was also enrichment of AMRs observed at subtelomeric regions of many chromosomes. Five out of 11 large AMRs localised to the protocadherin cluster on chromosome 5.

Conclusions

This study shows ASM extends far beyond genomic imprinting in humans and that gene regulatory regions are hotspots for ASM. Future studies of ASM in pedigrees should help to clarify transgenerational inheritance patterns in relation to genotype and disease phenotypes.

Genome-Wide Association Study Identifies a Susceptibility Locus for Comitant Esotropia and Suggests a Parent-of-Origin Effect

Purpose

To identify genetic variants conferring susceptibility to esotropia. Esotropia is the most common form of comitant strabismus, has its highest incidence in European ancestry populations, and is believed to be inherited as a complex trait.

Methods

White European American discovery cohorts with nonaccommodative (826 cases and 2991 controls) or accommodative (224 cases and 749 controls) esotropia were investigated. White European Australian and United Kingdom cohorts with nonaccommodative (689 cases and 1448 controls) or accommodative (66 cases and 264 controls) esotropia were tested for replication. We performed a genome-wide case-control association study using a mixed linear additive model. Meta-analyses of discovery and replication cohorts were then conducted.

Results

A significant association with nonaccommodative esotropia was discovered (odds ratio [OR] = 1.41, P = 2.84 × 10-09) and replicated (OR = 1.23, P = 0.01) at rs2244352 [T] located within intron 1 of the WRB (tryptophan rich basic protein) gene on chromosome 21 (meta-analysis OR = 1.33, P = 9.58 × 10-11). This single nucleotide polymorphism (SNP) is differentially methylated, and there is a statistically significant skew toward paternal inheritance in the discovery cohort. Meta-analysis of the accommodative discovery and replication cohorts identified an association with rs912759 [T] (OR = 0.59, P = 1.89 × 10-08), an intergenic SNP on chromosome 1p31.1.

Conclusions

This is the first genome-wide association study (GWAS) to identify significant associations in esotropia and suggests a parent-of-origin effect. Additional cohorts will permit replication and extension of these findings. Future studies of rs2244352 and WRB should provide insight into pathophysiological mechanisms underlying comitant strabismus.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Advances in understanding the association between Down syndrome and Hirschsprung disease (DS-HSCR)

The clinical association between Trisomy 21 (Down syndrome) and aganglionosis (Hirschsprung disease; DS-HSCR) is well-established, being of the order of 5% and remains the most common congenital association with Hirschsprung disease. However, little consensus exists as to the possible etiologic and genetic factors influencing this association. Recent research has identified a number of levels at which development of the enteric nervous system is potentially affected in Trisomy 21. These include a decreased central pool of available neuroblasts for migration into the enteric nervous system, abnormal neuroblast type, poor synaptic nerve function and early germline gene-related influences on the migrating neuroblasts due to genetic mutations of a number of important developmental genes, and possible somatic mutations resulting from alterations in the local tissue microenvironment. In this paper, we review available evidence for this association. In addition, we provide evidence of both germline and somatic gene mutations suggesting causation. Although the picture is complex, recent associations between specific RET proto-oncogene variations have been shown to be significant in Down syndrome patients with Hirschsprung disease, as they probably interfere with vital RET functions in the development of the autonomic and enteric nervous systems, increasing the risk of disturbed normal function. In addition, we explore potential role of other facilitatory influence of other susceptibility genes as well as potential other chromosome 21 gene actions and the microenvironment on the Down syndrome gastro-intestinal tract. The various ways in which trisomy of chromosome influences the enteric nervous system are becoming clearer. The sum of these effects influences the outcome of surgery in Down syndrome patients with Hirschsprung Disease.

Widespread domain-like perturbations of DNA methylation in whole blood of Down syndrome neonates

INTRODUCTION

Down syndrome (DS) is the most frequent genetic cause of intellectual disability. Despite the fact that more than 50 years have passed since the discovery of its genetic aberrations, the exact pathogenesis of the DS phenotype has remained largely unexplained. It was recently hypothesized that the DS pathogenesis involves complex (epi)genetic, molecular and cellular determinants. To date, many reports have addressed epigenetic aberrations associated with DS at different developmental stages/ages and tissue types, but to our best knowledge not in DS newborns. This study aimed to investigate genome-wide methylation patterns in DS newborns compared to non-trisomic newborns.

METHOD

We analyzed blood samples obtained from ten newborns with DS and five age-matched non-trisomic newborns. Epigenetic profiles were obtained from extracted DNA using the Illumina Infinium 450K array. Since aberrant blood cell distribution is known to be present in DS, we applied two distinct models: with and without correction for estimated blood cell distribution.

RESULTS

Differentially methylated position (DMP) analysis of the uncorrected model detected 19525 significant hits (51,2% hypomethylated). In the corrected model, we found 121953 significant DMPs (49,8% hypomethylated). Independent of the used model we observed a chromosome 21 dosage effect. Moreover, we detected 46 and 145 differentially methylated regions in the uncorrected and corrected model respectively, both showing hypomethylation overrepresentation. Replication analyses of DMPs and DMRs found by Bacalini et al. (2015) showed a large overlap.

CONCLUSION

In this study, we found methylation profile differences between DS newborns and controls reflecting a systematically affected epigenetic profile. The observed chromosome 21 dosage effect suggests the involvement of affected essential regulatory factors/regions or altered expression of chromatin modeling enzymes located on chromosome 21. Additional research is necessary to substantiate these hypotheses.

Maternal 5mCpG Imprints at the PARD6G-AS1 and GCSAML Differentially Methylated Regions Are Decoupled From Parent-of-Origin Expression Effects in Multiple Human Tissues

A hallmark of imprinted genes in mammals is the occurrence of parent-of-origin-dependent asymmetry of DNA cytosine methylation (5^mC) of alleles at CpG islands (CGIs) in their promoter regions. This 5^mCpG asymmetry between the parental alleles creates allele-specific imprinted differentially methylated regions (iDMRs). iDMRs are often coupled to the transcriptional repression of the methylated allele and the activation of the unmethylated allele in a tissue-specific, developmental-stage-specific and/or isoform-specific fashion. iDMRs function as regulatory platforms, built through the recruitment of chemical modifications to histones to achieve differential, parent-of-origin-dependent chromatin segmentation states. Here, we used a comparative computational data mining approach to identify 125 novel constitutive candidate iDMRs that integrate the maximal number of allele-specific methylation region records overlapping CGIs in human methylomes. Twenty-nine candidate iDMRs display gametic 5^mCpG asymmetry, and another 96 are candidate secondary iDMRs. We established the maternal origin of the 5^mCpG imprints of one gametic (PARD6G-AS1) and one secondary (GCSAML) iDMRs. We also found a constitutively hemimethylated, nonimprinted domain at the PWWP2AP1 promoter CGI with oocyte-derived methylation asymmetry. Given that the 5^mCpG level at the iDMRs is not a sufficient criterion to predict active or silent locus states and that iDMRs can regulate genes from a distance of more than 1 Mb, we used RNA-Seq experiments from the Genotype-Tissue Expression project and public archives to assess the transcriptional expression profiles of SNPs across 4.6 Mb spans around the novel maternal iDMRs. We showed that PARD6G-AS1 and GCSAML are expressed biallelically in multiple tissues. We found evidence of tissue-specific monoallelic expression of ZNF124 and OR2L13, located 363 kb upstream and 419 kb downstream, respectively, of the GCSAML iDMR. We hypothesize that the GCSAML iDMR regulates the tissue-specific, monoallelic expression of ZNF124 but not of OR2L13. We annotated the non-coding epigenomic marks in the two maternal iDMRs using data from the Roadmap Epigenomics project and showed that the PARD6G-AS1 and GCSAML iDMRs achieve contrasting activation and repression chromatin segmentations. Lastly, we found that the maternal 5^mCpG imprints are perturbed in several hematopoietic cancers. We conclude that the maternal 5^mCpG imprints at PARD6G-AS1 and GCSAML iDMRs are decoupled from parent-of-origin transcriptional expression effects in multiple tissues.

A Pair of Maternal Chromosomes Derived from Meiotic Nondisjunction in Trisomy 21 Affects Nuclear Architecture and Transcriptional Regulation

Eukaryotic genomes are organised into complex higher-order structures within the nucleus, and the three-dimensional arrangement of chromosomes is functionally important for global gene regulation. The existence of supernumerary chromosome 21 in Down syndrome may perturb the nuclear architecture at different levels, which is normally optimised to maintain the physiological balance of gene expression. However, it has not been clearly elucidated whether and how aberrant configuration of chromosomes affects gene activities. To investigate the effects of trisomy 21 on nuclear organisation and gene expression, we performed three-dimensional fluorescent imaging analysis of chromosome-edited human induced pluripotent stem cells (iPSCs), which enabled identification of the parental origin of the three copies of chromosome 21. We found that two copies of maternal chromosomes resulting from meiotic nondisjunction had a higher tendency to form an adjacent pair and were located relatively distant from the nuclear membrane, suggesting the conserved interaction between these homologous chromosomes. Transcriptional profiling of parental-origin-specific corrected disomy 21 iPSC lines indicated upregulated expression of the maternal alleles for a group of genes, which was accompanied by a fluctuating expression pattern. These results suggest the unique effects of a pair of maternal chromosomes in trisomy 21, which may contribute to the pathological phenotype.

Trans-acting epigenetic effects of chromosomal aneuploidies: lessons from Down syndrome and mouse models

An important line of postgenomic research seeks to understand how genetic factors can influence epigenetic patterning. Here we review epigenetic effects of chromosomal aneuploidies, focusing on findings in Down syndrome (DS, trisomy 21). Recent work in human DS and mouse models has shown that the extra chromosome 21 acts in trans to produce epigenetic changes, including differential CpG methylation (DS-DM), in specific sets of downstream target genes, mostly on other chromosomes. Mechanistic hypotheses emerging from these data include roles of chromosome 21-linked methylation pathway genes (DNMT3L and others) and transcription factor genes (RUNX1, OLIG2, GABPA, ERG and ETS2) in shaping the patterns of DS-DM. The findings may have broader implications for trans-acting epigenetic effects of chromosomal and subchromosomal aneuploidies in other human developmental and neuropsychiatric disorders, and in cancers.