- Research article
- Open access
- Published:
Deficits in mitochondrial TCA cycle and OXPHOS precede rod photoreceptor degeneration during chronic HIF activation
Molecular Neurodegeneration volume 18, Article number: 15 (2023)
Abstract
Background
Major retinal degenerative diseases, including age-related macular degeneration, diabetic retinopathy and retinal detachment, are associated with a local decrease in oxygen availability causing the formation of hypoxic areas affecting the photoreceptor (PR) cells. Here, we addressed the underlying pathological mechanisms of PR degeneration by focusing on energy metabolism during chronic activation of hypoxia-inducible factors (HIFs) in rod PR.
Methods
We used two-photon laser scanning microscopy (TPLSM) of genetically encoded biosensors delivered by adeno-associated viruses (AAV) to determine lactate and glucose dynamics in PR and inner retinal cells. Retinal layer-specific proteomics, in situ enzymatic assays and immunofluorescence studies were used to analyse mitochondrial metabolism in rod PRs during chronic HIF activation.
Results
PRs exhibited remarkably higher glycolytic flux through the hexokinases than neurons of the inner retina. Chronic HIF activation in rods did not cause overt change in glucose dynamics but an increase in lactate production nonetheless. Furthermore, dysregulation of the oxidative phosphorylation pathway (OXPHOS) and tricarboxylic acid (TCA) cycle in rods with an activated hypoxic response decelerated cellular anabolism causing shortening of rod photoreceptor outer segments (OS) before onset of cell degeneration. Interestingly, rods with deficient OXPHOS but an intact TCA cycle did not exhibit these early signs of anabolic dysregulation and showed a slower course of degeneration.
Conclusion
Together, these data indicate an exceeding high glycolytic flux in rods and highlight the importance of mitochondrial metabolism and especially of the TCA cycle for PR survival in conditions of increased HIF activity.
Background
The retina is an energetically highly active tissue and relies on a well-functioning blood supply for the delivery of sufficient oxygen and nutrients to meet its metabolic needs [1]. The photoreceptor (PR) cells, residing in the outer retina, are especially energy-demanding and have high oxygen and glucose consumption levels [2,3,4,5]. Oxygen reaching the outer retina is mainly used in the photoreceptor inner segments (IS), a specialized cell compartment tightly packed with mitochondria [6]. Here, ATP is produced providing the required energy to maintain ion gradients across the cell membrane by Na+/K+-ATPase and to sustain vision [7, 8]. Despite the high energy demand, isolated retina produces lactate in the presence of oxygen [9,10,11,12]. This incomplete oxidation of glucose under aerobic conditions, or aerobic glycolysis, and the associated high glucose consumption are metabolic adaptations characteristic of anabolic and proliferating cells facilitating the incorporation of nutrients into biomass [13]. Thus, retinal lactate production has been attributed to the highly anabolic PRs possibly supporting the diurnal renewal of the photoreceptor outer segments (OS) [12, 14,15,16]. A growing body of evidence suggests also metabolic coupling of PRs and the underlying retinal pigmented epithelium (RPE), in which the RPE uses lactate (produced by the PRs) as fuel, sparing glucose for the sake of PRs [16,17,18].
In the aging human eye, the RPE and the neighbouring choroid undergo a series of morphological changes, including accumulation of intracellular and extracellular deposits, thickening of the Bruch’s membrane, and reduction of choroidal blood flow [19,20,21,22]. These morphological alterations may hinder oxygen delivery to the RPE and the outer retina generating a tissue environment where oxygen demand may exceed oxygen availability (hypoxia). Thus, hypoxia is a factor in the development of retinal dysfunction and degeneration during aging [23, 24], as well as in other pathologies including diabetic retinopathy and retinal detachment [25, 26].
An evolutionarily conserved pathway mediated by the hypoxia-inducible transcription factors (HIFs) is crucial for the cellular adaptation to hypoxia [27]. HIFs have been associated with the development of age-related macular degeneration (AMD), providing potential therapeutic targets [28]. Central to the HIF pathway is the von Hippel-Lindau tumor suppressor protein (VHL), that forms an E3 ubiquitin ligase complex together with elongins B and C, cullin 2, and RING box protein 1 [29] and ubiquitinates hydroxylated HIF-alpha (HIFA) subunits leading to their proteasomal degradation in normoxia. Inactivation of VHL prevents HIFA degradation activating the heterodimeric transcription factors HIF1 and HIF2, and, thus, enables the activation of the cellular response to hypoxia even in normoxic conditions. We previously generated rod-specific Vhl knockout mice (\(rod^{\varDelta \ Vhl}\)) where chronic HIF1 (but not HIF2) activation led to a late-onset and age-dependent PR and RPE degeneration, mimicking aspects of retinal degenerative diseases associated with aging in humans [30, 31].
In this study, we addressed the relationship between energy metabolism and PR degeneration using chronic HIF activation as a model of age-related retinal degeneration. Using two-photon laser scanning microscopy (TPLSM) of genetically encoded biosensors, layer-specific proteomics and in situ enzymatic assays, we analyzed lactate and glucose dynamics at the cellular level, investigated oxidative phosphorylation pathway (OXPHOS) and the tricarboxylic acid (TCA) cycle, and addressed their contribution to PR survival.
Results
Chronic HIF activation in rods reduces outer segment length and increases lactate production
Age-related changes may impair oxygen flux from the choroid to the RPE and outer retina, generating a hypoxic environment in the retina. Hypoxia activates an intracellular molecular response that includes activation of HIF transcription factors and production of vascular endothelial growth factor (VEGF), the main driver of neovascularization in wet AMD [32]. However, chronic hypoxia not only drives neovascularization but also profoundly affects various cellular metabolic processes that may contribute to the development of dry AMD. To generate a mouse model that mimics the situation of the retina in the aged eye, we inactivated Vhl specifically in rods using the Vhlflox/flox;OpsinCre (\(rod^{\varDelta \ Vhl}\)) mice [30, 31] (Fig. 1a) . Vhl inactivation led to stabilization and chronic activation of HIF transcription factors in rods resulting in a late onset and slowly progressing HIF1-dependent PR degeneration [30, 31] (Fig. 1b).
To identify mechanisms leading to the degeneration in this disease model, we focused on early alterations. Already at 2.5 months of age, a time point before PR degeneration started, rod OS were significantly shorter in \(rod^{\varDelta \ Vhl}\) mice while outer nuclear layer (ONL) thickness and IS length were similar to controls (Fig. 1c). OS shortening progressed with time and was accompanied by reduced IS length and ONL thickness at 6 months of age (Fig. 1d). We hypothesized that chronic HIF activation led to early anabolic deficits affecting OS biosynthesis preceding the onset of PR degeneration. Since anabolism is closely linked to energy metabolism, we focused on glycolysis and mitochondrial metabolism, including OXPHOS and the TCA cycle, to unravel the mechanistic basis of the anabolic deficit (Fig. 1e).
To investigate glycolytic metabolites, we delivered FRET biosensors to PR using subretinal injections of adeno-associated viruses (AAV) and imaged the sensors in acute retinal slices by TPLSM (Fig. 2a,b). Since PRs are light sensitive cells and their metabolism may differ in light and darkness [7], we first investigated whether TPLSM activates phototransduction by monitoring Ca2+ levels during imaging using rod-specific expression of GCaMP6s (Fig. S1a,b). GCaMP6s is a genetically encoded single fluorophore calcium sensor consisting of a circularly permuted green fluorescent protein (cpGFP), the calcium-binding protein calmodulin (CaM), and the CaM-interacting M13 peptide [33]. Ca2+ binding to the CaM-M13 complex induces a conformational change in the protein leading to an increase in the brightness of cpGFP that is detectable by microscopy (Fig S1b).
TPLSM induced a drop in the fluorescence signal intensity in photoreceptor segments (PS) and synaptic terminals in the outer plexiform layer (OPL) of wild-type rods (Fig. S1c,d), indicating a decrease in intracellular Ca2+ concentrations during imaging. The stable GCaMP6s fluorescence signal in Gnat1a-/- mice with their light-insensitive rods [34] (Fig. S1e) verified that the drop observed in wild-type mice was due to the activation of the phototransduction cascade by the laser light and/or by the fluorescence emitted from the sensor. Considering that the emission spectra of all used biosensors overlap, our imaging data for lactate and glucose relate to light-exposed rods.
We used the FRET biosensor Laconic [35] to analyze lactate, the end product of glycolysis. Laconic is a fusion protein composed of the lactate binding bacterial transcription regulator LldR and the FRET pair mTFP and Venus. Binding of lactate decreases FRET efficiency leading to reduced Venus signal (Fig. 2b). To compare basal lactate levels in rods of \(rod^{\varDelta \ Vhl}\) and control mice, we depleted intracellular lactate by using the trans-acceleration property of the monocarboxylate transporters (MCTs) upon application of the non-metabolized MCT substrate oxamate (Fig. 2c) [36]. Due to the symporter type carrier characteristics of the MCTs, occupancy of the transporters in trans by oxamate induces import of oxamate into and extrusion of lactate out of the cell. Successful depletion of intracellular lactate by oxamate was confirmed by assessing lactate levels in aglycemic conditions, which reduced the Laconic FRET signal to a similar level (Fig. 2d). When normalized to this nominal zero lactate, there was no significant difference between rod lactate levels in \(rod^{\varDelta \ Vhl}\) and control mice (Fig. 2e,f). To test whether chronic HIF activation affected lactate production rather than lactate levels, we inhibited MCT mediated lactate transport across the cellular membrane with AR-C155858 [37] (Fig. 2g). This resulted in an increase of intracellular lactate (Fig. 2h) demonstrating that rods are net lactate producers in the presence of glucose as the exclusive fuel. Higher amplitude and higher rate of lactate accumulation suggested that PRs of \(rod^{\varDelta \ Vhl}\) mice are faster lactate producers than those of control mice (Fig. 2h,i).
Increased lactate production may be a consequence of increased glycolysis or decreased mitochondrial pyruvate consumption. We tested the first possibility by analysing intracellular glucose dynamics in retinal cells by TPLSM of the FRET glucose sensor FLIIP (Fig. 3a). FLIIP is composed of the glucose/galactose binding protein of Escherichia coli (MglB) and the FRET pair eCFP and Citrine and exhibits an increase in FRET efficiency upon glucose binding [38]. The sensor was delivered to PRs and to inner retinal cells by subretinal and intravitreal AAV injections, respectively (Fig. 3b). The rate of glycolysis at its entry point was estimated by pharmacological inhibition of glucose import using cytochalasin B, a competitive inhibitor of glucose transporters (Fig. 3c-e) [39,40,41]. Basal glucose levels were assessed by normalizing the FRET traces to nominal zero glucose using medium without glucose (aglycemia) at the end of the experiment (Fig. 3d,e). Surprisingly, baseline glucose levels and glucose consumption rate in PRs did not differ between \(rod^{\varDelta \ Vhl}\) and control mice (Fig. 3f,g), suggesting that chronic HIF activation did not accelerate glycolysis in PRs and that increased intracellular lactate levels may have resulted from reduced mitochondrial pyruvate utilization. The very sharp drop of the FRET signal after blocking glucose import indicated that PRs have a remarkably high glucose consumption rate that was significantly higher than that observed in cells of the inner nuclear layer (INL) and ganglion cell layer (GCL) (Fig. 3d,e,g). Similarly, basal glucose levels were higher in PRs than in cells of the inner retina (Fig. 3d-f). Although FRET sensors may be affected by the cellular environment [42], raw FRET signals from FLIIP in aglycemic conditions were comparable between PRs, and cells of the INL and GCL (Fig. 3h), suggesting a similar sensor behavior in the different retinal cell types.
Together, our data showed that even though glucose dynamics were not altered, lactate production was enhanced in rods of \(rod^{\varDelta \ Vhl}\) mice, suggesting that pyruvate is redirected to lactate upon chronic HIF activation. Since pyruvate is the main mitochondrial substrate and most of the mitochondria within PRs are localized in the IS, we established the retinal layer separation (ReLayS) method [43] and used proteomics (Fig. 4a) to investigate the protein landscapes in PS and soma (ONL) of \(rod^{\varDelta \ Vhl}\) and control mice. Principal component analysis (PCA) of PS and ONL proteomes showed a strong separation between the two groups (Fig. 4b), indicating their markedly different protein compositions. The PS samples were strongly enriched in proteins belonging to the phototransduction cascade, the cilium and mitochondria, while many of the most abundant proteins in the ONL samples were nuclear proteins (Fig. 4c and Table S1) validating the separation procedure. Several proteins involved in glycolysis were significantly upregulated in the ONL (Fig. 4d and Table S2) and PS (Fig. 4e and Table S3) of \(rod^{\varDelta \ Vhl}\) mice. Of particular interest were the increased levels of the glucose transporter GLUT3 (encoded by Slc2a3), since GLUT3 is not present in PRs under physiological conditions [44, 45]. This finding was verified by immunofluorescence studies showing a strong but patchy GLUT3 signal in the ONL and PS of \(rod^{\varDelta \ Vhl}\) but not control mice (Fig. 4f), which resembled the expression of Cre in the ONL of \(rod^{\varDelta \ Vhl}\) mice (Fig. S2). Importantly, we also detected increased levels of the lactate transporter MCT4 (encoded by Slc16a3) in both the ONL (Fig. 4d and Table S5) and PS (Fig. 4e and Table S3) of \(rod^{\varDelta \ Vhl}\) mice. This is of significance since MCT4 is not blocked by AR-C155858 [37], the inhibitor used to determine lactate production by TPLSM (Fig. 2g-i). Thus, our imaging results may have even somewhat underestimated the magnitude of the lactate production in PRs of \(rod^{\varDelta \ Vhl}\) mice, because lactate could still have been exported through the MCT4.
Finally, gene set enrichment analysis (GSEA) yielded an increased enrichment score of proteins belonging to the hallmark hypoxia and a reduced score of proteins involved in OXPHOS in both PS and ONL of \(rod^{\varDelta \ Vhl}\) mice (Fig. 4g,h). This indicated that the chronic activation of HIFs in rods caused a proteomic landscape characteristic for hypoxia and suggested that this intracellular environment profoundly affected OXPHOS.
Chronic HIF activation in rods reduces OXPHOS
Our TPLSM and proteomics data suggested a preferential conversion of pyruvate to lactate pointing to downregulation of OXPHOS in PRs of \(rod^{\varDelta \ Vhl}\) mice. To address OXPHOS in PRs directly, we examined the activity of respiratory complexes IV (cytochrome c oxidase (COX)) and II (succinate dehydrogenase (SDH)) by in situ enzymatic assays. Prominent COX and SDH signals in the IS, the OPL, and the inner plexiform layer (IPL) of control mice indicated strong OXPHOS activity in the retinal regions known to be enriched in mitochondria (Fig. 5a,b). In \(rod^{\varDelta \ Vhl}\) retinas, however, areas of reduced COX and SDH activities were evident in the IS already at 1.5 months of age, several weeks before the onset of degeneration. This phenotype deteriorated with age and depended on HIF1A but not HIF2A (black arrows in Fig. 5a,b) and, thus, on the same HIF isoform that drives retinal degeneration in \(rod^{\varDelta \ Vhl}\) mice [31]. A HIF1 dependent impairment of COX was further evidenced by the absence of COX4i1 (a catalytic subunit of COX) in parts of the IS in \(rod^{\varDelta \ Vhl}\) and \(rod^{\varDelta \ Vhl;Hif2\alpha }\) (white arrows in Fig. 5c) but not \(rod^{\varDelta \ Vhl;Hif1\alpha }\)or control retinas. The patchy pattern of COX4i1 expression was very similar to the pattern of COX and SDH activity (Fig. 5a,b) and was likely due to the mosaic Cre expression in OpsinCre mice [30, 46] (Fig. S2b). Indeed, only Cre-positive areas of the ONL were associated with reduced COX4i1 levels in the adjacent IS (white arrows in Fig. 5d). It is unlikely that this effect was specific for the mitochondria in the rod IS, however, mitochondria localized in the photoreceptor synaptic terminals could not be assessed due to their close proximity to dendritic mitochondria belonging to 2nd order neurons in the retina.
Since reduced activity of the OXPHOS complexes II and IV in rods of \(rod^{\varDelta \ Vhl}\) mice corroborated the GSEA analyses of the layer-specific proteomic data (Fig. 4h), we further examined proteins directly related to OXPHOS, such as mitochondrial complex subunits and assembly factors. The mouse MitoCarta3.0 (http://www.broadinstitute.org/mitocarta) lists 165 such proteins, 81 of which were identified with high confidence in the PS of \(rod^{\varDelta \ Vhl}\) and control mice. 28 of these proteins, belonging to complexes I, III, IV or V, including COX4i1, were significantly downregulated in \(rod^{\varDelta \ Vhl}\) PS (Fig. 5e). The only upregulated OXPHOS-related protein was the ATPase inhibitory factor 1 (ATPIF1), suggesting complex V inhibition in rods of \(rod^{\varDelta \ Vhl}\) mice. Interestingly, BNIP3 was detected in three of the six PS samples from \(rod^{\varDelta \ Vhl}\) but in none of the control mice (Table S5). Although its presence in \(rod^{\varDelta \ Vhl}\) was inconsistent, BNIP3 may be an interesting protein for further investigations due to its connection to cell death [47] and mitochondrial energetics [48]. Such future experiments will address the potential of BNIP3 to reduce rod survival in conditions of activated HIF proteins. Collectively, these data indicated that chronic HIF1 activation in rods caused OXPHOS dysregulation preceding the onset of cell degeneration.
Chronic HIF activation in rods affects mitochondrial protein composition
Downregulation of OXPHOS-related proteins in PRs of \(rod^{\varDelta \ Vhl}\) mice was unlikely due to loss of mitochondria. Only 14% of the mitochondrial proteins identified in the PS of \(rod^{\varDelta \ Vhl}\) mice were significantly downregulated (Fig. 6a and Table S4) and well-known mitochondrial marker proteins, such as mitochondrial apoptosis-inducing factor 1 (AIFM1), mitochondrial import receptor subunit TOM70 (TOM70) and dihydrolipoamide dehydrogenase (DLD) had comparable levels in PS of \(rod^{\varDelta \ Vhl}\) and control mice (Fig. 6b). Immunofluorescence labeling for voltage dependent anion channel 1 (VDAC1) also did not indicate loss of mitochondria in rods lacking COX4i1 in \(rod^{\varDelta \ Vhl}\) retinas (white arrows point to areas with reduced COX4i1 staining in Fig. 6c). Vdac1 expression was even significantly upregulated in the PS of \(rod^{\varDelta \ Vhl}\) mice (Fig. S3a). To address a possible increase in transcriptional regulation of mitochondrial proteins, we assessed the relative expression levels of translocase of outer mitochondrial membrane 20 (Tomm20), another outer mitochondrial membrane protein. However, Tomm20 expression was significantly downregulated in ONL and PS samples from \(rod^{\varDelta \ Vhl}\) mice (Fig. S3a) suggesting that there was no general increase in transcription of mitochondrial genes. Together, these data indicate that chronic HIF activation in rods inhibited OXPHOS without affecting the overall mitochondrial biomass.
Although only one (cytochrome c oxidase subunit 2 (MT-CO2)) of the five mtDNA-encoded proteins identified in our proteomics data was downregulated in rods of \(rod^{\varDelta \ Vhl}\) mice (Fig. 5e, Fig. 6d and Table S5), we examined whether OXPHOS downregulation was a result of compromised mtDNA integrity. Relative quantification of mtDNA copy number showed no difference between PRs of \(rod^{\varDelta \ Vhl}\) and control mice (Fig. 6e), indicating that mtDNA integrity was not strongly affected by chronic HIF activation at 2.5 months of age.
Nevertheless, immunofluorescence labeling for the nuclear-encoded mitochondrial transcription factor A (TFAM) showed decreased immunoreactivity in patchy areas in the ONL of \(rod^{\varDelta \ Vhl}\) mice as early as at 2 months of age, a phenotype that became even more pronounced with increasing age (Fig. 6f). Consistently, proteomics data showed a significant downregulation of TFAM and three other proteins involved in mtDNA maintenance (EXOG, PolG-\(\upbeta\), ATAD3) in PRs of 2.5 months old \(rod^{\varDelta \ Vhl}\) mice (Fig. 6g). In addition, two proteins involved in mtRNA metabolism (LACTB2 and REXO2) and three in mitochondrial translation (MRP-L4, MRP-S2, and TrpRS) were differentially and mostly downregulated in the ONL and PS of \(rod^{\varDelta \ Vhl}\) mice (Fig. 6g). Similarly, expression of two mtDNA-encoded genes (mt-Co2 and NADH-ubiquinone oxidoreductase chain 1 (mt-Nd1)) was significantly downregulated in ONL and PS samples of \(rod^{\varDelta \ Vhl}\) mice (Fig. S3b). Together, our data suggest that chronic HIF1 activation in rods led to an early severe downregulation of OXPHOS-related proteins without a significant loss of mitochondria despite an affected mtDNA maintenance machinery. Furthermore, rods seemed to survive despite OXPHOS dysfunction, at least for some time.
Inhibition of OXPHOS in rods does not lead to early anabolic deficiency
To further investigate the dependency of rod OS biogenesis and survival on OXPHOS, we generated \(rod^{\varDelta \ Cox10}\) mice that have impaired OXPHOS activity due to compromised complex IV assembly [49] in rods (Fig. 7a). Similar to chronic HIF activation, Cox10 deletion (indicated by the presence of Cre) led to the loss of the complex IV subunit COX4i1 (Fig. 7b) but not the outer mitochondrial membrane protein VDAC1 (Fig. 7c; white arrows point to areas with reduced COX4i1 staining) in rod IS. Furthermore, rods showed remarkable resilience against impaired OXPHOS activity and survived for several months (Fig. 7b,d-f) as indicated by the presence of Cre-positive rods in 6-month-old \(rod^{\varDelta \ Cox10}\) mice (Fig. 7b). At this time point, ONL thickness was only moderately affected (Fig. 7d,f) and scotopic and photopic electroretinograms were in the normal range (Fig. S4). Significant loss of rods and reduced scotopic function were observed only at 12 months of age (Fig. 7d,g and Fig. S4). In strong contrast to \(rod^{\varDelta \ Vhl}\) mice, OS length in \(rod^{\varDelta \ Cox10}\) mice was not affected at the early age (Fig. 7d,e), indicating that OS biogenesis did not depend on OXPHOS activity.
Chronic HIF activation in rods downregulates the citric acid cycle independently of OXPHOS
Since enzymatic staining showed a downregulation of SDH in \(rod^{\varDelta \ Vhl}\) mice, an enzyme complex serving not only OXPHOS but also the TCA cycle, we hypothesised that chronic HIF activation affects the TCA cycle and, thus, generation of TCA metabolites important for anabolism and OS biogenesis. Detailed analysis of the proteomics data revealed a significant downregulation not only of SDH subunit A (SDHA) but also of six additional TCA associated enzymes in \(rod^{\varDelta \ Vhl}\) mice (Fig. 8a-c). Immunofluorescence labeling for SDHA, citrate synthase (CS), and aconitase 2 (ACO2) validated the proteomics data and indicated that the TCA cycle was severely affected in \(rod^{\varDelta \ Vhl}\) but not \(rod^{\varDelta \ Cox10}\) mice (Fig. 8d-f). Together, these results suggest that chronic HIF activation concomitantly affected the TCA cycle and OXPHOS in rods, leading to anabolic dysfunction and ultimately to PR degeneration (Fig. 8g).
Discussion
We investigated PR energy metabolism in a mouse model of slow progressing retinal degeneration [30, 31]. Our results show that PRs have a remarkably high glucose consumption rate, and that chronic HIF activation induces an increase in lactate production while severely affecting OXPHOS and the TCA cycle in rods. We also show that OXPHOS deficiency on its own is surprisingly well tolerated by rods, while additional impairment of the TCA cycle leads to an early dysregulation of rod OS biogenesis and retinal degeneration.
Since metabolism is fundamental factor for cell survival, tissue integrity and disease development, the investigation of glucose, the main metabolic precursor, becomes central for the understanding of pathophysiological processes. Using TPLSM and FRET biosensors, we provide, to our knowledge, the first data on lactate (Fig. 2) and glucose (Fig. 3) dynamics at cellular resolution in the intact retinal tissue. Results indicate that PRs are not only net lactate producers in the presence of glucose as the only fuel (Fig. 2h) but that they also have markedly higher basal glucose levels and glucose consumption rate than cells of the inner retina (Fig. 3). These results confirm the long-standing hypothesis of PRs being the main site of aerobic glycolysis in the retina [4, 12, 16, 18], a phenomenon that has been observed on tissue level already a century ago by Otto Warburg [9].
It is of interest that rods express GLUT3 upon chronic HIF activation (Fig. 4f), a glucose transporter with higher glucose affinity and greater transport capacity than GLUT1 [50]. Since GLUT3 is not present in rods under normoxic conditions [44, 45, 51], it is conceivable that rods may attempt to increase glucose uptake as a response to reduced OXPHOS activity during HIF activation. However, there was no measurable increase in basal glucose levels or glucose consumption rate in rods of \(rod^{\varDelta \ Vhl}\) mice (Fig. 3). The absence of glycolysis stimulation in the face of OXPHOS failure was unexpected. One possibility is that in these cells the glycolytic rate is already near a maximum under normoxic conditions. Alternatively, chronic hypoxia may redirect anabolic branches like the pentose phosphate pathway towards the main glycolytic pathway for the production of ATP and lactate release, as observed in stressed neurons [52]. This redeployment of glucose would explain the combination of higher lactate production with conserved glucose consumption.
The increase of lactate production in rods upon chronic HIF activation (Fig. 2i) is in accordance to previous studies reporting oxygen level dependent lactate production in isolated retinas [53]. Notably, our data rather underestimated the magnitude of this increase for two reasons: First, the mosaic Cre expression (Fig. S2b) inevitably led to the inclusion of both Vhl knockout and wild type cells in our TPLSM measurements. Second, AR-C155858 potently inhibits MCT1 and MCT2 but not MCT4 [37] that is even significantly upregulated in the ONL and PS of \(rod^{\varDelta \ Vhl}\) mice (Fig. 4d,e) and, thus, lactate could be exported from the cell through the high affinity MCT4 [54], escaping detection by the biosensor.
Since the glycolytic flux through the hexokinases was not altered, increased lactate production in rods of \(rod^{\varDelta \ Vhl}\) mice was likely the result of a reduced pyruvate utilization by the mitochondria. In addition, feedback mechanisms and allosteric regulation of TCA cycle enzymes [55, 56] may also have affected function of mitochondrial metabolic pathways and thus contributed to the effect. Proteomics data from isolated PS and ONL (Fig. 4h, Fig. 5e, Fig. 8a,b) as well as results from enzymatic in situ experiments (Fig. 5a,b) and immunofluorescence studies (Fig. 5e, Fig. 8d-f) showed that chronic HIF1 (but not HIF2) activation in rods severely downregulated protein levels of several key enzymes involved not only in OXPHOS but also in the TCA cycle (Fig. 5, Fig. 8). Despite the downregulation of these proteins and the severely reduced activity of SDH and COX in photoreceptor segments, mitochondria remained present (Fig. 6c) and photoreceptors survived for several weeks (Fig. 1b,c). Survival of rods with OXPHOS deficiency was further corroborated in \(rod^{\varDelta \ Cox10}\) mice, suggesting that the high glucose consumption rate of rods may generate sufficient energy for the cells to function and survive in such conditions. This extends the list of cells (including astrocytes [57], oligodendrocytes [58], brain neurons [59, 60], and skeletal muscle cells [49]) that are reported to survive as glycolytic cells in the absence of Cox10. Additionally, it has been suggested that cells in the retina can also use fumarate as electron acceptor instead of oxygen by reversing the succinate dehydrogenase reaction [61]. The resulting succinate can then be transported to the neighbouring RPE, oxidised to malate, and shuttled back to the retina. The existence of such malate-succinate shuttle could also help to explain the resilience of COX deficient rod PRs. Although PRs eventually die in both \(rod^{\varDelta \ Vhl}\) (Fig. 1b-d) and \(rod^{\varDelta \ Cox10}\) mice (Fig. 7d-g), the degenerative phenotype was more severe and appeared earlier in \(rod^{\varDelta \ Vhl}\) mice. For survival, rods may thus not so much depend on OXPHOS but more heavily so on functional TCA cycle that is seemingly unaffected in \(rod^{\varDelta \ Cox10}\) but not \(rod^{\varDelta \ Vhl}\) mice. This underlines the significance of the TCA cycle with its pivotal role in anabolic and homeostatic processes that are essential, among others, for OS biogenesis. Reduced TCA cycle activity as well as shorter rod OS length have been reported in mouse models of retinitis pigmentosa (RP), an inherited retinal degenerative disorder with high genetic heterogeneity [62,63,64,65]. Shorter rod OS length was also described in RP patients [66, 67] as well as in the context of hypoxia in mice [68] and human [69]. Thus, TCA cycle dysregulation might represent a common molecular feature preceding PR degeneration. As suggested by others [64], supporting the TCA cycle may thus constitute a promising strategy to improve PR survival in a variety of patients.
Conclusions
Our study demonstrates the consequences of an activated hypoxic response for glucose flux and lactate levels in rods and highlights the importance of the TCA cycle for rod survival. Despite the need for high ATP levels for maintenance of the visual function [7], rods survive in the absence of OXPHOS, indicating that they can satisfy their energy needs by alternative pathways such as aerobic glycolysis and/or beta-oxidation of fatty acids. In fact, it has been proposed that only about 20% of glucose is completely oxidized in the outer retina while the rest is used for lactate production [4] and that the long-chain fatty acid palmitate may serve as a fuel for mitochondrial respiration in the isolated retina [70, 71]. However, no sufficiently effective alternative pathway to the TCA cycle may exist to provide satisfactory amounts of intermediate metabolites needed for the constant biosynthesis of complex cellular structures such as the outer segments of PRs.
Methods
Animals
All experimental procedures were approved by the local veterinary authorities, conforming to the guidelines of the Swiss Animal Protection Law, Veterinary Office, Canton Zurich (Act of Animal Protection 16 December 2005 and Animal Protection Ordinance 23 April 2008) and were performed in accordance with the respective national, federal and institutional regulations. Mice were maintained as breeding colonies at the Laboratory Animal Services Center of the University of Zurich in a 14h:10h light-dark cycle. Mice had access to food and water ad libitum. The average light intensity at cage levels was 60-150 lux, depending on the position in the rack.
Vhlflox/flox [72], Hif1\(\alpha\) flox/flox [73], Hif2\(\alpha\) flox/flox [74], Cox10flox/flox [49], and OpsinCre (LMOPC1 [46]) mice were intercrossed to obtain \(rod^{\varDelta \ Vhl}\) (Vhlflox/flox;OpsinCre), \(rod^{\varDelta \ Vhl;Hif1\alpha }\)(Vhlflox/flox;Hif1\(\alpha\) flox/flox;OpsinCre), \(rod^{\varDelta \ Vhl;Hif2\alpha }\) (Vhlflox/flox;Hif2\(\alpha\) flox/flox; OpsinCre), and \(rod^{\varDelta \ Cox10}\) (Cox10flox/flox;OpsinCre) mice. All breeding pairs were heterozygous for OpsinCre and littermates without OpsinCre served as controls. Rod-specific Cre expression in OpsinCre mice starts around postnatal day 7 and increases up to 6 weeks of age [46]. Homozygous rod transducin \(\alpha\)-subunit knockout mice [34] (Gnat1a-/-) were used as controls in TPLSM calcium imaging. All mice were homozygous for the Rpe65450Leu variant [75, 76]. For tissue sample preparations, mice were euthanized with CO2 followed by decapitation.
Retinal morphology analysis
Eyes were enucleated, fixed in glutaraldehyde (2.5% in cacodylate buffer) for 12-24 h at 4\(^{\circ }\)C, trimmed, post-fixed in 1% osmium tetroxide, and embedded in Epon 812 as described [77]. Retinal cross-sections of 0.5 \(\upmu\)m were cut through the optic nerve head, stained with toluidine blue, and analyzed by light microscopy (Zeiss, Axioplan). OS and IS length, and the thickness of the ONL were measured at indicated distances from the optic nerve head using the Adobe Photoshop CS6 ruler tool (Adobe). Mann-Whitney nonparametric test was used to compare the overall OS and IS length, as well as ONL thickness of \(rod^{\varDelta \ Vhl}\) and \(rod^{\varDelta \ Cox10}\) mice to their respective controls. P-values < 0.05 were considered to show significant differences.
Cloning and virus production
The GCaMP6s gene [78] was cloned into an AAV plasmid backbone containing a rod PR-specific mouse opsin promoter (mOP) (a gift from Sanford L. Boye) [79]. The laconic AAV plasmid was a gift from Luis Felipe Barros (Addgene plasmid #44238) [35]. The FLIIP AAV plasmid [38] was codon-diversified to avoid recombination during AAV production [80, 81]. All AAVs were produced by the Viral Vector Facility of the Neuroscience Center Zurich and packaged in the AAV2(QuadYF+TV; 7m8)/2 capsid [82].
Intraocular injections of AAVs carrying fluorescent sensors
For intraocular AAV injections, the pupils of mice were dilated with Cyclogyl 1% (Alcon Pharmaceuticals) and Neosynephrine 5% (Ursapharm Schweiz). Mice were anesthetized by a subcutaneous injection of ketamine (85 mg/kg, Parke-Davis) and xylazine (10 mg/kg, Bayer AG). Viscotears (Bausch & Lomb Swiss AG) were applied to keep the eyes moist. 1 x 109 total viral genome in 1 \(\upmu\)l total volume were injected either intravitreally or subretinally to transduce the inner or the outer retinal cells, respectively. Anesthesia was reversed by a subcutaneous injection of antisedan (2 mg/kg, Atipazemole). Two to three weeks after injection, mice were subjected to fluorescent funduscopy and OCT to test for transgene expression and detect potential injection-inflicted tissue damage including bleeding or persistent retinal detachment. Such eyes were excluded from the study. \(rod^{\varDelta \ Vhl}\) and ctrl mice were injected at 5-6 weeks of age and TPLSM experiments was performed at 10-14 weeks of age.
Acute retinal preparations for TPLSM ex vivo imaging
Mice were housed in a facility with a reversed 12/12 h light/dark cycle for at least 5 days before the experiments. Prior to dissection, mice were dark adapted for at least 2 h, and all steps of the dissection procedure were carried out in dim red light. Retinas were rapidly dissected through a slit in the cornea and placed in freshly prepared artificial cerebrospinal fluid (ACSF) containing: 126 mM NaCl, 3 mM KCl, 2 mM CaCl2, 1.25 mM NaH2PO4, 26 mM NaHCO3, 2 mM MgSO4, 6 mM D-glucose), bubbled with oxycarbon (95% O2/ 5% CO2, PanGas, 820, and kept at RT. PRs were imaged in the cross-section of acute retinal slices as previously described [83]. In brief, half retinas were flatmounted with the PRs side facing up on a nitrocellulose filter membrane and cut in 250 \(\upmu\)m-thick slices using a razor blade attached to a custom-made tissue chopper. For imaging inner retinal cells, flatmounted retinas with the ganglion cells facing up were used without slicing the tissue. The tissue was left to rest while slowly warming the ACSF to 35\(^{\circ }\)C for 30 min before imaging.
Two-photon excitation laser scanning microscopy
A custom-made two-photon laser scanning microscope equipped with a two-photon laser with < 120 fs temporal pulse width (InSight DeepSee Dual; Spectra-Physics) and a 16x water immersion objective (Nikon N16XLWD-PF, 0.8 NA, 3 mm WD) was used for image acquisition [84]. Excitation and emission beam paths were separated by a dichroic mirror (F73-825; AHF Analysentechnik). Dichroic mirrors at 506nm (F38-506; AHF Analysentechnik) and 560 nm (F38-560; AHF Analysentechnik) further separated the emission light into colored components, that were then focused on GaAsP photomultipliers (H10770PA-40sel; Hamamatsu Photonics) equipped with filters for cyan (475/50; AHF Analysentechnik), and green/yellow (542/50; AHF Analysentechnik) wavelenghts. ScanImage 3.8 [85] and custom-written LabVIEW software (Version 2012; National Instruments) were used for image control and data acquisition.
For imaging, constantly oxygenated ACSF was gravity-fed into the recording chamber at a flow rate of \(\approx\) 2 ml/min and temperature was maintained at 35\(^{\circ }\)C using a temperature controller (V TC05, Luigs & Neumann GmbH) with an in-bath temperature sensor. All solutions used during imaging were infused using the same gravity flow system and osmolarity was adjusted at \(\approx\) 300 mOsm. Data was acquired from acute retinal slices at 50-70 \(\upmu\)m below the cross-section surface, from cells in the GCL at 0-10 \(\upmu\)m and from cells in the INL at 50-60 \(\upmu\)m below the surface.
Calcium imaging parameters and analysis
GCaMP6s was excited at 940 nm wavelength and images were acquired at 11.84 Hz with a 128x256-pixel resolution. Analyses were performed with ImageJ (National Institutes of Health) and regions of interest were selected manually. Prism software (GraphPad) was used for data visualization.
Lactate and glucose imaging parameters and analysis
Laconic and FLIIP were excited at 870 nm wavelength and images were acquired at 5.94 Hz with a 256x256-pixel resolution. ACSF containing 20 mM oxamate was used in the trans-acceleration imaging protocol. ACSF containing 5 \(\upmu\)M AR-C155858 (BioTechne Tocris) was used in the lactate transport inhibition imaging protocol. ACSF containing 20 \(\upmu\)M cytochalasin B (BioTechne Tocris, 5474) was used in the glucose transport inhibition imaging protocol. ACSF containing 0 mM glucose was used in the aglycemia imaging protocol.
Whole-frame image analysis was carried out using custom-made code for MATLAB 2017b (MathWorks) and ImageJ (National Institutes of Health). For each experiment, images were aligned using a 2D convolution engine to account for xy drift in time. To optimize signal to noise ratio, time smoothing of 11 frames for Laconic and 5 frames for FLIIP was applied, as well as automatic thresholding using the Li’s minimum cross entropy method [86]. Ratiometric data were visualized with Prism software (GraphPad). Mann-Whitney nonparametric test was used to compare \(rod^{\varDelta \ Vhl}\) and ctrl mice. Kuskal-Wallis nonparametric test and Dunnett’s multiple comparison test were used to compare PRs, GCL and INL cells. Because AAVs were injected subretinally and biosensor expression was CMV driven, both rod and cone PRs may have been transduced. However, given that rods outnumber cone PRs by a factor of 30 in the mouse retina, all imaging data represent primarily rods.
Label-free quantitative proteome analysis
Sample preparation with the ReLayS method
Samples were prepared using the ReLayS method [43]. In brief, retinas from 2.5-months-old \(rod^{\varDelta \ Vhl}\) (N = 6) and control (N = 6) mice were rapidly dissected through a slit in the cornea and placed in PBS on ice. Vitreous was carefully removed using a pair of forceps and retinas were halved through the optic nerve head. The half-retinas were flattened on a small piece of nitrocellulose filter membrane with the photoreceptor side facing up and frozen on a metal platform placed on dry ice. Frozen retinal flatmounts were stored at -80\(^{\circ }\)C until further use. PS and ONL samples were consecutively separated from frozen half-retinal flatmounts using the adherence of the cellular structures to a nitrocellulose filter membrane placed on top of the photoreceptors or the ONL, respectively. The PS and ONL membranes were frozen on a metal platform placed on dry ice and transferred to 1.5 ml Protein LoBind\(^{\circledR }\) Tubes (Eppendorf; 0030108116). For protein isolation, samples were homogenized by sonication with an ultrasonic homogenizer in ice-cold Tris-HCl (100 mM, pH 8.0) containing protease inhibitors (Sigma-Aldrich, P2714). After sonication, 10% SDS (final concentration 1% SDS in Tris-HCl, pH 8.0) was added and samples were incubated at 75\(^{\circ }\)C for 10 min. Protein samples were stored at -20\(^{\circ }\)C until further use.
Mass spectrometry
The mass spectrometry proteomic data have been deposited to the ProteomeXchange Consortium via the PRIDE [87] partner repository with the dataset identifier PXD034057.
From each sample, total proteins were proteolyzed with LysC (Wako Chemicals, Neuss, Germany) and trypsin (Promega) using a suspension trapping protocol (S-Trap, Protifi) to remove SDS according to the manufacturer’s instructions. Briefly, samples were reduced and carbamidomethylated, followed by acidification with phosphoric acid and addition of methanol to a final concentration of >70% before loading to the trap columns. Proteins were washed while trapped on column, then digested on column with LysC (2 hours at RT) followed by trypsin (overnight, at 37\(^{\circ }\)C). Peptides were collected by centrifugation and acidified. Eluted peptides were analyzed on a Q Exactive HF mass spectrometer (Thermo Fisher Scientific) in the data dependent mode. Approximately 0.5 \(\upmu\)g peptides per sample were automatically loaded to the online coupled ultra-high-performance liquid chromatography (UHPLC) system (Ultimate 3000, Thermo Fisher Scientific). A nano trap column was used (300-\(\upmu\)m ID X 5mm, packed with Acclaim PepMap100 C18, 5 \(\upmu\)m, 100 Å; LC Packings) before separation by reversed phase chromatography (Acquity UHPLC M-Class HSS T3 Column 75 \(\upmu\)m ID X 250 mm, 1.8 \(\upmu\)m; Waters) at 40\(^{\circ }\)C. Peptides were eluted from the column at 250 nL/min using increasing ACN concentrations (in 0.1% formic acid) from 3% to 41% over a linear 95-min gradient. MS spectra were recorded at a resolution of 60 000 with an AGC target of 3 \(\times\) 106 and a maximum injection time of 50 ms from 300 to 1’500 m/z. From the MS scan, the 10 most abundant peptide ions were selected for fragmentation via HCD with a normalized collision energy of 28, an isolation window of 1.6 m/z, and a dynamic exclusion of 30 s. MS/MS spectra were recorded at a resolution of 15’000 with an AGC target of 105 and a maximum injection time of 50 ms. Unassigned charges and charges of +1 and above +8 were excluded from precursor selection.
Data processing - protein identification and label-free quantification
Proteome Discoverer 2.4 software (Thermo Fisher Scientific; version 2.4.1.15) was used for peptide and protein identification via a database search (Sequest HT search engine) against SwissProt mouse database (release 2020_02, 17061 sequences), considering full tryptic specificity, allowing for up to two missed tryptic cleavage sites, precursor mass tolerance 10 ppm, and fragment mass tolerance 0.02 Da. Carbamidomethylation of Cys was set as a static modification. Dynamic modifications included deamidation of Asn and Gln, oxidation of Met; and a combination of Met loss with acetylation on protein N-terminus. Percolator [88] was used for validating peptide spectrum matches and peptides, accepting only the top-scoring hit for each spectrum, and satisfying the cut-off values for FDR < 1%, and posterior error probability < 0.01. The final list of proteins complied with the strict parsimony principle.
The quantification of proteins, after precursor recalibration, was based on abundance values (intensity) for unique peptides. Abundance values were normalized to the total peptide amount to account for sample load errors. The protein abundances were calculated by summing the abundance values for admissible peptides. The final protein ratio was calculated using median abundance values of six replicate analyses each. The statistical significance of the ratio change was ascertained employing the approach described in [89], which is based on the presumption that we look for expression changes for proteins that are just a few in comparison to the number of total proteins being quantified. The quantification variability of the non-changing “background” proteins can be used to infer which proteins change their expression in a statistically significant manner.
Data analysis and visualization
PCA was performed on normalized abundance protein levels using ClustVis [90] (http://biit.cs.ut.ee/clustvis/). Heatmap of the top 25 proteins in ONL and PS was generated from the scaled normalized abundance values using the ComplexHeatmap package + [91] and the viridis color palette with the R package circlize [92]. The volcano plots were generated using VolcaNoseR [93] (https://huygens.science.uva.nl/VolcaNoseR2/). GSEA was performed using GSEA (Version 4.1.0) software provided by Broad Institute of Massachusetts Institute of Technology and Harvard University [94, 95]. Analysis was conducted on a pre-ranked protein list based on their relative abundance in \(rod^{\varDelta \ Vhl}\) against control samples for the PS and ONL, respectively. The hallmark gene set from Molecular Signatures Database [96, 97] (MSigDB, Version 7.2) was used for comparison. Normalized protein abundances and abundance ratios were visualized with Prism software (GraphPad).
mtDNA copy number assay
PS and ONL samples were prepared from 2.5-months-old \(rod^{\varDelta \ Vhl}\) (N = 6) and control (N = 6) mice using the ReLayS method [43] as described above. Tissue was lysed by proteinase K (Sigma-Aldrich, 03115879001) treatment at 56\(^{\circ }\)C for 45 min with shacking at 300 rpm and RNase A (100 mg/ml; Thermo Fisher Scientific, 12091021) was added to the samples. DNA isolation was performed with QIAamp DNA Blood Mini Kit (Qiagen, 51104) according to manufacturer instructions. DNA from PS and ONL samples was pooled together to obtain DNA samples from photoreceptors. Relative quantification of mtDNA levels was determined by the ratio of the mitochondrial ND1 (mt-Nd1) gene to the nuclear-encoded 18S rRNA gene using real-time PCR. 16 ng DNA was used as template and the genes of interest amplified using the PowerUp SYBR Green Master Mix (ThermoFisher Scientific) in the ABI QuantStudio 3 system (ThermoFisher Scientific) and specific primer pairs: ND1 fwd 5’-3’: CTAGCAGAAACAAACCGGGC and ND1 rev 5’-3’: CCGGCTGCGTATTCTACGTT; 18S rRNA fwd 5’-3’: CGCGGTTCTATTTTGTTGGT and 18S rRNA rev 5’-3’: AGTCGGCATCGTTTATGGTC. Data analysis was carried out using the 2\(^{\wedge }\)-ddCt method [98].
Immunohistochemistry
Eyes were marked at the nasal limbus, enucleated, and fixed in 4% PFA in PBS for a total of 4h at 4\(^{\circ }\)C. For cryosectioning, eyes were placed in 30% sucrose in PBS for 2h at 4\(^{\circ }\)C and embedded with tissue freezing medium (Leica Biosystems, 81-0771-00). 12-\(\upmu\)m-thick sections were prepared using a Leica cryostat (Biosystems Switzerland AG, CM1860). For paraffin sections, PFA fixed eyes were dehydrated in a series of ethanol solutions of increasing concentrations up to 100% and immersed in xylene prior to paraffin embedding. 5-\(\upmu\)m-thick sections were prepared using a Zeiss microtome (Microm HM 440E). Paraffin sections were deparaffinized, rehydrated, and epitopes were unmasked by heat-induced antigen retrieval at pH 6. For immunofluorescence staining, paraffin- and cryosections were blocked with PBS containing 3% normal goat serum and 1% Triton X-100. The following primary antibodies were incubated overnight at 4\(^{\circ }\)C: rabbit anti-GLUT3 (1:250; Sigma-Aldrich, ab1344), mouse anti-COX4i1 (1:100; Abcam, ab14744), rabbit anti-CRE (1:300; Merck Millipore, 69050-3), rabbit anti-VDAC1 (1:250; Abcam, ab15895), rabbit anti-TFAM (1:250; Abcam, ab131607), rabbit anti-CS (1:250; GeneTex, GTX110624), rabbit anti-ACO2 (1:250; Cell Signaling Technology, 6571), rabbit anti-SDHA (1:250; Cell Signaling Technology, 11998) or rabbit anti-cARR (1:1000; Merck Millipore Chemicals, AB15282). Secondary antibodies (Alexa Fluor 488- or 568-conjugated, 1:500) were applied for 1 h at RT. Fluorescence images were acquired with a Zeiss microscope (Axioplan 2, Zeiss), processed, and analyzed with ImageJ (National Institutes of Health). Pixel intensity profiles through the PS were obtained using the profile plots tool controlled by the line tool in ImageJ.
COX and SDH enzymatic histochemistry
Eyes were marked at the nasal limbus, enucleated, briefly washed in ice cold PBS, and immediately embedded in tissue freezing medium (Leica Biosystems, 81-0771-00). 12-\(\upmu\)m-thick sections were prepared using a Leica cryostat. For COX enzymatic histochemistry, sections were incubated with COX detection solution (1.6 mM DAB, 1.25 mg/ml cytochrom C, 10 \(\upmu\)l catalase in PBS, pH 7.0) for 45 min at 37\(^{\circ }\)C. For SDH enzymatic histochemistry, sections were incubated with SDH detection solution (1.55 mM nitrotetrazolium blue chloride, 0.13 mM sodium succinate, 0.2 mM phenazine methosulfate, 0.1 mM sodium azide in PBS, pH 7.0) for 25 minutes at 37\(^{\circ }\)C. Sections were dehydrated in a series of ethanol solutions of increasing concentrations up to 100% and immersed in xylene. Coverslips were mounted with Cytoseal (Thermo Scientific, 8312-4). Bright-field microscopy images were acquired with a Zeiss microscope (Axioplan 2), processed and analyzed with ImageJ (National Institutes of Health). Pixel intensity profiles through the PS were obtained using the profile plots tool controlled by the line tool in ImageJ.
Scotopic and photopic electroretinography (ERG)
Mice were processed for ERG as described [99]. Briefly, after overnight dark-adaptation and pupil dilation with 1% Cyclogyl (Alcon) and 5% Neosynephrin-POS, mice were anesthetized with ketamin (85 mg/kg) and xylazine (10 mg/kg). Electroretinograms were recorded simultaneously from both eyes using a Diagnosys Celeris rodent ERG device (Diagnosys). Ten flash intensities ranging from 8x10-6 cd*s/m2 to 3 cd*s/m2 and six flash intensities ranging from 1 cd*s/m2 to 200 cd*s/m2 were used for dark- (scotopic) and light-adapted (photopic) single-flash intensity ERG series, respectively. Five sweeps per intensity were averaged for the scotopic and ten sweeps per intensity for the photopic ERGs. The standard rod-suppressive background light (30 cd/m2) was used prior (5 min) and during recordings in photopic conditions. Statistical analysis was performed using two-way ANOVA and Bonferroni’s multiple comparison test (spidergrams) or nested t-test (bar graphs) with Prism software (GraphPad).
Semi-quantitative real-time PCR
PS and ONL samples were prepared from 2.5-months-old \(rod^{\varDelta \ Vhl}\) (N = 6) and control (N = 6) mice using the ReLayS method [43] as described above. Total RNA was isolated from the tissue on the membrane with an RNA isolation kit (Thermo Fisher PicoPure RNA Isolation Kit, KIT0204) including an on-column DNaseI treatment. cDNA synthesis was carried out with oligo-(d)T primers and M-MLV reverse transcriptase (Promega). For semi-quantitative real-time PCR, 10 ng cDNA was amplified using the PowerUp SYBR Green Master Mix (Thermo Fisher Scientific) and specific primer pairs (Actb fwd 5’-3’: CAACGGCTCCGGCATGTGC and rev 5’-3’: CTCTTGCTCTGGGCCTCG; Vdac1 fwd 5’-3’: CAAGGTCACACTGAACATGG and rev 5’-3’: TCACTTTGGTGGTTTCCGT; Tomm20 fwd 5’-3’: TGCATCTACTTCGACCGCAAA and rev 5’-3’: GTCCACACCCTTCTCGTAGTC; mt-Co2 fwd 5’-3’: CCTCCACTCATGAGCAGTCC and rev 5’-3’: AATAACCCTGGTCGGTTTG; mt-Nd1 fwd 5’-3’: CTAGCAGAAACAAACCGGGC and rev 5’-3’: CCGGCTGCGTATTCTACGTT) in the ABI QuantStudio 3 system (Thermo Fisher Scientific). Actin-beta (Actb) was used as a reference housekeeping gene. Data analysis was carried out using the 2\(^{\wedge }\)-ddCt method [98]. Data were visualized with Prism software (GraphPad) and Mann-Whitney nonparametric test was used to compare \(rod^{\varDelta \ Vhl}\) and ctrl mice.
Availability of data and materials
The mass spectrometry proteomic data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD034057.
Abbreviations
- AAV:
-
Adeno-associated viruses
- ACSF:
-
Artificial cerebrospinal fluid
- AMD:
-
Age-related macular degeneration
- COX:
-
Cytochrome c oxidase
- ERG:
-
Electroretinography
- GCL:
-
Ganglion cell layer
- GSEA:
-
Gene set enrichment analysis
- HIF:
-
Hypoxia-inducible factor
- HIFA:
-
HIF-alpha
- INL:
-
Inner nuclear layer
- IPL:
-
Inner plexiform layer
- IS:
-
Photoreceptor inner segments
- MCT:
-
Monocarboxylate transporters
- mOP:
-
Mouse opsin promoter
- MT-CO2:
-
Cytochrome c oxidase subunit 2
- mt-Nd1:
-
NADH-ubiquinone oxidoreductase chain 1
- ONL:
-
Outer nuclear layer
- OPL:
-
Outer plexiform layer
- OS:
-
Photoreceptor outer segments
- OXPHOS:
-
Oxidative phosphorylation pathway
- PCA:
-
Principal component analysis
- PR:
-
Photoreceptor
- PS:
-
Photoreceptor segments
- RPE:
-
Retinal pigmented epithelium
- SDH:
-
Succinate dehydrogenase
- TCA:
-
Tricarboxylic acid
- TFAM:
-
Mitochondrial transcription factor A
- Tomm20:
-
Translocase of outer mitochondrial membrane 20
- TPLSM:
-
Two-photon laser scanning microscopy
- VDAC1:
-
Voltage dependent anion channel 1
- VEGF:
-
Vascular endothelial growth factor
- VHL:
-
Von Hippel-Lindau tumor suppressor protein
References
Wong-Riley M. Energy metabolism of the visual system. Eye Brain. 2010;2:99–116. https://doi.org/10.2147/EB.S9078.
Medrano CJ, Fox DA. Oxygen consumption in the rat outer and inner retina: light- and pharmacologically-induced inhibition. Exp Eye Res. 1995;61(3):273–84. https://doi.org/10.1016/s0014-4835(05)80122-8.
Wang L, Kondo M, Bill A. Glucose metabolism in cat outer retina. Effects of light and hyperoxia. Investig Ophthalmol Vis Sci. 1997;38(1):48–55.
Wang L, Törnquist P, Bill A. Glucose metabolism in pig outer retina in light and darkness. Acta Physiol Scand. 1997;160(1):75–81. https://doi.org/10.1046/j.1365-201X.1997.00030.x.
Yu DY, Cringle SJ. Retinal degeneration and local oxygen metabolism. Exp Eye Res. 2005;80(6):745–51. https://doi.org/10.1016/j.exer.2005.01.018.
Meschede IP, Ovenden NC, Seabra MC, Futter CE, Votruba M, Cheetham ME, et al. Symmetric arrangement of mitochondria:plasma membrane contacts between adjacent photoreceptor cells regulated by Opa1. Proc Natl Acad Sci. 2020;117(27):15684–93. https://doi.org/10.1073/pnas.2000304117.
Okawa H, Sampath AP, Laughlin SB, Fain GL. ATP Consumption by Mammalian Rod Photoreceptors in Darkness and in Light. Current Biology. 2008;18(24):1917–21. https://doi.org/10.1016/j.cub.2008.10.029.
Ingram NT, Fain GL, Sampath AP. Elevated energy requirement of cone photoreceptors. Proc Natl Acad Sci. 2020;117(32):19599–603. https://doi.org/10.1073/pnas.2001776117.
Warburg O. The Metabolism of Carcinoma Cells. The Journal of Cancer Research. 1925;9(1):148–63.
Cohen LH, Noell WK. Glucose Catabolism of Rabbit Retina Before and After Development of Visual Function. J Neurochem. 1960;5(3):253–76. https://doi.org/10.1111/j.1471-4159.1960.tb13363.x.
Winkler BS. Glycolytic and oxidative metabolism in relation to retinal function. J Gen Physiol. 1981;77(6):667–92. https://doi.org/10.1085/jgp.77.6.667.
Chinchore Y, Begaj T, Wu D, Drokhlyansky E, Cepko CL. Glycolytic reliance promotes anabolism in photoreceptors. eLife. 2017;6:e25946. https://doi.org/10.7554/eLife.25946.
Heiden MGV, Cantley LC, Thompson CB. Understanding the Warburg Effect: The Metabolic Requirements of Cell Proliferation. Science. 2009;324(5930):1029–33. https://doi.org/10.1126/science.1160809.
Young RW. The renewal of photoreceptor cell outer segments. J Cell Biol. 1967;33(1):61–72. https://doi.org/10.1083/jcb.33.1.61.
Young RW. The renewal of rod and cone outer segments in the rhesus monkey. J Cell Biol. 1971;49(2):303–18. https://doi.org/10.1083/jcb.49.2.303.
Hurley JB. Retina Metabolism and Metabolism in the Pigmented Epithelium: A Busy Intersection. Ann Rev Vis Sci. 2021;7(1):665–92. https://doi.org/10.1146/annurev-vision-100419-115156.
Nolan ND, Caruso SM, Cui X, Tsang SH. Renormalization of metabolic coupling treats age-related degenerative disorders: an oxidative RPE niche fuels the more glycolytic photoreceptors. Eye. 2022;36(2):278–83. https://doi.org/10.1038/s41433-021-01726-4.
Kanow MA, Giarmarco MM, Jankowski CS, Tsantilas K, Engel AL, Du J, et al. Biochemical adaptations of the retina and retinal pigment epithelium support a metabolic ecosystem in the vertebrate eye. eLife. 2017;6:e28899. https://doi.org/10.7554/eLife.28899.
Lin JB, Tsubota K, Apte RS. A glimpse at the aging eye. NPJ Aging Mech Dis. 2016;2(1):1–7. https://doi.org/10.1038/npjamd.2016.3.
Lam KCA, Chan ST, Chan HLH, Chan B. The effect of age on ocular blood supply determined by pulsatile ocular blood flow and color Doppler ultrasonography. Am J Optom Physiol Optic. 2003;80(4):305–11. https://doi.org/10.1097/00006324-200304000-00008.
Curcio CA, Millican CL, Bailey T, Kruth HS. Accumulation of cholesterol with age in human Bruch’s membrane. Investig Ophthalmol Vis Sci. 2001;42(1):265–74.
Wassell J, Davies S, Bardsley W, Boulton M. The photoreactivity of the retinal age pigment lipofuscin. J Biol Chem. 1999;274(34):23828–32. https://doi.org/10.1074/jbc.274.34.23828.
Caprara C, Grimm C. From oxygen to erythropoietin: Relevance of hypoxia for retinal development, health and disease. Prog Retin Eye Res. 2012;31(1):89–119. https://doi.org/10.1016/j.preteyeres.2011.11.003.
van Vliet T, Casciaro F, Demaria M. To breathe or not to breathe: Understanding how oxygen sensing contributes to age-related phenotypes. Ageing Res Rev. 2021;67:101267. https://doi.org/10.1016/j.arr.2021.101267.
Dimopoulos IS, Freund PR, Redel T, Dornstauder B, Gilmour G, Sauvé Y. Changes in rod and cone-driven oscillatory potentials in the aging human retina. Investig Ophthalmol Vis Sci. 2014;55(8):5058–73. https://doi.org/10.1167/iovs.14-14219.
Weleber RG. The effect of age on human cone and rod ganzfeld electroretinograms. Investig Ophthalmol Vis Sci. 1981;20(3):392–9.
Kaelin WG, Ratcliffe PJ. Oxygen Sensing by Metazoans: The Central Role of the HIF Hydroxylase Pathway. Mol Cell. 2008;30(4):393–402. https://doi.org/10.1016/j.molcel.2008.04.009.
Mammadzada P, Corredoira PM, André H. The role of hypoxia-inducible factors in neovascular age-related macular degeneration: a gene therapy perspective. Cell Mol Life Sci. 2020;77(5):819–33. https://doi.org/10.1007/s00018-019-03422-9.
Maynard MA, Ohh M. Von Hippel-Lindau tumor suppressor protein and hypoxia-inducible factor in kidney cancer. Am J Nephrol. 2004;24(1):1–13. https://doi.org/10.1159/000075346.
Lange C, Heynen SR, Tanimoto N, Thiersch M, Le YZ, Meneau I, et al. Normoxic Activation of Hypoxia-Inducible Factors in Photoreceptors Provides Transient Protection against Light-Induced Retinal Degeneration. Investig Ophthalmol Vis Sci. 2011;52(8):5872–80. https://doi.org/10.1167/iovs.11-7204.
Barben M, Ail D, Storti F, Klee K, Schori C, Samardzija M, et al. Hif1a inactivation rescues photoreceptor degeneration induced by a chronic hypoxia-like stress. Cell Death Differ. 2018;25(12):2071–85. https://doi.org/10.1038/s41418-018-0094-7.
Mitchell P, Liew G, Gopinath B, Wong TY. Age-related macular degeneration. Lancet. 2018;392(10153):1147–59. https://doi.org/10.1016/S0140-6736(18)31550-2.
Chen TW, Wardill TJ, Sun Y, Pulver SR, Renninger SL, Baohan A, et al. Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature. 2013;499(7458):295–300. https://doi.org/10.1038/nature12354.
Calvert PD, Krasnoperova NV, Lyubarsky AL, Isayama T, Nicoló M, Kosaras B, et al. Phototransduction in transgenic mice after targeted deletion of the rod transducin a-subunit. Proc Natl Acad Sci. 2000;97(25):13913–8. https://doi.org/10.1073/pnas.250478897.
San Martín A, Ceballo S, Ruminot I, Lerchundi R, Frommer WB, Barros LF. A Genetically Encoded FRET Lactate Sensor and Its Use To Detect the Warburg Effect in Single Cancer Cells. PLoS ONE. 2013;8(2). https://doi.org/10.1371/journal.pone.0057712.
Delgado MG, Oliva C, López E, Ibacache A, Galaz A, Delgado R, et al. Chaski, a novel Drosophila lactate/pyruvate transporter required in glia cells for survival under nutritional stress. Sci Rep. 2018;8(1):1186. https://doi.org/10.1038/s41598-018-19595-5.
Ovens MJ, Davies AJ, Wilson MC, Murray CM, Halestrap AP. AR-C155858 is a potent inhibitor of monocarboxylate transporters MCT1 and MCT2 that binds to an intracellular site involving transmembrane helices 7–10. Biochem J. 2010;425(3):523–30. https://doi.org/10.1042/BJ20091515.
Takanaga H, Chaudhuri B, Frommer WB. GLUT1 and GLUT9 as major contributors to glucose influx in HepG2 cells identified by a high sensitivity intramolecular FRET glucose sensor. Biochim Biophys Acta (BBA) Biomembr. 2008;1778(4):1091–9. https://doi.org/10.1016/j.bbamem.2007.11.015.
Baeza-Lehnert F, Saab AS, Gutiérrez R, Larenas V, Díaz E, Horn M, et al. Non-Canonical Control of Neuronal Energy Status by the Na+ Pump. Cell Metabolism. 2019;29(3):668-680.e4. https://doi.org/10.1016/j.cmet.2018.11.005.
Bittner CX, Valdebenito R, Ruminot I, Loaiza A, Larenas V, Sotelo-Hitschfeld T, et al. Fast and reversible stimulation of astrocytic glycolysis by K+ and a delayed and persistent effect of glutamate. J Neurosci. 2011;31(12):4709–13. https://doi.org/10.1523/JNEUROSCI.5311-10.2011.
Bittner CX, Loaiza A, Ruminot I, Larenas V, Sotelo-Hitschfe T, Gutiérrez R, et al. High Resolution Measurement of the Glycolytic Rate. Front Neuroenerg. 2010;2. https://doi.org/10.3389/fnene.2010.00026.
Koveal D, Díaz-García CM, Yellen G. Fluorescent Biosensors for Neuronal Metabolism and the Challenges of Quantitation. Curr Opin Neurobiol. 2020;63:111–21. https://doi.org/10.1016/j.conb.2020.02.011.
Todorova V, Merolla L, Karademir D, Wögenstein GM, Behr J, Ebner LJA, et al. Retinal Layer Separation (ReLayS) method enables the molecular analysis of photoreceptor segments and cell bodies, as well as the inner retina. Sci Rep. 2022;12(1):20195. https://doi.org/10.1038/s41598-022-24586-8.
Mantych GJ, Hageman GS, Devaskar SU. Characterization of glucose transporter isoforms in the adult and developing human eye. Endocrinology. 1993;133(2):600–7. https://doi.org/10.1210/en.133.2.600.
Watanabe T, Matsushima S, Okazaki M, Nagamatsu S, Hirosawa K, Uchimura H, et al. Localization and ontogeny of GLUT3 expression in the rat retina. Dev Brain Res. 1996;94(1):60–6. https://doi.org/10.1016/0165-3806(96)00044-2.
Le YZ, Zheng L, Zheng W, Ash JD, Agbaga MP, Zhu M, et al. Mouse opsin promoter-directed Cre recombinase expression in transgenic mice. Mol Vis. 2006;12:389–98.
Kunimi H, Lee D, Ibuki M, Katada Y, Negishi K, Tsubota K, et al. Inhibition of the HIF-1a/BNIP3 pathway has a retinal neuroprotective effect. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 2021;35(8):e21829. https://doi.org/10.1096/fj.202100572R.
Rikka S, Quinsay MN, Thomas RL, Kubli DA, Zhang X, Murphy AN, et al. Bnip3 impairs mitochondrial bioenergetics and stimulates mitochondrial turnover. Cell Death Differ. 2011;18(4):721–31. https://doi.org/10.1038/cdd.2010.146.
Diaz F, Thomas CK, Garcia S, Hernandez D, Moraes CT. Mice lacking COX10 in skeletal muscle recapitulate the phenotype of progressive mitochondrial myopathies associated with cytochrome c oxidase deficiency. Hum Mol Genet. 2005;14(18):2737–48. https://doi.org/10.1093/hmg/ddi307.
Custódio TF, Paulsen PA, Frain KM, Pedersen BP. Structural comparison of GLUT1 to GLUT3 reveal transport regulation mechanism in sugar porter family. Life Sci Alliance. 2021;4(4). https://doi.org/10.26508/lsa.202000858.
Watanabe T, Nagamatsu S, Matsushima S, Kirino T, Uchimura H. Colocalization of GLUT3 and Choline Acetyltransferase Immunoreactivity in the Rat Retina. Biochem Biophys Res Commun. 1999;256(3):505–11. https://doi.org/10.1006/bbrc.1999.0369.
Rodriguez-Rodriguez P, Fernandez E, Almeida A, Bolaños JP. Excitotoxic stimulus stabilizes PFKFB3 causing pentose-phosphate pathway to glycolysis switch and neurodegeneration. Cell Death Differ. 2012;19(10):1582–9. https://doi.org/10.1038/cdd.2012.33.
Ames A, Li YY, Heher EC, Kimble CR. Energy metabolism of rabbit retina as related to function: high cost of Na+ transport. J Neurosci. 1992;12(3):840–53. https://doi.org/10.1523/JNEUROSCI.12-03-00840.1992.
Contreras-Baeza Y, Sandoval PY, Alarcón R, Galaz A, Cortés-Molina F, Alegría K, et al. Monocarboxylate transporter 4 (MCT4) is a high affinity transporter capable of exporting lactate in high-lactate microenvironments. J Biol Chem. 2019;294(52):20135–47. https://doi.org/10.1074/jbc.RA119.009093.
Martínez-Reyes I, Chandel NS. Mitochondrial TCA cycle metabolites control physiology and disease. Nat Commun. 2020;11(1):102. https://doi.org/10.1038/s41467-019-13668-3.
Arnold PK, Finley LWS. Regulation and function of the mammalian tricarboxylic acid cycle. J Biol Chem. 2022;102838. https://doi.org/10.1016/j.jbc.2022.102838.
Supplie LM, Düking T, Campbell G, Diaz F, Moraes CT, Götz M, et al. Respiration-Deficient Astrocytes Survive As Glycolytic Cells In Vivo. J Neurosci. 2017;37(16):4231–42. https://doi.org/10.1523/JNEUROSCI.0756-16.2017.
Fünfschilling U, Supplie LM, Mahad D, Boretius S, Saab AS, Edgar J, et al. Glycolytic oligodendrocytes maintain myelin and long-term axonal integrity. Nature. 2012;485(7399):517–21. https://doi.org/10.1038/nature11007.
Diaz F, Garcia S, Padgett KR, Moraes CT. A defect in the mitochondrial complex III, but not complex IV, triggers early ROS-dependent damage in defined brain regions. Hum Mol Genet. 2012;21(23):5066–77. https://doi.org/10.1093/hmg/dds350.
Fukui H, Diaz F, Garcia S, Moraes CT. Cytochrome c oxidase deficiency in neurons decreases both oxidative stress and amyloid formation in a mouse model of Alzheimer’s disease. Proc Natl Acad Sci. 2007;104(35):14163–8. https://doi.org/10.1073/pnas.0705738104.
Bisbach CM, Hass DT, Robbings BM, Rountree AM, Sadilek M, Sweet IR, et al. Succinate Can Shuttle Reducing Power from the Hypoxic Retina to the O2-Rich Pigment Epithelium. Cell Rep. 2020;31(5):107606. https://doi.org/10.1016/j.celrep.2020.107606.
Gao J, Cheon K, Nusinowitz S, Liu Q, Bei D, Atkins K, et al. Progressive photoreceptor degeneration, outer segment dysplasia, and rhodopsin mislocalization in mice with targeted disruption of the retinitis pigmentosa-1 (Rp1) gene. Proc Natl Acad Sci. 2002;99(8):5698–703. https://doi.org/10.1073/pnas.042122399.
Meschede IP, Burgoyne T, Tolmachova T, Seabra MC, Futter CE. Chronically shortened rod outer segments accompany photoreceptor cell death in Choroideremia. PLoS ONE. 2020;15(11):e0242284. https://doi.org/10.1371/journal.pone.0242284.
Rowe AA, Patel PD, Gordillo R, Wert KJ. Replenishment of TCA cycle intermediates provides photoreceptor resilience against neurodegeneration during progression of retinitis pigmentosa. JCI Insight. 2021;6(17). https://doi.org/10.1172/jci.insight.150898.
Jiang K, Mondal AK, Adlakha YK, Gumerson J, Aponte A, Gieser L, et al. Multi-omics analyses reveal early metabolic imbalance and mitochondrial stress in neonatal photoreceptors leading to cell death in Pde6brd1/rd1 mouse model of retinal degeneration. Hum Mol Genet. 2022;31(13):2137–54. https://doi.org/10.1093/hmg/ddac013.
Cideciyan AV, Jacobson SG, Roman AJ, Sumaroka A, Wu V, Charng J, et al. Rod function deficit in retained photoreceptors of patients with class B Rhodopsin mutations. Sci Rep. 2020;10(1):12552. https://doi.org/10.1038/s41598-020-69456-3.
Milam AH, Li ZY, Fariss RN. Histopathology of the human retina in retinitis pigmentosa. Prog Retin Eye Res. 1998;17(2):175–205. https://doi.org/10.1016/s1350-9462(97)00012-8.
Ebner LJA, Samardzija M, Storti F, Todorova V, Karademir D, Behr J, et al. Transcriptomic analysis of the mouse retina after acute and chronic normobaric and hypobaric hypoxia. Sci Rep. 2021;11(1):16666. https://doi.org/10.1038/s41598-021-96150-9.
Nork TM, Millecchia LL, Strickland BD, Linberg JV, Chao GM. Selective Loss of Blue Cones and Rods in Human Retinal Detachment. Arch Ophthalmol. 1995;113(8):1066–73. https://doi.org/10.1001/archopht.1995.01100080118039.
Joyal JS, Sun Y, Gantner ML, Shao Z, Evans LP, Saba N, et al. Retinal lipid and glucose metabolism dictates angiogenesis through the lipid sensor Ffar1. Nat Med. 2016;22(4):439–45. https://doi.org/10.1038/nm.4059.
Pearsall EA, Cheng R, Zhou K, Takahashi Y, Matlock HG, Vadvalkar SS, et al. PPARa is essential for retinal lipid metabolism and neuronal survival. BMC Biol. 2017;15(1):113. https://doi.org/10.1186/s12915-017-0451-x.
Haase VH, Glickman JN, Socolovsky M, Jaenisch R. Vascular tumors in livers with targeted inactivation of the von Hippel-Lindau tumor suppressor. Proc Natl Acad Sci U S A. 2001;98(4):1583–8. https://doi.org/10.1073/pnas.98.4.1583.
Ryan HE, Poloni M, McNulty W, Elson D, Gassmann M, Arbeit JM, et al. Hypoxia-inducible factor-1alpha is a positive factor in solid tumor growth. Cancer Res. 2000;60(15):4010–5.
Gruber M, Hu CJ, Johnson RS, Brown EJ, Keith B, Simon MC. Acute postnatal ablation of Hif-2alpha results in anemia. Proc Natl Acad Sci U S A. 2007;104(7):2301–6. https://doi.org/10.1073/pnas.0608382104.
Samardzija M, Wenzel A, Naash M, Remé CE, Grimm C. Rpe65 as a modifier gene for inherited retinal degeneration. Eur J Neurosci. 2006;23(4):1028–34. https://doi.org/10.1111/j.1460-9568.2006.04639.x.
Danciger M, Matthes MT, Yasamura D, Akhmedov NB, Rickabaugh T, Gentleman S, et al. A QTL on distal chromosome 3 that influences the severity of light-induced damage to mouse photoreceptors. Mamm Genome. 2000;11(6):422–7. https://doi.org/10.1007/s003350010081.
Samardzija M, Wenzel A, Aufenberg S, Thiersch M, Remé C, Grimm C. Differential role of Jak-STAT signaling in retinal degenerations. FASEB J. 2006;20(13):2411–3. https://doi.org/10.1096/fj.06-5895fje.
Stobart JL, Ferrari KD, Barrett MJP, Stobart MJ, Looser ZJ, Saab AS, et al. Long-term In Vivo Calcium Imaging of Astrocytes Reveals Distinct Cellular Compartment Responses to Sensory Stimulation. Cereb Cortex. 2018;28(1):184–98. https://doi.org/10.1093/cercor/bhw366.
Beltran WA, Boye SL, Boye SE, Chiodo VA, Lewin AS, Hauswirth WW, et al. rAAV2/5 gene-targeting to rods:dose-dependent efficiency and complications associated with different promoters. Gene Ther. 2010;17(9):1162–74. https://doi.org/10.1038/gt.2010.56.
Komatsubara AT, Matsuda M, Aoki K. Quantitative analysis of recombination between YFP and CFP genes of FRET biosensors introduced by lentiviral or retroviral gene transfer. Sci Rep. 2015;5(1):13283. https://doi.org/10.1038/srep13283.
Kim M, Eleftheriou A, Ravotto L, Weber B, Rivlin M, Navon G, et al. What do we know about dynamic glucose-enhanced (DGE) MRI and how close is it to the clinics? Horizon 2020 GLINT consortium report. Magn Reson Mater Phys Biol Med. 2022;35(1):87–104. https://doi.org/10.1007/s10334-021-00994-1.
Reid CA, Ertel KJ, Lipinski DM. Improvement of Photoreceptor Targeting via Intravitreal Delivery in Mouse and Human Retina Using Combinatory rAAV2 Capsid Mutant Vectors. Investig Ophthalmol Vis Sci. 2017;58(14):6429–39. https://doi.org/10.1167/iovs.17-22281.
Kulkarni M, Schubert T, Baden T, Wissinger B, Euler T, Paquet-Durand F. Imaging Ca2+ Dynamics in Cone Photoreceptor Axon Terminals of the Mouse Retina. J Visualized Exp. 2015;99:e52588. https://doi.org/10.3791/52588.
Mayrhofer JM, Haiss F, Haenni D, Weber S, Zuend M, Barrett MJP, et al. Design and performance of an ultra-flexible two-photon microscope for in vivo research. Biomed Optic Express. 2015;6(11):4228–37. https://doi.org/10.1364/BOE.6.004228.
Pologruto TA, Sabatini BL, Svoboda K. ScanImage: Flexible software for operating laser scanning microscopes. Biomed Eng OnLine. 2003;2(1):13. https://doi.org/10.1186/1475-925X-2-13.
Li CH, Lee CK. Minimum cross entropy thresholding. Pattern Recog. 1993;26(4):617–25. https://doi.org/10.1016/0031-3203(93)90115-D.
Perez-Riverol Y, Bai J, Bandla C, García-Seisdedos D, Hewapathirana S, Kamatchinathan S, et al. The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 2021;50(D1):D543–52. https://doi.org/10.1093/nar/gkab1038.
Käll L, Canterbury JD, Weston J, Noble WS, MacCoss MJ. Semi-supervised learning for peptide identification from shotgun proteomics datasets. Nat Methods. 2007;4(11):923–5. https://doi.org/10.1038/nmeth1113.
Navarro P, Trevisan-Herraz M, Bonzon-Kulichenko E, Núñez E, Martínez-Acedo P, Pérez-Hernández D, et al. General statistical framework for quantitative proteomics by stable isotope labeling. J Proteome Res. 2014;13(3):1234–47. https://doi.org/10.1021/pr4006958.
Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015;43(Web Server issue):W566–70. https://doi.org/10.1093/nar/gkv468.
Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–9. https://doi.org/10.1093/bioinformatics/btw313.
Gu Z, Gu L, Eils R, Schlesner M, Brors B. circlize Implements and enhances circular visualization in R. Bioinformatics. 2014;30(19):2811–2. https://doi.org/10.1093/bioinformatics/btu393.
Goedhart J, Luijsterburg MS. VolcaNoseR is a web app for creating, exploring, labeling and sharing volcano plots. Sci Rep. 2020;10(1):20560. https://doi.org/10.1038/s41598-020-76603-3.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73. https://doi.org/10.1038/ng1180.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25. https://doi.org/10.1016/j.cels.2015.12.004.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics (Oxford, England). 2011;27(12):1739–40. https://doi.org/10.1093/bioinformatics/btr260.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Kast B, Schori C, Grimm C. Hypoxic preconditioning protects photoreceptors against light damage independently of hypoxia inducible transcription factors in rods. Exp Eye Res. 2016;146:60–71. https://doi.org/10.1016/j.exer.2015.12.008.
Acknowledgements
We thank the Viral Vector Facility (VVF) of the Neuroscience Center Zurich (ZNZ) and the Laboratory Animal Services Center (LASC) of the University of Zurich.
Funding
This work was supported by Swiss National Science Foundation grants 31003A_173008 and 310030_200798.
Author information
Authors and Affiliations
Contributions
V.T., M.S., B.W. and C.G. conceived the study. V.T. performed the TPLSM experiments. V.T. and L.R. analysed and visualized the TPLSM data. V.T. and S.M.H. performed the proteomics analysis, V.T., D.K., L.J.A.E. and S.M.H. analysed and visualised the proteomics data. V.T., M.F.S. and C.I. performed the morphological analysis. V.T. and M.F.S. performed and analysed the immuno- and enzymatic histochemistry. V.T. and S.N. performed the mtDNA copy number assay. V.T. and L.M. performed the semi-quantitative real-time PCR analysis. V.T. performed and analysed the ERG experiments. All authors contributed to project design, organization of the experiments, interpretation of results and writing of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Figure S1.
Light-induced calcium response in rods during TPLSM. a, Expression of the calcium sensor GCaMP6s (green) in cone arrestin (cARR, red) negative photoreceptors after subretinal AAV application. White arrows: cARR-positive but GCaMP6s negative cells. b, Preparation of acute retinal slices from flat mounted half retinas and imaging of GCaMP6s by TPLSM. c, Representative TPLSM micrographs illustrating the drop of intracellular calcium levels during imaging. Left: Max intensity of whole image series. Middle and right: intensity weighted images at 0 s (middle) and 4 s (right) of imaging. White arrows indicate regions with reduced GCaMP6s signal after 4 s of imaging. Scale bars, 50 μm. d, Representative GCaMP6s traces in the outer plexiform layer (OPL), outer nuclear layer (ONL), and photoreceptor segments (PS) of wild type mice during 6 s of TPLSM. e, GCaMP6s signal (mean \(\pm\) SD) in the OPL of Gnat1a-/- mice (12 retinal slices from 4 mice) at 2, 4, and 6 s during TPLSM.
Additional file 2: Figure S2.
Variability and mosaicism of Cre levels in \(rod^{\varDelta\ Vhl}\) mice. a, Relative normalized abundance of Cre, SAG, and GAPDH in the outer nuclear layer (ONL) of \(rod^{\varDelta\ Vhl}\) mice at 2.5 months of age. Shown are means (normalized to 1) \(\pm\) SD and individual data points. The samples with the highest (red) and lowest (blue) Cre level are indicated. Statistics: Bartlett's test for homoscedasticity. b, Immunofluorescence labeling for Cre in \(rod^{\varDelta\ Vhl}\) mice at 2.5 months of age. White box: magnification of the ONL region. Scale bar, 50 μm. Panorama: Scale bar, 500 μm.
Additional file 3: Figure S3.
Expression of mitochondrial genes in the ONL and PS of \(rod^{\varDelta\ Vhl}\) mice. a, Relative expression of Vdac1 and Tomm20 in ONL and PS samples from 2.5 months old \(rod^{\varDelta\ Vhl}\) and ctrl mice. b, Relative expression of the mt-DNA encoded genes mt-Co2 and mt-Nd1 in ONL and PS samples from 2.5 months old \(rod^{\varDelta\ Vhl}\) and ctrl mice. N = 6 mice per genotype. Statistics: Mann-Whitney nonparametric test.
Additional file 4: Figure S4.
Retinal function of \(rod^{\varDelta\ Cox10}\) mice. a, Scotopic single flash ERG responses to light stimuli with increasing intensities (top to bottom) of \(rod^{\varDelta\ Cox10}\) (orange) and ctrl (black) mice at 6 (left) and 12 (right) months of age. Shown are averaged traces. b,c, Scotopic a-wave (b) and b-wave (c) amplitudes as a function of stimulus intensity derived from (a). Shown are means \(\pm\) SD. **: p \(\leq\) 0.01. ****: p \(\leq\) 0.0001. d, Photopic single flash ERG responses to light stimuli with increasing intensities (top to bottom) of \(rod^{\varDelta\ Cox10}\) (orange) and ctrl (black) mice at 6 (left) and 12 (right) months of age. Shown are averaged traces. e, Photopic b-wave amplitudes as a function of stimulus intensity derived from (d). Shown are means \(\pm\) SD. N = 5 (6-month-old \(rod^{\varDelta\ Cox10}\) mice); N = 6 (6- and 12-month-old ctrl mice; 12-month-old \(rod^{\varDelta\ Cox10}\) mice).
Additional file 5: Table S1.
Top 25 upregulated proteins in the PS and ONL.
Additional file 6: Table S2.
Top 50 differentially regulated proteins in the ONL of \(rod^{\varDelta\ Vhl}\) mice.
Additional file 7: Table S3.
Top 50 differentially regulated proteins in the PS of \(rod^{\varDelta\ Vhl}\) mice.
Additional file 8: Table S4.
Differentially regulated mitochondrial proteins in the PS of \(rod^{\varDelta\ Vhl}\) mice.
Additional file 9: Table S5.
Proteomic analysis of ONL and PS from \(rod^{\varDelta\ Vhl}\) and ctrl mice; includes UniProt IDs, number of unique peptides, ensembl IDs, gene symbols and short description, abundance ratios, adj-p values and raw data.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Todorova, V., Stauffacher, M.F., Ravotto, L. et al. Deficits in mitochondrial TCA cycle and OXPHOS precede rod photoreceptor degeneration during chronic HIF activation. Mol Neurodegeneration 18, 15 (2023). https://doi.org/10.1186/s13024-023-00602-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13024-023-00602-x