1 Introduction
Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease characterized by progressive degeneration of upper and lower motor neurons, leading to muscular weakness and atrophy [
1]. The incidence of ALS varies among different populations, ranging from 0.8–1.2 to 2.1–3.8 per 100 000 person–years in China and Europe, respectively [
2]. Accumulating evidence indicates that genetic background, environmental exposure, and aging are risk factors implicated in the pathogenesis of ALS [
3,
4]. Genetically, approximately 10% of ALS cases are familial, and over 40 genes have been directly linked to ALS through next-generation sequencing [
5]. In China,
SOD1 has been identified as the most frequent ALS pathogenic gene, followed by
TARDBP and
FUS [
6].
ALS typically occurs in individuals aged between 42 and 65 years, invariably leading to death due to respiratory failure 3–4 years after the onset [
1,
7]. Juvenile ALS (JALS), a rare form of ALS, is defined as age at onset less than 25 years [
8]. Some patients of familial ALS (fALS) subtypes, such as ALS2, ALS4, ALS5, ALS6, and ALS16, are juvenile-onset, and most of these cases are inherited as an autosomal recessive pattern. JALS is more frequently recognized to have a genetic origin than adult-onset ALS, suggesting that studies on JALS may provide an opportunity to identify more novel causative genes [
8]. Till now, several genes, including
SPTLC1 [
9],
FUS [
10],
ALS2 [
11],
SETX [
12],
SPG11 [
13], and
SIGMAR1 [
14] have been identified in JALS. These disease causative genes underlying JALS also play a role in adult-onset ALS [
9]. The identification of these genes contributed to disclosing the molecular mechanisms underlying ALS. However, the etiology is not fully illustrated, and further genetic research may discover novel disease-associated genes and provide new clues to the disease pathogenesis.
Mitochondrial malfunctions play a pivotal role in the development of ALS [
15–
17]. The mitochondria produce energy for cellular biochemical reactions, which are essential for the survival of motor neurons and are implicated in the process of generating reactive oxygen species (ROS). The changes occurring in the mitochondrial respiratory chain enzymes have been implicated in ALS pathophysiology [
18]. Haploinsufficiency of
C9ORF72, an ALS pathogenic gene, causes impairment of respiratory complex I assembly in patients with ALS [
19].
In this study, a novel ALS-associated gene, TRMT2B, was identified in a JALS family, and additional TRMT2B variants were discovered in large cohorts of patients with ALS. Functional analysis revealed that TRMT2B may participate in the pathogenesis of ALS by impairing the functional activity of the mitochondria.
2 Materials and methods
2.1 Patients and clinical analysis
A family (XY003) that included two members who had ALS was recruited in Hunan Province, China (Fig.1). The two patients and their one unaffected sister and parents underwent full physical and clinical examinations, including electromyography, electrocardiography, brain magnetic resonance imaging (MRI), cognitive function evaluation, and blood biochemistry tests.
In this study, 910 cases with ALS from China were enrolled for mutation screening of the candidate genes [
20]. This ALS cohort included 753 patients with ALS, which were reported previously [
7]. Among these patients, pathogenic polynucleotide expansions in
C9ORF72 and
ATXN2 were excluded. The patients’ demographic characteristics are listed in Tab.1.
TRMT2B gene variants were also analyzed in an independent publicly available ALS cohort (ALSdb) comprising 3317 sporadic ALS cases subjected to whole-exome sequencing (WES). For controls, 104 068 exomes of non-neuro individuals in the gnomAD version 2 database were analyzed. The ALS-associated genes excluded for all patients are shown in Table S1.
This study was approved by the Ethics Committee of Xiangya Hospital, Central South University, China, under approval number 202103191. All participants provided written informed consent to participate in the study.
2.2 Karyotype analysis and whole-genome genotyping
Karyotype analysis was performed on peripheral blood lymphocytes by using the G-banding method as described before [
21]. Genomic DNA was extracted from the peripheral blood of all participants [
22]. Whole-genome genotyping was performed using an Illumina ASA 750K Chip. Genotype analysis was carried out using the Illumina GenomeStudio Genotyping Module (version 2011.1). The CNV gap value was set as 100 M.
2.3 WES
Samples of the two patients and their unaffected parents were subjected to WES as previously described [
7,
23]. Different filtering pipelines were used as the disease may be inherited by the family in an X-linked or autosomal recessive pattern. Rare damage variants (RDVs) that fulfilled the following criteria were included for further analysis: (1) rarity: heterozygous variants with a minor allele frequency (MAF) of less than 0.1% in the 1000 Genome Project and gnomAD; (2) nonsynonymous substitutions, indels, and putative splice site variants; (3) variant pathogenicity predicted by at least two of 11
in-silico tools; (4) for the X-linked pattern, all consensus variations in the X chromosome were extracted from both patients in the family, and for the autosomal recessive inheritance pattern, all consensus homozygous or compound heterozygous variations were extracted from both patients in the family. The samples of the healthy sister and parents were included for co-segregation analysis.
Further screening of candidate variants was performed on the samples of the 910 patients with ALS in the cohort who underwent WES [
7].
2.4 Mutant mRNA expression and qPCR
Lymphoblastoid cells were isolated from blood samples and immortalized using the Epstein–Barr virus [
22]. Total RNA was extracted as described before [
22]. PCR was performed with primers TRMT2B-F/R (Table S2) by using Premix (Takara). The amplified product was inserted into the T-vector (Takara) and sequenced. qPCR was performed with actin and ND1 primers (Table S2) by using TB Green Premix (Takara) on QuantStudio (Invitrogen).
2.5 Plasmid construction and cell transfection
TRMT2B cDNA (NM_001167972.2) was amplified by PCR using primers TRM-F/R (Table S2). Premix (Takara) was used to amplify the target sequence and inserted into the PCDNA3.1(+) plasmid by using T4 ligase (Thermo Scientific). The mutated TRMT2B cDNA was amplified using specific primers (Table S2) and cloned as described above. HEK293 cells were transfected with jetPRIME in accordance with the manufacturer’s instructions.
2.6 Gene interfering
TRMT2B siRNA with target sequence of 5′-CTGGTCAAGCAGAGAAGATTT-3′ and 5′-GCACAGTATGTAAGGGAGATT-3′ was purchased from Sangon. HEK293 cells were transfected with jetPRIME in accordance with the manufacturer’s instructions. RNAi lentivirus particles, with identical target sequences, were purchased from Genechem. Samples were tested 14 days after the lentivirus infection.
2.7 Immunoblotting and immunofluorescence
For immunoblotting, the lymphoblastoid cells from the patients and normal participants were lysed in RIPA buffer. Protein electrophoresis was conducted as described before [
22]. The expression of TRMT2B and p62 was analyzed using antibodies. Anti-β-actin antibodies were used for loading control.
For immunofluorescence, HEK293 cells cultured in DMEM with 10% FBS (Invitrogen) were plated on collagen-coated coverslips and fixed with 4% ice-cold paraformaldehyde. The procedures of immunofluorescence were identical as described before [
24]. Images were analyzed by confocal microscopy (Zeiss). The antibodies used in this study are listed in Table S3. Co-localization analysis was performed using the ImageJ plugin colocalization-finder.
2.8 Activities of mitochondrial complexes I and III
In each sample, 5 × 106 lymphoblastoid cells were used for mitochondrial extraction. The activity of mitochondrial complex I was assessed by the NADH oxidation rate measured on the basis of absorbance at 340 nm. The activity of mitochondrial complex III was assessed by the rate of cytochrome C increase on the basis of absorbance at 550 nm. The assays on mitochondrial complexes I and III were performed using the mitochondrial complexes I and III activity kits, respectively (Solarbio).
2.9 Mitochondrial aerobic respiration
In each sample, 10
5 cells were seeded on a seahorse assay plate and cultured overnight. The medium was then replaced with the Seahorse XF base medium, and the cells were incubated in a 37 °C incubator without CO
2 for 1 h. After the baseline oxygen consumption rate (OCR) was measured, oligomycin (1 µmol/L, complex V inhibitor), carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone (FCCP, 2 µmol/L, uncoupler), and rotenone/antimycin A (0.5 µmol/L each, complexes I and III inhibitors) were added sequentially. The OCR measurements, including basal, maximum respiration, spare-respiratory capacity, ATP synthesis, proton leak, and non-mitochondrial respiration, were obtained and calculated as previously described [
25]. The OCR values were normalized to the sample protein concentration.
2.10 Flow cytometry detection of ROS
Cells were incubated with dihydroethidium (Applygen) diluted in PBS (1:1000) for 30 min at 37 °C, washed in PBS three times to remove the residual probe, and analyzed by flow cytometry using FACSCanto II (BD). The data were evaluated using the FlowJo software.
2.11 Transmission electron microscopy (TEM)
Cells were incubated in iced 2.5% glutaraldehyde solution for 1 h and then washed with iced PB buffer five times. Next, 2% osmic acid was added to each sample, and the samples were washed with iced ddH2O five times. They were then replaced with ethanol and embedded with SPI-Pon 812. Leica EM UC7 ultramicrotome was used to obtain 70 nm slices, and classical staining with 2% uranyl acetate and lead citrate was performed. All the TEM images were captured using Hitachi 7700.
2.12 Statistical analysis
Statistical analysis was performed using SPSS 25.0. Association analysis of RDVs across TRMT2B was performed in the ALS cohort by using Fisher’s exact test for each allele. An unpaired t-test was conducted using Prism. A P value of 0.05 indicated statistical significance.
2.13 Data availability statement
The datasets analyzed in this study are not publicly available. Further information about these datasets is available from the senior authors (Jungling Wang and Zhengmao Hu) upon reasonable request.
3 Results
3.1 Clinical characteristics of patients with JALS
The proband of the XY003 family (Fig.1) was a patient in his 30s who initially developed weakness of the lower limbs at childhood age. Before the age of 10 years, weakness of the upper limbs occurred. The patient also showed delayed cognitive and motor developmental milestones, and his school performance was worse than average. At the age of 20s, he suffered from a decrease in visual acuity of the left eye. He complained of progressive exacerbation of limb weakness but not of dysphagia or dysarthria. A detailed medical record was not available. On admission, the patient had a small stature, which was below the genetic expectation, and deformity of hands and feet (Fig. S1A–S1C). He also showed low activity endurance and ability to exercise. A neurological examination revealed obvious muscle weakness and atrophy of the four extremities, facial muscle bundle tremor, and normal eye movement but incomplete eye closure. Furthermore, the proband showed a slight increase in muscle tension in the lower limbs and stiffness of ankle joints, which induced ankle clonus. He also showed cognitive decline (executive dysfunction and memory impairment). The proband’s younger brother (XY003.P2) showed the same lower limb weakness at childhood age and similar symptoms at the age of 20s (Fig.1 and S1D–S1F). The detailed clinical information is described in Tab.2. Their sister and parents had no neurological symptoms or signs.
Nerve conduction studies revealed an obvious reduction in compound action potential in the lower and upper limb nerves of the proband, but nerve conduction velocity (NCV) was normal. Meanwhile, sensory NCV studies did not detect any abnormalities. The electromyogram of the proband showed chronic neurogenic reinnervation, enlarged long-duration motor unit potentials in the cervical and lumbar spinal cord regions, and ongoing denervation (moderate spontaneous potentials in lower limb muscles). Chronic neurogenic reinnervation and ongoing denervation were not found in the sternocleidomastoid muscle. Left visual evoked potential (VEP) indicated an obvious reduction in the amplitude of the P100 component. The proband’s younger brother also showed chronic neurogenic reinnervation in the cervical and lumbar regions and ongoing denervation in the thoracic region (moderate spontaneous potentials). He had the same VEP reduction as the proband. The proband and his younger brother had increased rest levels of creatine kinase (349 and 1099 U/L, respectively, vs. normal level of 50–310 U/L). In both patients, blood lactate was normal before exercise and slightly increased (less than twice) after exercise. The blood levels of ammonia, glucose, and lipids and thyroid function were normal. No obvious abnormalities were revealed by brain MRI (Fig. S2), ECG, and echocardiography.
3.2 Identification of TRMT2B variants associated with JALS
All coding and splice variants that were identified in the X chromosomes of the two affected patients were analyzed. After functional filtering, a hemizygous variation c.1356G>T (p.K452N, NM_001167972.2) was identified in TRMT2B. Sanger sequencing revealed segregation of the c.1356G>T mutation with the phenotypes in the XY003 family (Fig.1). The mother of the proband was heterozygous for this variant (Fig.1). The c.1356G>T mutation in TRMT2B has not been reported in 1000 Genomes Project or gnomAD databases, and it is highly conserved among species (Fig.1).
As TRMT2B is located on the X chromosome and the unaffected mother of the proband was a heterozygous carrier, TA cloning was conducted to determine the presence of the mutated TRMT2B mRNA in the lymphoblastoid cells of the mother. The results revealed only 8.5% (6/71) of mutant mRNA in the lymphoblastoid cells of the proband’s mother, less than the expected rate of 50%.
The karyotypes of the patients and their mother were normal (Fig.1). CNV analysis on the two patients showed no consensus abnormities (Fig. S3). An analysis based on an autosomal recessive hereditary pattern was performed; however, only two genes fulfilling the analysis criteria were found, and none of them was associated with potentially neurodegenerative processes (Table S4).
3.3 Screening of TRMT2B variants in patients with ALS
The occurrence of TRMT2B variants in a cohort of 910 patients with ALS was evaluated to explore the variation role of TRMT2B in ALS. Among 10 patients, two missense variants (c.250C>G, L84V and c. 1344T>G, F448L) and one splice variant (c.539-3T>C) of TRMT2B were identified, which were rare or absent in gnomAD. Both missense variants were predicted to be damaging (Tab.3 and Fig. S4). Notably, the L84V variant was identified in seven unrelated patients with sporadic ALS who had typical clinical features of adult-onset ALS without vision impairment. The mean age at onset was 57.22 ± 7.96 years (range of 42–66 years). A notable detail that three of the six patients who underwent cognitive tests had cognitive impairment (50%).
Among the 104 068 non-neuro individuals in the gnomAD version 2 database (47 831 females and 56 237 males), a total of 494 rare putative pathogenic TRMT2B variants from the WES data were identified in TRMT2B that fulfilled the same RDV analysis criteria. The results of gene burden testing revealed TRMT2B variations as a risk factor for the disease (P = 0.003 by Fisher’s exact test). For variation L84V, the results of gene burden test showed significance in patients (P < 0.0001) compared with GnomAD. For variation F448L, the results of gene burden did not show significance in patients (P = 0.0961 by Fisher’s exact test) compared with GnomAD or ExAC (P > 0.99 by Fisher’s exact test). However, the results of gene burden showed significance in patients compared within GnomAD male section (P = 0.0460 by Fisher’s exact test).
In an independent ALS cohort (ALSdb), three additional rare deleterious TRMT2B variants, namely, R334W (c.1000C>T), R315L (c.944G>T), and R162G (c.484C>G), were identified (Tab.3 and Fig. S4). However, no clinical data were available online.
3.4 Patients with TRMT2B variations showed impaired mitochondrial function
Lymphoblastoid cell lines of the two patients and controls, including their heterozygous mother, the healthy sister of the XY003 family, and one unrelated man, were constructed. The unrelated man aged 24 years at the collection of the sample. He was mentally and physically healthy, and he received his master’s degree 1 year after sample recruitment. He had no family history of neurodegeneration diseases, such as Alzheimer’s disease (AD), Parkinson’s disease, or ALS.
Compared with the controls, decreasing numbers of mitochondria (Fig.2 and 2B, P < 0.0001 by unpaired t-test) and swollen mitochondria (Fig.2 and 2D, P < 0.0001 by unpaired t-test) were found on TEM. Considering the TEM results showed morphological alteration in the mitochondria, putative functional alteration in the mitochondria was further detected among the patients. The patients had lower expression levels of ND1 at protein (Fig.3 and 3D, P = 0.0465 by unpaired t-test) and mRNA levels (Fig.3, P = 0.0331 by unpaired t-test) than the controls. However, the expression of TRMT2B seemed unchanged between the patients and the controls (Fig.3 and 3B). The activity of mitochondrial complexes I and III in lymphoblastoid cells was also examined from the two patients and the controls. The activities of complex I in the patients were lower than those in the other tested individuals (P = 0.0251 by unpaired t-test), including the carrier mother (P= 0.0356 by unpaired t-test). Although no significant difference was found among the three controls (carrier mother, healthy sister, and the unrelated man), the average value of the carrier mother was lower than that of the other controls. Meanwhile, no significant differences were observed in the activities of complex III (Fig.3 and S5). The level of ROS in the lymphoblastoid cells of the patients also increased (Fig.3 and S6; P < 0.0001 by unpaired t-test).
3.5 TRMT2B contribution to maintaining mitochondrial function
Two TRMT2B-interfering HEK293 cell lines were constructed using lentivirus to further explore the function of TRMT2B (Fig. S7A). A reduced number of mitochondria (Fig.4 and 4B) and swollen mitochondria (Fig.4 and 4D) were observed in the TRMT2B-interfering cell lines.
By interfering with the expression of TRMT2B in HEK293, the expression of ND1 decreased (Fig.5 and 5B) in patients with TRMT2B variations. Further, no alteration was detected in the expression of MTCO1, MTCO2, or TFAM (Fig.5). By overexpression of TRMT2B and TRMT2B mutated proteins in HEK293, a the lower expression of ND1 was found in the mutated groups (Fig.5 and 5E). Compared with those in NC, the basal respiration, ATP production, maximal respiration, and spare respiratory capacity in the two TRMT2B-interfering HEK293 cell lines decreased according to the Seahorse XF Pro Analyzer, which indicated decreasing mitochondrial aerobic respiration induced by the down-expression of TRMT2B (Fig. S8A–S8G). Increased ROS levels were also detected in patients (Fig.5 and S7B).
Immunocytochemistry analysis of HEK293 cells indicated that TRMT2B was co-localized with the mitochondrial marker TOMM20, and the mutated TRMT2B proteins were also located within the mitochondria (Fig. S9).
Thus, TRMT2B may play an important role in maintaining normal mitochondrial function.
4 Discussion
In this study, a novel ALS-associated gene called
TRMT2B was identified by analyzing patients with familial and sporadic ALS. The proband of the family exhibited progressive muscular weakness and atrophy, cognitive impairment, and deformities in all four extremities, as well as impaired vision. The neurological examination revealed signs of upper and lower motor neuron lesions. The neurological electrophysiology analysis showed extensive chronic and ongoing neurogenic damage, indicating damage to the anterior horn cells in the spinal cord and pyramidal tract. After other neurological diseases were ruled out, the patient was diagnosed with fALS and ALS-Plus syndrome in accordance with the revised El Escorial criteria from 2000 [
20,
26]. ALS-Plus is used to classify patients with ALS who exhibit non-motor symptoms, such as oculomotor disorders, sensory disturbances, and cerebellar and extrapyramidal disorders [
26,
27].
In this study, other diseases were excluded from the differential diagnosis. Multisystem proteinopathy (MSP) is a pleiotropic degenerative disorder characterized by a combination of at least two of the following: body myopathy, Paget disease of the bone, and ALS/frontotemporal dementia [
28]. Several ALS causative genes, such as
VCP,
SQSTM1,
HNRNPA2B1,
HNRNPA1, and
MATR3, have been identified to be correlated with MSP phenotypes [
28]. In the present study, MSP was also a possibility given the bone deformity and motor neuron dysfunction in the proband. However, on the basis of the phenotypes and ECG results of the patients, they were diagnosed with ALS-Plus. In addition, spinal muscular atrophy (SMA)-Plus and mitochondrial encephalomyopathy were not the appropriate diagnoses. The UMN signs did not support a diagnosis of SMA-Plus [
29], and mitochondrial encephalomyopathy was ruled out because of the un-elevated level of lactate, the electrophysiological chronic and ongoing neurogenic changes instead of myogenic changes, and the normal brain MRI. Therefore, ALS-Plus was the more appropriate diagnosis, which may extend the phenotypic spectrum of the ALS-Plus syndrome.
TRMT2B is a gene located on the X chromosome that encodes tRNA methyltransferase 2 homolog B. Although the function of TRMT2B is not fully understood, it is believed to be responsible for the m
5U methylation in mitochondrial tRNA and rRNA, which is crucial for tRNA stability and maturation [
30]. Inactivation of
TRMT2B has been shown to decrease the activity of respiratory chain complexes I, III, and IV [
31]. While one study suggested that
TRMT2B may not be essential in human cells [
32], dysregulation of tRNA post-transcriptional methylation has been linked to neural developments [
33] and neurodegenerative diseases [
30]. Other members of the TRMT gene family, such as
TRMT61B [
34],
TRMT1, and
TRMT10A [
35,
36], were supported to contribute to the pathogenesis of neurological disorders, such as AD and intellectual disability. Meanwhile, a previous study indicated that mice with a mutated
Trmt2b exhibited abnormal physical strength [
37]. These findings and the screening results of TRMT2B in the cohort may provide independent, albeit limited, evidence for the pathogenicity of
TRMT2B in ALS.
Some ALS-associated genes, including
TARDBP,
C9ORF72, and
SOD1, have been reported to play a role in the etiology of ALS by causing mitochondrial dysfunction [
38–
40]. An increasing body of research data suggested that imbalanced ROS plays a remarkable role in the onset and progression of various neurodegenerative disorders, such as ALS, AD, and PD [
41]. Herein, in the patients’ lymphocytes and
TRMT2B-interfering cell lines, mitochondrial morphology alteration, reduced ND1 expression, lower mitochondrial complex I activity, decreased mitochondrial aerobic respiration, and increased ROS levels were observed, indicating mitochondrial dysfunction induced by TRMT2B loss of function.
TA cloning and sequencing showed that only 8.5% of the TRMT2B mRNA in the unaffected carrier mother had the c.1356G > T substitution, suggesting that the majority of the protein could be unaffected. The mitochondrial complex I activity in the carrier mother was normal and higher than that in patients, although the average value was lower than that in other normal controls, similar to the ROS level (Fig. S6). These results indicated that the accumulative effect of the lower activity did not reach the threshold for ALS onset. Consequently, the TRMT2B functional alterations identified in patients with ALS, especially women, may contribute to pathogenesis through dose- and time-dependent accumulation.
Due to technical limitations, mitochondrial tRNA methylation could not be directly examined using isotope labeling of tRNA. Alternative methods should be utilized in future studies. In addition, the lack of lymphoblastoid cell lines prevented the examination of mitochondrial complex activity and ROS production in patients with sporadic ALS with mutated TRMT2B. Future studies could concentrate on investigating the functions of mutated TRMT2B and validating accumulated effects.
In conclusion, a novel JALS-associated gene called TRMT2B was identified, which may contribute to the pathogenesis of ALS through the impairment of mitochondrial function.