Abstract
Endothelin-1 (ET-1) is a vasoactive peptide that is elevated in aqueous humor as well as circulation of primary open angle glaucoma (POAG) patients. ET-1 has been shown to promote degeneration of optic nerve axons and apoptosis of retinal ganglion cells (RGCs), however, the precise mechanisms are still largely unknown. In this study, RNA-seq analysis was used to assess changes in ET-1 mediated gene expression in primary RGCs, which revealed that 23 out of 156 differentially expressed genes (DEGs) had known or predicted mitochondrial function, of which oxidative phosphorylation emerged as the top-most enriched pathway. ET-1 treatment significantly decreased protein expression of key mitochondrial genes including cytochrome C oxidase copper chaperone (COX17) and ATP Synthase, H+ transporting, Mitochondrial Fo Complex (ATP5H) in primary RGCs and in vivo following intravitreal ET-1 injection in rats. A Seahorse ATP rate assay revealed a significant decrease in the rate of mitochondrial ATP production following ET-1 treatment. IOP elevation in Brown Norway rats showed a trend towards decreased expression of ATP5H. Our results demonstrate that ET-1 produced a decrease in expression of vital components of mitochondrial electron transport chain, which compromise bioenergetics and suggest a mechanism by which ET-1 promotes neurodegeneration of RGCs in glaucoma.
Similar content being viewed by others
Introduction
Glaucoma is an optic neuropathy with an approximate prevalence of 60.5 million people worldwide and is projected to reach 111 million by 20401,2.The disease is commonly associated with elevated intraocular pressure (IOP), accompanied by optic nerve degeneration and loss of retinal ganglion cells (RGCs)3,4. RGC death via apoptosis is a culminating event in the pathophysiology of glaucoma, stemming from optic nerve axonal injury, leading to visual field loss. Elevated IOP is a major risk factor in primary open-angle-glaucoma and current therapeutic approaches are aimed at lowering IOP with medications, laser treatment, or surgery5,6. However in some patients, the progression of the disease continues to occur slowly7 despite lowering IOP, hence there is a compelling need for neuroprotection of RGCs and optic nerve axons and as an additional therapeutic modality. The molecular changes occurring specifically in the RGCs, during the progression of glaucoma, contributing to neurodegeneration, are still poorly understood; hence identifying new therapeutic targets could provide more efficacious neuroprotective treatments.
ET-1 is a 21 amino acid vasoactive peptide that acts through two G-protein coupled receptors namely, ETA and ETB receptors, to produce diverse effects in various ocular tissues8,9,10,11. A growing body of evidence suggests that endothelins and their receptors are major contributors to neuronal damage in glaucoma11,12,13,14,15,16. The role of endothelins in glaucomatous neurodegeneration has been the subject of several review articles12,17,18. However, the detailed cellular and molecular mechanisms that contribute to ET-1 mediated-neurodegeneration are not completely understood. ET-1 concentrations were also found to be increased in several animal models of glaucoma19,20,21. Both intravitreal and peribulbar administration of ET-1 has been shown to produce axon loss, RGC loss and disruption of nerve fiber layer11,22,23.
Several studies point to the role of mitochondrial dysfunction and oxidative stress as contributors to glaucomatous damage in animal models of glaucoma24,25,26,27. There is an abundance of mitochondria in retinal ganglion cells (RGCs) within the soma region and the initial unmyelinated axons (anterior to the lamina cribrosa) of the optic nerve and synaptic terminals as opposed to the post-laminar and myelinated region. Hence, RGCs are typically more sensitive to mitochondrial dysfunction than other neuronal populations28. Retrograde and anterograde transport of mitochondria which occurs continuously along the axons of neurons, is vital for synaptic function. Hence, any disturbance of axonal transport due to mitochondrial abnormalities and a compromise in the mitochondrial function can lead to severe metabolic crises because of energy depletion29. The impairment in axonal transport can be observed commonly in diseases such as Amyotrophic Lateral Sclerosis and Leber’s Hereditary Optic Neuropathy30. Mitochondrial function decreases with age and as RGCs are largely dependent on their mitochondria, the prevalence of primary open angle glaucoma (POAG) increases with age possibly due to a decline in mitochondrial function31. Impaired mitochondrial dysfunction is a phenomenon observed in most, if not all, neurodegenerative disorders32,33,34,35. ET-1 elevates reactive oxygen species by acting through the ETB receptor, which could be one mechanism by which ET-1 alters mitochondrial function36.
The effect of ET-1 on mitochondrial gene expression and function and their potential impact on neurodegeneration of RGCs has not been adequately explored. The purpose of this study was to use RNA-seq to investigate ET-1 mediated changes in gene expression of mitochondrial genes and proteins in primary RGCs treated with ET-1, as well as in vivo following either ET-1 administration or IOP elevation in rat eyes.
Results
Characterization, purity and yield of isolated primary RGCs
By sequential double immunopanning with the Thy1.1 antibody, we isolated approximately 20,000 to 40,000 RGCs per retina, with a yield of 40% to 50% of the RGCs in post-natal day 4–7 (P4-P7) rat retinas35. To confirm the purity of the RGCs, immunostaining using RGC specific markers, including, RBPMS (Fig. 1a) and Brn3a (Fig. 1b) was carried out and GFAP, an astrocyte marker was used as a negative control (Fig. 1c). It was found that approximately 90 to 95% of the cells were immunopositive for RBPMS and Brn3a and hence were characterized as RGCs, whereas 2% of the cells were immunoreactive for GFAP. Thus, the immunopanning technique yielded a highly enriched RGC culture.
Translatome analysis of primary RGCs following ET-1 treatment (100 nM)
Primary cultures of RGCs were isolated and treated with ET-1 for 24 h, following which polysomal RNA was isolated. RNA-seq analysis of polysomal RNA from untreated (control) and ET-1-treated RGCs was performed and bioinformatics analysis was carried out. From the STRING network analysis (Fig. 2), several mitochondrially-relevant genes which had a significant change in expression were identified. Gene Ontology (GO) classification with statistical overrepresentation tests indicate a significant enrichment in the differentially expressed genes (DEG) set for cellular components, including mitochondrion, endocytic vesicle membrane and cytoplasmic stress granules. The DEG set was also enriched for specific biological processes, including cation transport and metabolic processes. Twenty-three genes out of the 156 DEGs were identified to have a known or predicted mitochondrial function. Of the 23 genes, 14 of which are represented in the STRING network (indicated by the mitochondrial symbol on the node) are defined by some connected nodes. Four of the connected, mitochondrially-related genes were in the top 20 up-regulated and down-regulated genes in the list of 156 DEGs.
From the STRING network, mitochondrially-relevant genes were identified within the network and the direction of fold regulation was overlaid using red (up-regulated) and green (down-regulated) borders. It was found that 156 genes were differentially expressed following ET-1 treatment in RGCs. The thickness of the color indicator denotes the ranking of the gene: thick indicators are ranked in the top or bottom 20 of the ranked DEG list; medium borders are in the top / bottom 20 genes based only on magnitude of fold change; thin lines just indicate direction of fold change. Analysis of differential gene expression revealed a significant change in expression of several genes involved in mitochondrial function, oxidative metabolism and cell survival pathways. An increase in an expression of mitochondrial proteins cytochrome c oxidase copper chaperone (COX17) (3- fold), and ATP synthase, H+ transporting mitochondrial F0 complex (ATP5H) (3-fold) was found. On the other hand, a decrease in expression of glucose-6-phosphate dehydrogenase (5-fold), Forkhead box O1 (FOXO1) (986-fold), mitogen- activated protein kinase kinase kinase 11(Map3k11) (40-fold) and modulator of apoptosis (Moap1) (6255-fold) was observed.
Translatome changes in primary RGCs following ET-1 treatment
To explore the translatome differences in primary retinal RGCs between the control and ET-1 treated samples, a heatmap was generated showing log2 (fold-change)> 1.5. Hierarchical clustering of the Pearson correlation coefficients was performed using the 156 DEGs to determine if the gene expression profiles were different between the two groups (Fig. 3a). The results showed that the replicates cluster together based on similarities in gene expression, although there was some variability observed in the top genes between the replicates. Further, two clear sets of genes with increased expression as well as genes with decreased expression following ET-1 treatment emerged. In addition, principal component analysis (PCA) of the normalized gene expression values for the top 156 DEGs showed that axes 1 and 2 accounted for 76.3% and 19.4% of the total variation, respectively (Fig. 3b). The sign (+/-) of PC1 distinguishes controls from ET-1 treated replicates, where PC1 is positive for ET-1 treatment and negative for the controls. Ingenuity Pathway Analysis (IPA) canonical pathway analysis of the gene set revealed that oxidative phosphorylation pathway is the top-most enriched pathway within the DEGs following ET-1 treatment (Fig. 3c).
Validation of DEGs with qPCR
To validate RNA-seq findings, qPCR analysis was performed on several selected genes associated with either mitochondrial function, neurodegeneration or cell survival. A significant decrease in gene expression of Atp5h (30%) and Moap-1 (17%) was observed following treatment of primary RGCs with ET-1 (Fig. 4). On the other hand, an increasing trend in the expression of FOXO1 (31%) and Map3K11 (179%) genes and a decreasing trend in the expression of Cox17 (17%) and G6PD (15%) genes were observed in ET-1 treated RGCs (Fig. 4). However, these were not statistically significant. The expression pattern observed by conventional qPCR was different in comparison to the RNA-seq data for some DEGs including Atp5h, Cox17 FoxO1, and Map3k11. There was no appreciable change in the expression of mitochondrial ribosomal protein L30 (Mrpl30) gene following ET-1 treatment when analysed by qPCR (Fig. 4).
ET-1 treatment causes a decline in key mitochondrial proteins COX17 and ATP5H in primary RGCs
Mitochondrial dysfunction and oxidative damage are associated with many neurodegenerative disorders including glaucoma. Hence, we focused our attention on the expression of some mitochondrial genes that are important components of the mitochondrial oxidative phosphorylation pathway. To observe changes at the protein level, immunostaining for ATP5H and COX17 was carried out (Fig. 5a,b). Staining for both ATP5H and COX17 was observed both in the soma and axonal processes of RGCs.
Immunocytochemical analysis revealed that ET-1 treatment produced a decrease in immunostaining of key mitochondrial proteins ATP5H and COX17, indicating that ET-1 has the potential to compromise the bioenergetics of RGCs. Quantitation of the immunofluorescence of both ATP5H and COX17 revealed a decrease in integrated density/100 Brn3a positive cells following ET-1 treatment. (Fig. 5c,d) However, immunostaining of only COX17 showed a significant (p<0.05) decrease in integrated density/100 Brn3a positive cells.
A significant decrease in mitochondrial rate of ATP production was observed following ET-1 treatment in primary RGCs
Since a decrease in immunostaining for ATP5H and COX17 was found in RGCs treated with ET-1, a Seahorse XF ATP rate assay was carried out to determine if a functional decline in mitochondrial ATP production was observed. As seen in Fig. 6a, the OCR differed at the 4 h and 24 h timepoints of ET-1 treatment. At the 4 hour time point, RGCs had a higher basal level of OCR (35 pmol/min), compared to the 24 hour time point (20 pmol/min). At both 4 h and 24 h time points of ET-1 treatment, the OCR was reduced, compared to the corresponding untreated controls. ET-1 treatment produced a significant decrease in the rate of mitochondrial ATP production at both 4 h and 24 h time points following treatment with 100 nM ET-1 (Fig. 6b). The glycolytic ATP production rate showed a significant decrease following 4 h and a decreasing trend at 24 h following ET-1 treatment in RGCs. When the ratio of the mitochondrial to glycolytic ATP production was calculated, a significant decrease of the ratio was found at both 4 h and 24 h time points of ET-1 treatment (Fig. 6c).
Decrease in expression of mitochondrial proteins COX17 and ATP5H following intravitreal administration of ET-1
To determine if ET-1-mediated changes in ATP5H and COX17 were also observed in vivo, intravitreal injection of ET-1 was carried out in one eye of Brown Norway rats. The contralateral eye was injected with the vehicle used to dissolve ET-1. Rat retina section from vehicle-injected rat eyes showed immunoreactivity for ATP5H and COX17, which was assessed in multiple retinal layers including the nerve fiber layer (NFL), ganglion cell layer (GCL) and inner plexiform layer (IPL).
In the vehicle treated eyes, robust immunostaining for ATP5H was detected in the retinal ganglion cells, inner plexiform layer and outer plexiform layer (Fig. 7a). On the other hand, immunostaining for COX17 was found predominantly in the nerve fiber layer, retinal ganglion cells and inner plexiform layer in vehicle treated rat eyes (Fig. 7a). A semi-quantitative analysis measuring the change in integrated intensity was also performed using ImageJ software. Projections from confocal image z-stacks were analyzed for integrated density of ATP5H and COX17 expression in ET-1 treated and vehicle-treated eyes to determine fold change in expression. Compared to the vehicle-treated controls, retinas obtained from ET-1 treated eyes showed a statistically significant decrease in ATP5H expression mainly in the GCL (RFI = 0.4818 ± 0.1328, n = 6, p<0.01) and IPL (RFI = 0.5633 ± 0.1333, n = 6, p < 0.01) (Fig. 7b). Expression of ATP5H was co-localized with Brn3a (a RGC specific marker) in RGCs (Fig. 7b). Additionally, a modest decrease in immunostaining for ATP5H was also observed in the NFL compared to the vehicle-control. (RFI = 0.645 ± 0.1814, n = 6).
Similar to the pattern of expression observed for ATP5H, ET-1 treatment (24 h) showed decreased immunoreactivity for COX17 in NFL, GCL and IPL in rat retinas (Fig. 7c). A semi-quantitative analysis performed indicated a statistically significant decrease in COX17 expression in the NFL (RFI = 0.4914 ± 0.1738, n = 6, p <0.05) and GCs (RFI = 0.6483 ± 0.1518, n = 6, p <0.05) (Fig. 7c). A statistically significant difference was not found in the IPL (RFI = 0.6867 ± 0.1593, n = 6), however, the same trend of decrease as observed in ATP5H expression was still present. Taken together, the data suggest that there was a downregulation in the expression of the key mitochondrial proteins ATP5H and COX17, 24 hours following ET-1 treatment.
IOP elevation produced a modest decrease in immunoreactive levels of ATP5H
Retina sections were obtained from four retired breeder Brown Norway rats, which were IOP elevated in one eye, using the other eye as the contralateral control (Fig. 8a). Immunostaining was carried out using antibodies to ATP5H and COX17 and images were taken on a Zeiss LSM580 confocal microscope. A semi-quantitative analysis measuring the change in integrated intensity was also performed using ImageJ software. Projections from the confocal image z-stacks were analyzed for integrated density of ATP5H and COX17 immunostaining and compared between the IOP-elevated eyes and their corresponding contralateral eyes (Fig. 8b). It was found that IOP elevation produced a decreasing trend in immunostaining for ATP5H, compared to that of the corresponding contralateral eyes (Fig. 8d). There was no appreciable change in immunostaining for COX17 between IOP elevated eyes and contralateral eyes (Fig. 8c).
Discussion
ET-1 is a potent vasoactive peptide that has been shown to be elevated in both the aqueous humor and circulation of primary open angle glaucoma patients37,38,39,40. An elevation of both ET-1 and homocysteine was found in the plasma of primary open angle glaucoma patients, compared to normal tension glaucoma as well as control subjects, suggesting that elevated oxidative stress and ET-1 are involved in the early stages of glaucoma pathology41. Single nucleotide polymorphisms (SNPs) in both the ET-1 and ETA receptor genes that affected ET-1 plasma concentrations and higher ET-1 levels correlated with increased visual field damage42. A recent meta-analysis that included seven studies (212 cases, 164 controls) for normal tension glaucoma and six studies (160 cases, 174 controls) for the POAG analysis demonstrated increased circulating ET-1 levels in glaucoma patients, compared to control subjects43.
Several studies have pointed to ET-1’s key role in neurodegeneration, however, the precise mechanisms by which ET-1 produces neurodegenerative effects is still an active area of investigation. Previous studies from our laboratory have shown that intravitreal ET-1 injection in rats decreased anterograde axonal transport of vesicles associated with the mitochondrial subcomponent9. Axonal transport is an energy dependent process, mediated by motor proteins which use the energy of ATP hydrolysis to produce force and movement along the microtubules. The decrease in fast axonal transport could be due to a decline ATP production, disruption of the cytoskeleton or a combination of these factors. Most of the ATP needed for cellular metabolism/activities is generated by oxidative phosphorylation in the mitochondria. This raises the possibility that ET-1 could alter mitochondrial dynamics as well as oxidative metabolism in the optic nerve head during glaucoma. It has been shown that ETB receptor stimulation by administration of sarafotoxin (ETB receptor agonist) in rats increases the production of superoxide anions, in sympathetic neurons36. Increased oxidative stress by endothelin receptor stimulation could have damaging effects on the mitochondria, thereby providing a strong rationale to assess expression of mitochondrial genes and their contribution to neurodegeneration.
Dysregulation of gene expression contributes to neurodegeneration in glaucoma, by acting on several cellular pathways, ultimately enhancing pro-apoptotic effects while suppressing anti-apoptotic gene expression44,45. Numerous studies have been carried out using total cellular RNA from either the retina or cultured RGCs to evaluate changes in gene expression in various models of optic neuropathy46,47. While an assessment of changes in gene expression at the mRNA level is informative, it does not provide an indication of the corresponding changes in the protein levels. In the polysomal profiles, there is a drastic shift of distribution of mRNA and its association with multiple ribosomes, following metabolic insults before it is detected on the proteome level48. Hence, assessment of polysomal RNA will provide a reflection of the translatome, offering a glimpse into de novo protein synthesis which is an important manifestation of changes in gene expression at the protein level.
In the current study, using RNA-seq we analyzed the changes in gene expression of polysomal RNA which represent the pool of actively translated mRNAs. This provides a unique opportunity to identify novel genes whose expression are altered at the protein level during pathogenesis. Identification of novel genes that contribute to glaucomatous neurodegeneration will provide the basic information needed to develop new therapies aimed at neuroprotection in glaucoma. To our knowledge, this is the first study to use RNA-seq analysis to investigate the RGC translatome profile following ET-1 treatment. Although other researchers have conducted microarray analyses of gene expression under glaucomatous conditions, most of these studies focused on mRNA expression which may not be reflected at the level of the protein expression49,50,51.
The results of the RNA seq analysis were confirmed by qPCR analysis of some of the candidate genes. Some of the data derived from the qPCR data were different from that found in the RNA-seq analysis. For instance, while an increase in Atp5h and Cox17 was found in the RNA-seq analysis, a decrease in the expression of these genes were found by qPCR analysis. Some of these variations could be related to the template half-life of the individual transcripts. Less abundant transcripts in RNA samples with low concentrations could be rapidly degraded, particularly with longer handling times. The differences in the handling time between RNA-seq analysis and qPCR analysis could be one plausible reason for the discrepancies in the findings between RNA-seq and qPCR, as well as account for the big variations observed in RNA expression that contributed to the lack of statistical significance in many of the qPCR analyses. On the other hand, at the protein level, the data was very consistent and a decrease in levels of ATP5H was observed both in cultured RGCs treated with ET-1 and in retinas of rats intravitreally injected with ET-1.
A notable finding from the RNA-seq analysis was that 23 out of 156 DEG in the RNA-seq analysis were identified with known or predicted mitochondrial function. Ingenuity Pathway Analysis revealed oxidative phosphorylation as the top significantly enriched canonical pathway, which points to the involvement of mitochondrial bioenergetics in RGCs with ET-1 treatment. To specifically assess the nuclear encoded genes in the mitochondrial respiratory chain complex and oxidative phosphorylation, the mRNA and protein expression of Cox17, and Atp5h which are critical components of Complex IV and Complex V, respectively, were determined. COX17 is the terminal enzyme of the mitochondrial respiratory chain while mitochondrial ATP synthase catalyzes ATP synthesis, utilizing an electrochemical gradient of protons across the inner membrane during oxidative phosphorylation28,52. IOP elevation has previously shown to cause a reduction of Cytochrome c oxidase IV subunit 1 activity by damaging the mitochondria in response to mitochondrial fission53. There is a decline of mitochondrial function during early stages of POAG due to reduction in ATP synthesis, particularly making RGCs susceptible for oxidative phosphorylation impairment54. Mitochondrial mutation and mitochondrial genetic variation can also act as a potential risk factor in POAG pathogenesis55,56. Interestingly, a recent study showing association of mitochondrial proteins with the pathogenesis of POAG identified a number of key genes including Atp5h, Cox17, Ndufs7, Ndufb9, Uqcrc1 and Ak3 involved in oxidative phosphorylation in human POAG that were identical to some of the genes observed in our findings57. In the current study, we identified that there was a decline in two important mitochondrial proteins ATP5H and COX17 both in primary RGCs treated with ET-1 as well as in rats intravitreally injected with ET-1. In addition, a decreasing trend in immunostaining for ATP5H was found following IOP elevation in Brown Norway rats.
The IOP elevation model is a chronic model and the data presented is indicative of IOP-mediated changes that occur over longer period of time, compared to the acute effects that occur following a bolus intravitreal administration of 2 nmole of ET-1. Nevertheless, a trend towards a decrease in ATP5H was observed following IOP elevation for 2 weeks in Brown Norway rats. This suggests that a decline in mitochondrial function could compromise bioenergetics mediated by a decline in ATP levels. The Seahorse assay which exploits the coupling of ATP synthesis to oxygen consumption that occurs during oxidative phosphorylation, provided a clear evidence for the ability of ET-1 to significantly decline mitochondrial ATP production in RGCs. The use of oligomycin to inhibit mitochondrial ATP synthase and assess the decline in the oxygen consumption rate, provides a glimpse of ATP-linked respiration, which differed significantly between the control and the ET-1 treated RGCs. The data points to the ability of ET-1 to functionally impact ATP synthesis in RGCs. The findings have important implications for RGC death occurring through energy depletion, which could occur through opening of permeability transition pore, leading to mitochondrial depolarization and mitochondrial damage. It is well known that reactive oxygen species are a major contributor to mitochondrial outer membrane permeabilization. As mentioned earlier, ETB receptor activation has been shown to result in elevation of oxidative stress as reported in dorsal root ganglion neurons36. It is plausible that similar mechanisms could operate to generate oxidative stress through the actions of ET-1 on its receptors in RGCs.
Taken together, our findings have important implications for the role of mitochondria in the pathogenesis of glaucoma. Further work is required to understand the role of the genes identified in the RNA-seq analysis and elucidate the signaling pathways which contribute to the expression of these genes. The current study could provide an insight to advance our understanding of mitochondrial dynamics contributing to RGC loss in POAG.
Methods
Isolation of primary retinal ganglion cells
Primary rat RGCs isolation was performed according to the method of Barres et al.58. All protocols and procedures were in accordance with the policies of the Association for Research in Vision and Ophthalmology (ARVO) for use of animals in research and approved by the Institutional Animal Care and Use Committee (IACUC) (animal protocol #IACUC-2017-0024) at University of North Texas Health Science Center at Fort Worth, TX, USA. Post-natal day 4–7 rat pups were euthanized, and the isolated retinas were treated with 4.5 units/ml of papain solution to dissociate the cells. The cell suspension was subsequently incubated for 10 minutes with a rabbit anti-macrophage antibody and transferred to a petri dish coated with a goat anti-rabbit IgG (H + L chain) antibody for 35 minutes. Cells that were not attached to the coated anti-rabbit IgG were transferred to a petri-dish coated with anti-Thy1.1 antibody for 1 hour with intermittent shaking. Cells were then detached by trypsin treatment and collected by centrifugation. Cells were seeded either directly in a 6-well plate or on coverslips coated with poly-D lysine and mouse-laminin and cultured in serum-free Dulbecco’s modified Eagle’s medium containing brain-derived neurotrophic factor (BDNF) (50 ng/ml; Peprotech, Rocky Hill, NJ, USA), ciliary neurotrophic factor (10 ng/ml; Peprotech), and forskolin (5 ng/ml; Sigma-Aldrich Corp.). Cells were incubated at 37 °C in a humidified atmosphere of 10% CO2 and 90% air. One-third volume of the culture medium was changed every two days. The purity of the culture was determined by counting RBPMS and Brn3a positive cells per 100 cells from different regions in an image. The purity of the culture obtained is routinely between 90%-95% in most batches of isolated RGCs.
ET-1 treatment
ET-1 peptide was purchased from Bachem (Torrance, CA) and dissolved to a stock concentration of 500 μM and working stocks of 100 μM. The seeding density of the primary RGCs was approximately 350,000 / well in a 6-well plate. After 1 week in culture to promote neurite outgrowth, the isolated RGCs were treated with 100 nM ET-1 (final concentration) for 24 hours while untreated RGCs served as a control.
Polysomal RNA isolation
Following ET-1 treatment, a brief incubation with cycloheximide (10 μg/ml) to inhibit protein synthesis was done, total polysomes were isolated by magnesium precipitation by a modification of the method of Palmiter et al.59. Cells were washed, scraped in PBS and the cell pellet was suspended in 150 μl of homogenization buffer (25 mM Tris-HCl pH 7.5, 25 mM NaCl, 5 mM MgCl2, and 1% Triton X-100). After an incubation on ice for 5 min, the homogenate was centrifuged at 10,000xg in a microfuge for 10 min at 4 °C to sediment the nuclei and mitochondria. The pellet was discarded and the supernatant was treated with an equal volume of polysome precipitation buffer (25 mM Tris-HCl pH 7.5, 25 mM NaCl, 250 mM MgCl2 and 1 M sucrose). The content of the tubes were mixed by inverting and incubated on ice for one hour with intermittent mixing. The polysomal fraction was obtained by centrifugation at 40,000 × g for 1 hr in a Beckman ultracentrifuge at 4 °C using the TLA-55 fixed angle rotor. The supernatant was discarded and the polysomal pellet was extracted using the Trizol reagent to isolate polysomal RNA according to the manufacturer’s instructions.
Library preparation and RNA-seq analysis
Libraries for RNA-Seq were prepared in the UCLA core facility with Clonetech SMARTer Stranded Total RNA-Seq (Pico) Kit. The workflow consisted of first-strand synthesis, template switching, adaptor ligation, cleavage of ribosomal cDNA and PCR amplification. Different adaptors were used for multiplexing samples in one lane. Sequencing was performed on Illumina Hiseq. 3000 for a single read 50 bp run. Data quality check was done on Illumina SAV. The reads were first mapped to the latest UCSC transcript set using Bowtie2 version 2.1.060 and the gene expression level was estimated using RSEM v1.2.1561. TMM (trimmed mean of M-values) was used to normalize the gene expression. Differentially expressed genes were identified using the edgeR program62. Genes showing altered expression with p<0.05 and more than 1.5 fold changes were considered differentially expressed, since a fold increase of less than 1.5-fold would typically not produce appreciable phenotypic changes. Differentially expressed genes (DEGs) based on grouped comparison of untreated and treated cells were analyzed for biological and functional relevance using the Rattus norvegicus knowledgebase in PANTHER v.13.1; GO classifications based on cellular component, molecular function, biological process, and reactome pathways were analyzed. MitoMiner v.4.0 was used to identify the DEG subset that has a known or predicted mitochondrial function. A STRING network was constructed using the 156 genes to identify known and potential protein-protein interactions within the set of differentially expressed genes (DEGs). DEGs were ranked using the sign (fold change)*(1/p-value) method to identify genes that are both highly dysregulated (up/down fold change magnitude) and with a low probability of false discovery (low p-value).
Quantitative real-time PCR
Primary RGCs were isolated from post-natal day 4–6 rat pups (as described above) and total cellular RNA was extracted using the Trizol reagent (as described above). Complementary DNA was prepared from total cellular RNA using iScript Reverse Transcription Kit (BioRad, Hercules, CA), and Real-time PCR (SsoAdvanced™ SYBR® Green Supermix) was performed using cDNA as template to detect gene expression of some DEGs including Mitochondrial Ribosomal Protein L30 (Mrpl30), ATP Synthase, H+ Transporting, Mitochondrial F0 Complex, Subunit D (Atp5h), cytochrome C oxidase copper chaperone (Cox17), Glucose-6-phosphate dehydrogenase (G6pd), Forkhead box protein (Foxo1), Mitogen-Activated Protein Kinase Kinase Kinase 11 (Map3k11) and Modulator Of Apoptosis 1 (Moap1). PCR conditions were as follows: 95 °C for 3 minutes followed by 40 three-step cycles of 95 °C for 10 seconds, 60 °C for 30 seconds and 72 °C for 30 seconds. PCR primers were validated and the authenticity of the amplicons were confirmed by their melting temperatures. RNA samples from 3–4 biological replicates were used for real-time PCR analysis. Data was expressed as mean ± SEM and normalized to cyclophilin-1. The results were presented as relative fold change compared to untreated control. The statistical significance was calculated by Student’s t-test, and p values <0.05 were considered significant.
Immunocytochemistry
Purified RGCs were cultured for 7 days to promote neurite outgrowth. RGCs were then treated with 100 nM ET-1 for 24 h and fixed with 4% paraformaldehyde for 10 min. After permeabilization with 0.1% Triton X-100, the cells were blocked with 5% bovine serum albumin and 5% normal donkey serum at room temperature for 1 hour to prevent non-specific binding with the secondary antibody. Blocking was followed by overnight incubation at 4 °C with primary antibodies against mitochondrial proteins ATP5H and Cox17. Primary antibodies used were rabbit anti-COX 17 (1:100, #ab69611, Abcam), mouse anti-ATP5H (1:200, #ab110275, Abcam) and goat anti-Brn3a (1:200, #sc31984, Santa Cruz Biotechnology). Following primary antibody incubations, cells were washed with PBS and then incubated with the corresponding secondary antibody at room temperature (RT) in the dark. Secondary antibodies used were donkey anti-rabbit Alexa 647 (1:1000, Invitrogen), donkey anti-mouse Alexa 546 (1:1000, Invitrogen) and donkey anti-goat 488 (1:1000, Invitrogen). Immunostaining carried out with the exclusion of the primary antibody and using only the secondary antibody conjugated to the fluorophore served as the negative control, which was indicative of non-specific staining by the secondary antibody. Following the secondary antibody incubation, cells were washed with PBS and incubated with with 4′ 6 Diamidino-phenylindole dichloride (DAPI) to stain cell nuclei. The coverslips were mounted on slides with antifade medium (FluorSave; Calbiochem, La Jolla, CA, USA). Images were captured for each field of view using the Zeiss 510 meta confocal microscope.
Intravitreal injection of ET-1
ET-1 peptide was purchased from Bachem (Torrance, CA) and dissolved to a stock concentration of 500 μM. One group of male Brown Norway rats (n=6) were intravitreally injected in one eye with 4 μl of 500 μM ET-1 in 0.25% glacial acetic acid (adjusted to pH 7.0) while the contralateral eye served as control. The other group (n=6) received vehicle (0.25% glacial acetic acid, neutralized to pH 7.0) intravitreally in one eye, while the contralateral eye was control. Rats were euthanized 24 hours after injection by intraperitoneal pentobarbital injection (120 mg/kg body wt) followed with intracardiac pentobartital injection.
Morrison model of IOP elevation in rats
IOP elevation was carried out in one eye of four retired breeder female Brown Norway rats by injection of 100 μl of hypertonic saline through episcleral veins63. Briefly, rats were anesthetized by inhalation of isoflurane and ensured by lack of a toe pinch reflex. Following a conjunctival incision, hypertonic saline (1.8 M NaCl) was injected into an episcleral vein using a glass needle (TIP01TW1F, WPI). The companion eye was used as the contralateral control eye. A flow rate of 309 µL/min was maintained during the injection for 10 to 20 seconds, thereby injecting 50 to 100 μl of hypertonic saline into the episcleral veins. IOP was measured two times a week using a hand-held TonoLab tonometer (iCare, Finland). Six IOP readings were averaged from each IOP measurement and ten IOP measurements were obtained for each eye. IOP plots were generated from IOP values obtained from the surgically treated eye and contralateral control eye. The IOP exposure in each rat was computed by the integral product of the extent of IOP elevation and the number of days of IOP elevation (expressed as mm Hg-days).
Immunohistochemistry
Eyes were enucleated and a small incision (4 mm) was made posterior to the limbus and the eye was fixed in 4% paraformaldehyde (PFA) in phosphate buffered saline (PBS) for 30 min. Subsequently, an incision was made along the entire circumference of the limbus to remove the entire anterior segment. The optic cups were fixed in 4% PFA for an additional 3 hours. After fixation, the tissue was rinsed with 70% ethanol and later embedded in paraffin. Seven-micron sagittal retinal sections through the optic nerve head were obtained, de-paraffinized in xylene, and rehydrated with a graded series of ethanol. Following permeabilization with 0.1% Triton X-100 and blocking with 5% normal donkey serum containing 5% BSA in PBS, retinal sections were incubated with the appropriate primary antibodies overnight. Primary antibodies used were rabbit anti-Cox 17 (1:100, #ab69611, Abcam), mouse anti-ATP5H (1:200, #ab110275, Abcam)and goat anti-Brn3a (1:200, #sc31984, Santa Cruz Biotechnology). Secondary antibody incubation for 1 hr was carried out at room temperature. Secondary antibodies used were either donkey anti-rabbit Alexa 546 (1:1000, Invitrogen), donkey anti-mouse Alexa 647 (1:1000, Invitrogen) or donkey anti-goat 488 (1:1000, Invitrogen). Retinal sections incubated with no primary antibody served as the blank to assess non-specific staining by the secondary antibodies.
Seahorse XF real-time ATP rate assay
Agilent Seahorse XFe96 analyzer was used to measure multiple parameters including oxygen consumption rate (OCR), total ATP production by quantifying ATP production rate from both glycolytic and mitochondrial pathways and ATP rate index. Briefly, primary RGCs were isolated from post-natal day 4–6 rat pups. The RGCs were seeded RGC medium in a poly-D-lysine coated seahorse plate and maintained in culture for one week to allow neurite outgrowth. The RGCs were either untreated or treated with ET-1 for 4 h or 24 h in trophic factor-free medium. A day prior to the experiment seahorse sensor cartridge was hydrated with water and incubated in a non CO2 incubator at 37 °C. Two hours before the experiment, the water was replaced with a seahorse calibrant. On the day of the experiment, seahorse XF DMEM (pH=7.4) was supplemented with 1 mM pyruvate, 2 mM glutamine and 10 mM glucose (Agilent, USA). The RGC medium was then replaced with 180 µl of Seahorse XF base medium and the cells incubated in a non CO2 incubator at 37 °C for 1 hour. During this incubation period, oligomycin and Rotenone/antimycin (Agilent, USA) were prepared in the seahorse medium to achieve final concentrations of 1.5 μM and 0.5 μM respectively when injected. From these stock solutions, 20 μl and 22 μl of Oligomycin and rotenone/antimycin respectively were then loaded into the drug delivery ports A and B of the hydrated sensor cartridge and loaded into the seahorse XF analyzer to calibrate for 30 minutes. The calibration plate was replaced with the cell culture plate and oxygen consumption (OCR) and extracellular acidification rate (ECAR) was monitored following the sequential injection of Oligomycin and rotenone/antimycin with each cycle set as 3 min mix, 2 min delay and measure for 3 minutes. According to the manufacturer’s instructions (Agilent Technologies, Santa Clara, CA), the following equation was used to calculated mitochondrial ATP production:
The ATP production rate in ET-1 treated RGCs was compared to that of untreated RGCs for two time points of ET-1 treatment (4 h and 24 h). Data was normalized to the cell number of each well by Calcien AM assay.
Statistical analysis
GraphPad Prism 7 (La Jolla, CA, USA) was used to perform statistical analysis. Student’s t-test was used while comparing two experimental groups (control versus treated). Statistical significance of the experimental data was described as *p < 0.05; **p < 0.01. Data are presented as mean ±SEM.
References
Quigley, H. A. & Broman, A. T. The number of people with glaucoma worldwide in 2010 and 2020. Br. J. Ophthalmol. 90, 262–267, https://doi.org/10.1136/bjo.2005.081224 (2006).
Tham, Y. C. et al. Global prevalence of glaucoma and projections of glaucoma burden through 2040: a systematic review and meta-analysis. Ophthalmol. 121, 2081–2090, https://doi.org/10.1016/j.ophtha.2014.05.013 (2014).
Quigley, H. A. Neuronal death in glaucoma. Prog. retinal eye Res. 18, 39–57 (1999).
Quigley, H. A. 21st century glaucoma care. Eye 33, 254–260, https://doi.org/10.1038/s41433-018-0227-8 (2019).
The Advanced Glaucoma Intervention Study (AGIS): 7. The relationship between control of intraocular pressure and visual field deterioration.The AGIS Investigators. Am. J. Ophthalmol. 130, 429–440 (2000).
Vass, C. et al. Medical interventions for primary open angle glaucoma and ocular hypertension. The Cochrane database of systematic reviews, Cd003167, https://doi.org/10.1002/14651858.CD003167.pub3 (2007).
Leske, M. C. et al. Factors for glaucoma progression and the effect of treatment: the early manifest glaucoma trial. Arch. Ophthalmol. 121, 48–56 (2003).
Prasanna, G., Krishnamoorthy, R., Clark, A. F., Wordinger, R. J. & Yorio, T. Human Optic Nerve Head Astrocytes as a Target for Endothelin-1. Investig. Ophthalmol. Vis. Sci. 43, 2704–2713 (2002).
Stokely, M. E., Brady, S. T. & Yorio, T. Effects of endothelin-1 on components of anterograde axonal transport in optic nerve. Invest. Ophthalmol. Vis. Sci. 43, 3223–3230 (2002).
Narayan, S., Prasanna, G., Krishnamoorthy, R. R., Zhang, X. & Yorio, T. Endothelin-1 Synthesis and Secretion in Human Retinal Pigment Epithelial Cells (ARPE-19): Differential Regulation by Cholinergics and TNF-α. Investig. Ophthalmol. Vis. Sci. 44, 4885–4894, https://doi.org/10.1167/iovs.03-0387 (2003).
Krishnamoorthy, R. R. et al. Role of the ETB receptor in retinal ganglion cell death in glaucomaThis article is one of a selection of papers published in the special issue (part 1 of 2) on Forefronts in Endothelin. Can. J. Physiol. Pharmacology 86, 380–393, https://doi.org/10.1139/Y08-040 (2008).
Yorio, T., Krishnamoorthy, R. & Prasanna, G. Endothelin: is it a contributor to glaucoma pathophysiology? J. glaucoma 11, 259–270 (2002).
Howell, G. R. et al. Molecular clustering identifies complement and endothelin induction as early events in a mouse model of glaucoma. J. Clin. investigation 121, 1429–1444, https://doi.org/10.1172/jci44646 (2011).
Howell, G. R. et al. Combinatorial targeting of early pathways profoundly inhibits neurodegeneration in a mouse model of glaucoma. Neurobiol. Dis. 71, 44–52, https://doi.org/10.1016/j.nbd.2014.07.016 (2014).
McGrady, N. R. et al. Upregulation of the endothelin A (ETA) receptor and its association with neurodegeneration in a rodent model of glaucoma. BMC Neurosci. 18, 27, https://doi.org/10.1186/s12868-017-0346-3 (2017).
Minton, A. Z. et al. Endothelin B receptors contribute to retinal ganglion cell loss in a rat model of glaucoma. PLoS one 7, e43199, https://doi.org/10.1371/journal.pone.0043199 (2012).
Chauhan, B. C. Endothelin and its potential role in glaucoma. Can. J. Ophthalmol. 43, 356–360, https://doi.org/10.3129/i08-060 (2008).
Rosenthal, R. & Fromm, M. Endothelin antagonism as an active principle for glaucoma therapy. Br. J. pharmacology 162, 806–816, https://doi.org/10.1111/j.1476-5381.2010.01103.x (2011).
Kallberg, M. E. et al. Endothelin 1 levels in the aqueous humor of dogs with glaucoma. J. glaucoma 11, 105–109 (2002).
Prasanna, G. et al. Effect of elevated intraocular pressure on endothelin-1 in a rat model of glaucoma. Pharmacol. Res. 51, 41–50, https://doi.org/10.1016/j.phrs.2004.04.006 (2005).
Thanos, S. & Naskar, R. Correlation between retinal ganglion cell death and chronically developing inherited glaucoma in a new rat mutant. Exp. Eye Res. 79, 119–129, https://doi.org/10.1016/j.exer.2004.02.005 (2004).
Chauhan, B. C. et al. Model of Endothelin-1–Induced Chronic Optic Neuropathy in Rat. Investig. Ophthalmol. Vis. Sci. 45, 144–152, https://doi.org/10.1167/iovs.03-0687 (2004).
Lau, J., Dang, M., Hockmann, K. & Ball, A. K. Effects of acute delivery of endothelin-1 on retinal ganglion cell loss in the rat. Exp. Eye Res. 82, 132–145, https://doi.org/10.1016/j.exer.2005.06.002 (2006).
Cuenca, N. et al. Cellular responses following retinal injuries and therapeutic approaches for neurodegenerative diseases. Prog. retinal eye Res. 43, 17–75, https://doi.org/10.1016/j.preteyeres.2014.07.001 (2014).
Pinazo-Duran, M. D., Zanon-Moreno, V., Garcia-Medina, J. J. & Gallego-Pinazo, R. Evaluation of presumptive biomarkers of oxidative stress, immune response and apoptosis in primary open-angle glaucoma. Curr. Opin. pharmacology 13, 98–107, https://doi.org/10.1016/j.coph.2012.10.007 (2013).
Tezel, G. TNF-alpha signaling in glaucomatous neurodegeneration. Prog Brain Res. 173, 409–421, https://doi.org/10.1016/s0079-6123(08)01128-x (2008).
Zanon-Moreno, V., Garcia-Medina, J. J., Zanon-Viguer, V., Moreno-Nadal, M. A. & Pinazo-Duran, M. D. Smoking, an additional risk factor in elder women with primary open-angle glaucoma. Mol. Vis. 15, 2953–2959 (2009).
Carelli, V., Ross-Cisneros, F. N. & Sadun, A. A. Mitochondrial dysfunction as a cause of optic neuropathies. Prog. retinal eye Res. 23, 53–89, https://doi.org/10.1016/j.preteyeres.2003.10.003 (2004).
Osborne, N. N. Mitochondria: Their role in ganglion cell death and survival in primary open angle glaucoma. Exp. Eye Res. 90, 750–757, https://doi.org/10.1016/j.exer.2010.03.008 (2010).
Osborne, N. N. & del Olmo-Aguado, S. Maintenance of retinal ganglion cell mitochondrial functions as a neuroprotective strategy in glaucoma. Curr. Opin. pharmacology 13, 16–22, https://doi.org/10.1016/j.coph.2012.09.002 (2013).
Kong, G. Y., Van Bergen, N. J., Trounce, I. A. & Crowston, J. G. Mitochondrial dysfunction and glaucoma. J. glaucoma 18, 93–100, https://doi.org/10.1097/IJG.0b013e318181284f (2009).
Yu-Wai-Man, P. Mitochondrial Dysfunction in Glaucoma—Closing the LoopMitochondrial Dysfunction in Glaucoma. Investig. Ophthalmol. Vis. Sci. 53, 2438–2438, https://doi.org/10.1167/iovs.12-9815 (2012).
Ji, T. et al. Does perturbation in the mitochondrial protein folding pave the way for neurodegeneration diseases? Ageing Res. Rev. 57, 100997, https://doi.org/10.1016/j.arr.2019.100997 (2019).
Boyman, L., Karbowski, M. & Lederer, W. J. Regulation of Mitochondrial ATP Production: Ca(2+) Signaling and Quality Control. Trends Mol Med, https://doi.org/10.1016/j.molmed.2019.10.007 (2019).
Williams, P. A., Harder, J. M. & John, S. W. M. Glaucoma as a Metabolic Optic Neuropathy: Making the Case for Nicotinamide Treatment in Glaucoma. J. glaucoma 26, 1161–1168, https://doi.org/10.1097/ijg.0000000000000767 (2017).
Lau, Y. E., Galligan, J. J., Kreulen, D. L. & Fink, G. D. Activation of ETB receptors increases superoxide levels in sympathetic ganglia in vivo. Am. J. Physiol. Regulatory, Integr. Comp. physiology 290, R90–95, https://doi.org/10.1152/ajpregu.00505.2005 (2006).
Iwabe, S., Lamas, M., Vasquez Pelaez, C. G. & Carrasco, F. G. Aqueous humor endothelin-1 (Et-1), vascular endothelial growth factor (VEGF) and cyclooxygenase-2 (COX-2) levels in Mexican glaucomatous patients. Curr. Eye Res. 35, 287–294, https://doi.org/10.3109/02713680903545315 (2010).
Choritz, L., Machert, M. & Thieme, H. Correlation of Endothelin-1 Concentration in Aqueous Humor with Intraocular Pressure in Primary Open Angle and Pseudoexfoliation Glaucoma. Investig. Ophthalmol. Vis. Sci. 53, 7336–7342, https://doi.org/10.1167/iovs.12-10216 (2012).
Cellini, M. et al. Endothelin-1 plasma levels and vascular endothelial dysfunction in primary open angle glaucoma. Life Sci. 91, 699–702, https://doi.org/10.1016/j.lfs.2012.02.013 (2012).
Chen, H. Y., Chang, Y. C., Chen, W. C. & Lane, H. Y. Association between plasma endothelin-1 and severity of different types of glaucoma. J. glaucoma 22, 117–122, https://doi.org/10.1097/IJG.0b013e31822e8c65 (2013).
Lopez-Riquelme, N. et al. Endothelin-1 levels and biomarkers of oxidative stress in glaucoma patients. Int. Ophthalmol. 35, 527–532, https://doi.org/10.1007/s10792-014-9979-8 (2015).
Kosior-Jarecka, E. et al. Plasma endothelin-1 and single nucleotide polymorphisms of endothelin-1 and endothelin type A receptor genes as risk factors for normal tension glaucoma. Mol. Vis. 22, 1256–1266 (2016).
Li, S., Zhang, A., Cao, W. & Sun, X. Elevated Plasma Endothelin-1 Levels in Normal Tension Glaucoma and Primary Open-Angle Glaucoma: A Meta-Analysis. J. Ophthalmol. 2016, 2678017, https://doi.org/10.1155/2016/2678017 (2016).
Nickells, R. W., Howell, G. R., Soto, I. & John, S. W. Under pressure: cellular and molecular responses during glaucoma, a common neurodegeneration with axonopathy. Annu. Rev. Neurosci. 35, 153–179, https://doi.org/10.1146/annurev.neuro.051508.135728 (2012).
Nickells, R. W. & Pelzel, H. R. Tools and resources for analyzing gene expression changes in glaucomatous neurodegeneration. Exp. Eye Res. 141, 99–110, https://doi.org/10.1016/j.exer.2015.05.009 (2015).
Ahmed, F. et al. Microarray Analysis of Changes in mRNA Levels in the Rat Retina after Experimental Elevation of Intraocular Pressure. Investig. Ophthalmol. Vis. Sci. 45, 1247–1258, https://doi.org/10.1167/iovs.03-1123 (2004).
White, A. J., Heller, J. P., Leung, J., Tassoni, A. & Martin, K. R. Retinal ganglion cell neuroprotection by an angiotensin II blocker in an ex vivo retinal explant model. J. Renin-Angiotensin-Aldosterone Syst. 16, 1193–1201, https://doi.org/10.1177/1470320314566018 (2015).
Panda, A. C., Martindale, J. L. & Gorospe, M. Polysome Fractionation to Analyze mRNA Distribution Profiles. Bio-protocol 7, https://doi.org/10.21769/BioProtoc.2126 (2017).
He, S., Park, Y. H., Yorio, T. & Krishnamoorthy, R. R. Endothelin-Mediated Changes in Gene Expression in Isolated Purified Rat Retinal Ganglion CellsEndothelin-Mediated Changes in Gene Expression. Investig. Ophthalmol. Vis. Sci. 56, 6144–6161, https://doi.org/10.1167/iovs.15-16569 (2015).
Yang, Z. et al. Changes in gene expression in experimental glaucoma and optic nerve transection: the equilibrium between protective and detrimental mechanisms. Invest. Ophthalmol. Vis. Sci. 48, 5539–5548, https://doi.org/10.1167/iovs.07-0542 (2007).
Johnson, E. C. et al. Cell proliferation and interleukin-6-type cytokine signaling are implicated by gene expression responses in early optic nerve head injury in rat glaucoma. Invest. Ophthalmol. Vis. Sci. 52, 504–518, https://doi.org/10.1167/iovs.10-5317 (2011).
Srinivasan, S. & Avadhani, N. G. Cytochrome c oxidase dysfunction in oxidative stress. Free. Radic. Biol. Med. 53, 1252–1263, https://doi.org/10.1016/j.freeradbiomed.2012.07.021 (2012).
Ju, W. K. et al. Intraocular pressure elevation induces mitochondrial fission and triggers OPA1 release in glaucomatous optic nerve. Invest. Ophthalmol. Vis. Sci. 49, 4903–4911, https://doi.org/10.1167/iovs.07-1661 (2008).
Lee, S. et al. Impaired complex-I-linked respiration and ATP synthesis in primary open-angle glaucoma patient lymphoblasts. Invest. Ophthalmol. Vis. Sci. 53, 2431–2437, https://doi.org/10.1167/iovs.12-9596 (2012).
Abu-Amero, K. K., Morales, J. & Bosley, T. M. Mitochondrial Abnormalities in Patients with Primary Open-Angle Glaucoma. Investig. Ophthalmol. Vis. Sci. 47, 2533–2541, https://doi.org/10.1167/iovs.05-1639 (2006).
Collins, D. W. et al. Association of primary open-angle glaucoma with mitochondrial variants and haplogroups common in African Americans. Mol. Vis. 22, 454–471 (2016).
Khawaja, A. P. et al. Assessing the Association of Mitochondrial Genetic Variation With Primary Open-Angle Glaucoma Using Gene-Set AnalysesMitochondrial Genetic Variation and POAG. Investig. Ophthalmol. Vis. Sci. 57, 5046–5052, https://doi.org/10.1167/iovs.16-20017 (2016).
Barres, B. A., Silverstein, B. E., Corey, D. P. & Chun, L. L. Immunological, morphological, and electrophysiological variation among retinal ganglion cells purified by panning. Neuron 1, 791–803 (1988).
Palmiter, R. D. Magnesium precipitation of ribonucleoprotein complexes. Exped. Tech. isolation undergraded polysomes messenger ribonucleic acid. Biochem. 13, 3606–3615 (1974).
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. methods 9, 357–359, https://doi.org/10.1038/nmeth.1923 (2012).
Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 12, 323, https://doi.org/10.1186/1471-2105-12-323 (2011).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinforma. 26, 139–140, https://doi.org/10.1093/bioinformatics/btp616 (2010).
Morrison, J. C. et al. A rat model of chronic pressure-induced optic nerve damage. Exp. Eye Res. 64, 85–96, https://doi.org/10.1006/exer.1996.0184 (1997).
Acknowledgements
This work was supported by an extramural grant (EY028179) to RRK from the National Eye Institute.
Author information
Authors and Affiliations
Contributions
R.M.C. designed, performed experiments, analyzed data and wrote the manuscript. D.L.S. and S.H. designed and performed experiments and interpreted data. B.K. performed some experiments. N.P. carried out bioinformatics analysis and interpretation of the data. S.Y. and J.P. aided in the design, performance and interpretation of the Seahorse assay. R.K. was involved in conceptualization, data analysis and writing the manuscript. The authors thank Dr. Denise Inman for useful discussions in interpretation of the seahorse data.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chaphalkar, R.M., Stankowska, D.L., He, S. et al. Endothelin-1 Mediated Decrease in Mitochondrial Gene Expression and Bioenergetics Contribute to Neurodegeneration of Retinal Ganglion Cells. Sci Rep 10, 3571 (2020). https://doi.org/10.1038/s41598-020-60558-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-60558-6
This article is cited by
-
The role of oxidative stress in the pathogenesis of ocular diseases: an overview
Molecular Biology Reports (2024)
-
Flavinated SDHA underlies the change in intrinsic optical properties of oral cancers
Communications Biology (2023)
-
The role of the adaptive immune system and T cell dysfunction in neurodegenerative diseases
Journal of Neuroinflammation (2022)
-
Peptains block retinal ganglion cell death in animal models of ocular hypertension: implications for neuroprotection in glaucoma
Cell Death & Disease (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.