PMD presents a wide range of symptoms and complications caused by a deficit of myelin deposits in the white matter. Hypotonia, nystagmus, spastic quadriparesis, ataxia, extrapyramidal signs, and hypomyelination are commonly manifested on MRI [15]. The clinical diagnosis of PMD is based on clinical characteristics and MRI brain scans. All three patients in this study met the criteria for PMD, with characteristics including delayed motor milestones, hypotonia, nystagmus, and hypomyelination on MRI. According to the phenotypic classification (PCS) score (Table 1), Pt2 was classified as the connatal type (PCS = 1); correspondingly, Pt1 (PCS = 4) and Pt3 (PCS = 3) were the classical type. Genetic testing showed that Pt1-Pt3 carried a PLP1 duplication. Research has shown that the PLP1 duplication can increase the transcript level of PLP1, resulting in its cellular overexpression [8]. Our study confirmed that PLP1 expression was higher in patient than control fibroblasts (Fig. S2). However, the cellular and molecular mechanisms underlying the phenotype in PMD patients with a PLP1 duplication was unclear.
A study by Ruiz et al. associated PLP1 overexpression with the dysfunction of mitochondrial respiratory functions [9]. Mitochondria play important roles in bioenergy and cell biology pathways, and their predominant physiological function is the generation of ATP through oxidative phosphorylation. The mitochondrial proton circuit across the inner membrane is crucial to mitochondrial respiration. To evaluate mitochondrial respiratory function, mitochondrial proton currents were assessed by the oxygen consumption rate (OCR), which was measured using oxygen-dependent fluorescence quenching [35]. A lower maximum OCR may result from compromised mitochondrial integrity. A higher maximum OCR correlates with acute impairment [36–38]. Brand et al. suggested that a decrease in maximum OCR may be a strong indicator of potential mitochondrial dysfunction [35]. Our results showed that the maximal OCR was lower in fibroblast mitochondria of the PLP1 duplication mutants than the control group (Fig. 1d), which indicated that PLP1 duplication mutants might have a compromised mitochondrial integrity and impaired mitochondrial respiratory function. A previous study obtained similar results in COS-7 cells [39]; they found that compared with nontransfected COS-7 cells, the maximal and basal OCRs were lower in COS-7 cells expressing PLP, but the exact mechanism of mitochondrial dysfunction resulting from the PLP1 duplication remains unclear.
An organelle’s structure determines its function to some degree. To understand mitochondrial dysfunction, we linked dysfunction to changes in morphology. A normal mitochondrial morphology can be classified into tubular, short tubular, fragmented, and large spherical mitochondria (LASMs) [34], which were observed in our experiments (Fig. 4a). The tubular conformation has been reported to be maintained by the mitochondrial contact site and cristae organizing system (MICOS) complex, a core component in cristae and crista junctions. Subunits of the MICOS complex, Sam50, Mic19 and Mic60, constitute a subcellular structure known as the Sam50-Mic19-Mic60 axis, which sustains the structure of cristae. Disruption of this axis results in loss of crista junctions and abnormal mitochondrial morphology, consequently leading to mitochondrial dysfunction [40, 41].
The depletion of Mic60 or Sam50 knockout causes the formation of LASMs. In our study, there were more LASMs in the mitochondria of the PLP1 duplication mutant than the control (Fig. 4b), demonstrating that the mitochondrial structure could be disrupted. In addition to the abnormal formation of LASMs, we also found that the average length of mitochondria was shorter compared with the control (Fig. 4c). The length of mitochondria may be related to mitochondrial fragmentation, which is the first step in selective autophagosome (mitophagosome) formation or the process of mitochondrial fission [42]. A study has shown that in response to hypoxic stress, the outer mitochondrial membrane protein FUNDC1 recruits DNM1L/DRP1 to target the final scission of mitochondria. As a possible result of mitophagosomes, mitochondrial fragmentation could explain our observation of the shorter average length of mitochondria compared with the control (Fig. 4c).
Taken together, we demonstrated that respiratory functions and mitochondrial structural integrity were impaired in PLP1-duplicated fibroblasts, as evidenced by decreases in OCRs, more LASMs, and shorter average mitochondrial lengths. These results corroborated previous findings regarding the role of mitochondrial respiratory functions in PLP1-related diseases. However, previous studies have indicated that PLP1 is expressed in the endoplasmic reticulum (ER). Therefore, we propose a novel hypothesis that overexpression of PLP1 in the ER affects the mitochondria-associated membrane (MAM) structure, resulting in mitochondrial dysfunction and subsequent cytotoxicity of PLP1. To test our hypothesis, we needed to link these results to the morphology and functions of the ER.
First, we compared the ER morphology in control and PLP1-duplicated fibroblasts. The ER consists of the perinuclear ER and the peripheral ER [43–45]. ER sheets are mainly distributed in perinuclear sites, and only a low percentage of ER sheets can be observed in the periphery. Moreover, most ER tubules are located in the periphery. The ratio of sheets to tubules varies in different processes. In yeasts and COS7 cells, ER sheet rates are higher in translocation protein machinery processes, they may be the preferred site of protein translocation [46, 33, 47].
In comparison to the control, we found abnormal distributions of ER tubules and sheets in PLP1 duplication fibroblasts. As shown in Fig. 2b, the perinuclear ER in control subject fibroblasts had closed membrane sheets (or cisterna), and the peripheral ER consisted of a network of sheets and tubules. In contrast, in the PLP1 duplication mutant fibroblasts, ER sheets were observed in both the perinuclear and peripheral regions, and more ER sheets were observed in the perinuclear region (Fig. 2c). In addition, the ER sheets exhibited an ‘expanded’ morphology (Fig. 2d). Furthermore, PLP1 mainly colocalized with the peripheral ER (Fig. 2a), where more ER sheets were located. A study also demonstrated that the peripheral ER generates expanded ER sheets to withstand ER stress [48]. One possible explanation for our observations is that the increasing ER sheet size may provide proteins with more time to fold by avoiding aggregate formation. The increased number of ER sheets in our results might indicate their additional need to address the overexpression of PLP1. Moreover, their expanded morphology suggested that they might be in an abnormal circumstance, for example, protein folding itself. Therefore, the ER sheet distribution and morphology were changed in PLP1 duplication fibroblasts.
To link abnormal ER distributions to mitochondrial dysfunction, we investigated the dynamic structure of MAMs, which serve as a bridge between the ER and mitochondria. ER morphology and location might affect MAM function. A critical function of MAMs is to coordinate Ca2+ leakage from the ER: Ca2+ is first released through IP3R from the ER section that is in close contact with mitochondria and then transported through the mitochondrial Ca2+ uniporter (MCU) into mitochondria (Fig. 5) [35]. Therefore, MAM function could be affected by MAM distance, which reflects the expansion of the ER-mitochondrion interfaces. The interface explains the amount of overlap between the ER and mitochondria: the higher the level of the interface, the closer is the distance of the MAMs.
For the abnormal ‘expanded ER sheets’ found in the peripheral ER, we investigated the interface between expanded ER sheets and mitochondria. We used the colocalization coefficient to calculate the extent of the interface between the ER and mitochondria. Greater colocalization coefficients (Fig. 3b) were found near the expanded ER sheets and mitochondria in the PLP1 duplication fibroblasts, which indicated a closer distance between expanded ER sheets and mitochondria. The close distances between the ER and mitochondrial membranes might lead to continuous mitochondrial Ca2+ uptake during background Ca2+ release, which in turn could cause mitochondrial Ca2+ overloading and membrane permeabilization. Moreover, this pathway might trigger IP3-induced Ca2+ mobilization from the ER to mitochondria, which has a function in apoptosis [49–51].
Somayajulu et al. demonstrated that increased insertion of PLP into the mitochondria affects oxidative phosphorylation and leads to mitochondrial dysfunction [39]. By demonstrating the overexpression of PLP1 in COS-7 cells and in Plp1 transgenic mice, the authors showed that PLP1 translocates to mitochondria by specific cysteine motifs in the N-terminus. However, they did not explore the role of ER in mitochondrial dysfunction. As the ER is a major organelle for protein translation and processing, our study provides novel insight into the association of the ER and mitochondria through MAM function. In Fig. 5, we summarize the pathways in both healthy and PLP1 duplication fibroblasts: overexpression of PLP1 mainly colocalizes with the ER, leading to abnormal ER structures; next, MAM dysfunction may trigger mitochondrial Ca2+ overloading and membrane permeabilization, ultimately leading to mitochondrial dysfunction.
The novel pathway was derived by comparing control and patient fibroblasts, and the differences between the patients’ cells were not statistically significant, while the patients were representative of the two PMD subtypes. For the mitochondrial respiration function analysis, an abnormal ER and mitochondria morphology or mitochondrial length, there were no significant differences among the three patients. Fibroblasts of Pt2, who was diagnosed with the connatal type (the most severe type on clinical examination), had the lowest level of mitochondrial maximum respiration (Fig. 1d), the shortest mitochondrial length (Fig. 4c), and the highest ER sheet proportion in the periphery among the three patients’ cells (Fig. 2c). Although there were no significant differences among these three patients (Fig. 1d, 2c, 2d, 3b, 4c), these trends may be indicative of mitochondrial dysfunction and PMD clinical severity. The goal of this paper was to investigate the ER-mitochondrion interface in healthy and PLP1-duplicated fibroblasts; the connection between the risk of mitochondrial dysfunction and clinical severity will need to be explored in future studies with larger sample sizes.
As a major limitation, our control samples were from a 27-year-old healthy adult, which were not age- and number-matched to our experimental samples of the three children PMD patients. The primary concern was that sampling fibroblasts was an invasive procedure. While the sampling procedure was approved by our hospital’s ethics committee, it required sampling from the upper limb, possibly leading to bleeding and scar formation. While a study demonstrated that mitotic activity of fibroblasts was independent of the subject age [52] by measuring ratio of in vivo proportion of progenitor fibroblasts (MF) and functional postmitotic fibrocytes (PMF), the age-related cell stress can affect the ER and mitochondrial morphology, further altering the MAMs function. Consequently, age-related cell stress can affect the distributions of mitochondrial morphology (Fig. 4a). Therefore, for future work, biomarkers that measure the growth, differentiation, and apoptosis of fibroblasts such as the pH 6 activity of \(\beta\)-galactosidase or TGF- \(\beta 1\) will be needed [53]. For another, in this study, we demonstrated that in PLP1 overexpressed cells, closer ER-mitochondria interfaces led to pathological mitochondrial morphology, as mediated by MAMs’ structure change (Fig. 3). To detail the roles of MAM proteins, future work will be needed to study MAMs-related proteins such as FUNDC1 [42]. In addition, we proposed that mitochondrial Ca2+ overloading could underlie mitophage as a result of PLP1 overexpression (Fig. 5). To investigate this pathway, the formation of ER–mitochondria Ca2+ tunnel proteins [22] will need to be examined. An example protein is IP3,which interacts with the voltage-dependent anion-selective channel across the ER and mitochondria.