Whole-exome sequencing and digital PCR identified a novel compound heterozygous mutation in the NPHP1 gene in a case of Joubert syndrome and related disorders
© The Author(s). 2017
Received: 21 June 2016
Accepted: 14 March 2017
Published: 27 March 2017
Joubert syndrome and related disorders (JSRD) is a clinically and genetically heterogeneous condition with autosomal recessive or X-linked inheritance, which share a distinctive neuroradiological hallmark, the so-called molar tooth sign. JSRD is classified into six clinical subtypes based on associated variable multiorgan involvement. To date, 21 causative genes have been identified in JSRD, which makes genetic diagnosis difficult.
We report here a case of a 28-year-old Japanese woman diagnosed with JS with oculorenal defects with a novel compound heterozygous mutation (p.Ser219*/deletion) in the NPHP1 gene. Whole-exome sequencing (WES) of the patient identified the novel nonsense mutation in an apparently homozygous state. However, it was absent in her mother and heterozygous in her father. A read depth-based copy number variation (CNV) detection algorithm using WES data of the family predicted a large heterozygous deletion mutation in the patient and her mother, which was validated by digital polymerase chain reaction, indicating that the patient was compound heterozygous for the paternal nonsense mutation and the maternal deletion mutation spanning the site of the single nucleotide change.
It should be noted that analytical pipelines that focus purely on sequence information cannot distinguish homozygosity from hemizygosity because of its inability to detect large deletions. The ability to detect CNVs in addition to single nucleotide variants and small insertion/deletions makes WES an attractive diagnostic tool for genetically heterogeneous disorders.
KeywordsCase report Joubert syndrome NPHP1 Deletion Hidden Markov model Whole exome sequencing Digital PCR
Joubert syndrome and related disorders (JSRD) (OMIM 213300) is a rare neurological disorder with autosomal recessive or X-linked inheritance, characterized by neurological symptoms including hypotonia, ataxia, global developmental delay, abnormal eye movements, and breathing dysregulation [1–3]. JSRD is also often associated with visceral involvements and is classified into six clinical subtypes . The neuroradiological hallmark of JSRD is the so-called molar tooth sign, which reflects a complex cerebellar and brainstem malformation [4, 5]. To date, 21 genes have been reported to be responsible for JSRD [1, 6]. JSRD has clinical and genetic heterogeneity, which makes it difficult to identify the causative mutations in individual cases by Sanger sequencing alone . Whole-exome sequencing (WES) is becoming widely adopted as an efficient strategy to identify disease-causing mutations in genetically heterogeneous diseases. In addition to single nucleotide variants (SNVs) and small insertion/deletions, structural changes such as copy number variations (CNVs) contribute to human diseases [8, 9]; however, the accurate detection of CNVs using WES data remains challenging. In this study, WES and a read depth-based CNV detection algorithm using WES data of the family identified not only a novel nonsense mutation but also a large heterozygous deletion mutation in NPHP1, leading to the precise molecular diagnosis of JSRD.
Genomic DNA was extracted from the peripheral blood of each subject using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany). WES of the patient and her family members (her parents and her healthy older sister) was performed using amplicon-based next-generation sequencing. Briefly, libraries were constructed using an Ion AmpliSeq Library Kit v2.0 (Life Technologies, Carlsbad, CA) according to the manufacturer’s instructions. Quantification of the libraries was performed by 2200 TapeStation Instrument using High Sensitivity D1000 Reagents and High Sensitivity D1000 ScreenTape (Agilent Technologies, Santa Clara, CA). Amplified libraries were submitted to emulsion polymerase chain reaction (PCR) using an Ion OneTouch™ 2 Instrument with an Ion PI™ Template OT2 200 Kit v3. Ion sphere particles were enriched using Ion OneTouch ES and were loaded on an Ion PI Chip v2. Sequencing was performed by Ion Proton™ with an Ion PI Sequencing 200 Kit v3. Read sequence files were then run through two independent pipelines constructed with Bowtie2 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml)-GATK (the Genome Analysis Toolkit; https://www.broadinstitute.org/gatk/index.php) and BWA (Burrows Wheeler Alignment; http://bio-bwa.sourceforge.net/bwa.shtml)-Platypus (http://www.well.ox.ac.uk/platypus) to obtain variant call format files. Then, concordant genetic variants detected by two pipelines were annotated using ANNOVAR program (http://annovar.openbioinformatics.org). Then, variants in 21 known causative genes of JSRD were picked up. To calculate the depth of coverage, we used the GenomicRanges and IRanges packages from Bioconductor (http://www.bioconductor.org), with a sex-considered correction using the average read counts of the X chromosome in each sample, in consideration of the number of X chromosomes. Read counts of the patient were also compared with those of her father as a reference sample in the family. Then, we ran a hidden Markov model algorithm to call CNVs according to the programs described in the eXome Hidden Markov Model (XHMM) .
To confirm the mutation identified by WES, Sanger sequencing was performed by amplifying NPHP1 exon 7 with the following primers: 5′-AATTCCCATCTATGCTCAAGATTT-3′ (forward) and 5′-TTCCCACTTTTGTACCTTTGC-3′ (reverse). PCR products were sequenced using a BigDyeV3.1 terminator Kit on an ABI 310 automated sequencer (Life Technologies).
Digital PCR was performed on a QuantStudio™ 3D Digital PCR System (Life Technologies) according to the manufacturer’s instructions. Briefly, reaction mixtures were prepared with the QuantStudio 3D Digital PCR Master Mix and custom TaqMan probes. Primers for NPHP1 exon 7 (Hs00138437_cn) and BUB1 exon 8 (Hs02420021_cn) were purchased from Life Technologies. BUB1 is located 500 kb upstream of NPHP1 and was used as a control. Reaction mixtures were loaded onto a QuantStudio 3D Digital PCR Chip, and PCR was performed using recommended conditions. The gene copy number was analyzed using a QuantStudio 3D Digital PCR Instrument.
A deletion of NPHP1 had previously been reported in JSRD [11, 12], so we next performed CNV assessment using a read depth-based CNV detection algorithm. The CNV-calling algorithm predicted the presence of an approximately 470 kb heterozygous deletion including the entire NPHP1 gene in both the patient and her mother (Additional file 2: Figure S1). To confirm this deletion mutation, digital PCR targeted to NPHP1 exon 7 and BUB1 exon 8 was performed. The copy numbers of NPHP1 (copies/μl) in individuals I-1, I-2, II-1, and II-2 were 245.2, 133.6, 243.7, and 122.6, respectively (Fig. 2b, left panel). This compares with those of BUB1: 235.7, 231.0, 223.4, and 221.2, respectively (Fig. 2b, middle panel). The copy number ratios of NPHP1 to BUB1 of I-1, I-2, II-1, and II-2 were 1.04, 0.58, 1.09, and 0.55, respectively (Fig. 2b, right panel), indicating the presence of a heterozygous deletion in individuals I-2 (the patient’s mother) and II-2 (the patient). These findings indicated that the patient was compound heterozygous for the paternal nonsense mutation and the maternal deletion mutation.
JSRD is a genetically heterogeneous disorder; therefore, the identification of the causative mutation for each patient by Sanger sequencing is a daunting task even after narrowing down the candidates based on clinical phenotype. As shown here, WES is also useful to detect not only single nucleotide changes but also CNVs that cannot be identified by Sanger sequencing. CNV assessment using WES data of the family led to the detection of a large deletion mutation in combination with a novel nonsense mutation of the remaining allele. In the present case, WES and Sanger sequencing of the patient detected the nonsense mutation in an apparently homozygous state; however, this was refuted by genetic testing to assess CNV. Therefore, it should be noted that analytical pipelines that focus purely on sequence information cannot distinguish homozygosity from hemizygosity because of its inability to detect large deletions. Thus, in this case, comprehensive analysis of the family was necessary to avoid imprecise genetic diagnosis.
CNV detection is an important issue in genetic diagnosis that might be overlooked by conventional sequencing strategies. Although several algorithms capable of identifying CNVs from WES data have been proposed, there is no WES-based CNV-calling method that can detect all CNVs in a large size range [13–15]. Because non-allelic homologous recombination between low-copy repeats on chromosome 2q13 have been shown to be a major cause of large deletions involving the NPHP1 gene [11, 16, 17], we chose to use XHMM based CNV-calling algorithm, which is suitable for predicting CNVs ranging from 10 kb to 1 Mb . Although we have not determined the precise breakpoint of the deletion, the CNV-calling algorithm was nevertheless successful at predicting the heterozygous large deletion including the entire NPHP1 gene. It should be noted that WES data could identify not only SNVs but also clinically relevant CNVs without additional experiments such as genome-wide array comparative genomic hybridization or single nucleotide polymorphism array.
CNV validation methods include real-time quantitative PCR and multiplex ligation-dependent probe amplification. In the present case, we performed chip-based digital PCR to confirm the NPHP1 deletion. This new PCR technique allows absolute quantification of copy numbers without reference to a standard curve, and produces precise results compared with conventional real-time quantitative PCR [18, 19]. We show that digital PCR is an alternative efficient method to confirm clinical relevant CNVs in genetic diagnosis.
In conclusion, we identified a novel compound heterozygous mutation (p.Ser219*/deletion) in the NPHP1 gene in a case of JSRD using WES. The present study showed a clinical usefulness of WES and digital PCR to detect clinically relevant CNVs in patients with genetically heterogeneous conditions. It should be noted that patients with genetic disorders might harbor CNVs that cannot be detected by conventional sequencing strategies.
Copy number variation
Joubert syndrome and related disorders
Magnetic resonance imaging
Online mendelian inheritance in man
Polymerase chain reaction
Single nucleotide variant
eXome hidden markov model
The authors thank the patient and her family members who participated in this study. We are grateful to Ms. T. Seino for her technical assistance.
This report did not require any funding source.
Availability of data and materials
All data supporting our findings are included in the manuscript.
SK acquired the clinical data. SK and HS carried out the molecular genetic studies. All authors analyzed data. SK wrote the initial draft of the manuscript which was edited by ME and TK. All authors have read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Written informed consent was obtained from the patient’s mother for publication of this Case report and any accompanying images. A copy of the written consent is available for review by the Editor of this journal.
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of the Yamagata University Faculty of Medicine. Written informed consent was obtained from the patient and her family members.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Romani M, Micalizzi A, Valente EM. Joubert syndrome: congenital cerebellar ataxia with the molar tooth. Lancet Neurol. 2013;12(9):894–905.View ArticlePubMedGoogle Scholar
- Brancati F, Dallapiccola B, Valente EM. Joubert Syndrome and related disorders. Orphanet J Rare Dis. 2010;5:20.View ArticlePubMedPubMed CentralGoogle Scholar
- Parisi MA. Clinical and molecular features of Joubert syndrome and related disorders. Am J Med Genet C Semin Med Genet. 2009;151C(4):326–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Maria BL, et al. "Joubert syndrome" revisited: key ocular motor signs with magnetic resonance imaging correlation. J. Child Neurol. 1997;12(7):423–30.View ArticleGoogle Scholar
- Poretti A, et al. Joubert syndrome and related disorders: spectrum of neuroimaging findings in 75 patients. AJNR Am J Neuroradiol. 2011;32(8):1459–63.View ArticlePubMedGoogle Scholar
- Valente EM, et al. Clinical utility gene card for: Joubert syndrome--update 2013. Eur J Hum Genet. 2013;21:10.View ArticleGoogle Scholar
- Tsurusaki Y, et al. The diagnostic utility of exome sequencing in Joubert syndrome and related disorders. J Hum Genet. 2013;58(2):113–5.View ArticlePubMedGoogle Scholar
- Cooper GM, et al. A copy number variation morbidity map of developmental delay. Nat Genet. 2011;43(9):838–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Pinto D, et al. Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010;466(7304):368–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Fromer M, et al. Discovery and statistical genotyping of copy-number variation from whole-exome sequencing depth. Am J Hum Genet. 2012;91(4):597–607.View ArticlePubMedPubMed CentralGoogle Scholar
- Parisi MA, et al. The NPHP1 gene deletion associated with juvenile nephronophthisis is present in a subset of individuals with Joubert syndrome. Am J Hum Genet. 2004;75(1):82–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Castori M, et al. NPHP1 gene deletion is a rare cause of Joubert syndrome related disorders. J Med Genet. 2005;42(2):e9.View ArticlePubMedPubMed CentralGoogle Scholar
- Samarakoon PS, et al. Identification of copy number variants from exome sequence data. BMC Genomics. 2014;15:661.View ArticlePubMedPubMed CentralGoogle Scholar
- Tan R, et al. An evaluation of copy number variation detection tools from whole-exome sequencing data. Hum Mutat. 2014;35(7):899–907.View ArticlePubMedGoogle Scholar
- Miyatake S, et al. Detecting copy-number variations in whole-exome sequencing data using the eXome Hidden Markov Model: an 'exome-first' approach. J Hum Genet. 2015;60(4):175–82.View ArticlePubMedGoogle Scholar
- Stankiewicz P, Lupski JR. Genome architecture, rearrangements and genomic disorders. Trends Genet. 2002;18(2):74–82.View ArticlePubMedGoogle Scholar
- Konrad M, et al. Large homozygous deletions of the 2q13 region are a major cause of juvenile nephronophthisis. Hum Mol Genet. 1996;5(3):367–71.View ArticlePubMedGoogle Scholar
- Huggett JF, Whale A. Digital PCR as a novel technology and its potential implications for molecular diagnostics. Clin Chem. 2013;59(12):1691–3.View ArticlePubMedGoogle Scholar
- Nuzzo F, et al. Identification of a novel large deletion in a patient with severe factor V deficiency using an in-house F5 MLPA assay. Haemophilia. 2015;21(1):140–7.View ArticlePubMedGoogle Scholar