- Case report
- Open Access
Novel manifestations of Warburg micro syndrome type 1 caused by a new splicing variant of RAB3GAP1: a case report
BMC Neurology volume 21, Article number: 180 (2021)
The present study aimed to determine the underlying genetic factors causing the possible Warburg micro syndrome (WARBM) phenotype in two Iranian patients.
A 5-year-old female and a 4.5-year-old male were referred due to microcephaly, global developmental delay, and dysmorphic features. After doing neuroimaging and clinical examinations, due to the heterogeneity of neurodevelopmental disorders, we subjected 7 family members to whole-exome sequencing. Three candidate variants were confirmed by Sanger sequencing and allele frequency of each variant was also determined in 300 healthy ethnically matched people using the tetra-primer amplification refractory mutation system-PCR and PCR-restriction fragment length polymorphism. To show the splicing effects, reverse transcription-PCR (RT-PCR) and RT-qPCR were performed, followed by Sanger sequencing. A novel homozygous variant—NM_012233.2: c.151-5 T > G; p.(Gly51IlefsTer15)—in the RAB3GAP1 gene was identified as the most likely disease-causing variant. RT-PCR/RT-qPCR showed that this variant can activate a cryptic site of splicing in intron 3, changing the splicing and gene expression processes. We also identified some novel manifestations in association with WARBM type 1 to touch upon abnormal philtrum, prominent antitragus, downturned corners of the mouth, malaligned teeth, scrotal hypoplasia, low anterior hairline, hypertrichosis of upper back, spastic diplegia to quadriplegia, and cerebral white matter signal changes.
Due to the common phenotypes between WARBMs and Martsolf syndrome (MIM: 212720), we suggest using the “RABopathies” term that can in turn cover a broad range of manifestations. This study can per se increase the genotype-phenotype spectrum of WARBM type 1.
Warburg Micro syndrome (WARBM)—also known as Micro Syndrome—is a rare autosomal recessive disorder with unknown true incidence. This syndrome is characterized by neurodevelopmental abnormalities such as congenital or postnatal microcephaly, severe intellectual disability, pachy- or polymicrogyria, and hypoplasia/agenesis of the corpus callosum in addition to the ocular manifestations including congenital cataract, microcornea, microphthalmia, and optic atrophy [1, 2]. This syndrome embraces four subtypes as in WARBM type 1 (MIM: 600118) which is caused by biallelic mutations of the RAB3 GTPase-activating protein 1 gene (RAB3GAP1; MIM: 602536) , WARBM type 2 (MIM: 614225) by mutations of RAB3GAP2 (MIM: 609275) , WARBM type 3 (MIM: 614222) by mutations of the Ras-related protein Rab-18 (RAB18; MIM: 602207) , and WARBM type 4 (MIM: 615663) by mutations in the TBC1 Domain Family Member 20 (TBC1D20; MIM: 611663) . Of these, RAB3GAP1 mutations are by far the most common and approximately account for 40% of WARBM cases .
RAB3GAP1 contains 24 exons and encodes the catalytic subunit of a Rab GTPase activating protein (GAP) that has specificity for the Rab3 subfamily (i.e. RAB3A, RAB3B, RAB3C, and RAB3D) . The catalytic activity of RAB3GAP1 is attributed to its C-terminus (i.e. codons 601–981) . Rab3 proteins are required for normal eye and brain development , proper exocytosis of neurotransmitters  and hormones , and also are imperative for other neurodevelopmental processes such as proliferation, migration, and differentiation of neuronal cells [11, 12]. Rab3A is the most common isoform in the human brain playing a pivotal role in neurotransmitter release and synaptic plasticity of RAB3GAP, which, in turn, converts active Rab3-GTP to inactive Rab3-GDP .
One of the major obstacles for clinicians presented with patients with a combination of different symptoms—e.g. postnatal growth retardation, microcephaly, microphthalmia, and cataracts—is to make a differential diagnosis ; however, using some advanced technologies such as whole-exome sequencing (WES), to some extent, has addressed these concerns and helped to make as precisely as diagnosis.
Although different studies have been conducted to investigate into the genotype-phenotype correlation in WARBM patients, yet this still remains blanketed in mystery in different populations. To fill this gap, we describe the clinical characteristics of two Iranian patients with WARBM type 1. This study may provide insight into the pathogenesis and clinical spectrum and facilitate early diagnosis of this syndrome.
Two patients—a 5-year-old female and a 4.5-year-old male—were referred to our department due to postnatal microcephaly and ocular manifestations including congenital cataract, microphthalmia, and optic atrophy. Pregnancies were uneventful; the anthropometric data (e.g. weight, length, and head circumference) in their first month of life were reported normal. No contractures were detected in both patients, while mild facial hypertrichosis was evident and limited to the temporal areas of their faces. Their parents were phenotypically normal and also consanguineous. The mothers (II.1 and II.3) had not any confrontation with infection, radiation, or even drug consumption during their pregnancies.
All of the patients’ clinical information and medical history were collected at the Department of Medical Genetics, Tarbiat Modares University, Tehran, Iran. Detailed patients’ history including disease onset, symptoms, progression, and also family history were gathered at this department. In addition to the traditional biochemistry tests, the patients were subjected to karyotyping using standard procedures to exclude any detectable chromosomal abnormalities. Neuroimaging was carried out using non-contrast brain magnetic resonance imaging (MRI).
The proband (III.1), a 5-year-old female, was born through normal vaginal delivery at 34 weeks of gestation without any complications. Her birth weight, length, and head circumference (HC) were reported normal. The motor developmental delay was firstly noted after the age of 6 months. After this age, the patient was not able to roll over or even sit on her own. On top of that, she had a severe speech and cognition delay. No seizure was reported by the parents until examination time. Cataract surgery was performed on the patient at the age of 4 years.
Physical and neurologic examination at the age of 5 years revealed a severe head lag, axial hypotonia, sparse voluntary movements, and peripheral spasticity in addition to severe speech and intellectual disability. Head circumference was measured 46 ± 0.2 cm (<− 3 SD) showing postnatal progressive microcephaly. From skeletal points of view, mild thoracolumbar scoliosis (Fig. 1a) and cortical thumb (Fig. 1b) were evident. Some craniofacial features—e.g. large low-set ears, congenital bilateral cataract, microphthalmia, and ptosis (Fig. 1b, c, d—were also apparent. No obvious abnormality was found regarding the patient’s feet (Fig. 1e).
Brain MRI without contrast including T1-W, T2-W, and FLAIR sequences was performed on the proband at the age of 5 years and revealed agenesis of the corpus callosum, cerebral white matter hypomyelination as diffuse bilateral abnormal signal changes in subcortical, deep, and periventricular white matter areas, cortical thickening especially in the frontal area, subcortical band heterotopia in the anterior part of the brain, and hypomyelination in the pre-dentate white matter of the cerebellum (Fig. 1f-l). Similar findings were obtained for the III.4. Further examinations showed no other abnormality using echocardiography, abdominal ultrasonography, and further skeletal investigations—i.e. examination for osteopetrosis, kyphoscoliosis, osteopenia, and growth hormone deficiency. Besides, using auditory brainstem response (ABR) and pure tone otoacoustic emissions (OAEs) testing, no hearing impairment was detected in the patient.
Hematological examination, thyroid, liver, and renal function tests, serum calcium, ammonia, lactate, pyruvate, and TORCH studies (serum anti-toxoplasmosis, rubella, cytomegalovirus, and herpes virus IgG and IgM antibodies) were all reported normal. Karyotyping did not reveal any detectable chromosomal abnormalities (Additional file 1). In Table 1, we summarized some of the important clinical findings in this patient, compared to different types of WARBM disorder and Martsolf syndrome (MIM: 212720).
The second case (III.4) was a 4.5-year-old male born through a normal vaginal delivery without any complications. He had normal birth weight and length and his birth HC was measured as 31 ± 0.2 cm (7th centile). No seizures were recorded in the patient’s medical history. According to the medical history, the motor developmental delay was first suspected after the age of 8 months. After this period, the patient was not able to roll over or sit down without any assistant. Furthermore, speech delay and cognitive disability were also noted after 1 year of age. Until the age of 4.5 years, no serious medical evaluation had been performed for a better diagnosis.
Clinical examinations at the age of 4.5 years showed quite similar findings in this patient such as noticeably delayed gross motor development with no head control, axial hypotonia, profound intellectual and speech abilities, and four limb spasticity. HC was measured 49 ± 0.2 cm (<− 3 SD) at the age of examination, indicating postnatal microcephaly. The general examination demonstrated a soft cleft palate, downturned corners of the mouth, wide nasal bridge, relatively short nose, mild micrognathia, a mild prominent forehead with bitemporal hollowing, bilateral low-set prominent ears, hypogenitalism with micropenis and bilateral cryptorchidism, and mild scoliosis (Fig. 1m-q). As a prominent clinical feature, the patient showed hypertrichosis of the upper back (Fig. 1r). Ocular abnormalities such as cataract, nasolacrimal duct obstruction, and bilateral microphthalmia were also detected at the time of examination.
Brain MRI without contrast—including T1-W, T2-W, and FLAIR sequences—performed at this age and showed the hypoplasia of the corpus callosum, cerebral white matter hypomyelination, cortical thickening (especially in the frontoparietal area), widened Sylvain fissure, enlarged ventricles, and hypomyelination in the white matter of the cerebellum. Electrocardiogram, bone scan, neonatal TORCH studies, and basic metabolic panel tests were all normal. G-banded chromosome analysis revealed a normal male karyotype as well (Additional file 1).
Genetic findings and determination of allele frequency of the candidate variants
Due to the heterogeneity of neurodevelopmental disorders , we subjected the patients, parents, and also a healthy sibling to paired-end WES screening to detect the possible underlying genetic factor. WES was performed based on the previous studies [16,17,18]. The mean depth of coverage was around 100× and approximately 98% of targeted regions were covered (Additional file 2). The overlapped filtering strategy was used to detect the common variants among the patients and the parents (Fig. 2a, b). The filtering steps were carried out according to the previous studies . In total, 780,385 common SNV and Indel variants among all the 7 genotyped individuals (are shown by an asterisk (*) in Fig. 2a) were detected. By doing the overlapped filter strategy, 4297 variants were detected. By excluding the synonymous, low-functional profile, and variants with mutation allele frequencies (MAF) greater than 1% in publicly available databases—e.g dbSNP150 , 1000 Genomes Project , Exome Sequencing Project , and gnomAD  or ExAC databases —only three variants were identified in MAP 3 K19, RAB3GAP1, and XIRP2 genes, respectively. Different tools, e.g. SIFT , Polyphen-2 , MutationTaster , and Provean  were also executed to evaluate the pathogenicity of the identified variants (Table 2). To investigate the co-segregation analysis, the primers were designed using Primer3.0 online tool  (Additional file 3) and the variants were subjected to Sanger sequencing (Additional file 4; Fig. 2c). The reference sequences of NM_001199144.1 for XIRP2, NM_001018044.2 for MAP3K19, and NM_012233.2 for RAB3GAP1 were used based on the human genome assembly GRCh37 (hg19).
In the nucleotide level, conservational analyses using ConSurf  and ‘2-Way Pseudogene Annotation Set’ of UCSC  showed that c.151 T is located in a highly conserved region (Fig. 2d). On the other hand, c.9835 T > C, p.(Ser3279Pro) in XIRP2 and c.332 T > C, p.(Ile111Thr) in MAP3K19 were categorized as non-conserved variants. MetaDome  was also used to visualize and confirm the data (Additional file 4).
To exclude the variants with low detrimental effects, we determined the allele frequency in 300 healthy and ethnically matched controls. For this purpose, we used PCR-Restriction Fragment Length Polymorphism (PCR-RFLP) and also tetra-primer Amplification Refractory Mutation System-PCR (tetra-primer ARMS-PCR). Using PCR-RFLP (Additional file 2), the allele frequency of c.9835 T > C variant in XIRP2 was determined 0.031, moreover, a homozygous individual was detected in healthy controls (Additional file 5). c.9835 T > C variant causes a loss of the XapI restriction site. The PCR products were digested with XapI and run on 2% agarose gel to assess whether the variant allele existed or not. This data was in line with the findings of Iranome  (allele frequency: 0.01063), suggesting that this variant cannot be considered pathogenic, so we excluded the variant from further investigations.
To determine the allele frequency of c.332 T > C in MAP3K19 and c.151-5 T > G in RAG3GAP1 in 300 healthy ethnicity-matched controls, tetra-primer ARMS-PCR was performed (Additional File 6; Fig. 3a). No homozygote/heterozygote was identified among 300 samples; consistently, Iranome and gnomAD verified these findings. The primers and size of expected bands in each reaction are depicted in Additional file 3. The procedure was carried out according to the previous study  (Additional file 2).
c.151-5T> G affects splicing
Using ‘MaxEnt Acceptor site algorithms’ of Human Splicing Finder V.3.1 (HSF; https://hsf.genomnis.com/) , we showed that c.151-5 T > G may alter the wild-type acceptor site, thus can most probably affect splicing. This finding was confirmed by Berkeley Drosophila Genome Project (BDGP; https://www.fruitfly.org)  and NetGene2 (http://www.cbs.dtu.dk/services/NetGene2)  tools. To predict if the variant can affect the cryptic splice-site activation and exon skipping, CRYP-SKIP (https://cryp-skip.img.cas.cz/) program  was used; these data showed that the c.151-5 T > G variant can increase the probability of cryptic splice site activation. To confirm, we used RT-PCR on the samples of the proband (III.1) and her parents (II.1 and II.2). The PCR products were then sequenced using Sanger sequencing with same primers (Additional file 3). The results revealed that the novel variant can activate a cryptic splice site within intron 3 of RAB3GAP1. This can change the splicing process of the RAB3GAP1 gene by adding 4 nucleotides between exons 3 and 4. In sum, this leads to the losing 930 amino acids (Fig. 3b-d).
The relative expression levels of RAB3GAP1 in the peripheral blood cells of the index case (III.1), parents (II.1 and II.2), and also a control/normal individual (III.3) were compared. The normal individual (control) was initially screened for the c.151-5 T > G variant and he was negative for this variant. RT-qPCR was performed using SYBR® Premix Taq™ II (TAKARA, Tokyo, Japan) on an ABI StepOne Sequence Detection System (Applied Biosystems, Foster City, USA). Each experiment was performed in triplicate and expression data were normalized to ALB (albumin) as an internal reference gene. Primer sequences are put forward in Additional file 3. The 2−ΔΔCt method was used for quantification by comparing the Ct values [37, 38]. The RT-qPCR analysis showed that the RAB3GAP1 expression was significantly decreased in the proband (Fig. 3e, f) and also demonstrated that only a trivial amount of RAB3GAP1 expression remained as a result of homozygous c.151-5 T > G splice site variant. This underscores the potential pathogenic effect of this novel variant (Fig. 3b, c). On the other hand, the decreased expression of RAB3GAP1 was determined in parents who were heterozygous for the novel variant (Fig. 3e, f).
To show whether the impaired genes make a contribution to brain development/function, we explored the regional expression of MAP3K19 and RAB3GAP1 in the normal adult human brain. For this purpose, we used microarray data (Affymetrix Exon 1.0 ST) from human post-mortem brain tissue collected by the UK Human Brain Expression Consortium . To confirm, we also used the data available in the Allen Mouse Brain Atlas [40, 41]. To show how the expression of Map3k19 (MGI: 1203481) and Rab3gap1 (MGI: 2445001) changes during brain development, we utilized the ‘Gene Expression Database (GXD) Project’  of MGI which is an international database resource for the laboratory mouse. The results showed that Map3k19 has a subtle expression during brain development that per se cannot justify a broad range of brain abnormalities in the patients. On the other hand, the aberrant expression of the RAB3GAP1 was significantly expected to be associated with the diverse brain impairments in WARBM disorders (Additional file 7).
Discussion and conclusions
WARBM is a heterogeneous group of autosomal recessive disorders that can be resulted from loss-of-function mutations in RAB3GAP1/2, RAB18, and TBC1D20 [43, 44]. These genes encode the ‘RAB proteins’ that mainly function as molecular switches ; they can also regulate membrane trafficking in a spatially and temporally restricted manner . GAP proteins can regulate RAB protein cycling , e.g., RAB3GAP1 (catalytic subunit) and RAB3GAP2 (non-catalytic subunit) make a heterodimeric enzyme complex, which has GAP activity that is specific for the RAB3 family of protein. This heterodimeric enzyme complex hydrolyzes GTP into GDP and also modulates the Ca2+ mediated exocytosis of hormones and neurotransmitters .
Most of the mutations in RAB3GAP1—in WARBM patients—are predicted to result in stop-codons and are distributed throughout the coding region from codons 89 to 934 . In a nutshell, all kinds of these mutations can lead to a truncated protein either before or within the regions that are important for catalytic activity  (Fig. 3c). This also can explain why there are some phenotypic variations between the WARBM with different mutations in RAB3GAP1. Rab3gap1 mutant and knockout mice models have previously been generated, while none of them embraced all major clinical features of WARBM . In mice, disruption of Rab3gap1 results in increased short-term depression and also paired-pulse facilitation at CA1 hippocampal neurons . On the other hand, human case reports have successfully covered a broad range of phenotypes caused by mutations in RAB3GAP1.
The majority of WARBM cases, similar to our case, are resulted from consanguineous marriages . Here, after performing WES, as a final candidate, we identified a homozygous pathogenic variant c.151-5 T > G; p.(Gly51IlefsTer15) in the intron 3 of the RAB3GAP1 gene, justifying the phenotypes that were in association with autosomal recessive WARBM type 1. Data from in silico analyses and also RT-PCR showed that c.151-5 T > G can activate a cryptic splice site in intron 3. This led to the addition of four nucleotides upstream of exon 4, which in turn changes the mRNA reading frame; this causes a loss of 971 codons, due to the premature stop codon. Moreover, RT-qPCR data showed that the expression of RAB3GAP1 harboring c.151-5 T > G variant was significantly decreased in the patients than the normal individual (III.3). Although the parents showed lower mRNA levels than the control group, they did not show any symptoms of the disease. This can per se confirm that the low amount of RAB3GAP1 protein is also sufficiently functional. There are two fates for mRNAs containing premature termination codons: nonsense-mediated mRNA decay (NMD)  or translation to truncated proteins . NMD is an evolutionarily conserved quality control pathway in eukaryotic cells that is responsible for inspecting mRNA for any possible errors, so eliminating any error-containing transcripts and controlling the amount of non-mutated transcript in the transcriptome (reviewed in ). Therefore, NMD leads to loss-of-function allele ; herein, RT-qPCR data supported the NMD hypothesis. Axiomatically, RAB3GAP1 protein that embraces p.(Gly51IlefsTer15) will be destroyed or be malfunctioned. It was previously suggested that pathogenic mutations in RAB3GAP1 with WARBM type 1 were mostly frameshift, nonsense, and splicing mutations causing complete loss of function ; Our data underscore this hypothesis as well. We also reclassified the variant based on the American College of Medical Genetics and Genomics (ACMG) variant interpretation guideline  into the ‘Pathogenic’ variant.
Homozygous and compound heterozygous mutations in human RAB3GAP1 cause WARBM type 1 and also a somehow milder phenotype called Martsolf syndrome . These two disorders are genetically heterogeneous and comprise a phenotypic spectrum described as “RAB18 deficiency” based on molecular etiology . Mutations in RAB3GAP2 can cause Martsolf syndrome. To exclude this syndrome, we meticulously compared the phenotypes of the patients in Table 1. Indeed, patients with Martsolf syndrome manifest different degrees of ID—from mild to severe—while WARBM patients show often severe-to-profound ID with no or even very limited communication skills. Spastic quadriplegia has been reported in WARBM, whereas in Martsolf syndrome, spasticity is limited to lower limbs . Molecularly, mutations reported in association with WARBM are mostly predicted to result in NMD and/or loss-of-protein-function . On top of that, we also found that some clinical features, e.g. hypoplasia of the corpus callosum, microcephaly, and developmental delay, can be observed in other types of WARBM in addition to Martsolf syndrome. Besides, based on some common clinical features among these patients and Martsolf syndrome, e.g. prominent antitragus, malaligned/misaligned teeth, and small penis, we suggest using the “RABopathies” term that can efficiently cover all clinical features between different types of WARBM and Martsolf syndrome.
We also introduced some novel manifestations of WARBM type 1 including abnormal philtrum, prominent antitragus, downturned corners of the mouth, malaligned teeth, scrotal hypoplasia, low anterior hairline in temporal areas, hypertrichosis of upper back, spastic diplegia to quadriplegia, and signal changes in deep and subcortical white matter areas. Perhaps, one of the most recognizable features of WARBM is the brain phenotype. Herein, clinical examinations showed that the patients (III.1 and III.4) manifested postnatal microcephaly which has been observed in most cases of the Martsolf syndrome and the WARBM. Developmental milestones in two patients were compatible with severe ID (Table 1).
The present study demonstrated that a novel homozygous variant in RAB3GAP1 is the main cause of WARBM type 1 disorder characterized by abnormal philtrum, prominent antitragus, downturned corners of the mouth, cerebral white matter signal changes, malaligned teeth, scrotal hypoplasia, low anterior hairline, hypertrichosis of upper back, spastic diplegia to quadriplegia as novel manifestations in WARBM type 1. We also coined the term “RABopathies” that can cover the overlapped clinical features among WARBM and Martsolf syndrome. Moreover, we reconfirmed that the WES can properly detect underlying genetic factors contributing to the neurodevelopmental disorder that show a great phenotype overlapping. This study also provides important information to guide genetic counseling and expands the genotype-phenotype spectrum of WARBM type 1.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. The variant and pertinent phenotypes caused by a mutation in RAB3GAP1 are accessible at ClinVar (SCV001426394) and Leiden Open Variation Database (LOVD; individual number: 279559).
American College of Medical Genetics and Genomics
Gene Expression Database
Magnetic resonance imaging
PCR-Restriction Fragment Length Polymorphism
Nonsense-mediated mRNA decay
Quantitative reverse transcription-PCR
- Tetra-primer ARMS-PCR:
Tetra-primer Amplification Refractory Mutation System-PCR
Warburg micro syndrome
Bem D, Yoshimura S-I, Nunes-Bastos R, Bond FF, Kurian MA, Rahman F, et al. Loss-of-function mutations in RAB18 cause Warburg micro syndrome. Am J Hum Genet. 2011;88(4):499–507. https://doi.org/10.1016/j.ajhg.2011.03.012.
Dursun F, Güven A, Morris-Rosendahl D. Warburg micro syndrome. J Pediatr Endocrinol Metab. 2012;25(3–4):379–82. https://doi.org/10.1515/jpem-2011-0459.
Morris-Rosendahl DJ, Segel R, Born AP, Conrad C, Loeys B, Brooks SS, et al. New RAB3GAP1 mutations in patients with Warburg micro syndrome from different ethnic backgrounds and a possible founder effect in the Danish. Eur J Hum Genet. 2010;18(10):1100–6. https://doi.org/10.1038/ejhg.2010.79.
Borck G, Wunram H, Steiert A, Volk AE, Körber F, Roters S, et al. A homozygous RAB3GAP2 mutation causes Warburg micro syndrome. Hum Genet. 2011;129(1):45–50. https://doi.org/10.1007/s00439-010-0896-2.
Liegel RP, Handley MT, Ronchetti A, Brown S, Langemeyer L, Linford A, et al. Loss-of-function mutations in TBC1D20 cause cataracts and male infertility in blind sterile mice and Warburg micro syndrome in humans. Am J Hum Genet. 2013;93(6):1001–14. https://doi.org/10.1016/j.ajhg.2013.10.011.
Aligianis IA, Morgan NV, Mione M, Johnson CA, Rosser E, Hennekam RC, et al. Mutation in Rab3 GTPase-activating protein (RAB3GAP) noncatalytic subunit in a kindred with Martsolf syndrome. Am J Hum Genet. 2006;78(4):702–7. https://doi.org/10.1086/502681.
Fukui K, Sasaki T, Imazumi K, Matsuura Y, Nakanishi H, Takai Y. Isolation and characterization of a GTPase activating protein specific for the Rab3 subfamily of small G proteins. J Biol Chem. 1997;272(8):4655–8. https://doi.org/10.1074/jbc.272.8.4655.
Ng EL, Tang BL. Rab GTPases and their roles in brain neurons and glia. Brain Res Rev. 2008;58(1):236–46. https://doi.org/10.1016/j.brainresrev.2008.04.006.
Coppola T, Perret-Menoud V, Lüthi S, Farnsworth CC, Glomset JA, Regazzi R. Disruption of Rab3–calmodulin interaction, but not other effector interactions, prevents Rab3 inhibition of exocytosis. EMBO J. 1999;18(21):5885–91. https://doi.org/10.1093/emboj/18.21.5885.
Tahara S, Sanno N, Teramoto A, Osamura RY. Expression of Rab3, a Ras-related GTP-binding protein, in human nontumorous pituitaries and pituitary adenomas. Mod Pathol. 1999;12(6):627–34.
Iwasaki K, Toyonaga R. The Rab3 GDP/GTP exchange factor homolog AEX-3 has a dual function in synaptic transmission. EMBO J. 2000;19(17):4806–16. https://doi.org/10.1093/emboj/19.17.4806.
Khvotchev MV, Ren M, Takamori S, Jahn R, Südhof TC. Divergent functions of neuronal Rab11b in Ca2+−regulated versus constitutive exocytosis. J Neurosci. 2003;23(33):10531–9. https://doi.org/10.1523/jneurosci.23-33-10531.2003.
Handley MT, Aligianis IA. RAB3GAP1, RAB3GAP2 and RAB18: disease genes in Micro and Martsolf syndromes. Biochem Soc Trans. 2012;40(6). https://doi.org/10.1042/bst20120169.
Shotelersuk V, Desudchit T, Suwanwela N. Postnatal growth failure, microcephaly, mental retardation, cataracts, large joint contractures, osteoporosis, cortical dysplasia, and cerebellar atrophy. Am J Med Genet A. 2003;116(2):164–9. https://doi.org/10.1002/ajmg.a.10067.
Dajani DR, Llabre MM, Nebel MB, Mostofsky SH, Uddin LQ. Heterogeneity of executive functions among comorbid neurodevelopmental disorders. Sci Rep. 2016;6(1):36566. https://doi.org/10.1038/srep36566.
Gilani N, Razmara E, Ozaslan M, Abdulzahra IK, Arzhang S, Tavasoli AR, et al. A novel deletion variant in CLN3 with highly variable expressivity is responsible for juvenile neuronal ceroid lipofuscinoses. Acta Neurol Belg. 2021:1–12. https://doi.org/10.1007/s13760-021-01655-9.
Soleimanipour F, Razmara E, Rahbarizadeh F, Fallahi E, Khodaeian M, Tavasoli AR, et al. A novel missense variant in the LMNB2 gene causes progressive myoclonus epilepsy. Acta Neurol Belg. 2021:1–9. https://doi.org/10.1007/s13760-021-01650-0.
Razmara E, Azimi H, Tavasoli AR, Fallahi E, Sheida SV, Eidi M, et al. Novel neuroclinical findings of autosomal recessive primary microcephaly 15 in a consanguineous Iranian family. Eur J Med Genet. 2020;63(12):104096. https://doi.org/10.1016/j.ejmg.2020.104096.
Garshasbi M, Mahmoudi M, Razmara E, Vojdanian M, Aslani S, Farhadi E, et al. Identification of RELN variant p.(Ser2486Gly) in an Iranian family with ankylosing spondylitis; the first association of RELN and AS. Eur J Hum Genet. 2020;28(6):1–9. https://doi.org/10.1038/s41431-020-0573-4.
Sherry ST, Ward M-H, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11. https://doi.org/10.1093/nar/29.1.308.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. 1000 genome project data processing Subgroup. 2009. The sequence alignment/map format and samtools. Bioinformatics. 2009;25(16):2078–9. https://doi.org/10.1093/bioinformatics/btp352.
Auer PL, Johnsen JM, Johnson AD, Logsdon BA, Lange LA, Nalls MA, et al. Imputation of exome sequence variants into population-based samples and blood-cell-trait-associated loci in African Americans: NHLBI GO exome sequencing project. Am J Hum Genet. 2012;91(5):794–808. https://doi.org/10.1016/j.ajhg.2012.08.031.
Wang Q, Pierce-Hoffman E, Cummings BB, Alföldi J, Francioli LC, Gauthier LD, et al. Landscape of multi-nucleotide variants in 125,748 human exomes and 15,708 genomes. Nat Commun. 2020;11(1):1–13. https://doi.org/10.1038/s41467-019-12438-5.
Karczewski KJ, Weisburd B, Thomas B, Solomonson M, Ruderfer DM, Kavanagh D, et al. The ExAC browser: displaying reference data information from over 60 000 exomes. Nucleic Acids Res. 2017;45(D1):D840–D5. https://doi.org/10.1093/nar/gkw971.
Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31(13):3812–4. https://doi.org/10.1093/nar/gkg509.
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;76(1):7.20. 1–7. 41. https://doi.org/10.1002/0471142905.hg0720s76.
Schwarz JM, Rödelsperger C, Schuelke M, Seelow D. MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods. 2010;7(8):575–6. https://doi.org/10.1038/nmeth0810-575.
Choi Y, Chan AP. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015;31(16):2745–7. https://doi.org/10.1093/bioinformatics/btv195.
Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JA. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35(suppl_2):W71–W4. https://doi.org/10.1093/nar/gkm306.
Glaser F, Pupko T, Paz I, Bell RE, Bechor-Shental D, Martz E, et al. ConSurf: identification of functional regions in proteins by surface-mapping of phylogenetic information. Bioinformatics. 2003;19(1):163–4. https://doi.org/10.1093/bioinformatics/19.1.163.
Karolchik D, Baertsch R, Diekhans M, Furey TS, Hinrichs A, Lu Y, et al. The UCSC genome browser database. Nucleic Acids Res. 2003;31(1):51–4. https://doi.org/10.1093/nar/gkg129.
Wiel L, Baakman C, Gilissen D, Veltman JA, Vriend G, Gilissen C. MetaDome: pathogenicity analysis of genetic variants through aggregation of homologous human protein domains. Hum Mutat. 2019;40(8):1030–8. https://doi.org/10.1002/humu.23798.
Fattahi Z, Beheshtian M, Mohseni M, Poustchi H, Sellars E, Nezhadi SH, et al. Iranome: a catalog of genomic variations in the Iranian population. Hum Mutat. 2019;40(11):1968–84. https://doi.org/10.1002/humu.23880.
Daneshjoo O, Garshasbi M. Novel compound heterozygote mutations in the ATP7B gene in an Iranian family with Wilson disease: a case report. J Med Case Rep. 2018;12(1):1–7. https://doi.org/10.1186/s13256-018-1608-0.
Desmet F-O, Hamroun D, Lalande M, Collod-Béroud G, Claustres M, Béroud C. Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009;37(9):e67–e. https://doi.org/10.1093/nar/gkp215.
Spradling AC, Stern D, Beaton A, Rhem EJ, Laverty T, Mozden N, et al. The Berkeley Drosophila genome project gene disruption project: single P-element insertions mutating 25% of vital Drosophila genes. Genetics. 1999;153(1):135–77.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Hayat Nosaeid M, Mahdian R, Jamali S, Maryami F, Babashah S, Maryami F, et al. Validation and comparison of two quantitative real-time PCR assays for direct detection of DMD/BMD carriers. Clin Biochem. 2009;42(12):1291–9. https://doi.org/10.1016/j.clinbiochem.2009.04.016.
Soreq L, Rose J, Soreq E, Hardy J, Trabzuni D, Cookson MR, et al. Major shifts in glial regional identity are a transcriptional hallmark of human brain aging. Cell Rep. 2017;18(2):557–70. https://doi.org/10.1016/j.celrep.2016.12.011.
Wang Q, Ding S-L, Li Y, Royall J, Feng D, Lesnar P, et al. The Allen mouse brain common coordinate framework: a 3D reference atlas. Cell. 2020;181(4):936–953.e20. https://doi.org/10.1016/j.cell.2020.04.007.
Xiong J, Ren J, Luo L, Horowitz M. Mapping histological slice sequences to the Allen mouse brain atlas without 3D reconstruction. Front Neuroinform. 2018;12:93. https://doi.org/10.3389/fninf.2018.00093.
Hill DP, Begley DA, Finger JH, Hayamizu TF, McCright IJ, Smith CM, et al. The mouse Gene Expression Database (GXD): updates and enhancements. Nucleic Acids Res. 2004;32(suppl_1):D568–D71. https://doi.org/10.1093/nar/gkh069.
Kabzinska D, Mierzewska H, Senderek J, Kochanski A. Warburg micro syndrome type 1 associated with peripheral neuropathy and cardiomyopathy. Folia Neuropathol. 2016;3:273–81. https://doi.org/10.5114/fn.2016.62537.
Hozhabri H, Talebi M, Mehrjardi MY, De Luca A, Dehghani M. Martsolf syndrome with novel mutation in the TBC1D20 gene in a family from Iran. Am J Med Genet A. 2020;182(5):957–61. https://doi.org/10.1002/ajmg.a.61543.
Geyer M, Wittinghofer A. GEFs, GAPs, GDIs and effectors: taking a closer (3D) look at the regulation of Ras-related GTP-binding proteins. Curr Opin Struct Biol. 1997;7(6):786–92. https://doi.org/10.1016/s0959-440x(97)80147-9.
Pfeffer SR. Rab GTPases: master regulators of membrane trafficking. Curr Opin Cell Biol. 1994;6(4):522–6. https://doi.org/10.1016/0955-0674(94)90071-x.
Ingmundson A, Delprato A, Lambright DG, Roy CR. Legionella pneumophila proteins that regulate Rab1 membrane cycling. Nature. 2007;450(7168):365–9. https://doi.org/10.1038/nature06336.
Handley MT, Aligianis IA. RAB3GAP1, RAB3GAP2 and RAB18: disease genes in Micro and Martsolf syndromes. Biochem Soc Trans. 2012;40(6). https://doi.org/10.1042/bst20120169.
Sakane A, Manabe S, Ishizaki H, Tanaka-Okamoto M, Kiyokage E, Toida K, et al. Rab3 GTPase-activating protein regulates synaptic transmission and plasticity through the inactivation of Rab3. Proc Natl Acad Sci. 2006;103(26):10029–34. https://doi.org/10.1073/pnas.0600304103.
Tenawi S, Al Khudari R, Alasmar D. Novel mutation in the RAB3GAP1 gene, the first diagnosed Warburg Micro syndrome case in Syria. Oxford Med Case Rep. 2020;2020(4):omaa031. https://doi.org/10.1093/omcr/omaa031.
Brogna S, Wen J. Nonsense-mediated mRNA decay (NMD) mechanisms. Nat Struct Mol Biol. 2009;16(2):107–13. https://doi.org/10.1038/nsmb.1550.
Roest PA, Roberts RG, Sugino S, van Ommen G-JB, den Dunnen JT. Protein truncation test (PTT) for rapid detection of translation-terminating mutations. Hum Mol Genet. 1993;2(10):1719–21. https://doi.org/10.1093/hmg/2.10.1719.
Chang Y-F, Imam JS, Wilkinson MF. The nonsense-mediated decay RNA surveillance pathway. Annu Rev Biochem. 2007;76(1):51–74. https://doi.org/10.1146/annurev.biochem.76.050106.093909.
Asahina M, Endoh Y, Matsubayashi T, Fukuda T, Ogata T. Novel RAB3GAP1 compound heterozygous mutations in Japanese siblings with Warburg micro syndrome. Brain Dev. 2016;38(3):337–40. https://doi.org/10.1016/j.braindev.2015.09.006.
Biesecker LG, Harrison SM. The ACMG/AMP reputable source criteria for the interpretation of sequence variants. Genet Med. 2018;20(12):1687–8. https://doi.org/10.1038/gim.2018.42.
Gumus E. Case report of four siblings in Southeast Turkey with a novel RAB3GAP2 splice site mutation: Warburg micro syndrome or Martsolf syndrome? Ophthalmic Genet. 2018;39(3):391–5. https://doi.org/10.1080/13816810.2018.1432065.
Koparir A, Karatas OF, Yilmaz SS, Suer I, Ozer B, Yuceturk B, et al. Revealing the functions of novel mutations in RAB3GAP1 in Martsolf and Warburg micro syndromes. Am J Med Genet A. 2019;179(4):579–87. https://doi.org/10.1002/ajmg.a.61065.
We thank the family members for their willingness and contribution.
The funders provided roles in study design, execution, and data interpretation. This research received specific grants from de recherché du Québec–Santé, Québec, Canada (Grant Number: 269552), and Tarbiat Modares University, Tehran, Iran (Grant Number: TMU.8198.52D).
Ethics approval and consent to participate
The research protocol was approved by the ethics committee of Tarbiat Modares University, Tehran, Iran under the ethical code of “TMU.8198.52D.” The written informed consent was received from each guardian and they also provided a signed written consent to publish all personal and medical details included in this study in addition to the accompanying pictures. Beyond that, all methods were performed in accordance with the relevant guidelines and regulations of the Declaration of Helsinki.
Consent for publication
The informed consent to publish was obtained from each guardian.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The karyotypes of the proband (a) and the male affected individual or III.4 (b). The karyotypes did not show any obvious chromosomal changes
Detailed procedures of Whole-Exome Sequencing and Tetraplex-Amplification refractory mutation system (Tetraple-ARMS-PCR)
A list of primers that was used in this study
Sanger sequencing and conservational analysis for the candidate variants in MAP3K19 and XIRP2 genes. a) Chromatograms shed light on the cosegregation of c.332 T > C of MAP3K19 in the family members. b) Chromatograms showing nucleotide sequences of XIRP2 in the regions of c.9835 T > C. c) MetaDome shows that the c.332 T is located in a region with an average conservational profile in MAP3K19 protein. d) In this XIRP2 tolerance landscape, the region harboring the novel missense variant is partially tolerant in comparison with other parts in this protein. The affected region is located in a somehow variable region
2% agarose gel electrophoresis of tetra primer ARMS-PCR test products for c.332 T > C in the MAP3K19 gene. In this figure: well 1: 100 bp ladder, well 2: negative PCR test control (no template control; NTC), well 3: Homozygous for the normal allele, well 4: Heterozygous for the wild-type allele. The allele frequency was estimated at 0.031 for the c.332 T > C variant
A simple and rapid PCR-RFLP assay was used to detect the c.9835 T > C variant in the XIRP2 gene. a) This schematic figure shows the generated fragments after digestion with the XapI restriction enzyme. The variant causes losing the restriction enzyme site. Agarose gel (2.0%) electrophoresis with ethidium bromide staining following the XapI digestion of the PCR products is shown. PCR-RFLP results in normal control showing 340, 249, and 130 bp (T/T: wild-type allele); after XapI digestion, a heterozygous sample will show four distinct bands including 470, 340, 249, and 130 bp. The homozygote reveals three distinct bands consisting of 470, 249, and 130 bp. b) PCR-RFLP assay was used on the samples in order to show the genotyping. In this figure, 1: heterozygote, 2: patient, 3–14, and 16–22 are wild-type alleles
Gene expression of the RAB3GAP1 is changing in different brain regions as age increases. b) Gene expression of RAB3GAP1 in different brain regions and sexes. c) Boxplots of RAB3GAP1 mRNA expression levels in ten adult brain regions. The expression levels are based on exon array experiments and are plotted on a log2 scale (y-axis). This dataset was generated with Affymetrix Exon 1.0 ST arrays and brain tissue originating from 134 control individuals, collected by the Medical Research Council Sudden Death. The plot also shows significant variation in RAB3GAP1 expression across the ten brain regions analyzed, such that expression is higher in the temporal cortex than in any other region. In this figure, TCTX: temporal cortex, FCTX: frontal cortex, OCTX: occipital cortex, SNIG: substantia nigra, THAL: thalamus, PUTM: putamen, HIPP: hippocampus, CRBL: cerebellar cortex, MEDU: medulla (specifically the inferior olivary nucleus), WHMT: intralobular white matter, and N: number of samples analyzed for each brain region. These figures shed light on RAB3GAP1 expression in the different parts of the brain
About this article
Cite this article
Khalesi, R., Razmara, E., Asgaritarghi, G. et al. Novel manifestations of Warburg micro syndrome type 1 caused by a new splicing variant of RAB3GAP1: a case report. BMC Neurol 21, 180 (2021). https://doi.org/10.1186/s12883-021-02204-w
- Tetra-primer ARMS-PCR
- WARBM type 1
- Whole-exome sequencing