Kabuki syndrome (KS), first described in 1981 by Niikawa et al and Kuroki et al, is a rare congenital mental retardation condition with an incidence ranging from 1/32,000 to 1/86,000. Patients with KS typically present with characteristic facial features including long palpebral fissures with ectropion of eversion of the lateral third of the lower eyelids, arched eyebrows, a depressed nasal tip, and prominent ears.
KS is also known as Kabuki make-up syndrome due to the resemblance of the typical facial dysmorphia to actors making up in Japanese traditional theater. Additional manifestations that are frequently associated with KS include a high incidence of intellectual disability, cardiovascular diseases, short stature, and musculoskeletal abnormalities.
Both KMT2D and KDM6A pathogenic variants were associated with KS, resulting in KS type 1 (KS1, OMIM #147920) and KS type 2 (KS2, OMIM #300867), respectively. Approximately 55% to 80% KS patients were attributed to de novo or autosomal dominant KMT2D (12q13.12) variants. Loss-of-function variants of KDM6A (Xq11.3) have been identified in several KMT2D variant-negative patients since Lederer et al reported 3 KS cases in 2012. The KDM6A contains 29 exons and encodes the lysine-specific demethylase 6A (KDM6A), which demethylates the H3K27me3/me2 mark specifically.[7–9] The H3K27me3/me2 mark was notified as a DNA methylation-independent imprinting mark recently.[7–9] Moreover, KDM6A and KMT2D encode methyltransferase activating signal cointegrator-2-containing complex that contains various unique components, contributing to equilibrating the active/repressive histone methylations.[10–13]
Since deletions of KDM6A were initially identified in 3 KMT2D-variant-negative KS patients by array comparative genomic hybridization analysis, 45 X-linked KS2 patients have been reported in various countries during the last 6 years,[3–6,12,14–20] but only 1 case was reported in Chinese population. In terms of reported KDM6A variants, KS2 seems more common in Europe (51.1%, 23/45), where the first 3 KS2 patients were identified.
The present study identified a novel de novo variant of the KDM6A gene in a Chinese patient with characteristic facial features by a combination of facial dysmorphology analysis and NGS-based multigene panel approach. The consequence of the variant was analyzed through RNA assay.
Participants and methods
Participants and facial dysmorphology analysis
Genomic DNA of all the participants, including the patient, her parents and 215 non-relative healthy controls matched for ethnicity, were extracted from their peripheral blood samples according to standard procedures. Face2Gene (http://www.face2gene.com), a suite of phenotyping application that contributed to genetic disease evaluation, was used for facial dysmorphology analysis. A cluster of syndromes was suggested after information of the proband (her photo, date of visit, date of birth, gender, height, weight, etc.) were filled out online.
All the participants were recruited from the outpatient department of the International Peace Maternity and Child Health Hospital (IPMCH), Shanghai Jiao Tong University School of Medicine, China. This study was approved by the Ethics Review Committee of IPMCH (No. GKLW2017–149) and performed according to the Declaration of Helsinki. Written informed consent was provided by legal guardians of the KS patient and other participants.
Targeted sequencing and data analysis
Targeted capture of 2181 genes and NGS had been described in detail previously.[22,23] The genomic DNA was isolated from peripheral blood sample and sheared into segments in size from 200 to 250 bp. Sequencing of the captured target regions was performed on Illumina HiSeq2500 Analyzers (Illumina, San Diego, CA, USA). Alignment against the human reference genome hg19/GRCh37 was performed with Burrows–Wheeler Aligner (BWA; version 0.7.12) (http://bio-bwa.sourceforge.net/). Variant calling and annotation were achieved using the Genome Analysis Toolkit (GATK; version 3.5) (https://software.broadinstitute.org/gatk/) and ANNOVAR, respectively. The resulting variants were prioritized with our in-house developed pipeline (MultiOmics One). The interpretation and classification of variants were based on the American College of Medical Genetics and Genomics and the Association for Molecular Pathology Guideline.
A total of 3.88 Gb raw data were obtained with an average coverage of 189.04× and 97.3% of bases in the defined target reached the coverage of at least 10×. Subsequent variants were filtered using our in-house MultiOmics One pipeline and intersected with genes associated with the suggested syndromes according to Face2Gene analysis.
Reverse-transcription PCR (RT-PCR), Sanger sequencing and identity testing
RNA extracted from peripheral blood mononuclear cells of the proband and a healthy female control were reversed to cDNA with a PrimeScript™ RT reagent Kit with gDNA Eraser (NO. RR047Q, Takara, Japan). Nonsense-mediated mRNA decay (NMD) was subsequently detected by polymerase chain reaction (PCR) amplification and Sanger sequencing following the standard methods. Primer sequences for verifying KDM6A variant were designed with the Primer 3 Web (http://primer3.ut.ee/). Sequence variant was named based on the current Human Genome Variation Society nomenclature (version 15.11) (http://varnomen.hgvs.org/) and checked with the Mutalyzer (https://mutalyzer.nl/). The short tandem repeat (STR) typing method was applied for identity testing, as described previously. All primers were listed in supplemental materials (Additional Table 1, http://links.lww.com/JR9/A1).
Clinical features of the patient
The patient, a Chinese 24-year-old married woman, came to the outpatient department of IPMCH for pre-pregnancy consultation (Fig. 1A, patient II-1). She had moderate learning difficulties, short stature (140 cm, <1st centile) and skeletal malformations including brachydactyly and hip dysplasia, leading to unsteady walking with unequal legs. She presented with characteristic facial features, including arched and laterally sparse eyebrows, long palpebral fissure, left eye exotropia, eversion of lateral third lower eyelid, long eyelashes, depressed tip of nose with short columella, and dental anomalies (large upper incisors, irregular teeth and malocclusion) (Fig. 1B–E). Her parents and husband, showed in the pedigree chart, were normal in these physical examinations and they denied similar patients in their family.
Candidate syndromes revealed by facial dysmorphology analysis
Considering characteristic facial features of the patient, the phenotyping application Face2Gene was applied to match the molecularly confirmed syndromes by analyzing facial features of the proband. The top four syndromes suggested by Face2Gene were KS, Mucopolysaccharidoses, Prader-Willi syndrome and odontoonychodermal dysplasia (Fig. 2A). Remarkably, the proband has a high gestalt similarity to composite public photos of KS patients (Fig. 2B).
Identification of a novel KDM6A pathogenic variant using NGS-based targeted sequencing
The proband was heterozygous for a novel nonsense variant NM_021140.2(KDM6A):c.3521G>A, predicting to introduce a premature stop codon (p.Trp1174∗) of the KDM6A protein (Fig. 3A).
Sanger sequencing and variant classification
The result of Sanger sequencing was consistent with NGS, validating the presence of NM_021140.2(KDM6A):c.3521G>A heterozygous variant in the proband (Fig. 3B). The variant was not detected in her parents. Furthermore, identity testing, using short tandem repeat (STR) loci and amelogenin, confirmed the de novo characteristics (Additional Table 2, http://links.lww.com/JR9/A1). The variant was also excluded in 1000 Genomes Project (http://browser.1000genomes.org/), Exome Aggregation Consortium (ExAC; http://exac.broadinstitute.org/), database of human single-nucleotide polymorphisms (dbSNP; http://www.ncbi.nlm.nih.gov/projects/SNP/), as well as 215 unrelated controls.
The nonsense variant in exon 24 was predicted to lead to a premature stop codon (p.Trp1174∗), which might result in the disruption of Jumonji C domain, the most functionally important part of KDM6A in catalyzing the demethylation of H3K27. More importantly, the variant was not located in the last two exons, suggesting the corresponding variant transcript could undergo the NMD.
RT-PCR analysis of KDM6A mRNA
To validate our hypothesis, RT-PCR analysis of RNA extracted from peripheral blood mononuclear cells of the proband and a healthy female control was performed with the primers in exon 22 and exon 27 (Additional Table 1, http://links.lww.com/JR9/A1). cDNA sequencing only showed the wild-type KDM6A transcripts of the proband (Fig. 3C), suggesting monoallelic expression of the wild-type allele and NMD of the variant KDM6A transcripts. Together, this variant can be categorized to pathogenic variant according to the American College of Medical Genetics and Genomics and the Association for Molecular Pathology guideline using pathogenicity very strong and pathogenicity strong evidence.
The integrated phenotype and genotype-based prioritization strategy has been shown to be highly effective for the final diagnosis of KS. However, performance of this strategy also requires detailed phenotypes and clinical experience. Facial dysmorphology technology is becoming more accessible and may help physicians narrow down candidate genetic diagnosis behind disorders associated with characteristic facial features. The Face2Gene has been limited to be applied in an Asian population. Recently, a study indicated ethnic background of the patient could influence the evaluation. In this study, Face2Gene was used to match the genetically confirmed syndromes by analyzing facial features of the proband. Notably, KS was prioritized at the top of suggested syndromes list, consistent with the final molecular diagnosis, which might reflect the face of the KS patient is significantly different from the average, general population face.
KMT2D and KDM6A are the only two genes known to cause KS, among which KDM6A variants account for 3% to 8% of KS cases. A recent study demonstrated that clinical manifestations of two subtypes of KS, KMT2D associated KS1 and KDM6A associated KS2, were so similar that it was hard to distinguish between them. The most evident distinction between KS1 and KS2 is that patients with KDM6A pathogenic variants tend to have sex-specific phenotypic differences, that is, male patients showed much more severe intellectual disability and developmental delay than females.[5,12,16] It is especially evident in the three-generation KS2 family, two inherited sons with severe phenotypes, and their mother and maternal grandmother with attenuated phenotypes. For our patient, the 24-year-old female presented with mild learning difficulties besides a typical face, short stature, and hip dysplasia. She was not aware of the aforementioned features until she considered having a baby, consistent with the previous findings that female patients were mild in intellectual manifestation. Besides, the long halluces, large central incisors and hypertrichosis features in adults and hypoglycemia in children have been more frequently described in KS2 and served as good indicators of KDM6A variants.[5,12] However, only the large central incisors were found in our patient.
To date, only 45 KS2 cases have been reported all over the world, of them 23 (51.1%) patients are from Europe, 11 (24.4%) from Asia, 5 (11.1%) from America, 1 (2.2%) from Australia and 5 (11.1%) cases lack nationality information (Additional Table 3, http://links.lww.com/JR9/A1). With regard to updated KS2 cases, a large proportion of KS type 2 patients were from Europe. Nevertheless, disparity in techniques and awareness of KS2 between countries should be taken into account for further prediction of the ethnic or racial predilection of KS2.
The evidence for genotype-phenotype correlations between KDM6A variants and clinical severity is ambiguous with respect to the limited reported variants of KDM6A. Recently, Bögershausen et al indicated that KDM6A splice-site variants are the most common variant type, whereas the updated published variants combined with our identified case showed an equality between the number of nonsense variants and splice-site variants in total (Fig. 4A). It is found that a large proportion of KDM6A variants distributed across exons 16 to 29, which is the same as the observation of Bögershausen et al. Variants in exon 6 are the most frequently observed among these single-exon-involved variants, leaving out variants in exon 17, which contained a three-generation family with four patients (Fig. 4B), suggesting a variant hot spot may locate in exon 6. However, approximately 30% of patients with KS-like phenotypes are negative for KDM6A and KMT2D, implying that there may be other candidate genes associated with KS. Genes similar to KDM6A and/or KMT2D in function or biological processes would be more likely to be new candidate genes for further consideration.
KDM6A is able to escape X-inactivation, and haploinsufficiency of KDM6A is likely to be the pathogenic mechanism of KS2. Thus, a pathogenic variant of KDM6A in females, on either the inactive or the active X chromosome, would be disease-causing. For our patient, the risk for an affected son or daughter is 50%. The successful genetic diagnosis of the proband thus could provide an opportunity for prenatal diagnosis of the family.
In conclusion, we reported a novel de novo KDM6A pathogenic variant associated with KS via integrated facial analysis and targeted sequencing. The novel pathogenic variant is the first KDM6A point variant identified in Chinese patients. Moreover, our results provide evidence that the integrated facial dysmorphology analysis and multigene panel sequencing approach is valuable for the diagnosis of the clinical and genetic variability in KS.
WS, JZ, SC, HF, YC, SL, CC and LZ carried out the experiments; JZ and HH analyzed the data; JZ made the figures; WS and JZ drafted and CX revised the paper. All authors approved the final version of the paper.
This work was supported by the National Natural Science Foundation of China (No. 81471506, 81401219, 81501276, and 81771638), the National Key Research and Development Program of China (No. 2016YFC0905103), the Shanghai Municipal Commission of Science and Technology Program, China (No. 15411966700, 15411964000, and 17411972900), the Municipal Human Resources Development Program for Outstanding Young Talents in Medical and Health Sciences in Shanghai, China (No. 2018YQ39), the Shanghai Municipal Commission of Health and Family Planning Program, China (No. 20154Y0039 and 15GWZK0701), the Shanghai Jiao Tong University Program, China (No. YG2017MS39), and the Innovation Foundation of Translational Medicine of Shanghai Jiao Tong University School of Medicine, Shanghai SJTUSM Biobank, China (No. 15ZH4011).
Conflicts of interest
The authors declare that they have no conflicts of interest.
Institutional review board statement
This study was approved by the Ethics Review Committee of IPMCH (No. GKLW2017–149) and performed according to the Declaration of Helsinki.
1. Niikawa N, Matsuura N, Fukushima Y, et al. Kabuki make-up syndrome: a syndrome of mental retardation, unusual facies, large and protruding ears, and postnatal growth deficiency. J Pediatr
2. Kuroki Y, Suzuki Y, Chyo H, et al. A new malformation syndrome of long palpebral fissures, large ears, depressed nasal tip, and skeletal anomalies associated with postnatal dwarfism and mental retardation. J Pediatr
3. Gole H, Chuk R, Coman D. Persistent hyperinsulinism in Kabuki syndrome 2: case report and literature review. Clin Pract
4. Van Laarhoven PM, Neitzel LR, Quintana AM, et al. Kabuki syndrome genes KMT2D and KDM6A: functional analyses demonstrate critical roles in craniofacial, heart and brain development. Hum Mol Genet
5. Banka S, Lederer D, Benoit V, et al. Novel KDM6A (UTX) mutations and a clinical and molecular review of the X-linked Kabuki syndrome (KS2). Clin Genet
6. Lederer D, Grisart B, Digilio MC, et al. Deletion of KDM6A, a histone demethylase interacting with MLL2, in three patients with Kabuki syndrome. Am J Hum Genet
7. Martin C, Zhang Y. The diverse functions of histone lysine methylation. Nat Rev Mol Cell Biol
8. Sengoku T, Yokoyama S. Structural basis for histone H3 Lys 27 demethylation by UTX/KDM6A. Genes Dev
9. Inoue A, Jiang L, Lu F, et al. Maternal H3K27me3 controls DNA methylation-independent imprinting. Nature
10. Goo YH, Sohn YC, Kim DH, et al. Activating signal cointegrator 2 belongs to a novel steady-state complex that contains a subset of trithorax group proteins. Mol Cell Biol
11. Agger K, Cloos PA, Christensen J, et al. UTX and JMJD3 are histone H3K27 demethylases involved in HOX gene regulation and development. Nature
12. Bögershausen N, Gatinois V, Riehmer V, et al. Mutation update for Kabuki syndrome genes KMT2D and KDM6A and further delineation of X-linked Kabuki syndrome subtype 2. Hum Mutat
13. Issaeva I, Zonis Y, Rozovskaia T, et al. Knockdown of ALR (MLL2) reveals ALR target genes and leads to alterations in cell adhesion and growth. Mol Cell Biol
14. Lederer D, Shears D, Benoit V, et al. A three generation X-linked family with Kabuki syndrome phenotype and a frameshift mutation in KDM6A. Am J Med Genet A
15. Miyake N, Koshimizu E, Okamoto N, et al. MLL2 and KDM6A mutations in patients with Kabuki syndrome. Am J Med Genet A
16. Miyake N, Mizuno S, Okamoto N, et al. KDM6A point mutations cause Kabuki syndrome. Hum Mutat
17. Lindgren AM, Hoyos T, Talkowski ME, et al. Haploinsufficiency of KDM6A is associated with severe psychomotor retardation, global growth restriction, seizures and cleft palate. Hum Genet
18. Cheon CK, Sohn YB, Ko JM, et al. Identification of KMT2D and KDM6A mutations by exome sequencing in Korean patients with Kabuki syndrome. J Hum Genet
19. Micale L, Augello B, Maffeo C, et al. Molecular analysis, pathogenic mechanisms, and readthrough therapy on a large cohort of Kabuki syndrome patients. Hum Mutat
20. Morgan AT, Mei C, Da Costa A, et al. Speech and language in a genotyped cohort of individuals with Kabuki syndrome. Am J Med Genet A
21. Yang P, Tan H, Xia Y, et al. De novo exonic deletion of KDM6A in a Chinese girl with Kabuki syndrome: a case report and brief literature review. Am J Med Genet A
22. Liu Y, Wei X, Kong X, et al. Targeted next-generation sequencing for clinical diagnosis of 561 Mendelian diseases. PLoS One
23. Zhang JY, Chen SC, Chen YY, et al. Targeted sequencing
identifies a novel SH2D1A pathogenic variant in a Chinese family: carrier screening and prenatal genetic testing. PLoS One
24. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics
25. Yang H, Wang K. Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR. Nat Protoc
26. Richards S, Aziz N, Bale S, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med
27. Untergasser A, Cutcutache I, Koressaar T, et al. Primer3–new capabilities and interfaces. Nucleic Acids Res
28. den Dunnen JT, Dalgleish R, Maglott DR, et al. HGVS recommendations for the description of sequence variants: 2016 update. Hum Mutat
29. Wildeman M, van Ophuizen E, den Dunnen JT, et al. Improving sequence variant descriptions in mutation databases and literature using the Mutalyzer sequence variation nomenclature checker. Hum Mutat
30. Lee MG, Villa R, Trojer P, et al. Demethylation of H3K27 regulates polycomb recruitment and H2A ubiquitination. Science
31. Khajavi M, Inoue K, Lupski JR. Nonsense-mediated mRNA decay modulates clinical outcome of genetic disease. Eur J Hum Genet
32. Bone WP, Washington NL, Buske OJ, et al. Computational evaluation of exome sequence data using human and model organism phenotypes improves diagnostic efficiency. Genet Med
33. Lumaka A, Cosemans N, Lulebo Mampasi A, et al. Facial dysmorphism is influenced by ethnic background of the patient and of the evaluator. Clin Genet
34. Lintas C, Persico AM. Unraveling molecular pathways shared by Kabuki and Kabuki-like syndromes. Clin Genet