Potential biomarker identification for Friedreich’s ataxia using overlapping gene expression patterns in patient cells and mouse dorsal root ganglion
Authors:
Marissa Z. McMackin aff001; Blythe Durbin-Johnson aff002; Marek Napierala aff003; Jill S. Napierala aff003; Luis Ruiz aff001; Eleonora Napoli aff001; Susan Perlman aff004; Cecilia Giulivi aff001; Gino A. Cortopassi aff001
Authors place of work:
Department of Molecular Biosciences, University of California, Davis, Davis, California, United States of America
aff001; Bioinformatics, University of California, Davis, Davis, California, United States of America
aff002; Department of Biochemistry and Molecular Genetics, University of Alabama at Birmingham, Birmingham, Alabama, United States of America
aff003; Department of Neurology, University of California, Los Angeles, Los Angeles, California, United States of America
aff004
Published in the journal:
PLoS ONE 14(10)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0223209
Summary
Friedreich’s ataxia (FA) is a neurodegenerative disease with no approved therapy that is the result of frataxin deficiency. The identification of human FA blood biomarkers related to disease severity and neuro-pathomechanism could support clinical trials of drug efficacy. To try to identify human biomarkers of neuro-pathomechanistic relevance, we compared the overlapping gene expression changes of primary blood and skin cells of FA patients with changes in the Dorsal Root Ganglion (DRG) of the KIKO FA mouse model. As DRG is the primary site of neurodegeneration in FA, our goal was to identify which changes in blood and skin of FA patients provide a 'window' into the FA neuropathomechanism inside the nervous system. In addition, gene expression in frataxin-deficient neuroglial cells and FA mouse hearts were compared for a total of 5 data sets. The overlap of these changes strongly supports mitochondrial changes, apoptosis and alterations of selenium metabolism. Consistent biomarkers were observed, including three genes of mitochondrial stress (MTIF2, ENO2), apoptosis (DDIT3/CHOP), oxidative stress (PREX1), and selenometabolism (SEPW1). These results prompted our investigation of the GPX1 activity as a marker of selenium and oxidative stress, in which we observed a significant change in FA patients. We believe these lead biomarkers that could be assayed in FA patient blood as indicators of disease severity and progression, and also support the involvement of mitochondria, apoptosis and selenium in the neurodegenerative process.
Keywords:
Gene expression – Blood – Mouse models – mitochondria – apoptosis – biomarkers – Antioxidants – Selenium
Introduction
Friedreich’s ataxia (FA) is a lethal neurodegenerative disease with a pediatric onset, and the most frequently occurring autosomal recessive inherited ataxia. FA is caused by a trinucleotide (GAAn) repeat expansion in intron 1 of the nuclear encoded gene frataxin (FXN) [1], which results in a 90% decrease in mitochondrial frataxin’s gene expression [2]. The lower the frataxin level, the earlier the age of onset and the worse the severity of the disease [3]. Frataxin protein is thought to be important for Fe-S cluster biogenesis and heme synthesis [4], thus supporting metabolically active cells through mitochondrial functions. However how this Fe-S defect leads to neurodegeneration is less clear and could potentially be addressed through high-dimensional gene expression analysis as in the current study.
The primary site of neurodegeneration in FA patients is the dorsal root ganglion (DRG), and progressive neural tissue damage ascends the dorsal spine including the spinocerebellar tract to degenerate the cerebellum and leads to a loss of voluntary motor function [5, 6]. Results of autopsies from FA patients have clearly indicated that the DRGs are significantly implicated in the progressive pathology of the disease [7]. Because DRG tissue is unavailable from living patients, finding overlap of blood-based or skin-based human FA biomarkers with the FA pathoneurophysiological process could be important for understanding the disease progression. Accessible peripheral cells could provide a window into the neuropathophysiological process inside the nervous system.
The first symptom of the disease is usually motor coordination loss in lower extremities [8] making the lumbar spinal tissue of particular interest. Adolescence is the average age of onset in FA, when symptoms present most often with lower extremity loss of coordination [9]. By choosing to analyze gene expression changes in the mid-lumbar spinal cord mouse DRGs, we assessed the target tissues innervating the lower extremities where clinical signs first present in humans [10].
Although there is support that frataxin defects affects iron-sulfur biogenesis, how this leads to DRG neurodegeneration is not clear, so there is still urgent need to understand pathomechanism as well as investigate therapeutic avenues for patients. The KIKO mouse model is an established FA model that results in reduced FXN transcription as well as reduced frataxin protein expression [11]. The neurobehavioral phenotype of these mice has also recently been clarified by our lab [12]. There is a defect in mitochondrial biogenesis in KIKO mice which also occurs in FA patient skin and white blood cells [13].
Using adolescent KIKO mice we isolated the primary target tissue of lumbar (L1-L4) DRG and spinal cord, and conducted an RNA sequencing experiment with WT and KIKO littermates. The DRG transcriptome genes were then correlated to FXN gene expression within each sample because FXN gene expression is associated with disease severity [14, 15]. The genes that were found to be significantly correlated with FXN expression in the mouse DRG tissues were then compared collaboratively with independent data sets from patient samples for comparative gene expression profiles in FA human tissues.
Lymphocytes were collected from patient blood and were analyzed for gene expression changes in a comparison of control, carrier, and patient cells [16]. Gene expression data was collected on fibroblasts derived from FA patient or unaffected control skin biopsies [17] in a comparison of FA to control [18]. For additional independent confirmation of frataxin-dependent gene expression changes mouse glial and mouse heart array data were compared as fourth and fifth data sets for confirmation of genes likely to be significantly altered in FA.
There was substantial overlap in significantly altered gene expression across the five data sets. Those genes that overlapped in greater than 3 data sets were evaluated for their relative strength as potential biomarkers of FA. The list of overlapping altered genes was also analyzed for their functional association to each other using STRING interaction network mapping and resulted in a distinct groups of antioxidant genes, apoptosis regulation, translation regulation, and mitochondrial genes.
Materials and methods
Mouse DRG
Mice were bred in the UC Davis laboratory colony until tissue extraction at 4 months of age. All research involving mouse tissues has been approved by The Institutional Animal Care and Use Committee (IACUC) protocol #18070. Mice were maintained on a (12h light/12hr dark cycle), and given food and water ad libitum. A total of 7 knock-in knock-out (KIKO) mice (C57BL6/j; fxn^GAA230/-), and 10 littermate wild type mice were used for these data. Mice were euthanized using inhaled isoflurane, followed by cervical dislocation before rapid tissue removal and storage in RNAlater, Vilnius, Lithuania). Mouse DRG tissue was removed using a modified protocol for dorsal approach DRG removal of the lumbar spinal cord [19]. Total RNA was extracted from DRG lumbar tissue lysate using RNeasy Plus Mini Kit (Qiagen, Valencia, CA) following manufacturer’s instructions. The mRNA quantity was measured by NanoDrop 2000c Spectrophotometer (Thermo Scientific, Waltham, MA), and quality was measured by the 2100 Bioanalyzer using the Agilent RNA 6000 Nano Kit (Agilent Technologies).
RNA Sequencing
The DNA Technologies Core prepared barcoded mRNA-Seq libraries for each sample, pooled them, and ran the pool on one lane of HiSeq (PE100). After evaluation of the first lane for quality and number of reads per sample, additional lanes were run to increase the number of reads (or pairs of reads). For RNAseq studies in vertebrates, 10 million fragments should detect ~80% of annotated genes, while 30 million fragments should detect ~90% of annotated genes (ENCODE Consortium Standards, Guidelines and Best Practices for RNA-Seq, V1.0, June 2011;[20]). Illumina read quality assessment was performed using FastQC. Scythe and Sickle were used for Illumina adapter and quality trimming. Trimmed reads were aligned to the Mus musculus mm10 genome using Tophat2 [21]. Cufflinks2 [22] were used on the paired-end reads to identify potential novel splice variants. The raw counts were derived from the alignments using the HTSeq-count python script [23]. Statistical analyses (tests of differential expression and tests of correlation were conducted using EdgeR or Limma-Voom).
Genes with expression less than 0.2 counts per million reads were filtered prior to analysis, leaving 20,444 genes. Analyses of the correlation of frataxin expression with the expression of other genes were conducted by fitting linear regression models using the Limma-Voom Bioconductor pipeline. Raw p-value was calculated for the test that the logFC for frataxin is different from 0. Benjamini-Hochberg false discovery rate (FDR) is used to calculate the adjusted P-value.
Lymphocytes
The second independent set of data was derived from human blood samples collected from patients, carriers, and control individuals in a clinical setting at The UCLA Department of Neurology Program in Neurogenetics, and the UCLA Ataxia Center. Blood samples were collected from ~760 different individuals over ~ 5 years and results were also correlated to GAA repeat length to analyze peripheral gene expression data from patients with FA. All raw gene expression data is available for download in NCBI Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/gds) under accession number GSE102008.
All research involving patient/patient lymphocyte tissues was been approved by the Institutional Review Boards (IRB) protocol IRB#10–000833. The lymphocytes gene expression data was overlapped with the KIKO DRG gene expression data in search of commonly altered genes.
Fibroblast
The third separate data set was generated from RNA sequencing reactions performed on 18 FA patient and 17 control fibroblast cell lines [17]. Fibroblast repositories were created from more than 50 patient skin biopsies, and numerous controls, and transformed in to induced pluripotent stem cells (iPSCs) All research involving patient/patient tissues has been approved by the Institutional Review Boards (IRB) at The University of Alabama at Birmingham (IRB Protocols: N160923005 and N160922011) and Children’s Hospital of Philadelphia (IRB Protocol 10–007864). The fibroblast gene expression results were then analyzed for overlap with mouse DRG transcriptome results.
Mouse heart, glia
The fourth and fifth data sets were included as microarrays of Glial knockdowns of frataxin, and microarrays of mouse hearts lacking frataxin. Reverse transfection involving simultaneously transfecting and plating cells was performed using lipofectamine 2000 according to the manufacture protocol (Invitrogen). Briefly, transfection mix including siRNA (30nM final concentration, except 40nM for ND7/23) and transfection reagent was made and added to the wells in 6-well plates. The cells were harvested by trypsinization and about 0.3 million cells were plated in the wells with the transfection mix. siRNAs used in the study include: human frataxin siRNA described previously[20], AAC GUG GCC UCA ACC AGA UUU, and scrambled siRNA, CAG UCG CGU UUG CGA CUG GdTdT, human frataxin siRNA (ON-TARGETplus duplex #J-006691-07) and non-targeting siRNA (D-001810-0X), Rat frataxin siRNA (ON-TARGETplus SMARTpool, L-104901-01), mouse frataxin siRNA (ON-TARGETplus SMARTpool, #L-045500-00), and control siRNA (Non-Targeting Pool, #D-001810-10) from Dharmacon [24, 25]. Data was analyzed individually for each sample type using dChip v1.2. Within each group, samples were normalized to the median intensity chip and fluorescence values were generated using the perfect match-only model. Probe sets with a pCall > 40% and P < .05 were considered significantly altered. FDR testing was not performed.
GPX activity
Glutathione peroxidase activity in FA patient whole blood was evaluated by using a commercially available kit from Abcam (#ab102530) and according to the manufacturer’s protocol. Hemoglobin was evaluated by using the ferricyanide-cyanide reagent. Results were expressed as IU per g Hb.
Data analysis
All experiments were analyzed as a log fold change of FA cells compared to control, except for the mouse DRGs. Mouse DRGs results were analyzed as a regression analysis function of FXN expression. LogFC of DRG samples represents the change of a gene compared to FXN (+ r value represents the increase in a gene as FXN also increases). String Association Network Version 10.0 [26] tools were applied to overlap analysis of the 87 genes that were significant in 4/5 datasets. All settings were automated to default selections including a minimum required interaction score of 0.400, no clustering parameters specified, and with active interaction sources including: Textmining, Experiments, Databases, Co-expression, Neighborhood, Gene Fusion, and Co-occurrence. Network nodes represent all of “the proteins produced by a single protein-coding gene locus”. Edges represent protein-protein associations, “where proteins jointly contribute to a shared function” to determine the interactions shown. Further clarification of clusters associations and gene pathways were analyzed using DAVID Bioinformatics Resources [27], and using ENRICHR Comprehensive Gene Set Enrichment Analysis [28].
Results
There is substantial overlap of gene profile in frataxin-deficient cells and tissues. In the DRG transcriptome sequencing experiment of KIKO mice, 3961 genes were significantly correlated with FXN expression (FDR<0.05). In the human lymphocyte data set 2617 genes were significantly different in the comparison of patient to carriers (p<0.05). When DRG and lymphocyte data sets were compared to find common genes, 443 genes were found to be linked to frataxin expression on both lists. Next we incorporated expression data from FA patient fibroblasts, and found that 100 genes overlapped between the 3 experiments (Fig 1). We also tested the overlap of genes affected by frataxin depletion in FA mouse heart and FA glial cells (Fig 1). 87 most 'frataxin responsive' or 'consistently altered' transcripts differentially expressed in at least 4 experiments were identified (S1 Table).
Although there was substantial overlap both of gene expression changes, and sometimes of the direction of change, there are also cases in which transcript direction was different in patients vs. controls in different tissues. For the remainder of our analysis we disregarded those differences, picking those 87 genes as the most 'frataxin-responsive', irrespective of sign change positive or negative in the patient vs. control comparison (Fig 2 red = negative value, green = positive value).
STRING network analysis suggests multiple functional clusters. We then used STRING Network Association version 10.0 [26] to analyze the 87 frataxin-responsive genes. Three distinct clusters emerged (Fig 3). The functional enrichment results for our network (p-value = 0.00783) included biological processes and cellular components (FDR < .05) The three distinct clusters were: biological processes related to apoptosis regulation in Cluster 1, mitochondrial translation and transcription in Cluster 2, and antioxidant genes in Cluster 3 (Fig 4). 18 of the 87 shared genes were mitochondrial. The STRING Association Network Cellular Compartment mitochondrial genes are: MTIF2, TOMM40, TXNRD1, LPIN1, PTEN, FYN, KIF1B, CERK, RPL10A, ADSL, DCAKD, CYB5B, ECHDC2, SEPW1, SLC11A2, TTC3, HAX1, and SNN (FDR = 0.01) (Fig 5). Thus 20% of genes altered in frataxin-deficient cells across 4 tissues are mitochondrial. Frataxin is a nuclear-encoded gene that is translated on cytosolic ribosomes and translocated to mitochondria [15]. One interesting gene from the mitochondrial cluster is Translocase of Outer Mitochondrial Membrane 40 (TOMM40). TOMM40 was significantly lower in patient lymphocytes as compared to carriers, and in frataxin-knock down glial cells. These data imply that frataxin deficiency causes many mitochondrial changes, and that mitochondrial biomarkers are altered in multiple FA patient tissues, and multiple frataxin-deficient cell lines. The results on (mitochondrial) translation and mitochondria in general suggest a fundamental involvement of mitochondrial biogenesis and metabolism in Friedreich’s, which has recently been confirmed by identifying a mitochondrial biogenesis defect in FA mouse model KIKO, patient cells and patient blood [13]
The composition of cluster 1 includes 10 genes from the biological process pathway of Positive Regulation of Apoptotic Process (FDR = 0.0324). Genes that induce apoptosis also scored highly, they include DDIT3, JUN, PTEN, SMAD3, RHOB, SLC11A2, ADRB2, TSC22D1, RPS27A, and RNPS1. These 10 apoptosis related genes could potentially play a role downstream of frataxin loss and upstream of neurodegeneration. A predisposition to apoptosis in FA (Wong et al., 1999) may be induced as a result of the intrinsic pathway of apoptosis due to DNA damage [29]. Notably among the cluster 1 genes, DNA Damage Inducible Transcript 3 (DDIT3 = CHOP = GADD153) is significantly upregulated in 4/5 of the experiments (S3b Fig), and significantly correlated with FXN expression in KIKO mouse DRG (r = 0.56). Stimulation of DDIT3/CHOP/GADD153 stress response [30–33] has been shown to occur as the immediate result of mitochondrial inhibition and the mitochondrial stress response [30, 31, 34]. Thus DDIT3/CHOP shows great promise as a biomarker of the mitochondrial stress that occurs in the pathophysiological mechanism of FA. DDIT3/CHOP induction also occurs as a result of cellular stress, and multiple genes in our Cluster 2 results relate to translational stress.
The genes in cluster 2 are composed of the genes related to mitochondrial translation and transcription. They are translation-related genes MTIF2 (involved in mitochondrial translation), EIF4b, EEF2, RPL8, RPS27A, RPS2, and RPL10A; and transcription genes BTF3, and RNPS1 (Enrichr). In fact MTIF2 was conserved in 5/5 overlapping datasets, supporting a defect/alteration in mitochondrial protein translation (Fig 2, S1a and S1b Fig). We recently demonstrated a defect in mitochondrial biogenesis in Friedreich’s ataxia in FA cells, FA KIKO mouse tissues, and FA human patients[13], which could result from this consistent defect in MTIF, involved in mitochondrial protein translation. Eukaryotic Translation Initiation Factor 4b (EIF4B) are cluster 2 genes that are significantly altered in all 5 data sets (S1c Fig). The cluster is also composed of RPS27A and RPS2, which encode for ribosomal protein components of the 40s ribosomal subunit. RPL10A and RPL8 encode for ribosomal 60s subunit components. EEF2 promotes the GTP-dependent translocation of the nascent protein chain from the A-site to the P-site of the ribosome [35]. BTF3 basic transcription factor 3 is required for transcriptional initiation. RNPS1 detects incomplete translation as an mRNA nuclear transport and mRNA surveillance protein, detecting truncated mRNA and initiates nonsense-mediated mRNA decay. Thus these high dimensional studies suggest a (mitochondrial) protein translation defect, and in fact other studies have recently confirmed a mitochondrial biogenesis defect in multiple FA mouse and human tissues.
The genes in cluster 3 are primarily antioxidant genes related to selenium and glutathione. Multiple antioxidant response genes altered are related to selenium metabolism and glutathione peroxidase activity: SEPW1, GPX7, and TXNRD1 (Fig 4), and if one reduces the stringency to 3/5 experiments, selenoprotein M (SELM) was also significant (p<0.05). SEPW1 and SELM are both biomarkers of bioavailable selenium [36–42], and both genes are down regulated when FXN is reduced in the DRGs. Because selenocysteine is at the active site of many selenium-dependent antioxidant enzymes, deficiency in selenium results in reduced antioxidant activity, including glutathione peroxidase activity [43–45], and the gene expression changes observed in our data may be important to understanding the involvement of frataxin as an antioxidant. Alanyl aminopeptidase (ANPEP) has been shown to positively regulate glutathione synthase [46], and its expression was also altered in 4/5 of the experiments (Fig 4e). Thioredoxin reductase 1 (TXNRD1) is a selenoenzyme that protects against oxidative stress [47], [48] which was significant in 4/5 of the experiments. Thus, frataxin-deficiency is associated with multiple parameters related to oxidative stress and selenometabolism. If frataxin deficiency alters selenium bioavailability, then reduced selenium bioavailability could potentially decrease the selenium-dependent antioxidant functions of thioredoxin reductase and other glutathione peroxidases[49], underlying the intrinsic sensitivity of FA patient cells to oxidant stress [50, 51].
For direct confirmation of the glutathione peroxidase antioxidant pathway in frataxin-deficient cells, a 6x6 study of FA patient whole blood was compared to control whole blood samples. There was no difference between controls and patients in hemoglobin (mg/mL) (p = 0.400), but GPx activity was significantly increased (IU/g Hb) in patient cells (p = 0.047) (Fig 5).
Search for overlapping genes with available blood tests. Since the underlying goal of this work is to identify new biomarkers that could be used clinically, we decreased screening stringency to altered in 3/5 of the experiments, which increased the total gene list, and then we crossed this expanded gene list with clinical blood tests available and orderable at the Mayo clinic diagnostic testing labs. The genes altered in 3/5 experiments and also with an orderable blood test from Mayo include: ENO2, LPIN1, SLC22A5, ATP6AP2, PTEN, SMAD3, TMEM127, ADSL, AARS, HAX1, and KIF1B. ENO2 (Mayo Clinic), and lipin 1 (LPIN1) (S4a Fig), are potentially interesting in terms of clinical monitoring of FA disease progression. ENO2 is found in mature neurons, and is linked to brain iron accumulation associated neurodegeneration [52]. Lipin 1 (LPIN1) was significantly altered in 4/5 of the data sets. Abnormal LPIN1 is associated with metabolic syndromes, vacuole regulation, and diabetes [53, 54]. Integral aspects of FA patient pathologies include increased likelihood of developing diabetes, metabolic abnormalities, and abnormal vacuole accumulation in the DRGs.
Discussion
Here we conducted differential expression analyses of clinical and preclinical models of Friedreich’s ataxia. Our results suggest defects in mitochondrial translation and selenium metabolism are part of the pathophysiological process that leads to apoptosis and neurodegeneration. Since these clusters were identified in both clinical and preclinical samples, we suggest that they represent promising targets for further pathomechanistic studies. These results suggest that decreased FXN may alter selenometabolism which could explain the translational deficits (cluster 2), and may decrease the activity of antioxidant selenoenzymes (cluster 3), increasing the burden of oxidative stress, and leading to increased apoptosis (cluster 1). Thus, our results further support a more specific formulation of the oxidative stress hypothesis for FA, which has long been hypothesized to be a key mechanism of FA [55].
Transcription and Translation Cluster: The apoptosis genes may connect cluster 1 to cluster 2. This implies that the regulation of cell death is inextricably linked to the transcription and translation changes seen in frataxin-deficient cells. The GAA repeat expansion in FXN causes R-Loop formation [56], which leads to transcriptional repression. Basic transcription factor 3 (BTF3), which was significantly altered in 4/5 of the experiments, is required for transcriptional initiation, and may be responding to the need for increased transcription of FXN. Additionally, RNPS1 was a significantly altered gene in 4/5 of the experiments. RNPSI is integral to mRNA nuclear transport and mRNA surveillance, which may be detecting truncated frataxin mRNA and initiating nonsense-mediated mRNA decay.
Expression of translation genes was also altered, suggesting that cellular compensation for frataxin deficiency has an impact on translation efficiency of mitochondrial and cytoplasmic translation. Most mitochondrial proteins are translated on the cytoplasmic ribosomes, including Frataxin and other Fe-S cluster proteins [57]. Our results included significantly altered cytoplasmic translation genes: EEF2, RPL8, RPS27A, RPS2, and RPL10A. MTIF2 is a mitochondrial gene that is altered in 5/5 of the experiments. MTIF2 transfers mt-mRNAs into the mitoribosome [58] and translocation of frataxin to the mitochondria appears to be integral to oxidative phosphorylation [59, 60].
A selenium-dependent cluster of genes: Relationship to glutathione and oxidative stress and potential use as a biomarker
Selenometabolism cluster: Selenium metabolism is intimately related with oxidative stress, as many important antioxidants enzymes (thioredoxin reductase, glutathione peroxidase) require the abnormal amino acid selenocysteine in order to carry out their activities. Selenium is a constitutive component of glutathione peroxidase [61]. If frataxin deficiency reduces the bioavailability of selenocysteine, causing selenocysteine deficiency, then one might predict multiple consequences. First, since selenocysteine is a rare but used amino acid with its own selenocysteine-tRNA, one might expect translational defects, and translation was a major functional cluster (Fig 3). The deficiency of SEPW1 and SELM in FA patient blood is consistent with a deficiency of bioavailable selenocysteine in FA patients.
Second, a defect in selenocysteine bioavailability would be predicted to reduce the antioxidant activity of multiple antioxidant selenoenzymes, including TXNRD1 and GPX7 which were both significantly altered in 4/5 of the experiments. Glutathione peroxidases and thioredoxin have also previously been shown to be significantly altered in the blood samples of FA patients [62]. These selenoenzymes were also found to be reduced in DRG of another FA mouse model, YG8 [49]. In support of an alteration in selenometabolism in FA, we observed a significant increase in GPX1 activity (Fig 5). These data tend to support an alteration of selenometabolism, possibly GPX1 activity is induced in response to a deficiency in the antioxidant selenium.
Alterations in selenoenzymes have been reported by others, including Helveston et al, 1996 [63]. Altered bioavailability of selenium results in reduced glutathione peroxidase activity, which has selenocysteine at its active site [45, 64, 65]. A deficiency in bioavailable selenium or selenocysteine could be an important consequence of frataxin deficiency, because many selenoenzymes are antioxidants. Additionally, selenium supplementation has been shown to increase the viability of FA fibroblast’s viability [66]. SEPW1, and SELM to some extent, are biomarkers of selenium status [67], and TXNRD1, thioredoxin reductase 1, is a selenium containing enzyme. Many glutathione peroxidases require selenium [64], GPX7 does not [68]. The selenoenzymes that are thought to provide antioxidant protection against ROS induced cellular damage are GPX—GPX1, GPX3, GPX4, GPX5, and GPX6 [69], and there is substantial support for antioxidant deficiency in FA [34, 49, 70].
Several markers of selenium status and selenium-related enzymes were significantly overlapping in 3/5 experiments. Those have not been included in these results to remove speculation. However, animals deficient in selenium or vitamin E develop white muscle disease, a myodegenerative condition considered to result from defects in antioxidant activity [71, 72]. Selenoprotein W was altered in 4/5 experiments. Selenoprotein W has been shown by [73] Flohe, and Steinbrenner) as the most responsive blood biomarker to selenium status. Based on the observed data we propose that frataxin deficiency may be associated with a change in bioavailable selenium or selenocysteine, which results in a deficiency in selenoenzymes and an antioxidant enzyme deficiency. Selenium level in blood and serum is available to be ordered through Mayo Clinic laboratories.
Apoptosis Regulation Cluster: From the STRING analysis, multiple apoptosis regulation genes were altered. These are consistent with the idea that frataxin deficiency may cause oxidative stress to the cells, result in DNA damage, and increase sensitivity to apoptosis [50]. DDIT3 and PTEN are notable transcripts from our cluster 1 results. DDIT3 (also known as CHOP), which we observed to be significantly altered in 4/5 of our data has already been shown to be activated downstream of mitochondrial stress [4, 34]. DDIT3 has also been linked to global shut down type of cellular response to stress [32].
PTEN was significantly altered in 4/5 of the data sets. PTEN has been linked to PREX1 activation [74](S4c Fig). PREX1 = Phosphatidylinositol-3,4,5-Trisphosphate Dependent Rac Exchange Factor 1 was also significantly changed in 3/5 of the experiments. PREX1 has been shown to bind to Rac1, which is a NADPH oxidase [75], and we have previously observed PREX1 to be altered in a smaller screen for FA biomarkers [76]. PREX1 induces PI3K inhibition-induced apoptosis and also has potential as a stress response biomarker. PREX1 was significant in 3/5 of the experiments. PREX1 has been shown to bind to Rac1, which is a NADPH oxidase [75], and we have also previously observed PREX1 to be altered in a smaller screen for FA biomarkers [76].
Other potential blood based biomarkers of FA that already exist as blood tests is ENO2. We observed significantly altered expression of ENO2 in 3/5 experiments. Enolase 2 = Neuron-specific enolase is expressed in neurons and lymphocytes, and was one of the most altered genes in FA blood. It has also been shown that increased ENO2 corresponds to demyelination [77]. ENO2 is known to play a role in at least one other neuronal disease that results in motor coordination loss [78]. Recently, serum neuron specific enolase (NSE) levels have been reported as a biomarker for multiple sclerosis disease progression [79]. ENO2 is a clinically-run blood test (Mayo Clinic, Quest Diagnostics), and could be used to assess markers of mitochondrial function and/or myelination status in FA patients. If we relax constraints to 3/5 experiments and then cross these by blood tests available at Mayo clinic the list of blood tests include: ENO2, LPIN1, SLC22A5, ATP6AP2, PTEN, SMAD3, TMEM127, ADSL, AARS, HAX1, and KIF1B. Thus these are potential biomarkers to be explored in FA for which commercially available clinical tests exist.
List of potential biomarkers with support from blood to mice. Our results suggest a group of biomarkers to be tested in FA patients vs. controls, including these four seleno-related: SEPW1, SELM1, TXRD1, ANPEP, and the antioxidant-related PREX1. As a marker of apoptosis and mitochondrial defects DDIT3 also shows potential. ENO2 is currently available from clinical labs, was quite altered in blood, and shows promise. MTIF2 was a significant biomarker in 5/5 of the experiments and could potentially act to clarify ongoing mitochondrial changes as a result of disease progression or therapy.
Conclusions
Overlapping patterns of gene expression in 5 independent frataxin-deficient tissues led to a limited set of 5 potential biomarkers. Cluster analysis identified a mitochondrial transcription translation cluster, and a selenium metabolism cluster, and an apoptotic cluster. Multiple genes that could be considered as biomarkers of frataxin deficiency in blood, that agreed with other tissues, were identified, and these included DDIT3, MTIF2, SEPW1, ENO2, and PREX1. Of these, blood tests for ENO2, Lpin1, Serum Selenium, and Blood Selenium are clinically available blood tests from Mayo Clinic laboratories, and thus could be included in new biomarker clinical trials for Friedreich’s. These 'most consistent' biomarkers appear to be ready for testing in humans, and could be tested in blood relative to repeat length, extent of frataxin deficiency, age of onset, and clinical assessment of motor score. Furthermore, these data suggest pathomechanistically that frataxin deficiency perhaps through selenium metabolism alters mitochondrial biogenesis and consequently triggers apoptosis and neurodegeneration.
Supporting information
S2 Fig [tif]
A substantial number of mitochondrial genes are significantly associated with FXN expression in 4/5 gene sets.
S1 Table [tif]
The 87 genes that were significantly changed in 4/5 experiments are shown as raw values.
Zdroje
1. Pandolfo M. Friedreich’s ataxia: clinical aspects and pathogenesis. Semin Neurol. 1999;19(3):311–21. Epub 2002/08/27. doi: 10.1055/s-2008-1040847 12194387.
2. Campuzano V, Montermini L, Lutz Y, Cova L, Hindelang C, Jiralerspong S, et al. Frataxin is Reduced in Friedreich Ataxia Patients and is Associated with Mitochondrial Membranes. Human Molecular Genetics. 1997;6(11):1771–80. doi: 10.1093/hmg/6.11.1771 9302253
3. Durr A, Cossee M, Agid Y, Campuzano V, Mignard C, Penet C, et al. Clinical and genetic abnormalities in patients with Friedreich’s ataxia. N Engl J Med. 1996;335(16):1169–75. Epub 1996/10/17. doi: 10.1056/NEJM199610173351601 8815938.
4. Napoli E, Taroni F, Cortopassi GA. Frataxin, iron-sulfur clusters, heme, ROS, and aging. Antioxid Redox Signal. 2006;8(3–4):506–16. Epub 2006/05/09. doi: 10.1089/ars.2006.8.506 16677095.
5. Filla A, De Michele G, Cavalcanti F, Pianese L, Monticelli A, Campanella G, et al. The relationship between trinucleotide (GAA) repeat length and clinical features in Friedreich ataxia. Am J Hum Genet. 1996;59(3):554–60. 8751856.
6. Marmolino D. Friedreich’s ataxia: past, present and future. Brain Res Rev. 2011;67(1–2):311–30. Epub 2011/05/10. doi: 10.1016/j.brainresrev.2011.04.001 21550666.
7. De Biase I, Rasmussen A, Endres D, Al-Mahdawi S, Monticelli A, Cocozza S, et al. Progressive gaa expansions in dorsal root ganglia of Friedreich’s ataxia patients. Annals of Neurology. 2007;61(1):55–60. doi: 10.1002/ana.21052 17262846
8. Pandolfo M. Friedreich ataxia: the clinical picture. Journal of neurology. 2009;256 Suppl 1:3–8. Epub 2009/04/11. doi: 10.1007/s00415-009-1002-3 19283344.
9. Delatycki MB, Paris DB, Gardner RJ, Nicholson GA, Nassif N, Storey E, et al. Clinical and genetic study of Friedreich ataxia in an Australian population. Am J Med Genet. 1999;87(2):168–74. Epub 1999/10/26. doi: 10.1002/(sici)1096-8628(19991119)87:2<168::aid-ajmg8>3.0.co;2-2 10533031.
10. Rigaud M, Gemes G, Barabas ME, Chernoff DI, Abram SE, Stucky CL, et al. Species and strain differences in rodent sciatic nerve anatomy: Implications for studies of neuropathic pain. Pain. 2008;136(1–2):188–201. doi: 10.1016/j.pain.2008.01.016 18316160.
11. Miranda CJ, Santos MM, Ohshima K, Smith J, Li L, Bunting M, et al. Frataxin knockin mouse. FEBS Lett. 2002;512(1–3):291–7. Epub 2002/02/20. doi: 10.1016/s0014-5793(02)02251-2 11852098.
12. McMackin MZ, Henderson CK, Cortopassi GA. Neurobehavioral deficits in the KIKO mouse model of Friedreich’s ataxia. Behav Brain Res. 2017;316:183–8. Epub 2016/08/31. doi: 10.1016/j.bbr.2016.08.053 27575947.
13. Jasoliya MJ, McMackin MZ, Henderson CK, Perlman SL, Cortopassi GA. Frataxin Deficiency Impairs Mitochondrial Biogenesis in Cells, Mice and Humans. Hum Mol Genet. 2017. Epub 2017/04/27. doi: 10.1093/hmg/ddx141 28444186.
14. Evans-Galea MV, Lockhart PJ, Galea CA, Hannan AJ, Delatycki MB. Beyond loss of frataxin: the complex molecular pathology of Friedreich ataxia. Discov Med. 2014;17(91):25–35. Epub 2014/01/15. 24411698.
15. Santos R, Lefevre S, Sliwa D, Seguin A, Camadro JM, Lesuisse E. Friedreich Ataxia: Molecular Mechanisms, Redox Considerations, and Therapeutic Opportunities. Antioxid Redox Signal. 2010;13(5):651–90. doi: 10.1089/ars.2009.3015 20156111.
16. Nachun D, Gao F, Isaacs C, Strawser C, Yang Z, Dokuru D, et al. Peripheral blood gene expression reveals an inflammatory transcriptomic signature in Friedreich’s ataxia patients. Hum Mol Genet. 2018. Epub 2018/05/24. doi: 10.1093/hmg/ddy198 29790959.
17. Li Y, Lu Y, Polak U, Lin K, Shen J, Farmer J, et al. Expanded GAA repeats impede transcription elongation through the FXN gene and induce transcriptional silencing that is restricted to the FXN locus. Hum Mol Genet. 2015;24(24):6932–43. Epub 2015/09/25. doi: 10.1093/hmg/ddv397 26401053.
18. Napierala JS, Li Y, Lu Y, Lin K, Hauser LA, Lynch DR, et al. Comprehensive analysis of gene expression patterns in Friedreich’s ataxia fibroblasts by RNA sequencing reveals altered levels of protein synthesis factors and solute carriers. 2017;10(11):1353–69. doi: 10.1242/dmm.030536 29125828.
19. Sleigh JN, Weir GA, Schiavo G. A simple, step-by-step dissection protocol for the rapid isolation of mouse dorsal root ganglia. BMC Res Notes. 2016;9. doi: 10.1186/s13104-016-1915-8 26864470.
20. Wang Y, Ghaffari N, Johnson CD, Braga-Neto UM, Wang H, Chen R, et al. Evaluation of the coverage and depth of transcriptome by RNA-Seq in chickens. BMC Bioinformatics. 2011;12(10):S5. doi: 10.1186/1471-2105-12-S10-S5 22165852
21. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biology. 2013;14(4):R36. doi: 10.1186/gb-2013-14-4-r36 23618408
22. Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011;27(17):2325–9. Epub 2011/06/24. doi: 10.1093/bioinformatics/btr355 21697122.
23. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. Epub 2014/09/28. doi: 10.1093/bioinformatics/btu638 25260700.
24. Lu C, Schoenfeld R, Shan Y, Tsai C, Hammock B, Cortopassi G. Frataxin Deficiency Induces Schwann Cell Inflammation and Death. Biochim Biophys Acta. 2009;1792(11):1052–61. doi: 10.1016/j.bbadis.2009.07.011 19679182.
25. Schoenfeld RA, Napoli E, Wong A, Zhan S, Reutenauer L, Morin D, et al. Frataxin deficiency alters heme pathway transcripts and decreases mitochondrial heme metabolites in mammalian cells. Hum Mol Genet. 2005;14(24):3787–99. Epub 2005/10/22. doi: 10.1093/hmg/ddi393 16239244.
26. Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41(Database issue):D808–15. Epub 2012/12/04. doi: 10.1093/nar/gks1094 23203871.
27. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57. Epub 2009/01/10. doi: 10.1038/nprot.2008.211 19131956.
28. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Research. 2016. doi: 10.1093/nar/gkw377 27141961
29. Igoillo-Esteve M, Gurgul-Convey E, Hu A, Romagueira Bichara Dos Santos L, Abdulkarim B, Chintawar S, et al. Unveiling a common mechanism of apoptosis in β-cells and neurons in Friedreich’s ataxia. Human Molecular Genetics. 2015;24(8):2274–86. doi: 10.1093/hmg/ddu745 25552656
30. Harding HP, Novoa I, Zhang Y, Zeng H, Wek R, Schapira M, et al. Regulated translation initiation controls stress-induced gene expression in mammalian cells. Mol Cell. 2000;6(5):1099–108. Epub 2000/12/07. doi: 10.1016/s1097-2765(00)00108-8 11106749.
31. Oyadomari S, Mori M. Roles of CHOP//GADD153 in endoplasmic reticulum stress. Cell Death Differ. 2003;11(4):381–9.
32. Paschen W, Proud CG, Mies G. Shut-down of translation, a global neuronal stress response: mechanisms and pathological relevance. Curr Pharm Des. 2007;13(18):1887–902. Epub 2007/06/23. doi: 10.2174/138161207780858401 17584115.
33. Sok J, Wang XZ, Batchvarova N, Kuroda M, Harding H, Ron D. CHOP-Dependent Stress-Inducible Expression of a Novel Form of Carbonic Anhydrase VI. Mol Cell Biol. 1999;19(1):495–504. doi: 10.1128/mcb.19.1.495 9858573.
34. Silva JM, Wong A, Carelli V, Cortopassi GA. Inhibition of mitochondrial function induces an integrated stress response in oligodendroglia. Neurobiol Dis. 2009;34(2):357–65. Epub 2009/02/24. doi: 10.1016/j.nbd.2009.02.005 19233273.
35. Kaul G, Pattan G, Rafeequi T. Eukaryotic elongation factor-2 (eEF2): its regulation and peptide chain elongation. Cell biochemistry and function. 2011;29(3):227–34. Epub 2011/03/12. doi: 10.1002/cbf.1740 21394738.
36. Steinbrenner H, Speckmann B, Klotz LO. Selenoproteins: Antioxidant selenoenzymes and beyond. Arch Biochem Biophys. 2016;595:113–9. Epub 2016/04/21. doi: 10.1016/j.abb.2015.06.024 27095226.
37. De Spirt S, Eckers A, Wehrend C, Micoogullari M, Sies H, Stahl W, et al. Interplay between the chalcone cardamonin and selenium in the biosynthesis of Nrf2-regulated antioxidant enzymes in intestinal Caco-2 cells. Free Radic Biol Med. 2016;91:164–71. Epub 2015/12/25. doi: 10.1016/j.freeradbiomed.2015.12.011 26698667.
38. Combs G. Who Can Benefit from Selenium? Diversity of Selenium Functions in Health and Disease. Oxidative Stress and Disease: CRC Press; 2015. p. 3–16.
39. Gladyshev VN, Arner ES, Berry MJ, Brigelius-Flohe R, Bruford EA, Burk RF, et al. Selenoprotein Gene Nomenclature. J Biol Chem. 2016;291(46):24036–40. Epub 2016/09/21. doi: 10.1074/jbc.M116.756155 27645994.
40. Brigelius-Flohe R, Flohe L. Selenium and redox signaling. Arch Biochem Biophys. 2016. Epub 2016/08/09. doi: 10.1016/j.abb.2016.08.003 27495740.
41. Brigelius-Flohe R. The evolving versatility of selenium in biology. Antioxid Redox Signal. 2015;23(10):757–60. Epub 2015/09/26. doi: 10.1089/ars.2015.6469 26406357.
42. Steinbrenner H. Interference of selenium and selenoproteins with the insulin-regulated carbohydrate and lipid metabolism. Free radical biology & medicine. 2013;65:1538–47. Epub 2013/07/23. doi: 10.1016/j.freeradbiomed.2013.07.016 23872396.
43. Auchere F, Santos R, Planamente S, Lesuisse E, Camadro JM. Glutathione-dependent redox status of frataxin-deficient cells in a yeast model of Friedreich’s ataxia. Hum Mol Genet. 2008;17(18):2790–802. Epub 2008/06/20. doi: 10.1093/hmg/ddn178 18562474.
44. Abeti R, Parkinson MH, Hargreaves IP, Angelova PR, Sandi C, Pook MA, et al. 'Mitochondrial energy imbalance and lipid peroxidation cause cell death in Friedreich/’s ataxia'. Cell Death Dis. 2016;7:e2237. doi: 10.1038/cddis.2016.111 27228352
45. Fryer MJ. Rationale for clinical trials of selenium as an antioxidant for the treatment of the cardiomyopathy of Friedreich’s ataxia. Medical hypotheses. 2002;58(2):127–32. Epub 2002/01/29. doi: 10.1054/mehy.2001.1474 11812188.
46. Yim H-S, Cho YS, Guang X, Kang SG, Jeong J-Y, Cha S-S, et al. Minke whale genome and aquatic adaptation in cetaceans. Nat Genet. 2014;46(1):88–92. http://www.nature.com/ng/journal/v46/n1/abs/ng.2835.html#supplementary-information. 24270359
47. Dato S, De Rango F. Antioxidants and Quality of Aging: Further Evidences for a Major Role of TXNRD1 Gene Variability on Physical Performance at Old Age. 2015;2015:926067. doi: 10.1155/2015/926067 26064428.
48. Kiermayer C, Michalke B, Schmidt J, Brielmeier M. Effect of selenium on thioredoxin reductase activity in Txnrd1 or Txnrd2 hemizygous mice. Biol Chem. 2007;388(10):1091–7. Epub 2007/10/17. doi: 10.1515/BC.2007.133 17937623.
49. Shan Y, Schoenfeld RA, Hayashi G, Napoli E, Akiyama T, Iodi Carstens M, et al. Frataxin Deficiency Leads to Defects in Expression of Antioxidants and Nrf2 Expression in Dorsal Root Ganglia of the Friedreich’s Ataxia YG8R Mouse Model. Antioxidants & redox signaling. 2013;19(13):1481–93. doi: 10.1089/ars.2012.4537 23350650.
50. Wong A, Yang J, Cavadini P, Gellera C, Lonnerdal B, Taroni F, et al. The Friedreich’s ataxia mutation confers cellular sensitivity to oxidant stress which is rescued by chelators of iron and calcium and inhibitors of apoptosis. Hum Mol Genet. 1999;8(3):425–30. Epub 1999/02/09. doi: 10.1093/hmg/8.3.425 9949201.
51. Chantrel-Groussard K, Geromel V, Puccio H, Koenig M, Munnich A, Rotig A, et al. Disabled early recruitment of antioxidant defenses in Friedreich’s ataxia. Hum Mol Genet. 2001;10(19):2061–7. Epub 2001/10/09. doi: 10.1093/hmg/10.19.2061 11590123.
52. Takano K, Shiba N, Wakui K, Yamaguchi T, Aida N, Inaba Y, et al. Elevation of neuron specific enolase and brain iron deposition on susceptibility-weighted imaging as diagnostic clues for beta-propeller protein-associated neurodegeneration in early childhood: Additional case report and review of the literature. Am J Med Genet A. 2016;170a(2):322–8. Epub 2015/10/21. doi: 10.1002/ajmg.a.37432 26481852.
53. Sasser T, Qiu QS, Karunakaran S, Padolina M, Reyes A, Flood B, et al. Yeast lipin 1 orthologue pah1p regulates vacuole homeostasis and membrane fusion. The Journal of biological chemistry. 2012;287(3):2221–36. Epub 2011/11/29. doi: 10.1074/jbc.M111.317420 22121197.
54. Bego T, Dujic T, Mlinar B, Semiz S, Malenica M, Prnjavorac B, et al. Association of LPIN1 gene variations with markers of metabolic syndrome in population from Bosnia and Herzegovina. Med Glas (Zenica). 2015;12(2):113–21. Epub 2015/08/16. 26276647.
55. Sanchez N, Chapdelaine P, Rousseau J, Raymond F, Corbeil J, Tremblay JP. Characterization of frataxin gene network in Friedreich’s ataxia fibroblasts using the RNA-Seq technique. Mitochondrion. 2016;30:59–66. http://dx.doi.org/10.1016/j.mito.2016.06.003. 27350085
56. Groh M, Lufino MM, Wade-Martins R, Gromak N. R-loops associated with triplet repeat expansions promote gene silencing in Friedreich ataxia and fragile X syndrome. PLoS Genet. 2014;10(5):e1004318. Epub 2014/05/03. doi: 10.1371/journal.pgen.1004318 24787137.
57. Stemmler TL, Lesuisse E, Pain D, Dancis A. Frataxin and Mitochondrial FeS Cluster Biogenesis. J Biol Chem. 2010;285(35):26737–43. doi: 10.1074/jbc.R110.118679 20522547.
58. Audano M, Ferrari A, Fiorino E, Kuenzl M, Caruso D, Mitro N, et al. Energizing Genetics and Epi-genetics: Role in the Regulation of Mitochondrial Function. Curr Genomics. 2014;15(6):436–56. doi: 10.2174/138920291506150106151119 25646072.
59. Smits P, Smeitink J, van den Heuvel L. Mitochondrial Translation and Beyond: Processes Implicated in Combined Oxidative Phosphorylation Deficiencies. Journal of Biomedicine and Biotechnology. 2010;2010:24. doi: 10.1155/2010/737385 20396601
60. Smeitink J, van den Heuvel L, DiMauro S. The genetics and pathology of oxidative phosphorylation. Nat Rev Genet. 2001;2(5):342–52. http://www.nature.com/nrg/journal/v2/n5/suppinfo/nrg0501_342a_S1.html. 11331900
61. Rotruck JT, Pope AL, Ganther HE, Swanson AB, Hafeman DG, Hoekstra WG. Selenium: biochemical role as a component of glutathione peroxidase. Science. 1973;179(4073):588–90. Epub 1973/02/09. doi: 10.1126/science.179.4073.588 4686466.
62. Piemonte F, Pastore A, Tozzi G, Tagliacozzi D, Santorelli FM, Carrozzo R, et al. Glutathione in blood of patients with Friedreich’s ataxia. Eur J Clin Invest. 2001;31(11):1007–11. Epub 2001/12/12. doi: 10.1046/j.1365-2362.2001.00922.x 11737244.
63. Helveston W, Hurd R, Uthman B, Wilder BJ. Abnormalities of glutathione peroxidase and glutathione reductase in four patients with Friedreich’s disease. Mov Disord. 1996;11(1):106–7. Epub 1996/01/01. doi: 10.1002/mds.870110124 8771081.
64. Baker RD, Baker SS, LaRosa K, Whitney C, Newburger PE. Selenium regulation of glutathione peroxidase in human hepatoma cell line Hep3B. Arch Biochem Biophys. 1993;304(1):53–7. Epub 1993/07/01. doi: 10.1006/abbi.1993.1320 8391784.
65. Kirches E, Andrae N, Hoefer A, Kehler B, Zarse K, Leverkus M, et al. Dual role of the mitochondrial protein frataxin in astrocytic tumors. Lab Invest. 2011;91(12):1766–76. doi: 10.1038/labinvest.2011.130 21863062
66. Jauslin ML, Wirth T, Meier T, Schoumacher F. A cellular model for Friedreich Ataxia reveals small-molecule glutathione peroxidase mimetics as novel treatment strategy. Hum Mol Genet. 2002;11(24):3055–63. Epub 2002/11/06. doi: 10.1093/hmg/11.24.3055 12417527.
67. Kohrl J, Brigelius-Flohe R, Bock A, Gartner R, Meyer O, Flohe L. Selenium in biology: facts and medical perspectives. Biol Chem. 2000;381(9–10):849–64. Epub 2000/11/15. doi: 10.1515/BC.2000.107 11076017.
68. Brigelius-Flohe R, Maiorino M. Glutathione peroxidases. Biochim Biophys Acta. 2013;1830(5):3289–303. Epub 2012/12/04. doi: 10.1016/j.bbagen.2012.11.020 23201771.
69. Tinggi U. Selenium: its role as antioxidant in human health. Environ Health Prev Med. 2008;13(2):102–8. doi: 10.1007/s12199-007-0019-4 19568888.
70. Carletti B, Piemonte F. Friedreich’s Ataxia: A Neuronal Point of View on the Oxidative Stress Hypothesis. Antioxidants. 2014;3(3):592–603. doi: 10.3390/antiox3030592 26785073.
71. Ishihara H, Kanda F, Matsushita T, Chihara K, Itoh K. White muscle disease in. J Neurol Neurosurg Psychiatry. 1999;67(6):829–30. doi: 10.1136/jnnp.67.6.829 10617385.
72. Koller LD, Exon JH. The two faces of selenium-deficiency and toxicity—are similar in animals and man. Can J Vet Res. 1986;50(3):297–306. Epub 1986/07/01. 3527390.
73. Kipp AP, Frombach J, Deubel S, Brigelius-Flohe R. Selenoprotein W as biomarker for the efficacy of selenium compounds to act as source for selenoprotein biosynthesis. Methods in enzymology. 2013;527:87–112. Epub 2013/07/09. doi: 10.1016/B978-0-12-405882-8.00005-2 23830627.
74. Dillon LM, Bean JR, Yang W, Shee K, Symonds LK, Balko JM, et al. P-REX1 creates a positive feedback loop to activate growth factor receptor, PI3K/AKT, and MEK/ERK signaling in breast cancer. Oncogene. 2015;34(30):3968–76. doi: 10.1038/onc.2014.328 25284585.
75. Damoulakis G, Gambardella L, Rossman KL, Lawson CD, Anderson KE, Fukui Y, et al. P-Rex1 directly activates RhoG to regulate GPCR-driven Rac signalling and actin polarity in neutrophils. Journal of Cell Science. 2014;127(11):2589–600. doi: 10.1242/jcs.153049 24659802
76. Hayashi G, Cortopassi G. Lymphoblast Oxidative Stress Genes as Potential Biomarkers of Disease Severity and Drug Effect in Friedreich’s Ataxia. PLoS ONE. 2016;11(4). doi: 10.1371/journal.pone.0153574 27078885.
77. Mathur D, López-Rodas G, Casanova B, Marti MB. Perturbed Glucose Metabolism: Insights into Multiple Sclerosis Pathogenesis. Frontiers in Neurology. 2014;5(250). doi: 10.3389/fneur.2014.00250 25520698
78. Nguyen TT, Oh SS, Weaver D, Lewandowska A, Maxfield D, Schuler M-H, et al. Loss of Miro1-directed mitochondrial movement results in a novel murine model for neuron disease. Proceedings of the National Academy of Sciences. 2014;111(35):E3631–E40. doi: 10.1073/pnas.1402449111 25136135
79. Koch MW, George S, Wall W, Wee Yong V, Metz LM. Serum NSE level and disability progression in multiple sclerosis. Journal of the neurological sciences. 2015;350(1–2):46–50. Epub 2015/02/18. doi: 10.1016/j.jns.2015.02.009 25686504.
Článok vyšiel v časopise
PLOS One
2019 Číslo 10
- Metamizol jako analgetikum první volby: kdy, pro koho, jak a proč?
- Parazitičtí červi v terapii Crohnovy choroby a dalších zánětlivých autoimunitních onemocnění
- Nejasný stín na plicích – kazuistika
- Masturbační chování žen v ČR − dotazníková studie
- Eliquis (apixaban) nově hrazen ze zdravotního pojištění
Najčítanejšie v tomto čísle
- Correction: Low dose naltrexone: Effects on medication in rheumatoid and seropositive arthritis. A nationwide register-based controlled quasi-experimental before-after study
- Combining CDK4/6 inhibitors ribociclib and palbociclib with cytotoxic agents does not enhance cytotoxicity
- Risk factors associated with IgA vasculitis with nephritis (Henoch–Schönlein purpura nephritis) progressing to unfavorable outcomes: A meta-analysis
- Prevalence of pectus excavatum (PE), pectus carinatum (PC), tracheal hypoplasia, thoracic spine deformities and lateral heart displacement in thoracic radiographs of screw-tailed brachycephalic dogs