Skip to main content

Gaucher disease: clinical phenotypes and refining GBA mutational spectrum in Thai patients



Gaucher disease (GD) is a rare lysosomal storage disorder, characterized by hepatosplenomegaly and pancytopenia, with or without neurologic involvement. The disorder is categorized into three phenotypes: GD type 1 or nonneuronopathic GD; GD type 2 or acute neuronopathic GD; and GD type 3 or chronic neuronopathic GD. The purposes of this study were to describe clinical characteristics of Thai GD in patients diagnosed and/or followed up during 2010–2018 and to perform re-genotyping including analysis of GBA recombinant alleles which had not been investigated in Thai patients before.


There were 27 patients from seven medical centers, enrolled in the study. All the cases had pediatric onset. GD3 (44.5%) was the most common phenotype, followed by GD2 (40.7%) and GD1 (14.8%), with one case of neonatal GD. The median age of onset for GD1, GD2, and GD3 was 72, 4 and 12 months, respectively, suggesting relatively earlier onset of GD1 and GD3 in Thai patients. All patients with GD1 and most patients with GD3 received ERT. Four patients with GD3 had ERT followed by HSCT. Patients with GD3 who received no or late ERT showed unfavorable outcomes. We identified 14 variants including two novel (p.S384F and p.W533*) and 12 reported pathogenic variants: p.L483P, p.N409S, p.R159W, p.P305A, p.A175G, p.D448H, p.V414L, IVS2+1G>A, IVS6-1G>C, IVS7+1G>C, IVS9-3C>G, and Rec1a. The p.L483P was the most prevalent allele found in this study, at 66% (33/50 alleles), followed by IVS2+1G>A, Rec1a, and IVS6-1G>C. Twenty-four percent of patients were reassigned with validated genotypes, most of whom (4 of 6) were patients with GD2. The [p.S384F + p.W533*] being compounded with p.L483P, was found in the patient with neonatal GD, suggesting that the p.S384F could potentiate the deleterious effect of the p.W533*, and/or vice versa.


Neuronopathic GD was strikingly prevalent among Thai affected population. Homozygous p.L483P was the most common genotype identified in Thai patients. Recombinant allele Rec1a and splicing mutations were associated with GD2 and severe cases of GD3. Mutation spectrum could be useful for designing stepwise molecular analysis, genetic screenings in population, and new therapeutic research for neuronopathic GD.


Gaucher disease (GD) is a lysosomal storage disorder, caused by deficiency of glucocerebrosidase enzyme, resulting in multi-systemic manifestations including hepatosplenomegaly, pancytopenia, with or without neurologic involvement. The disorder is categorized into three clinical phenotypes [1]. Non-neuronopathic or GD type 1 (GD1; OMIM 230800) has predominant hematologic and visceral organ involvement with additional bone symptoms in some patients. Acute neuronopathic GD type 2 (GD2; OMIM 230900) has onset before age two years and rapid progression leading to death before age four years. Subacute neuronopathic GD type 3 (GD3; OMIM 231000) may have onset before or after two years of age, with later onset and slower progression of neurological symptoms. GD3 can be further divided into 3a: mild visceral symptoms but progressive neurologic decline; 3b: severe visceral involvement; 3c: fibrosis and calcification of aorta and cardiac valves; and Swedish L444P-related Norrbottnian phenotype which is characterized by progressive kyphosis and early onset visceral and neurological manifestations [2]. There is also a very rare distinct subtype, perinatal-lethal form associated with severe hepatosplenomegaly and ichthyosiform, nonimmune hydrops, and early death [1]. Moreover, there was a reported patient of an unclassified subtype of Gaucher disease, who presented with a neonatal-onset of hepatosplenomegaly and pancytopenia, without profound neurological disease [3].

GD1 accounts for 95% of GD found among Caucasian population while GD2 and GD3 are more prevalent among Asian affected populations [1, 4,5,6,7].

The glucocerebrosidase enzyme is encoded by acid-beta glucosidase (GBA) gene located on chromosome 1q21 and spanning 7 kb with 11 coding exons. GBAP, a non-functional pseudogene of GBA, shares 96% homology of the overall sequence and is located 16 kb downstream [8]. The presence of GBAP complicates mutation analysis of GBA and can lead to recombinant allele which is much more difficult to detect [6,7,8,9].

There have been > 449 damaging GBA variants identified, including 357 missense/nonsense/splicing variants and 23 complex rearrangement (the Human Gene Mutation Database, Most GBA mutations resulting in GD are point mutations with majority clustering within exons 8–11. N370S (p.Asn409Ser) and L444P (p.Leu483Pro) are the most prevalent variants described among Jewish and non-Jewish affected populations, respectively.

The complex rearrangement is resulted from homologous recombination between the GBA and GBAP. It could be either non-reciprocal (gene conversion) or reciprocal recombination (gene fusion and partial gene duplication). These alleles are likely underrepresented because it is difficult to analyze and often left out from the routine analysis due to its complexity. It has been shown that recombination accounts for up to 24% in some affected populations, especially in neuronopathic GD [5, 7, 10] and low at 3% in some others [11, 12].

Among GD with complex mutations, gene conversion represents 58% of the recombination. A specific allele, namely Rec1a, is the most prevalent (91%) variant amongst the gene conversion alleles identified [7, 8]. Rec1a has been previously reported in Asian GD [4, 7, 8, 13, 14] but never been detected in Thai patients [15,16,17]. The possible reason that Rec1a had not been reported in Thai patients, likely because the complex rearrangement/conversion was not detected by the methods used in those studies.

To date, reports of GD with mutation data from Thailand are limited. We previously described four patients with neuronopathic GD and their molecular defects. Herein, we present 27 unreported patients and their genotypes including the first study of recombinant alleles in Thai GD.



This is a cross-sectional and descriptive study of seven major teaching hospitals: Ramathibodi Hospital, Queen Sirikit National Institute of Child Health, Srinagarind Hospital (Khon Kaen University), Siriraj Hospital, King Chulalongkorn Memorial Hospital, Chiang Mai University Hospital, and Surin Hospital. These institutes cover most of the patients diagnosed with GD in the country. Patients with GD who were newly diagnosed and/or followed during 2010 and 2018 (9 years) were enrolled in the study.

Clinical data including age at onset, age at diagnosis, clinical manifestations, type of GD, enzyme activity, specific treatment received including enzyme replacement therapy (ERT) and hematopoietic stem cell transplantation (HSCT) were retrospectively reviewed. Type of GD was categorized by the treating physicians, according to clinical ground, age of onset, and the presence of neurological involvement. The neurological manifestations referred to bulbar signs (stridor, strabismus, feeding difficulty), pyramidal signs (opisthotonus, spasticity, trismus) and oculomotor apraxia (failure of saccadic eye movement) [1].

GD1 is assigned for patients with the absence of neurological symptoms related to GD throughout the time course. The presence of spasticity, seizure, and/or developmental regression before 2 years of age was assigned as GD2. Patients with any neurological manifestation and/or the ocular signs other than those of GD2 would be categorized as GD3.

Re-genotyping of the GBA gene and analysis of recombinant alleles

Patient’s mutation data prior to this study were collected, and later compared with the result from the re-genotyping. Archived/left over specimens were sent to Ramathibodi Hospital for re-genotyping. We used amino acid numbering system following Human Genome Variation Society ( nomenclature with preceding p., and also provided traditional number without the proceeding p., for example p.L483P and L444P representing the same amino acid variant.

To avoid amplification of non-functional pseudogene, GBAP, a nested PCR approach was employed to selectively amplify the GBA gene, following established protocols [6, 9]. In brief, the first amplification (10 cycles) using specific primers sets (GBA-f and MTX1P-r) resulted in DNA fragment containing the entire GBA gene and part of MTX1P (pseudo-MTX1 gene). This was followed by nested PCR (35 cycles) to amplify the entire GBA gene, excluding the MTX1P, using different primer pair (GBA-f and GBA-r; Additional file 1: Table S1 and Figure S1A). Both PCR reactions shared similar constituents as follow: 25 ul of total volume, 250 ng of genomic DNA, 2.5 µM of each primer, 1 × PCR buffer, 2 mM MgCl2, and 2.5 U of Taq polymerase. (TaKaRa LA Taq™). The PCR cycling conditions were 95 °C, 1 min; then 10 (or 35) cycles of denaturation at 95 °C, 30 s and annealing at 60 °C, 30 s with extension at 72 °C, 5 min; then final extension at 68 °C for 10 min.

The expected final 6.5 kb-PCR product was then subjected to Sanger sequencing of all 11 GBA exons, using ABI 3730xl DNA Analyzer (Macrogen®, Korea). Details of sequencing primers were provided in the Supplementary data (Additional file 1: Figure S1B and Table S1). Gene conversion allele could be detected by careful examination of the sequence, focusing on every site known to be different between the GBA and GBAP (Additional file 1: Figure S2). The presence of three or more point mutations in cis, distinctive sequence of GBAP, and/or detection of the cross over site between GBA and GBAP (as a result of GBAP being converted into the GBA gene) would indicate gene conversion. The assignments of gene conversion allele were as follows: Rec1a − [L444P + A456P + V460V]; Rec2a −  [D409H + L444P + A456P + V460V]; Rec3a −  [55-bp deletion in exon 9 + D409H + L444P + A456P + V460V], as previously established [6,7,8,9, 18]. As for Rec4a, Rec5a, and Rec8a, the final PCR product’s size would be several-hundred-bp-smaller than 6.5 kb which could be seen on the agarose gel, in addition to the [55-bp deletion in exon 9 + D409H + L444P + A456P + V460V] plus the distinctive variants of GBAP (Additional file 1: Figure S2). Rec6a and Rec7a alleles could not be detected by the method used in this study.

The GenBank reference sequences used were ENST00000368373.8 and NM_000157.4. Variant pathogenicity was determined in accordance with the guidelines of the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) [19]. Rec1a positive control was kindly given by Prof. Yin-Hsiu Chien, National Taiwan University Hospital.


Patients’ demographic and clinical data

There were 27 patients (12 male and 15 female) enrolled in the study, including 4 GD1, 10 GD2,12 GD3, and 1 neonatal GD. To the best of our knowledge, this report covers about 80% of GD patients in Thailand during the time of study. Eleven of 27 (40.7%) patients were recruited from Ramathibodi Hospital and the remainders (59.3%) were enrolled via the other six participating hospitals (3 in the Central region, 1 in the Northern and 2 in the Northereastern area). Of 27 patients, 10 originated from the Central region, 10 from the Northern area, and the other 7 from Northeastern region.

The age at onset ranged from 2 weeks (wk) to 10 years (y), and age at diagnosis ranged between 4 months (m) and 26 y 11 m (Tables 1, 2). The median and mean age at onset, at diagnosis, and time to diagnosis for each clinical type was as detailed in Table 1.

Table 1 Demographic data of Gaucher disease patients in the present study (n = 27)
Table 2 Clinical characteristics and specific treatment received of each patient (n = 27)

Abdominal distension is the most consistent finding in all patients (100%), followed by anemia (59.3% or 16/27), thrombocytopenia (55.6% or 15/27), and neurological symptoms including developmental delay/regression with/without seizure disorders (51.9% or 14/27). Oculomotor apraxia (supranuclear gaze palsy) and/or esotropia was noted at the first presentation or during the time course in 44.4% (12/27).

All (4) patients with GD1 received ERT and are still alive. Ten of 12 patients with GD3 received ERT. Four of the 10 patients subsequently underwent HSCT and had ERT discontinued soon after and they were alive at the time of this writing. One patient with GD3 (GD-26) had ERT initiated at 11 years of age, she later developed seizures and neurological deterioration, leading to death at the age of 13 years. Two patients with GD3 were not treated with ERT due to unavailability of the treatment at that time (year 2010). As for GD2, all (10) patients did not receive ERT and died at the age 10 months–2 years 6 months. One patient with neonatal-typed GD was attempted with ERT for 6 months which led to an improvement of hematologic and visceral symptoms, but ongoing neurological decline and eventually death at the age of 2 years. There was no apparent correlation between the residual enzyme activity and type of GD.

As for why some (4) patients underwent HSCT, there were multiple reasons for this. The four patients had GD3 and genotype of homozygous p.L483P, at risk for progressive neurological deterioration. Before the year 2013, ERT was not reimbursable in the country, two patients (GD-4 and GD-6) had received ERT under a compassionate use program with unclear future of a long-term supply. Therefore, HSCT was offered to the patients, following 1–2 years of ERT. After 2013, additional families (GD-7 and GD-8) were interested in having HSCT for a cure instead of life-long ERT and for that it might be helpful in mitigating the progression of neurological manifestations. One patient (GD-4) received HLA-matched unrelated donor transplantaion and the other three patients (GD-6, GD-7, and GD-8) received HLA-haploidentical transplantation. The details of patients’ clinical and laboratory data including outcomes are being prepared for a separate manuscript. It should be mentioned that HSCT for severe thalassemia, a common genetic disorder in Thailand, has been widely accepted in our country which may influence the decision of the families to accept HSCT for GD3.

GBA mutational spectrum and recombinant allele, Rec1a

Prior to the study, mutation data of the enrolled patients were determined in four different laboratories including three in the country and one abroad. These tests did not cover analysis for recombinant alleles.

The re-genotyping in this study identified 14 variants including two novel (p.S384F and p.W533*) and 12 reported pathogenic variants: p.L483P, p.N409S, p.R159W, p.P305A, p.A175G, p.D448H, p.V414L, IVS2+1G>A, IVS6-1G>C, IVS7+1G>C, IVS9-3C>G, and Rec1a (RecNciI) (Table 3, Fig. 1 and Additional file 1: Figure S3).

Table 3 Genotype according to phenotype of each GD patient (n = 27)
Fig. 1
figure 1

Prevalence of Gaucher disease phenotypes and GBA mutant allele in Thai patients. A Data of the present cohort (n = 27 for phenotype analysis; 25 individuals of 50 alleles for genotype analysis). B Data combined with three previous reports (total n = 36 for phenotype analysis)

The p.S384F was classified as likely pathogenic (PM1, PM2, PP2, PP5) and the p.W533* as pathogenic (PVS1, PM2, PP3), according to the ACMG/AMP guideline. These two variants were not found in the gnomAD database (MAF = 0) and Thai Reference Exome Database (MAF 0/2000; T-REx, Multiple lines of computational evidence (DEOGEN2, FATHMM-MKL, M-CAP and MutationTaster) suggested pathologic impact of the variants. The p.S384F and p.W533* were submitted to ClinVar Database and were assigned references#, VCV000917863 and VCV000917862, respectively ( Interestingly, the two variants were present in cis [p.S384F + p.W533*] and genetically compounded with p.L483P in a single patient with neonatal GD (GD-28; Table 3). The [p.S384F + p.W533*] allele was inherited from the patient’s asymptomatic mother and the p.L483P from the father.

The prevalence of variants identified according to phenotypes was shown in Fig. 1A and Additional file 1: Table S2. Genotype of each patient was summarized in Table 3. Allele frequency of the GBA variants in affected populations was compared between this study and other studies (Table 4).

Table 4 Frequencies of GBA pathogenic variants in Thai Gaucher disease and other ethnicities


To our knowledge, this the first multicenter study and analysis of GBA recombinant allele in Thai patients, representing the landscape of GD in Thailand in recent years. In the present cohort, all the cases had pediatric onset. GD3 (44.5%) was the most common phenotype in the Thai population, followed by GD2 (40.7%) and GD1 (14.8%), making neuronopathic GD accounting for the majority of cases. When combining our data to previous reports with available information on types of GD and mutations, it still supported that GD3 was the most prevalent phenotype (41.7%) in Thai patients, followed by GD2 (33.3%), and higher proportion of GD1 (25.0%) (Fig. 1B and Additional file 1: Table S2) [15,16,17].

Notably, our patient GD-28 had neonatal-onset hepatosplenomegaly, anemia, and thrombocytopenia and skin bruises. There was no evidence of neurological disease indicating GD2 or any ichthyosiform or collodion skin abnormalities found in perinatal-lethal form GD. We performed literature serach and found only one patient who was described with homozygous p.D448H and the phenotypes similar to our patient [3].

In the present cohort, 85.2% (23/27) of the patients developed first symptoms at age ≤ 5 years and 74.1% (20/27) at age ≤ 1 year of age, indicating relatively early and severe phenotype of Thai GD as compared to those of Western population. All GD3 patients manifested their first symptoms at age ≤ 3 years and two-thirds at age ≤ 1 year. Nine of 12 GD3 patients had subtype 3b, while the others had phenotypes overlapping between 3a and 3b. Patients with GD3 who did not receive ERT died at very young age while those with comparable onset who had ERT were still living. Early diagnosis and timely initiation of ERT appeared to result in clinical improvement and slowing down and/or preventing the progression of the neurological manifestations. Delayed diagnosis and treatment with ERT led to inevitable neurological progression and death, as seen in one of our patients (GD-26).

The median age of onset of GD1 and GD3 was at 72 and 12 months, respectively, indicating earlier age of onset as compared to Western affected population. As for, GD2, the onset of symptom was at 4.0 months and natural course did not seem to be different from those described in the literature [1]. The average time to diagnosis was shortest in GD2, at 5.4 months, and longer in GD3 at 19.6 months., and GD1 at 45.3 months. This likely reflected the nature of rapid disease progression of GD2, demanding for immediate investigation.

The p.L483P was the most prevalent allele found in this study, at 66% (33/50 alleles), followed by IVS2+1G>A (8%), Rec1a (4%), and IVS6-1G>C (4%) (Table 2). These four alleles together accounted for 82% (41/50) of the mutant alleles identified in this cohort or 75% (51/68) of total Thai GD alleles (Additional file 1: Figure S2) [15,16,17].

Rec1a allele was found in two patients with neuronopathic GD in the present study, supporting its association with severe disease. Rec1a has been described to be associated with poor clinical outcomes and sometimes mistaken as p.L483P (L444P) [7, 8]. Previous studies have shown that a significant number of GD patients were initially assigned with point mutation(s), but in fact they carried recombinant allele(s) [6,7,8,9]. All other variants were found, each in one patient. The alleles associated with GD2 in the present study included p.R159W, p.D488H, p.L483P, Rec1a, IVS2+1G>A, IVS6-1G>C, IVS7+1G>A, and IVS9-3C>G.

In this cohort, the most common genotype was homozygous p.L483P, found in 11 of 12 individuals with GD3, suggesting more homogeneous population of GD3 in Thai affected individuals. Of the 11 patients with homozygous p.L483P, 5 were from the Northern area, 4 from the Northeastern area, and 2 from the Central region, implying a nationwide distribution of this genotype. The second common genotype was genetic compound between p.L483P and other allele which was found in 11 patients: 1 GD3, 6 GD2, 1 neonatal GD, and 3 GD1. We observed that genetic compounds between p.L483P and the recombinant allele, splicing variants, or known severe missense alleles (p.D488H and p.R159W) would lead to neuronopathic GD while the genetic compounds between p.L483P and p.N409S or p.P305A would result in GD1 with slow clinical progression, as described previously [20]

Twenty-four percent (6 of 25 studied) of patients in this cohort were reassigned with validated genotypes, most of whom (4 of 6) were patients with GD2. Therefore, we suggest inclusion of recombinant allele analysis for Thai GD patients especially those with suspected GD2 and GD3 with rapid progression. Plausible explanations for the discrepancy between the results from the re-genotyping and the prior mutations were using different detailed method of mutation screening, such as using vs not using long ranged PCR, excluding vs not excluding the pseudogene, and limited mutation screening (only for common mutations). Our data support previous observation that the detection of recombinant allele and correct variants of GBA depends on the methodology of genetic analysis [12]. The correct identification of mutant alleles could lead to accurate genetic counseling, reproductive planning, and therapeutic decision making.

As for why mutations in some patients (GD-11, GD-17, and GD-27) were drastically different from its initial reports, our guess was that the orginal sequencing was just incorrect, as the method used for the GBA analysis did not excluded the GBA pseudogene (GBAP). The IVS2+1G>A is an allele present in the GBAP (Additional file 1: Figure S2), while it could be a pathogenic mutation present in the GBA gene of some patients. The p.R159W is a known pathogenic variant associated with severe phenotypes of GD. The p.A175G and p.L483P represent common alleles in GBAP (Additional file 1: Figure S2) whereas the p.P305A and c.IVS6-1G>C are not of the GBAP and have been described as pathogenic variants in GD.

Given the [p.S384F + p.W533*] being compounded with p.L483P, resulted in an extreme phenotype of neonatal GD, we hypothesized that the p.S384F could potentiate the deleterious effect of the nonsense mutation, p.W533*, and/or vice versa. Further investigations are required to confirm the pathogenic mechanism of these intriguing alleles. The presence of two alleles in cis has been shown associated with a more severe phenotype in GD, homozygous D409H leading to GD3 with cardiac phenotype in most cases whilst homozygous [D409H + H255Q] resulting in GD2 [21].

Our result suggested that mutation spectrum of GBA in Thai patients was more homogeneous than ever believed and was different from those seen in Caucasian populations, and more similar to those of East Asian descendants [4, 6, 11, 13, 15,16,17, 22].

The older report of GD in Thailand, during 1966–1998, consisting of 25 GD patients diagnosed over 33 year-period or 0.8 case per year [23], When combining our data with other recent reports of Thai GD, a total number of patients was 33 (27 + 5 + 1) over the 9 years of study (2010–2018) or 3.7 new cases per year [15, 17]. These data could reflect the better awareness of GD among Thai physicians, resource availability including clinical geneticists, diagnostic facilities, and reimbursement for ERT (since 2013) and HSCT (since 2015), which have been escalating lately. We estimate the prevalence of GD in the Thai population at least 1 in 156,818 and carrier frequency of 1 in 198 as ascertained by using Hardy–Weinberg equation and 6,449,073 live births during the study years, and based on the assupmption that our study population accounts for 80% of total patients across the country ( We believe that these numbers could be rather underestimated due to some hospitals not included in the study and patients escaping diagnosis especially those with later onset and mild GD (possible GD1).

As it is known that ERT does not cross blood brain barrier and that patients with GD3 could exhibit more advanced neurological manifestation some time in their life; therefore, new therapeutic approach and/or adjunct treatment targeting to protect neurological involvement may be required to overcome these challenges. Long-term outcome of GD3 with quite homogeneous genotype of homozygous p.L483P deserves further investigation..

There are a number of limitations of this work including missing data due to the natures of the retrospective study and varying experiences of clinicians in detecting specific subtle signs (i.e. oculomotor apraxia) which may perturb the accuracy of clinical data collected and the assignment of GD phenotype in particular GD1 and GD3. Lastly, analysis for reciprocal recombinant alleles and some gene conversion alleles was not included in this study, which might affect the variant distribution, though the effect was believed to be small due to its rarity.


GD3 followed by GD2 was strikingly prevalent among the Thai affected population. Homozygous p.L483P was the most common genotype identified in Thai GD. Recombinant allele Rec1a and splicing mutations were associated with GD2 and severe cases of GD3. Mutation spectrum and its frequency could be useful for designing stepwise molecular analysis, genetic screenings in population, and new therapeutic research for neuronopathic GD.

Availability of data and materials

Additional data are provided in Supplementary materials and ClinVar Database (



Enzyme replacement therapy


Acid-beta glucosidase gene


Acid-beta glucosidase pseudogene gene


Gaucher disease


Gaucher disease type 1


Gaucher disease type 2


Gaucher disease type 3


Hematopoietic stem cell transplantation


  1. Pastores GM, Hughes DA. Gaucher Disease. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Mirzaa G, et al., editors. GeneReviews((R)). Seattle (WA). 1993.

  2. Germain DP. Gaucher’s disease: a paradigm for interventional genetics. Clin Genet. 2004;65(2):77–86.

    Article  CAS  Google Scholar 

  3. Roth P, Sklower Brooks S, Potaznik D, Cooma R, Sahdev S. Neonatal Gaucher disease presenting as persistent thrombocytopenia. J Perinatol. 2005;25(5):356–8.

    Article  Google Scholar 

  4. Eto Y, Ida H. Clinical and molecular characteristics of Japanese Gaucher disease. Neurochem Res. 1999;24(2):207–11.

    Article  CAS  Google Scholar 

  5. Ida H, Iwasawa K, Kawame H, Rennert OM, Maekawa K, Eto Y. Characteristics of gene mutations among 32 unrelated Japanese Gaucher disease patients: absence of the common Jewish 84GG and 1226G mutations. Hum Genet. 1995;95(6):717–20.

    Article  CAS  Google Scholar 

  6. Jeong SY, Park SJ, Kim HJ. Clinical and genetic characteristics of Korean patients with Gaucher disease. Blood Cells Mol Dis. 2011;46(1):11–4.

    Article  CAS  Google Scholar 

  7. Wan L, Hsu CM, Tsai CH, Lee CC, Hwu WL, Tsai FJ. Mutation analysis of Gaucher disease patients in Taiwan: high prevalence of the RecNciI and L444P mutations. Blood Cells Mol Dis. 2006;36(3):422–5.

    Article  CAS  Google Scholar 

  8. Tayebi N, Stubblefield BK, Park JK, Orvisky E, Walker JM, LaMarca ME, et al. Reciprocal and nonreciprocal recombination at the glucocerebrosidase gene region: implications for complexity in Gaucher disease. Am J Hum Genet. 2003;72(3):519–34.

    Article  CAS  Google Scholar 

  9. Jeong SY, Kim SJ, Yang JA, Hong JH, Lee SJ, Kim HJ. Identification of a novel recombinant mutation in Korean patients with Gaucher disease using a long-range PCR approach. J Hum Genet. 2011;56(6):469–71.

    Article  CAS  Google Scholar 

  10. Koprivica V, Stone DL, Park JK, Callahan M, Frisch A, Cohen IJ, et al. Analysis and classification of 304 mutant alleles in patients with type 1 and type 3 Gaucher disease. Am J Hum Genet. 2000;66(6):1777–86.

    Article  CAS  Google Scholar 

  11. Horowitz M, Pasmanik-Chor M, Borochowitz Z, Falik-Zaccai T, Heldmann K, Carmi R, et al. Prevalence of glucocerebrosidase mutations in the Israeli Ashkenazi Jewish population. Hum Mutat. 1998;12(4):240–4.

    Article  CAS  Google Scholar 

  12. Tajima A, Yokoi T, Ariga M, Ito T, Kaneshiro E, Eto Y, et al. Clinical and genetic study of Japanese patients with type 3 Gaucher disease. Mol Genet Metab. 2009;97(4):272–7.

    Article  CAS  Google Scholar 

  13. Feng Y, Huang Y, Tang C, Hu H, Zhao X, Sheng H, et al. Clinical and molecular characteristics of patients with Gaucher disease in Southern China. Blood Cells Mol Dis. 2018;68:30–4.

    Article  CAS  Google Scholar 

  14. Zhang W, Oehrle M, Prada CE, Schwartz IVD, Chutipongtanate S, Wattanasirichaigoon D, et al. A convenient approach to facilitate monitoring Gaucher disease progression and therapeutic response. Analyst. 2017;142(18):3380–7.

    Article  CAS  Google Scholar 

  15. Chavananon S, Sripornsawan P, Songthawee N, Chotsampancharoen T. Successful treatment of gaucher disease with matched sibling hematopoietic stem cell transplantation: a case report and literature review. J Pediatr Hematol Oncol. 2021;43(8):e1153–5.

    Article  Google Scholar 

  16. Suwannarat P, Keeratichamroen S, Wattanasirichaigoon D, Ngiwsara L, Cairns JR, Svasti J, et al. Molecular characterization of type 3 (neuronopathic) Gaucher disease in Thai patients. Blood Cells Mol Dis. 2007;39(3):348–52.

    Article  CAS  Google Scholar 

  17. Tammachote R, Tongkobpetch S, Srichomthong C, Phipatthanananti K, Pungkanon S, Wattanasirichaigoon D, et al. A common and two novel GBA mutations in Thai patients with Gaucher disease. J Hum Genet. 2013;58(9):594–9.

    Article  CAS  Google Scholar 

  18. Hruska KS, LaMarca ME, Scott CR, Sidransky E. Gaucher disease: mutation and polymorphism spectrum in the glucocerebrosidase gene (GBA). Hum Mutat. 2008;29(5):567–83.

    Article  CAS  Google Scholar 

  19. 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(5):405–24.

    Article  Google Scholar 

  20. Filocamo M, Mazzotti R, Stroppiano M, Seri M, Giona F, Parenti G, et al. Analysis of the glucocerebrosidase gene and mutation profile in 144 Italian gaucher patients. Hum Mutat. 2002;20(3):234–5.

    Article  Google Scholar 

  21. Michelakakis H, Moraitou M, Dimitriou E, Santamaria R, Sanchez G, Gort L, et al. Homozygosity for the double D409H+H255Q allele in type II Gaucher disease. J Inherit Metab Dis. 2006;29(4):591.

    Article  CAS  Google Scholar 

  22. Barney AM, Danda S, Abraham A, Fouzia NA, Gowdra A, Abraham SSC, et al. Clinicogenetic profile, treatment modalities, and mortality predictors of gaucher disease: a 15-year retrospective study. Public Health Genomics. 2021;2021:1–10.

    Google Scholar 

  23. Tanphaichitr VS, Suvatte V, Mahasandana C, Sachapong P, Veerakul G, Kankirawatana S, et al. Gaucher’s disease;thirty-two years experience at Siriraj Hospital. Southeast Asian J Trop Med Public Health. 1999;30(Suppl 2):143–7.

    PubMed  Google Scholar 

Download references


We would like to thank Mahidol University and Ramathibodi Foundation for partial funding of this project, Dr.Ponghatai Damrongphol for assistance in data collection, and Prof. Yin-Hsiu Chien, National Taiwan University Hospital for kind generosity in giving positive control of the recombinant allele.. DW and TT are recipients of the Research Career Development Awards from the Faculty of Medicine Ramathibodi Hospital, Mahidol University.


Mahidol University and Ramathibodi Foundation.

Author information

Authors and Affiliations



TP performed molecular analysis and prepared the manuscript draft. TT interpreted biochemical and molecular data. AK and SN performed and assisted in biochemical and/or molecular analysis. TP, TT, CK, KS, KW, PC, AS, AT, NN, and DW provided patient care, clinical information, biochemical and prior genotype data, and/or left over specimens for re-genotyping. DW designed the study concept, supervised the study and data interpretation, and critically revised the manuscript. All authors read and approved the final manuscript.

Author’s information

TP, TT, CK, KS, KW, AS and DW are pediatric geneticists. AK and SN are laboratory technologists. PC is a pediatric hematologist. AT is an adult hematologist. NN is a neonatologist.

Corresponding author

Correspondence to Duangrurdee Wattanasirichaigoon.

Ethics declarations

Ethics approval and consent to participate

The study was conducted following an approval by the Ramathibodi Hospital Human Research Ethics Committee (protocol ID 467; MURA2019/587). Consent to participate was not applicable because of retrospective study and using left over/archived specimens.

Consent for publication

Not applicable.

Competing interest


Additional information

Publisher's Note

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

Tim Phetthong started the study at Faculty of Medicine Ramathibodi Hospital and continued the work till completion at Phramongkutklao College of Medicine

Supplementary Information

Additional file 1.

Supplemental methods and results.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Phetthong, T., Tim-Aroon, T., Khongkraparn, A. et al. Gaucher disease: clinical phenotypes and refining GBA mutational spectrum in Thai patients. Orphanet J Rare Dis 16, 519 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: