- Micro report
- Open Access
- Published:
Loss of park7 activity has differential effects on expression of iron responsive element (IRE) gene sets in the brain transcriptome in a zebrafish model of Parkinson’s disease
Molecular Brain volume 14, Article number: 83 (2021)
Abstract
Mutation of the gene PARK7 (DJ1) causes monogenic autosomal recessive Parkinson’s disease (PD) in humans. Subsequent alterations of PARK7 protein function lead to mitochondrial dysfunction, a major element in PD pathology. Homozygous mutants for the PARK7-orthologous genes in zebrafish, park7, show changes to gene expression in the oxidative phosphorylation pathway, supporting that disruption of energy production is a key feature of neurodegeneration in PD. Iron is critical for normal mitochondrial function, and we have previously used bioinformatic analysis of IRE-bearing transcripts in brain transcriptomes to find evidence supporting the existence of iron dyshomeostasis in Alzheimer’s disease. Here, we analysed IRE-bearing transcripts in the transcriptome data from homozygous park7−/− mutant zebrafish brains. We found that the set of genes with “high quality” IREs in their 5′ untranslated regions (UTRs, the HQ5′IRE gene set) was significantly altered in these 4-month-old park7−/− brains. However, sets of genes with IREs in their 3′ UTRs appeared unaffected. The effects on HQ5′IRE genes are possibly driven by iron dyshomeostasis and/or oxidative stress, but illuminate the existence of currently unknown mechanisms with differential overall effects on 5′ and 3′ IREs.
Introduction
Parkinson’s disease (PD) is the second most common neurodegenerative disease, affecting approximately 1% of the population aged over 60 years. Most cases of PD are idiopathic, but a clear genetic link has been identified in approximately 5–10% of PD cases [1]. One gene, PARK7, implicated in autosomal recessive early-onset PD, encodes Parkinson disease protein 7 (PARK7). PARK7 protein has been suggested to act as a GSH-independent glyoxalase for detoxification of methylglyoxal and as a protein glycase responsible for restoring the function of proteins damaged by oxidative stress. However, these activities are disputed in PD (reviewed and analysed in [2]). PARK7 also has critical roles in maintaining mitochondrial function, sensing and responding to reactive oxygen species (ROS), and ultimately, acts in neuroprotection (reviewed in [3]).
PD is characterized by the specific depletion of substantia nigra pars compacta dopaminergic (SNc DA) neurons. These neurons make large numbers of synapses in the basal ganglia. Consequently, their high energy demands may make SNc DA neurons sensitive to energy deficiency [4]. Many factors can affect energy production by the process of oxidative phosphorylation. In particular, ferrous iron (Fe2+) is incorporated into Fe-S clusters central to the function of the electron transport chain (ETC) in oxidative phosphorylation [5]. ETC dysfunction causes oxidative stress that may lead to the (primarily) cytosolic PARK7 protein translocating into mitochondria to regulate the effects of reactive oxygen species (ROS) [5]. This process is possibly altered in individuals with mutation of PARK7, leading to ETC dysfunction and dyshomeostasis of iron (via effects on Iron Regulatory Proteins, IRP1 & IRP2) and damaging dopaminergic neurons.
IRP1 and IRP2 bind IREs in the mRNAs of genes involved in iron homeostasis to regulate their translation and stability (reviewed in [6]). IRPs are regulated both by cellular ferrous iron status and oxidative stress [6]. Previously, we defined sets of genes bearing IREs in either the 5′ or 3′UTRs of their transcripts (at lower or higher similarity to an IRE consensus sequence) in humans, mice and zebrafish [7]. Using these, we found evidence supporting iron dyshomeostasis in Alzheimer’s disease (AD) brains, and in animal models of AD [7].
Orthologues of PD genes have previously been identified and manipulated in zebrafish. For example, Hughes et al. [8] developed a novel zebrafish model to examine the function of PARK7. Zebrafish which lack the PARK7-orthologous gene: park7 (park7−/−) display a movement phenotype at three months of age, and show changes to gene expression in their brain transcriptome suggestive of disrupted mitochondrial metabolism (i.e., upregulation of genes associated with oxidative phosphorylation) at four months of age [8]. Therefore, we hypothesized that oxidative stress and/or iron dyshomeostasis in park7−/− zebrafish brains would alter binding of IRPs to transcripts containing IREs, thereby altering transcript stability. To explore this, we reanalysed the zebrafish brain transcriptome data from Hughes et al. to test for changes in the representation of IRE-containing gene sets in park7−/− zebrafish brains. We found that the HQ5′IRE gene set is significantly altered in 4-month-old park7−/− brains, while the HQ3′IRE gene set is not.
Results
We previously defined sets of zebrafish genes according to whether these possess IRE-like motifs in either the 5′ or 3′ UTRs of their mRNAs, and whether their IREs match a canonical (high quality, HQ) or non-canonical IRE sequence (all) [7]. We found that only transcripts of the HQ5’IRE gene set show statistically significant changes to gene expression as a group (Fig. 1a) in 4-month-old park7−/− brains. Interestingly, the most upregulated gene of the HQ5’IRE gene set is alas2.
Enrichment analysis of IRE gene sets in park7−/− zebrafish brains. a Bar plot showing the significance of enrichment of the four IRE gene sets in 4-month-old park7−/− zebrafish brains. The green bar shows that genes with high quality 5’IREs are significantly altered as a group. The vertical line indicates when the FDR-adjusted p-value = 0.05. b Heatmap indicating the log2 fold change (logFC) of all HQ5’IRE genes detectable in the RNA-seq experiment
Discussion
Using our method to detect evidence for iron dyshomeostasis in RNA-seq data, we found highly significant alteration of the expression of genes with IREs in the 5’UTRs of their mRNAs in 4-month-old park7−/− brains.
Iron homeostasis is maintained by regulation of gene expression at several levels, including via transcription, mRNA stability, and mRNA translation (reviewed in [10]). The latter two phenomena are modulated by the iron regulatory proteins IRP1 and IRP2 when these bind to IREs. The highly significant alteration of expression of the HQ5′IRE gene set in 4-month-old park7−/− brains is likely due to changes in binding of IRP1 and/or IRP2 to IREs in the transcripts of these genes. However, since mutation of PARK7 is known to cause oxidative stress [11] and oxidative stress can also affect IRP formation (reviewed in [6]), it is difficult to differentiate between oxidative stress and iron dyshomeostasis as contributing to changes in HQ5′IRE gene transcript abundance. Indeed, since iron is so important for mitochondrial function, iron dyshomeostasis and oxidative stress often co-occur [12]. Interestingly, the HQ3′IRE gene set appeared unaffected in park7−/− brains and we currently have no explanation for why this should be so. However, if the effects on HQ5’IRE transcript abundance are, indeed, due to binding of IRPs, this points to the existence of mechanisms that can discriminate in binding of IRPs to IREs (or cause differences in the effects of such binding), depending on whether an IRE resides in the 5′ or 3′ UTR of a transcript.
Among members of the HQ5’IRE gene set, alas2 transcript levels were observed to be increased in park7−/− brains. The relationship between alas2 activity and IRPs is discussed in Additional file 1.
Taken together, our analysis of transcriptome data from park7−/− zebrafish brains supports the possibility of iron dyshomeostasis and/or oxidative stress as early preclinical events in PD. In future studies, we will explore the nature of the effects of PD-linked genes on iron homeostasis and mitochondrial function in zebrafish. Understanding these effects can provide mechanistic insight into PD for the development of therapeutics.
Availability of data and materials
The R code used to re-analyse the raw transcript counts of Hughes et al. can be found at https://github.com/karissa-b/dj1KO-RNAseq-IRE. The raw data from Hughes et al. is available from the Gene Expression Omnibus(GEO) database GSE135271 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE135271). The list of genes which contain an iron responsive element (IRE) in the untranslated regions of their mRNAs in zebrafish can be found at https://github.com/nhihin/ire.
Abbreviations
- PD:
-
Parkinson’s disease
- IRE:
-
Iron responsive element
- HQ:
-
High quality
- UTRs:
-
Untranslated
- GSH:
-
Glutathione
- SNc DA:
-
Substantia nigra pars compacta dopaminergic neurons
- ETC:
-
Electron transport chain
- ROS:
-
Reactive oxygen species
- IRP1:
-
Iron regulatory proteins 1
- IRP2:
-
Iron regulatory proteins 2
- mRNA:
-
Messenger ribonucleic acid
- AD:
-
Alzheimer’s disease
- RNA-seq:
-
RNA sequencing
- GSEA:
-
Gene set enrichment analysis
- alas2 :
-
Delta-aminolevulinate synthase 2
- FDR:
-
False discovery rate
References
Lee JH, Han J-h, Kim H, Park SM, Joe E-h, Jou I. Parkinson’s disease-associated LRRK2-G2019S mutant acts through regulation of SERCA activity to control ER stress in astrocytes. Acta Neuropathol Commun. 2019;7(1):68.
Andreeva A, Bekkhozhin Z, Omertassova N, Baizhumanov T, Yeltay G, Akhmetali M, et al. The apparent deglycase activity of DJ-1 results from the conversion of free methylglyoxal present in fast equilibrium with hemithioacetals and hemiaminals. J Biol Chem. 2019;294(49):18863–72.
Dolgacheva LP, Berezhnov AV, Fedotova EI, Zinchenko VP, Abramov AY. Role of DJ-1 in the mechanism of pathogenesis of Parkinson’s disease. J Bioenerg Biomembr. 2019;51(3):175–88.
Pacelli C, Giguère N, Bourque MJ, Lévesque M, Slack RS, Trudeau L. Elevated mitochondrial bioenergetics and axonal arborization size are key contributors to the vulnerability of dopamine neurons. Curr Biol. 2015;25(18):2349–60.
Junn E, Jang WH, Zhao X, Jeong BS, Mouradian MM. Mitochondrial localization of DJ-1 leads to enhanced neuroprotection. J Neurosci Res. 2009;87(1):123–9.
Urrutia PJ, Bórquez DA, Núñez MT. Inflaming the brain with iron. Antioxidants. 2021;10(1).
Hin N, Newman M, Pederson S, Lardelli M. Iron responsive element (IRE)-mediated responses to iron dyshomeostasis in Alzheimer’s disease. bioRxiv. 2021; 2020.05.01.071498.
Hughes GL, Lones MA, Bedder M, Currie PD, Smith SL, Pownall ME. Machine learning discriminates a movement disorder in a zebrafish model of Parkinson’s disease. Dis Models Mech. 2020;13(10):d045815.
Wu D, Lim E, Vaillant F, Asselin-Labat M-L, Visvader JE, Smyth GK. ROAST: rotation gene set tests for complex microarray experiments. Bioinformatics. 2010;26(17):2176–82.
Zhou ZD, Tan E-K. Iron regulatory protein (IRP)-iron responsive element (IRE) signaling pathway in human neurodegenerative diseases. Mol Neurodegener. 2017;12(1):75.
Canet-Avilés RM, Wilson MA, Miller DW, Ahmad R, McLendon C, Bandyopadhyay S, et al. The Parkinson’s disease protein DJ-1 is neuroprotective due to cysteine-sulfinic acid-driven mitochondrial localization. Proc Natl Acad Sci U S A. 2004;101(24):9103–8.
Yambire KF, Rostosky C, Watanabe T, Pacheu-Grau D, Torres-Odio S, Sanchez-Guerrero A, et al. Impaired lysosomal acidification triggers iron deficiency and inflammation in vivo. Elife. 2019;8.
Acknowledgements
We would like to thank Mary Elizabeth Pownall and Katherine Newling for making available the raw transcript counts of [8], and Nhi Hin for providing the sets of zebrafish genes containing iron-responsive elements.
Funding
KB is supported by an Australian Government Research Training Program Scholarship and by funds from the Carthew Family Charity Trust. ML is an academic employee of the University of Adelaide. LCP is an academic employee of the University of Adelaide and is further supported by a Barbara Kidman Fellowship.
Author information
Authors and Affiliations
Contributions
HYC drafted the manuscript. KB wrote the methodology, generated the diagrams and performed the bioinformatic analysis. ML, LCP and KB supervised and provided advice on the analysis of this paper. ML, KB, LCP and HYC edited 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.
Detailed Methodology and Additional Discussion.
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
Chin, H.Y., Lardelli, M., Collins-Praino, L. et al. Loss of park7 activity has differential effects on expression of iron responsive element (IRE) gene sets in the brain transcriptome in a zebrafish model of Parkinson’s disease. Mol Brain 14, 83 (2021). https://doi.org/10.1186/s13041-021-00792-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13041-021-00792-9
Keywords
- DJ-1
- PARK7
- Parkinson disease
- Zebrafish
- RNA-seq
- Enrichment analysis
- Iron dyshomeostasis