- 1Department of Transfusion, Shaoxing People’s Hospital, Shaoxing, China
- 2Department of Medical Genetics, School of Basic Medical Sciences, Southern Medical University, Guangzhou, China
- 3Experimental Education/Administration Center, School of Basic Medical Science, Southern Medical University, Guangzhou, China
- 4Department of Hematology, Shaoxing People’s Hospital, Shaoxing, China
Gaucher disease (GD, ORPHA355) is a rare autosomal recessive genetic disease caused by mutations in GBA1, which encodes the lysosomal enzyme glucocerebrosidase (GCase). Here, we report a patient with GD who carried the heterozygous c.1240G > C (p.Val414Leu) mutation and the heterozygous pathogenic c.1342G > C (p.Asp448His) mutation in GBA1. Bioinformatics analysis suggested that the two mutations are pathogenic. Functional studies showed that GBA1 mRNA and GCase protein levels of mutant types were significantly less than the wild-type. In the cell lysates, the two mutations of GBA1 c.1240G > C and c.1342G > C caused a decreased GCase concentration, while the two mutations did not change the distribution in the cell. The pathogenicity of the compound heterozygous mutations was verified. Early diagnosis and treatment can improve the quality of life and prevent unnecessary procedures in patients with GD.
Introduction
Gaucher disease (GD, ORPHA355) is a rare autosomal recessive genetic disease, affecting 1/50,000 to 1/100,000 individuals in the general population. However, its frequency is as high as 1:800 in Ashkenazi Jews (1). GD is caused by mutations in GBA1 (OMIM #606463), which encodes the lysosomal enzyme glucocerebrosidase (GCase) (2). To date, more than 500 mutations have been described in the GBA1 gene, including single base changes, splicing alterations, partial and total deletions, insertions, and rearrangements (3, 4). GD is categorized into three clinical phenotypes based on the presence or absence of neurological involvement (5). A diagnosis of GD is established by the detection of GCase or by molecular genetics testing (6).
In this study, we report a patient with GD who carried the c.1240G > C (p.Val414Leu) mutation and c.1342G > C (p.Asp448His) mutation in GBA1. The pathogenicity of the compound heterozygous mutation was verified with cell experiments.
Materials and methods
Patient
The diagnosis was based on clinical data as well as reduced GCase activity and molecular analysis. This study conformed to human research guidelines stated in the Declaration of Helsinki and written informed consent was obtained from all subjects. The female patient was a 62-year-old who first visited Shaoxing People's Hospital because of massive splenomegaly and thrombocytopenia with no neurological impairment.
Bioinformatics analysis
The three-dimensional (3D) protein structures of the wild-type and mutant GCase were predicted using I-TASSER (http://zhanglab.ccmb.med.umich.edu/I-TASSER) and Pymol. ClustalX was used to analyze evolutionary conservation among different paralogs and orthologs.
DNA extraction, purification, and single-gene sequencing
DNA was extracted from fresh blood using a standard phenol-chloroform extraction method after digestion with proteinase K to ensure that the DNA quality and quantity met the experimental requirements. Specific primers were used to amplify the 8–11 exon of GBA1 gene. The forward and reverse sequences were as follows: forward primer (5′-TGTGTGCAAGGTCCAGGATCAG-3′), reverse primer (5′-ACCACCTAGAGGGGAAAGTG-3′). Sanger sequencing was performed by Sangon Biotech (Shanghai, China). The sequences obtained were aligned to the available reference sequence in The National Centre for Biotechnology Information GeneBank database to detect variants.
Plasmid construct
The length of the GBA1 cDNA is approximately 1.6 kb. Total RNA from peripheral blood lymphocytes of healthy family members was extracted with Trizol reagent (Invitrogen, Carlsbad, CA, United States) and was reverse transcribed into first-strand cDNA with the HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, Jiangsu, China) according to the manufacturer's instructions. The wild-type GBA1 cDNA was cloned using primers containing the XhoI and EcoRI restriction sites. The forward and reverse sequences were as follows: forward primer (5′-CCGCTCGAGAGATGGAGTTTTCAAG-3′),reverse primer (5′-GGAATTCCACTGGCGACGCCACAG-3′). To clone the verified PCR products into the pEGFP-N1 vector by double digestion with XhoI and EcoRI. GBA1 mutations were introduced into the wild-type plasmid by PCR-based site-directed mutagenesis to obtain the pEGFP-N1-GBA1-MUT vector (c.1240G > C:forward primer: 5′-AACCTCCTGTACCATCTGGTCGGCTGGACC-3′; reverse primer: 5′-GGTCCAGCCGACCAGATGGTACAGGAGGTT-3′; c.1342G > C: forward primer: 5′-GTAGACATCACC AAGCACACGTTTTACAAA-3′; reverse primer: 5′-TTTGTAAAACGTGTGCTTGGTGATGTCTAC-3′) with the Phanta Max Super-Fidelity DNA polymerase and confirmed by sequencing. The plasmids were extracted with the endotoxin-free plasmid mini-extraction kit (Tiangen, Beijing, China).
HEK 293 T cell line culture and transient transfection
HEK 293 T cells were grown in Dulbecco's Modified Eagle Medium (Gibco, NY, United States) supplemented with 10% fetal bovine serum (Gibco), 100 U/ml penicillin, and 100 U/ml streptomycin at 37 °C in a humidified atmosphere enriched with 5% CO2. Cell transfection was performed according to the manufacturer's instructions using Lipofectamine 3,000 transfection reagent (Invitrogen) for the pEGFP-N1 vector, and the GBA1 wild-type (WT) and mutant-type (MT), respectively. After 48 h transfection, cells were collected for GCase enzyme assays and Western blotting. All transfections were performed in duplicate in three individual experiments.
Immunofluorescence detection and subcellular localization
The cells were grown on Petri dishes and fixed with cold 4% formaldehyde for 30 min at −20 °C and blocked in 0.2% Triton X-100 for 5 min at room temperature. The fixed cells were separately stained with DAPI to visualize the nuclear DNA. Fluorescence microscopy was carried out with an LSM880 confocal microscope equipped with appropriate filter sets.
ELISA
When HEK 293 T cells reached 70% confluence in 6-well plates, 2.5 µg pEGFP-N1-GBA1-WT and MT were transfected into the cells using Lipofectamine 3,000 transfection reagent (Invitrogen). The cells were collected by centrifugation and washed with pre-chilled PBS after 48 h transfection, then 500 µl PBS was added to resuspend the cells, and the total cellular protein was extracted by repeated freezing and thawing. GCase concentrations were detected using the Human Glucosylceramidase (GCase) ELISA Kit (CUSABIO).
Gcase enzyme activity assay
When HEK 293 T cells reached 50% confluence in 6-well plates, 2.5 µg pEGFP-N1-GBA1-WT and MT were transfected into the cells using Lipofectamine 3,000 transfection reagent (Invitrogen). 48 h after transfection, cells were collected and washed with cold PBS (Sigma–Aldrich), and the total cell protein was extracted by repeated freezing and thawing before resuspension of the cells in 150 µl PBS containing 1% PMSF. Cell debris was removed by centrifugation, and cell lysates were quantified using the Pierce™ Rapid Gold BCA Protein Assay Kit (Thermo Fisher Scientific). 15 μg protein was mixed with 20 ul 4-methylumbelliferyl-β-glucosidase substrate, incubated at 37 °C for 2 h, and 200 ul glycine buffer was added to terminate the reaction. Sample 4-methylumbelliferone fluorescence intensity was measured using a microplate reader at an excitation wavelength of 355 nm and an emission wavelength of 460 nm.
RNA expression
Total RNA was extracted with TRIzol reagent (Invitrogen) from the transfected cells. RT-PCR was performed using the HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, Nanjing, China). GBA1 mRNA expression levels were analyzed using the ChamQ SYBR qPCR Master Mix (Vazyme) with the specific primers as follows: forward primer5′-TCCCAAGCCTTTGAGTAGGG-3′, reverse primer 5′-GCCGAAGCTTTTAGGGATGC-3′. Gene expression levels were calculated using the 2−ΔΔCT method. Expression levels of the target gene was normalized to that of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) with the specific primers: forward primer 5′-GTGAAGGTCGGAGTCAACG-3′, reverse primer 5′-TGAGGTCAATGAAGGGGTC-3′.
Western blotting
The transfected cells were lysed using cell lysis buffer for Western blot (Beyotime, Beijing, China) containing 1% PMSF. Protein loading buffer was added to the cell lysates before boiling for 5 min. Protein (25 µg) was loaded for electrophoresis on 10% SDS-PAGE gels and transferred to PVDF membranes. Then, the membranes were blocked with 5% non-fat milk in TBS containing 0.1% Tween 20 for 2 h at room temperature and incubated with antibodies against EGFP (1:2,000, Proteintech, Wuhan, China) at 4 °C overnight. Finally, the membranes were incubated with a secondary goat anti-mouse IgG (1:5,000) for 2 h at room temperature. The signals were detected using the ECL Western blotting detection reagent.
Statistical analysis
The results are presented as means ± standard deviation of at least three independent biological replicates. At least three technical replicates were used per experiment. The differences between the WT and MT enzyme activities were analyzed using a t-test with GraphPad Prism 5 (GraphPad Software Inc., CA, United States). A p-value of ≤0.05 was considered statistically significant.
Results
Clinical phenotype
The patient was diagnosed with Gaucher disease type I. She was conscious and had no anemic appearance and no ecchymosis or petechiae on the skin of limbs (Figure 1). There was no significant superficial lymphadenopathy. Her breath sounds were clear, without dry or wet sounds. The heart rate was 65 beats/min and was rhythmical. The abdomen was flat and soft, without pressure pain or rebound pain. There was no percussion pain in either kidney and no swelling in the lower limbs. The neurological examination was normal. Routine blood tests showed normal hemoglobin and mild thrombocytopenia. The GCase activity was lower than the reference value. The pathological results showed that the proportion of hematopoietic tissue and fat in the bone and bone marrow was essentially normal. The proportion of granulocytes and erythrocytes in the hematopoietic tissue was approximately normal, and megakaryocytes were occasionally seen. The bone density was generally decreased in the lower and middle femur, knee, pelvis, and both hips. Hyperosteogeny was seen in the knee and both hips. There were no other skeletal abnormalities, fractures, pain and bruising. The physical examination results of the patient are shown in Table 1.
Identification of GBA1 mutations
Guangzhou Women's and Children's Hospital (Guangzhou, China) had previously performed sanger sequencing for 11 exons and the exon-intron junction region of the GBA1 gene. DNA sequencing revealed a mutation c.1240G > C (p.Val414Leu) and a known pathogenic mutation c.1342G > C (p. Asp448His) (Figure 2A). The mutations were predicted to alter both the secondary structure and the three-dimensional structures of GCase. In addition, the three-dimensional structures of the WT and MT proteins were visually distinct in both the overall and partial views. The c.1240G > C mutation might affect the catalytic function by disturbing the glycosylation process, and the substitution of the bulky phenylalanine with a small leucine residue might bring in a larger side chain that could form steric effects in the catalytic region (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1092645/fped-11-1092645-HTML/image_m/fped-11-1092645-g002.jpg)
Figure 2. Sequencing chromatogram and bioinformatics analysis. (A) Sanger sequencing results of the c.1240G > C (p.Val414Leu) and c.1342G > C (p.Asp448His) mutations. (B) Protein structure predicted by I-TASSER.
Evolutionary conservation analysis of multiple sequence alignments of the GCase protein and its homologs showed that the amino acid at residues 414 and 448 are evolutionarily conserved (Figure 3A), suggesting that these amino acids may play an important role in the function of GCase. We verified the constructed plasmids with sequencing, and the results are shown in Figure 3B.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1092645/fped-11-1092645-HTML/image_m/fped-11-1092645-g003.jpg)
Figure 3. (A) Evolutionary conservation analysis of a multiple sequence alignment of the GCase. The comparison between different eukaryotic species, showing that the base pair involving the nucleotide at positions 1,240 and 1,342 in GCase are highly conserved. (B) Sequencing validation of the constructed plasmids.
Functional analysis after plasmid transfection
To further study the functional consequences of both GBA1 mutation in vitro, we examined the expression of GBA1 using qPCR and Western blotting. Quantitative detection revealed that the WT GBA1 was more highly expressed than the MT at both the mRNA and protein levels (p < 0.05) (Figures 4A–C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1092645/fped-11-1092645-HTML/image_m/fped-11-1092645-g004.jpg)
Figure 4. Functional analysis after plasmid transfection. (A) mRNA expression level of GBA1 in HEK293T cells. There was a significant difference in the mRNA level of the wild-type and the mutant. (B) Western blotting analysis of GCase expression. (C) Protein expression level of GCase in HEK293T cells. Quantitative detection found that the WT GBA1 was more highly expressed than the MT at both the mRNA and protein levels (p < 0.05). (D) The subcellular localization analysis showed that both the MT and WT were localized at the cytoplasm. (E) The GCase concentration was significantly different between the WT and the V414L mutation or V414L + D448H. (F) The GCase enzyme activity assay was significantly different between the WT and the V414L mutation or V414L + D448H.
The subcellular localization analysis showed that both the MT and WT were localized at the cytoplasm, suggesting there is no difference in subcellular localization (Figure 4D). We also detected the GCase concentration, and the results showed that there was a significant difference between the WT and the MT with the V414L mutation, or that with V414L + D448H. However, no significant difference was found between the WT and D448H (Figure 4E). The results of GCase enzyme activity assay showed that there also was a significant difference between the WT and the MT with the V414L mutation, or that with the V414L + D448H mutation (Figure 4F).
Discussion
GD is a rare hereditary lysosomal storage disorder caused by defects in the GBA1 gene (7). Splenomegaly, bone disease, and thrombocytopenia are the most common characteristics of GD (8). The GBA1 gene is located in a gene-rich region on chromosome 1q21 and has a highly homologous pseudogene GBAP1 sequence located 16 kb downstream that shares 96% exonic homology with the functional gene. It is composed of 11 exons. A useful distinction for molecular diagnostic applications is that the pseudogene GBAP1 has a 55 bp deletion in exon 9. By Sanger sequencing of the exon 8–11 region of DNA from the patient's peripheral blood, we not only confirmed those two heterozygous mutations but also confirmed that both mutations are located in GBA1 rather than GBAP1. Many mutations associated with GD have been reported (4, 9). There are different types of GBA1 gene mutations and different pathogenic modes, so it is necessary to use appropriate methods to evaluate the pathogenicity of GBA1 gene mutations (3). Large fragment deletion-insertion mutations can directly affect the process of transcription and translation and are often considered pathogenic. In nonsense mutations, termination codons appear early, potentially leading to premature mRNA degradation or producing truncated proteins that affect its stability and normal function, so these are often considered pathogenic. Missense mutations are more difficult to determine and must be evaluated using biological information and functional studies.
In this study, we identified a heterozygous c.1240G > C (p.Val414Leu) mutation, on which functional studies, including subcellular localization, mRNA level, and protein expression, were performed to confirm its pathogenicity. This mutation has not been previously reported, and the possibility of polymorphism was excluded. This c.1240G > C mutation is predicted to be possibly damaging, and the amino acids are highly conserved in different species, which supports that this mutation is likely to be pathogenic. Another heterozygous pathogenic c.1342G > C (p.Asp448His) mutation, also known as D409H, contributes to α-synuclein accumulation and leads to GCase deficiency. This has been widely reported in at least 33 individuals (10). The allele frequency of this mutation in a cohort of 436 GD patients in Portugal and Spain was 3.3% (11). This variant typically causes a severe form of the disease, and could cause different severity of GD (3).
In vitro functional studies of GBA1 gene mutations can help us understand the molecular pathogenesis of these mutations and evaluate the pathogenicity of mutations. We found that the decrease in mRNA caused by compound heterozygous mutation was between two simple independent mutations. However, the decrease in protein level was not as obvious as that of the two simple independent mutations. This inconsistency may be due to protein stability and protein-protein interactions. There was no difference in subcellular localization, and the protein seemed to be located more in the nucleus after co-transformation, which may be due to protein degradation slowing down after co-transformation, resulting in slightly more protein. In this study, all results were consistent. At the molecular function level, the V414L mutation was more significantly altered than the V414L + D448H mutaition. Clinically, we observed that people carrying only the V414L heterozygous mutation did not have a severe phenotype. This may be because human have two chromosomes, one from the mother and one from the father. When one of them is mutated, the other can still normally translate proteins to meet body needs. However, when both chromosomes are mutated and the body is unable to translate completely the correct protein, it can lead to severe disease. From the molecular function results presented in this study, we can predict that patients will present with a more distinct disease phenotype when they carry the V414L homozygous mutation.
Genetic analysis can be used to confirm the diagnosis of GD, to help analyze the relationship between clinical phenotype and genotype, to evaluate the prognosis of the disease, and to aid in genetic counseling and prenatal diagnosis. Early diagnosis of GD can be achieved by combining clinical manifestations, laboratory testing, and genetic analyses (12, 13).
Our research provides a reference for the same mutant loci discovered in the future and to avoid duplication of work. On the other hand, With the continuous application of new technology and the deepening of the research, new mutant loci are continuously discovered, which may form a network of GBA1 gene mutations associated with Gaucher disease and provide new perspectives for further research on the pathogenesis of Gaucher disease.
Multiple studies have implicated heterozygous GBA1 mutations as a major genetic risk factor for Parkinson's disease (14, 15). Therefore, it is necessary to identify the diagnosis with a combination of clinical and biochemical examinations.
Conclusion
The mutation spectrum of GBA1 exhibits ethnic and regional disparity in Asian patients. Our study provides evidence that V414L and D448H in GBA1 gene can be pathogenic and lead to GD. The findings expand the phenotypic and mutation spectra, and will benefit further genetic studies of the GBA1 gene in patients with GD.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
QL: conceptualization; methodology; formal analysis. ZS: performed experiments; writing—original draft. HP: software, investigation, funding acquisition. SM: data curation. FX: visualization, supervision, project administration. FH: validation; writing—review and editing; project administration. All authors contributed to the article and approved the submitted version.
Funding
The work was supported by the National Key S&T Special Projects (grant no. 2021YFC100530), Department of Health of Zhejiang Province (grant no. 2021KY1137), National Natural Science Foundation of China (grant nos. 32170617, 31970558), Natural Science Foundation of Guangdong Province of China (grant no. 2022A1515012621).
Acknowledgments
The authors thank the patients and control subjects for their participation in this study.
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.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Horowitz M, Braunstein H, Zimran A, Revel-Vilk S, Goker-Alpan O. Lysosomal functions and dysfunctions: molecular and cellular mechanisms underlying Gaucher disease and its association with Parkinson disease. Adv Drug Deliv Rev. (2022) 187:114402. doi: 10.1016/j.addr.2022.114402
2. Arevalo NB, Lamaizon CM, Cavieres VA, Burgos PV, Álvarez AR, Yañez MJ, et al. Neuronopathic Gaucher disease: beyond lysosomal dysfunction. Front Mol Neurosci. (2022) 15:934820. doi: 10.3389/fnmol.2022.934820
3. Dimitriou E, Moraitou M, Cozar M, Serra-Vinardell J, Vilageliu L, Grinberg D, et al. Gaucher disease: biochemical and molecular findings in 141 patients diagnosed in Greece. Mol Genet Metab Rep. (2020) 24:100614. doi: 10.1016/j.ymgmr.2020.100614
4. Guo JN, Guan M, Jiang N, Li N, Li YJ, Zhang J, et al. Establishment and phenotypic analysis of the novel Gaucher disease mouse model with the partially humanized Gba1 gene and F213I mutation. Front Genet. (2022) 13:892457. doi: 10.3389/fgene.2022.892457
5. Phetthong T, Tim-Aroon T, Khongkraparn A, Noojarern S, Kuptanon C, Wichajarn K, et al. Gaucher disease: clinical phenotypes and refining GBA mutational spectrum in Thai patients. Orphanet J Rare Dis. (2021) 16:519. doi: 10.1186/s13023-021-02151-2
6. Huang Y, Jia X, Tang C, Liu S, Sheng H, Zhao X, et al. High risk screening for Gaucher disease in patients with splenomegaly and/or thrombocytopenia in China: 55 cases identified. Clin Chim Acta. (2020) 506:22–7. doi: 10.1016/j.cca.2020.03.016
7. Roh J, Subramanian S, Weinreb NJ, Kartha RV. Gaucher disease - more than just a rare lipid storage disease. J Mol Med. (2022) 100:499–518. doi: 10.1007/s00109-021-02174-z
8. Lepe-Balsalobre E, Santotoribio JD, Nunez-Vazquez R, García-Morillo S, Jiménez-Arriscado P, Hernández-Arévalo P, et al. Genotype/phenotype relationship in Gaucher disease patients. Novel mutation in glucocerebrosidase gene. Clin Chem Lab Med. (2020) 58:2017–24. doi: 10.1515/cclm-2020-0306
9. Garcia RS, de Frutos LL, Arreguin EA, González CC, Ortiz JEG, Ornelas SF, et al. Gaucher disease: identification and novel variants in Mexican and Spanish patients. Arch Med Res. (2021) 52:731–7. doi: 10.1016/j.arcmed.2021.05.001
10. Kim D, Hwang H, Choi S, Kwon SH, Lee S, Park JH, et al. D409h GBA1 mutation accelerates the progression of pathology in A53T alpha-synuclein transgenic mouse model. Acta Neuropathol Commun. (2018) 6:32. doi: 10.1186/s40478-018-0538-9
11. Liou B, Kazimierczuk A, Zhang M, Scott CR, Hegde RS, Grabowski GA. Analyses of variant acid beta-glucosidases: effects of Gaucher disease mutations. J Biol Chem. (2006) 281:4242–53. doi: 10.1074/jbc.M511110200
12. Kilavuz S, Basaranoglu M, Epcacan S, Bako D, Ozer A, Donmez YN, et al. A rare cause of hydrops fetalis in two Gaucher disease type 2 patients with a novel mutation. Metab Brain Dis. (2022) 37:1283–7. doi: 10.1007/s11011-022-00942-5
13. Erdal I, Yildiz Y, Onal G, Aktepe OH, Düzgün SA, Sağlam A, et al. Splenic gaucheroma leading to incidental diagnosis of Gaucher disease in a 46-year-old man with a rare GBA mutation: a case report. Endocr Metab Immune Disord Drug Targets. (2022) 23:230–4. doi: 10.2174/1871530322666220420134935
14. Toffoli M, Chen X, Sedlazeck FJ, Lee C-Y, Mullin S, Higgins A, et al. Comprehensive short and long read sequencing analysis for the Gaucher and Parkinson’s disease-associated GBA gene. Commun Biol. (2022) 5:670. doi: 10.1038/s42003-022-03610-7
Keywords: Gaucher disease, GBA1, c.1240G > C, c.1342G > C, cell experiment
Citation: Liu Q, Shen Z, Pan H, Ma S, Xiong F and He F (2023) The molecular mechanism of Gaucher disease caused by compound heterozygous mutations in GBA1 gene. Front. Pediatr. 11:1092645. doi: 10.3389/fped.2023.1092645
Received: 8 November 2022; Accepted: 6 January 2023;
Published: 26 January 2023.
Edited by:
Nahid Tayebi, Hematopoiesis Section, Genetics and Molecular Biology Branch (NHGRI), United StatesReviewed by:
Khalid Al-Thihli, Sultan Qaboos University, OmanNicolina Cristina Sorrentino, Telethon Institute of Genetics and Medicine (TIGEM), Italy
© 2023 Liu, Shen, Pan, Ma, Xiong and He. 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: Fu Xiong xiongfu@smu.edu.cn Fei He 15218809216@163.com
†These authors have contributed equally to this work
Specialty Section: This article was submitted to Genetics of Common and Rare Diseases, a section of the journal Frontiers in Pediatrics