Re-analysis of whole-exome sequencing data reveals a novel splicing variant in the SLC2A1 in a patient with GLUT1 Deficiency Syndrome 1 accompanied by hemangioma: a case report
BMC Medical Genomics volume 14, Article number: 197 (2021)
GLUT1 Deficiency Syndrome 1 (GLUT1DS1) is a neurological disorder caused by either heterozygous or homozygous mutations in the Solute Carrier Family 2, Member 1 (SLC2A1) gene. SLC2A1 encodes Glucose transporter type 1 (GLUT1) protein, which is the primary glucose transporter at the blood–brain barrier. A ketogenic diet (KD) provides an alternative fuel for brain metabolism to treat impaired glucose transport. By reanalyzing exome data, we identified a de novo heterozygous SLC2A1 variant in a girl with epilepsy. After reversed phenotyping with neurometabolic tests, she was diagnosed with GLUT1DS1 and started on a KD. The patient's symptoms responded to the diet. Here, we report a patient with GLUT1DS1 with a novel SLC2A1 mutation. She also has a hemangioma which has not been reported in association with this syndrome before.
A 5-year 8-month girl with global developmental delay, spasticity, intellectual disability, dysarthric speech, abnormal eye movements, and hemangioma. The electroencephalography (EEG) result revealed that she had epilepsy. Magnetic resonance imaging (MRI) showed that non-specific white matter abnormalities. Whole Exome Sequencing (WES) was previously performed, but the case remained unsolved. The re-analysis of WES data revealed a heterozygous splicing variant in the SLC2A1 gene. Segregation analysis with parental DNA samples indicated that the variant occurred de novo. Lumbar puncture (LP) confirmed the diagnosis, and the patient started on a KD. Her seizures responded to the KD. She has been seizure-free since shortly after the initiation of the diet. She also had decreased involuntary movements, her speech became more understandable, and her vocabulary increased after the diet.
We identified a novel de novo variant in the SLC2A1 gene in a patient who previously had a negative WES result. The patient has been diagnosed with GLUT1DS1. The syndrome is a treatable condition, but the differential diagnosis is not an easy process due to showing a wide range of phenotypic spectrum and the overlapping symptoms with other neurological diseases. The diagnosis necessitates a genomic testing approach. Our findings also highlight the importance of re-analysis to undiagnosed cases after initial WES to reveal disease-causing variants.
The Solute Carrier Family 2, Member 1 gene (SLC2A1, OMIM *138140) is located on the 1p34.2 and encodes Glucose transporter type 1 (GLUT1) composed of 492 amino acids [1, 2]. GLUT1 is the primary membrane transporter responsible for the delivery of glucose in human erythrocytes and blood-tissue barriers [3,4,5]. Impaired glucose transport at the blood–brain barrier caused by mutations in the SLC2A1 result in GLUT1 Deficiency Syndrome 1 (GLUT1DS1, OMIM #606777), which is a rare neurometabolic disorder [6, 7]. In 1991, De Vivo et al. firstly identified GLUT1DS1 in the two patients with low glucose concentrations in cerebrospinal fluid (CSF)-termed hypoglycorrhachia, seizures, delayed development, dysarthria, acquired microcephaly, and movement disorder including spasticity, ataxia, and dystonia. Seven years after this report, Seidner et al. revealed that mutations in the SLC2A1 gene are the genetic factor behind the GLUT1DS1 .
The current diagnosis approaches are generally based on the molecular analysis of the SLC2A1 gene and lumbar puncture for CSF evaluation  since hypoglycorrhachia is the hallmark of the syndrome . The CSF/blood glucose ratio is usually less than 0.4 and CSF glucose level should be less than 60 mg/dl for the diagnosis . Another sign for the syndrome is the SLC2A1 mutations which predominantly occur de novo with the autosomal dominant condition . In familial cases, mutations are usually inherited by autosomal dominant pattern; yet, autosomal recessive variants have also been shown [12,13,14,15]. Thus far, it is estimated that more than 140 different pathogenic variants have been reported in the SLC2A1 . Despite these reports, the relationship between phenotype and genotype for the syndrome could not be clarified thoroughly . However, missense mutations can be associated with the whole symptom spectrum of GLUT1DS1, while deletions in the SLC2A1 are related to more severe forms of the syndrome. . Furthermore, the phenotypic range of GLUT1DS has been broadened in recent years . All these facts make a precise diagnosis for this syndrome a challenging task for clinicians. On the other hand, it is vital to identify GLUT1DS1 in the early stages since the ketogenic diet (KD) therapy has the potential to improve symptoms .
The patient was a 5-year 8-month-old girl at the time of presentation, with epilepsy, intellectual disability, and a movement disorder. She was born to a 29-year-old woman, G3P1, at term, with c- section. She had extreme crying at the age of 2.5 months. She also had left lower eyelid hemangioma at birth. She was diagnosed with epilepsy at the age of 4 months. She had head control at the age of 1 year; she sat at 3 years; she never walked without assistance. She had delayed speech. She had microcephaly (46.5 cm, < 3 p); she has a dysarthric speech. Her tonus was increased in all four extremities; she also had involuntary movements (choreoathetosis and dystonia). Her deep tendon reflexes were increased in the lower extremities. She could only walk with assistance. Her 3 T magnetic resonance imaging (MRI) showed non-specific millimetric white matter hyperintensities. Her EEG showed left temporoparietal and right posterior temporal interictal epileptiform discharges that also generalized. The patient was on Levetiracetam for seizures. Lumbar puncture (LP) was performed to confirm the genetic finding. LP showed no cells, protein 20,40 mg/dl, glucose 32 mg/dl, simultaneous blood glucose 87 mg/dl (before the lumbar puncture, after a 4-h fasting), and a negative CSF culture.
Whole exome sequencing data analysis
The quality of raw reads was checked using FastQC (v0.11.8) . Skewer (v0.2.2)  was used to remove adapters and low quality (Phred quality score < 20) bases. Trimmed reads were aligned to the reference human genome (hg19) using the Burrows-Wheeler Aligner (BWA mem, v0.7.17) . Post-alignment processing was performed following the Genome Analysis Toolkit’s (GATK, v126.96.36.199)  best practice recommendations. Picard tool package (v1.141)  was used to mark duplicate reads and fix the information of mate reads. Haplotyping and joint genotyping were also performed following GATK recommendation. Annovar (v2018-12-04)  was used to annotate each variant by getting biologically significant information from several databases.
Variant filtration strategy
After WES data analysis, 10,800 heterozygous and 5,500 homozygous variants passed the depth and quality filter. A phenotype-driven prioritization strategy was employed to detected variants. According to the mode of inheritance, variant filtering was performed by either absence-approach that assumes that disease-causing variants are not found in healthy population databases or setting the minor allele frequency threshold at < 0.1% . The variants that are more than ten bases away from the exon–intron boundaries were excluded to filter the splicing variants. After the population frequency and splice site filtering, the number of variants decreased to 402 and 15, respectively. Several gene-level assessments, including the genic intolerance score  and gnomAD metrics , and variant level assessments such as CADD , REVEL , M-CAP  were utilized to interpret the pathogenicity. Then, the variants were prioritized based on the symptoms of the affected individual by collecting evidence from various sources; Online Mendelian Inheritance in Man (OMIM) , Mouse Genome Informatics (MGI) , and literature search. The variants were visualized via the Integrative Genomics Viewer (IGV)  to coverage check. The whole bioinformatics workflow was shown in Fig. 1.
Combining all of these computational evaluations provided strong evidence that the novel variant in the SLC2A1 gene causes the observed phenotype. The heterozygous mutation in SLC2A1 (NM_006516), c.275 + 1del classified as “Uncertain significance” based on the 2015 guidelines of the American College of Medical Genetics and Genomics (ACMG) . The mutation was detected with high-depth reading in IGV (Fig. 2).
Parental segregation analysis
Since the patient was a proband-only WES case, Sanger sequencing was carried out for parental segregation analysis. Peripheral blood samples were collected from the patient and her biological parents following the standard procedures. Genomic DNA was isolated from peripheral blood samples using the Quick-DNA™ Miniprep Plus Kit (Zymo Research, USA) following the manufacturer’s protocol. DNA quality was checked by agarose gel electrophoresis and NanoDrop 2000c Spectrophotometer (Thermo Fisher Scientific, USA). Sanger sequencing was carried out using standard protocols. PCR conditions and primer sequences can be provided upon request. Sequence electropherograms were visualized using the 4Peaks (Mekentosj, Amsterdam).
Sanger results indicate that the variant, c.275 + 1del, occurred de novo. Biological parents of the patient have reference alleles in the same position. In addition to the parental verification, mixed and low-quality traces from the heterozygous deletion point in both forward and reverse sequence direction confirm the IGV visualization (Fig. 3).
Discussion and conclusions
Clinical exome sequencing is considered as a powerful approach for identifying disease-causing variants, even though it has a 25–30% diagnostic yield for patients . Some of these undiagnosed patients already have a pathogenic variant in WES, but it may not be identified in the initial analysis. A recent study showed that 10% of the undiagnosed patients could get a precise diagnosis through the re-analysis of the same WES data, by applying different workflows and with the help of growing knowledge in the literature .
Here, we report a patient who had previously unsolved WES data. The in-house variant prioritization workflow was employed to the raw data of the patient. The workflow revealed the novel heterozygous variant c.275 + 1del in the SLC2A1 gene. To our knowledge, this variant has not been reported in the literature or the available databases. Sanger results indicate that the variant has occurred de novo. After extensive neurometabolic and genetic screen, the patient was diagnosed with GLUT1DS1. Then the patient started on the KD, and her seizures responded to the KD. She has been seizure-free shortly after the initiation of the diet. She also had decreased involuntary movements, her speech became more understandable, and her vocabulary increased after the diet. Thus, this study has also shown that the importance of dietary therapeutic approaches after the investigation of genetic causes behind a disease.
Interestingly, the missense splicing variant c.275 + 1G > A, which is in the same position as our patient, was reported before by Ismayilova et al.  (Table 1). The nucleotide sequences on the exon–intron boundaries are usually highly conserved and crucial for the regulation of gene expression [40, 41]. Mutations at these splicing sites may affect the gene expression drastically, which can explain the disease phenotype occurred in both individuals. Although both patients had some phenotypic overlap, including global developmental delay, seizure, signal abnormality in white matter, our patient showed some additional symptoms including microcephaly and several complex movement anomalies such as dystonia and choreoathetosis. However, the previously reported patient has a lower CSF glucose level and CSF/serum ratio. The symptomatic differences between these cases may arise from alternate modifying genetic features and epigenetic factors. The causative mechanisms have to be further investigated to understand the dissimilarity.
In addition to typical symptoms of GLUT1DS1, our patient also has a hemangioma as an unusual feature. Although the hemangioma pathogenesis is not completely clarified, GLUT1 protein has been used as a selective marker to differentiate hemangioma from other vascular malformations [42, 43]. It is supposed that the expression of GLUT1 in hemangioma tissue is a potential sign of the pathogenicity because the protein is not expressed in the healthy skin tissue . Even though GLUT1 plays a role in both conditions, this occurrence can be coincidental. Further molecular and bioinformatic analyses are needed to establish whether there is a relationship between hemangioma and GLUT1DS1, but, based on our knowledge, these two conditions were not previously reported together.
Availability of data and materials
The novel variant revealed during the study has been submitted to ClinVar repository with the accession number SCV001548248. The datasets analyzed during the study are available in the NCBI Sequence Read Archive (SRA) with the bioproject accession number PRJNA706171. The Sanger Sequencing data generated in the study has been submitted to NCBI GenBank BankIt with the accession numbers MW883607 and MW883608. Reference sequences used in this study for SLC2A1 are available in the following link (https://www.ncbi.nlm.nih.gov/gene/6513#reference-sequences). Databases used in this study were human hg19 reference genome assembly (http://hgdownload.soe.ucsc.edu/goldenPath/hg19/bigZips/hg19.fa.gz), ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar), dbSNP (https://www.ncbi.nlm.nih.gov/snp/), gnomAD Browser (https://gnomad.broadinstitute.org/).
American College of Medical Genetics and Genomics
Combined Annotation Dependent Depletion
Genome Analysis Toolkit
Glucose Transporter Type 1
GLUT1 Deficiency Syndrome 1
The Genome Aggregation Database
The Integrative Genomics Viewer
Mendelian Clinically Applicable Pathogenicity;
Mouse Genome Informatics
Magnetic Resonance Imaging
Online Mendelian Inheritance in Man
Rare Exome Variant Ensemble Learner
- SLC2A1 :
Solute Carrier Family 2, Member 1
Whole Exome Sequencing
Wang D, Pascual JM, Yang H, Engelstad K, Jhung S, Sun RP, et al. Glut-1 deficiency syndrome: Clinical, genetic, and therapeutic aspects. Ann Neurol. 2005;57:111–8.
Baldwin SA. Mammalian passive glucose transporters: members of an ubiquitous family of active and passive transport proteins. Biochimica et Biophysica Acta (BBA) - Reviews on Biomembranes. 1993;1154:17–49.
Gerhart DZ, Levasseur RJ, Broderius MA, Drewes LR. Glucose transporter localization in brain using light and electron immunocytochemistry. J Neurosci Res. 1989;22:464–72.
Harik SI, Kalaria N, Andersson L, Lundahl P. lmmunocytochemical localization of the erythroid glucose transporter: abundance in tissues with barrier functions. J Neurosci. 1990;10:3862–72.
Pardridges M, Boado J, Farrell R. Brain-type glucose transporter (GLUT-1) is selectively localized to the blood-brain barrier. J Biol. 1990;265:18035–40.
De Vivo DC, Trifiletti RR, Jacobson RI, Ronen GM, Behmand RA, Harik SI. Defective glucose transport across the blood-brain barrier as a cause of persistent hypoglycorrhachia, seizures, and developmental delay. N Engl J Med. 1991;325:703–9.
Larsen J, Johannesen KM, Ek J, Tang S, Marini C, Blichfeldt S, et al. The role of SLC2A1 mutations in myoclonic astatic epilepsy and absence epilepsy, and the estimated frequency of GLUT1 deficiency syndrome. Epilepsia. 2015;56:e203–8.
Seidner G, Alvarez MG, Yeh JI, O’Driscoll KR, Klepper J, Stump TS, et al. GLUT-1 deficiency syndrome caused by haploinsufficiency of the blood-brain barrier hexose carrier. Nat Genet. 1998;18:188–91.
De Giorgis V, Veggiotti P. GLUT1 deficiency syndrome 2013: Current state of the art. Seizure. 2013;22:803–11.
Pearson TS, Akman C, Hinton VJ, Engelstad K, De Vivo DC. Phenotypic spectrum of glucose transporter type 1 deficiency syndrome (Glut1 DS). Curr Neurol Neurosci Rep. 2013;13:342.
Verrotti A, D’Egidio C, Agostinelli S, Gobbi G. Glut1 deficiency: When to suspect and how to diagnose? Eur J Paediatr Neurol. 2012;16:3–9.
Brockmann K, Wang D, Korenke CG, Von Moers A, Ho Y-Y, Pascual JM, et al. Autosomal dominant Glut-1 deficiency syndrome and familial epilepsy. Ann Neurol. 2001;50:476–85.
Klepper J. Autosomal dominant transmission of GLUT1 deficiency. Hum Mol Genet. 2001;10:63–8.
Rotstein M, Engelstad K, Yang H, Wang D, Levy B, Chung WK, et al. Glut1 deficiency: inheritance pattern determined by haploinsufficiency. Ann Neurol. 2010;68:955–8.
Klepper J, Scheffer H, Elsaid MF, Kamsteeg E-J, Leferink M, Ben-Omran T. Autosomal recessive inheritance of GLUT1 deficiency syndrome. Neuropediatrics. 2009;40:207–10.
Di Vito L, Licchetta L, Pippucci T, Baldassari S, Stipa C, Mostacci B, et al. Phenotype variability of GLUT1 deficiency syndrome: description of a case series with novel SLC2A1 gene mutations. Epilepsy Behav. 2018;79:169–73.
Koch H, Weber YG. The glucose transporter type 1 (Glut1) syndromes. Epilepsy Behav. 2019;91:90–3.
Zaman SM, Mullen SA, Petrovski S, Maljevic S, Gazina EV, Phillips AM, et al. Development of a rapid functional assay that predicts GLUT1 disease severity. Neurol Genet. 2018;4:e297.
Castellotti B, Ragona F, Freri E, Solazzi R, Ciardullo S, Tricomi G, et al. Screening of SLC2A1 in a large cohort of patients suspected for Glut1 deficiency syndrome: identification of novel variants and associated phenotypes. J Neurol. 2019;266:1439–48.
Klepper J, Diefenbach S, Kohlschütter A, Voit T. Effects of the ketogenic diet in the glucose transporter 1 deficiency syndrome. Prostaglandins Leukot Essent Fatty Acids. 2004;70:321–7.
The Babraham Bioinformatics group. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Jiang H, Lei R, Ding S-W, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15:182.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
DePristo MA, Banks E, Poplin RE, Garimella KV, Maguire JR, Hartl C. A framework for variation discovery and genotyping using next- generation DNA sequencing data. Nat Genet. 2011;43:491–8.
Broad Institute. Picard. 2020. http://broadinstitute.github.io/picard/.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164–e164.
Bamshad MJ, Ng SB, Bigham AW, Tabor HK, Emond MJ, Nickerson DA, et al. Exome sequencing as a tool for Mendelian disease gene discovery. Nat Rev Genet. 2011;12:745–55.
Petrovski S, Wang Q, Heinzen EL, Allen AS, Goldstein DB. Genic intolerance to functional variation and the interpretation of personal genomes. PLoS Genet. 2013;9:e1003709.
Genome Aggregation Database Consortium, Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581:434–43.
Kircher M, Witten DM, Jain P, O’Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014;46:310–5.
Ioannidis NM, Rothstein JH, Pejaver V, Middha S, McDonnell SK, Baheti S, et al. REVEL: an ensemble method for predicting the pathogenicity of rare missense variants. Am J Hum Genet. 2016;99:877–85.
Jagadeesh KA, Wenger AM, Berger MJ, Guturu H, Stenson PD, Cooper DN, et al. M-CAP eliminates a majority of variants of uncertain significance in clinical exomes at high sensitivity. Nat Genet. 2016;48:1581–6.
Amberger JS, Bocchini CA, Scott AF, Hamosh A. OMIM.org: leveraging knowledge across phenotype–gene relationships. Nucleic Acids Res. 2019;47:D1038–43.
Bult CJ, Blake JA, Smith CL, Kadin JA, Richardson JE, The Mouse Genome Database Group, et al. Mouse Genome Database (MGD) 2019. Nucleic Acids Res. 2019;47:D801–6.
Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.
Richards S, Aziz N, Bale S, Bick D, Das S, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–23.
Yang Y, Muzny DM, Reid JG, Bainbridge MN, Willis A, Ward PA, et al. Clinical whole-exome sequencing for the diagnosis of mendelian disorders. N Engl J Med. 2013;369:1502–11.
Wenger AM, Guturu H, Bernstein JA, Bejerano G. Systematic reanalysis of clinical exome data yields additional diagnoses: implications for providers. Genet Med. 2017;19:209–14.
Ismayilova N, Hacohen Y, MacKinnon AD, Elmslie F, Clarke A. GLUT-1 deficiency presenting with seizures and reversible leukoencephalopathy on MRI imaging. Eur J Paediatr Neurol. 2018;22:1161–4.
Chang W-C. Alternative splicing of U12-type introns. Front Biosci. 2008;13:1681.
Parada GE, Munita R, Cerda CA, Gysling K. A comprehensive survey of non-canonical splice sites in the human transcriptome. Nucleic Acids Res. 2014;42(16):10564–78.
North PE, Waner M, Mizeracki A, Mihm MC Jr. GLUT1: a newly discovered immunohistochemical marker for juvenile hemangiomas. Hum Pathol. 2000;31(1):11–22.
Leon-Villapalos J, Wolfe K, Kangesu L. GLUT-1: an extra diagnostic tool to differentiate between haemangiomas and vascular malformations. Br J Plast Surg. 2005;58:348–52.
We thank the family for participating in this study and contributing to the literature.
There was no funding received for the conduction of the study.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board (IRB) of Acibadem Mehmet Ali Aydinlar University (ATADEK-2019-17/36). The parent provided written consent on behalf of the child participant. The written consent was also received from the patient’s parents for their own participation in the study.
Consent for publication
All family members and the parent on behalf of the child provided written informed consent for publication. The institutional ethics committee approved this procedure.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Bozkurt, T., Alanay, Y., Isik, U. et al. Re-analysis of whole-exome sequencing data reveals a novel splicing variant in the SLC2A1 in a patient with GLUT1 Deficiency Syndrome 1 accompanied by hemangioma: a case report. BMC Med Genomics 14, 197 (2021). https://doi.org/10.1186/s12920-021-01045-3