2. Materials and Methods
2.1. Materials
H2DCFDA, MitoSOX Red, MitoTracker Red CM-H2XRos were purchased from Invitrogen; DMEM, FBS, Glutamax, Sodium Pyruvate and TrypLE™ Express were purchased from Gibco, USA; HBSS was purchased from PanEco, Russia.
2.2. Cell Lines
All cell cultures were kindly provided by prof. Andrey Y. Abramov (UCL Institute of Neurology, London, Great Britain). Cell cultures of human skin fibroblasts carrying mutations in the genes encoding alpha-synuclein (a-syn, A53T het), PINK1 (homozygous p.Try90Leufsx12), PINK1 and Parkin (PINK1/Parkin, PARK2 R275W/WT+PINK1 p.Try90Leu fs*12/WT), and G2019S mutation in LRRK2 gene, as well as control lines of fibroblasts were used as objects in this study (
Table 1). We used the same cell lines as in the previous study [
14].
Cells were cultured on 25 cm2 culture flasks in DMEM containing 10% FBS (Invitrogen, USA), 2 mM glutamine, 1 mM pyruvate at 37 °C, 5% CO2, and 100% humidity. Upon reaching 80–85% confluence, the cells were splitted to maintain the culture, or seeded on round 25 mm glass coverslips for the experiment. All cells used for experiments were not older than 18 passage.
2.3. Hydrogen Peroxide-Induced Stress
Hydrogen peroxide was added to the incubation culture medium at 37 degrees for 1 hour, after which it was washed three times and left to rest for another hour, after which the experiment was started. The concentration of hydrogen peroxide was selected from the calculation of the dose leading to the death of 20% of cells for control fibroblast lines.
2.4. ROS Measurements
To assess cytosolic ROS production, the cells were loaded with the H2DCFDA probe (10 μM, 40 min) followed by washing with HBSS. Fluorescence registration was carried out using a Spark10M tablet reader (Tecan Group, Switzerland) or an imaging station based on a Leica DMI6000 B inverted microscope (Leica Microsystems, Germany) using a standard FITC filter set (excitation: 494 ± 10 nm, emission: 535 ± 10 nm) with a 20x objective lens. Recording speed 1 frame/5 sec. To assess ROS production in mitochondria, cells were loaded with a MitoSOX Red probe (5 μM, 15 min) or MitoTracker Red CM-H2XRos (100 nM, 30 min) followed by washing with HBSS. Registration was performed using a standard Texas Red filter set (excitation: 575 ± 10 nm, registration: 624 ± 20 nm) and a 20x objective lens. Recording speed 1 frame/5 sec.
2.5. Mitochondrial Membrane Potential Measurements
The mitochondrial membrane potential (Δψm) was measured by incubating cells with 25 nM tetramethylrhodamine methyl ester (TMRM) fluorescent dye in a buffered saline solution (HBSS) for 40 minutes at room temperature. Fluorescent images were obtained on Zeiss 900 CLSM (Carl Zeiss Microscopy GmbH, Jena, Germany) confocal microscope equipped with a × 63 oil immersion objective and during the measurements, the TMRM remained in the HBSS solution. Illumination intensity was kept to a minimum (0.1–0.2 % of laser output) to avoid phototoxicity. The cells were excited with a laser at 561 nm, and the fluorescence was detected above 580 nm. FCCP, which is uncouplers of oxidative phosphorylation, were used to assess the mitochondrial function. When FCCP is added, TMRM fluorescence is quenched. The mitochondrial membrane potential was estimated as the difference in TMRM fluorescence (maximum signal minus the signal with the addition of 2 μM FCCP) and was taken as 100% in control fibroblasts without treatment.
2.6. Mitochondrial Network Morphology Analysis
To analyze the MN morphology, images of cells loaded with TMRM (25 nM), which constantly remained in the working solution, were taken on Zeiss 900 CLSM (Carl Zeiss Microscopy GmbH, Jena, Germany) confocal microscope equipped with a × 63 oil immersion objective. Resolution is 23.4 pixels per micron. The mitochondrial area may depend on the cell volume. In order to take into account this effect on the MN area, we have previously estimated the volume of cells using Calcein AM fluorescence. Calcein at a concentration of 5 μM was loaded into cells for 30 minutes at room temperature, after which Z-stacks images were obtained on a confocal microscope and the fluorescence area was calculated. It was found that there is no difference in this parameter in the studied cell lines. Next, we used a self-written Fiji plugin for automatic batch processing. The protocol for assessing the MN morphology was based on the approaches described earlier [
12,
15,
16]. Detailed analysis steps are shown in
Figure 2.
1. At the first stage, the image is preprocessed, which included background subtraction (rolling ball radius set on 50 pixels), median filter and local contrast enhance (CLACHE) functions. 3. Next, the image was binarized using threshold by “Otsu” method. 3.1 Next, the image was copied and the mitochondrial footprint was calculated by the sum of the mitochondrial areas in the “Particle analyzer” Fiji plugin. 4. Further, the image was skeletonized using the “Skeletonize” function and processed using the built-in “Skeleton Analysis” function. 5. Next, the table of the received data was transmitted to the R-programming environment. The data table was filtered using the “filter” function from “dplyr” package into two data subsets: “mitochondrial network” and “individual mitochondria” according to the described scheme (mitochondrial network were considered as the object what consist of more one branch and non-zero junction, whereas individual mitochondria consist of only one branch and zero junction).
2.7. Statistical Analysis
Computer processing and data images analysis of cell cultures were carried out using Fiji and R-Studio. OriginPro2018 was used for plotting. For experiments with the assessment of mitochondrial membrane potential and the rate of cytosolic and mitochondrial ROS production, two control lines of fibroblasts were combined and taken as 100%. When processing data on the measurement of ROS production, curves were obtained, then were linearly approximated, and the rate of the fluorescent signal increase was calculated. Statistical analysis was performed using “dplyr” package and OriginPro2018 programs using one-way ANOVA (post Tukey test with Bonferroni correction) parametric analysis. The results were presented as mean ± standard deviation of the mean (SD). Differences were considered statistically significant at: * p < 0.05, ** p < 0.01.
4. Discussion
This work is a quantitative analysis of the mitochondrial network (MN) morphology (parameters such as: mitochondrial footprint, length of mitochondria in NW, the ratio of individual mitochondria and their length, as well as mitochondrial connectivity) in fibroblasts obtained from patients with hereditary forms of PD and established mutations, and also their response to hydrogen peroxide-induced stress. Two fibroblast lines studied in this work are associated with a mitophagy-related genes PINK1 and PARK2 malfunction and one line contains a mutation in a Leucine-rich repeat kinase 2 (LRRK2), a protein involved in autophagy. Another fibroblast line contains a point mutation in the gene encoding α-synuclein (A53T, SNCA), the aggregates of which are found in the brains of PD patients. As a control, we used fibroblasts from conditionally healthy individuals comparable in age to patients with PD, between which in most cases there were no significant differences. This fact is important for such an analysis, since cells from donors are highly heterogeneous. In addition to the MN morphology, we evaluated the mitochondrial membrane potential and the rate of cytosolic and mitochondrial ROS production in these fibroblasts. A change in the stress response may indicate pathology. In our study, hydrogen peroxide was used as a stress model, since it is known that it causes ROS-dependent ROS production and in some cases [
17], it has revealed differences in the cells parameters that were not present under normal conditions. It is believed that MN morphology may depend on the external stress on the cell, substrate deprivation, etc., and the rearrangement of MN morphology is the most important for cell adaptation under various conditions [
3,
10,
18].
The MN may be a potential target for the therapy of neurodegenerative diseases, including PD [
19]. It is shown that with aging, mitochondria undergo fragmentation and lose area [
20], which is also may be associated with the development of neurodegeneration [
21]. However, in this work we showed a significant change in the mitochondrial footprint only for the Pink/Parkin mutation fibroblasts, which, on the contrary, led to an increase in this area of MN (
Figure 3b), probably, due to a violation of mitophagy. Interestingly, when exposed to stress, there was indeed a decrease in the mitochondrial footprint of fibroblasts, but for the most part it was not significant and mitochondrial footprint in Pink/Parkin mutation fibroblasts remained the largest. The change in the mitochondria length may occur due to increased fusion or disruption of fission. An increase in the mitochondria length may indicate, for instance, an intensification of the ATP production processes [
22]. In our work, we found an increase in the median mitochondria length of MN (mitochondrial branch) in cells with Pink/Parkin and A53T mutations (
Figure 3c). Interestingly, although it was previously shown that the expression of α-synuclein with A53T mutation in mammalian cells leads to fragmentation of mitochondria, we found the opposite effect [
23]. Since Mfn1 and Mfn2 involved in mitochondrial fusion are known to be substrates for Parkin ubiquitination [
24], it is not surprising that impaired Parkin function leads to an increase in the length of mitochondria due intense fusion. Moreover, fusion of damaged mitochondria with a network has been shown to separate intact mDNA between these mitochondria [
7] and possibly this process took place in the fibroblast in cells with Pink/Parkin and A53T mutations. At the same time, it is worth noting that the peroxide-induced stress led to a significant decrease in the mitochondria length in all the studied cells, except for fibroblasts with the A53T mutation. It is possible that α-synuclein somehow involved in the mitochondrial fission in fibroblasts, which remains to be elucidated in future studies.
The accumulation of individual mitochondria may also indicate a violation of mitochondrial dynamics and mitophagy [
16,
15]. We found that for fibroblasts with Pink/Parkin and PINK1 mutations, the ratio of individual mitochondria/to network ones is higher compared to the control (
Figure 4a), which may indicate increased division or/and inefficiency of mitophagy. Interestingly, despite mutations in Parkin and PINK1 involved in the recruitment of mitochondria into mitophagy 15 [
15,
25,
26], early reports of these cells with mutations showed the levels of mitophagy in fibroblast similar to control cells, studied by double staining of lysosomes and mitochondria [
14]. Perhaps from our results it can be concluded that when analyzing mitophagy, it is necessary to look not only at the colocalized lysosomes and mitochondria, but also the ratio of individual mitochondria in cells. It is noteworthy that under the action of hydrogen peroxide, the ratio of individual mitochondria to mitochondrial networks in fibroblasts with PINK1 cell mutation decreased, whereas for the double Pink/Parkin mutation it did not change. A decrease in the ratio of individual mitochondria may be associated with the activation of transcription factor, which is known to be involved in the launch of mitophagy and is dependent on stress [
27,
28,
29]. Perhaps such stress can cause mitophagy activation, but only in the case of a single PINK1 mutation. This is consistent with data where it was previously shown mild and sustained hydrogen peroxide (H
2O
2) stimulation induces Parkin-dependent mitophagy accompanied by downregulation of the mitophagy-associated proteins OPTN, NDP52, and MFN2 [
31]. This protective effect of moderate stress during exercise or starvation may also play a role in increasing mitochondrial biogenesis via PGC-1α and DRP1 [
22,
31,
32].
In addition to ratio of individual mitochondria we made an attempt to study which size of individual mitochondria mainly represents in fibroblasts. To do this, we estimated the median length of individual mitochondria (
Figure 4b). It was found that in control fibroblasts and fibroblasts with mutations, individual mitochondria are predominantly up to 1 µm in length. Surprisingly, the PINK1 mutation is characterized by shorter individual mitochondria, while the A53T mutation, on the contrary, is longer compared to the control. The impact of stress had an ambiguous effect on the size of individual mitochondria: in control1 and in cells with the Pink/Parkin mutation, they became shorter, although there were no significant changes in other fibroblasts.
The MN connectivity can also change in pathology and as a response to physiological conditions, for example, starvation [
16]. We have shown, that the MN connectivity does not differ in fibroblasts with mutations and control cells under normal conditions (
Figure 5). Interestingly, after hydrogen peroxide addition, the number of junctions in the MN of controls fibroblasts and fibroblasts with the LLRK2 mutation significantly decreases, while for PINK1, A53T and Pink/Parkin fibroblasts it does not significantly change in response to stress. Possibly this may be due to a violation in these cells of the mitochondrial fission machinery, because fission requires the coordinated work of the GTPase separating the inner and outer mitochondrial membranes [
33]. When only the inner membrane separates, the mitochondria can still be connected to the rest of the MN [
13].
According to modern concepts, mitochondrial dysfunction is a key factor in the pathogenesis of hereditary and sporadic PD. The mitochondrial function can be indicated by the potential of the mitochondrial membrane. We evaluated this parameter by TMRM fluorescence (
Figure 6). It was found that fibroblasts with Pink/Parkin and A53T mutations are characterized by a significant increase – hyperpolarization of the mitochondrial membrane, and such a trend is observed for PINK1 mutation. Hyperpolarization of the mitochondrial membrane has been shown for aging cells and may be associated with oxidative stress and also as consequence of increased energy requirements, substrate deprivation or disruption of the ETC [
34,
35]. Interestingly, the cells with hyperpolarized membranes turned out to be sensitive to stress – in the control lines, the addition of hydrogen peroxide did not lead to any significant differences, whereas for A53T and Pink/Parkin and PINK1 cells, a significant decrease in mitochondrial potential occurred.
Oxidative stress is a hallmark of neurodegenerative diseases, in addition, mitochondrial dysfunction also causes oxidative stress [
1]. Though it is believed that oxidative stress in PD is characteristic of brain neurons, we found an increase in the rate of cytosolic ROS production in almost all fibroblasts with mutations (
Figure 7a). Despite we cannot estimate the contribution of any particular form of ROS, using mitochondrial-oriented probes, it is possible to measure whether mitochondria are sources of ROS overproduction in fibroblasts with mutations. Since under normal conditions it was not possible to detect the difference in mitochondrial ROS production between control cells and fibroblasts with mutations, then we applied a model with peroxide-induced stress (
Figure 7b). After the hydrogen peroxide addition, a significant increase in mitochondrial ROS production was detected in fibroblasts with PINK1 and A53T mutations, that may also indicate ETC disruption and mitochondrial dysfunction [
1]. It should also be noted that we did not find an increase in mitochondrial ROS production in fibroblasts with Pink/Parkin mutation, possibly due to a decrease in mitochondrial potential after the stress, which can disrupt the accumulation of fluorescent probe near the mitochondrial membrane.
To summarize, in fibroblasts from patients with PD and established mutations, we found various changes of MN morphology compared to control cells from conditionally healthy donors. Fibroblasts with a double Pink/Parkin mutation have the most different morphology, whose MN is larger and contains longer mitochondria. Despite the fact that the A53T mutation in SNCA gene is not directly related to mitochondrial dynamics, the fibroblasts that contain it also showed significant differences in the MN morphology, similar to Pink/Parkin cells. It is worth noting that a new parameter was introduced in this work – the ratio of individual mitochondria and MN, using which the accumulation of individual mitochondria in cells with the PINK1 and Pink/Parkin mutation was revealed. Remodeling of the MN in response to stress also differed in conditionally healthy control cells and fibroblasts with mutations and, to the greatest extent, in cells with a double Pink/Parkin mutation.
It is important to note that the morphology of the mitochondrial network does not always correlate with function, but may be a consequence of cell adaptation to conditions. In addition, it is necessary to investigate not only the MN morphology, but also its relationship with metabolic signatures and the state of the cell [
36]. Impairment in the mitochondrial membrane potential and the rate of mitochondrial and cytosolic ROS production were also estimated in this work, but we have not established a connection whether this is the cause or consequence of changes in the MN morphology in these fibroblasts.
Surprisingly, cells with the c mutation in the gene encoding LLRK2 kinase showed no significant differences in the studied parameters, although it is believed that this protein is also associated with autophagy and mitochondrial function [
37,
38].
We believe that such violation in fibroblasts with mutations do not have critical consequences as in neurons, since fibroblasts are relatively rapidly dividing cells, whereas neurons live throughout the life of an individual. However, as we have shown, these cells can be a promising and convenient model for studying the relationship between MN morphology and pathology associated with neurodegeneration [
39]. Furthermore, a recent article has shown that the MN can be a target for regulation of stem cell differentiation [
40]. This is of a great interest in light of the fact that neurons can be obtained from the fibroblasts of patients with PD to replace dead cells in the brain.