- 1Department of Neurology and Institute of Neurology, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China
- 2Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
- 3Department of Neuromuscular Disorders, The Third Hospital of Hebei Medical University, Shijiazhuang, China
- 4Department of Neurology, Luoyang Central Hospital Affiliated to Zhengzhou University, Luoyang, China
- 5Fujian Key Laboratory of Molecular Neurology, Fujian Medical University, Fuzhou, China
Background: Hereditary spastic paraplegia (HSP) caused by mutations in ALDH18A1 have been reported as spastic paraplegia 9 (SPG9), with autosomal dominant and autosomal recessive transmission (SPG9A and SPG9B). SPG9 is rare and has shown phenotypic and genotypic heterogeneity in previous reports.
Methods: This study screened ALDH18A1 mutations in autosomal recessive HSP patients using combined whole exome sequencing and RNA splicing analysis. We conducted in silico investigations, co-segregation analysis, and ELISA-based analysis of P5CS (Δ1-pyrroline-5-carboxylate synthetase; encoded by ALDH18A1) concentration to validate the pathogenicity of the detected ALDH18A1 variants. All previously reported bi-allelic ALDH18A1 mutations and cases were reviewed to summarize the genetic and clinical features of ALDH18A1-related HSP.
Results: A novel missense mutation c.880T>C, p.S294P and an intronic splicing mutation c.-28-13A>G were both detected in ALDH18A1 in an autosomal recessive family presenting with a complicated form HSP. ELISA assays revealed significantly decreased P5CS concentration in the proband's plasma compared with that in the healthy controls. Moreover, review of previously reported recessive cases showed that SPG9B patients in our cohort presented with milder symptoms, i.e., later age at onset and without cognitive impairment.
Conclusion: The present study expands the genetic and clinical spectrum of SPG9B caused by ALDH18A1 mutation. Our work defines new genetic variants to facilitate future diagnoses, in addition to demonstrating the highly informative value of splicing mutation prediction in the characterization of disease-related intronic variants.
Introduction
Hereditary spastic paraplegia (HSP) is an umbrella group comprised of neurodegenerative diseases characterized by progressive weakness and spasticity of the lower extremities (1). HSP caused by mutations in ALDH18A1 (MIM#616586), encoding delta-1-pyrroline-5-carboxylate synthase (P5CS), has been reported as hereditary spastic paraplegia type 9 (SPG9), i.e., SPG9A and SPG9B, which are distinguished by autosomal dominant and autosomal recessive modes of inheritance, respectively (2, 3). The P5CS enzyme is known to catalyze proline and ornithine biosynthesis, and disease-related mutations in ALDH18A1 have been correlated with loss or decrease of P5CS function (4). Patients with SPG9 present with either a pure or complicated form of HSP (3). Pure forms are characterized by lower limb spasticity, without prominent additional clinical findings other than mild urinary symptoms and impaired distal vibratory sensation (5). Complicated forms are associated with additional, often severe, clinical features such as cognitive impairment, seizures, neuropathy, amyotrophy, short stature, and vision abnormalities, among others (6).
HSP exhibits considerable genetic heterogeneity with more than 80 different HSP (or SPG) loci reported, including ALDH18A1 (2, 7). Massively parallel high throughput sequencing strategies, such as targeted gene-panel sequencing, whole exome sequencing (WES), and whole genome sequencing (WGS), have facilitated genetic diagnosis of HSP (8–11). As a first-line test, WES can significantly reduce the time and cost of the diagnostic process through accelerated identification of the causative variants associated with HSP. However, the rate of successful diagnoses may only reach ~50% in carefully selected cohorts of patients (12, 13). Unbiased genome-wide association studies have revealed that non-coding regions can play a significant role in gene regulation, and may account for 90% of causal disease loci in human complicated diseases (13–15). Penetrant non-coding variants that disrupt the normal mRNA splicing patterns, despite lying outside the canonical GC-AG splice sites, have long been recognized to play a significant role in rare genetic diseases (16). However, intronic splice mutations are often overlooked in WES sequence data due to the difficulty in accurate identification of these splice-altering variants outside the canonical GT-AG splice sites.
In this study, we report the identification of a novel splicing mutation c.-28-13A>G in intron 1 and a missense mutation c.880T>C, p.S294P in the ALDH18A1 gene in an autosomal recessive family presenting with a complicated form HSP. This is the first time that a hidden intronic mutation has been described for SPG9B. The mutation disrupts the normal accepter site, leading to the creation of an aberrant transcript that lacks exon 2, in which the original, correct AUG start codon is located. Thus, internal translation may initiate at downstream AUG start sites, thus leading to hypomorphic protein, and potentially resulting in mRNA degradation due to non-sense-mediated mRNA decay. Our results demonstrate that the genetic causes of Mendelian-inheritance diseases can be determined by focusing on non-coding DNA sequences.
Materials and Methods
Patients and Samples
The proband was recruited and examined at the Department of Neurology of the First Affiliated Hospital of Fujian Medical University. Detailed neurological examinations were performed on all members of the family by at least two senior neurologists. Blood samples were collected according to standard methods. Written informed consent was obtained from all participants in this study. The research project was approved by the Ethics Committee of the First Affiliated Hospital of Fujian Medical University.
Whole Exome Sequencing
Genomic DNA was extracted from the peripheral blood cells of patients and normal controls using Qiagen kits. Whole exome DNA was captured using a SureSelect Human All Exon V6 kit (Agilent) and sequenced using the Illumina HiSeq 3000 platform as described in our previous study (17). Fragment sequences were aligned to the consensus sequence (UCSC hg38). Variant calling was performed by Genome Analysis Toolkit (GATK) (18) and annotated with ANNOVAR (19). We filtered the SNP sets obtained from the total sequencing data using the parameters of GQ>20 and DP>8 for genotype quality (GQ) and read depth (DP), respectively. Variants were excluded based on the following criteria: (a) the variants did not affect the amino acid sequence; (b) the allele frequency was >1% in the 1,000 Genomes Project, ESP database, or gnomAD. Then we screened potential disease-causing genes according to their patterns of inheritance.
RNA Splicing Analysis
Total RNA was isolated from the peripheral blood leukocytes of family members. RNA was extracted using a Trizol extraction kit (Invitrogen, Carlsbad, CA, USA) and then synthesized to cDNA with a PrimeScript RT reagent kit (TAKARA BIO, Kusatsu, Shiga, Japan) according to the manufacturer's protocols. Primers were designed to target ALDH18A1 exons 1–3 to confirm aberrant splicing (forward:5′-ACGGAAGAAAAAAGAGAGTGAG-3′; reverse: 5′-AAGGACTTGCCATGTGTACG-3′). The products amplified by the above primers were separated by agarose gel electrophoresis. The DNA fragments were then purified from the gels and sent for Sanger sequencing.
High-Throughput Sequencing
The products amplified by the above primers were used for high-throughput sequencing (HTS). The detailed methods were described in the previous literature (20). In brief, the sequencing library was constructed (New England Biolabs), quality checked, and sequenced with paired-end 150-bp reads using the Illumina HiSeq X-Ten platform (Novogene). Paired end reads of 150 bp were pre-processed to remove adaptors and low quality reads using fastp (v 0.19.6) (21). The following criteria were used to remove the low-quality reads: (i) reads containing more than 15 N bases, (ii) a quality score of < 20, and (iii) more than 50% of the sequenced bases have low quality scores. Clean reads were subsequently aligned to the appropriate reference sequence using the Smith-Waterman algorithm. Mutation types were extracted from the aligned results and the number of all types of mutations were also counted. Subsequently, percentage of diverse transcripts was calculated as the number of reads align to mutant sequence over the total number of reads covering the same locus.
Enzyme-Linked Immunosorbent Assay
The blood samples from the patient carrying ALDH18A1 mutations and three gender matched healthy controls were centrifuged and the plasma was fractionated prior to storage at −80°C. The concentration of P5CS, encoded by ALDH18A1, was measured by human enzyme linked immunosorbent assay (ELISA) (Shanghai Jianglai Biotech, JL15386). The principle of ELISA kit was based on a double-antibody sandwich enzyme-linked immunosorbent assay to measure the concentration of ALDH18A1 by one step. Briefly, the plasma and the horseradish peroxidase (HRP)-labeled detection antibody were added to the microwells, which were pre-coated with ALDH18A1 antibody. Then, the mixture was incubate at 37°C for an hour. After washing with 1 × PBS, 3,3',5,5'-Tetramethylbenzidine was added as the chromogenic substrate. The wavelength at 450 nm excitation was used to assess the concentration of ALDH18A1 in the subject samples. The exact values were determined to the standard curve.
Statistical Analysis
Data are presented as mean ± SEM. One-way ANOVA with Dunnett's multiple comparisons test was used and probability values P < 0.05 were considered statistically significant. The statistical analyses were performed with GraphPad Prism 5.0.
Results
Clinical Description
The patient II-2 was a 46-year-old woman in ongoing follow-up at the neurology outpatient clinic of the First Affiliated Hospital of Fujian Medical University due to a 20-year-long history of altered gait (Table 1). She was born in a non-consanguineous family and showed normal motor and psychological milestones. She experienced unstable walking and abnormal walking posture at the age of 26 and symptoms worsened after pregnancy. As gait abnormality progressed, she developed dysarthria, although no muscular atrophy or autonomic symptoms were present. No symptoms of cutis laxa were observed. Physical examination revealed decreased muscle strength and increased muscle tone of the lower limbs, hyperreflexia in all limbs, and bilaterally extensor plantar responses. Hoffmann sign was observed in her upper limbs. Spinal MRI presented atrophy of thoracic spinal cord (Supplementary Figure 1). Moreover, the plasma levels of amino acids, including proline, arginine, ornithine, and citrulline, were within the normal range (Supplementary Table 1). The proband had a younger brother, subject II-3, who carried the same mutations. He presented with difficulty in walking in his 30's. Walking abnormalities progressed gradually, but no other symptoms emerged. He can still walk independently at the age of 44. Their parents and elder sister carried the heterozygous mutation without detectable symptoms (Figure 1A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627531/fneur-12-627531-HTML-r1/image_m/fneur-12-627531-g001.jpg)
Figure 1. ALDH18A1 variants in the family and molecular studies. (A) Family pedigree and genotype data for the ALDH18A1 variants. Filled black squares and circles indicate affected individuals. The proband in the family is indicated by an arrow. The genotypes of all available family members were determined. Dash symbol indicates reference allele. (B) Sequence analysis of the patients' genomic DNA revealed a splice site variant c.-28-13A>G in intron1 and a missense mutation c.880T>C. (C) Agarose gel electrophoresis of RT-PCR products showed an additional fragment (193 bp) in patients II-2 and II-3, but not in healthy family members. Sanger sequencing showed two populations of mRNA, one corresponding to skipped transcription of exon 2. (D) Transcriptional consequences of the c.-28-13A> G variant. The schematic diagram above shows the structure of the wild-type ALDH18A1 transcript (exons 1–3). The following schematic diagram shows the structure of the ALDH18A1 transcript isoform generated by the c.-28-13A> G variant (exon 2 skipping). (E) The plasma P5CS concentration of the proband decreased as compared with that of three, gender matched, healthy controls. Data represent mean ± SEM, *p < 0.05.
Molecular Studies
We sequenced the exomes of the proband and her parents to 100X average depth in target regions, with 96% of the target base pairs on average having at least 20X coverage. WES analysis identified a heterozygous missense variant in ALDH18A1 (chr10: 95628421A>G hg38, NM_002860: exon8: c.880T>C: p.S294P), which was inherited from her mother (Figure 1B). This variant was completely absent from the ExAC, 1,000 Genomes Project, and Genome Aggregation Database (gnomAD) genome and exome databases. It also showed high evolutionary conservation (GERP score: 5.79), was predicted to be deleterious by at least three different algorithms (PolyPhen-2, MutationTaster, CADD), and was found to be located in the G5K domain of P5CS. Based on these results, we classified the variant as “Likely Pathogenic” according to ACMG/AMP guidelines (PM1, PM2, PM3, and PP1) (26).
As the clinical presentation of the proband was compatible with spastic paraplegia, and we did not find any other rare candidate variants in either the exon or conventional splicing region in the exome sequence, we then re-examined the exome data and found a potentially deleterious variant in the intron 1. This variant (chr10: 95653418T>C hg38, NM_002860: c.-28-13A>G) was intronic, completely absent from ExAC and 1,000 Genomes Project, occurred at extremely low frequencies among East Asian populations according to the gnomAD database (1/5,196), and has not been previously reported. This variant was inherited from her father and segregation in the family was compatible with pathogenicity of the variant. In silico analysis using SpliceAI, a novel tool for cryptic splice site prediction based on artificial neural network deep learning (15), predicted the disruption of an acceptor site located 13 bp downstream of the c.-28-13A>G variant. The loss of the original acceptor splice site likely led to the transcriptional skipping of exon 2.
We therefore sought to validate the SpliceAI prediction through analysis of cDNA reverse transcribed from RNA extracted from the peripheral blood samples of patients and family members. Two transcripts were observed in gel electrophoresis and Sanger sequencing confirmed the prediction that exon 2 was absent in the truncated transcript due to transcriptional skipping caused by the c.-28-13A>G variant (Figure 1C). The percentage of truncated transcript was 14%, which was quantified using HTS (Supplementary Table 2). This novel mRNA isoform lacked the original start codon necessary for initiation of full length ALDH18A1 transcription, instead initiating mRNA transcription at the second possible start codon, in a different downstream reading frame (Figure 1D). Each of these reading frames, i.e., for the correct and aberrant start codons, resulted in a different amino acid sequence, and hence one of which was non-functional P5CS protein.
Functional Analysis of ALDH18A1 Mutations
The P5CS enzyme is known to catalyze proline and ornithine biosynthesis. The effects of disease-related mutations in ALDH18A1 are correlated with loss or decrease of P5CS function (2, 24, 27). To validate the functional involvement of ALDH18A1 in SPG9 in this family, we next performed ELISA analysis of the proband's plasma. The results showed that plasma P5CS concentration in patient II-2 was significantly decreased compared to that of the three, gender matched healthy controls (p < 0.05) (Figure 1E). Based on this functional validation, we proposed that compound heterozygosity for these two mutations most likely explains the phenotype of the patient.
Discussion
It has been estimated that genetic variants that affect splicing comprise roughly 9–10% of the pathogenic variants in rare genetic disorders (15, 28, 29), although this number appears very likely to be an underestimation. To date, disease-causing intronic splicing mutations have been reported in only a few genetic disorders such as ataxia-telangiectasia and neurofibromatosis type I (30, 31). Recently, a deep intronic splice mutation causing SPG7 was identified through detection by WGS, which led to a diagnosis of HSP (9). This discovery shed light on the importance of studying mutations at both the DNA and RNA levels. A number of in silico prediction tools have been developed to assess the effects of DNA sequence variations on splicing, which have substantially improved diagnosis of several human genetic diseases (32). In this work, the pathogenicity of intronic splice mutation was predicted through Illumina SpliceAI (15). Based on a deep learning neural network, SpliceAI can predict splice junctions from random pre-mRNA transcript sequences and generate anticipated start and end positions of aberrant exons with very high accuracy. The SpliceAI-predicted Delta scores and the position of the lost accepter for the variant (-28-13A>G) were 0.41 and −13, respectively, which indicated that exon2 was skipped (Figures 1C,D). The predicted length of the aberrant exon matched its observed length in the patient-derived cDNA sequence data.
In this study, we report an SPG9B family in which WES coupled with cDNA sequencing led to the identification of causative mutations. To our knowledge, this is the first description of a hidden intronic mutation in ALDH18A1, outside the essential GT and AG splice site dinucleotides, and moreover, only eight families with bi-allelic mutations in ALDH18A1 have been described to date. Clinical details for these patients and the previously reported SPG9B patients are presented in Table 1. Mutations in ALDH18A1 were initially identified in autosomal recessive or dominant neurocutaneous syndrome, characterized by severe developmental delay with marked cognitive impairment, associated with progeroid features, cutis laxa, joint hyperlaxity, short stature, cataract, and frequent microcephaly (2, 33, 34). Pyramidal signs have been found in patients with neurocutaneous syndrome, although not cardinal. Subsequently, it became clear that ALDH18A1 was the causal gene in some families with autosomal dominant or recessive HSP (SPG9A and SPG9B) (2). Pyramidal syndrome was reported as severe in SPG9 patients, while cutis laxa has never been noted, although it was previously considered to be an obligate sign of neurocutaneous syndrome (Table 1). Moreover, the clinical features of patients with SPG9B largely differ from those with SPG9A: they generally exhibit earlier onset and higher severity than SPG9A patients (27). In addition, unlike SPG9A, SPG9B patients share common features with those previously reported for neurocutaneous syndromes: developmental delay (7/14), intellectual deficiency (13/14), short stature (7/14), and facial dysmorphism (9/14) (2). Interestingly, these syndromes were demonstrated in 9/10 Caucasian patients, but none of which were observed in our patients. Notably, the clinical presentation was not homogeneous among the respective families, regardless of differences in their ethnicity. Age, disease duration, lifestyle, and genetic differences among ethnic groups are all potential modifiers underlying this phenotypic variability.
A spectrum of 13 mutations in the ALDH18A1 gene have been identified in 17 patients with autosomal recessive HSP, including 10 missenses, 2 splicing mutations, and a non-sense mutation (Table 1). Most notably, the previously reported p.R128H mutation was found in seven patients (2, 3, 25), making it the most frequently occurring ALDH18A1 mutation thus far identified in SPG9B patients. The patients reported to carry the p.R128H mutation invariably exhibited an earlier age of onset. However, truncating mutations, unlike other mutation types, have not been associated with an early age at onset or severe presentation. Moreover, genotype–phenotype correlations have been difficult to establish due to the complexity of both the phenotype and the function of the gene. Variability in the age at onset ranges from <1 to 32 years in SPG9B patients and no significant correlation has been established between phenotype and mutational class. It is therefore possible that other genotype–phenotype correlations will emerge as larger numbers of patients are studied.
Conclusion
With this report, we expand the genotype and phenotype spectrum of SPG9B and provide further evidence that RNA splicing assays, in conjunction with clinical and family inheritance analyses, can determine the clinical significance of intronic variants. In addition, the identification of this splicing mutation within intronic sequence should prompt closer scrutiny of ALDH18A1 non-coding genomic sequence in HSP patients for whom coding mutations in HSP-causative genes have not been identified.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession numbers can be found at: NCBI GenBank; PRJNA699107 and PRJNA699100.
Ethics Statement
The studies involving human participants were reviewed and approved by the First Affiliated Hospital of Fujian Medical University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Conception and drafting of the work were performed by Y-JC and XL. Data acquisition and data interpretation were performed by Y-JC, Z-QZ, M-WW, Y-SQ, R-YY, E-LD, ZZ, and H-TZ. Revision of the manuscript for intellectual content was conducted by XL, NW, and W-JC. All authors have read and approved the final manuscript.
Funding
This work was funded by National Natural Science Foundation of China Grants 81771230, 81801130, and U1905210; Joint Funds for the Innovation of Science and Technology of Fujian Province Grants 2017Y9094 and 2018Y9082; Key Clinical Specialty Discipline Construction Program of Fujian grant; Natural Science Foundation of Fujian Province Grant 2019J05076; Startup Fund for scientific research of Fujian Medical University Grants 2019QH2028.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We sincerely thank the patients and relatives for their participation in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.627531/full#supplementary-material
References
1. Hedera P. Hereditary spastic paraplegia overview. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Mirzaa G, Amemiya A, editors. GeneReviewsⓇ [Internet]. Seattle, WA: University of Washington (2020).
2. Coutelier M, Goizet C, Durr A, Habarou F, Morais S, Dionne-Laporte A, et al. Alteration of ornithine metabolism leads to dominant and recessive hereditary spastic paraplegia. Brain. (2015) 138:2191–205. doi: 10.1093/brain/awv143
3. Koh K, Ishiura H, Beppu M, Shimazaki H, Ichinose Y, Mitsui J, et al. Novel mutations in the ALDH18A1 gene in complicated hereditary spastic paraplegia with cerebellar ataxia and cognitive impairment. J Hum Genet. (2018) 63:1009–13. doi: 10.1038/s10038-018-0477-0
4. Pérez-Arellano I, Carmona-Alvarez F, Martínez AI, Rodríguez-Díaz J, Cervera J. Pyrroline-5-carboxylate synthase and proline biosynthesis: from osmotolerance to rare metabolic disease. Protein Sci. (2010) 19:372–82. doi: 10.1002/pro.340
5. Lee S, Park H, Zhu PP, Jung SY, Blackstone C, Chang J. Hereditary spastic paraplegia SPG8 mutations impair CAV1-dependent, integrin-mediated cell adhesion. Sci Signal. (2020) 13:eaau7500. doi: 10.1126/scisignal.aau7500
6. Fink JK. Hereditary spastic paraplegia: clinico-pathologic features and emerging molecular mechanisms. Acta Neuropathol. (2013) 126:307–28. doi: 10.1007/s00401-013-1115-8
7. Shribman S, Reid E, Crosby AH, Houlden H, Warner TT. Hereditary spastic paraplegia: from diagnosis to emerging therapeutic approaches. Lancet Neurol. (2019) 18:1136–46. doi: 10.1016/S1474-4422(19)30235-2
8. Dong EL, Wang C, Wu S, Lu YQ, Lin XH, Su HZ, et al. Clinical spectrum and genetic landscape for hereditary spastic paraplegias in China. Mol Neurodegener. (2018) 13:36. doi: 10.1186/s13024-018-0269-1
9. Verdura E, Schlüter A, Fernández-Eulate G, Ramos-Martín R, Zulaica M, Planas-Serra L, et al. A deep intronic splice variant advises reexamination of presumably dominant SPG7 Cases. Ann Clin Transl Neurol. (2020) 7:105–11. doi: 10.1002/acn3.50967
10. Schüle R, Wiethoff S, Martus P, Karle KN, Otto S, Klebe S, et al. Hereditary Spastic paraplegia: clinicogenetic lessons from 608 patients. Ann Neurol. (2016) 79:646–58. doi: 10.1002/ana.24611
11. Lin X, Su HZ, Dong EL, Lin XH, Zhao M, Yang C, et al. Stop-gain mutations in UBAP1 cause pure autosomal-dominant spastic paraplegia. Brain. (2019) 142:2238–52. doi: 10.1093/brain/awz158
12. Kara E, Tucci A, Manzoni C, Lynch DS, Elpidorou M, Bettencourt C, et al. Genetic and phenotypic characterization of complex hereditary spastic paraplegia. Brain. (2016) 139:1904–18. doi: 10.1093/brain/aww111
13. Farh KK, Marson A, Zhu J, Kleinewietfeld M, Housley WJ, Beik S, et al. Genetic and epigenetic fine mapping of causal autoimmune disease variants. Nature. (2015) 518:337–43. doi: 10.1038/nature13835
14. Maurano MT, Humbert R, Rynes E, Thurman RE, Haugen E, Wang H, et al. Systematic localization of common disease-associated variation in regulatory DNA. Science. (2012) 337:1190–5. doi: 10.1126/science.1222794
15. Jaganathan K, Kyriazopoulou Panagiotopoulou S, McRae JF, Darbandi SF, Knowles D, Li YI, et al. Predicting splicing from primary sequence with deep learning. Cell. (2019) 176:535–48.e24. doi: 10.1016/j.cell.2018.12.015
16. Cooper TA, Wan L, Dreyfuss G. RNA and disease. Cell. (2009) 136:777–93. doi: 10.1016/j.cell.2009.02.011
17. Chen YJ, Wang MW, Dong EL, Lin XH, Wang N, Zhang ZQ, et al. Chinese patients with adrenoleukodystrophy and Zellweger spectrum disorder presenting with hereditary spastic paraplegia. Parkinsonism Relat Disord. (2019) 65:256–60. doi: 10.1016/j.parkreldis.2019.06.008
18. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. (2010) 20:1297–303. doi: 10.1101/gr.107524.110
19. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. (2010) 38:e164. doi: 10.1093/nar/gkq603
20. Lin X, Chen H, Lu YQ, Hong S, Hu X, Gao Y, et al. Base editing-mediated splicing correction therapy for spinal muscular atrophy. Cell Res. (2020) 30:548–50. doi: 10.1038/s41422-020-0304-y
21. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560
22. Steenhof M, Kibæk M, Larsen MJ, Christensen M, Lund AM, Brusgaard K, et al. Compound heterozygous mutations in two different domains of ALDH18A1 do not affect the amino acid levels in a patient with hereditary spastic paraplegia. Neurogenetics. (2018) 19:145–9. doi: 10.1007/s10048-018-0547-7
23. Wei Q, Dong HL, Pan LY, Chen CX, Yan YT, Wang RM, et al. Clinical features and genetic spectrum in Chinese patients with recessive hereditary spastic paraplegia. Transl Neurodegener. (2019) 8:19. doi: 10.1186/s40035-019-0157-9
24. Magini P, Marco-Marin C, Escamilla-Honrubia JM, Martinelli D, Dionisi-Vici C, Faravelli F, et al. P5CS expression study in a new family with ALDH18A1-associated hereditary spastic paraplegia SPG9. Ann Clin Transl Neurol. (2019) 6:1533–40. doi: 10.1002/acn3.50821
25. Kalmár T, Maróti Z, Zimmermann A, Sztriha L. Tremor as an early sign of hereditary spastic paraplegia dueto mutations in ALDH18A1. Brain Dev. (2020) 43:144–51. doi: 10.1016/j.braindev.2020.07.015
26. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, 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. (2015) 17:405–24. doi: 10.1038/gim.2015.30
27. Panza E, Escamilla-Honrubia JM, Marco-Marín C, Gougeard N, De Michele G, Morra VB, et al. ALDH18A1 gene mutations cause dominant spastic paraplegia SPG9: loss of function effect and plausibility of a dominant negative mechanism. Brain. (2016) 139:e3. doi: 10.1093/brain/awv247
28. Anna A, Monika G. Splicing mutations in human genetic disorders: examples, detection, and confirmation. J Appl Genet. (2018) 59:253–68. doi: 10.1007/s13353-018-0444-7
29. Hiraide T, Nakashima M, Ikeda T, Tanaka D, Osaka H, Saitsu H. Identification of a deep intronic POLR3A variant causing inclusion of a pseudoexon derived from an Alu element in Pol III-related leukodystrophy. J Hum Genet. (2020) 65:921–5. doi: 10.1038/s10038-020-0786-y
30. Teraoka SN, Telatar M, Becker-Catania S, Liang T, Onengüt S, Tolun A, et al. Splicing defects in the ataxia-telangiectasia gene, ATM: underlying mutations and consequences. Am J Hum Genet. (1999) 64:1617–31. doi: 10.1086/302418
31. Ars E, Serra E, García J, Kruyer H, Gaona A, Lázaro C, et al. Mutations affecting mRNA splicing are the most common molecular defects in patients with neurofibromatosis type 1. Hum Mol Genet. (2000) 9:237–47. doi: 10.1093/hmg/9.2.237
32. Zhao M, Chen YJ, Wang MW, Lin XH, Dong EL, Chen WJ, et al. Genetic and clinical profile of Chinese patients with autosomal dominant spastic paraplegia. Mol Diagn Ther. (2019) 23:781–9. doi: 10.1007/s40291-019-00426-w
33. Fischer-Zirnsak B, Escande-Beillard N, Ganesh J, Tan YX, Al Bughaili M, Lin AE, et al. Recurrent de novo mutations affecting residue Arg138 of pyrroline-5-carboxylate synthase cause a progeroid form of autosomal-dominant cutis laxa. Am J Hum Genet. (2015) 97:483–92. doi: 10.1016/j.ajhg.2015.08.001
34. Baumgartner MR, Hu CA, Almashanu S, Steel G, Obie C, Aral B, et al. Hyperammonemia with reduced ornithine, citrulline, arginine and proline: a new inborn error caused by a mutation in the gene encoding delta(1)-pyrroline-5-carboxylate synthase. Hum Mol Genet. (2000) 9:2853–8. doi: 10.1093/hmg/9.19.2853
Keywords: spastic paraplegia 9, ALDH18A1, intronic splicing mutation, SpliceAI, RNA splicing assay
Citation: Chen Y-J, Zhang Z-Q, Wang M-W, Qiu Y-S, Yuan R-Y, Dong E-L, Zhao Z, Zhou H-T, Wang N, Chen W-J and Lin X (2021) Novel Compound Missense and Intronic Splicing Mutation in ALDH18A1 Causes Autosomal Recessive Spastic Paraplegia. Front. Neurol. 12:627531. doi: 10.3389/fneur.2021.627531
Received: 09 November 2020; Accepted: 07 April 2021;
Published: 19 May 2021.
Edited by:
Guy Rouleau, McGill University, CanadaReviewed by:
Giovanni Stevanin, INSERM U1127 Institut du Cerveau et de la Moelle épinière (ICM), FranceNicolas Dupre, Laval University, Canada
Ziv Gan-Or, McGill University, Canada
Junling Wang, Central South University, China
Copyright © 2021 Chen, Zhang, Wang, Qiu, Yuan, Dong, Zhao, Zhou, Wang, Chen and Lin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiang Lin, linxiang1988@fjmu.edu.cn; Wan-Jin Chen, wanjinchen75@fjmu.edu.cn
†These authors have contributed equally to this work