Mutation spectrum of MLL2 in a cohort of kabuki syndrome patients
- Lucia Micale†1,
- Bartolomeo Augello†1,
- Carmela Fusco†1,
- Angelo Selicorni2,
- Maria N Loviglio1,
- Margherita Cirillo Silengo3,
- Alexandre Reymond4,
- Barbara Gumiero2,
- Federica Zucchetti2,
- Ester V D'Addetta1,
- Elga Belligni3,
- Alessia Calcagnì1,
- Maria C Digilio5,
- Bruno Dallapiccola5,
- Francesca Faravelli6,
- Francesca Forzano6,
- Maria Accadia7,
- Aldo Bonfante8,
- Maurizio Clementi9,
- Cecilia Daolio9,
- Sofia Douzgou10,
- Paola Ferrari11,
- Rita Fischetto12,
- Livia Garavelli13,
- Elisabetta Lapi14,
- Teresa Mattina15,
- Daniela Melis16,
- Maria G Patricelli17,
- Manuela Priolo18,
- Paolo Prontera19,
- Alessandra Renieri20,
- Maria A Mencarelli20,
- Gioacchino Scarano21,
- Matteo della Monica21,
- Benedetta Toschi22,
- Licia Turolla23,
- Alessandra Vancini24,
- Adriana Zatterale25,
- Orazio Gabrielli26,
- Leopoldo Zelante†1 and
- Giuseppe Merla†1Email author
© Micale et al; licensee BioMed Central Ltd. 2011
Received: 11 March 2011
Accepted: 9 June 2011
Published: 9 June 2011
Kabuki syndrome (Niikawa-Kuroki syndrome) is a rare, multiple congenital anomalies/mental retardation syndrome characterized by a peculiar face, short stature, skeletal, visceral and dermatoglyphic abnormalities, cardiac anomalies, and immunological defects. Recently mutations in the histone methyl transferase MLL2 gene have been identified as its underlying cause.
Genomic DNAs were extracted from 62 index patients clinically diagnosed as affected by Kabuki syndrome. Sanger sequencing was performed to analyze the whole coding region of the MLL2 gene including intron-exon junctions. The putative causal and possible functional effect of each nucleotide variant identified was estimated by in silico prediction tools.
We identified 45 patients with MLL2 nucleotide variants. 38 out of the 42 variants were never described before. Consistently with previous reports, the majority are nonsense or frameshift mutations predicted to generate a truncated polypeptide. We also identified 3 indel, 7 missense and 3 splice site.
This study emphasizes the relevance of mutational screening of the MLL2 gene among patients diagnosed with Kabuki syndrome. The identification of a large spectrum of MLL2 mutations possibly offers the opportunity to improve the actual knowledge on the clinical basis of this multiple congenital anomalies/mental retardation syndrome, design functional studies to understand the molecular mechanisms underlying this disease, establish genotype-phenotype correlations and improve clinical management.
Kabuki syndrome (KS, MIM #147920), also known as Niikawa-Kuroki syndrome, is a rare, multiple congenital anomalies/mental retardation syndrome characterized by a peculiar face, which is defined by long palpebral fissures with eversion of the lateral third of the lower eyelids, short columella with a broad and depressed nasal tip, prominent ears, and a cleft or high-arched palate. Additional features include short stature, skeletal, visceral and dermatoglyphic abnormalities, cardiac anomalies, and immunological defects [1, 2]. Kabuki syndrome has an incidence of 1 in 32,000, likely largely underestimated . The vast majority of reported cases are sporadic. After initial and controversial data that associated this condition to chromosomal rearrangement [4, 5], mutations in the MLL2 gene identified the underlying cause of Kabuki syndrome in approximately 72% of affected individuals [6, 7]. The encoded MLL2 protein is a member of the Mixed Lineage Leukemia (MLL) family of histone methyl transferases (HMT). The MLL proteins (MLLs) are part of the SET (Su(var)3-9, Enhancer-of-zeste, Trithorax) family of proteins . The highly conserved SET domain of MLLs confers histone methyltransferase activity, which is the core function of HMTs. MLLs are important in the epigenetic control of active chromatin states . They act as transcriptional co-activators and are involved in embryogenesis and development through, for example, regulation of the expression of the HOX genes and their interaction with nuclear receptors [10, 11].
The MLL2 gene encodes a multi-domain-containing protein of 5,537 amino acid residues that can methylate the Lys-4 position of histone H3 (H3K4), an epigenetic mark correlated with transcriptional active chromatin [12, 13]. MLL2 is involved in estrogen receptor α (ERα)-mediated signal transduction, acting as a coactivator of a complex that includes ASH2, RBQ3, and WDR5 .
In the present study, by direct sequencing of DNA samples from 62 Kabuki patients we identified 42 MLL2 variants, 38 of which are novel.
Subjects and Clinical Data
Clinical features of our cohort of Kabuki syndrome patients
36/62 (58.1%) Male
26/62 (41.9) Female
Long palpebral fissures
Everted lower eyelids
Large dysplastic ears
Arched eyebrows, sparse lateral one third
Flat nasal tip
Broad nasal root
Thin upper and full lower lip
Persistent fetal pads
Other clinical features (most recurrent)
leftearly breast development
leftagenesis/dysgenesis corpus callosum
Genomic DNAs were extracted from fresh and/or frozen peripheral blood leukocytes of the probands and their available family members using an automated DNA extractor and commercial DNA extraction Kits (EZ1, Qiagen, Hilden, Germany).
PCR-based sequencing of MLL2
Primers were designed using the Primer 3 Output program (http://frodo.wi.mit.edu/primer3/) to amplify the 54 coding exons of MLL2 (RefSeq NM_003482.3) gene including the intronic flanking sequences. Amplicons and primers were checked both by BLAST and BLAT against the human genome to ensure specificity. A complete list of primers is reported in Additional file 1, Table S1. The amplified products were subsequently purified and sequenced with a ready reaction kit (BigDye Terminator v1.1 Cycle, Applied Biosystems). The fragments obtained were purified using DyeEx plates (Qiagen) and resolved on an automated sequencer (3130xl Genetyc analyzer DNA Analyzer, ABI Prism). Sequences were analyzed using the Sequencer software (Gene Codes, Ann Arbor, Michigan). Whenever possible the mutations identified were confirmed on a second independent blood sample. The issue of whether the novel MLL2 missense alterations were causative mutations or neutral polymorphisms was addressed by searching dbSNP (http://www.ncbi.nlm.nih.gov/SNP) for their presence; the screening of 100 alleles from healthy unrelated control subjects and from the 1000 Genomes database  were used to assess their presence/absence in the general population. All existing and new mutations were described following the recommendations of the Human Genome Variation Society (http://www.hgvs.org/mutnomen).
In silico analysis
The putative causal and functional effect of each identified nucleotide variant was estimated by using the following in silico prediction tools: Polyphen http://genetics.bwh.harvard.edu/pph, Align GVGD http://agvgd.iarc.fr/agvgd_input.php, and MutPred http://mutdb.org/profile. Splice sites variants were evaluated for putative alteration of regulatory process at the transcriptional or splicing level with NetGene2 http://www.cbs.dtu.dk/services/NetGene2 and NNSPLICE http://www.fruitfly.org/seq_tools/splice.html. RESCUE-ESE http://genes.mit.edu/burgelab/rescue-ese/ and Fas-ESS http://genes.mit.edu/fas-ess/ online tools were used to predict exon-splicing enhancer and silencer, respectively. RepeatMasker http://www.repeatmasker.org/ was used to screen DNA sequences for the presence of direct repeats. The Coils program http://www.ch.embnet.org/software/COILS_form.html was employed to calculate the probability that the variant induces a conformational change in the coiled-coil domains. Finally multiple species alignment of MLL2 protein was made with the ClustalW software http://www.ebi.ac.uk/Tools/clustalw2/index.html using the following orthologs sequences obtained through the Ensembl genome browser http://www.ensembl.org: P. troglodytes (ENSPTRP00000041051), M. musculus (ENSMUSP00000023741), C. familiaris (gENSCAFP00000012833), B. taurus (ENSBTAP00000019193), X. troplicalis (ENSXETP00000024427), and D. rerio (ENSDARP00000053862).
Results and Discussion
Nonsense and frameshift mutations
In agreement with previous reports we identified a majority of truncating mutations (70%, 29/42), three of which were reported previously (Figure 2, Additional file 1, Table S2) [6, 7]. Most of the variants are predicted, if translated, to encode shorter MLL2 proteins either by loss of the entire C-terminal region or parts of it (Figure 2). This region harbors highly conserved domains that are found in a variety of chromatin-associated proteins [17–19]: (i) the helical LXXLL regions involved in the recruitment of the MLL2 complex to the promoters of ERα target genes (Figure 2); (ii) FYRN and FYRC sequence motifs, two poorly characterized phenylalanine/tyrosine-rich regions of around 50 and 100 amino acids , respectively; and (iii) a catalytic SET motif that confers histone methyltransferases activity.
Although it has not been yet experimentally verified for the MLL2 gene, the prevalence of premature termination mutations may result in the partial transcripts degradation through nonsense-mediated mRNA decay (NMD). NMD is an evolutionarily conserved process that typically degrades transcripts containing premature termination codons (PTCs) to prevent translation of unnecessary or aberrant and possible transcripts . The NMD process takes place when PTCs are located more than 50-55 nucleotides upstream of an exon-exon junction . As 86% (25/29) of such detected MLL2 mutations follow this rule it is likely that the consequent MLL2 haploinsufficiency could be the driving force for the onset of the Kabuki syndrome.
Our screen revealed three not yet described indel variants located in the C-terminal region of the protein (see samples KB71, KB77, and KB53 in Additional file 1, Table S2). They might have resulted from slipped mispairing between direct repeats or through the insertion or deletion of a single base within a mononucleotide tract (Additional file 1, Table S3), as already reported .
COILS algorithm predicted the amplification of one of the five coiled-coil putative domains for the c.11819_11836dupTTCAACAACAGCAGCAGC (p.Lys3940_Gln3945dup) (Figure 2 and data not shown), a domain involved in protein-protein interaction. This variant was inherited from the apparently asymptomatic mother. It is thus impossible to conclusively determine the pathogenic nature of the resulting protein.
Splice site variants
We detected 3 variants located at the splice site junctions, two of which are novel  (Additional file 1, Table S2); the in silico modeling predicted complete or partial abrogation of the junction formation with a pathogenic impact. The c.400+1G>A, occurring within the invariant GT donor splice site in intron 3-4, results in the disruption of the canonical splice site and it is expected to produce an aberrant protein of only 135 residues. The c.401-3A>G, occurring 3 bp away from the next intron-exon junction, is predicted to create a new acceptor splice site at position -3 within intron 3-4 that could lead to a frameshift encoding a mutant protein with a premature stop at 84 codons downstream. Finally, c.13999+5G>A decreases the donor site score prediction, possibly resulting in a less efficient intron splicing. Unfortunately, RNA from these patients was unavailable preventing further investigation of the effect of these variants.
Missense de novo variants have already been found in Kabuki patients. Ng and colleagues reported 8 pathogenic missense variants, two of which were recurrent in affected patients. As these were all mapping within the last exon of MLL2 that encodes the different conserved C-terminal domains of the protein (see above), the authors suggested that such mutations are tolerated, while mutations elsewhere are lethal. By in silico analysis, Paulussen et al. proposed the pathogenicity of two missense variants located within that C-terminal region [6, 7]. We detected seven patients with a single or two missense variants (KB28 and KB38 patients; Additional file 1, Table S2). PCR amplification, cloning and sequencing showed that both sets of the two sequence changes in KB28 and KB38 patients are located on one allele. From a phenotypic point of view, the two patients with pairs of missense variants do not appear to be more severely affected than affected individuals with single variants. Sequencing of the corresponding exons in the KB38 parents demonstrated that both variants arose de novo, while the KB28 patient inherited the variant from the apparently asymptomatic mother. We had also accessed to DNA of parents of carriers of missense variants (Additional file 1, Table S2). Yet, in both cases the MLL2 variant was inherited from the apparently asymptomatic father.
In silico prediction of pathogenic effect of MLL2 missense and splice site variants
Prediction of damaging effect at the protein level
Probability of deleterious mutation
Confident in silico hypothesis
loss of loop (p = 0.0288, loss of catalytic residue at R1258 (p = 0.0301); gain of helix (p = 0.0349)
gain of one ESS
loss of one ESE
gain of phoshorylation P2841 (p = 0.028)
gain of two ESEs
loss of one ESS
gain of three ESEs
loss of three ESSs
gain of ubiquitination at K5344 (p = 0.0396)
gain of two ESEs
Splice site modification prediction
Loss of DS
Loss of DS
135 AA with novel 2 AA
New AS at-3
New AS at-3
217 AA with novel 84 AA
c. 13999+5 G>A
Lower confidence of DS prediction
Lower confidence of DS prediction
Finally, we employed RESCUE-ESE and Fas-ESS tools on missense variants and frameshift mutations to predict associated splicing phenotypes by identifying sequence changes that disrupt or alter predicted Exonic Splicing Enhancers (ESE) and Exonic Splicing Silencers (ESS). ESE and ESS are short oligonucleotides that can enhance or inhibit pre-mRNA splicing when present in exons, playing important roles in constitutive and alternative splicing. A variation that disrupts an ESE, for instance, could cause exon skipping which would result in the exclusion of an entire exon from the mRNA transcript. Conversely, a substitution in the ESS sequence promotes the use of adjacent splice sites, often contributing to alternative splicing. As reported in Table 2 and Additional file 1, Table S3, we found that some of the MLL2 mutations lead to creation of new ESEs and/or to disruption of predicted wild type ESEs/ESSs. As secondary structures or adjacent negative elements also participate to the modulation of the splicing event mediated by ESE and ESS, we retain that association to functional studies will enable to better understand the role of the reported cases of ESEs/ESSs disruption or alteration in the complex phenotypic spectrum observed in the Kabuki patients.
Our study increases the number of identified MLL2 mutations and variants, and emphasizes some characteristics of the spectrum of MLL2 mutations associated with this pathology, further providing insight into its etiology. The in silico analysis predicts that the identified MLL2 missense, splice-site and indel variants might be pathogenic. Other studies reported the presence of such MLL2 variants predicted to be associated with the disease. However, their biological significance and pathogenicity were not unambiguously demonstrated; therefore further and more functionally oriented studies are needed to understand the nature of these variants and their possible role in the disease. Solving these issues is relevant to avoid any incorrect interpretation and diagnosis of other Kabuki cases carrying such MLL2 variants.
We were unable to found any detectable point mutation and/or small del/dup in the coding region of MLL2 gene in 27% (17/62) of the Kabuki syndrome patients. Mutations in MLL2 regulatory regions, exon microduplications and/or microdeletions, as well as genetic heterogeneity of the syndrome may account for these negative results. Alternatively, some of these patients might have been misdiagnosed as a result of the complex clinical spectrum covered by this pathology, thus possibly highlighting the need to more accurately select Kabuki cases before conducting the analysis.
In summary, this study underlines the relevance of mutational screening of the MLL2 gene among patients with Kabuki syndrome. The identification of a large spectrum of MLL2 mutations will offer the opportunity to improve the actual knowledge on the clinical basis of this multiple congenital mental retardation syndrome, to design functional studies to understand the molecular mechanisms underlying the disease, to establish genotype-phenotype correlations, to improve the clinical management, and to identify potential targets for therapy.
Acknowledgements and funding
We thank all the families that agree to participate and get possible this study. This work was in part supported by grants from the Italian Ministry of Health (Ricerca Corrente 2008-10), the Fondazione Banca del Monte di Foggia "Domenico Siniscalco Ceci", the Jerome Lejeune Foundation, and with the contribution of Ministero degli Affari Esteri, Direzione Generale per la Promozione e la Cooperazione Culturale (2009-2010) to GM. We thank Mariani Foundation, Milan for the support to the activity of Ambulatorio di Genetica Clinica Pediatrica Fondazione MBBM to AS. AR is supported by grants from the Swiss National Science Foundation. We thank the Galliera Genetic Bank - Network of Telethon Genetic Biobanks (project GTB07001) for providing us with some Kabuki samples.
The founders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Written consent for publication was obtained from the patients or their relatives.
- Kuroki Y, Suzuki Y, Chyo H, Hata A, Matsui I: 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. 1981, 99 (4): 570-573. 10.1016/S0022-3476(81)80256-9.View ArticlePubMedGoogle Scholar
- Niikawa N, Matsuura N, Fukushima Y, Ohsawa T, Kajii T: Kabuki make-up syndrome: a syndrome of mental retardation, unusual facies, large and protruding ears, and postnatal growth deficiency. J Pediatr. 1981, 99 (4): 565-569. 10.1016/S0022-3476(81)80255-7.View ArticlePubMedGoogle Scholar
- Niikawa N, Kuroki Y, Kajii T, Matsuura N, Ishikiriyama S, Tonoki H, Ishikawa N, Yamada Y, Fujita M, Umemoto H: Kabuki make-up (Niikawa-Kuroki) syndrome: a study of 62 patients. Am J Med Genet. 1988, 31 (3): 565-589. 10.1002/ajmg.1320310312.View ArticlePubMedGoogle Scholar
- Cusco I, del Campo M, Vilardell M, Gonzalez E, Gener B, Galan E, Toledo L, Perez-Jurado LA: Array-CGH in patients with Kabuki-like phenotype: identification of two patients with complex rearrangements including 2q37 deletions and no other recurrent aberration. BMC Med Genet. 2008, 9: 27.PubMed CentralView ArticlePubMedGoogle Scholar
- Milunsky JM, Huang XL: Unmasking Kabuki syndrome: chromosome 8p22-8p23.1 duplication revealed by comparative genomic hybridization and BAC-FISH. Clin Genet. 2003, 64 (6): 509-516. 10.1046/j.1399-0004.2003.00189.x.View ArticlePubMedGoogle Scholar
- Ng SB, Bigham AW, Buckingham KJ, Hannibal MC, McMillin MJ, Gildersleeve HI, Beck AE, Tabor HK, Cooper GM, Mefford HC, Lee C, Turner EH, Smith JD, Rieder MJ, Yoshiura K, Matsumoto N, Ohta T, Niikawa N, Nickerson DA, Bamshad MJ, Shendure J: Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome. Nat Genet. 2010, 42 (9): 790-793. 10.1038/ng.646.PubMed CentralView ArticlePubMedGoogle Scholar
- Paulussen AD, Stegmann AP, Blok MJ, Tserpelis D, Posma-Velter C, Detisch Y, Smeets EE, Wagemans A, Schrander JJ, van den Boogaard MJ, van der Smagt J, van Haeringen A, Stolte-Dijkstra I, Kerstjens-Frederikse WS, Mancini GM, Wessels MW, Hennekam RC, Vreeburg M, Geraedts J, de Ravel T, Fryns JP, Smeets HJ, Devriendt K, Schrander-Stumpel CT: MLL2 mutation spectrum in 45 patients with Kabuki syndrome. Hum Mutat. 2010.Google Scholar
- Dillon SC, Zhang X, Trievel RC, Cheng X: The SET-domain protein superfamily: protein lysine methyltransferases. Genome Biol. 2005, 6 (8): 227. 10.1186/gb-2005-6-8-227.PubMed CentralView ArticlePubMedGoogle Scholar
- Issaeva I, Zonis Y, Rozovskaia T, Orlovsky K, Croce CM, Nakamura T, Mazo A, Eisenbach L, Canaani E: Knockdown of ALR (MLL2) reveals ALR target genes and leads to alterations in cell adhesion and growth. Mol Cell Biol. 2007, 27 (5): 1889-1903. 10.1128/MCB.01506-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Ansari KI, Mandal SS: Mixed lineage leukemia: roles in gene expression, hormone signaling and mRNA processing. FEBS J. 2010, 277 (8): 1790-1804. 10.1111/j.1742-4658.2010.07606.x.View ArticlePubMedGoogle Scholar
- Eissenberg JC, Shilatifard A: Histone H3 lysine 4 (H3K4) methylation in development and differentiation. Dev Biol. 2010, 339 (2): 240-249. 10.1016/j.ydbio.2009.08.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Malik S, Bhaumik SR: Mixed lineage leukemia: histone H3 lysine 4 methyltransferases from yeast to human. FEBS J. 2010, 277 (8): 1805-1821. 10.1111/j.1742-4658.2010.07607.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Strahl BD, Ohba R, Cook RG, Allis CD: Methylation of histone H3 at lysine 4 is highly conserved and correlates with transcriptionally active nuclei in Tetrahymena. Proc Natl Acad Sci USA. 1999, 96 (26): 14967-14972. 10.1073/pnas.96.26.14967.PubMed CentralView ArticlePubMedGoogle Scholar
- Mo R, Rao SM, Zhu YJ: Identification of the MLL2 complex as a coactivator for estrogen receptor alpha. J Biol Chem. 2006, 281 (23): 15714-15720. 10.1074/jbc.M513245200.View ArticlePubMedGoogle Scholar
- Kawame H, Hannibal MC, Hudgins L, Pagon RA: Phenotypic spectrum and management issues in Kabuki syndrome. J Pediatr. 1999, 134 (4): 480-485. 10.1016/S0022-3476(99)70207-6.View ArticlePubMedGoogle Scholar
- Durbin RM, Abecasis GR, Altshuler DL, Auton A, Brooks LD, Gibbs RA, Hurles ME, McVean GA: A map of human genome variation from population-scale sequencing. Nature. 2010, 467 (7319): 1061-1073. 10.1038/nature09534.View ArticlePubMedGoogle Scholar
- Balciunas D, Ronne H: Evidence of domain swapping within the jumonji family of transcription factors. Trends Biochem Sci. 2000, 25 (6): 274-276. 10.1016/S0968-0004(00)01593-0.View ArticlePubMedGoogle Scholar
- Doerks T, Copley RR, Schultz J, Ponting CP, Bork P: Systematic identification of novel protein domain families associated with nuclear functions. Genome Res. 2002, 12 (1): 47-56. 10.1101/gr.203201.PubMed CentralView ArticlePubMedGoogle Scholar
- Prasad R, Zhadanov AB, Sedkov Y, Bullrich F, Druck T, Rallapalli R, Yano T, Alder H, Croce CM, Huebner K, Mazo A, Canaani E: Structure and expression pattern of human ALR, a novel gene with strong homology to ALL-1 involved in acute leukemia and to Drosophila trithorax. Oncogene. 1997, 15 (5): 549-560. 10.1038/sj.onc.1201211.View ArticlePubMedGoogle Scholar
- Garcia-Alai MM, Allen MD, Joerger AC, Bycroft M: The structure of the FYR domain of transforming growth factor beta regulator 1. Protein Sci. 2010, 19 (7): 1432-1438. 10.1002/pro.404.PubMed CentralView ArticlePubMedGoogle Scholar
- Maquat LE: When cells stop making sense: effects of nonsense codons on RNA metabolism in vertebrate cells. RNA. 1995, 1 (5): 453-465.PubMed CentralPubMedGoogle Scholar
- Nagy E, Maquat LE: A rule for termination-codon position within intron-containing genes: when nonsense affects RNA abundance. Trends Biochem Sci. 1998, 23 (6): 198-199. 10.1016/S0968-0004(98)01208-0.View ArticlePubMedGoogle Scholar
- Tappino B, Biancheri R, Mort M, Regis S, Corsolini F, Rossi A, Stroppiano M, Lualdi S, Fiumara A, Bembi B, Di Rocco M, Cooper DN, Filocamo M: Identification and characterization of 15 novel GALC gene mutations causing Krabbe disease. Hum Mutat. 2010, 31 (12): E1894-1914. 10.1002/humu.21367.PubMed CentralView ArticlePubMedGoogle Scholar
- Ensembl Genome Browser. [http://www.ensembl.org].
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.