key: cord-0070736-qbgkal2j
authors: Wen, Xinmei; Zhu, Wenjia; Xia, Nan L.; Li, Qianwen; Di, Li; Zhang, Shu; Chen, Hai; Lu, Yan; Wang, Min; Xu, Min; Wang, Suobin; Shen, Xin-Ming; Lu, Jie; Da, Yuwei
title: Missense Mutations of Codon 116 in the SOD1 Gene Cause Rapid Progressive Familial ALS and Predict Short Viability With PMA Phenotype
date: 2021-11-17
journal: Front Genet
DOI: 10.3389/fgene.2021.776831
sha: a80a700c74ee82fa5442ccd94a535984e113b30b
doc_id: 70736
cord_uid: qbgkal2j

Amyotrophic lateral sclerosis (ALS) is the most common form of motor neuron disease, characterized by a great variety of both clinical presentations and genetic causes. Previous studies had identified two different missense mutations in SOD1 (p.R116C and p.R116G) causing familial ALS. In this study, we report a novel heterozygous missense mutation in the SOD1 gene (p.R116S) in a family with inherited ALS manifested as fast-deteriorating pure lower motor neuron symptoms. The patient displayed similar clinical picture and prognostic value to previous reported cases with different R116 substitution mutations. Modeling of all R116 substitutions in the resolved SOD1 protein structure revealed a shared mechanism with destroyed hydrogen bonds between R116 and other two residues, which might lead to protein unfolding and oligomer formation, ultimately conferring neurotoxicity.

Amyotrophic lateral sclerosis (ALS) is a relentless neurodegenerative disorder and the most common adult-onset motor neuron disease (MND), which causes muscle weakness and atrophy. It is progressively fatal, and the majority of patients die within 3 years after onset of symptom (Turner et al., 2013) . It is considered highly heterogeneous in nature, which displays a great variety of both clinical presentations and underlying pathogenic mechanisms. About 5%-10% of ALS cases are familial (fALS). The genetic spectrums of ALS vary among different populations. In European and American populations, the most common mutations are C9orf72 intronic repeat expansions, whereas in Chinese and other Asian populations, SOD1 mutations are most frequent (Zou et al., 2017; Liu et al., 2018) . So far, more than 180 distinct SOD1 mutations have been identified as genetic causes for ALS (Mathis et al., 2019) (https://alsod.ac.uk/). Mutant SOD1 proteins cause neurodegeneration mainly through a gain of toxicity pathway. Other possible mechanisms also include increased protein instability, protein aggregation, and probabilities of fibrillization (Wright et al., 2019) .

The genotype-phenotype correlation has been described for a few SOD1 mutations for disease durations but not for clinical features. p.A5V and p.G42S (previously denoted as p.A4V and p.G41S, respectively) SOD1 substitutions have been consistently associated with a fast progressive phenotype, while patients with a p.H47R (previously denoted as p.H46R) mutation showed a more benign phenotype and a much longer life expectancy (Battistini et al., 2010; Takazawa et al., 2010) . However, these mutations are marked by symptoms representing a varying degree of lower and/ or upper motor neuron involvement of spinal, bulbar, and cortical regions. The genetic and phenotypic heterogeneities have been suggested to play critical roles in determining disease prognosis and in designing clinical trials (Manera et al., 2019) . Whether there is a correlation between mutation status and clinical traits remains to be answered.

Here, we report a novel missense mutation in the SOD1 gene (p.R116S), causing rapid deteriorating lower motor neuron symptoms in a fALS. The function of the mutant protein was explored to try to find a correlation between perturbed protein functions and clinical severity.

The proband was from a Chinese family of Han ethnicity with 28 known family members from three generations ( Figure 1A) . Overall, four affected family members have been diagnosed with MND. Two members from this family (proband and patient II-6) were seen at the ALS/MND Clinic, Xuanwu Hospital, Beijing, China, and further information on additional family members (II-7, III-7) was provided by the proband afterward. Disease onset was defined as first reported symptoms of weakness, dysarthria, or dysphagia. Patients underwent extensive laboratory studies to rule out other causes of neuropathy. The study was approved by the local medical ethical committee, with all participants providing written informed consent.

DNA was extracted from peripheral blood samples using a standard phenol-chloroform method. Whole-exome sequencing was used to detect gene mutations. DNA libraries were prepared with KAPA Library Preparation Kit and enriched using SeqCap hybridization probes following the manufacturer's instructions (Roche, Basel, Switzerland). Captured DNA libraries were subsequently sequenced on the Illumina NovaSeq platform. Low-quality reads and adapter sequences were first excluded. Single-nucleotide variants (SNVs) and insertions and deletions (INDELs) were filtered using GATK. All called variants were annotated based on several public databases, including the 1,000 Genomes Project, Human Gene Mutation Database Professional, gnomAD, and China National Gene Bank (CNGB). The pathogenicity of the identified variants was analyzed mainly based on criteria following the American College of Medical Genetics and Genomics (ACMG) guidelines.

Candidate mutation detected in exome sequencing was confirmed by Sanger sequencing.

The sequence-based assessment of the potential pathogenicity of each variant was assessed by the SIFT web server (https://sift.bii. a-star.edu.sg), Provean Protein web server (http://provean.jcvi. org/genome_submit_2.php, PolyPhen-2 web server (http:// genetics.bwh.harvard.edu/pph2/), CamSol web server (https:// www-cohsoftware.ch.cam.ac.uk), and PopMuSiC web server (https://soft.dezyme.com). The protein crystal structure of the wild-type human SOD1 was retrieved from the Protein Data Bank repository (PDB code 2C9V), in silico mutagenesis was performed for the three mutations discussed in this study using Pymol (The PyMOL Molecular Graphics System, Schrödinger, LLC, http://www.pymol.org), and the Find function was used to identify potential hydrogen bonds.

The proband was a 50-year-old female presented with a 2.5month history of progressive muscle weakness. She first experienced a heavy feeling in her right leg when climbing stairs. It worsened over the next 2 weeks that she had trouble standing up from a squatting position without support. The weakness continued to progress and spread to her left leg in the following 2 weeks. By then, she noticed muscle wasting in her right thigh and experienced significant difficulty climbing stairs as well as trouble getting up from a sitting even with support. She soon showed problems reaching upward using her right arm. Neurological examination revealed mild to moderate weakness of proximal muscle groups in the lower limbs of both sides and the right upper limb. Muscle atrophy of the right thigh was noted. Deep tendon reflexes were decreased to normal with down-going plantar reflex on both sides. Other pathological reflexes were not elicited. There were no abnormal findings on sensory system examination. There were no cognitive or mood changes. Pulmonary function test indicated restrictive deficit. Electromyography showed in both upper and lower limbs as well as thoracic paraspinal muscles a neurogenic pattern. Sensory and motor nerve conduction velocities were within normal ranges. 18 F-FDG PET did not detect any hypometabolism at the motor and premotor cortices or frontal-temporal regions. 3 Tesla magnetic resonance imaging (MRI) of the brain and spinal cord were unremarkable. Diffusion tensor imaging (DTI) reconstruction of the bilateral cortical-spinal tracts (CSTs) was performed, and no overt deficits were observed ( Figure 1C ). Together, a suspected diagnosis of progressive muscular atrophy (PMA) was made, and based on the familial history of MND, her blood sample was also sent out for WES and genetic analysis.

The proband received standard treatment with Edaravone and Riluzole and was evaluated every 3-6 months. Unfortunately, due to the COVID-19 pandemic, her follow-ups were only conducted via Telemedicine; therefore, detailed neurological check-ups were not possible. Her initial ALSFRS-R score was 39/48, which declined to 22/48 at her 3-month follow-up when she developed bulbar symptoms. Her ALSFRS-R score further decreased to 16/48 at her 5-month follow-up. One year after the disease onset, she was completely bedridden and required fulltime non-invasive positive pressure ventilation (NIPVV). By then, her ALSFRS was 8/48. Six months later, she was admitted to the hospital and intubated due to respiratory failure and died 1 month after. The proband's mother (patient II-6) was diagnosed of MND at age 67 and died 1 year later. Her symptoms started with weakness and muscle wasting of the lower limbs. Later, she experienced weakness in upper limbs and complained of shortness of breath. The symptoms were more pronounced on the right side. Upon neurological examination, substantial proximal > distal bilateral muscle weakness of her lower limbs and mild weakness of her right upper limb were found. Atrophies of multiple muscle groups were also observed. Deep tendon reflexes were diminished as well. No pathological reflexes were detected, and both cranial nerves and sensory modalities were intact. There was no cerebellar sign. The pulmonary function test indicated restrictive deficit. Needle EMG showed chronic and active denervation in bulbar and spinal regions.

Patient III-7 was diagnosed of MND at age 52 and died of respiratory failure 1 year later. Her initial symptom included weakness in the upper limb. Another patient (II-7) in this family was also diagnosed of MND with lower-limb onset when she was 62 years old and died at age 63. Both underwent detailed laboratory studies that did not reveal any other cause of diffuse denervation. Both the mother (II-3) and the grandmother (I-2) of patient III-7 were suspected mutation carriers who died of other causes over the age of 70 years.

The DNA sample of patient III-15 was subjected to whole-exome sequencing (WES) to get a comprehensive profile of genetic variants. A novel heterozygous c.346C>A mutation (NM_000454.5, p.R116S, previously denoted as p.R115S) in the SOD1 gene, was revealed ( Figure 1B) . Due to the death of all other affected members, segregation of the mutations could not be confirmed in the family presented here. However, this mutation was absent in the databases of healthy controls, including Genome Aggregation Database, the Exome Sequencing Project (ESP), and Exome Aggregation Consortium (ExAC) databases. The mutation altered a highly conserved residue and was predicted to be deleterious by both SIFT and PROVEAN and ranked as Probably Damaging and Disease Causing by PolyPhen-2 and Mutation Taster, respectively ( Figure 1E ).

SOD1 R116S mutation-induced alterations in protein functions were explored using a battery of algorithms. Mutation at this site resulted in a polar to non-polar amino acid change, implying a significant change in amino acid chemistry (Grantham score of 110) and affecting a highly conserved residue ( Figures 1D,E) .

The impact of R116S mutation on protein solubility was predicted using the CamSol Intrinsic web server, which suggested that the R116S SOD1 mutant was less soluble compared to the wild-type protein. Additionally, free energy change upon mutation was predicted using PopMuSiC, and an arginine-to-serine change at codon 116 was predicted to have a calculated ΔΔG of 1.62, suggesting a destabilizing effect.

The potential effects of the mutation on SOD1 protein stabilization were further investigated by protein modeling. SOD1 protein is a homodimer composed of two β-barrel monomers with copper and zinc in the active site. Codon 116 is located at the dimer contact site, being part of the interlocking Greek key loop connecting two four-stranded anti-parallel β-sheets facing each other (Figure 2A ) (Wright et al., 2019) . Interatomic interactions between the side chain of mutated residue and residues in the vicinity were analyzed, and hydrogen bonds were represented in the form of pseudo bonds using Pymol. Arginine at position 116 was shown to form hydrogen bonds with neighboring amino acids, particularly glutamic acid 50 and cysteine 112. R116S was found to disrupt hydrogen bonding with E50 and C112, suggesting alterations in native protein conformation and structure destabilization, thereby serving as a basis for neurotoxicity ( Figures 2B,D) .

We found a Chinese fALS family carrying a novel R116S change in SOD1. All four patients displayed an aggressive disease course with disease duration less than 2 years from initial symptom to death. Molecular modeling and bioinformatics analysis on mutant SOD1 protein suggested that loss of hydrogen bonds at the dimer interface, together with decreases in protein solubility and stability, which might explain the aggressive nature associated with the mutation.

The mean age at onset was 58 years in the family. The proband and her mother were characterized by a rapid disease progression with prominent lower motor neuron symptoms. Physical examinations showed no evidence of upper motor neuron signs, and there were no CSTs deficits on DTI in the proband, suggesting a progressive spinal muscular atrophy (PMA) phenotype of MND. It was suggested that around 10% of MND patients clinically presented as PMA (Van den Berg-Vos et al., 2009). Genetic causes might partially determine the motor neuron/clinical phenotypes. Some specific SOD1 missense mutations have been associated with a consistent clinical presentation. Here, we summarized a group of SOD1 mutations that have been found to cause a predominant lower motor neuron syndrome ( Table 1) . Most of the mutations belong to exon 4 and are located at β-strand structures of SOD1 protein ( Figure 2A ; Table 1 ). Regardless of the mutation status, these patients were more likely to have a limb onset with bulbar involvement during disease course and were also associated with a rapid disease progression. Consistent with clinical features, pathological findings showed severe spinal lower motor neuronal loss with minimal involvement of upper Frontiers in Genetics | www.frontiersin.org November 2021 | Volume 12 | Article 776831 motor neurons of motor cortex and cortical spinal tracts ( Table 1) . Of notice, different missense mutations of the same site at the SOD1 gene could lead to different disease phenotypes (i.e., A5V, A5T, previously denoted as A4V, A4T) ( Table 1) . It is reasonable to speculate that the relative position, property of the side chain of the mutated residue, and nature of the adjacent amino acid of the mutated residues were likely to affect the stability and conformation of SOD1 protein, thus conferring toxicity on mutant proteins. To our knowledge, the R116S mutation has never been reported so far (http://alsod.iop.kcl.ac.uk). However, R116G was found to be the most common SOD1 mutation in German fALS patients, which accounts for up to 44% of SOD1 associated ALS patients and was suggested to arise from a common founder (Niemann et al., 2004; Rabe et al., 2010) . A different point mutation R116C was identified in one Italian sporadic ALS patient (Tortelli et al., 2013) . Consistently, these two mutations are also associated with rapid disease progression and predominant lower motor neuron symptoms ( Table 2) . R116 mutations of SOD1 in three distinct ethnic groups caused the similar severe phenotype, indicating that R116 is critical for the protein function. Additionally, it has been suggested that protein aggregation is favored by mutations that bring the net charge of the protein closer to neutrality (Lindberg et al., 2005) . Thus, the substitution of a positively charged amino acid (arginine) with a neutral one (glycine, cysteine, or serine) would probably promote protein aggregation and lead to neurotoxicity. Within the resolved protein structure of SOD1, R116 forms hydrogen bonds with E50 and C112, respectively, both of which were reported to play important roles in SOD1 protein stability and function. Specifically, E50 locates in a loop structure, whose alterations could cause β-sheet partial unfolding, rendering it more prone to aggregation (Ding and Dokholyan, 2008) . On the other hand, C112 is a primary target for redox attacks, the oxidation of which results in conformational changes of SOD1 protein and formation of oligomers, ultimately leading to axon transport deficits and motor neuron death (Bosco et al., 2010; Nagano et al., 2015) . Notably, all R116 substitution mutations reported to date, including R116S, R116C, and R116G, are predicted to abolish the two hydrogen bonds and likely to affect the relative positions of E50 and C112 ( Figures 2C-E) , which in turn might cause β-sheet partial unfolding and might make C112 more accessible for oxidation, leading to mutant SOD1 misfolding and toxicity.

Taken together, we report a SOD1 mutation site causing the most consistent phenotype: ubiquitous limb onset characterized by prominent LMN presentation with rapid progression, adding to the list of ALS causing mutations bearing a genotype-phenotype correlation. Additionally, it displays a prognostic value of SOD1 R116 missense mutation. 

The datasets for this article are not publicly available due to concerns regarding participant/patient anonymity. Requests to access the datasets should be directed to the corresponding author.

The studies involving human participants were reviewed and approved by the Xuanwu Hospital Medical Ethical Committee. The patients/participants provided their written informed consent to participate in this study.

Conceptualization: XW and YD; patient evaluation: XW, WZ, LD, SZ, and YD; electrophysiology evaluation: YL and HC; protein structure modeling: XW and NX; imaging acquisition and analysis: QL and JL; writing: XW and YD; resources: MW, MX, SW; and supervision: X-MS, JL, and YD. All authors read, discussed, and approved the manuscript as submitted.

Clinical Characteristics of Familial Amyotrophic Lateral Sclerosis with CuZn Superoxide Dismutase Gene Mutations

A4T Mutation in the SOD1 Gene Causing Familial Amyotrophic Lateral Sclerosis

Phenotypic Heterogeneity in Motor Neuron Disease Patients with CuZn-Superoxide Dismutase Mutations in Scandinavia

Sixteen Novel Mutations in the Cu/Zn Superoxide Dismutase Gene in Amyotrophic Lateral Sclerosis: a Decade of Discoveries, Defects and Disputes

Variance of Age at Onset in a Japanese Family with Amyotrophic Lateral Sclerosis Associated with a Novel Cu/Zn Superoxide Dismutase Mutation

G41SSOD1mutation: A Common Ancestor for Six ALS Italian Families with an Aggressive Phenotype

Wild-type and Mutant SOD1 Share an Aberrant Conformation and a Common Pathogenic Pathway in ALS

Clustering of ALS Patients in central Italy Due to the Occurrence of the L84F SOD1 Gene Mutation

Progressive Muscular Atrophy Variant of Familial Amyotrophic Lateral Sclerosis (PMA/ALS)

Limited Corticospinal Tract Involvement in Amyotrophic Lateral Sclerosis Subjects with the A4V Mutation in the Copper/zinc Superoxide Dismutase Gene

A Novel L67P SOD1 Mutation in an Italian ALS Patient

Dynamical Roles of Metal Ions and the Disulfide Bond in Cu, Zn Superoxide Dismutase Folding and Aggregation

Two Novel Mutations in the Gene for Copper Zinc Superoxide Dismutase in UK Families with

Anticipation and Phenotypic Heterogeneity in Korean Familial Amyotrophic Lateral Sclerosis with Superoxide Dismutase 1 Gene Mutation

Autosomal Dominant Amyotrophic Lateral Sclerosis: a Novel Mutation in the Cu/Zn Superoxide Dismutase-1 Gene

Systematically Perturbed Folding Patterns of Amyotrophic Lateral Sclerosis (ALS)-associated SOD1 Mutants

The Epidemiology and Genetics of Amyotrophic Lateral Sclerosis in China

Regional Spreading of Symptoms at Diagnosis as a Prognostic Marker in Amyotrophic Lateral Sclerosis: a Population-Based Study

Genetics of Amyotrophic Lateral Sclerosis: A Review

Comprehensive Analysis of the Mutation Spectrum in 301 German ALS Families

A Cysteine Residue Affects the Conformational State and Neuronal Toxicity of Mutant SOD1 in Mice: Relevance to the Pathogenesis of ALS

Familial ALS in Germany: Origin of the R115G SOD1 Mutation by a Founder Effect

A SOD1 Gene Mutation in a Patient with Slowly Progressing Familial ALS

The Epidemiology of CuZn-SOD Mutations in Germany: a Study of 217 Families

SOD1 Missense Mutation in an Italian Family with ALS

The G93C Mutation in Superoxide Dismutase 1: Clinicopathologic Phenotype and Prognosis

Genotype-phenotype Correlation and Evidence for a Common Ancestor in Two Italian ALS Patients with the D124G SOD1 Mutation

SOD1 Mutation Is Associated with Accumulation of Neurofilaments in Amyotrophic Lateral Sclerosis

Rapid Disease Progression Correlates with Instability of Mutant SOD1 in Familial ALS

Clinicopathological Phenotype of ALS with a Novel G72C SOD1 Gene Mutation Mimicking a Myopathy

Confirmation of the Severe Phenotypic Effect of Serine at Codon 41 of the Superoxide Dismutase 1 Gene

Familial Amyotrophic Lateral Sclerosis with a Mutation in the Cu/Zn Superoxide Dismutase Gene

Familial Amyotrophic Lateral Sclerosis with a Novel G85S Mutation of Superoxide Dismutase 1 Gene: Clinical Features of Lower Motor Neuron Disease

Amyotrophic Lateral Sclerosis: a New Missense Mutation in the SOD1 Gene

Controversies and Priorities in Amyotrophic Lateral Sclerosis

A Long-Term Prospective Study of the Natural Course of Sporadic Adult-Onset Lower Motor Neuron Syndromes

The Biophysics of Superoxide Dismutase-1 and Amyotrophic Lateral Sclerosis

Genetic Epidemiology of Amyotrophic Lateral Sclerosis: a Systematic Review and Meta-Analysis

We thank our patients and their family members for participating in this study. We thank members of the Division of Neuromuscular Disease at Department of Neurology, Xuanwu Hospital, for the suggestions and critical evaluation of this work.

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.