1. Introduction
The advent of next generation sequencing techniques has greatly eased the identification of novel genes causative of Mendelian disorders while disclosing a previously unsuspected degree of pleiotropy. Indeed, variants in the same gene are often found to cause a broad spectrum of phenotypes, with variable severity and clinical presentation, possibly encompassing distinct disorders. The reasons underlying such pleiotropy are still mostly unknown, and genotype-phenotype correlates remain scarce.
One such example is the peptidase mitochondrial processing alpha (
PMPCA) gene, encoding the alpha subunit of mitochondrial processing peptidase (α-MPP), a heterodimeric enzyme responsible for the cleavage of nuclear-encoded mitochondrial precursor proteins after import in mitochondria. α-MPP is characterized by the glycine-rich loop (GRL, residues 345–358), a key structural element of MPP that recognizes and binds the precursor protein substrate toward the active site [
1,
2,
3].
Biallelic pathogenic missense variants in the
PMPCA gene have been associated with two distinct conditions. The first is an early onset, non-progressive cerebellar ataxia (NPCA) with or without developmental delay and intellectual disability, termed “recessive spinocerebellar ataxia type 2” (SCAR2), which should be kept distinct from SCA2 (spinocerebellar ataxia type 2), an autosomal dominant, adult-onset spinocerebellar ataxia due to a pathological trinucleotide expansion in the ataxin 2 gene. Nineteen patients from five families have been reported, all presenting a typical pattern of non-progressive cerebellar atrophy at brain imaging. Muscle biopsy was normal, although mitochondrial complexes activity was mildly decreased in some cases [
4,
5]. Patients from three families were homozygous for the p.Ala377Thr variant, the proband in the fourth family was compound heterozygous for the p.Ser96Leu and p.Gly515Arg variants, while the two affected siblings in the fifth family were homozygous for the p.Val256Met variant.
The second condition is a much more severe progressive mitochondrial encephalopathy, presenting with profound developmental delay, ptosis, ophthalmoplegia, weakness, respiratory insufficiency, blindness, cardiomyopathy, liver failure, and lactic acidemia, reported in two siblings who were compound heterozygous for the recurrent p.Ala377Thr and the novel p.Gly356Ser variant. Brain MRI showed diffuse cerebellar and cerebral atrophy, with a lactate peak at MRI spectroscopy, while muscle biopsy demonstrated enlarged, structurally abnormal mitochondria [
3].
Recently, we reported a 7-year-old Italian boy (here termed “Ita-2”) homozygous for a novel
PMPCA missense variant (p.Arg185Trp), presenting a phenotype of intermediate severity characterized by a progressive encephalopathy without extra-neurological involvement. He had psychomotor delay with regression, intellectual disability, a spastic-ataxic gait, and upper limb dystonia. Mitochondrial involvement was demonstrated by a slight elevation of blood lactate, as well as multiple defects of oxidative metabolism and signs of mitochondrial proliferation at muscle biopsy. In addition to cerebellar atrophy, brain MRI showed bilateral symmetric hyperintensity in the striatum, resembling the pattern of bilateral striatal necrosis typical of Leigh syndrome [
6]. A similar presentation was very recently reported in a 15-year-old Japanese patient [
7].
Here, we further define the PMPCA-associated phenotypic spectrum by reporting a 12-year-old Italian girl with a similar “intermediate” phenotype of progressive encephalopathy lacking extra-neurological signs. Although she did not present biochemical signs of mitochondrial involvement, brain imaging showed the same peculiar combination of cerebellar atrophy and bilateral “Leigh-like” striatal hyperintensity, suggesting this neuroradiological feature may represent a hallmark of the PMPCA-associated intermediate phenotype.
2. Patients and Methods
2.1. Genetic Analysis
The proband (NG2880) was ascertained through a research project aimed at detecting the underlying genetic defect in a large cohort of patients with pediatric ataxias. The local ethics committee approved the project, and parents signed a written informed consent. Clinical, laboratory, and imaging data were obtained through routine diagnostic examinations. Parents consented to a skin biopsy for research purposes; cultured fibroblasts of the Ita-2 patient were already available.
Clinical exome sequencing was performed using SureSelect QXT Focused Exome (Agilent, Santa Clara, CA, United States) and run on NetSeq500 sequencer (Illumina, San Diego, CA, United States), as per manufacturer’s protocol. Bioinformatic analysis was carried out as reported [
8]. Variants were filtered based on effect (retained only variants affecting protein sequence or splicing), frequency (minor allele frequency <1% in general population), and quality (reads >10 and positive quality passive filter). Retained variants were classified according to American college of Medical Genetics (ACMG) criteria, using Varsome [
9]. Candidate variants were confirmed by Sanger sequencing using the Big Dye Terminator chemistry (Applied Biosystems, Foster City, CA, USA).
2.2. Cell Cultures
Skin-derived primary fibroblasts from both Italian patients (NG2880 and Ita-2) and three healthy controls were cultured in Roswell Park Memorial Institute Medium (RPMI 1640 w/o L-glutamine) (Carlo Erba Reagents, Cornaredo, Italy) supplemented with inactivated 20% fetal bovine serum (FBS), 1% L-glutamine 200 mM, and 1% penicillin/streptomycin 100X (Euroclone, Pero, Italy), at 5% CO2 and 37 °C. The medium was replaced three times per week until fibroblasts reached confluence. Fibroblasts were detached with trypsin (Trypsin-EDTA with Phenol Red, Euroclone) and either replated, pelleted or cryopreserved.
2.3. Real-Time PCR Analysis
Total RNA was extracted from fibroblasts with the combination of Trizol and column kit activity (Direct-zolTM RNA MiniPrep Plus w/Zymo-SpinTM IIICG Columns, Euroclone). A total amount of 1000 ng was reverse transcribed using High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific, Waltham, MA, USA).
After cDNA synthesis, 24 ng for each sample (1.6 ng/μL), together with 14.4 μL PowerUp™ SYBR™ Green Master Mix (Applied Biosystems) and 0.6 μL primer mix (40 μM each) were loaded in duplicate into hard-shell 96-well PCR plates and run on CFX Connect detection system (Biorad, Hercules, CA, United States). Primer sequences are available upon request. Data were analyzed with the ΔΔCT method and glyceraldehyde-3-phopaspate dehydrogenase (GAPDH) was used as internal control. Statistical analysis was performed using two-tailed Student’s t-test.
2.4. Cell Extracts and Western Blot
Protein lysates were obtained by resuspending fibroblast pellets in 100 µL of CellLytic lysis buffer containing protease and phosphatase inhibitors (Sigma-Aldrich, St. Louis, MI, United States). Samples were kept on ice and placed on a shaker for 30 min, then centrifuged 16,000× g for 20 min at 4 °C. Proteins quantification was performed by the bicinchoninic acid method (Sigma). The absorbance was measured at 562 nm wavelength on CLARIOstarPlus microplate reader (BMG LABTECH, Ortenberg, Germany).
Protein samples (20 μg) were resolved by electrophoretic run on a Mini-Protean® TGX™ gel and transferred to nitrocellulose membranes through the Trans-Blot Turbo Transfer System (Biorad). Membranes were blocked with a commercial buffer (Li-Cor Biosciences, Lincoln, NE, USA) and immunoblotted with primary (anti PMPCA, 1:500—Santa Cruz Biotechnology, Dallas, TX, United States; anti β-actin, 1:10,000—St. John’s Laboratory, London, UK) and secondary antibodies (IRDye 800 donkey anti goat and IRDye 700 donkey anti-mouse, Li-Cor). Blots were imaged with the Odyssey System (Li-Cor) and the fluorescent signal obtained from each protein was normalized with the corresponding β-actin signal.
Data were expressed as percentage compared to the average of control samples, in each membrane. Statistical analysis was performed using two-tailed Student’s t-test.
2.5. Molecular Modeling and Dynamics Simulation
Atomic coordinates of the human MPP heterodimer (Mitochondrial Processing Peptidase alpha and beta subunits, PMPCA and PMPCB) were obtained through homology modeling with SwissModel [
10] using the Yeast Mitochondrial Processing Peptidase structure (PDB id: 1hr6, chain A, sequence identity = ~36%) as a template for PMPCA and the Bovine Cytochrome bc1 complex subunit 1 (PDB id:
5okd, chain A, sequence identity = ~57%) as a template for PMPCB. Then, the PMPCA (residues 64 to 510) and PMPCB (residues 45 to 488) models were superimposed to the Yeast Mitochondrial Processing Peptidase structure using the MatchMaker extension of UCSF Chimera to recreate the heterodimer. The wild-type structure of the heterodimer was mutated in silico to introduce the following missense variants: p.Ala377Thr, p.Gly356Ser, p.Arg185Trp, p.Trp278Leu, p.Gly356Ser, and p.Val256Met.
These were simulated by the Gaussian accelerated Molecular Dynamics (GaMD) technique [
11]. A time step of 2 femtoseconds (fs) was set, and then 250 nanoseconds (ns) of GaMD simulations were performed. Each system was simulated three times using AMBER 18 [
12] running on 1 NVIDIA Quadro P6000 and 3 NVIDIA RTX 2080Ti [
13]. Root-mean-square deviation (RMSD), to measure the average distance between the positions of the Cα atomic coordinates and those of the reference X-ray structure over time, and dynamic cross-correlation maps (DCCMs), to study the long-range interactions between all pairs of atoms, were computed [
14].
4. Discussion
An increasing number of genes are progressively associated with a broad spectrum of phenotypes, even encompassing different disorders. In these cases, an accurate mapping of gene-associated phenotypes is mandatory, as it would help the diagnostic lab to select the appropriate gene panels (or virtual gene panels in case of exome sequencing), and therefore increase the diagnostic yield.
Here we contribute to delineate the phenotypic spectrum associated with recessive mutations in the PMPCA gene. So far, two well-distinct conditions have been delineated: (1) a severe progressive mitochondrial encephalopathy with multiorgan involvement, characterized by profound psychomotor delay, ptosis, ophthalmoplegia, hypotonia and severe generalized weakness, optic atrophy, cardiomyopathy, liver failure, global brain atrophy, and evidence of mitochondrial involvement (lactic acidemia, lactate peak at MRI spectroscopy, abnormal mitochondria at muscle biopsy); (2) a much milder neurological disorder characterized by non-progressive cerebellar atrophy and NPCA with or without intellectual disability, lacking any obvious mitochondrial sign.
We now confirm the existence of a third PMPCA-associated phenotype of intermediate severity, presenting as a progressive encephalopathy with psychomotor regression, intellectual disability, spastic ataxia, and a variable hyperkinetic movement disorder, but without extra-neurological signs of mitochondriopathy. In addition to cerebellar atrophy, patients present a peculiar neuroradiological hallmark, represented by bilateral symmetric T2-weighted hyperintensities of the basal ganglia, resembling those seen in Leigh syndrome. However, it is important to note that our patient lacked obvious signs of mitochondrial involvement (blood lactate, muscle biopsy, and MRI spectroscopy were normal), which instead were present in the first patient. Just before submission of this manuscript, a 15-year-old Japanese girl was reported with a similar phenotype of severe and progressive developmental delay, cerebellar ataxia, and extrapyramidal symptoms. Brain imaging showed cerebellar atrophy and abnormal signal in the basal ganglia, which was interpreted as iron accumulation; however, it is likely that this girl also fell in the “
PMPCA intermediate phenotype” category highlighted here, characterized by peculiar Leigh-like basal ganglia hyperintensities in association to cerebellar atrophy [
7].
All pathogenic variants in the
PMPCA gene are missense variants variably distributed throughout the protein (
Figure 3) that are highly conserved among vertebrates and consistently predicted as damaging by prediction software. Previous studies have attempted to define the functional impact of specific amino acid substitutions on the stability and function of α-MPP protein and correlate these effects with the observed phenotype. In particular, mutant α-MPP was found to impair processing of the mitochondrial protein frataxin, causing an accumulation of the intermediate forms and reduced levels of the processed protein [
15]. Interestingly, abnormal frataxin processing was demonstrated in lymphoblasts and/or fibroblasts from patients presenting either with the severe (compound heterozygous for p.Ala377Thr and p.Gly356Ser variants) or the mild phenotype (homozygous for p.Ala377Ser or for p.Val256Met variants) [
3,
4,
5]. Moreover, α-MPP levels were measured in lymphoblasts from patients presenting the mild phenotype homozygous for two distinct mutations: while cells with the p.Ala377Thr variant showed significantly reduced levels of α-MPP, these were comparable to controls in cells carrying the p.Val256Met variant [
4,
5]. Taken together, these findings suggest a loss of function mechanisms underlying
PMPCA missense variants (as also pointed out by the frameshift mutation recently reported in the Japanese patient) but fail to disclose any specific genotype-phenotype correlates at the functional level. This is further confirmed by our findings on fibroblasts from the two Italian patients with intermediate phenotypes, as we unexpectedly observed opposite trends regarding α-MPP expression levels, which were halved than controls in cells from the patient Ita-2 (homozygous for the p.Arg185Trp variant) but were significantly increased compared to controls in patient NG2880 (compound heterozygous for p.Trp278Leu and p.Arg362Gly variants). Since RNA expression levels were comparable in the two patients and similar to controls, the altered levels of α-MPP are likely related to the variable impact of missense variants on the protein stability, influencing its physiological turnover within the cell.
To further attempt a dissection of the functional effect of distinct missense variants, we performed molecular modeling of the human MPP heterodimer, and a systematic dynamic simulation of several missense mutants compared to wild type (
Figure 4). Such an approach was already proposed by Jobling and collaborators. They suggested that the substitution of a Serine with the bulkier Threonine at position 377, lying in close proximity to the GRL, would interfere with the loop flexibility and, therefore, with the ability to bind substrates stably [
4]. Similarly, Joshi and coworkers performed computer modeling for the same p.Ala377Thr and for p.Gly356Ser, which were compound heterozygous in their patients with severe phenotype. Since p.Gly356Ser falls within the GRL, the authors proposed that both mutations would restrict the dynamic movements of the catalytic loop, therefore affecting substrate binding or catalytic activity of the protein [
3].
Our widespread dynamic simulations showed that all analyzed mutants could variably impact on the GRL flexibility and some also influenced the overall stability of the entire MPP heterodimer. In particular, the variants p.Arg185Trp, p.Arg362Gly, and p.Ala377Thr exhibited lower RMSD values than the wt for the simulations’ entire duration, reflecting an increased rigidity and significant impairment of the GRL mobility. Conversely, the impact of p.Val256Met, p.Trp278Leu, and p.Gly356Ser on the GRL seemed less evident, although flexibility is always lower than wt. However, these three variants (as well as p.Arg362Gly) showed a more relevant impact on the overall stability of the heterodimer, as they induced a strong increase in anticorrelation motions between the two proteins. These data confirm that missense variants of PMPCA variably affect either the degree of rigidity of the GRL and/or the global stability of the protein, likely resulting in impaired substrate binding or catalytic activity of the complex. However, our findings argue against the existence of reliable genotype-phenotype correlations, suggesting that the phenotypic manifestations of PMPCA-related disease are likely influenced by modifier factors, which remain to be disclosed.