Fragmentation of the mitochondrial network in skin in vivo

Abstract

Mitochondria form dynamic networks which adapt to the environmental requirements of the cell. We investigated the aging process of these networks in human skin cells in vivo by multiphoton microscopy. A study on the age-dependency of the mitochondrial network in young and old volunteers revealed that keratinocytes in old skin establish a significantly more fragmented network with smaller and more compact mitochondrial clusters than keratinocytes in young skin. Furthermore, we investigated the mitochondrial network during differentiation processes of keratinocytes within the epidermis of volunteers. We observe a fragmentation similar to the age-dependent study in almost all parameters. These parallels raise questions about the dynamics of biophysical network structures during aging processes.

About the Authors:

Daniel Mellem

Affiliations Center for Free-Electron Laser Science (CFEL), University of Hamburg, Hamburg, Germany, Beiersdorf AG, Applied Biophysics, Hamburg, Germany

Martin Sattler

Affiliation: Beiersdorf AG, Applied Biophysics, Hamburg, Germany

Sonja Pagel-Wolff

Affiliation: Beiersdorf AG, Applied Biophysics, Hamburg, Germany

Sören Jaspers

Affiliation: Beiersdorf AG, Applied Biophysics, Hamburg, Germany

Horst Wenck

Affiliation: Beiersdorf AG, Applied Biophysics, Hamburg, Germany

Michael Alexander Rübhausen

Affiliation: Center for Free-Electron Laser Science (CFEL), University of Hamburg, Hamburg, Germany

Frank Fischer

* E-mail: [email protected]

Affiliation: Beiersdorf AG, Applied Biophysics, Hamburg, Germany

Introduction

Cells in human tissue produce their energy in form of adenosine triphosphate (ATP) basically by oxidative phosphorylation. This process takes place in mitochondria, that form highly dynamic networks adapting fast to the cells’ environment and their metabolic requirements. [1] In recent years, numerous experiments have been performed to investigate correlations between mitochondrial network states and corresponding metabolic processes within cells. Cancer cells lacking of a glycolytic medium establish fusion states of the mitochondrial network correlating with a change in energy production going from glycolysis to oxidative phosphorylation. [2] In contrast, a transfer from respiratory to glycolytic conditions leads to a fragmentation of the mitochondrial network. [3] During moderate stress mitochondrial networks form hyperfusion states which coincide with increased ATP production. [4] On the contrary, high stress levels induce to a fragmentation of networks. [5] Mitochondrial fission and fusion states are considered to be quality saving mechanisms of the cell. Fission states support repelling of heavily damaged mitochondria from the network, e.g. during autophagy and mitophagy. [6, 7] In contrast, fusion states help to compensate for defect protein complexes or rare metabolites among mitochondria. [8, 9] Recently, a fragmentation of the mitochondrial network with age was observed in vitro. [10] The interplay of mitochondrial dynamics was simulated biophysically in a probabilistic quality model which also revealed a significant fragmentation of the mitochondrial network during aging. [11]

In this paper, we present investigations of mitochondrial morphologies during aging and differentiation of keratinocytes in human skin in vivo. We performed two studies: The first study explored variations of the mitochondrial network in epidermal keratinocytes of young and old volunteers. The second study extended the results by investigating the role of the mitochondrial network during the epidermal turnover.

Materials and methods

Tissue can be investigated non-invasively in vivo by imaging endogenous fluorophores using multiphoton microscopy. [12–14] An important metabolic fluorophore in the epidermis is nicotinamide adenine dinucleotide (NADH) [15, 16] which serves as an electron carrier from the Krebs cycle to respiratory chain. Thus, NADH is almost exclusively located in mitochondria. We tested this premise by in vitro investigations using a confocal scanning microscope (Leica TCS SP5, Leica Microsystems, Mannheim, Germany). Neonatal human epidermal keratinocytes were purchased from Lonza Group AG (Basel, Switzerland) and cultured in Dulbecco’s Modified Eagle Medium (DMEM) (Life Technologies, Carlsbad, USA), supplemented with 10% fetal calf serum (PAA Laboratories, Pasching, Austria), L-glutamine and penicillin/streptomycin (both: Life Technologies, Carlsbad, USA). During culturing cells were maintained in a humidifying incubator with a 5% CO₂ atmosphere at 37°C. Mitochondria in the cells were marked with MitoTracker Red CMXRos (Life Technologies, Carlsbad, USA). We correlated the fluorescence of the dye from 570nm to 650nm to the autofluorescence of keratinocytes in the NADH emission spectrum from 410 nm to 540 nm (Fig 1). The measurements were performed with a 63x objective (HCX PL APO lambda blue 63.0x1.20 WATER UV, Leica). Single images of 2048*2048 pixels were generated with a bidirectional scan either at zoom 3.0 (82μm x 82 μm, step size: 0.04 μm) with scanning rate of 100Hz or at zoom 7.9 (31.2μm x 31.2μm, step size: 0.04μm) with scanning rate of 200Hz.

[Figure omitted. See PDF.]

Fig 1.

Fluorescence of dyed mitochondria (a) and autofluorescence of NADH (b) in keratinocytes (zoom 3.0). Fluorescence of dyed mitochondria (c) and autofluorescence of NADH (d) in a single keratinocyte (zoom 7.9).

https://doi.org/10.1371/journal.pone.0174469.g001

We found a nearly perfect agreement of morphologies of stained mitochondria with the autofluorescence signal of NADH. Aberrations are possibly caused by an overstaining of the dye, leading to very high intensities and consequently to agglomerations of mitochondrial morphologies in the micrograph. Intensity deviations among autofluorescent mitochondrial structures are assumed to be fluctuations in the metabolic activity. Previous studies indicate that the concentration of NADH is five times higher [17] and the quantum yield is 1.5 to 2.5 greater than of coenzyme NADPH [18] that is fluorescing in the same spectral range. So, potential interference of NADH with the coenzyme NADPH can be neglected.

Imaging mitochondrial network morphologies in vivo by the autofluorescence of NADH was performed with a multiphoton microscope (Dermainspect, developed by Jenlab GmbH (Jena, Germany) in collaboration with Beiersdorf AG (Hamburg, Germany) (S1 Fig). The device includes a femtosecond laser with a repetition rate of 80MHz and a pulse width of 150fs. Laser pulses excite the autofluorescence of NADH in the epidermis at a wavelength of 750nm. The autofluorescence is separated from the laser light by a 548nm ± 150nm band-pass filter (HQ 548/305m-2P, Schott AG, Mainz, Germany) and detected by a photomultiplier tube (PMT). In order to image a single cell layer in the epidermis, the focus of the laser through an objective (40x magnification, 1.3 numerical aperture) is scanned parallel to the surface of the skin with an image acquisition time of about 12,5s. During the measurement the z-position is controlled by a piezo element. To prevent disturbance in the images by the movement of the volunteer, the objectives adheres to the skin of the volunteers. Each image has the dimensions of 110μm x 110μm and a resolution of 512pixels x 512pixels. In all measurements the laser power after the objective was set to 20 mW. For the quantitative analyses of mitochondrial morphologies we developed an image binarization algorithm that extracted mitochondrial structures from multiphotonic measurements: Firstly, the cytoplasm of each keratinocyte was taken from the measurements by manually annotating the cell borders and the corresponding nuclei in all images. Then, an Otsu-based [19] image algorithm extracted mitochondrial morphologies using binarization (S1 Fig). A pixel in one annotations had to exceed two thresholds in order to be categorized as NADH-autofluorescence during the binarization process: A global threshold which factors the gray values of all pixels in the annotation and a local threshold which is defined by the gray value of the surrounding pixels. Both of these thresholds were calculated via the Otsu-method. The algorithm calculates gradually all possible binarizations and selects the gray value as the optimal threshold that minimizes the Otsu variance:(1)

In the Eq (1) T represents the threshold of the gray value for the binarization, σ_ω(T) its Otsu variance, the variances of the gray values of black and white pixels and ω_i(T) their corresponding weights. After binarization procedures, the algorithm assigns all signal related pixels to common signal clusters, so that they fulfill the definition of 4-neighbor-connectivity. The number of pixels per cluster defines its area A, the number of edges which border on black pixels (the background of the measurement) defines its circumference U and the parameter C = U²/4πA defines its circularity. A represents the volume of a mitochondrial structure, C its circularity. Mann-Whitney-Tests were used for statistical analysis (p≤0.05:*; p≤0.01:**; p≤0.001:***). The analysis was performed with the software “Statistica” by StatSoft (Europe) GmbH (Hamburg, Germany). All corresponding plots were depicted with the software OriginPro 8 (OriginLab Corporation, Northampton, USA).

Two clinical in vivo studies were performed. Written informed consent was obtained from each volunteer in both studies. All volunteers provided skin without pathologies and possessed skin types II and III of the Fitzpatrick phototyping scale. [20] Prior and during the studies all volunteers were required to desist from visits of solariums and intensive sun exposure. Treatment of skin in the investigated areas with cosmetic substances and medicals was prohibited during the studies. Measurements in four areas at the inner side of the forearm were performed. In each area multiphotonic images were acquired at three positions. The results of all positions were averaged in each area. The examinations were performed by trained personnel at standard atmospheric conditions (23°C ± 1°C and 43% ± 2% relative humidity). All volunteers were provided with fifteen minutes of acclimatization prior to the measurements. The studies were conducted according to the recommendations of the current version of the Declaration of Helsinki and the Guideline of the International Conference on Harmonization Good Clinical Practice, (ICH GCP). In addition, this study was approved and cleared by the institutional ethics review board (Beiersdorf AG, Hamburg, Germany). Written informed consent was obtained from each volunteer in both studies.

Results

In the first study the mitochondrial morphology in epidermal keratinocytes of twelve young (mean ± SD: 23.75 ± 1.67 years) and twelve old (72.17 ± 4.15 years) volunteers was analysed by statistically comparing 48 areas in each age group. We analysed the number, area and circularity of the mitochondrial clusters per keratinocyte in the stratum granulosum in a depth of about 15 μm. The results of all examined parameters were averaged in each area before statistical comparison. Mitochondrial clusters in the stratum granulosum of young volunteers are significantly larger (p = 0.029) and have a significantly higher circularity (p = 0.014) (Fig 2). In contrast, the number of mitochondrial clusters normalized to the size of each keratinocyte is significantly higher (p = 0.006) in old skin (S2 Fig).

[Figure omitted. See PDF.]

Fig 2.

(a) Comparison of mitochondrial networks in young volunteers and old volunteers before and after binarization. Age-related comparison of size (b) and circularity (c) of mitochondrial clusters in keratinocytes. The p-value represents the probability if the two compared probability distributions derive from the same population. The square in the box represents the mean value, the horizontal line in the center the median. 25% percent of the measured values are lower than the bottom of the box, 75% of them are lower than the top of the box. The antennas (whiskers) are the corresponding limits for 2.5% and 97.5%. Crosses represent the highest and the lowest value in the distribution.

https://doi.org/10.1371/journal.pone.0174469.g002

In the second study the morphology of the mitochondrial networks in keratinocytes of the stratum spinosum (depth of about 25 μm) and the stratum granulosum (depth of about 15 μm) were analysed by statistically comparing 32 areas in each layer. For determination of the layers, a prescan through the epidermis was performed during which the mosaic structure of the spinolar cells and the larger and rounder granular cells were identified by trained personnel. All eight volunteers were of the same age group (62 ± 1.31 years). Again, investigated parameters were averaged in each area. The average size per mitochondrial cluster is significantly larger (p ≤ 0,001) and its circularity significantly higher (p ≤ 0,001) in the stratum spinosum (Fig 3). Contrary, the number of mitochondrial clusters normalized to the cell size tends to be higher (p ≤ 0,086) in the stratum granulosum (S3 Fig).

[Figure omitted. See PDF.]

Fig 3.

(a) Comparison of mitochondrial networks in stratum spinosum and stratum granulosum representing different states during the differentiation process of keratinocytes. Differentiation-related comparison of size (b) and circularity (c) of mitochondrial clusters in keratinocytes. The p-value represents the probability if the two compared probability distributions derive from the same population. The square in the box represents the mean value, the horizontal line in the center the median. 25% percent of the measured values are lower than the bottom of the box, 75% of them are lower than the top of the box. The antennas (whiskers) are the corresponding limits for 2.5% and 97.5%. Crosses represent the highest and the lowest value in the distribution.

https://doi.org/10.1371/journal.pone.0174469.g003

Discussion

Our study of the mitochondrial network during aging reveals, that granular epidermal keratinocytes in young skin have less clusters in total than keratinocytes in old skin. Additionally, clusters in granular keratinocytes of young volunteers are bigger in size and less circular in shape. These findings reveal a high connectivity among mitochondria in the stratum granulosum in young skin and a fragmented mitochondrial network in the stratum granulosum of old skin. Granular keratinocytes of young volunteers seem to prefer a fusion state of the mitochondrial network which is not maintained during aging.

In our second study we investigated the changes of mitochondrial morphologies during epidermal differentiation. At the beginning of the epidermal turnover proliferation of keratinocytes occurs in the deepest layer of the epidermis, the stratum basale. During their differentation process keratinocytes migrate to upper layers, i.e. the stratum spinosum and the stratum granulosum, and finally transform into dead corneocytes in the stratum corneum. [21] For an analysis of mitochondrial network states during differentiation, we performed examinations in two layers above the stratum basale: in the stratum spinosum and the stratum granulosum. The stratum spinosum is the layer next to the stratum basale, so that keratinocytes in the stratum spinosum are less differentiated than in the stratum granulosum which is the last living epidermal layer before cornification occurs. We observed that the number of mitochondrial clusters tends to increase during differentiation. Additionally, clusters significantly shrink to smaller volumes and more compact shapes. These parameters indicate that mitochondrial networks in the stratum spinosum establish fused states which fragment during differentiation processes.

Comparing the results of the two studies the dynamics in mitochondrial network morphologies during differentiation agree well with changes during skin aging in volunteers. In the following we refer to several processes which play a role in the fragmentation of the mitochondrial network.

Reactive oxygen species (ROS) are produced due to long-term solar UV radiation [22] as well as intrinsic generation by mitochondria and other cell organelles. [23] Oxidative stress causes long-term damage to intramitochondrial structures [24] and leads to mitochondrial fission states during aging. [25] Numerous experiments point to the importance of oxidative stress for differentiation processes of keratinocytes. [26] The exclusion of heavily ROS-damaged mitochondria from the mitochondrial network could thereby lead to a fragmentation of the network’s morphology during aging and differentiation.

Mitochondrial aging processes lead to a quality decay of mitochondria [27] which correlates with an age-related decrease of mitochondrial membrane potentials (MMP). [28] The decrease of the MMP also plays a role during the differentiation of keratinocytes. [21] Experiments revealed that mitochondria with low MMPs are excluded from the network by autophagic processes. [29]

Mitophagic and autophagic processes exclude heavily damaged mitochondria from the network. Mitochondrial autophagy increases during the aging of cells [30], thereby fragmenting intermitochondrial connections. Autophagic activity is similarly involved in the differentiation of keratinocytes. [31]

In cultivated keratinocytes an increased glycolytic activity in old volunteers was observed. [32] This shift from oxidative phosphorylation to glycolysis during aging correlates with a fragmentation of the mitochondrial network. Recent publications disclose the change of mitochondrial metabolisms during the differentiation process. [33]

There are several alternatives to explain the similarities of mitochondrial morphologies during aging and differentiation. For instance, keratinocytes could exploit natural aging processes to catalyze differentiation procedures which consequently leads to a fragmentation of the mitochondrial network. Another possible explanation is that mitochondrial morphologies are independent from differentiation and purely influenced by aging processes during the epidermal turnover. Further investigations are necessary to clarify the physiological connections between aging and differentiation regarding mitochondrial morphologies in cells. The influence of aging processes on mitochondrial quality was investigated in computational quality models. [11, 34] A fragmentation of the mitochondrial network during aging of cells and a deceleration of fission and fusion cycles as a quality saving process were observed. The models reveal quality saving mechanisms for mitochondria such as the interplay of mitochondrial networking and recycling and the quality saving benefit of a decreased mitochondrial repair ability in old cells. Hence, our findings are qualitatively in good agreement with the model.

Conclusion

In conclusion, we have measured and analyzed the morphology of mitochondrial networks in human skin in vivo for the first time. We found a significant fragmentation of mitochondrial morphologies in granular keratinocytes during aging. Our results are qualitatively in good agreement with age-dependent investigations in vitro and with simulations of biophysical models. Moreover, we observed a fragmentation of the network during the differentiation process of keratinocytes in the epidermis. The parallels of both studies raise questions about the linkage of aging and differentiation concerning mitochondrial morphologies.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. (a) setup of multiphoton microscope (b) measurement of autofluorescence of NADH in stratum granulosum (c) binarization of measured keratinocyte with otsu algorithm.

https://doi.org/10.1371/journal.pone.0174469.s001

(TIF)

S2 Fig. Comparison of number of mitochondrial clusters in young and old volunteers normed to the size of keratinocytes.

https://doi.org/10.1371/journal.pone.0174469.s002

(TIF)

S3 Fig. Comparison of number of mitochondrial clusters in stratum spinosum and stratum granulosum normed to the size of keratinocytes.

https://doi.org/10.1371/journal.pone.0174469.s003

(TIF)

S1 Dataset. The .csv-file provides the data of mitochondrial cluster numbers, cluster size and cluster circularity in all 40 areas that were investigated in the skin of young and old volunteers and in the stratum spinosum and the stratum granulosum.

https://doi.org/10.1371/journal.pone.0174469.s004

(CSV)

Author Contributions

1. Conceptualization: FF DM.

2. Data curation: DM SPW MS FF.

3. Formal analysis: DM SPW MS FF.

4. Investigation: DM SPW MS.

5. Methodology: DM MS SPW FF.

6. Project administration: FF.

7. Resources: DM SPW MS FF.

8. Software: DM.

9. Supervision: FF HW SJ MR.

10. Validation: DM SPW MS FF.

11. Visualization: DM.

12. Writing – original draft: DM.

13. Writing – review & editing: DM FF MR HW SJ.

Citation: Mellem D, Sattler M, Pagel-Wolff S, Jaspers S, Wenck H, Rübhausen MA, et al. (2017) Fragmentation of the mitochondrial network in skin in vivo. PLoS ONE 12(6): e0174469. https://doi.org/10.1371/journal.pone.0174469

References

1. Youle RJ, van der Bliek AM. Mitochondrial fission, fusion, and stress. Science. 2012;337(6098):1062–5. pmid:22936770

2. Rossignol R, Gilkerson R, Aggeler R, Yamagata K, Remington SJ, Capaldi RA. Energy substrate modulates mitochondrial structure and oxidative capacity in cancer cells. Cancer Res. 2004;64(3):985–93. pmid:14871829

3. Jakobs S, Martini N, Schauss AC, Egner A, Westermann B, Hell SW. Spatial and temporal dynamics of budding yeast mitochondria lacking the division component Fis1p. J Cell Sci. 2003;116(Pt 10):2005–14. pmid:12679388

4. Tondera D, Grandemange S, Jourdain A, Karbowski M, Mattenberger Y, Herzig S, et al. SLP-2 is required for stress-induced mitochondrial hyperfusion. EMBO J. 2009;28(11):1589–600. pmid:19360003

5. van der Bliek AM, Shen Q, Kawajiri S. Mechanisms of mitochondrial fission and fusion. Cold Spring Harb Perspect Biol. 2013;5(6). pmid:23732471

6. Singh SB, Ornatowski W, Vergne I, Naylor J, Delgado M, Roberts E, et al. Human IRGM regulates autophagy and cell-autonomous immunity functions through mitochondria. Nat Cell Biol. 2010;12(12):1154–65. pmid:21102437

7. Kim I, Rodriguez-Enriquez S, Lemasters JJ. Selective degradation of mitochondria by mitophagy. Arch Biochem Biophys. 2007;462(2):245–53. pmid:17475204

8. Nakada K, Inoue K, Ono T, Isobe K, Ogura A, Goto YI, et al. Inter-mitochondrial complementation: Mitochondria-specific system preventing mice from expression of disease phenotypes by mutant mtDNA. Nat Med. 2001;7(8):934–40. pmid:11479626

9. Rolland SG, Motori E, Memar N, Hench J, Frank S, Winklhofer KF, et al. Impaired complex IV activity in response to loss of LRPPRC function can be compensated by mitochondrial hyperfusion. Proc Natl Acad Sci U S A. 2013;110(32):E2967–76. pmid:23878239

10. Regmi SG, Rolland SG, Conradt B. Age-dependent changes in mitochondrial morphology and volume are not predictors of lifespan. Aging (Albany NY). 2014;6(2):118–30. pmid:24642473

11. Figge MT, Reichert AS, Meyer-Hermann M, Osiewacz HD. Deceleration of fusion-fission cycles improves mitochondrial quality control during aging. PLoS Comput Biol. 2012;8(6):e1002576. pmid:22761564

12. Lakowicz JR, Szmacinski H, Nowaczyk K, Johnson ML. Fluorescence lifetime imaging of free and protein-bound NADH. Proc Natl Acad Sci U S A. 1992;89(4):1271–5. pmid:1741380

13. Skala MC, Riching KM, Gendron-Fitzpatrick A, Eickhoff J, Eliceiri KW, White JG, et al. In vivo multiphoton microscopy of NADH and FAD redox states, fluorescence lifetimes, and cellular morphology in precancerous epithelia. Proc Natl Acad Sci U S A. 2007;104(49):19494–9. pmid:18042710

14. Puschmann S, Rahn CD, Wenck H, Gallinat S, Fischer F. Approach to quantify human dermal skin aging using multiphoton laser scanning microscopy. J Biomed Opt. 2012;17(3):036005. pmid:22502563

15. Mayevsky A, Rogatsky GG. Mitochondrial function in vivo evaluated by NADH fluorescence: from animal models to human studies. Am J Physiol Cell Physiol. 2007;292(2):C615–40. pmid:16943239

16. Balu M, Mazhar A, Hayakawa CK, Mittal R, Krasieva RB. In vivo multiphoton NADH fluorescence reveals depth-dependent keratinocyte metabolism in human skin. Biophys J. 2013;104(1):258–67. pmid:23332078

17. Klaidman LK, Leung AC, Adams J J D. High-performance liquid chromatography analysis of oxidized and reduced pyridine dinucleotides in specific brain regions. Anal Biochem. 1995;228(2):312–7. pmid:8572312

18. Avi-Dor Y, Olson JM, Doherty MD, Kaplan NO. Fluorescence of Pyridine Nucleotides in Mitochondria. J Biol Chem. 1962;237:2377–2383.

19. Otsu N. A Threshold Selection Method from Gray-Level Histograms. Systems, Man and Cybernetics, IEEE Transactions. 1979;9(1):5.

20. Fitzpatrick TB. Soleil et peau. Journal de Médecine Esthétique. 1975;2:2.

21. Allombert-Blaise C, Tamiji S, Mortier L, Fauvel H, Tual M, Delaporte E, et al. Terminal differentiation of human epidermal keratinocytes involves mitochondria- and caspase-dependent cell death pathway. Cell Death Differ. 2003;10(7):850–2. pmid:12815468

22. D’Orazio J, Jarrett S, Amaro-Ortiz A, Scott T. UV Radiation and the Skin. Int J Mol Sci. 2013;14(6):12222–48. pmid:23749111

23. Brown GC, Borutaite V. There is no evidence that mitochondria are the main source of reactive oxygen species in mammalian cells. Mitochondrion. 2012;12(1):1–4. pmid:21303703

24. Yakes FM, Van Houten B. Mitochondrial DNA damage is more extensive and persists longer than nuclear DNA damage in human cells following oxidative stress. Proc Natl Acad Sci U S A. 1997;94(2):514–9. pmid:9012815

25. Wu S, Zhou F, Zhang Z, Xing D. Mitochondrial oxidative stress causes mitochondrial fragmentation via differential modulation of mitochondrial fission-fusion proteins. FEBS J. 2011;278(6):941–54. pmid:21232014

26. Hamanaka RB, Chandel NS. Mitochondrial metabolism as a regulator of keratinocyte differentiation. Cell Logist. 2013;3(1):e25456. pmid:24475371

27. Rooyackers OE, Adey DB, Ades PA, Nair KS. Effect of age on in vivo rates of mitochondrial protein synthesis in human skeletal muscle. Proc Natl Acad Sci U S A. 1996;93(26):15364–9. pmid:8986817

28. Hagen TM, Yowe DL, Bartholomew JC, Wehr CM, Do KL, Park JY, et al. Mitochondrial decay in hepatocytes from old rats: membrane potential declines, heterogeneity and oxidants increase. Proc Natl Acad Sci U S A. 1997;94(7):3064–9. pmid:9096346

29. Twig G, Shirihai OS. The interplay between mitochondrial dynamics and mitophagy. Antioxid Redox Signal. 2011;14(10):1939–51. pmid:21128700

30. Oberley TD, Swanlund JM, Zhang HJ, Kregel KC. Aging results in increased autophagy of mitochondria and protein nitration in rat hepatocytes following heat stress. J Histochem Cytochem. 2008;56(6):615–27. pmid:18379016

31. Aymard E, Barruche V, Naves T, Bordes S, Closs B, Verdier M, et al. Autophagy in human keratinocytes: an early step of the differentiation? Exp Dermatol. 2011;20(3):263–8. pmid:21166723

32. Prahl S, Kueper T, Biernoth T, Wohrmann Y, Munster A, Furstenau M, et al. Aging skin is functionally anaerobic: importance of coenzyme Q10 for anti aging skin care. Biofactors. 2008;32(1–4):245–55. pmid:19096122

33. Hamanaka RB, Glasauer A, Hoover P, Yang S, Blatt H, Mullen AR, et al. Mitochondrial reactive oxygen species promote epidermal differentiation and hair follicle development. Sci Signal. 2013;6(261):ra8. pmid:23386745

34. Mellem D, Jaspers S, Wenck H, Rübhausen MA, Fischer F. Quality Saving Mechanisms of Mitochondria during Aging in a Fully Time-Dependent Computational Biophysical Model. PLoS ONE. 2016;11(1):1–22. pmid:26771181