Clinical summary
A 9-month-old girl was born after normal labor with a birth weight of 3.30 kg and height of 50 cm at 39 weeks gestation. Her parents were inconsanguineous. She is the second child of her family and has an elder sister. Her parents and sister are unaffected. The girl was hospitalized because of severe respiratory problems at birth, and then for delayed global development at 9 months. She can raise her head until 4 months, and can turn over until 8 months, cannot sat all by herself, and cannot stand when supporting station. She can grasp objects on her own initiative. She can identify her mother and speak with only monosyllables “mama, baba, et al.” She often stuck out her tongue unconsciously. The girl had mild retardation for language and fine movements, moderate retardation for adaptability and social abilities, and severe retardation for large motor abilities, as evaluated by the Chinese Developmental Scale for children aged 0–6 years (WS/T 580–2017). The girl was fed with milk powder. The girl had a height of 74 cm (+ 1 SD), weight of 9 kg (median ~ + 1 SD), and occipitofrontal circumference of 45 cm (+ 1 SD) at 9 months old. Physical examination showed strabismus, hypotonia and lower myodynamia. No abnormalities were mentioned for the heart, liver, gallbladder, spleen, kidneys, ureter, or bladder. Magnetic resonance imaging of her brain revealed that the superior cerebellar peduncle was thickened and lengthened (molar tooth sign), and the shape of the cerebellar vermis was irregular and partially missing (Fig. 1). There was a link near the midline of the bilateral cerebellar hemispheres. The upper part of the fourth ventricle seemed to be bat wings. Thus, the girl was diagnosed with Joubert syndrome.
Genetic investigations
A 2 mL EDTA anticoagulant venous blood sample was obtained from the proband, her parents and her elder sister. Genomic DNA was extracted from whole blood using QIAamp® DNA Blood mini Kit (QIAGEN, Germany) according to the manufacturer’s protocol. Whole-Exome Sequencing (WES) and SNVs/InDel analysis were performed for the proband as previously described[18]. Splice variants were evaluated using the online tool SpliceAI (https://spliceailookup.broadinstitute.org). Identification and annotation of genetically mobile domains and the analysis of domain architectures were performed by SMART (a Simple Modular Architecture Research Tool) (http://smart.embl-heidelberg.de/).
The WES data were further analyzed for copy number variants (CNVs). We made a virtual gene panel for JBTS and other ciliopathies using Phenolyzer[19] and OMIM (https://omim.org/) (There were a total of 169 related genes, data not shown). We obtained the chromosomal position information of exons of the classical transcripts from UCSC (https://genome.ucsc.edu/) for all candidate genes, and generated a bed file for all exons. Based on the hypothesis that exon deletion/duplication rarely occurs in JBTS, we performed CNV assessment using a read depth-based CNV detection algorithm[20]. Three WES data at the same sequencing batch, including the proband and two unaffected individuals were involved in CNV analysis. The mean depth of each exon was obtained with Samtools (https://anaconda.org/bioconda/samtools). Using Integrative Genomics Viewer (IGV, http://software.broadinstitute.org/software/igv/), the called exon deletions/duplications were confirmed manually.
SNV verification
PCR was performed using specific primer pairs for candidate SNV validation screened from WES in the proband and her family members. The Supplementary Table 1 lists the primer sequences.
CNV verification and precise breakpoint determination
To verify the CNV, we firstly carried out a relative quantitative PCR to analyze genome relative copy number of exon 15 of KIAA0586 (NM_001244189.1) in the family. Triplicate quantitative PCR for genomic DNA was performed using SYBR Green qPCR Master Mix (Life Technologies, USA) on ABI7500 (Life Technologies, USA). ALB (albumin) was used as an internal control. The delta-delta CT value analysis method was used to evaluate the relative copy number of genome KIAA0586 exon 15 for all samples. Target and internal control gene primer pairs are listed in Supplementary Table 1. Then, we performed a long PCR to determine the precise breakpoint of the deletion. The long PCR primer sequences are shown in Supplementary Table 1.
RNA extraction and reverse transcription
The Tempus™ Spin RNA Isolation Kit (Invitrogen, USA) was used for RNA extraction from whole blood cells of the proband and her parents and sister. One microgram RNA was reverse transcribed into cDNA using a SuperScript™ IV First-Strand Synthesis System Kit (Invitrogen, USA). The primer sequences for investigating the effect of deletion in the transcript are listed in Supplementary Table 1.
The PCR amplification products from the verification of SNVs and CNVs were analyzed by agarose gel electrophoresis, and then Sanger sequencing was performed on an ABI3730xl Genetic Analyzer (Life Technologies, USA) following the manufacturer's protocol.