Skip to main content

A novel CHD7 variant disrupting acceptor splice site in a patient with mild features of CHARGE syndrome: a case report



CHARGE syndrome (MIM# 214800)—which is characterised by a number of congenital anomalies including coloboma, ear anomalies, deafness, facial anomalies, heart defects, atresia choanae, genital hypoplasia, growth retardation, and developmental delay—is caused by a heterozygous variant in the CHD7 (MIM# 608892) gene located on chromosome 8q12. We report the identification of a novel c.5535-1G > A variant in CHD7 and provide the evaluation of its effect on pre-mRNA splicing.

Case presentation

In this study, we report on a female presenting features of CHARGE syndrome. A novel heterozygous CHD7 variant c.5535-1G > A located in the acceptor splice site of intron 26 was identified in the proband’s DNA sample after analysis of whole exome sequencing data. In silico predictions indicating that the variant is probably pathogenic by affecting pre-mRNA splicing were verified by genetic analysis based on reverse transcription of the patient’s RNA followed by PCR amplifications performed on synthesised cDNA and Sanger sequencing. Sanger sequencing of cDNA revealed that the c.5535-1G > A variant disrupts the original acceptor splice site and activates a cryptic splice site only one nucleotide downstream of the pathogenic variant site. This change causes the omission of the first nucleotide of exon 27, leading to a frameshift in the mRNA of the CHD7 gene. Our results suggest that the alteration induces the premature truncation of the CHD7 protein (UniProtKB: Q9P2D1), thus resulting in CHARGE syndrome.


Genetic analysis of novel splice site variant underlines its importance for studying the pathogenic splicing mechanism as well as for confirming a diagnosis.

Peer Review reports


Rapid advancement of high-throughput next generation sequencing (NGS) technologies has provided the necessary conditions for the identification of disease-causing variants in the human genome [1]. Even though the technique for detecting variants is now becoming routine, the challenge of performing a detailed analysis and determining the function of the altered gene is still open [2].

According to data from the 1000 Genomes Project, a typical human genome contains 149–182 sites with protein truncating variants, more than 10,000 sites with peptide sequence altering variants, and about 459,000–565,000 variant sites affecting known regulatory regions [3]. Previous studies have shown that many genetic diseases have been caused by coding and non-coding variants disrupting splice sites [4, 5]. Variants altering pre-mRNA splicing account for at least 10% of disease-causing gene alterations [6].

The CHD7 (Chromodomain Helicase DNA Binding Protein 7; MIM #608892; NM_017780.4) gene, which is located on chromosome 8q12, encodes protein that contains several helicase family domains [7]. Heterozygous variants within this gene are the major cause of CHARGE syndrome (MIM #214800; ORPHA: 138). This syndrome is a rare autosomal dominant disorder estimated to occur in about 1:10,000 births worldwide [8]. CHARGE syndrome is characterised by a number of different congenital anomalies, including Coloboma, Heart defects, Atresia of the choanae, Retardation of growth and development, Genital hypoplasia, Ear anomalies, and deafness [9, 10]. CHARGE syndrome shares several major findings and phenotypic overlap with Mandibulofacial dysostosis, Guion-Almeida type (MFDGA; MIM #610536; ORPHA: 79113) [11, 12] and Kabuki syndrome 1 (KABUK1; MIM #147920; ORPHA: 2322) [13], which are caused by heterozygous variants in EFTUD2 (MIM #603892) and KMT2D (MIM #602113), respectively. There is some evidence that the phenotype can be caused by a variant in the SEMA3E (Semaphorin-3E; MIM #608166) gene [14], but almost 60–70% of patients have a pathogenic variant in the CHD7 gene [15, 16]. In the CHD7 database [17, 18], 801 pathogenic variants in CHD7 are currently recorded. According to the Human Gene Mutation Database (HGMD) [19, 20], the most frequent variants in the CHD7 gene are nonsense, small deletions, and insertions. Moreover, a number of different pathogenic variants in CHD7 splice sites have also been recorded.

This study aims to reveal the consequences and the pathogenicity of a novel heterozygous splice site c.5535-1G > A variant in the CHD7 gene in a female patient with mild features of CHARGE syndrome.

Case presentation

Clinical findings

A female patient was referred for genetic assessment at 16 years of age for hearing impairment, scoliosis, and left side facial palsy. The proband was the second child of healthy unrelated Lithuanian parents aged 29 at the proband’s birth. The proband had neither a history of prenatal exposure to teratogenic agents nor any family history of congenital malformations. During the pregnancy, polyhydramnios was diagnosed. Congenital heart defect of the foetus was detected at 35 weeks of gestation. The proband was born at 38 weeks by Caesarian section due to congenital heart defect of the foetus. At birth, her weight was 2800 g (3rd centile), her length was 49 cm (10th centile), the circumference of her head was 35 cm (50th centile), and her Apgar score was 6 at 1 min and 8 at 5 min. Complete atrioventricular communication with a small ventricular septal defect, a moderately sized atrial septal defect with an aneurysm in the interatrial septum, and insufficiency of the mitral valve were diagnosed. Additionally, left side choanal atresia, left side facial palsy, a short neck, dysplasia of the sacrum, and dysplasia of the hip joint were observed. As a neonate, the patient had feeding difficulties and received a nasogastric tube for 1 month. Poor weight gain and cardiac insufficiency manifested at the age of 7 months. Auditory Brainstem Response Audiometry was performed and unilateral deafness in the left was diagnosed during infancy. Bilateral aplasia of the semicircular canals and cochlear nerve canal atresia on the left were found on temporal bone CT scans (Fig. 1). At the age of 1 year, radical surgical correction of atrioventricular communication was performed. Pyelonephritis and grade IV vesicoureteral reflux were first diagnosed at 1 year of age, and recurrent episodes of pyelonephritis continued until the age of 7 years. At age of 4 years, surgical repair of choanal atresia on the left was performed and a tympanostomy tube was placed in the same side due to recurrent otitis media. The choanal atresia operation was performed again at the age of 5–6 years. Bone-anchored hearing aid (BAHA) surgery was performed at 7 years of age. Due to a skin reaction at the location of the implant, another operation to implant a new BAHA attract system was performed after the device had been used for 4 years. Adenoidectomy was performed at 11 years of age. The proband had hypermetropia, for which glasses were prescribed. Scoliosis of 17° was evaluated at the age of 16 years.

Fig. 1

Temporal bone CT image showed aplasia of the semi-circular canals (a) and cochlear nerve canal atresia (b)

Motor development of the proband was delayed; she walked alone at the age of 27 months. According to the DISC scale, at the age of 29 months her psychomotor development was evaluated at the range of 21 (gross motor) to 29 (fine motor, visual attention and memory) months. Eruption of permanent teeth was also mildly delayed. The girl attended a regular school and had no behavioural problems. Her psychological evaluation at the age of 16 years yielded a verbal score of 90, a performance score of 97, and a full score of 94 on the Wechsler Intelligence Scale for Children (WISC-III).

During an examination at the age of 16 years, her height was 162 cm (25th centile), her weight was 48.5 kg (3rd–10th centile), and the circumference of her head was 53.5 cm (3rd–10th centile). Her phenotype was remarkable for left side facial palsy, low forehead, deformed nasal bridge, synophrys, protruding ears, auricular pits, low posterior hairline, short neck, limited movements of the neck, pectus excavatum, hirsutism, and Tanner PH-4, B-3. Additionally, primary amenorrhoea and premature ovarian insufficiency were diagnosed at 18 years of age. She had low levels of estradiol and near normal concentrations of prolactin and follicle-stimulating hormone. A pelvic ultrasound visualised a uterus (49x15x31 mm) and ovaries (26-27 × 12 mm).

Whole exome sequencing

Whole genomic DNA (gDNA) was extracted from peripheral blood of our proband and her healthy parents following the standard phenol-chloroform extraction protocol.

Whole exome sequencing (WES) using the high throughput next generation sequencing (NGS) technique was used to sequence the sample of the proband (Illumina, Inc., USA). DNA libraries generated using TruSeq Rapid Exome Library Prep kit (8x3plex) (Illumina, Inc., USA). In order to precisely measure the concentration of DNA libraries, Qubit dsDNA BR Assay kit (ThermoFisher Scientific, USA) and Qubit fluorimeter (ThermoFisher Scientific, USA) were used. Clusters amplificated using cBot system (Illumina, Inc., USA), TruSeq PE Cluster Kit v3-HS (Illumina, Inc., USA) and TruSeq Dual Index Sequencing Primer Box – Paired End (Illumina, Inc., USA). WES was performed using TruSeq SBS Kit v3-HS (Illumina, Inc., USA) and by employing the HiScanSQ (Illumina, Inc., USA) genetic analyzer.

Analyses of high throughput sequencing data started with alignment against The Human NCBI Build GRCh37 (hg19/2009) reference genome. The annotation of NGS data was made using the ANNOVAR v.2018Apr16 program. The pathogenicity of variants was assessed using ACMG criteria, taking into account the data provided by the ANNOVAR program, the available databases (ExAC Browser, 1000 Genome Project, NCBI dbSNP, NCBI dbVar, HGMD, NCBI OMIM, NCBI ClinVAR, Leiden Open Variation Database, NCBI Genome, Deafness Variation Database), and the relevant scientific literature. The pathogenic or probably pathogenic sequence variants were checked by analyzing proband’s BAM files using the visualization tool Integrative Genomics Viewer (IGV).

Reverse transcription reaction

Total RNA of the proband and her healthy parents’ samples was extracted from whole blood using a Tempus™ Blood RNA- Tube and Tempus™ Spin RNA Isolation Kit (Thermo Fisher Scientific, USA) according to the optimized manufacturer’s protocols. Total RNA concentration and quality were checked using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). Complementary DNA (cDNA) was reverse-transcribed (RT) from a total RNA using a High-Capacity RNA-to-cDNA Kit (Thermo Fisher Scientific, USA) and a ProFlex PCR system (Thermo Fisher Scientific, USA) following manufacturer’s protocol and recommended conditions.

Polymerase chain reactions and Sanger sequencing

Polymerase chain reactions (PCR) of gDNA and cDNA sequences flanking splice site c.5535-1G > A variant of CHD7 were performed using specific primers designed with Primer Blast tool [21, 22]. gDNA was PCR amplified using primer pair designed on exon 26 (forward primer) and on intron 27 (reverse primer), while PCR amplifications of cDNA sequence were performed using the same forward primer, but different reverse primers, which were designed on the 28–29 exon junction as well as on the 29–30 exon junction (Additional file 1: Table S1). PCR products were fractioned according to standard agarose gel electrophoresis.

gDNA samples was used for the segregation analysis. In order to elucidate the pathogenicity of detected variant, the cDNA samples were analysed. The PCR products were sequenced with BigDye® Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific, USA) and ABI 3130xL Genetic Analyser (Thermo Fisher Scientific, USA). The sequences were aligned with the reference sequence of the CHD7 (NCBI: NM_017780.4) gene.

In silico analysis

Mutation Taster [23, 24] and Human Splicing Finder [25] databases were used for predicting splice site alterations. Sequences of evolutionary distinct species were obtained from the Ensembl genome browser [26, 27], while a sequence alignment was produced using Clustal Omega tool [28]. Possible splice site c.5535-1G > A variant’s effect on the CHD7 protein (UniProtKB: Q9P2D1) was predicted using different tools and databases, such as ExPASy Bioinformatics Resource Portal [29, 30] and Pfam 32.0 database [31, 32].

Genetic findings

Heterozygous CHD7 variant NC_000008.11(NM_017780.4):c.5535-1G > A, located in the acceptor splice site of intron 26, was identified in the proband’s gDNA sample after analysis of WES data.

The acceptor splice site variant is novel and not recorded in the CHD7 database, HGMD, or other databases. The variant was confirmed by Sanger sequencing on the proband’s gDNA sample. Segregation analysis by Sanger sequencing on the gDNA samples of her parents confirmed the de novo origin of the variant. The splice site variant was evaluated by in silico analysis.

A variant effect prediction method, Mutation Taster, predicted the variant to be likely pathogenic and disease causing. Moreover, the alteration of the acceptor splice site was predicted by the Human Splicing Finder database. The variant was predicted to affect pre-mRNA splicing by the altering acceptor splice site at − 1 position and the possible use of a cryptic splice site. Sequence alignment of the CHD7 protein in seven species revealed that the region containing the splice site variant is highly conserved (Fig. 2b).

Fig. 2

a The chromatogram and schematic representation of wild-type and mutated sequences illustrate the disrupted acceptor splice site at position − 1 of intron 26, which lead to a frameshift in the cDNA sample. b A comparative sequence alignment produced by ClustalO of the CHD7 protein across seven evolutionarily distant species. Frameshift of 23 new amino acids and therefore a truncated 1869 amino acid sequence is highlighted in red. Below the protein sequences is a key denoting conserved sequence (*), conservative variants (:), semi-conservative variants (.), and non-conservative variants (). c 3D model of the BRK domain, which is absent in the patient. Picture produced by the LiteMol program [33]. d A schematic representation of the CHD7 protein, which is 2,997 amino acids in length, with the arrangement of the main domains: two chromodomains (green), SNF2 domain (red), helicase domain (blue), and two BRK domains (yellow). A schematic view was adopted from the Pfam protein domain database

To study the potential pathogenicity of the CHD7 gene c.5535-1G > A variant, the analysis of the cDNA sequence of the patient and healthy controls was performed. Sanger sequencing of the patient’s cDNA sample revealed that the c.5535-1G > A variant causes the loss of an original acceptor splice site at position − 1 in the intron 26 of CHD7 and consequently activates a cryptic splice site only one nucleotide downstream of the pathogenic variant site (Fig. 2a). Meanwhile, Sanger sequencing of cDNA samples from healthy controls did not expose any alteration in the splicing mechanism. In silico, this variant leads to the frameshift of 23 new amino acids and therefore a truncated 1869 amino acid protein NP_060250.2:p.(Gly1846Glufs*23) (Fig. 2d). The truncated CHD7 protein contains two functionally important chromodomains, the SNF2 domain and helicase domain, but there is a loss of two BRK domains (Fig. 2c).

Discussion and conclusions

Variants that affect pre-mRNA splicing and eventually translated proteins have been shown to be one of the most frequent causes of hereditary disorders [4]. The characterisation of splice site variants of unknown consequence and clinical significance is critical for appropriate and comprehensive genetic testing [34]. The most common consequence of splicing variants is skipping of one or more exons, but in certain cases the activation of cryptic 5′ (donor) or 3′ (acceptor) splice site and partial or full retention of specific introns may occur [35]. However, the effects of these variants on splicing are often not pursued as a cause of disease. Since there are different consequences of these variants, they are especially difficult to predict without mRNA or cDNA analysis [5].

In this study, we presented the characterisation of a previously unreported de novo heterozygous acceptor splice site c.5535-1G > A variant in the CHD7 gene. Earlier studies of variants in the 3′ splice site have shown that transitions, as well as in this case guanine (G) substitution for adenine (A), occur much more frequently than transversions, and the most frequently mutated base of an aberrant 3′ splice site is G [36]. The patient carrying this acceptor splice site variant fulfils major diagnostic characteristics of CHARGE syndrome outlined by Blanke et al. [37] and updated by Verloes [38]: the proband has three major characteristics (unilateral choanal atresia, facial palsy, temporal bone abnormalities) and three minor characteristics (delayed puberty, delayed mylestones, cardiovascular malformation). The haploinsufficiency of the CHD7 gene is assumed to be the pathogenic mechanism of CHARGE syndrome [10]. This mechanism is supported by studies in which a heterozygous knock-out mouse model was engineered. In these studies, heterozygous Chd7 mice presented features similar to those associated with CHARGE syndrome, providing strong support to the importance of the hemizygosity of CHD7 in this syndrome [39,40,41]. Heterozygous variants of the CHD7 gene play important roles in a number of vitally important processes such as embryonic development, cell cycle arrest, apoptosis, transcriptional regulation, DNA repair, replication, and recombination. Furthermore, the CHD7 protein is essential for regulation of distal genomic elements that control the expression of critical neural crest transcriptional factors. For this reason, it is hypothesised that the loss of 50% of the CHD7 protein is insufficient to maintain normal function of this gene, thus resulting in clinical manifestation [10, 16, 42].

Our further investigation of the acceptor splice site c.5535-1G > A variant detected by WES aimed to elucidate the pathogenicity of the novel variant in the CHD7 gene. In silico analysis predicted this variant to be likely pathogenic. Moreover, the sequence alignment of the CHD7 protein across seven evolutionarily distant species revealed that the region containing the variant at the acceptor splice site is highly conserved, thus suggesting that this region is fundamentally important in vertebrates and alterations could be pathogenic (Fig. 2b). Our following research was based on cDNA analysis of the CHD7 gene transcript. Unexpectedly, the data obtained on the patient’s cDNA showed that the c.5535-1G > A variant does not lead to adjacent exon skipping or intron retention; instead it activates a cryptic 3′ splice site only one nucleotide downstream from the pathogenic variant site. Eng et al. [43] presented a similar study of the acceptor splice site c.4437-1G > A variant in the ATM gene, which caused the same consequences to the cDNA sequence. In both cases, the first nucleotide of a particular exon is a G. The cryptic 3′ splice site has therefore been created using the mutated A of the intron and the G of the exon. These findings could suggest that if the original 3′ splice site is mutated the subunits of U2 spliceosomal snRNAs recognises the next AG site, which is located downstream of the original acceptor splice site. However, Moon et al. [44] reported two abnormal transcripts of the HRPT2 gene, also known as the CDC73 gene, resulting from the same c.238-1G > A variant of intron 2. The first transcript eliminated the whole exon 3, while the another lacked the first 23 nucleotides of exon 3 due to the use of a cryptic acceptor splice site in exon 3. For this reason, the exact mechanism for how one cryptic splice site is selected over another within the exon remains unclear. According to literature that we reviewed and the results of our experiment, the variant at − 1 position in the acceptor splice site and the first G nucleotide in the adjacent exon sequence is quite frequent [45, 46]. However, such variants usually disturb the reading frame by exon skipping [35] instead of the slippage of the splicing mechanism only 1 nt downstream as in our case. Moreover, the possible effect of altered splice sites is usually not based on mRNA or cDNA analysis and is predicted by computational algorithms [45, 47, 48] though the exact effect of the specific variant should be verified in functional studies [49].

Further in silico analysis showed that the activation of the cryptic 3′ splice site, including the first G nucleotide of exon 27, consequently disrupted the reading frame of the patient’s cDNA sequence, resulting in a truncated 1,869 amino acid protein NP_060250.2:p.(Gly1846Glufs*23). The original CHD7 protein (UniProtKB: Q9P2D1), which is 2,997 amino acids in length, has two N-terminal chromodomains and an SNF2 family domain that are involved in a variety of cellular processes but mainly affect the regulation of transcription. It also contains a C-terminal conserved helicase domain as well as two C-terminal BRK domains. The function of the BRK domain is still unknown (Fig. 2), but it is hypothesised that this domain may interact with chromatin components that are important and unique to higher eukaryotes [50, 51]. In our study, the truncated CHD7 protein is predicted to result in the loss of two BRK domains.

Based on these findings, our hypothesis is that the heterozygous acceptor splice site c.5535-1G > A variant in the CHD7 gene results in partial haploinsufficiency of the CHD7 gene in a patient with suspected CHARGE syndrome, because the variant causes premature termination of the translation of the mutated allele. The partial haploinsuficiency of the CHD7 gene could also clarify the mild features of CHARGE syndrome that are characteristic of our patient. The results reported here explain the pathogenicity of this novel variant and predict the possible effect on the CHD7 protein, but further proteomic experiments are needed to clarify and validate the molecular effects caused by the c.5535-1G > A splice variant on the CHD7 protein.

In conclusion, despite a growing number of reported splicing variants and associated phenotypes, the effect on the genotype and the phenotype persists difficult to demonstrate, and therefore it remains one of the principal challenges in human molecular genetics. Our study showed that the novel splice acceptor variant c.5535-1G > A has a complete impact on the splicing of the CHD7 gene, and the effect is predicted to be a truncated protein, thus resulting in CHARGE syndrome. These results enhance the importance of studying the pathogenicity of the splicing variant at the cDNA level and provide an opportunity to validate the suspected diagnosis and tailor clinical management for this family. Finally, the cDNA analysis provides a unique possibility to understand the etiology and pathophysiology and identify the molecular basis of many hereditary diseases and conditions, including those associated with congenital anomalies.

Availability of data and materials

The main data generated and analysed during this study are included in this article and its supplementary information file. Any additional information is available from the authors upon request.





Gene symbol for gene that encodes the protein ATM Serine/Threonine kinase


Bone-anchored hearing aid

CDC73 :

Gene symbol for gene that encodes the protein cell division cycle 73


Complementary DNA

CHD7 :

Gene symbol for gene that encodes the protein chromodomain helicase DNA binding


Computed tomography


Gene symbol for gene that encodes the protein elongation factor Tu GTP-binding domain-containing 2




Genomic DNA


Human gene mutation database


Gene symbol for gene that encodes the protein hyperparathyroidism 2


Integrative genomics viewer


Kabuki syndrome 1


Keratitis-ichthyosis-deafness syndrome


Gene symbol for gene that encodes the protein lysine-specific methyltransferase 2d


Mandibulofacial dysostosis, Guion-Almeida type


Next generation sequencing


Polymerase chain reaction


Reference sequence database in NCBI


reverse transcription


Gene symbol for gene that encodes the protein semaphorin-3E


Whole exome sequencing


Wechsler intelligence scale for children


  1. 1.

    Kircher M, Kelso J. High-throughput DNA sequencing - concepts and limitations. BioEssays. 2010;32:524–36.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Cooper DN, Krawczak M, Polychronakos C, Tyler-Smith C, Kehrer-Sawatzki H. Where genotype is not predictive of phenotype: towards an understanding of the molecular basis of reduced penetrance in human inherited disease. Hum Genet. 2013;132:1077–130.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al. A global reference for human genetic variation. Nature. 2015;526:68–74.

    CAS  Article  Google Scholar 

  4. 4.

    López-Bigas N, Audit B, Ouzounis C, Parra G, Guigó R. Are splicing mutations the most frequent cause of hereditary disease? FEBS Lett. 2005;579:1900–3.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Wang G-S, Cooper TA. Splicing in disease: disruption of the splicing code and the decoding machinery. Nat Rev Genet. 2007;8:749–61.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Krawczak M, Thomas NST, Hundrieser B, Mort M, Wittig M, Hampe J, et al. Single base-pair substitutions in exon-intron junctions of human genes: nature, distribution, and consequences for mRNA splicing. Hum Mutat. 2007;28:150–8.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Nagase T, Kikuno R, Ishikawa KI, Hirosawa M, Ohara O, Nakayama M, et al. Prediction of the coding sequences of unidentified human genes. XVI. The complete sequences of 150 new cDNA clones from brain which code for large proteins in vitro. DNA Res. 2000;7:65–73.

    CAS  Article  Google Scholar 

  8. 8.

    Issekutz KA, Graham JM, Prasad C, Smith IM, Blake KD. An epidemiological analysis of CHARGE syndrome: preliminary results from a Canadian study. Am J Med Genet. 2005;133(A):309–17.

    Article  Google Scholar 

  9. 9.

    Zentner GE, Hurd EA, Schnetz MP, Handoko L, Wang C, Wang Z, et al. CHD7 functions in the nucleolus as a positive regulator of ribosomal RNA biogenesis. Hum Mol Genet. 2010;19:3491–501.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Vissers LELM, van Ravenswaaij CMA, Admiraal R, Hurst JA, de Vries BBA, Janssen IM, et al. Mutations in a new member of the chromodomain gene family cause CHARGE syndrome. Nat Genet. 2004;36:955–7.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Gordon CT, Petit F, Oufadem M, Decaestecker C, Jourdain A-S, Andrieux J, et al. EFTUD2 haploinsufficiency leads to syndromic oesophageal atresia. J Med Genet. 2012;49:737–46.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Luquetti DV, Hing AV, Rieder MJ, Nickerson DA, Turner EH, Smith J, et al. “Mandibulofacial dysostosis with microcephaly” caused by EFTUD2 mutations: expanding the phenotype. Am J Med Genet Part A. 2013;161:108–13.

    CAS  Article  Google Scholar 

  13. 13.

    Sakata S, Okada S, Aoyama K, Hara K, Tani C, Kagawa R, et al. Individual clinically diagnosed with CHARGE syndrome but with a mutation in KMT2D, a gene associated with kabuki syndrome: a case report. Front Genet. 2017;8:1–6.

    CAS  Article  Google Scholar 

  14. 14.

    Lalani SR. SEMA3E mutation in a patient with CHARGE syndrome. J Med Genet. 2004;41:e94–4.

    CAS  Article  Google Scholar 

  15. 15.

    Martin DM, Sheldon S, Gorski JL. CHARGE association With choanal atresia and inner ear hypoplasia in a child with a de novo chromosome translocation t(2;7)(p14;q21.11). Am J Med Genet. 2001;99:115–9.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Zentner GE, Layman WS, Martin DM, Scacheri PC. Molecular and phenotypic aspects of CHD7 mutation in CHARGE syndrome. Am J Med Genet Part A. 2010;152A:674–86.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Janssen N, Bergman JEH, Swertz MA, Tranebjaerg L, Lodahl M, Schoots J, et al. Mutation update on the CHD7 gene involved in CHARGE syndrome. Hum Mutat. 2012;33:1149–60.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    CHD7 database. Accessed 3 Feb 2019.

  19. 19.

    Stenson PD, Mort M, Ball EV, Evans K, Hayden M, Heywood S, et al. The human gene mutation database: towards a comprehensive repository of inherited mutation data for medical research, genetic diagnosis and next-generation sequencing studies. Hum Genet. 2017;136:665–77.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    The Human Gene Mutation Database. Cardiff University. 2007. Accessed 3 Feb 2019.

  21. 21.

    Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:134.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    National Center for Biotechnology Information. U.S. National Library of Medicine. Accessed 25 Mar 2018.

  23. 23.

    Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11:361–2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    MutationTaster. Accessed 30 Apr 2018.

  25. 25.

    Human Splicing Finder. Aix Marseille Université. 2013. Accessed 30 Apr 2018.

  26. 26.

    Hubbard T, Barker D, Birney E, Cameron G, Chen Y, Clark L, et al. The Ensembl genome database project. Nucleic Acids Res. 2002;30:38–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Ensembl. Accessed 2 May 2018.

  28. 28.

    EMBL-EBI. Wellcome Genome Campus, Hinxton, Cambridgeshire, UK. Accessed 2 May 2018.

  29. 29.

    ExPASy Bioformatics Resources Portal. Accessed 2 May 2018.

  30. 30.

    Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40:W597–603.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47:D427–32.

    Article  PubMed  Google Scholar 

  32. 32.

    Pfam. 2018. Accessed 2 May 2018.

  33. 33.

    Sehnal D, Deshpande M, Vařeková RS, Mir S, Berka K, Midlik A, et al. LiteMol suite: interactive web-based visualization of large-scale macromolecular structure data. Nat Methods. 2017;14:1121–2.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Rodenburg RJ. The functional genomics laboratory: functional validation of genetic variants. J Inherit Metab Dis. 2018;41:297–307.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Krawczak M, Reiss J, Cooper D. The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum Genet. 1992;90:41–54.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Vořechovský I. Aberrant 3′ splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization. Nucleic Acids Res. 2006;34:4630–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Blake KD, Davenport SLH, Hall BD, Hefner MA, Pagon RA, Williams MS, et al. CHARGE association: an update and review for the primary pediatrician. Clin Pediatr (Phila). 1998;37:159–73.

    CAS  Article  Google Scholar 

  38. 38.

    Verloes A. Updated diagnostic criteria for CHARGE syndrome: a proposal. Am J Med Genet Part A. 2005;133A:306–8.

    Article  PubMed  Google Scholar 

  39. 39.

    Bosman EA, Penn AC, Ambrose JC, Kettleborough R, Stemple DL, Steel KP. Multiple mutations in mouse Chd7 provide models for CHARGE syndrome. Hum Mol Genet. 2005;14:3463–76.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Hurd EA, Capers PL, Blauwkamp MN, Adams ME, Raphael Y, Poucher HK, et al. Loss of Chd7 function in gene-trapped reporter mice is embryonic lethal and associated with severe defects in multiple developing tissues. Mamm Genome. 2007;18:94–104.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Ogier JM, Arhatari BD, Carpinelli MR, McColl BK, Wilson MA, Burt RA. An intronic mutation in Chd7 creates a cryptic splice site, causing aberrant splicing in a mouse model of CHARGE syndrome. Sci Rep. 2018;8:5482.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Bajpai R, Chen DA, Rada-Iglesias A, Zhang J, Xiong Y, Helms J, et al. CHD7 cooperates with PBAF to control multipotent neural crest formation. Nature. 2010;463:958–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Eng L, Coutinho G, Nahas S, Yeo G, Tanouye R, Babaei M, et al. Nonclassical splicing mutations in the coding and noncoding regions of the ATM gene: maximum entropy estimates of splice junction strengths. Hum Mutat. 2004;23:67–76.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Moon SD, Park JH, Kim EM, Kim JH, Han JH, Yoo SJ, et al. A novel IVS2-1G> a mutation causes aberrant splicing of the HRPT2 gene in a family with hyperparathyroidism-jaw tumor syndrome. J Clin Endocrinol Metab. 2005;90:878–83.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Simpson MA, Hsu R, Keir LS, Hao J, Sivapalan G, Ernst LM, et al. Mutations in FAM20C are associated with lethal Osteosclerotic bone dysplasia (Raine syndrome), highlighting a crucial molecule in bone development. Am J Hum Genet. 2007;81:906–12.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Gelincik A, Demir S, Olgaç M, Karaman V, Toksoy G, Çolakoğlu B, et al. Idiopathic angioedema with F12 mutation: is it a new entity? Ann Allergy Asthma Immunol. 2015;114:154–6.

    Article  PubMed  Google Scholar 

  47. 47.

    Müller CR, Fregin A, Srsen S, Srsnova K, Halliger-Keller B, Felbor U, et al. Allelic heterogeneity of alkaptonuria in Central Europe. Eur J Hum Genet. 1999;7:645–51.

    Article  PubMed  Google Scholar 

  48. 48.

    Bernard G, Chouery E, Putorti ML, Tétreault M, Takanohashi A, Carosso G, et al. Mutations of POLR3A encoding a catalytic subunit of RNA polymerase pol III cause a recessive Hypomyelinating Leukodystrophy. Am J Hum Genet. 2011;89:415–23.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Villate O, Ibarluzea N, Fraile-Bethencourt E, Valenzuela A, Velasco EA, Grozeva D, et al. Functional analyses of a novel splice variant in the CHD7 gene, found by next generation sequencing, confirm its pathogenicity in a Spanish patient and diagnose him with CHARGE syndrome. Front Genet. 2018;9:1–6.

    CAS  Article  Google Scholar 

  50. 50.

    Daubresse G, Deuring R, Moore L, Papoulas O, Zakrajsek I, Waldrip WR, et al. The Drosophila kismet gene is related to chromatin-remodeling factors and is required for both segmentation and segment identity. Development. 1999;126:1175–87.

    CAS  PubMed  Google Scholar 

  51. 51.

    Brehm A, Tufteland KR, Aasland R, Becker PB. The many colours of chromodomains. BioEssays. 2004;26:133–40.

    CAS  Article  PubMed  Google Scholar 

Download references


The authors are very thankful to the family for their contribution to this study.


This work was funded by the Research Council of Lithuania [No. S-MIP-17-19/LSS-150000-1179, Ingenes project]. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information




ES wrote the manuscript with consultation and support from EP and VM. ES extracted RNA of the proband and her healthy parents’ samples and performed reverse transcription reaction. ES and GP carried out genetic experiments and analysed obtained results. ES and TR performed bioinformatic analysis. ŽM and VM was involved in planning the experiment. VM and AM performed the WES experiment. JB provided primary medical counselling to the proband and was responsible for the temporal bone computed tomography procedure. AU and VK supervised the findings of this work. EP provided a genetic counselling to the family, conceived the original idea, and is a leader of Ingenes project, which funded this work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Evelina Siavrienė.

Ethics declarations

Ethics approval and consent to participate

The blood has been collected from patient and her parents, after obtaining written informed consent at the Centre for Medical Genetics, Vilnius University Hospital Santaros Klinikos. All procedures performed in studies involving human participants were in accordance with the 1964 Helsinki declaration and its later amendments. The study was approved by the Vilnius Regional Biomedical Research Ethics Committee of Lithuania.

Consent for publication

The patient’s parents provided the written consent for the case report to be published. Informed consent from parents regarding photographs for publication has been obtained.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. The conditions for PCR amplification of CHD7 gene in gDNA and cDNA samples. (DOCX 55 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Siavrienė, E., Petraitytė, G., Mikštienė, V. et al. A novel CHD7 variant disrupting acceptor splice site in a patient with mild features of CHARGE syndrome: a case report. BMC Med Genet 20, 127 (2019).

Download citation


  • CHD7
  • C.5535-1G > A splice site variant
  • CHARGE syndrome
  • Congenital anomalies
  • cDNA analysis