Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice.

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically shows phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption-the root of the pathogenesis-is similar in the different disease-relevant cell types. This is possible in principle, because all these cell types are subject to effects from the same causative gene, which has the same kind of function (e.g., methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2 and find that the chromatin accessibility changes in neurons are mostly distinct from changes in B or T cells. This is not because the neuronal accessibility changes occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that regulatory elements disrupted in B/T cells do show chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators and suggest that blood-derived episignatures, although useful diagnostically, may not be well suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Growth deficiency in a mouse model of Kabuki syndrome 2 bears mechanistic similarities to Kabuki syndrome 1.

Growth deficiency is a characteristic feature of both Kabuki syndrome 1 (KS1) and Kabuki syndrome 2 (KS2), Mendelian disorders of the epigenetic machinery with similar phenotypes but distinct genetic etiologies. We previously described skeletal growth deficiency in a mouse model of KS1 and further established that a Kmt2d-/- chondrocyte model of KS1 exhibits precocious differentiation. Here we characterized growth deficiency in a mouse model of KS2, Kdm6atm1d/+. We show that Kdm6atm1d/+ mice have decreased femur and tibia length compared to controls and exhibit abnormalities in cortical and trabecular bone structure. Kdm6atm1d/+ growth plates are also shorter, due to decreases in hypertrophic chondrocyte size and hypertrophic zone height. Given these disturbances in the growth plate, we generated Kdm6a-/- chondrogenic cell lines. Similar to our prior in vitro model of KS1, we found that Kdm6a-/- cells undergo premature, enhanced differentiation towards chondrocytes compared to Kdm6a+/+ controls. RNA-seq showed that Kdm6a-/- cells have a distinct transcriptomic profile that indicates dysregulation of cartilage development. Finally, we performed RNA-seq simultaneously on Kmt2d-/-, Kdm6a-/-, and control lines at Days 7 and 14 of differentiation. This revealed surprising resemblance in gene expression between Kmt2d-/- and Kdm6a-/- at both time points and indicates that the similarity in phenotype between KS1 and KS2 also exists at the transcriptional level.

Missense variants causing Wiedemann-Steiner syndrome preferentially occur in the KMT2A-CXXC domain and are accurately classified using AlphaFold2.

Wiedemann-Steiner syndrome (WDSTS) is a neurodevelopmental disorder caused by de novo variants in KMT2A, which encodes a multi-domain histone methyltransferase. To gain insight into the currently unknown pathogenesis of WDSTS, we examined the spatial distribution of likely WDSTS-causing variants across the 15 different domains of KMT2A. Compared to variants in healthy controls, WDSTS variants exhibit a 61.9-fold overrepresentation within the CXXC domain-which mediates binding to unmethylated CpGs-suggesting a major role for this domain in mediating the phenotype. In contrast, we find no significant overrepresentation within the catalytic SET domain. Corroborating these results, we find that hippocampal neurons from Kmt2a-deficient mice demonstrate disrupted histone methylation (H3K4me1 and H3K4me3) preferentially at CpG-rich regions, but this has no systematic impact on gene expression. Motivated by these results, we combine accurate prediction of the CXXC domain structure by AlphaFold2 with prior biological knowledge to develop a classification scheme for missense variants in the CXXC domain. Our classifier achieved 92.6% positive and 92.9% negative predictive value on a hold-out test set. This classification performance enabled us to subsequently perform an in silico saturation mutagenesis and classify a total of 445 variants according to their functional effects. Our results yield a novel insight into the mechanistic basis of WDSTS and provide an example of how AlphaFold2 can contribute to the in silico characterization of variant effects with very high accuracy, suggesting a paradigm potentially applicable to many other Mendelian disorders.

Promoter CpG Density Predicts Downstream Gene Loss-of-Function Intolerance.

The aggregation and joint analysis of large numbers of exome sequences has recently made it possible to derive estimates of intolerance to loss-of-function (LoF) variation for human genes. Here, we demonstrate strong and widespread coupling between genic LoF intolerance and promoter CpG density across the human genome. Genes downstream of the most CpG-rich promoters (top 10% CpG density) have a 67.2% probability of being highly LoF intolerant, using the LOEUF metric from gnomAD. This is in contrast to 7.4% of genes downstream of the most CpG-poor (bottom 10% CpG density) promoters. Combining promoter CpG density with exonic and promoter conservation explains 33.4% of the variation in LOEUF, and the contribution of CpG density exceeds the individual contributions of exonic and promoter conservation. We leverage this to train a simple and easily interpretable predictive model that outperforms other existing predictors and allows us to classify 1,760 genes-which are currently unascertained in gnomAD-as highly LoF intolerant or not. These predictions have the potential to aid in the interpretation of novel variants in the clinical setting. Moreover, our results reveal that high CpG density is not merely a generic feature of human promoters but is preferentially encountered at the promoters of the most selectively constrained genes, calling into question the prevailing view that CpG islands are not subject to selection.

Coexpression patterns define epigenetic regulators associated with neurological dysfunction.

Coding variants in epigenetic regulators are emerging as causes of neurological dysfunction and cancer. However, a comprehensive effort to identify disease candidates within the human epigenetic machinery (EM) has not been performed; it is unclear whether features exist that distinguish between variation-intolerant and variation-tolerant EM genes, and between EM genes associated with neurological dysfunction versus cancer. Here, we rigorously define 295 genes with a direct role in epigenetic regulation (writers, erasers, remodelers, readers). Systematic exploration of these genes reveals that although individual enzymatic functions are always mutually exclusive, readers often also exhibit enzymatic activity (dual-function EM genes). We find that the majority of EM genes are very intolerant to loss-of-function variation, even when compared to the dosage sensitive transcription factors, and we identify 102 novel EM disease candidates. We show that this variation intolerance is driven by the protein domains encoding the epigenetic function, suggesting that disease is caused by a perturbed chromatin state. We then describe a large subset of EM genes that are coexpressed within multiple tissues. This subset is almost exclusively populated by extremely variation-intolerant genes and shows enrichment for dual-function EM genes. It is also highly enriched for genes associated with neurological dysfunction, even when accounting for dosage sensitivity, but not for cancer-associated EM genes. Finally, we show that regulatory regions near epigenetic regulators are genetically important for common neurological traits. These findings prioritize novel disease candidate EM genes and suggest that this coexpression plays a functional role in normal neurological homeostasis.

A screen of 1,049 schizophrenia and 30 Alzheimer's-associated variants for regulatory potential.

Recent genome-wide association studies (GWAS) identified numerous schizophrenia (SZ) and Alzheimer's disease (AD) associated loci, most outside protein-coding regions and hypothesized to affect gene transcription. We used a massively parallel reporter assay to screen, 1,049 SZ and 30 AD variants in 64 and nine loci, respectively for allele differences in driving reporter gene expression. A library of synthetic oligonucleotides assaying each allele five times was transfected into K562 chronic myelogenous leukemia lymphoblasts and SK-SY5Y human neuroblastoma cells. One hundred forty eight variants showed allelic differences in K562 and 53 in SK-SY5Y cells, on average 2.6 variants per locus. Nine showed significant differences in both lines, a modest overlap reflecting different regulatory landscapes of these lines that also differ significantly in chromatin marks. Eight of nine were in the same direction. We observe no preference for risk alleles to increase or decrease expression. We find a positive correlation between the number of SNPs in linkage disequilibrium and the proportion of functional SNPs supporting combinatorial effects that may lead to haplotype selection. Our results prioritize future functional follow up of disease associated SNPs to determine the driver GWAS variant(s), at each locus and enhance our understanding of gene regulation dynamics.

Missense variants in the chromatin remodeler CHD1 are associated with neurodevelopmental disability.

BACKGROUND: The list of Mendelian disorders of the epigenetic machinery has expanded rapidly during the last 5 years. A few missense variants in the chromatin remodeler CHD1 have been found in several large-scale sequencing efforts focused on uncovering the genetic aetiology of autism. OBJECTIVES: To explore whether variants in CHD1 are associated with a human phenotype. METHODS: We used GeneMatcher to identify other physicians caring for patients with variants in CHD1. We also explored the epigenetic consequences of one of these variants in cultured fibroblasts. RESULTS: Here we describe six CHD1 heterozygous missense variants in a cohort of patients with autism, speech apraxia, developmental delay and facial dysmorphic features. Importantly, three of these variants occurred de novo. We also report on a subject with a de novo deletion covering a large fraction of the CHD1 gene without any obvious neurological phenotype. Finally, we demonstrate increased levels of the closed chromatin modification H3K27me3 in fibroblasts from a subject carrying a de novo variant in CHD1. CONCLUSIONS: Our results suggest that variants in CHD1 can lead to diverse phenotypic outcomes; however, the neurodevelopmental phenotype appears to be limited to patients with missense variants, which is compatible with a dominant negative mechanism of disease.

A mouse model of Weaver syndrome displays overgrowth and excess osteogenesis reversible with KDM6A/6B inhibition.

Weaver syndrome is a Mendelian disorder of the epigenetic machinery (MDEM) caused by germline pathogenic variants in EZH2, which encodes the predominant H3K27 methyltransferase and key enzymatic component of Polycomb repressive complex 2 (PRC2). Weaver syndrome is characterized by striking overgrowth and advanced bone age, intellectual disability, and distinctive facies. We generated a mouse model for the most common Weaver syndrome missense variant, EZH2 p.R684C. Ezh2R684C/R684C mouse embryonic fibroblasts (MEFs) showed global depletion of H3K27me3. Ezh2R684C/+ mice had abnormal bone parameters, indicative of skeletal overgrowth, and Ezh2R684C/+ osteoblasts showed increased osteogenic activity. RNA-Seq comparing osteoblasts differentiated from Ezh2R684C/+, and Ezh2+/+ BM-mesenchymal stem cells (BM-MSCs) indicated collective dysregulation of the BMP pathway and osteoblast differentiation. Inhibition of the opposing H3K27 demethylases KDM6A and KDM6B substantially reversed the excessive osteogenesis in Ezh2R684C/+ cells both at the transcriptional and phenotypic levels. This supports both the ideas that writers and erasers of histone marks exist in a fine balance to maintain epigenome state and that epigenetic modulating agents have therapeutic potential for the treatment of MDEMs.

KMT2D regulates activation, localization, and integrin expression by T-cells.

Individuals with Kabuki syndrome present with immunodeficiency; however, how pathogenic variants in the gene encoding the histone-modifying enzyme lysine methyltransferase 2D (KMT2D) lead to immune alterations remain poorly understood. Following up on our prior report of KMT2D-altered integrin expression in B-cells, we performed targeted analyses of KMT2D's influence on integrin expression in T-cells throughout development (thymocytes through peripheral T-cells) in murine cells with constitutive- and conditional-targeted Kmt2d deletion. Using high-throughput RNA-sequencing and flow cytometry, we reveal decreased expression (both at the transcriptional and translational levels) of a cluster of leukocyte-specific integrins, which perturb aspects of T-cell activation, maturation, adhesion/localization, and effector function. H3K4me3 ChIP-PCR suggests that these evolutionary similar integrins are under direct control of KMT2D. KMT2D loss also alters multiple downstream programming/signaling pathways, including integrin-based localization, which can influence T-cell populations. We further demonstrated that KMT2D deficiency is associated with the accumulation of murine CD8+ single-positive (SP) thymocytes and shifts in both human and murine peripheral T-cell populations, including the reduction of the CD4+ recent thymic emigrant (RTE) population. Together, these data show that the targeted loss of Kmt2d in the T-cell lineage recapitulates several distinct features of Kabuki syndrome-associated immune deficiency and implicates epigenetic mechanisms in the regulation of integrin signaling.

Pumping the brakes on RNA velocity by understanding and interpreting RNA velocity estimates.

BACKGROUND: RNA velocity analysis of single cells offers the potential to predict temporal dynamics from gene expression. In many systems, RNA velocity has been observed to produce a vector field that qualitatively reflects known features of the system. However, the limitations of RNA velocity estimates are still not well understood. RESULTS: We analyze the impact of different steps in the RNA velocity workflow on direction and speed. We consider both high-dimensional velocity estimates and low-dimensional velocity vector fields mapped onto an embedding. We conclude the transition probability method for mapping velocity estimates onto an embedding is effectively interpolating in the embedding space. Our findings reveal a significant dependence of the RNA velocity workflow on smoothing via the k-nearest-neighbors (k-NN) graph of the observed data. This reliance results in considerable estimation errors for both direction and speed in both high- and low-dimensional settings when the k-NN graph fails to accurately represent the true data structure; this is an unknown feature of real data. RNA velocity performs poorly at estimating speed in both low- and high-dimensional spaces, except in very low noise settings. We introduce a novel quality measure that can identify when RNA velocity should not be used. CONCLUSIONS: Our findings emphasize the importance of choices in the RNA velocity workflow and highlight critical limitations of data analysis. We advise against over-interpreting expression dynamics using RNA velocity, particularly in terms of speed. Finally, we emphasize that the use of RNA velocity in assessing the correctness of a low-dimensional embedding is circular.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice.

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically exhibits phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption - the root of the pathogenesis - is similar in the different disease-relevant cell types. This is possible in principle, since all these cell-types are subject to effects from the same causative gene, that has the same kind of function (e.g. methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2, and find that the chromatin accessibility abnormalities in neurons are mostly distinct from those in B or T cells. This is not because the neuronal abnormalities occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that the regions disrupted in B/T cells do exhibit chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators, and suggest that blood-derived "episignatures" may not be well-suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Novel mouse model of Weaver syndrome displays overgrowth and excess osteogenesis reversible with KDM6A/6B inhibition.

Weaver syndrome is a Mendelian disorder of the epigenetic machinery (MDEM) caused by germline pathogenic variants in EZH2, which encodes the predominant H3K27 methyltransferase and key enzymatic component of Polycomb repressive complex 2 (PRC2). Weaver syndrome is characterized by striking overgrowth and advanced bone age, intellectual disability, and distinctive facies. We generated a mouse model for the most common Weaver syndrome missense variant, EZH2 p.R684C. Ezh2R684C/R684C mouse embryonic fibroblasts (MEFs) showed global depletion of H3K27me3. Ezh2R684C/+ mice had abnormal bone parameters indicative of skeletal overgrowth, and Ezh2R684C/+ osteoblasts showed increased osteogenic activity. RNA-seq comparing osteoblasts differentiated from Ezh2R684C/+ and Ezh2+/+ bone marrow mesenchymal stem cells (BM-MSCs) indicated collective dysregulation of the BMP pathway and osteoblast differentiation. Inhibition of the opposing H3K27 demethylases Kdm6a/6b substantially reversed the excessive osteogenesis in Ezh2R684C/+ cells both at the transcriptional and phenotypic levels. This supports both the ideas that writers and erasers of histone marks exist in a fine balance to maintain epigenome state, and that epigenetic modulating agents have therapeutic potential for the treatment of MDEMs.

Leveraging the Mendelian disorders of the epigenetic machinery to systematically map functional epigenetic variation.

Although each Mendelian Disorder of the Epigenetic Machinery (MDEM) has a different causative gene, there are shared disease manifestations. We hypothesize that this phenotypic convergence is a consequence of shared epigenetic alterations. To identify such shared alterations, we interrogate chromatin (ATAC-seq) and expression (RNA-seq) states in B cells from three MDEM mouse models (Kabuki [KS] type 1 and 2 and Rubinstein-Taybi type 1 [RT1] syndromes). We develop a new approach for the overlap analysis and find extensive overlap primarily localized in gene promoters. We show that disruption of chromatin accessibility at promoters often disrupts downstream gene expression, and identify 587 loci and 264 genes with shared disruption across all three MDEMs. Subtle expression alterations of multiple, IgA-relevant genes, collectively contribute to IgA deficiency in KS1 and RT1, but not in KS2. We propose that the joint study of MDEMs offers a principled approach for systematically mapping functional epigenetic variation in mammals.

Deficiency of TET3 leads to a genome-wide DNA hypermethylation episignature in human whole blood.

TET3 encodes an essential dioxygenase involved in epigenetic regulation through DNA demethylation. TET3 deficiency, or Beck-Fahrner syndrome (BEFAHRS; MIM: 618798), is a recently described neurodevelopmental disorder of the DNA demethylation machinery with a nonspecific phenotype resembling other chromatin-modifying disorders, but inconsistent variant types and inheritance patterns pose diagnostic challenges. Given TET3's direct role in regulating 5-methylcytosine and recent identification of syndrome-specific DNA methylation profiles, we analyzed genome-wide DNA methylation in whole blood of TET3-deficient individuals and identified an episignature that distinguishes affected and unaffected individuals and those with mono-allelic and bi-allelic pathogenic variants. Validation and testing of the episignature correctly categorized known TET3 variants and determined pathogenicity of variants of uncertain significance. Clinical utility was demonstrated when the episignature alone identified an affected individual from over 1000 undiagnosed cases and was confirmed upon distinguishing TET3-deficient individuals from those with 46 other disorders. The TET3-deficient signature - and the signature resulting from activating mutations in DNMT1 which normally opposes TET3 - are characterized by hypermethylation, which for BEFAHRS involves CpG sites that may be biologically relevant. This work expands the role of epi-phenotyping in molecular diagnosis and reveals genome-wide DNA methylation profiling as a quantitative, functional readout for characterization of this new biochemical category of disease.

Precocious neuronal differentiation and disrupted oxygen responses in Kabuki syndrome.

Chromatin modifiers act to coordinate gene expression changes critical to neuronal differentiation from neural stem/progenitor cells (NSPCs). Lysine-specific methyltransferase 2D (KMT2D) encodes a histone methyltransferase that promotes transcriptional activation and is frequently mutated in cancers and in the majority (>70%) of patients diagnosed with the congenital, multisystem intellectual disability disorder Kabuki syndrome 1 (KS1). Critical roles for KMT2D are established in various non-neural tissues, but the effects of KMT2D loss in brain cell development have not been described. We conducted parallel studies of proliferation, differentiation, transcription, and chromatin profiling in KMT2D-deficient human and mouse models to define KMT2D-regulated functions in neurodevelopmental contexts, including adult-born hippocampal NSPCs in vivo and in vitro. We report cell-autonomous defects in proliferation, cell cycle, and survival, accompanied by early NSPC maturation in several KMT2D-deficient model systems. Transcriptional suppression in KMT2D-deficient cells indicated strong perturbation of hypoxia-responsive metabolism pathways. Functional experiments confirmed abnormalities of cellular hypoxia responses in KMT2D-deficient neural cells and accelerated NSPC maturation in vivo. Together, our findings support a model in which loss of KMT2D function suppresses expression of oxygen-responsive gene programs important to neural progenitor maintenance, resulting in precocious neuronal differentiation in a mouse model of KS1.

Precocious chondrocyte differentiation disrupts skeletal growth in Kabuki syndrome mice.

Kabuki syndrome 1 (KS1) is a Mendelian disorder of the epigenetic machinery caused by mutations in the gene encoding KMT2D, which methylates lysine 4 on histone H3 (H3K4). KS1 is characterized by intellectual disability, postnatal growth retardation, and distinct craniofacial dysmorphisms. A mouse model (Kmt2d+/ßGeo) exhibits features of the human disorder and has provided insight into other phenotypes; however, the mechanistic basis of skeletal abnormalities and growth retardation remains elusive. Using high-resolution micro-CT, we show that Kmt2d+/ßGeo mice have shortened long bones and ventral bowing of skulls. In vivo expansion of growth plates within skulls and long bones suggests disrupted endochondral ossification as a common disease mechanism. Stable chondrocyte cell lines harboring inactivating mutations in Kmt2d exhibit precocious differentiation, further supporting this mechanism. A known inducer of chondrogenesis, SOX9, and its targets show markedly increased expression in Kmt2d-/- chondrocytes. By transcriptome profiling, we identify Shox2 as a putative KMT2D target. We propose that decreased KMT2D-mediated H3K4me3 at Shox2 releases Sox9 inhibition and thereby leads to enhanced chondrogenesis, providing a potentially novel and plausible explanation for precocious chondrocyte differentiation. Our findings provide insight into the pathogenesis of growth retardation in KS1 and suggest therapeutic approaches for this and related disorders.