Recessive LAMA5 Variants Associated With Partial Epilepsy and Spasms in Infancy

Objective: The LAMA5 gene encodes the laminin subunit α5, the most abundant laminin α subunit in the human brain. It forms heterotrimers with the subunit β1/β2 and γ1/γ3 and regulates neurodevelopmental processes. Genes encoding subunits of the laminin heterotrimers containing subunit α5 have been reported to be associated with human diseases. Among LAMAs encoding the laminin α subunit, LAMA1-4 have also been reported to be associated with human disease. In this study, we investigated the association between LAMA5 and epilepsy.

Methods: Trios-based whole-exome sequencing was performed in a cohort of 118 infants suffering from focal seizures with or without spasms. Protein modeling was used to assess the damaging effects of variations. The LAMAs expression was analyzed with data from the GTEX and VarCards databases.

Results: Six pairs of compound heterozygous missense variants in LAMA5 were identified in six unrelated patients. All affected individuals suffered from focal seizures with mild developmental delay, and three patients presented also spasms. These variants had no or low allele frequencies in controls and presented statistically higher frequency in the case cohort than in controls. The recessive burden analysis showed that recessive LAMA5 variants identified in this cohort were significantly more than the expected number in the East Asian population. Protein modeling showed that at least one variant in each pair of biallelic variants affected hydrogen bonds with surrounding amino acids. Among the biallelic variants in cases with only focal seizures, two variants of each pair were located in different structural domains or domains/links, whereas in the cases with spasms, the biallelic variants were constituted by two variants in the identical functional domains or both with hydrogen bond changes.

Conclusion: Recessive LAMA5 variants were potentially associated with infant epilepsy. The establishment of the association between LAMA5 and epilepsy will facilitate the genetic diagnosis and management in patients with infant epilepsy.

Introduction

The LAMA5 gene (OMIM* 601033) encodes the laminin subunit α5, expressed in the human brain, especially in the cortex and during the early stages of life (Durkin et al., 1997). The laminin subunit α5 forms heterotrimers with subunits β1/β2 and γ1/γ3 and regulates neurodevelopmental biological processes, including epiblast polarization, neurite outgrowth, neuronal migration, synaptic stability, and cell adhesion, differentiation, migration, and signaling (Liang and Crutcher, 1992; Luckenbill-Edds, 1997; Domogatskaya et al., 2012; Mukherjee et al., 2020). In mice, homozygous knock-out of LAMA5 caused lethality throughout fetal growth and development, and led to exencephaly, megalencephaly, and other neural tube defects (Miner et al., 1998; Kikkawa and Miner, 2006). Genes encoding subunits of the laminin heterotrimers containing subunit α5 have been reported to be associated with human diseases. Among LAMAs encoding the laminin α subunit, LAMA1-4 have also been reported to be associated with human disease. The association between the LAMA5 gene, which encodes the most abundant laminin α subunit in the human brain, and human diseases, has not been determined.

Epilepsy is one of the most common neurological disorders in children with an estimated prevalence of 4–5 per 10,000. Infancy is the critical period of brain development, and epilepsy presents the highest incidence in infancy (Wirrell et al., 2011; Eltze et al., 2013; Hirose, 2013). Multiple seizures may appear in infancy epilepsy, such as focal (partial), myoclonic, spasms seizures, and spasms typically occur in infancy (Bayat et al., 2021). Clinically, both spasms and focal seizures are common in infancy. A proportion of infants with epilepsy have acquired causes, such as trauma, infection, and immune, but the etiologies in the majority are unknown. Previous studies have shown that genetic factors play an important role in the etiology of infant epilepsy (Gan et al., 2019; Kang et al., 2019). The established causative genes include PRRT2, KCNQ2, SCN1A, SCN2A, STXBP1, CDKL5, and ARX, which contribute to approximately 19% of patients with infant-onset epilepsy (Song et al., 2021).

In this study, trios-based whole-exome sequencing (WES) was performed in a cohort of infants with epilepsy. Six pairs of compound heterozygous missense variants in LAMA5 were detected in six unrelated cases. The present study suggests that recessive LAMA5 variants were potentially associated with infant epilepsy.

Materials and Methods

Subjects

One hundred and eighteen infants who suffered from focal seizures without any acquired causes were recruited from four hospitals, including the Second Affiliated Hospital of Guangzhou Medical University, Foshan Maternal and Child Health Hospital, the Second Affiliated Hospital of Shandong University, and Children’s Hospital of Shandong University, from June 2019 to July 2021. Clinical information of the affected individuals was collected, including age at onset, type and frequency of seizures, family history, systemic and neurological findings, and effective antiepileptic drugs. The structural abnormalities were detected by brain magnetic resonance imaging (MRI) scans. Long-term (24 h) video electroencephalography (EEG) monitoring records were performed with electrodes being arranged according to the international standard of 10–20 reduced montage system. The procedures of open-close eyes test, hyperventilation, intermittent photic stimulation, and sleeping recording were obtained. The EEG results were reviewed by at least two qualified electroencephalographers. The Chinese version of the Gesell development scales was utilized in the neurodevelopment evaluation of the participants according to their ages. Its scores were evaluated using tests for gross motor, fine motor, adaptive behavior, language, and personal-social behavior. Epileptic seizures and epilepsy syndromes were diagnosed according to the criteria of the Commission on Classification and Terminology of the International League Against Epilepsy (1981, 1989, 2001, 2010, and 2017). All patients were followed up for more than 1 year.

For the controls, WES was performed on 296 healthy Chinese volunteers who served as a normal control group as our previous report (Wang et al., 2018). Frequencies of the identified variants were also compared with that in the other control populations, including East Asian and general populations in the Genome Aggregation Database (gnomAD¹), and the 8,364 persons without known neuropsychiatric conditions in the Epi25 WES Browser² (Epi25 Collaborative, 2019; Karczewski et al., 2020).

This study adhered to the principles of the International Committee of Medical Journal Editors concerning patient consent for research or participation and received approval from the Ethics Committee of the Second Affiliated Hospital of Guangzhou Medical University. Written informed consents were provided by the patient’s legal guardians.

Whole Exon Sequencing

Blood samples of the probands, their parents, and other available family members were collected to identify the source of the genetic variation and to aid in the analysis of the pathogenicity of variants. Genomic DNAs were extracted from blood samples using the Qiagen Flexi Gene DNA kit (Qiagen, Hilden, Germany). WES was performed using a NextSeq500 sequencing instrument (Illumina, San Diego, CA, United States) following the standard procedures previously described (Wang et al., 2018). The sequencing data were generated by massively parallel sequencing with an average depth of >125× and >98% coverage of the capture region on the chip for obtaining high-quality reads that were mapped to the Genome Reference Consortium Human genome build 37 by Burrows–Wheeler alignment. Single-nucleotide point variants and indels were called with the Genome Analysis Toolkit.

Genetic Analysis

A case-by-case analytical approach was adopted to identify candidate causative variants in each trio. Primarily, the rare variants were prioritized with a minor allele frequency <0.005 in the Genome Aggregation Database (see Text Footnote 1). Then, potentially pathogenic variants were retained, including frameshift, nonsense, canonical splice site, initiation codon, and missense variants predicted as being damaging in silico tools. Lastly and importantly, the potentially disease-causing variants in each trio were analyzed with an individualized protocol. The variants in each trio were sorted according to the following five models: (1) epilepsy-associated gene model (Wang et al., 2017); (2) de novo variant dominant model; (3) autosomal recessive inheritance model, including homozygous and compound heterozygous variants; (4) X-linked model; (5) co-segregation analysis model. To identify novel epilepsy-associated genes, the known epilepsy-associated genes were put aside. Genes with de novo variants, bi-allelic variants, hemizygous variants, or variants with segregations, which represent the genetic difference between the patients and normal individuals in a family and potentially explain the occurrence of disease, were selected for further studies to define the gene-disease association. LAMA5 emerged as one of the candidate genes with recurrent bi-allelic variants in this cohort of infancy epilepsy. Sanger sequencing was used to validate the candidate variants. All LAMA5 variants identified in this study were annotated into the reference transcript NM_005560.4.

Bioinformatic Analyses

Protein modeling was performed by using the Iterative Threading ASSEmbly Refinement software (I-TASSER³) to evaluate the damaging effect of candidate variants (Yang and Zhang, 2015). The confidence of each model was quantitatively measured by a C-score in the range of [−5, 2]. PyMOL Molecular Graphics System (Version 2.3.2; Schrödinger, LLC; New York, NY, United States) was used for three-dimensional protein structure visualization and analysis. I-Mutant Suite server was used for the prediction of protein stability changes upon single-nucleotide variants that lead to changes in the amino acid⁴ (Capriotti et al., 2005). The alteration of the protein stability was assessed by the free energy change value (ΔΔG, kcal/mol). Values greater than 0.5 kacl/mol imply a large increase in protein stability, values less than −0.5 kacl/mol are considered to be a large decrease in protein stability, and others imply neutral stability. The consequences of all the missense variants were predicted by several common tools, including SIFT, PolyPhen2_HVAR, CADD, MutationTaster, GenoCanyon, fitCons, M_CAP, and GERP.

Statistical Analysis

R (version 4.0.3) was used for data processing. The frequencies of the LAMA5 variants between the epilepsy cohort and the controls were compared by a two-sided Fisher exact test (CONVERGE Consortium, 2015). The burden of recessive variants was analyzed according to the method recommended recently (Martin et al., 2018). P-value < 0.05 was considered statistically significant.

Results

Identification of LAMA5 Variants

Six pairs of compound heterozygous missense variants, i.e., (c.1337G > A/p.Arg446Gln and c.10699C > T/p.Pro3567Ser), (c.1418G > A/p.Pro473Leu and c.3608C > T/p.Arg1203Gln), (c.5426G > A/p.Arg1809His and c.7394G > A/p.Arg2465Gln), (c.3170C > T/p.Ser1057Leu and c.6388C > T/p.Arg2130Cys), (c.9448G > A/p.Gly3150Ser and c.10744C > T/p.Arg3582Trp), and (c.1963G > A/p.Gly655Ser and c.2192C > G/p.Ala731Gly), were identified in six unrelated cases with focal seizures, including three infants with also spasms (Figures 1A,B and Table 1). The compound heterozygous variants originated from their asymptomatic mothers and fathers, consistent with a classical recessive inheritance pattern.

FIGURE 1

Figure 1. Genetic data of cases with LAMA5 variants. (A) Pedigrees of the six cases with LAMA5 mutations and their corresponding phenotypes. (B) DNA sequence chromatogram of the LAMA5 mutations. Arrows indicate the positions of the variants.

TABLE 1

Table 1. Clinical features of the patients with LAMA5 variants.

These variants presented no or low frequencies (MAF < 0.005) in the gnomAD databases (Table 2). Nine of the variants did not present in the normal control of Epi25 WES Brower, and the other three variants presented extremely low frequencies (MAF < 0.0005). None of the variants, except p.Ala731Gly, were presented in the 296 normal controls.

TABLE 2

Table 2. Analysis of the aggregate frequency of LAMA5 variants identified in this study.

When the burden of recessive variants was analyzed, (Martin et al., 2018) the LAMA5 variants in the present cohort were significantly more than the expected number by chance in the East Asian population (MAF < 0.005, P = 6.4654 × 10^–6). Furthermore, the aggregate frequency of the variants in this cohort was significantly higher than that in the six controls (Table 2), including the controls of 296 normal individuals (12/236 vs. 1/592, p = 2.243 × 10^–6), the normal control of Epi25 WES Browser (vs. 6/16844, p < 2.2 × 10^–16), the gnomAD-all population (vs. 472/248,058; p = 8.306 × 10^–14), the controls of gnomAD-all population (vs. 218/100,060; p = 4.728 × 10^–13), the gnomAD-East Asian population (vs. 167/18250, p = 3.36 × 10^–6), and the controls of the gnomAD-East Asian population (vs. 71/8836, p = 1.749 × 10^–6).

All LAMA5 variants identified in this study were predicted to be damaging by at least two in silico tools (Supplementary Table 1). The probability that transcript falls into distribution of recessive genes (pRec) is 0.99353 for the LAMA5 gene, indicating that it is very likely intolerant to recessive loss-of-function variations (Karczewski et al., 2020). None of the six affected individuals had pathogenic or likely pathogenic variants in genes known to be associated with epilepsy (Wang et al., 2017).

Clinical Features of the Cases With LAMA5 Variants

The main clinical features of the six patients with LAMA5 variants were summarized in Table 1. All patients showed infant-onset epilepsy with onset ages ranging from 1 to 9 months. The patients of cases 1, case 2, and case 3 suffered from infrequent focal seizures or secondary tonic-clonic seizures (monthly or yearly) and became seizure-free after treatment with valproate monotherapy or valproate in combination with levetiracetam. Their EEGs showed bilateral, unilateral, and multiple discharges, predominantly at the frontal, central, and temporal lobe, mainly during sleep (Figure 2A). The patients of case 4 and case 5 had daily focal seizures with spasms. Seizure-free was achieved after treatment with a combination of adrenocorticotropic hormones and topiramate. The patient of case 6 also had focal seizures and spasms (daily). The spasms and focal seizures disappeared after treatment, but myoclonic seizures appeared, which were infrequent (3–4 times/month). Interictal multifocal and generalized discharges were recorded in the three cases with spasms (Figure 2B).

FIGURE 2

Figure 2. Ictal and interictal EEG in the cases with LAMA5 variants. (A) Interictal EEG in case 1 at the age of 8 months showed irregular spike-slow waves in the middle and posterior temporal region. (B) Interictal EEG in case 4 at the age of 3 months showed diffuse high amplitude spike and spike-slow waves, and focal spike-slow waves in the anterior region.

These patients were full-term and delivered without abnormalities. Brain MRI was normal in the six cases. All patients showed mild global developmental delays.

Molecular Alteration of Laminin Subunit α5

The laminin subunit α5 contains one signal peptide, one Laminin N terminal domain, 22 Laminin EGF-like domains, one Laminin IV type A domain, and five Laminin G-like domains (Uniport-id: O15230) (Figure 3A). Five of the variants were located in Laminin EGF-like domains, including p.Arg446Gln, p.Pro473Leu, p.Gly655Ser, p.Ala731Gly, and p.Arg2130Cys. Three variants, p.Glu3150Ser, p.Pro3567Ser, and p.Arg3582Trp, were located in the Laminin G-like domains. Variant p.Val1809His was located in the Laminin IV type A domain. The other three variants p.Ser1057Leu, p.Arg1203Gln, and p.Arg2465Gln are located between structural domains (Figure 3A).

FIGURE 3

Figure 3. Schematic illustration of LAMA5 variants. (A) Schematic illustration of the laminin subunit α5 and the location of the LAMA5 variants identified in this study. Two variants of the same height represent a pair of biallelic mutations. The variants of red-colored represent cases with both spasms and focal seizures; those of blue-colored represent cases with only focal seizures. (B) Hydrogen bond changes and free energy stability changes (ΔΔG, Kcal/mol) value of the variants from the present study.

The molecular effects of the variants were analyzed by using I-TASSER for protein modeling and PyMOL for visualization. Eight of the variants, including p.Arg446Gln, p.Pro3567Ser, p.Arg1203Gln, p.Arg2465Gln, p.Ser1057Leu, p.Arg2130Cys, p.Arg3582Trp, and p.Ala731Gly changed their hydrogen bonds with surrounding residues. The other four variants, p.Pro473Leu, p.Gly655Ser, p.Arg1809His, and p.Gly3150Ser, did not change their hydrogen bonds but were predicted to decrease the protein stability significantly with ΔΔG values less than −0.5 kcal/mol. In each pair of compound heterozygous variants, at least one variant had hydrogen bonds change (Figure 3B).

It is notable that among the compound heterozygous variants in the three cases with only focal seizures (case 1, case 2, and case 3), two variants of each pair were located in different structural domains or domains/links (Figure 3A). Case 1 had the compound heterozygous variants with two variants located furthest apart (p.Arg446Gln and p.Pro3567Ser); the patient showed a milder phenotype than the others, e.g., the latest onset age (8-month-old) and infrequent seizure frequency (yearly). In the three cases with spasms, two pairs of compound heterozygous variants (case 5 and case 6) were constituted by two variants in the identical functional domains. The two variants of case 4 were in non-identical structural domains but presented the most pronounced hydrogen bonding changes (a total of eight hydrogen bonds disrupted).

The Expression Profile of LAMAs

The functional laminins are cruciform heterotrimers that consist of three short arms formed by the N-terminal portion of α, β, and γ subunits, respectively, and a long arm polymerized by the C-terminal parts of the three subunits (Figure 4A). The laminin short arms (N-terminus) are involved in laminin’s ability to polymerize with those of other laminins to form a polygonal network, while the laminin long arm interacts with the extracellular matrix components to regulate biological processes (Hohenester, 2019). Laminin subunit α5 appears in laminin-511 (with β1 and γ1), laminin-521 (with β2 and γ1), and laminin-523 (with β2 and γ3). Previous studies show that laminins containing subunit α5 play an essential role in embryonic development and are intensely expressed on the surface of the ectoderm (Copp et al., 2011). Tissue-specific expression is the basis of gene function and subsequently the clinical phenotype. We thus compared the expression of LAMAs in the human brain from the data in VarCards and Genotype-Tissue Expression databases (GTEx Consortium, 2013). The LAMA5 gene presented the highest expression in the human brain (Figure 4B). Furthermore, it was more abundant in the cerebral cortex, substantia nigra, frontal cortex, hippocampus, and anterior cingulate cortex (Figure 4C).

FIGURE 4

Figure 4. Schematic illustration of laminin heterotrimers and the expression profile of LAMAs. (A) Schematic illustration of the heterotrimers containing laminin subunit α5. Laminin subunit α5 formed laminin-511 with subunit β1 and γ1, laminin-521 with subunit β2 and γ1, and laminin-523 with subunit β2 and γ3. (B) The overall expression of LAMAs in the human brain. RPKM, Reads Per Kilobase per Million mapped reads. (C) Heatmap and hierarchical clustering of LAMAs expression in the sub-regions of human brain. Columns represent individual sub-region and rows represent individual genes. The darker the color, the higher the expression. The lines represent cluster analysis. More details were presented in the GTEx database (www.gtexportal.org).

Discussion

The LAMA5 gene, predominantly expressed in the early stage of life, plays a vital role in neurodevelopmental biological processes such as neurite outgrowth, epiblast polarization, neuronal migration, synaptic stability, and cell adhesion, differentiation, migration, and signaling. A functional study suggested that the laminin subunit α5 was deposited at synapses in the brain and was essential for dendritic spine structural regulation and synapse stability (Omar et al., 2017). Experimentally, mutant protein of LAMA5 resulted in decreased binding to the synaptic vesicle protein encoded by SV2A (Maselli et al., 2018). Abnormalities in the central nervous system were monitored in the conditional allele knockout LAMA5 mice (Omar et al., 2017). In this study, six pairs of biallelic variants in LAMA5 were identified in patients with infant-onset epilepsy. All variants presented no or low allele frequencies in controls. The aggregate frequency of the LAMA5 variants identified in the case cohort was significantly higher than that in controls. The recessive burden analysis also showed that the LAMA5 variants in the present cohort were significantly more than the expected number in the East Asian population. These findings suggest a potential association between LAMA5 variants and epilepsy.

Previously, a de novo canonical splice site variant (c.10828 + 1G > A) was identified in a patient with developmental delays and epilepsy (Han et al., 2018). In a patient with cortical developmental malformations, a pair of compound heterozygous variants (p.Glu3567Lys and p.Asp2372Asn) were identified (Wiszniewski et al., 2018). One pair of homozygous variants in LAMA5 (p.Arg2659Trp) and one homozygous variant in LAMA1 were identified in a patient with the presynaptic congenital myasthenic syndrome who has non-classified facial tics or tics (Maselli et al., 2018). A heterozygous variant in LAMA5 (p.Val3140Met) was identified in a family with a complex multisystem syndrome, in which several family members with the LAMA5 variant had seizures (Sampaolo et al., 2017). These data provided clues for the possible association between LAMA5 and epilepsy, but the pathogenic role of LAMA5 variants in epilepsy could not be determined due to the variable phenotype, co-appearance of variants in other potentially pathogenic genes, or the single affected case. In the present study, we identified LAMA5 variants in six unrelated cases, and the other possible pathogenic genes were excluded. Therefore, this study provided direct evidence in supporting the association between LAMA5 and epilepsy.

In animals, the homozygous LAMA5 knockout mice model showed preweaning lethality with complete penetrance⁵, while the heterozygous knockout mice showed fewer abnormalities, suggesting a dose-effect. The pRec, a metric for intolerant of biallelic loss-of-function variants, is 0.99353 for the LAMA5 gene, indicating that it is highly intolerant for recessive loss-of-function variations (Karczewski et al., 2020). The LAMA5 variants identified in the present study were all biallelic variants. It was noted that among the compound heterozygous variants in the three cases with only focal seizures, two variants of each pair variant were located in different structural domains or domains/links. In the three cases with spasms, two pairs of variants were constituted by two variants in the identical functional domains, and another pair of variants presented the most pronounced hydrogen bonding changes. Taken together that homozygous LAMA5 knockout mice displayed preweaning lethality with complete penetrance, it was considered that homozygous variants in LAMA5 would cause the most severe damage effect and resulted in the most severe phenotype or even early death; the compound heterozygous variants with two variants located in identical functional domains potentially created relatively severe damage effect and caused severe epilepsy such as spasms; and the compound heterozygous variants with two variants located in non-identical functional domains potentially caused relatively mild damage effects and led to mild phenotypes. It is suggested that the location of variants in compound heterozygous variants was potentially associated with the phenotype severity, providing one of the explanations for the phenotype variation.

Structurally, the laminin subunit α5 forms heterotrimers with laminin subunits β1/β2 and γ1/γ3, which were encoded by LAMB1, LAMB2, LAMC1, and LAMC3 that were associated with neurodevelopmental diseases and epilepsy. The LAMB1 gene is the responsible gene of lissencephaly 5 (OMIM #615191) characterized by focal and spasmodic seizures and psychomotor development delay (Radmanesh et al., 2013; Tonduti et al., 2015). The LAMB2 gene is the causative gene of Pierson syndrome (OMIM #609049), with which the majority of patients died early and the survivors presented severe neurodevelopmental delays (Pierson et al., 1963; Zenker et al., 2004). The LAMC1 gene has been repeatedly reported in Dandy–Walker malformation with occipital cephalocele, and most of the affected individuals had an infant-onset intellectual disability with or without seizures (Carvalho et al., 2006; Darbro et al., 2013). Variants in LAMC3 were the causes of occipital cortical malformations (OMIM #614115), and the affected individuals experienced seizures (Barak et al., 2011). Experimentally, these genes caused neurological abnormalities through a common mechanism of disrupted function of laminin heterotrimers in regulating neuronal migration and other biological processes (Pierson et al., 1963; Zenker et al., 2004; Carvalho et al., 2006; Darbro et al., 2013; Radmanesh et al., 2013; Tonduti et al., 2015). However, the LAMA5 gene, which encodes one of the indispensable parts of the heterotrimers, has not been confirmed to be associated with neurological disorders. In this study, the patients with biallelic LAMA5 variants presented focal seizures and developmental delays, suggesting a potential role of LAMA5 in epilepsy with the involvement of neurodevelopment.

Regarding LAMAs that encode the laminin α subunits, the associations between LAMA1-4 and human diseases have been established. Variants in LAMA1 cause Poretti–Boltshauser syndrome (OMIM #615960), characterized by delayed motor development, speech delay, and cognitive function; and seizures, tics, and spasticity have also been observed (Aldinger et al., 2014; Elmas et al., 2020). LAMA2 is the causative gene of autosomal recessive limb-girdle muscular dystrophy-23 (OMIM #618138) and congenital merosin deficient or partially deficient muscular dystrophy (OMIM #607855), in which epilepsy was regarded as one of the core features (Chan et al., 2014; Xiong et al., 2015; Salvati et al., 2021). The LAMA3 gene is the responsible gene of Herlitz type junctional epidermolysis bullosa (OMIM #226700), generalized atrophic benign epidermolysis bullosa (OMIM #226650), and laryngo-onycho-cutaneous syndrome (OMIM #245660) (Kivirikko et al., 1995; McGrath et al., 1995; Vidal et al., 1995; Nakano et al., 2002; McLean et al., 2003). The LAMA4 gene is the causative gene of dilated cardiomyopathy 1JJ (OMIM #615235) (Knoll et al., 2007). Compared with the other LAMAs genes, the LAMA5 gene had the highest expression in the human brain and was abundant in the cerebral cortex, substantia nigra, frontal cortex, hippocampus, and anterior cingulate cortex. Generally, tissue-specific expression is the basis of gene function and subsequently the clinical phenotype. The highest expression of LAMA5 in the human brain provided an anatomical basis for the association between LAMA5 and neurological diseases.

The proteins encoded by LAMAs are of tissue-specific expression and independent function, while the protein encoded by LAMA5 interacts with that encoded by LAMBs and LAMCs. Therefore, the possibility of digenic or polygenic mechanisms could not be excluded. No pathogenetic or likely pathogenetic variants in LAMBs, LAMCs, or other epilepsy genes were identified in these patients, suggesting a potential monogenic role of LAMA5. However, the patients with LAMA5 variants presented a relatively moderate phenotype, and several variants presented in the population with low frequency. Thus, a modifier role of LAMA5 could not be excluded. Further studies with large cohorts are needed to validate the pathological nature of LAMA5 variants.

Among the disease-causing genes in humans (OMIM⁶), 1,008 genes were associated with disease in a dominant inheritance pattern, whereas 1,936 genes were in a recessive inheritance pattern. Considering that the genome in humans is diploid, it is possible that recessive variants were more common than dominant variants in the etiology of human diseases. Currently, most infantile spasms-related genes (36/47) (Bayat et al., 2021) and all focal epilepsy-associated genes were of autosomal dominant inheritance. Thus, more attention should be paid to recessive variants in epilepsy. This study revealed LAMA5 as a potential novel autosomal recessive gene in infant epilepsy, enriching the genetic etiology of epilepsy.

This study has several limitations. The direct functional effects of the variants were not examined. The consequences of the variants on the interactions of the alpha subunit with its partners warrant further studies. Besides point variants, whether CNV in the LAMA5 gene was pathogenic should also be considered.

Conclusion

The association between LAMA5 and epilepsy was supported by multiple pieces of evidence, such as common clinical features, unique gene functions, and statistical evidence. The establishment of the association between LAMA5 and epilepsy will facilitate the genetic diagnosis and management in patients with infant epilepsy.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI GenBank: OM994299 – OM994307 and OM994308 – OM994334.

Ethics Statement

This study adhered to the principles of the International Committee of Medical Journal Editors concerning patient consent for research or participation and received approval from the Ethics Committee of the Second Affiliated Hospital of Guangzhou Medical University. Written informed consent was provided by the patients’ legal guardians.

Author Contributions

SL, Z-GL, and W-PL designed the experiments. JuW, J-XL, X-GY, XL, Q-XZ, X-RL, JiW, L-DG, F-LL, B-ML, Z-FG, Q-HG, and Y-HY collected and analyzed the clinical data and patient samples. Z-LY and HL performed the computational modeling. SL, Z-GL, and W-PL wrote the manuscript, with contributions from all the authors.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 81971216 and 81871015), Natural Science Foundation of Guangdong Province (2020A1515010108), Science and Technology Project of Guangzhou (Grant Nos. 201904010292 and 201904020028), Science and Technology Project of Guangdong Province (Grant No. 2017B030314159), Shandong Medical and Health Science and Technology Development Plan (Grant No. 202106010271), and Multi-Center Clinical Research Fund Project of the Second Affiliated Hospital of Guangzhou Medical University (2020-LCYJ-DZX-03).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The help of patients and clinicians participating in this work are greatly appreciated.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2022.825390/full#supplementary-material

Footnotes

1. ^ gnomad.broadinstitute.org
2. ^ https://epi25.broadinstitute.org/
3. ^ https://zhanglab.ccmb.med.umich.edu/I-TASSER/
4. ^ https://folding.biofold.org/i-mutant/i-mutant2.0.html
5. ^ http://www.informatics.jax.org/allele/MGI:1934917
6. ^ www.omim.org

References

Aldinger, K. A., Mosca, S. J., Tetreault, M., Dempsey, J. C., Ishak, G. E., Hartley, T., et al. (2014). Mutations in LAMA1 cause cerebellar dysplasia and cysts with and without retinal dystrophy. Am. J. Hum. Genet. 95, 227–234. doi: 10.1016/j.ajhg.2014.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Barak, T., Kwan, K. Y., Louvi, A., Demirbilek, V., Saygi, S., Tuysuz, B., et al. (2011). Recessive LAMC3 mutations cause malformations of occipital cortical development. Nat. Genet. 43, 590–594. doi: 10.1038/ng.836

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayat, A., Bayat, M., Rubboli, G., and Moller, R. S. (2021). Epilepsy syndromes in the first year of life and usefulness of genetic testing for precision therapy. Genes (Basel) 12:1051. doi: 10.3390/genes12071051

PubMed Abstract | CrossRef Full Text | Google Scholar

Capriotti, E., Fariselli, P., and Casadio, R. (2005). I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 33, W306–W310. doi: 10.1093/nar/gki375

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, D. R., Giuliani, L. R., Simao, G. N., Santos, A. C., and Pina-Neto, J. M. (2006). Autosomal dominant atretic cephalocele with phenotype variability: report of a Brazilian family with six affected in four generations. Am. J. Med. Genet. A 140, 1458–1462. doi: 10.1002/ajmg.a.31255

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, S. H., Foley, A. R., Phadke, R., Mathew, A. A., Pitt, M., Sewry, C., et al. (2014). Limb girdle muscular dystrophy due to LAMA2 mutations: diagnostic difficulties due to associated peripheral neuropathy. Neuromuscul. Disord. 24, 677–683. doi: 10.1016/j.nmd.2014.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

GTEx Consortium (2013). The genotype-tissue expression (GTEx) project. Nat. Genet. 45, 580–585. doi: 10.1038/ng.2653

PubMed Abstract | CrossRef Full Text | Google Scholar

CONVERGE Consortium (2015). Sparse whole-genome sequencing identifies two loci for major depressive disorder. Nature 523, 588–591. doi: 10.1038/nature14659

PubMed Abstract | CrossRef Full Text | Google Scholar

Copp, A. J., Carvalho, R., Wallace, A., Sorokin, L., Sasaki, T., Greene, N. D., et al. (2011). Regional differences in the expression of laminin isoforms during mouse neural tube development. Matrix Biol. 30, 301–309. doi: 10.1016/j.matbio.2011.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Darbro, B. W., Mahajan, V. B., Gakhar, L., Skeie, J. M., Campbell, E., Wu, S., et al. (2013). Mutations in extracellular matrix genes NID1 and LAMC1 cause autosomal dominant Dandy-Walker malformation and occipital cephaloceles. Hum. Mutat. 34, 1075–1079. doi: 10.1002/humu.22351

PubMed Abstract | CrossRef Full Text | Google Scholar

Domogatskaya, A., Rodin, S., and Tryggvason, K. (2012). Functional diversity of laminins. Annu. Rev. Cell Dev. Biol. 28, 523–553. doi: 10.1146/annurev-cellbio-101011-155750

PubMed Abstract | CrossRef Full Text | Google Scholar

Durkin, M. E., Loechel, F., Mattei, M. G., Gilpin, B. J., Albrechtsen, R., and Wewer, U. M. (1997). Tissue-specific expression of the human laminin alpha5-chain, and mapping of the gene to human chromosome 20q13.2-13.3 and to distal mouse chromosome 2 near the locus for the ragged (Ra) mutation. FEBS Lett. 411, 296–300. doi: 10.1016/s0014-5793(97)00686-8

CrossRef Full Text | Google Scholar

Elmas, M., Gogus, B., and Solak, M. (2020). Understanding what you have found: a family with a mutation in the LAMA1 gene with literature review. Clin. Med. Insights Case Rep. 13:1179547620948666. doi: 10.1177/1179547620948666

PubMed Abstract | CrossRef Full Text | Google Scholar

Eltze, C. M., Chong, W. K., Cox, T., Whitney, A., Cortina-Borja, M., Chin, R. F., et al. (2013). A population-based study of newly diagnosed epilepsy in infants. Epilepsia 54, 437–445. doi: 10.1111/epi.12046

PubMed Abstract | CrossRef Full Text | Google Scholar

Epi25 Collaborative (2019). Ultra-rare genetic variation in the epilepsies: a whole-exome sequencing study of 17,606 individuals. Am. J. Hum. Genet. 105, 267–282. doi: 10.1016/j.ajhg.2019.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Gan, J., Cai, Q., Galer, P., Ma, D., Chen, X., Huang, J., et al. (2019). Mapping the knowledge structure and trends of epilepsy genetics over the past decade: a co-word analysis based on medical subject headings terms. Medicine (Baltimore) 98:e16782. doi: 10.1097/MD.0000000000016782

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, J. Y., Jang, J. H., Park, J., and Lee, I. G. (2018). Targeted next-generation sequencing of korean patients with developmental delay and/or intellectual disability. Front. Pediatr. 6:391. doi: 10.3389/fped.2018.00391

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirose, G. (2013). [An overview of epilepsy: its history, classification, pathophysiology and management]. Brain Nerve 65, 509–520.

PubMed Abstract | Google Scholar

Hohenester, E. (2019). Structural biology of laminins. Essays Biochem. 63, 285–295. doi: 10.1042/EBC20180075

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, K. W., Kim, W., Cho, Y. W., Lee, S. K., Jung, K. Y., Shin, W., et al. (2019). Genetic characteristics of non-familial epilepsy. PeerJ 7:e8278. doi: 10.7717/peerj.8278

PubMed Abstract | CrossRef Full Text | Google Scholar

Karczewski, K. J., Francioli, L. C., Tiao, G., Cummings, B. B., Alfoldi, J., Wang, Q., et al. (2020). The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443. doi: 10.1038/s41586-020-2308-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kikkawa, Y., and Miner, J. H. (2006). Molecular dissection of laminin alpha 5 in vivo reveals separable domain-specific roles in embryonic development and kidney function. Dev. Biol. 296, 265–277. doi: 10.1016/j.ydbio.2006.04.463

PubMed Abstract | CrossRef Full Text | Google Scholar

Kivirikko, S., McGrath, J. A., Baudoin, C., Aberdam, D., Ciatti, S., Dunnill, M. G., et al. (1995). A homozygous nonsense mutation in the alpha 3 chain gene of laminin 5 (LAMA3) in lethal (Herlitz) junctional epidermolysis bullosa. Hum. Mol. Genet. 4, 959–962. doi: 10.1093/hmg/4.5.959

PubMed Abstract | CrossRef Full Text | Google Scholar

Knoll, R., Postel, R., Wang, J., Kratzner, R., Hennecke, G., Vacaru, A. M., et al. (2007). Laminin-alpha4 and integrin-linked kinase mutations cause human cardiomyopathy via simultaneous defects in cardiomyocytes and endothelial cells. Circulation 116, 515–525. doi: 10.1161/CIRCULATIONAHA.107.689984

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, S., and Crutcher, K. A. (1992). Neuronal migration on laminin in vitro. Brain Res. Dev. Brain Res. 66, 127–132. doi: 10.1016/0165-3806(92)90148-p

CrossRef Full Text | Google Scholar

Luckenbill-Edds, L. (1997). Laminin and the mechanism of neuronal outgrowth. Brain Res. Brain Res. Rev. 23, 1–27. doi: 10.1016/s0165-0173(96)00013-6

CrossRef Full Text | Google Scholar

Martin, H. C., Jones, W. D., McIntyre, R., Sanchez-Andrade, G., Sanderson, M., Stephenson, J. D., et al. (2018). Quantifying the contribution of recessive coding variation to developmental disorders. Science 362, 1161–1164. doi: 10.1126/science.aar6731

PubMed Abstract | CrossRef Full Text | Google Scholar

Maselli, R. A., Arredondo, J., Vazquez, J., Chong, J. X., Bamshad, M. J., Nickerson, D. A., et al. (2018). A presynaptic congenital myasthenic syndrome attributed to a homozygous sequence variant in LAMA5. Ann. N Y Acad. Sci. 1413, 119–125. doi: 10.1111/nyas.13585

PubMed Abstract | CrossRef Full Text | Google Scholar

McGrath, J. A., Kivirikko, S., Ciatti, S., Moss, C., Dunnill, G. S., Eady, R. A., et al. (1995). A homozygous nonsense mutation in the alpha 3 chain gene of laminin 5 (LAMA3) in Herlitz junctional epidermolysis bullosa: prenatal exclusion in a fetus at risk. Genomics 29, 282–284. doi: 10.1006/geno.1995.1246

PubMed Abstract | CrossRef Full Text | Google Scholar

McLean, W. H., Irvine, A. D., Hamill, K. J., Whittock, N. V., Coleman-Campbell, C. M., Mellerio, J. E., et al. (2003). An unusual N-terminal deletion of the laminin alpha3a isoform leads to the chronic granulation tissue disorder laryngo-onycho-cutaneous syndrome. Hum. Mol. Genet. 12, 2395–2409. doi: 10.1093/hmg/ddg234

PubMed Abstract | CrossRef Full Text | Google Scholar

Miner, J. H., Cunningham, J., and Sanes, J. R. (1998). Roles for laminin in embryogenesis: exencephaly, syndactyly, and placentopathy in mice lacking the laminin alpha5 chain. J. Cell Biol. 143, 1713–1723. doi: 10.1083/jcb.143.6.1713

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukherjee, C., Saleem, S., Das, S., Biswas, S. C., and Bhattacharyya, D. (2020). Human placental laminin: role in neuronal differentiation, cell adhesion and proliferation. J. Biosci. 45:93.

PubMed Abstract | Google Scholar

Nakano, A., Chao, S. C., Pulkkinen, L., Murrell, D., Bruckner-Tuderman, L., Pfendner, E., et al. (2002). Laminin 5 mutations in junctional epidermolysis bullosa: molecular basis of Herlitz vs. non-Herlitz phenotypes. Hum. Genet. 110, 41–51. doi: 10.1007/s00439-001-0630-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Omar, M. H., Kerrisk Campbell, M., Xiao, X., Zhong, Q., Brunken, W. J., Miner, J. H., et al. (2017). CNS neurons deposit laminin alpha5 to stabilize synapses. Cell Rep. 21, 1281–1292. doi: 10.1016/j.celrep.2017.10.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Pierson, M., Cordier, J., Hervouuet, F., and Rauber, G. (1963). [An unusual congenital and familial congenital malformative combination involving the eye and kidney]. J. Genet. Hum. 12, 184–213.

PubMed Abstract | Google Scholar

Radmanesh, F., Caglayan, A. O., Silhavy, J. L., Yilmaz, C., Cantagrel, V., Omar, T., et al. (2013). Mutations in LAMB1 cause cobblestone brain malformation without muscular or ocular abnormalities. Am. J. Hum. Genet. 92, 468–474. doi: 10.1016/j.ajhg.2013.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Salvati, A., Bonaventura, E., Sesso, G., Pasquariello, R., and Sicca, F. (2021). Epilepsy in LAMA2-related muscular dystrophy: a systematic review of the literature. Seizure 91, 425–436. doi: 10.1016/j.seizure.2021.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Sampaolo, S., Napolitano, F., Tirozzi, A., Reccia, M. G., Lombardi, L., Farina, O., et al. (2017). Identification of the first dominant mutation of LAMA5 gene causing a complex multisystem syndrome due to dysfunction of the extracellular matrix. J. Med. Genet. 54, 710–720. doi: 10.1136/jmedgenet-2017-104555

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, T. Y., Deng, J., Fang, F., Chen, C. H., Wang, X. H., Wang, X., et al. (2021). [The etiology of 340 infants with early-onset epilepsy]. Zhonghua Er Ke Za Zhi 59, 387–392. doi: 10.3760/cma.j.cn112140-20201016-00947

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonduti, D., Dorboz, I., Renaldo, F., Masliah-Planchon, J., Elmaleh-Berges, M., Dalens, H., et al. (2015). Cystic leukoencephalopathy with cortical dysplasia related to LAMB1 mutations. Neurology 84, 2195–2197. doi: 10.1212/WNL.0000000000001607

PubMed Abstract | CrossRef Full Text | Google Scholar

Vidal, F., Baudoin, C., Miquel, C., Galliano, M. F., Christiano, A. M., Uitto, J., et al. (1995). Cloning of the laminin alpha 3 chain gene (LAMA3) and identification of a homozygous deletion in a patient with Herlitz junctional epidermolysis bullosa. Genomics 30, 273–280. doi: 10.1006/geno.1995.9877

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Lin, Z. J., Liu, L., Xu, H. Q., Shi, Y. W., Yi, Y. H., et al. (2017). Epilepsy-associated genes. Seizure 44, 11–20. doi: 10.1016/j.seizure.2016.11.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. Y., Zhou, P., Wang, J., Tang, B., Su, T., Liu, X. R., et al. (2018). ARHGEF9 mutations in epileptic encephalopathy/intellectual disability: toward understanding the mechanism underlying phenotypic variation. Neurogenetics 19, 9–16. doi: 10.1007/s10048-017-0528-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wirrell, E. C., Grossardt, B. R., Wong-Kisiel, L. C., and Nickels, K. C. (2011). Incidence and classification of new-onset epilepsy and epilepsy syndromes in children in Olmsted County, Minnesota from 1980 to 2004: a population-based study. Epilepsy Res. 95, 110–118. doi: 10.1016/j.eplepsyres.2011.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiszniewski, W., Gawlinski, P., Gambin, T., Bekiesinska-Figatowska, M., Obersztyn, E., Antczak-Marach, D., et al. (2018). Comprehensive genomic analysis of patients with disorders of cerebral cortical development. Eur. J. Hum. Genet. 26, 1121–1131. doi: 10.1038/s41431-018-0137-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, H., Tan, D., Wang, S., Song, S., Yang, H., Gao, K., et al. (2015). Genotype/phenotype analysis in Chinese laminin-alpha2 deficient congenital muscular dystrophy patients. Clin. Genet. 87, 233–243. doi: 10.1111/cge.12366

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., and Zhang, Y. (2015). I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res. 43, W174–W181. doi: 10.1093/nar/gkv342

PubMed Abstract | CrossRef Full Text | Google Scholar

Zenker, M., Tralau, T., Lennert, T., Pitz, S., Mark, K., Madlon, H., et al. (2004). Congenital nephrosis, mesangial sclerosis, and distinct eye abnormalities with microcoria: an autosomal recessive syndrome. Am. J. Med. Genet. A 130A, 138–145. doi: 10.1002/ajmg.a.30310

PubMed Abstract | CrossRef Full Text | Google Scholar