Loading ALS Journal
Article Sections

Edited by

Zia-ur-RahmanUniversity of Lahore, Pakistan

Reviewed by

Mohammad Zahid MustafaDepartment of Inserm, Unité 1193, Université Paris-Saclay, Villejuif 94800, France
Abdul QadirDepartment of Plant Pathology, Agriculture College Quetta, Pakistan

Figures

Tables

Discovery and analysis of a novel mutation (G774E) in the Lysine Demethylase 6A (KDM6A) gene causing congenital heart disease with various neurodevelopmental disorders
Mahmood Rasool1, Majed Alsulami2,3, Ayat Mohammed Shorbaji1,4, Absarul Haque5, Loubna Siraj Mira1, Mohammad Basabrain1, Sherin Bakhashab1,4, Mohamed Nabil Alama7, Sajjad Karim1, Isam M. Abu Zeid3,6, Hisham N. Altayb3, Peter Natesan Pushparaj1
  1. Institute of Genomic Medicine Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
  2. Department of Biology, College of Science, University of Hafr Al Batin, Hafr Al Batin, Saudi Arabia
  3. Department of Biological Sciences, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
  4. Department of Biochemistry, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
  5. King Fahd Medical Research Center, Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
  6. Centre of Excellence in Bionanoscience Research, King Abdulaziz University, Jeddah, Saudi Arabia
  7. Department of Cardiology, King Abdulaziz University Hospital, King Abdulaziz University Jeddah, Saudi Arabia

Abstract

Background: Congenital heart defects (CHD) are the most common birth defects, affecting approximately 0.8% of live births worldwide. CHD impairs the function and structure of the heart and blood vessels. Due to this damage, blood flow is impaired, which affects many major organs, including the brain, and causes various neurodevelopmental disorders.

Methods: In this study, we recruited a five-generation pedigree for analysis. The proband was born with congenital heart disease and subsequently developed various neurodevelopmental disorders. To understand the causes of the disease, we performed clinical whole-exome sequencing and applied various bioinformatics tools to determine the pathogenicity of the mutation.

Results: Molecular investigation revealed a novel lethal mutation (c.2321G>A) in KDM6A, causing the substitution of Glycine to Glutamic acid (Gly774Glu). The mutation was further confirmed using Sanger sequencing. Various bioinformatics tools were used to predict the lethality of the mutations. KDM6A disruption causes various diseases, among which Kabuki syndrome is the most prevalent.

Conclusion: Our findings may aid in the further development of genome-based medicines, leading to a reduction in mortality rates and improved healthcare in newborns.

Keywords

Congenital heart disease, Kabuki syndrome, Neurodevelopmental disorder, Whole exome sequencing, Novel mutation

Introduction

Congenital heart disease (CHD) is defined as the birth of a child with structural abnormalities of the heart or great vessels and is prevalent in almost 0.8% of newborns [1]. CHD is responsible for most infant deaths worldwide [2]. In Saudi Arabia, the incidence remains largely unexplored, but a recent study estimated that the incidence is much higher, that is 1.48% [3]. The presence of genetic underpinnings in CHD is strongly considered, as the disease shows high recurrence risks and familial forms, as well as the presence of chromosomal anomalies in many cases of CHD [4]. To date, studies have confirmed the presence of more than 400 genes associated with CHD pathogenesis. Most genetic mutations are present in transcription factors, cell signaling transducers, and chromatin modifiers that interact with specific cells during differentiation, which are important for structural and functional heart development [5].

Most genes encode proteins that work in coherence and are connected by interactive functional networks, leading to a broad network of proteins that may be associated with CHDs [6, 7]. However, to date, the majority of CHD cases remain to be elucidated (~60%), owing to the complexity of the genetic etiology of CHD in human cases and the presence of genetic heterogeneity associated with CHD [8]. This leads to variable expressivity, and similar variants show different phenotypic variability or variable penetrance in some individuals. Therefore, not all CHD cases follow Mendelian inheritance patterns and are better represented by complex genetics.

The emergence of new techniques, such as whole-exome or whole-genome sequencing, has made it far easier to study complex genetics, such as CHD. Currently, the most common technique used is trio studies (using patient and both parents’ samples), where patients are sequenced with their unaffected parents to identify the pathogenic variants that are missing in the parents (de novo mutations). In familial forms of CHDs, multiple members of the same family are sequenced to identify the common variants transmitted to family members. Moreover, in large cohort studies, a large number of patient samples were collected from different regions along with healthy control samples, and sequencing was performed to identify the presence of a single gene or group of genes that were enriched in variants in patients compared to healthy individuals. These studies revealed a high number of genes expressed in the heart or at some stage of heart development [8].

After the successful detection of genes associated with CHD, prenatal diagnosis is the best way to prevent the disease in the next generation. Prenatal diagnosis supplements the family with proper guidance on prognosis, planning adequate obstetric follow-up, better preparation for intrauterine management for complex and uncommon cases, and well before time refers to planning for pregnant women to transfer to tertiary care facilities (capable of providing therapeutic catheterization and, if needed, heart surgery), which will significantly reduce the morbidity and mortality in the newborn [9].

In this study, we examined the genetic causes of congenital heart diseases and neurodevelopmental disorders in a large, inbred family using molecular techniques such as whole exome sequencing and further studied their effects using various bioinformatics tools.

Methods

Pedigree ascertainment

In this study, we recruited a five-generation pedigree from King Abdulaziz University Hospital (Fig. 1). All clinical information and patient features were obtained by interviewing family members and reviewing the medical records of the hospital. Before the start of the study, ethical approval was obtained from the Ethics Committee of the Center of Excellence in Genomic Medicine Research (07-CEGMR-Bioeth-2020). Moreover, we obtained written informed consent from all the study participants, according to the Declaration of Helsinki. The proband (V-3) was a male, one year old, born full-term via normal vaginal delivery. The patient had various cardiac anomalies, including pulmonic valve stenosis, atrial septal defect, small patent ductus arteriosus, and a patent foramen ovale. The patient had difficulty breathing, showed signs of apnea and cyanosis, was admitted to the neonatal intensive care unit at the time of birth because of bronchopneumonia and poor respiration, and was intubated and linked to mechanical ventilation. Later in life, various issues were observed, such as hypoxic-ischemic encephalopathy, weak muscles, global developmental delay, poor feeding, mild to moderate intellectual disability, and speech problems. Unfortunately, the child died later at approximately 18 months of age because of multiple cardiac and other complications. The proband’s older brother (V-2) also exhibited overlapping symptoms and died due to similar complications.

Blood sampling and DNA extraction

To investigate the molecular causes of the disease, we collected blood samples from the patient and his parents in EDTA tubes, immediately transferred them to the Center of Excellence in Genomic Medicine Research, and stored them at 4 °C for further investigation. We extracted DNA using a DNeasy® Blood and Tissue Kit (QIAGEN) according to the manufacturer’s protocol. The quality and quantity of DNA were measured using a spectrophotometer (Nanodrop 2000; https://www.thermofisher.com/order/catalog/product/ND2000CLAPTOP) ensuring an A260/A280 ratio of 1.7–2.0. To further confirm DNA integrity, we performed horizontal gel electrophoresis on a 1% agarose gel using SYBR Safe dye (Thermo Fisher, USA).

Whole exome sequencing

After ensuring the quality and quantity of DNA, we proceeded to the next step of the whole exome sequencing (WES) experiment to identify pathogenic mutations. WES was performed using a NovaSeq6000 machine. Sequencing samples were processed according to the manufacturer’s protocol (SureSelect_v6 Agilent, USA). For the template we used 50 ng of genomic DNA for each sample. For purity and quality, we used 260/280 absorbance ratio values of 1.7–2.0 as the standard. To prepare the libraries, we used Illumina DNA Prep with the following enrichment protocol:

1) Tagment Genomic DNA. This step uses enrichment bead-linked transposomes (eBLT) to fragment and tag DNA with adapter sequences to target the DNA.

2) Post Tagmentation Clean Up. In this step, we washed the adapter-tagged DNA on the eBLT. 3) Amplified tagged DNA. We used PCR to amplify tagmented DNA.

4) Clean Up Libraries. In this study, we purified the amplified libraries using a double-sided bead purification system.

5) Pool Libraries. Subsequently, we pooled 12 DNA libraries into a single pool.

6) Hybridize Probes. In this step, DNA-targeted regions were bonded to capture probes.

7) Capture Hybridized Probes. Streptavidin Magnetic Beads (SMB) were used to capture probes for the hybridization of the targeted DNA regions.

8) The enriched library was amplified. PCR was performed to amplify the enriched libraries. 9) Clean Up Amplified Enriched Library. AMPure XP Beads were used to purify and remove unwanted products from the enriched libraries. To assess the quantity of the enriched libraries we ran 1 µl of the sample using the Qubit dsDNA BR Assay Kit, and library yields of ≥ 3 ng/µl were considered good. To assess the quality of the enriched libraries, we used an Agilent Technology 2100 Bioanalyzer with a High Sensitivity DNA kit, and if the fragment size was approximately 350 bp and the overall range from ~200 bp to ~1000 bp was found to be suitable for further processing. Finally, the libraries were diluted to 2 nM starting concentration and 175-185 pM for final loading concentrations and run on NovaSeq6000 to obtain sequencing data.

Raw sequencing data were analyzed using Aziz Supercomputer. Initially, the base calls from BCL files were converted to FASTQ format by “bcl2fastq” and then pre-processed FASTQ files, first removed adaptors, second, low-quality bases were trimmed “<20”, third, further filtering short reads “<20 bp”. The filtered generated FASTQ files were then mapped and aligned to the reference human genome sequence UCSC hg19 version using BWA mapping software (version 0.7.12) to generate the SAM file format, which was further converted to the BAM file format (compressed binary version of SAM). The resultant allele frequency variants were annotated using ANNOVAR tools and base quality score recalibration (BQSR) was performed. GATK-Haplotype Caller was used to generate variants in the vcf file format with approximately 40 thousand variants which were further reduced to approximately 3 thousand variants after further application of filtration criteria to eliminate common variants in the general population based on minor allele frequencies (MAF <0.001). The Variant Interpreter (Illumina) was used for annotation, interpretation, and detection of disease variants. The pathogenicity of the resultant variants and their disease associations were extracted from multiple databases, including ClinVar, OMIM, PubMed, and VarSome.

Sanger’s Sequencing

To confirm the pathogenic mutation in the patient and family, we performed Sanger sequencing on an ABI XL3500 sequencer using BIGDye Terminator® technology. Primers (forward and reverse) were designed using the primer3 software (https://primer3.ut.ee/). We used forward primer sequence 5’CTGGCTCTGGTATTCAGAATC3′ and reverse primer sequence 5’GAGCTGAGGATTGTCTGAACTTA3′. KDM6A gene sequence was retrieved using accession number NM_001291415. Sequencing was performed according to the manufacturer’s instructions, and the PCR products were purified using a master mix of 100% ethanol (2.5 volume), NaOAc (1:10 volume), and EDTA (1:10 volume), which were distributed and mixed in all sample tubes. The mixture was vortexed and centrifuged at 2146 x g for 45 min. Then 100 µl of 70% ethanol was added, and the mixture was vortexed and centrifuged at 2146 x g for 10 min. The supernatant was discarded and dried at 50°C for 20 min, nuclease-free water was added for elution, and the cycle sequencing step was initiated. We prepared a master mix of 10 µl as 5x sequencing buffer 2 µl, Primers (forward and reverse) 1 µl each, Big Dye Terminator 1 µl, template DNA 1 µl, and H2O 5 µl.  The samples were then run on the thermocycler for initial denaturation at 95°C for 5 min, followed by 25 cycles of 95°C for 10 s, 50°C for 5 s, 60°C for 4 min, and held at 4°C, followed by a second purification step. The second purification was the same as the first, and after drying the samples, we added 10-12 µl of Hi-Di formamide to each sample. The Hi-Di formamide samples were denatured at 95 °C for 5 min and immediately placed on ice.

To read the resultant sequencing peaks and for analysis, we used BioEdit Sequence Alignment Editor Version 7.2.5.

In Silico Predictions

To evaluate the specific genetic variant G774E in KDM6A, we used the Combined Annotation-Dependent Depletion (CADD) method. The genomic coordinates and alleles for the G774E variant in KDM6A were entered in CADD available online at https://cadd.gs.washington,  and selected human genome assembly UCSC hg19 and then submitted the input for further processing [10,11].

Because the complete three-dimensional structure of the KDM6A protein is not available, we used the SwissModel server to build homology models of both wild-type (WT) and mutant KDM6A proteins, as described previously [12] using the AlphaFold DB model of KDM6A_HUMAN (gene: KDM6A, organism: Homo sapiens (Human)) (https://swissmodel.expasy.org/templates/afdb_v4/O15550).

To evaluate the impact of G774E on the WT protein, we used the missense 3D tool (http://missense3d.bc.ic.ac.uk/missense3d/). As an additional confirmation, we used the mutant homology structure to check whether the replacement of E with G reversed the deleterious impact using the missense 3D tool [13, 14].

Furthermore, STRING software (version 12.0) was used to predict protein-protein interaction networks and perform functional enrichment analysis (https://string-db.org/). STRING provides both physical and functional associations generated through different computational predictions and interactions gathered from other databases [15].

Results

The proband was a child of one year old (12 months) and his two siblings were the result of a consanguineous marriage. The entire family showed multiple cousin marriages, suggesting the concentration of alleles and genes in the same family (Fig. 1). The proband’s elder brother also died in early childhood due to cardiac complications and multiple neurodevelopmental issues. From birth, the child had congenital heart disease, including pulmonic valve stenosis, atrial septal defect, patent ductus arteriosus, and patent foramen ovale. The child had bronchopneumonia and poor respiration, which may have resulted in hypoxic-ischemic encephalopathy due to lack of oxygen supply to the brain. With increasing age, multiple disease features have been noted, such as weak muscles, global developmental delay, poor feeding habits and intellectual disability. Unfortunately, the proband died within a few months after donating blood for genetic analysis due to cardiac complications.

Whole exome sequencing, based on filtration criteria and disease manifestations, revealed a pathogenic mutation in the KDM6A gene, where guanine at position c.2321 was substituted with adenine (c.2321G>A). Thus, in the protein, the amino acid glycine is replaced by glutamic acid at position 774 (Gly774Glu or G774E). Furthermore, the results were confirmed by Sanger sequencing (Fig. 2). The KDM6A gene is located on the X chromosome; therefore, one affected allele on the X chromosome is sufficient to cause disease in male patients. The mother was heterozygous for this allele, whereas the father was hemizygous for the normal allele. This confirmed the transfer of the disease allele from the mother to the patient.

CADD software was used to study the lethality of the substitution, and it was predicted to be “deleterious” based on the Sorting Intolerant From Tolerant (SIFT) category (SIFTcat) with a SIFT value (SIFTval) of 0.02 (Table 1). The results of this analysis provide insights into how the G774E mutation affects the KDM6A protein’s demethylase activity, its interactions with other proteins or DNA, and its overall role in epigenetic regulation. Such information is essential for elucidating the potential pathogenic effects of this mutation and its implications in related disorders. PHRED-scaled C-scores > 20 for G774E SNV in KDM6A predicted that it was in the top 1% of deleterious variants. Furthermore, we utilized other tools, such as ExAC and 1000 genome variants, but it was not found in either ExAC or 1000G.

The 3D analysis of the missense mutation (G774E) using the WT homology model of the KDM6A protein obtained using the SwissModel server (Fig. 3) indicated a clash, as the introduction of the mutant Glutamic Acid (E) at position 774 resulted in a significantly higher clash score (76.02), compared to the Wild Type Glycine (55.21). This is based on the criteria that the mutant KDM6A homology model has a MolProbity clash score ≥ 30, which is > 18 compared to the wild type (Table 1).

Furthermore, we used the STRING software to generate a protein-protein interaction network and performed functional enrichment analysis of the KDM6A gene. STRING resulted in an interaction map that clearly showed the close association between KDM6A and KMT2D (both appear to cause Kabuki syndrome, showing similar and overlapping features as our patients). Moreover, the KDM6A gene shows bonding with NKX2-5, GATA4, and TBX5, which are exclusively involved in cardiogenesis, and if this interaction is inversely affected by one of the genes, it may result in overall impaired functioning and result in reduced or non-functional heart and body organ development (Fig. 4).

Discussion

Kabuki syndrome is a very rare disorder that affects many organs and systems in newborns, and it is estimated that 1 in 32000 newborns have this disease [16,17]. To date, approximately four genes have been reported to cause Kabuki syndrome, and one gene has features similar to those of Kabuki syndrome [18]. However, pathogenic variants in the two genes, KDM6A and KMT2D, have been reported in the literature to be the most common cause of Kabuki syndrome, where 28-80% of patients suffer from CHDs [17-19]. The first gene discovered to cause Kabuki syndrome was KMT2D (MLLe), located on chromosome 12q13.12 [20]. The second most frequently involved gene is KDM6A, which is present on the X chromosome at location Xp11.3 [21]. KDM6A and KMT2D encode proteins involved in chromatin transcription; KDM6A demethylates H3K27 (lysine-27 histone-3), and KMT2D methylates H3K4 (lysine-4 histone-3) [2, 22]. Pathogenic mutations in these two genes lead to impaired histone methylation and result in the malfunctioning of multiple organ systems [23, 24]. KMT2D gene mutations are mostly transmitted in an autosomal dominant pattern and are responsible for most cases of Kabuki syndrome (KS). The KDM6A gene, located on the X-chromosome, is transmitted from the unaffected carrier female to the patient, where the father mostly shows a normal allele, and many patients have de novo mutations rather than those inherited from the parents [15, 16]. In our case, the c. 2321G > A mutation was transmitted through the carrier mother, who had an affected allele, while the father had a normal allele (Fig. 2).

Faundes et al. (2021), in one of the major studies that involves 80 patients, analyzed the molecular basis of heterozygous and hemizygous gene alterations in the KDM6A gene and found that patients with protein-altering variations (PAVs) had shorter birth spans than those with protein termination variants (PTVs). Patients with protein-altering variants (PAVs) have shorter birth lengths than those with protein-truncating variants (PTVs). They found that intellectual disability and impaired central nervous system were more common in PTVs than in PAVs (97.6% vs. 80% and 71.4% vs. 28.6%, respectively), although the differences were not statistically significant. Hence, they concluded that PTVs had more lethal effects on the patient’s phenotype, and patients with PAVs showed a more variable phenotype [25].

The third gene involved in Kabuki syndrome is RAP1A, which is responsible for a minority of patients [26] and is located in the region of the 12-MB duplication of chromosome 1p13.1-p22.1, that is known to be involved in Kabuki syndrome [27]. The fourth gene, RAP1B, is closely related to RAP1A [26]. The last gene that is considered to cause overlapping features, such as Kabuki syndrome with nodular heterotopia, is HNRNPK [28, 29].

The major cardiac features of Kabuki syndrome include atrial and ventricular septal defects, double-outlet right ventricle, aortic coarctation, pulmonary stenosis, dysplastic mitral valve, and tetralogy of Fallot [18, 30-32]. In our study, we found overlapping features, as described in the literature, in our studied patient, who also showed various cardiac anomalies, including pulmonic valve stenosis, atrial septal defect, small patent ductus arteriosus, and patent foramen ovale.

In addition to heart features, patients with Kabuki syndrome show developmental delay, poor and difficult feeding, various growth defects, renal impairment, skeletal deformities, and distinct facial features such as sparse eyebrows, eversion of the lateral third of the lower eyelids, long palpebral fissures, and large everted ears [18, 32-34].  Our patient also showed global developmental delay, poor and difficult feeding, mild to moderate intellectual disability, speech problems, and difficulty in breathing.

In our study using STRING analysis, we found that KDM6A interacts closely with KMT2D, while it is also associated with NKX2-5, GATA4, and TBX5, which are important CHD-associated transcription factors that have been identified in the majority of CHD patients and show multiple cardiac structural and functional impairments [5]. These transcription factors are responsible for the transcriptional activation and regulation of cardiac-specific genes encoding cardiac-specific regulatory or structural proteins [35, 36].

NKX2-5 is one of the earliest transcription factors expressed in cardiac cell differentiation, atrial compartment formation, atrioventricular outflow, and conduction system formation [37, 38]. NKX2-5 is located on chromosome 5q35.1, with two coding exons, and encodes a protein of 324 amino acids. Many studies using linkage analysis and single-candidate gene research have investigated the roles of NKX2-5 and GATA4 in heart disease in large cohorts [39]. GATA4 transcription factor is expressed in almost all types of cardiac cells and is required for cardiac differentiation, proliferation, morphogenesis, heart tube formation, and proepicardium generation [40-43]. GATA4 is located on chromosome 8p23.1-p22, has seven exons, and encodes a protein of 442 amino acids. GATA4 is essential for normal cardiogenesis because of its significance during embryonic development and its continued expression in the adult myocardium [44]. TBX5 (T-box transcription factor 5) is a member of the T-box transcription factor family and is involved in forelimb and cardiac development. Pathogenic mutations in this gene impair the cardiac septum and cardiac conduction system, resulting in Holt–Oram syndrome [45]. TBX5 gene is located on chromosome 12q24.21, contains nine exons, and encodes a protein of 518 amino acids with a T-box domain of 180 amino acids located between amino acid residues 56 and 236 [46].

KDM6A, also known as UTX, is a histone demethylase that plays a crucial role in epigenetic regulation and embryonic development [47]. This gene, located on the X chromosome, encodes a protein that specifically demethylates di- and trimethylated lysine 27 on histone H3 (H3K27me2/3), thereby promoting gene activation. The significance of KDM6A in human development and disease has been underscored by its association with Kabuki syndrome, a rare genetic disorder characterized by distinctive facial features, intellectual disability, and various congenital anomalies [48]. KDM6A belongs to the Jumonji C (JmjC) domain-containing histone demethylases, which are essential for maintaining the proper epigenetic landscape in cells. By removing methyl groups from H3K27, KDM6A counteracts the repressive effects of Polycomb Repressive Complex 2 (PRC2), which is responsible for depositing methyl marks [47]. This dynamic interplay between histone methylation and demethylation is critical for the regulation of gene expression during development and in adult tissue.  Kabuki syndrome is a heterogeneous condition caused by mutations in several genes, with KDM6A being one of the primary causative genes, along with KMT2D [49]. Pathogenic variants in KDM6A account for approximately 5-8% of Kabuki syndrome cases, particularly in females, owing to its X-linked inheritance pattern [50]. The identification of KDM6A mutations in patients with KS has provided valuable insights into the molecular mechanisms underlying this complex developmental disorder.  The clinical presentation of Kabuki syndrome is highly variable, with affected individuals displaying a wide range of symptoms and severity. Common features include distinctive facial characteristics (such as elongated palpebral fissures, arched eyebrows, and large ears), skeletal abnormalities, cardiac defects, immunological problems, and developmental delays [48, 51]. This variability in clinical presentation may be attributed to the diverse roles of KDM6A in the regulation of gene expression across different tissues and developmental stages.  The involvement of KDM6A in Kabuki syndrome emphasizes the critical role of epigenetic regulation in human development. As a histone demethylase, KDM6A influences gene expression patterns by modifying the chromatin structure, which is essential for proper embryonic development and tissue-specific gene regulation [52]. Disruption of KDM6A function can lead to aberrant gene expression, resulting in diverse clinical manifestations observed in patients with Kabuki syndrome. Recent studies have demonstrated that KDM6A is involved in various developmental processes, in addition to those affected by Kabuki syndrome. For example, it has been implicated in cardiac development, neurogenesis and hematopoiesis [48]. The wide-ranging effects of KDM6A on multiple organ systems highlight its importance as a master regulator of cellular development and differentiation.  Furthermore, KDM6A plays a role in the biology of cancer. Mutations or dysregulation of KDM6A have been associated with various cancers, including bladder cancer, renal cell carcinoma, and acute myeloid leukemia [52]. These findings suggest that KDM6A may function as a tumor suppressor in certain contexts, further emphasizing its significance in the maintenance of cellular homeostasis. Understanding the role of KDM6A in Kabuki syndrome has not only advanced our knowledge of the etiology of the disorder but has also opened new avenues for potential therapeutic interventions in Kabuki syndrome. Research on KDM6A function and its downstream targets may lead to the development of targeted treatments for Kabuki syndrome and other related developmental disorders [48]. Current therapeutic approaches to Kabuki syndrome primarily focus on managing individual symptoms and providing supportive care to patients. However, the identification of KDM6A as a causative gene has sparked interest in developing targeted therapies. One potential avenue of research involves exploring the use of epigenetic modifiers to compensate for the loss of KDM6A function. For example, inhibitors of histone deacetylases (HDACs) or DNA methyltransferases (DNMTs) could potentially be used to alter the epigenetic landscape and partially restore normal gene expression patterns in affected individuals [50]. Additionally, advances in gene therapy and genome editing techniques, such as CRISPR-Cas9, offer promising possibilities for correcting KDM6A mutations and modulating their expression in affected tissues. Although these approaches are still in the early stages of development, they represent exciting potential avenues for future therapeutic interventions in Kabuki syndrome and other disorders caused by epigenetic dysregulation [53].

In our study, the links between KDM6A and other cardiac-related genes can potentially address important questions regarding cardiac function and structural development, along with neurodevelopmental and other disorders [54-56]. The missense mutation at position 774 in the KDM6A protein, resulting in the substitution of glycine with glutamic acid, has significant implications for the development of Kabuki Syndrome. This mutation affects the catalytic domain of KDM6A, a histone demethylase that is crucial for epigenetic regulation and gene expression. The replacement of a small, non-polar amino acid (glycine) with a larger, negatively charged amino acid (glutamic acid) likely alters the protein structure and function in multiple ways.  Primarily, the introduction of a negatively charged residue in place of a neutral one may disrupt the local electrostatic environment within the protein. This alteration could potentially affect the protein’s ability to interact with its substrates or cofactors, possibly modifying its catalytic efficiency. Furthermore, the larger size of glutamic acid compared to that of glycine may introduce steric hindrance, further compromising protein function.  Structural modifications induced by this mutation may impair the enzyme’s capacity to remove methyl groups from histone proteins, a process critical for regulating gene expression. Histone demethylation is a key mechanism of epigenetic control that facilitates dynamic changes in chromatin structure and accessibility. Disruption of this process can lead to aberrant gene expression patterns throughout the genome.  The implications of this mutation extend beyond the immediate biochemical effects of KDM6A protein. The altered epigenetic landscape resulting from impaired histone demethylation can have far-reaching effects on cellular function and development. Genes that are normally tightly regulated during embryonic and postnatal development may be inappropriately expressed or silenced, leading to a cascade of developmental abnormalities.  These disruptions in normal developmental processes contribute to the characteristic features of Kabuki Syndrome. Hence, the missense mutation (G774E) in KDM6A serves as a critical example of how subtle molecular changes can have profound effects on human health and development. Continued research on this and similar mutations is essential for advancing our understanding of epigenetic regulation and developing targeted therapies for Kabuki Syndrome and related disorders.

In conclusion, we identified a novel mutation in the KDM6A gene (G774E) responsible for Kabuki syndrome. The patient exhibited various cardiac abnormalities, breathing issues, hypoxic-ischemic encephalopathy, weak muscles, global developmental delay, poor feeding, mild to moderate intellectual disability, and speech problems. Furthermore, in the future, large cohorts of patient recruitment are needed, in combination with next-generation sequencing technologies combined with animal and cell models, to shed light on molecular insights into this disease and enhance our awareness and knowledge for the development of molecular/genomic therapies and to improve clinical practice.

Conclusion

Statement & Declarations

Conflict of Interest

The authors declare that they have no competing interests

Author Contributions

Mahmood Rasool, Sherin Bakhashab and Isam M. Abu Zeid designed and supervised the study. Majed Alsulami, Ayat Mohammed Shorbaji, Loubna Siraj Mira and Mohammad Basabrain performed experiments, Mohamed Nabil Alama provided samples. Absarul Haque wrote and edited manuscript. Sajjad Karim, Hisam N. Altayb and Peter Natesan Pushparaj performed in silico analysis.

Acknowledgment

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through project number (574).

References

  1. Bouma BJ, Mulder BJ. Changing Landscape of Congenital Heart Disease.  Circulation research, (2017); 120(6): 908-922.
  2. Sen R, Lencer E, Geiger EA, Jones K, Shaikh TH, Artinger KB. The role of Kabuki Syndrome genes KMT2D and KDM6A in development: Analysis in Human sequencing data and compared to mice and zebrafish. bioRxiv, (2020); 03.024646: preprint.
  3. Majeed-Saidan MA, Atiyah M, Ammari AN, AlHashem AM, Rakaf MS, et al. Patterns, prevalence, risk factors, and survival of newborns with congenital heart defects in a Saudi population: a three-year, cohort case-control study. Journal of Congenital Cardiology, (2019); 3: 2.
  4. Pierpont ME, Brueckner M, Chung WK, Garg V, Lacro RV, et al. Genetic Basis for Congenital Heart Disease: Revisited: A Scientific Statement From the American Heart Association. Circulation, (2018); 138(21): e653-e711.
  5. Williams K, Carson J, Lo C. Genetics of Congenital Heart Disease. Biomolecules, (2019); 9(12): 879.
  6. Lage K, Greenway SC, Rosenfeld JA, Wakimoto H, Gorham JM, et al. Genetic and environmental risk factors in congenital heart disease functionally converge in protein networks driving heart development. Proceedings of the National Academy of Sciences of the United States of America, (2012); 109(35): 14035-40.
  7. Li Y, Klena NT, Gabriel GC, Liu X, Kim AJ, Lemke K, et al. Global genetic analysis in mice unveils central role for cilia in congenital heart disease. Nature, (2015); 521(7553): 520-4.
  8. Zaidi S, Brueckner M. Genetics and Genomics of Congenital Heart Disease. Circulation research, (2017); 120(6):923-940.
  9. Meller CH, Grinenco S, Aiello H, Córdoba A, Sáenz-Tejeira MM, et al. Congenital heart disease, prenatal diagnosis and management. Archivos argentinos de pediatria, (2020); 118(2): e149-e161.
  10. Rentzsch P, Witten D, Cooper GM, Shendure J, Kircher M. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Research, (2019); 47(D1): D886-D894.
  11. Schubach M, Maass T, Nazaretyan L, Röner S, Kircher M. CADD v1.7: using protein language models, regulatory CNNs and other nucleotide-level scores to improve genome-wide variant predictions. Nucleic Acids Research, (2024); 52(D1): D1143-D1154.
  12. Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Research, (2018); 46(W1): W296-W303.
  13. Ittisoponpisan S, Islam SA, Khanna T, Alhuzimi E, David A, Sternberg MJE. Can Predicted Protein 3D Structures Provide Reliable Insights into whether Missense Variants Are Disease Associated? Journal of Molecular Biology, (2019); 431(11):2197-2212.
  14. Pennica C, Hanna G, Islam SA, Sternberg MJE, David A. Missense3D-PPI: A Web Resource to Predict the Impact of Missense Variants at Protein Interfaces Using 3D Structural Data. Journal of Molecular Biology, (2023); 435(14): 168060.
  15. Szklarczyk D, Kirsch R, Koutrouli M, Nastou K, Mehryary F, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Research, (2023); 51(D1): D638-D646.
  16. Salguero MV, Chan K, Greeley SAW, Dyamenahalli U, Waggoner D, et al. Novel KDM6A Kabuki Syndrome Mutation with Hyperinsulinemic Hypoglycemia and Pulmonary Hypertension Requiring ECMO. Journal of the Endocrine Society, (2022); 6(4): bvac015.
  17. Adam MP, Banka S, Bjornsson HT, Bodamer O, Chudley AE, et al. Kabuki Syndrome Medical Advisory Board. Kabuki syndrome: international consensus diagnostic criteria. Journal of Medical Genetics, (2019); 56(2): 89-95.
  18. Digilio MC, Gnazzo M, Lepri F, Dentici ML, Pisaneschi E, et al. Congenital heart defects in molecularly proven Kabuki syndrome patients. American Journal of Medical Genetics Part A, (2017); 173(11): 2912-2922.
  19. Cuvertino S, Hartill V, Colyer A, Garner T, Nair N, et al. A restricted spectrum of missense KMT2D variants cause a multiple malformations disorder distinct from Kabuki syndrome. Genetics in Medicine, (2020); 22(5): 867-877.
  20. Ng SB, Bigham AW, Buckingham KJ, Hannibal MC, McMillin MJ, et al. Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome. Nature Genetics, (2010); 42(9): 790-3.
  21. Lederer D, Grisart B, Digilio MC, Benoit V, Crespin M, et al. Deletion of KDM6A, a histone demethylase interacting with MLL2, in three patients with Kabuki syndrome. American Journal of Human Genetics, (2012); 90(1): 119-24.
  22. Ali M, Hom RA, Blakeslee W, Ikenouye L, Kutateladze TG. Diverse functions of PHD fingers of the MLL/KMT2 subfamily. Biochimica et biophysica acta, (2014); 1843(2): 366-71.
  23. Yap KL, Johnson AEK, Fischer D, Kandikatla P, Deml J, et al. Congenital hyperinsulinism as the presenting feature of Kabuki syndrome: clinical and molecular characterization of 9 affected individuals. Genetics in Medicine, (2019); 21(1): 233-242.
  24. Micale L, Augello B, Maffeo C, Selicorni A, Zucchetti F, et al. Molecular analysis, pathogenic mechanisms, and readthrough therapy on a large cohort of Kabuki syndrome patients. Human Mutation, (2014); 35(7): 841-50.
  25. Faundes V, Goh S, Akilapa R, Bezuidenhout H, Bjornsson HT, et al. Clinical delineation, sex differences, and genotype-phenotype correlation in pathogenic KDM6A variants causing X-linked Kabuki syndrome type 2. Genetics in Medicine, (2021); 23(7): 1202-1210.
  26. Bögershausen N, Tsai IC, Pohl E, Kiper PÖ, Beleggia F, et al. RAP1-mediated MEK/ERK pathway defects in Kabuki syndrome. Journal of Clinical Investigation, (2015); 125(9): 3585-99.
  27. Lo IF, Cheung LY, Ng AY, Lam ST. Interstitial Dup(1p) with findings of Kabuki make-up syndrome. American Journal of Medical Genetics, (1998); 78(1): 55-7.
  28. Au PYB, You J, Caluseriu O, Schwartzentruber J, Majewski J, et al. GeneMatcher aids in the identification of a new malformation syndrome with intellectual disability, unique facial dysmorphisms, and skeletal and connective tissue abnormalities caused by de novo variants in HNRNPK. Human Mutation, (2015); 36(10): 1009-1014.
  29. Lange L, Pagnamenta AT, Lise S, Clasper S, Stewart H, et al. A de novo frameshift in HNRNPK causing a Kabuki-like syndrome with nodular heterotopia. Clinical Genetics, (2016); 90(3): 258-62.
  30. Cocciadiferro D, Augello B, De Nittis P, Zhang J, Mandriani B, et al. Dissecting KMT2D missense mutations in Kabuki syndrome patients. Human Molecular Genetics, (2018); 27(21): 3651-3668.
  31. Gažová I, Lengeling A, Summers KM. Lysine demethylases KDM6A and UTY: The X and Y of histone demethylation. Molecular genetics and metabolism, (2019); 127(1): 31-44.
  32. Yap CS, Jamuar SS, Lai AHM, Tan ES, Ng I, et al. Identification of KMT2D and KDM6A variants by targeted sequencing from patients with Kabuki syndrome and other congenital disorders. Gene, (2020); 731: 144360.
  33. Kawame H, Hannibal MC, Hudgins L, Pagon RA. Phenotypic spectrum and management issues in Kabuki syndrome. The Journal of pediatrics, (1999); 134(4): 480-5.
  34. Kluijt I, van Dorp DB, Kwee ML, Toutain A, Keppler-Noreuil K, et al. Kabuki syndrome – report of six cases and review of the literature with emphasis on ocular features. Ophthalmic genetics, (2000); 21(1): 51-61.
  35. Cao C, Li L, Zhang Q, Li H, Wang Z, et al. Nkx2.5: a crucial regulator of cardiac development, regeneration and diseases. Frontiers in cardiovascular medicine, (2023); 10: 1270951.
  36. Razzaq SS, Khan I, Naeem N, Salim A, Begum S, Haneef K. Overexpression of GATA binding protein 4 and myocyte enhancer factor 2C induces differentiation of mesenchymal stem cells into cardiac-like cells. World journal of stem cells, (2022); 14(9): 700-713.
  37. Peña-Martínez EG, Rivera-Madera A, Pomales-Matos DA, Sanabria-Alberto L, Rosario-Cañuelas BM, et al. Disease-associated non-coding variants alter NKX2-5 DNA-binding affinity. Biochimica et biophysica acta. Gene regulatory mechanisms, (2023); 1866(1): 194906.
  38. Henry S, Kokity L, Pirity MK. Polycomb protein RYBP activates transcription factor Plagl1 during in vitro cardiac differentiation of mouse embryonic stem cells. Open Biology, (2023); 13(2): 220305.
  39. Su W, Zhu P, Wang R, Wu Q, Wang M, et al. Congenital heart diseases and their association with the variant distribution features on susceptibility genes. Clinical Genetics, (2017); 91(3): 349-354.
  40. Fan D, Pang S, Chen J, Shan J, Cheng Q, Yan B. Identification and functional study of GATA4 gene regulatory variants in atrial septal defects. BMC Cardiovascular Disorders, (2021); 21(1): 321.
  41. Tremblay M, Sanchez-Ferras O, Bouchard M. GATA transcription factors in development and disease. Development, (2018); 145(20): dev164384.
  42. Zhou L, Liu J, Xiang M, Olson P, Guzzetta A, et al. Gata4 potentiates second heart field proliferation and Hedgehog signaling for cardiac septation. Proceedings of the National Academy of Sciences of the United States of America, (2017); 114(8): E1422-E1431.
  43. Watt AJ, Battle MA, Li J, Duncan SA. GATA4 is essential for formation of the proepicardium and regulates cardiogenesis. Proceedings of the National Academy of Sciences of the United States of America, (2004); 101(34): 12573-8.
  44. Jumppanen M, Kinnunen SM, Välimäki MJ, Talman V, Auno S, et al. Synthesis, Identification, and Structure-Activity Relationship Analysis of GATA4 and NKX2-5 Protein-Protein Interaction Modulators. Journal of medicinal chemistry, (2019); 62(17): 8284-8310.
  45. Steimle JD, Moskowitz IP. TBX5: A Key Regulator of Heart Development. Current topics in developmental biology, (2017); 122:195-221.
  46. Basson CT, Bachinsky DR, Lin RC, Levi T, Elkins JA, et al. Mutations in human TBX5 [corrected] cause limb and cardiac malformation in Holt-Oram syndrome. Nature Genetics, (1997); 15(1): 30-5.
  47. Agger K, Cloos PA, Christensen J, Pasini D, Rose S. et al. UTX and JMJD3 are histone H3K27 demethylases involved in HOX gene regulation and development. Nature, (2007); 449(7163): 731-4.
  48. Shpargel KB, Quickstad G. SETting up the genome: KMT2D and KDM6A genomic function in the Kabuki syndrome craniofacial developmental disorder. Birth Defects Research, (2023); 115(20): 1885-1898.
  49. Miyake N, Koshimizu E, Okamoto N, Mizuno S, Ogata T, et al. MLL2 and KDM6A mutations in patients with Kabuki syndrome. American Journal of Medical Genetics Part A, (2013); 161A(9): 2234-43.
  50. Banka S, Lederer D, Benoit V, Jenkins E, Howard E, et al. Novel KDM6A (UTX) mutations and a clinical and molecular review of the X-linked Kabuki syndrome (KS2). Clinical Genetics, (2015); 87(3): 252-8.
  51. Van Laarhoven PM, Neitzel LR, Quintana AM, Geiger EA, Zackai EH, et al. Kabuki syndrome genes KMT2D and KDM6A: functional analyses demonstrate critical roles in craniofacial, heart and brain development. Human Molecular Genetics, (2015); 24(15): 4443-53.
  52. Van der Meulen J, Speleman F, Van Vlierberghe P. The H3K27me3 demethylase UTX in normal development and disease. Epigenetics, (2014); 9(5): 658-68.
  53. Ma W, Fang H, Pease N, Filippova GN, Disteche CM, Berletch JB. Sex-biased and parental allele-specific gene regulation by KDM6A. Biol Sex Differ, (2022); 13(1): 40.
  54. Miyake N, Mizuno S, Okamoto N, Ohashi H, Shiina M, et al. KDM6A point mutations cause Kabuki syndrome. Human Mutation, (2013); 34(1): 108-10.
  55. Murabit FA, Etwebi KA. First Clinical Case report of a Libyan Girl with Kabuki Syndrome and Literature review. Libyan Journal of Medical Research, (2020); 14(2):70–80.
  56. Liu S, Hong X, Shen C, Shi Q, Wang J, et al. Kabuki syndrome: a Chinese case series and systematic review of the spectrum of mutations. BMC Medical Genetics, (2015); 16: 26.
Article Sections

Edited by

Zia-ur-RahmanUniversity of Lahore, Pakistan

Reviewed by

Mohammad Zahid MustafaDepartment of Inserm, Unité 1193, Université Paris-Saclay, Villejuif 94800, France
Abdul QadirDepartment of Plant Pathology, Agriculture College Quetta, Pakistan

Figures

Tables

Editors & Reviewers

Edited by

Zia-ur-RahmanUniversity of Lahore, Pakistan

Reviewed by

Mohammad Zahid MustafaDepartment of Inserm, Unité 1193, Université Paris-Saclay, Villejuif 94800, France
Abdul QadirDepartment of Plant Pathology, Agriculture College Quetta, Pakistan

Figures

Tables

Share this article: