- Case report
- Open Access
- Published:
Novel frameshift mutation in LIS1 gene is a probable cause of lissencephaly: a case report
BMC Pediatrics volume 22, Article number: 545 (2022)
Abstract
Background
Lissencephaly (LIS) is a cortical malformation, characterized by smooth or nearly smooth cerebral surface and a shortage of gyral and sulcal development, which is caused by deficient neuronal migration during embryogenesis. Neuronal migration involves many gene products, among which is the product of the PAFAH1B1 gene, associated with this disease. LIS is a rare disease, characterized by low population frequency, and with non-specific clinical symptoms such as early epilepsy, developmental delay or cerebral palsy-like motor problems. Given that high-throughput sequencing techniques have been improving diagnosis, we have chosen this technique for addressing this patient.
Case presentation
We present the case of a seven years old male patient with an undiagnosed rare disease, with non-specific clinical symptoms possibly compatible with lissencephaly.
The patient was enrolled in a study that included the sequencing of his whole genome. Sequence data was analyzed following a bioinformatic pipeline. The variants obtained were annotated and then subjected to different filters for prioritization. Also mitochondrial genome was analyzed. A novel candidate frameshift insertion in known PAFAH1B1 gene was found, explaining the index case phenotype. The assessment through in silico tools reported that it causes nonsense mediated mechanisms and that it is damaging with high confidence scores. The insertion causes a change in the reading frame, and produces a premature stop codon, severely affecting the protein function and probably the silencing of one allele. The healthy mother did not carry the mutation, and the unaffected father was not available for analysis.
Conclusions
Through this work we found a novel de novo mutation in LIS1/PAFAH1B1 gene, as a likely cause of a rare disease in a young boy with non-specific clinical symptoms. The mutation found correlates with the phenotype studied since the loss of function in the gene product has already been described in this condition. Since there are no other variants in the PAFAH1B1 gene with low population frequency and due to family history, a de novo disease mechanism is proposed.
Background
Lissencephaly (LIS) is a subtype of malformations of cortical development (MCD), which are a heterogenous group of disorders with diverse phenotypic and genotypic presentations. Patients with LIS may present different degrees of developmental delays, seizures, severe psychomotor impairment, muscle spasticity or hypotonia [1].
Lissencephaly is a disorder caused by a defect in neuronal migration, which occurs between 12 and 24 weeks of gestation and results in a lack of development of brain folds (gyri) and grooves (sulci) [2]. Neuronal migration is a complex process, which requires the coordination of many gene products.
LIS1 and DCX were the first genes that were associated with LIS, discovered in 1993 and 1998, respectively [3, 4]. In the past years with the advent of new molecular genomics technologies, many additional genes were found. These LIS-related genes include ACTB, ACTG1, ARX, CDK5, CRADD, DYNC1H1, KIF2A, KIF5C, NDE1/NDEL1, TUBA1A, TUBA8, TUBB, TUBB2B, TUBB3, TUBG1, RELN and VLDLR. Many of these 19 LIS-associated genes are related to microtubule structural proteins (tubulin) or microtubule-associated proteins [5].
The PAFAH1B1 gene (Genbank accession number: NM_000430), located at chromosome 17p13.3, encodes the alpha subunit of the 1B isoform of the platelet-activation factor acetylhydrolase regulatory, a highly conserved protein of 410 amino acids, known as LIS1 or PAFAH1B1 [6]. It has two protein coding transcripts and several non-coding ones. LIS1/PAFAH1B1 forms the non-catalytic subunit of the G protein-like heterotrimeric cytosolic platelet-activating factor acetylhydrolase (PAF-AH) brain isoform Ib (PAFAH1B1) [7]. Along with two other subunits, PAFAH1B2 and PAFAH1B3, LIS1/PAFAH1B1 forms a trimeric complex which regulates the level of platelet activating factor (PAF) in the brain, by catalyzing the removal of the acetyl group at the SN-2 position of platelet-activating factor [8, 9]. The regulation of optimal concentrations of PAF in the brain may be critical for correct neuronal migration, essential for normal brain development and function. LIS1/PAFAH1B1 has also been shown to play a central role in the organization of the cytoskeleton, which in turn affects neuronal proliferation and migration [6]. Mutations in this gene have previously been associated with cortical brain malformation in children (Table 1).
Here, we comment on the case of a seven years old boy with an undiagnosed rare disease, with non-specific symptoms that could be compatible with LIS, but with an unclear presentation. Whole genome sequencing (WGS) of the patient was performed in the context of a genomics project (urugenomes.org) and sequence data was analyzed following a bioinformatics pipeline which concluded with a set of annotated and prioritized variants. A novel candidate frameshift variant was found that fits with the boy’s phenotype. To support the pathogenicity of the variant we used computational prediction tools and made segregation analysis with Sanger sequencing. To the best of our knowledge this is the first time this variant is reported [11] and it is the most likely cause of the patient’s disease.
Case presentation
The index case is a seven years old boy with perinatal clinical records of poorly controlled pregnancy, homelessness and multiple drugs abuse through all gestation. Both parents are healthy and non-consanguineous. Early term delivery, low birthweight with microcephaly at birth with a head circumference of 32 cm (Z score -3).
The patient develops a spastic bilateral cerebral palsy, Gross Motor Function Classification System (GMFSC) IV, Bimanual fine motor function (BFMF) IVa, Communication Function Classification System (CFCS) IV. Associated with a profound intellectual disability, visual and auditory sensory deficit, pharmacoresistant epilepsy with generalized tonic clonic seizures, and congenital microcephaly with a head circumference growth always under -3 standard deviations. No dysmorphic signs were detected.
Cerebral magnetic resonance imaging (MRI) shows a diffuse lissencephaly-pachygyria spectrum with main affectation at posterior brain areas (Fig. 1).
Genetic and evolutionary analysis
Whole genome sequencing (WGS) of the patient was performed and sequence data was analyzed following a bioinformatic pipeline which included analysis of the quality of reads [12], mapping onto human reference genome (hg19) [13], mark of duplicates, sorting and variant calling [14]. The variants obtained were annotated [15] and then subjected to different sets of filters to detect potentially causative mutations (see Supplementary Material). After these filters were applied, we obtained 40 homozygous or hemizygous variants with population frequency less than 1% located at splicing sites or coding regions, 458 heterozygous variants with population frequency less than 1% and 439 heterozygous variants with population frequency less than 0.5% and located at splicing sites or coding regions.
Among these prioritized variants we found a potential causative mutation in heterozygous state in the LIS1/PAFAH1B1 gene. This gene was previously associated with the phenotype (LIS), especially with an autosomal dominant mechanism of inheritance. The potential causative mutation found is a frameshift insertion of a single nucleotide in exon 8 (PAFAH1B1:NM_000430:exon8:c.681dupG:p.(Lys228Glufs*28)), that lies between the first 23% to 55% of the protein depending on the transcript, according to the SIFT Indel tool [16]. The frameshift indel was reported as damaging with a confidence score of 0.858, and causing a nonsense mediated decay (NMD) response. It generates a premature stop codon 27 amino acids later, causing the loss of 156 amino acids.
If the mutated gene evaded NMD and led to a final product, this would be a 254 amino acids protein instead of the wild type 410 residues. A crystal structure has been described for LIS1 complexed to the brain cytosolic PAF-AH [17]. The complex shows that LIS1 folds into a beta propeller and interacts as a homodimer with a PAF-AH homodimer. From 14 reported surface interacting residues with PAF-AH, 8 are missing from our patient’s hypothetical protein. We predict that the mutated LIS1 could have self-aggregation tendencies, as the 27 new residues composing the shorter C-terminal region, not only would not allow the correct folding into a complete beta-propeller, but in addition would be highly disordered. As a qualitative indicator for this, the homology model in this C-terminal region has very low quality, in particular for the ‘HRTQRMGTY’ amino acid stretch, according to QMEANDiscO scoring function [18]. Even without aggregation and assuming the protein could fold into a ‘half propeller’, this protein would be unable to productively interact with PAF-AH as well as its additional molecular partners, notably dynein and a number of dynein-associated proteins. Indeed, LIS1 has been described as a molecular hub at a crossroad of several pathways, coupling PAF signaling to dynein regulation [17]. We expect all these functions to be hampered or inexistent in the protein product coded by this allelet.
The unaffected mother was sequenced at the proposed variant position and no mutation was detected. This is considered to support the hypothesis of a de novo mutation in the patient. Unfortunately the father of the patient (also unaffected) was not available for analysis.
Figure 2 A shows an IGV view of the candidate position. 37 reads are covering that location with a good quality. Additionally, the variant was confirmed with Sanger Sequencing in the patient (Fig. 2B, top) and was not seen in the mother (Fig. 2B, bottom)
Variant impact at genome and protein level. A IGV view of the reads mapping onto the gene PAFAH1B1. The gene’s structure is shown on top and the exon, where the mutation is located is marked with an asterisk. The reads are shown in gray and the Insertion is marked in violet. The inserted G is shown on top of the reference genome. The originated premature stop codon is also marked. B Sanger sequencing of the variant in mother (top) and patient (bottom). Mother presents no variant, while patient variant was confirmed. C Top view of the truncated LIS1 beta propeller in orange. The C-terminal region of 27 residues introduced between the frameshift and the stop codon is painted in red. For reference, the wild-type crystal structure of LIS1 (PDB access code 1VYH) is shown superimposed as a transparent grey cartoon. Also, eight crucial interacting residues are indicated in atom sticks, completely absent from the hypothetical truncated LIS1
Discussion and conclusion
We found a novel probably causative frameshift variant in a patient with a previously undiagnosed rare disease using WGS. Previous genetic tests (sequencing of MECP2 and ARX genes, and methylation analysis for Angelman syndrome) were performed with inconclusive results. This is expected since Lissencephaly and epileptic encephalopathy are highly heterogeneous genetic disorders in their etiology: ie. different genes are associated with several presentations of this pathology. For example, RELN gene is affected in the Norman-Roberts syndrome (LIS2) [19], heterozygous mutations in TUBA1A are responsible for the LIS3 syndrome [20], homozygous mutations in the NDE1 gene are associated to LIS4 [21], among many others (LAMB1 to LIS5 [22], KATNB1 to LIS6 [23], CDK5 to LIS7 [24], TMTC to LIS8 [25], MACF1 LIS9 [26], CEP85L LIS10 [27]). Additionally X-linked forms of Lissencephaly are caused by DCX and ARX genes [28]. Hence, usually WGS or WES are accurate strategies for assessing patients with epileptic encephalopathy. However, in this case if we had had the MRI results (fairly consistent with LIS) before we had done the NGS sequencing, we might have end up doing a targeted sequencing approach, such as the PAFAH1B1 gene or at least a subset of genes or WES, instead of doing the complete genome. This being a matter of costs and resources and not crucial for the patient’s diagnosis.
The variant we detected was an insertion of one nucleotide (G) in the coding sequence of LIS1 gene, causing a change in the reading frame. The localization of the variant corresponds to the first 23% to 55% of the protein (depending on the transcript) and as a consequence a premature stop codon is produced causing the loss of the last 156 amino acids of the protein. Therefore, a severe affectation of the protein function is expected and probably a silencing of this allele either by encoding a truncated protein or by the mechanism of degradation of messenger RNA mediated by terminator mutations (NMD).
This variant has not been previously described and does not appear in the population frequency databases. It corresponds to the phenotype of the patient and the loss of function in the gene product is a mechanism already described in this condition (truncating mutations were described in other patients and being the gene involved in the microdeletion of Miller-Dieker lissencephaly syndrome). Since there are no other variants in the LIS1/PAFAH1B1 gene with less than 1% population frequency and due to family history, we proposed a de novo mechanism for this case. This was (partially) confirmed by Sanger sequencing of the mother who doesn’t have the mutation. Father was unavailable for analysis, so this aspect remains unknown. Nevertheless, we consider that there is sufficient evidence that supports the pathogenic classification of the novel variant.
According to ACMG (American College of Medical Genetics and Genomics) variant interpretation guidelines [29] the frameshift variant found corresponds to the PVS1 (pathogenicity very strong) rule. It is a null variant (frameshift) in a gene where loss of function (LOF) is a known mechanism of disease (in ExAC database PAFAH1B1 gene has a maximal probability of being LOF intolerant, pLI = 1 [30]). The frameshift mutation is also classified as a PM2 (moderate evidence of pathogenicity) since it was absent in population databases (1000 Genomes Project, GnomAD, etc.). We also consider applying rule PP4 (Patient’s phenotype is highly specific for a disease with a single genetic etiology), since the MRI findings are very specific for PAFAH1B1-related LIS. According to ACMG rules, the variant is classified as pathogenic, as it belongs to one very strong (PVS1), one moderate (PM2) and one supporting category (PP4).
We could also consider applying PP3 (supporting evidence of pathogenicity) since the pathogenic computational verdict is based on one pathogenic prediction from SIFT Indel Tool, one pathogenic prediction from the conservation score GERP [31] and no benign predictions. However, some studies avoid [32] applying PP3 in LoF variants when PVS1 is valid.
Furthermore there are other disruptive (frameshift or stop codon) variants in the same gene region reported as pathogenic [33], supporting the importance of the region for proper gene product function.
Through this work we were able to find a molecular diagnosis of a rare disease in a seven years old boy with severe and heterogeneous neurological symptoms. We found a de novo novel frameshift mutation in the LIS1/PAFAH1B1 gene that most likely causes a silencing of one allele. This finding shows the benefit of the use of NGS as a diagnosis tool of rare diseases.
Availability of data and materials
The variant found in this study is available in ClinVar, ID: SUB9799916 [11]. The bam files generated and analyzed during the current study are only available from the corresponding author upon reasonable request.
Abbreviations
- LIS:
-
Lissencephaly
- MCD:
-
Malformations of cortical development
- PAFAH1B1:
-
Platelet-activating factor acetylhydrolase brain isoform Ib
- PAF:
-
Platelet activating factor
- PAF-AH:
-
Platelet-activating factor acetylhydrolase
- WGS:
-
Whole Genome Sequencing
- GMFSC:
-
Gross Motor Function Classification System
- BFMF:
-
Bimanual fine motor function
- CFCS:
-
Communication Function Classification System
- MRI:
-
Magnetic Resonance Imaging
- ARX:
-
Aristaless Related Homeobox
- MECP2:
-
Methyl-GpG binding protein 2
- NGS:
-
Next Generation Sequencing
- NMD:
-
Nonsense Mediated Decay
References
Jones KL. Smith’s recognizable patterns of human malformation. 6th ed. Philadelphia: Elsevier Saunders; 2006.
Dobyns WB. Developmental aspects of lissencephaly and the lissencephaly syndromes. Birth Defects Orig Artic Ser. 1987;23(1):225–41.
Gleeson JG, Allen KM, Fox JW, Lamperti ED, Berkovic S, et al. Doublecortin, a brain-specific gene mutated in human X-linked lissencephaly and double cortex syndrome, encodes a putative signaling protein. Cell. 1998;92(1):63–72.
Reiner O, Carrozzo R, Shen Y, Wehnert M, Faustinella F, et al. Isolation of a Miller-Dicker lissencephaly gene containing G protein β-subunit-like repeats. Nature. 1993;364(6439):717–21.
Di Donato N, Chiari S, Mirzaa GM, Aldinger K, Parrini E, et al. Lissencephaly: expanded imaging and clinical classification. Am J Med Genet A. 2017;173(6):1473–88.
Cardoso C, Leventer RJ, Dowling JJ, Ward HL, Chung J, et al. Clinical and molecular basis of classical lissencephaly: mutations in theLIS1 gene (PAFAH1B1). Hum Mutat. 2001;19(1):4–15.
Hattori M, Adachi H, Tsujimoto M, Arai H, Inoue K. Miller-Dieker lissencephaly gene encodes a subunit of brain platelet-activating factor acetylhydrolase. Nature. 1994;370:216–8.
Hattori M, Aoki J, Arai H, Inoue K. PAF and PAF acetylhydrolase in the nervous system. J Lipid Mediators Cell Signal. 1996;14:99–102.
Cardoso C, Leventer RJ, Matsumoto N, Kuc JA, Ramocki MB, et al. The location and type of mutation predict malformation severity in isolated lissencephaly caused by abnormalities within the LIS1 gene. Hum Mol Genet. 2000;9(20):3019–28.
Landrum MJ, Lee JM, Benson M, Brown G, Chao C, Chitipiralla S, Maglott DR. ClinVar: public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016;44(D1):D862–8.
VCV001162266.1 - ClinVar. https://www.ncbi.nlm.nih.gov/clinvar/variation/1162266/#id_first. Accessed Jun 2021.
Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164–e164.
Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31(13):3812–4.
Tarricone C, Perrina F, Monzani S, Massimiliano L, Kim MH, et al. Coupling PAF signaling to dynein regulation: structure of LIS1 in complex with PAF-acetylhydrolase. Neuron. 2004;44(5):809–21.
Studer G, Rempfer C, Waterhouse AM, Gumienny R, Haas J, Schwede T. QMEANDisCo—distance constraints applied on model quality estimation. Bioinformatics. 2020;36(6):1765–71.
Hong SE, Shugart YY, Huang DT, Al Shahwan S, Grant PE, et al. Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with human RELN mutations. Nat Genet. 2000;26(1):93–6.
Poirier K, Keays DA, Francis F, Saillour Y, Bahi N, et al. Large spectrum of lissencephaly and pachygyria phenotypes resulting from de novo missense mutations in tubulin alpha 1A (TUBA1A). Hum Mutat. 2007;28(11):1055–64.
Alkuraya FS, Cai X, Emery C, Mochida GH, Al-Dosari MS, et al. Human mutations in NDE1 cause extreme microcephaly with lissencephaly. Am J Hum Genet. 2011;88(5):536–47.
Radmanesh F, Caglayan AO, Silhavy JL, Yilmaz C, Cantagrel V, et al. Mutations in LAMB1 cause cobblestone brain malformation without muscular or ocular abnormalities. Am J Hum Genet. 2013;92(3):468–74.
Mishra-Gorur K, Çağlayan AO, Schaffer AE, Chabu C, Henegariu O, et al. Mutations in KATNB1 cause complex cerebral malformations by disrupting asymmetrically dividing neural progenitors. Neuron. 2014;84(6):1226–39.
Magen D, Ofir A, Berger L, Goldsher D, Eran A, et al. Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with a loss-of-function mutation in CDK5. Hum Genet. 2015;134(3):305–14.
Jerber J, Zaki MS, Al-Aama JY, Rosti RO, Ben-Omran T, et al. Biallelic mutations in TMTC3, encoding a transmembrane and TPR-containing protein, lead to cobblestone lissencephaly. Am J Hum Genet. 2016;99(5):1181–9.
Dobyns WB, Aldinger KA, Ishak GE, Mirzaa GM, Timms AE, et al. MACF1 mutations encoding highly conserved zinc-binding residues of the GAR domain cause defects in neuronal migration and axon guidance. Am J Hum Genet. 2018;103(6):1009–21.
Tsai MH, Muir AM, Wang WJ, Kang YN, Yang KC, et al. Pathogenic variants in CEP85L cause sporadic and familial posterior predominant lissencephaly. Neuron. 2020;106(2):237–45.
Kato M, Dobyns WB. Lissencephaly and the molecular basis of neuronal migration. Human Mol Genet. 2003;12(suppl_1):R89–96.
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(5):405–23.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536(7616):285–91.
Cooper GM, Stone EA, Asimenos G, Green ED, Batzoglou S, et al. Distribution and intensity of constraint in mammalian genomic sequence. Genome Res. 2005;15(7):901–13.
Abou Tayoun AN, Pesaran T, DiStefano MT, et al. Recommendations for interpreting the loss of function PVS1 ACMG/AMP variant criterion. Hum Mutat. 2018;39(11):1517–24. https://doi.org/10.1002/humu.23626.
VCV000021181.4 - ClinVar. https://www.ncbi.nlm.nih.gov/clinvar/variation/21181/#id_first. Accessed 20 Dec 2020.
Acknowledgements
We are very thankful to the patient and his family for kindly accepting to participate in the study.
Funding
This project was funded by URUGENOMES Project IP011-17/CEI/LC/MB. This funding supported the sequencing of the patients and the capacitation of human resources to analyze the data. The funding body was not involved in any decision regarding experimental design, patient selection nor analysis.
Author information
Authors and Affiliations
Contributions
CS: Bioinformatic Analysis (Methodology), data analysis and interpretation, wrote the manuscript. MG: Evolutionary and protein structure prediction analysis. Read and edit the manuscript. SR: Genetic diagnosis, data interpretation, read the manuscript. FBY: Clinical assessment of the patient, data interpretation, read the manuscript. AT: Genetic diagnosis, data interpretation, read the manuscript. ND: Genetic diagnosis, data interpretation, read the manuscript. HN: Supervision, funding, project administration, conceptualization, reviewed the manuscript. VR: Clinical evaluation and diagnosis, supervision of patients, data collection, data interpretation, reviewed the manuscript. LS: Bioinformatic Analysis (Methodology), supervision, conceptualization, coordination of data collection and interpretation, wrote the manuscript. All authors have read and approved the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The patient and his family accepted to participate in the study, and granted permission for the publication of the results involved. This study was approved by the ethics committee of Institut Pasteur de Montevideo (IP011-17/CEI/LC/MB).
Consent for publication
The patient and his family granted permission for the publication of the results involved. Written informed consent was obtained from the patient’s family to the submission of the case report.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
12887_2022_3595_MOESM1_ESM.pdf
Additional file 1.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Simoes, C., Graña, M., Rodriguez, S. et al. Novel frameshift mutation in LIS1 gene is a probable cause of lissencephaly: a case report. BMC Pediatr 22, 545 (2022). https://doi.org/10.1186/s12887-022-03595-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12887-022-03595-6
Keywords
- Lissencephaly
- PAFAH1B1
- Whole-genome sequencing
- Rare disease
- Novel mutation
- Case report