Manganese exposure in juvenile C57BL/6 mice increases glial inflammatory responses in the substantia nigra following infection with H1N1 influenza virus

About the Authors:

Collin M. Bantle

Roles Conceptualization, Formal analysis, Investigation, Writing – original draft

Affiliation: Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado, United States of America

C. Tenley French

Roles Data curation, Formal analysis, Writing – original draft

Affiliation: Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado, United States of America

Jason E. Cummings

Roles Data curation, Formal analysis, Writing – original draft

Affiliation: Department of Microbiology, Immunology and Pathology, Colorado State University, Fort Collins, Colorado, United States of America

ORCID logo https://orcid.org/0000-0001-6137-9163

Shankar Sadasivan

Roles Formal analysis, Methodology

Affiliation: Department of Neuroscience, Vickie & Jack Farber Institute for Neuroscience, Thomas Jefferson University, Philadelphia, Pennsylvania, United States of America

Kevin Tran

Roles Formal analysis, Investigation

Affiliation: Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado, United States of America

Richard A. Slayden

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft

Affiliation: Department of Microbiology, Immunology and Pathology, Colorado State University, Fort Collins, Colorado, United States of America

Richard J. Smeyne

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Writing – original draft

Affiliation: Department of Neuroscience, Vickie & Jack Farber Institute for Neuroscience, Thomas Jefferson University, Philadelphia, Pennsylvania, United States of America

Ronald B. Tjalkens

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado, United States of America

ORCID logo https://orcid.org/0000-0002-6848-4115

Introduction

Parkinson’s disease (PD) is characterized by the loss of voluntary motor control due to the degeneration of dopaminergic neurons in the substantia nigra pars compacta (SNpc) with associated α-synuclein protein-aggregation, neuroinflammatory activation of glial cells, mitochondrial dysfunction and oxidative stress [1]. Although genome-wide association studies (GWAS) have identified genetic variants in familial forms of the disease, these represent a very small percentage of individuals with PD, with the majority of PD cases thought to be sporadic or of an unknown etiology [2]. Epidemiological and experimental evidence suggests that environmental neurotoxin exposure and viral infections are possible risk factors for sporadic PD. Following the 1918 “Spanish Flu” pandemic, nearly every patient who had an acute episode of encephalitis lethargica (EL) from the H1N1 infection went on to develop postencephalitic parkinsonism and individuals who were born between 1888 and 1924 had a two to three-fold higher risk of developing Parkinson’s disease later in life than those born outside of that range [3–6].

Our lab and others have recently shown that exposure to certain classes of enveloped RNA viruses, Western equine encephalitis virus (WEEV) and H5N1 (strain, A/VN/1203/04) via intranasal infection can induce loss of dopaminergic neurons in the SNpc [7, 8]. Infection with viruses such as H5N1 avian influenza virus, WEEV and H1N1 induce neuronal loss in part through the activation of microglia and astrocytes and subsequent release of glial-derived neurotoxic inflammatory mediators [7–9]. Microglia and astrocytes express a variety of damage-associated molecular patterns (DAMPs) that promote inflammation and disease progression in postencephalitic parkinsonism [10]. Reactive microglia and astrocytes have increasingly become the focus of studies examining the pathophysiology of PD, suggesting that neuroinflammation may be a link between viral encephalitis and the development of parkinsonian neurological symptoms [11]. Additionally, astrocytes and microglia have innate immunological memory in the brain to facilitate a rapid inflammatory response to recurrent inflammatory stressors, and it has been postulated that this acute and exacerbated inflammatory response from glia may have the capacity to exacerbate neuronal injury following secondary insults [12].

How encephalitic infections from non-neurotropic H1N1 virus cause neurological dysfunction and neurodegeneration in certain individuals is not entirely clear but may be related to the severity of the neuroinflammatory response [7, 8, 13–17]. This could be due to prior exposure to environmental neurotoxins that activate innate immune inflammatory signaling in microglia and/or astrocytes, thereby sensitizing the tissue environment of the brain to greater inflammatory activation of glial cells during an encephalitic infection. Data from recent studies examining innate immune memory in microglia support this possibility, where multiple injections of LPS in mice resulted in immune training that amplified activation of microglia and astrocytes as well as levels of inflammatory cytokines upon later immunological challenge [12]. Similarly, mice infected with H1N1 prior to treatment with the neurotoxin, 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) had significantly greater activation of microglia and loss of dopamine neurons in the SNpc that mice treated only with MPTP [14]. Elevated levels of Mn during juvenile development in mice can exacerbate neuroinflammatory activation of glia and increase the production of reactive oxygen and nitrogen species and inflammatory cytokines following a subsequent exposure to Mn [18]. Rats exposed to Mn during juvenile development at concentrations comparable to low environmental exposure levels did not have overt neurotoxicity but showed evidence of deficits in mitochondrial respiration, as well as oxidative stress and chronic neuroinflammation [19]. Both Mn and viral infection induce inflammatory activation of microglia and astrocytes and mitochondrial stress in PD-relevant brain regions, suggesting that neuroinflammatory responses to Mn and encephalitic infections may evoke common pathophysiological signaling mechanisms in glial cells leading to altered innate immune responses in the brain [18, 20, 21].

Previous work suggests that glial inflammatory responses to Mn can exacerbate neuronal injury following secondary neurotoxic insults [10, 18, 22–25] but whether Mn can also enhance the severity of neurological damage from encephalitic viral infections with H1N1 is unknown. In the current study we investigated whether elevated levels of an Mn during juvenile development could enhance neuroinflammatory damage to dopaminergic neurons after infection with H1N1 influenza virus and increase one’s risk of developing neurological disease later in life. We tested this hypothesis by exposing juvenile C57BL/6 mice to MnCl2 in drinking water (50 mg/kg/day) for 30 days from days 21–51 PN, followed by intranasal infection with H1N1 at PN72. Control mice received only drinking water followed by either mock infection or infection with H1N1. Stereological counts of dopaminergic neurons and microglia in the SNpc were performed. We noted pronounced microglia activation following dual treatment with Mn and H1N1 relative to either treatment alone, as well as marked astrogliosis and increased in the number of reactive A1 astrocytes. RNA sequencing (RNAseq) analysis revealed activation of multiple stress response pathways involved in antioxidant activity, mitophagy, anti-viral activity and neurodegeneration in mice treated with Mn and subsequently exposed to H1N1. Collectively, these findings suggest that exposure to elevated levels of Mn during juvenile development increases neuroinflammatory activation of glia following encephalitic infection with H1N1 influenza virus later in life, likely through epigenetic modification of histones that increases secondary innate immune responses in microglia and astrocytes.

Materials and methods

Exposure protocol with manganese and H1N1

All procedures were approved by Colorado State University and St. Jude Children’s hospital Institutional Animal Care and Use Committee (IACUC) and were conducted in compliance of the National Institute of Health guidelines. Dosing was performed as previously published [26]. C57Bl/6 mice were obtained from the Jackson Laboratory and housed in a temperature-controlled room (maintained at 22–24°C on a 12 hr light/dark cycle) with ad libitum access to standard chow. At day P21, male and female C57BL/6 mice were administered MnCl2 (50mg/kg/day; Sigma) or normal drinking water. The dose of MnCl2 was calculated by monitoring water intake and weight gain for thirty days, with the concentration in drinking water adjusted to deliver 50 mg/Kg/day according to water consumption. At P51, Mn-treated mice were placed back on regular drinking water for a period of one month. Mice were then intranasally infected with A/California/04/2009 (CA/09) H1N1 or mock-infected with saline. Infection with H1N1 was performed as described previously [27]. Briefly, for infections, mice were lightly anesthetized with isofluorane and intranasally inoculated with either 103 TCID50 of CA/09 in 25μl of phosphate-buffered saline (PBS) or PBS alone and monitored daily for 21 days post-infection by assessing any possible neurobehavioral abnormalities or clinical signs of illness. A clinical scoring system was used to identify any animals with untoward morbidity, the treatments performed did not produce any overt morbidity and all animals maintained normal body weight and feeding behavior, relative to untreated/mock-infected controls. Clinical scoring was performed by laboratory staff and was supported by daily observation from dedicated laboratory animal veterinary personnel. After 21 days, infected or mock-infected control mice were euthanized for tissue collection. Euthanasia was performed under deep isofluorane anesthesia.

Preparation of A/California/04/2009 (CA/09) H1N1 inoculum

A detailed description of inoculum preparation and administration has been previously reported by our group [27]. In brief, A/California/04/2009 (CA/09) H1N1 virus was passaged in the allantoic cavity of 10-day-old specific pathogen-free embryonated chicken eggs. At 48 to 72 hours post-infection, allantoic fluid was harvested, clarified by centrifugation, and stored at −70°C. Tissue culture infectious dose 50% (TCID50) titers were determined using Madin-Darby canine kidney (MDCK) cells and evaluated by the method of Reed and Muench [28].

Tissue processing for immunohistochemistry and immunofluorescence

Mice were anesthetized with Avertin and transcardially perfused with 4% paraformaldehyde in PBS. Brains were dissected and processed for paraffin embedding. Brains were then sectioned on the microtome at 10μm thickness and mounted on polyionic slides (Superfrost-plus, Fisher Scientific). Deparaffinized SN sections were incubated with primary antibody for identification of dopaminergic neurons [mouse monoclonal anti-tyrosine hydroxylase (TH; Sigma-Aldrich; 1:500), microglia [rabbit polyclonal anti-IBA1 (Wako Chemicals; 1:500), astrocytes [rabbit polyclonal anti-S100beta (Abcam; 1:500)], complement C3 [rat monoclonal anti-C3 (Abcam; 1:100)], SerpinA3 [mouse monoclonal anti-SerpinA3 (thermos; 1:250)], IP-10 [mouse monoclonal anti-IP-10 (Santa Cruz; 1:500)], IP-10/Cxcl10 [mouse monoclonal anti-cxcl10 (Santa Cruz; 1:250)], CCL2 [mouse monoclonal anti-Ccl2 (Millipore; 1:500)], and acetylated Lysine residues [rabbit polyclonal anti-acetylated lysine (Cell Signaling; 1:500)]. For immunohistochemical analysis, the secondary antibodies included biotinylated mouse IgG (for TH, 1:1000) or biotinylated rabbit IgG (for IBA1, 1:1000). Diaminobenzidine (DAB) or a VIP kit (Vector labs) reaction was used to yield a brown (TH) or a purple (IBA1) color, respectively. For immunofluorescence, anti-mouse, anti-rabbit or anti-rat IgG alexa flour 555, alexa flour 488, or alexa flour 647 were diluted in TBS (2% Triton) at 1:500. Sections were washed 5X (5 min) and stained with DAPI in the final wash after an hour incubation period in the secondary antibody. Sections were then mounted with medium, coverslipped and stored at 4°C until imaged.

Quantification of TH+ dopaminergic neurons and Iba1+ microglia in the substantia nigra pars compacta

Quantitation of neurons and glial cells was performed as previously reported [27]. In brief, TH+ dopaminergic neurons and IBA1+ microglia in the SNpc were estimated using standard model-based stereological methods [16, 29]. Counts of total dopaminergic neurons and activated microglia were estimated using Microbrightfield StereoInvestigator (MBF Biosciences, Williston, VT) and the optical fractionator method using an Olympus BX-51 microscope and 100X objective [30, 31]. The identification of resting and activated microglia was based on defined morphological criteria, as previously reported [16, 30]. Resting microglia were defined as having a small, oval IBA1+ cell body that averaged 3 microns in diameter with long slender processes, while microglia were classified as activated when the cell body was slightly increased in size compared to resting microglia and had an irregular shape, with shorter and thickener processes. The investigator was blinded from all experimental groups during imaging and cell quantitation.

Quantification of astrocyte-specific inflammatory markers in the substantia nigra pars compacta

Formalin-fixed, paraffin-embedded 10 μm brain sections were immunofluorescently-labeled using a Leica Bond RXM automated robotic staining system. Sections were immunohistochemically stained on a Leica Bond-III IHC automated stainer. Antigen retrieval was performed with Bond Epitope Retrieval Solution 2 for 20 minutes. Sections were then incubated with primary antibodies for S100β+ (Abcam; rabbit, 1:500), complement C3 (Abcam; rat monoclonal, 1:100), SerpinA3 (Thermos; mouse monoclonal, 1:250), IP-10/CXCL10 (Santa Cruz; mouse monoclonal, 1:250) and CCL2 (Millipore; mouse monoclonal, 1:500). Secondary antibodies included AlexaFluor anti-rabbit IgG AlexaFluor 488, anti-rat IgG AlexaFluor 555 and anti-mouse AlexaFluor 647. Whole-brain immunofluorescence montage images of labeled tissue sections were imaged using an automated Olympus BX51 fluorescence microscope equipped with a Hamamatsu ORCA-flash 4.0 LT CCD camera and collected using Olympus Cellsens software (v 1.15). Quantitative analysis was performed on dual- or triple-labeled fluorescent images generated by montage imaging of an entire coronal mouse brain section compiled from individual images acquired using an Olympus Plan Apochromat 20X air objective (0.40 N.A.). All slides were scanned under the same conditions for acquisition time, magnification, exposure time, lamp intensity and camera gain. The substantia nigra was delineated by neuroanatomical landmarks and referenced to the Allen brain atlas, following application of an adaptive threshold with shape factor and area (μm2) object filters for automatic S100β+ astrocyte cell detection. The number of cells was divided over the area (μm2) of the region. To measure expression of inflammatory proteins within S100β+ astrocytes, mean intensities of complement C3, serpinA3, IP-10/CXCL10 and CCL2 were measured by generating automated individual ROIs around all S100β+ astrocytes within the SNpc. We assessed the presence of the inflammatory molecules on two coronal sections per animal that were 10 μm in thickness, spaced at 200 μm intervals within the SNpc, with an N = 6–8 mice for each treatment group. The investigator was blinded from all experimental groups during imaging and cell quantitation.

Quantification of histone acetylation in dopaminergic neurons, microglia and astrocytes in the substantia nigra par compacta

Formalin-fixed, paraffin-embedded 10 μm brain sections were immunofluorescently-labeled as described above using a Leica Bond RXM automated robotic staining system. Sections were then incubated with primary antibodies for tyrosine hydroxylase (WAKO; goat, 1:500), S100β (Abcam; rabbit, 1:500), IBA1 (Abcam; goat, 1:100) and lycine-acetylated histones (Cell Signaling; rabbit, 1:500). Secondary antibodies included anti-goat IgG AlexaFluor 647, anti-mouse IgG AlexaFluor 488, and anti-rabbit IgG AlexaFluor 555. To measure total lysine acetylation in dopaminergic neurons, astrocytes and microglia in the SNpc, labeled tissue sections were imaged using a Olympus Plan Apochromat 20X air objective (0.40 N.A.) and an automated Olympus BX51 fluorescence microscope equipped with a Hamamatsu ORCA-flash 4.0 LT CCD camera and collected using Olympus Cellsens software (v 1.15). Quantitative analysis was performed as described above based on triple-labeled fluorescent images montage images of an entire coronal mouse brain section. All slides were scanned under the same conditions for acquisition time, magnification, exposure time, lamp intensity and camera gain. The substantia nigra was delineated by neuroanatomical landmarks and referenced to the Allen brain atlas, following application of an adaptive threshold with shape factor and area (μm2) object filters for automatic TH+ dopaminergic neurons, IBA1+ microglia, and S100β+ astrocyte cell detection. Mean intensities of acetylated lysine residues were measured within TH+ dopaminergic neurons, IBA1+ microglia, and S100β+ astrocytes in the SNpc by generating automated ROIs around each cell type. We assessed total lysine histone acetylation in dopaminergic neurons, astrocytes and microglia in two coronal sections per animal that were 10 μm in thickness, spaced at 200 μm intervals within the SNpc, with an N = 6–8 mice for each treatment group. The investigator was blinded from all experimental groups during imaging and cell quantitation.

RNA sequencing of the substantia nigra

This method was adapted from a previous protocol [32]. A systems-based transcriptional analysis of substantia nigra brain tissues from each treatment group was performed, and the treatment conditions included control, H1N1, and Mn+H1N1. Tissue was anatomically dissected from the basal midbrain and included the substantia nigra but not the ventral tegmental area. Isolated tissue was immediately flash frozen in liquid nitrogen and later homogenized in Trizol reagent (Thermo Fisher) for purification. Samples were then treated with DNAse (Fermentas, Burlington, Ontario) for 30 minutes and purified by phenol/chloroform/isoamyl alcohol (25:24:1) (Fisher Scientific, Pitts- burgh, PA) extraction and ammonium acetate precipitation. Quality and integrity of total RNA was assessed using the 4200 Agilent Tapestation, and samples were confirmed to have RIN scores of >7. RNA [transcripts & non-coding RNA] was isolated from total host RNA followed by library construction and template preparation with the Ion Total RNA-Seq kit and Ion Chef system kit. Sample libraries were prepared using the Ion Total RNA-Seq kit v2 (Life Technologies) and multiplexed on a P1 chip using Ionxpress RNA-Seq 1–16 kit (Life Technologies). Whole mouse transcriptome sequencing was then performed using the Ion Proton Next Generation Sequencer (Life Technologies) through the core facility at the Infectious Disease Research Complex (IDRC) at Colorado State University. Following Next Generation RNA-seq, we used read count coverage (RPKM or FPKM values) to compare the differential gene expression between groups. Advanced RNA-seq analysis was done using the Tuxedo package in Linux command line, including alignments using Bowtie2 and differential gene expression analysis using Cufflinks. Local realignment and base quality score recalibration (BQSR) methods were used as needed to reduce false-positive base calls and improve alignments.

Analysis of next-generation sequencing data

Data analysis was performed as previously published [32]. In brief, FASTQ files were analyzed using Galaxy for quality trimming, with minimum PHRED quality threshold set at 20 and all read length greater than 20bp. Trimmed reads were then aligned to Mus musculus mm9 using Bowtie2 and gene expression determined using Cufflinks. Expression output was normalized in FPKM format (fragments per kilobase of exon per million reads). Replicate mean values were calculated, and the data was further reduced to FPKM values greater than two. Venn diagrams and Pie chats were generated by comparing the reduced FPKM transcript totals for each treatment group. PANTHER (Protein Analysis Through Evolutionary Relations) Classification System (http://pantherdb.org/publications.jsp#HowToCitePANTHER) was used to ascertain functional pathways driving differences in gene expression by analyzing the complete gene lists within each treatment group for their corresponding annotations, accessed from the Gene Ontology (GO) Consortium. Secondarily, the PANTHER statistical overrepresentation test enabled a comparison of those GO annotations across each treatment group to identify functionally related genes relevant to neurodegeneration. The PANTHER gene list analysis was used to perform a functional classification of all the GO-annotated transcripts within each treatment group. This process employed the GO Term of Molecular Function and GO-slim annotation data sets to analyze the expression profiles of each gene list. We charted any differences between treatment groups as the number of gene hits (the % of gene hits for a GO term /total number of annotated genes in that category). The only GO Term categories shown are those with demonstrable differences in the number of gene hits between the treatment groups. The exact genes returned within each GO classification were analyzed in Excel using MATCH function and dual comparisons to find the specific transcripts unique to each treatment category. To further classify RNAseq transcripts according to function, gene lists were analyzed using the PANTHER overrepresentation test. This yielded statistically over or under-represented annotations among gene lists relative to the GO Biological Process and Molecular Function Mus musculus datasets. Complete gene lists for each treatment group were imported into the analysis tool. Fisher’s Exact test for significance was used with FDR multiple test correction. Resulting GO Terms with an enrichment score >1.5 and a FDR<0.05 were considered significant. Genes relevant to neurodegeneration were selected from significant GO Term gene lists. The Search Tool for the Retrieval of Interacting Genes database (STRING) was used to create a network diagram of functional associations between protein products of unique genes within the H1N1+Mn treatment group. Each protein-protein association is weighted according to evidence channels delineated by differing colors. The confidence cutoff was set to 0.4.

Statistical analysis

All data was presented as mean +/- SEM, unless otherwise noted. Experimental values from each mean were analyzed with a Grubb’s (α = 0.05) test for exclusion of significant outliers. Differences between each experimental group were analyzed by a one-way ANOVA following a Tukey post hoc multiple comparisons test. Significance was identified as ++P < 0.01, *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. All statistical analysis was conducted using Prism (version 6.0; Graph Pad Software, San Diego, CA).

Results

Juvenile manganese exposure increases the number of activated microglial following encephalitic infection with H1N1 in the SNpc

To assess if Mn exposure during juvenile development would enhance the neurological effects of infection with H1N1 in the substania nigra, three week old C57BL/6 mice were administered MnCl2 (50mg/kg/day) or normal drinking water for a total of 30 days and then intranasally infected with H1N1 or mock-infected with saline at 3 months of age (Fig 1A). Stereological determination of TH+ dopaminergic neurons and morphological analysis of Iba1+ microglia at 21 DPI revealed that pre-treatment with MnCl2 during juvenile development induced persistent morphological changes in microglia consistent with an activated phenotype and increased their reactivity to a subsequent infection with H1N1, characterized by retraction of cytoplasmic processes and adoption of an amoeboid phenotype (Fig 1B–1I). We did not observe significant changes in the number of dopaminergic neurons, ⍺-synuclein protein aggregation, or the total number of resting microglia in the SNpc at 21 DPI (Fig 1F and 1G).

[Figure omitted. See PDF.]

Fig 1. Pre-treatment with manganese increases microglial activation in the substantia nigra following infection with H1N1.

(A) C57BL/6 mice were divided into groups with and without Mn in drinking water (50 mg/Kg) from day 21–51 PN and then exposed to H1N1 influenza virus three weeks later. Groups: Control, Mn, H1N1, H1N1+MnCl2. (B-E) IHC labeling of dopamine neurons (TH) and microglia (IBA1). (F) Stereological determination of the number of TH+ neurons in the substantia nigra pars compacta (SNpc). (G) Resting microglia in the SNpc. (H) Activated microglia in the SNpc. (I) Ration of activated/resting microglia in the SNpc. *P<0.05 **P<0.01 ***P<0.001****P<0.0001. n = 6 mice/group.

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

Juvenile manganese exposure increases neuroinflammatory activation of astrocytes following encephalitic infection with H1N1 in the SNpc

Given the increased number of reactive microglia in the substantia nigra and previous work showing that astrocytes play a significant role in microglial activation through glial-glial communication [7, 22, 23, 33, 34], we examined the extent and severity of astrocyte activation in the basal ganglia following treatment with MnCl2 (50mg/kg/day) and intranasal infection with H1N1 at 21 DPI (Fig 2). To determine if dual treatment with Mn enhanced the inflammatory phenotype of astrocytes and increased the number of A1 neurotoxic astrocytes in the SNpc, we measured the level of astrogliosis (Fig 2A–2E), as well as the expression of A1-specific astrocyte inflammatory markers (C3, SerpinA3, IP10, CCL2) in S100β+ astrocytes with immunofluorescence co-localization (Fig 2F–2Y). At 21 DPI, there was a significant increase in S100β+ astrocytes in the SNpc with dual treatment compared to control and Mn treatment alone (Fig 2E). Additionally, complement C3 and SerpinA3 were significantly increased following pre-treatment with Mn and H1N1 infection compared to control, Mn and H1N1 treatment alone (Fig 2J–2Y). Dual treatment did not induce any differences in expression of Interferon gamma-induced protein 10 (IP-10) or monocyte chemotactic protein (CCL2) (Fig 2P–2Y).

[Figure omitted. See PDF.]

Fig 2. Pre-treatment with manganese induces proliferation and inflammatory activation in astrocytes in the substantia nigra following infection with H1N1.

(A-D) IF labeling of astrocytes (s100β). Groups: Control, Mn, H1N1, H1N1+ MnCl2. (E) Cell counts of the s100β+ astrocytes in the substantia nigra pars compacta (SNpc). (F-I) IF colocalization of astrocytes (s100β) with complement C3 (C3). (J) C3 mean intensity measurements in s100β+ astrocytes in the SNpc. (K-N) IF colocalization of astrocytes (s100β) with SerpinA3. (O) SerpinA3 mean intensity measurements in S100β+ astrocytes in the SNpc. (P-S) IF colocalization of astrocytes (s100β) with complement IP-10/Cxcl10 (IP-10). (T) IP-10 mean intensity measurements in s100β+ astrocytes in the SNpc. (U-X) IF colocalization of astrocytes (s100β) with complement Ccl2. (J) Ccl2 mean intensity measurements in s100β+ astrocytes in the SNpc. *P<0.05 **P<0.01 ***P<0.001****P<0.0001. n = 6 mice/group.

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

Dual treatment with Mn and H1N1 alters of histone acetylation in the substantia nigra

To examine the basis for the heightened innate immune inflammatory response observed in glial cells in the basal ganglia following the two-hit exposure model with Mn and H1N1, we assessed histone acetylation in glia and neurons in the substantia nigra. To quantitate the level of cell type-specific histone acetylation in the basal midbrain, brain sections were stained for total acetylated histone lysine residues in TH+ dopaminergic neurons (Fig 3A–3E), S100β+ astrocytes (Fig 3F–3J), and IBA1+ microglia in the SNpc (Fig 3K–3O). Infection with H1N1 alone or following Mn pre-treatment significantly decreased histone lysine acetylation in dopaminergic neurons in the SNpc (Fig 3A–3E). Minimal differences were noted in histone acetylation in astrocytes (Fig 3F–3J) In contrast, mice infected with H1N1 alone or following pre-treatment with Mn had significantly increased histone acetylation in IBA1+ microglia in the SNpc (Fig 3F–3O), consistent with previous findings in human PD brains [35].

[Figure omitted. See PDF.]

Fig 3. Dual treatment with MnCl2 and H1N1 alters histone acetylation in dopaminergic neurons, microglia, and astrocytes in the SNpc.

(A-D) Immunofluorescence co-localization of dopaminergic neurons (TH) and total acetylated lysine residues (AL). Groups: Control, Mn, H1N1, H1N1+ MnCl2. (E) AL mean intensity measurements in TH+ dopaminergic neurons in the substantia nigra pars compacta (SNpc). (F-G) IF colocalization of astrocytes (s100β) with AL. (J) AL mean intensity measurements in s100β+ astrocytes in SNpc. (K-N) IF colocalization of microglia (Iba1) with AL. (O) AL mean intensity measurements in Iba1+ microglia in the SNpc. *P<0.05 **P<0.01 ***P<0.001****P<0.0001. n = 6 mice/group.

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

Dual treatment with Mn and H1N1 induces a unique transcriptional signature in the substantia nigra consistent with a neurodegenerative phenotype

Given the differences in patterns of histone acetylation in the SNpc following treatment with Mn and H1N1, we performed Next Generation RNA sequencing (RNAseq) of brain tissue from the SN to assess global transcriptional patterns in control mice and those infected with H1N1 as adults with and without Mn pre-treatment during juvenile development (Fig 3). Given our previous findings showing minimal neuropathology in adult mice following treatment with Mn during juvenile development [18, 24, 25], we directly compared H1N1 infected mice to those infected with H1N1 + Mn in the RNA-seq analysis, with both groups referenced to untreated control mice receiving only mock viral infection. Thus, the strategy was to determine how pretreatment with Mn would modulate the severity of infection with H1N1 with respect to changes in gene expression in the SN. RNA was isolated from the SN at 21 DPI for sequencing. The resulting FASTQ files were analyzed using Galaxy for quality trimming, with minimum PHRED quality threshold set at 20 and all read length greater than 20 bp. Trimmed reads were then aligned to the mouse genome using Bowtie2 and gene expression determined using Cufflinks. Local realignment and base quality score recalibration (BQSR) methods were used as needed to reduce false-positive base calls and improve alignments. The expression output was normalized in FPKM format. The replicate mean values were calculated, and the data was further reduced to FPKM values greater than two.

This unbiased global analysis of the transcriptional profile between Control, H1N1 and H1N1+Mn treatment groups revealed that of the 13,168 transcripts annotated from the SN, 351 where unique to control, 290 were unique to H1N1 and 951 transcripts were unique to the dual treatment group (Fig 4A). The complete list of annotated transcripts is provided in S1 Table. To assess whether the increased inflammatory phenotype of microglia and astrocytes in the dual treatment group was a gain of function or loss of function, we next determined the major biological pathways that were altered in each treatment group by using the Gene Ontology (GO) Consortium and PANTHER Classification System pathway and overrepresentation analyses. Resulting transcript annotations were quantitatively different between treatment groups in the following four Molecular Function GO Term categories: antioxidant activity (GO:0016209) (Fig 4A and 4B), catalytic hydrolase activity (GO:0016788,GO:0016462,GO:0003924) (Fig 4C), catalytic transferase activity (GO:0008168,GO:0016757,GO:0016301) (Fig 4D), and G-protein coupled receptor activity (GO:0004930) (Fig 4E). The exact genes returned within each GO Term classification were further analyzed to find the specific transcripts unique to each treatment type. Unique transcripts are listed in Fig 4B–4E. Of note, the dual-treatment group active transcript list includes Interferon-inducible GTPase1, Dynein heavy chain, homeobox protein, mitochondrial hydrodoxymethylglutaryl-CoA synthase (Fig 4B–4E), as well as DJ-1/Park 7, other interferon regulatory proteins, autophagy related proteins (Atg), amyloid-beta precursor bindings proteins, NFκB related inflammatory proteins and histone acetyltransferases (S2 Table). Additionally, control and H1N1 treated mice uniquely expressed glutathione peroxidase 8, Oasl3, glutathione-s-transferase, thioredoxin and colony-stimulating factor receptor 1 (CSF1R) (Fig 4B–4E), suggesting that specific stress response pathways are activated and protective mechanisms are lost following dual treatment with Mn and H1N1 when compared to Control and H1N1 treatment alone.

[Figure omitted. See PDF.]

Fig 4. RNA sequencing of the H1N1+MnCl2 treatment group highlights uniquely transcribed genes annotated to functionally relevant biological processes.

(A) Venn diagram depicting unique and overlapping transcripts. (B-E) GO Consortium & PANTHER Classification System pathway analysis: Pie charts describe the relative proportions of RNAseq transcripts annotated to each GO Term within the molecular function Gene Ontology. The exact gene lists unique or uniquely absent from the H1N1+ Mn treatment group are listed for each category. Genes were classified according to molecular function and biological process with Panther overrepresentation and Functional Enrichment analyses. Groups: Control, H1N1, H1N1+ Mn.

https://doi.org/10.1371/journal.pone.0245171.g004

GO PANTHER analyses were used to further classify and compare transcript lists according to functional pathways. To understand system-level changes in cellular processes induced by pre-exposure to Mn, we interrogated the protein–protein interaction networks in the dual treatment group using the Search Tool for the Retrieval of Interacting Genes database (STRING) (Fig 5A) [36, 37], wherein network nodes represent proteins and lines represent functional associations between those proteins. The color of each line represents the origin and/or type of evidence supporting that protein-protein interaction and the weights of each line correspond to a confidence score for that evidence type (Fig 5A). From a functional perspective, an association can mean direct physical binding, but STRING evaluates each protein-protein interaction according to evidence from seven different classification channels: neighborhood, co-occurrence, co-expression, experiments, textmining, database and fusion. The neighborhood category is related by conserved and co-transcribed operons; the co-ocurrence category is related by phylogenetic distribution of orthologs of all proteins in a given organism; the co-expression category is related by predicted association between genes based on observed patterns of simultaneous expression of genes; the experiments category is related by known experimental interactions; the textmining category conducts statistical co-citation analysis across a large number of scientific texts, including all PubMed abstracts and OMIM; the database category is expertly curated and imported from pathway databases; and the fusion category is related by proteins that are fused in some genomes and are most likely to be functionally linked [2, 36, 37]. The more lines shown between each protein-protein interaction represents a more likely biological interaction. The most highly represented interactions in the current study were centralized around Lingo 2, Pak6, Tbpl1 and Ctnnb1 (Fig 5A). Analyzing the RNAseq data using STRING methods, we identified a list of genes unique to the H1N1+Mn treatment group that are relevant to neurodegeneration (Table 1). Notable pathways include those contributing to multiple stress responses such as neuroinflammation, oxidative stress, protein misfolding and neurodegeneration. These included anti-oxidant genes such as Kelch-like ECH-associated protein 1 (Keap1), the familial PD gene, Park7, and mitogen activated protein kinase 4k4 (Map4k4), which is involved in the response to environmental stressors and inflammatory cytokines. Other genes at highly interactive nodes unique to the H1N1 + Mn group included Sox11 and 13, which are involved in neurogenesis, and the eukaryotic translation initiation factor 2C1 (Eif2c1), which is involved in protein synthesis and stress responses through post-trasncriptional silencing and damage repair pathways.

[Figure omitted. See PDF.]

Fig 5. Transcriptional connectome of dual treatment group with associated gene products.

(A) Search Tool for the Retrieval of Interacting Genes database (STRING) diagram of functional associations between protein products of unique genes within the H1N1+Mn treatment group. Nodes represent proteins and edges correspond to functional interactions. Edge colors differ according to a scored confidence scale based on the extent and type of evidence supporting that particular association. Evidence channels: Green Line: Neighborhood, Navy Line: Co-occurrence, Black line: Co-expression, Pink line: Experiments, Yellow line: Textmining, Light blue line: Databases, Red line: Fusion.

https://doi.org/10.1371/journal.pone.0245171.g005

[Figure omitted. See PDF.]

Table 1. H1N1+Mn unique genes relevant to neurodegeneration.

https://doi.org/10.1371/journal.pone.0245171.t001

Discussion

Although the majority of PD cases are sporadic and of unknown etiology, there is increasing evidence that environmental stressors such as neurotoxic chemicals and encephalitic infections can increase susceptibility to this and other neurological diseases [4, 38, 39]. Chronic inflammatory activation of glia in the nigro-striatal pathway is a well-established feature of both sporadic PD and post-encephalitic parkinsonism and could represent a cellular mechanism linking environmental exposures early to potentiation of neurodegeneration following subsequent insults. Using a two-hit model with H1N1 and MPTP, it was recently shown that influenza infections can enhance innate immune responses of microglia in the SNpc following exposure to the parkinsonian agent, MPTP [16]. It is not known whether chronic exposure to environmental metals such as Mn have a similar capacity to enhance neurological injury from viral infection. However, earlier studies in mice reported that a single dose of Mn or Cd enhanced the neurovirulence of infection with several alphaviruses and increased the severity of symptoms, neuroinflammation and mortality [40]. To test if exposure to Mn during juvenile development could similarly prime glial cells in the SNpc for a more severe neuroinflammatory response following encephalitic infection with H1N1, we utilized a two-hit exposure model with Mn and H1N1 and examined whether there were epigenetic and transcriptomic changes in the SNpc that could explain the heightened innate immune response to viral infection in microglia and astrocytes.

Our results suggest that Mn exposure during juvenile development induces lasting neuroinflammatory and epigenetic alterations in glia that exacerbate the neuroinflammatory response to H1N1 compared to mice that did not receive excess Mn as juveniles. Although we did not observe overt loss of dopaminergic neurons or protein aggregation of ⍺-synuclein at 21 DPI in any treatment group (Fig 1), there was a marked increase in activation of microglia that correlated with a downward trend in the number of dopaminergic neurons that would likely accelerate with aging, thus predisposing to neurological dysfunction. Therefore, these pathological findings may reflect early neuropathological changes corresponding the prodromal stage of PD or viral parkinsonism. This highlights the importance of neurotoxic environmental exposures that promote a reactive inflammatory phenotype in microglia that could predispose to injury within the nigro-striatal dopamine system. Recent studies support a critical role for Mn in modulating innate immunity in response to viral infection through enhancing the sensitivity of pattern recognition receptors that then stimulate anti-viral gene expression [41]. Release of Mn into the cytosol during infection with dsDNA viruses increases activation of the cGAS-STING pathway by elevating production of cGAMP that stimulates NFκB-dependent expression of the anti-viral interferon response. This is consistent with previous data from our laboratory demonstrating that Mn directly stimulates soluble guanylate cyclase and thereby enhances NFκB-induced expression of inflammatory genes through increases in cGMP [25].

Chronic neuroinflammatory activation of microglia and astrocytes is a central feature of aging, viral encephalitis and neurodegenerative disease. Increasing evidence implicates microglia-astrocyte signaling in neuroinflammatory responses that enhance neurodegeneration [42–45]. Our lab and others have demonstrated that microgliosis precedes astrogliosis and neuronal loss in models of PD [22, 46, 47], indicating that cytokine and chemokine signaling from neuroinflammatory activation of microglia likely plays a critical role in inflammatory activation of astrocytes in the SNpc. Given that dual treatment with Mn and H1N1 potentiated inflammatory activation of the A1 astrocyte markers, SerpinA3 and C3 (Fig 2), and also increased the number of reactive microglia (Fig 1), these results provide additional evidence that innate immune signaling in activated microglia plays a critical role in inflammatory activation of astrocytes. Mn directly stimulates expression of TNF in microglia that can enhance the production of inflammatory cytokines and chemokines in astrocytes [22], including CCL2, that promote neuronal injury [33]. Thus, inflammatory ‘priming’ of microglia during juvenile exposure to Mn could enhance activation of astrocytes into a neurotoxic A1 phenotype during subsequent exposure to H1N1. The data reported in the current study supports this model, with marked increases in expression of the A1 proteins SerpinA3 and C3 in reactive astrocytes in dual-treated mice (Fig 3). However, Mn pre-exposure did not increase levels of IP10 and CCL2 in astrocytes following H1N1 infection, suggesting that individual neurotoxic exposures or pathogens likely have a unique molecular signature associated with distinct inflammatory phenotypes in astrocytes.

The neurotoxic effects of Mn are mediated both through direct toxic effects in neurons, as well as through activation of inflammatory signaling pathways which further injure neurons through overproduction of reactive oxygen and nitrogen species and inflammatory cytokines [19]. Thus, the capacity of environmental exposure to sensitize neural tissue to additional damage from subsequent H1N1 encephalitic infection may involve persistent inflammatory changes in glial cells [14]. It has previously been shown that juvenile exposure to Mn can induce lasting effects on the neuroinflammatory status of astrocytes and microglia which continues into adulthood and this likely contributes increased susceptibility to secondary environmental insults and infections [18, 24, 25, 48, 49]. However, the mechanisms regulating this sustained inflammatory state after environmental insults in microglia and astrocytes are unclear. Previous studies have reported that manganese can alter histone acetylation and gene expression, chromatin remodeling, cell cycle progression, DNA repair and apoptosis in neurons and glia [12, 50]. Glial cells possess innate immune memory to environmental stimuli through histone acetylation that alters subsequent inflammatory responses [12]. Histone acetylation is also altered in PD, likely associated with microglial activation in the SNpc that increases dopaminergic neurodegeneration [35, 51, 52]. Additionally, decreased acetylation in TH neurons has been noted in clinical PD cases, and histone acetylation is thought to be modulated in glia following exposure to other environmental insults associated with PD [35, 51–53]. Here we observed that juvenile exposure to Mn followed by adult infection with H1N1 causes changes in histone acetylation in dopaminergic neurons, microglia and astrocytes in the SNpc, relative to infection with H1N1 alone (Fig 3). Also, dual treatment with Mn and H1N1 decreased histone acetylation in dopaminergic neurons, corresponding with data in clinical cases of PD [35, 51–53]. We surmise that these epigenetic changes in glia likely play a role in their increased inflammatory activation to infection with H1N1. However, whether these epigenetic changes persist throughout aging remains to be determined.

It is still unclear whether neuroinflammatory activation of glia increases susceptibility to neurodegeneration primarily through decreased release of neurotrophic factors or from excessive synthesis of neurotoxic inflammatory mediators [47, 54]. Glial activation is regulated through multiple pathways including mitogen-activated protein kinases (MAPKs), activator protein-1 (AP-1), Janus kinase (JAK)/signal transducer and activator of transcription (STAT), interferon regulatory factor families (IFN), as well as through the nuclear factor kappa B (NFκB) pathway. To address this question, we performed Next Generation RNA-sequencing to examine gene networks regulated by Mn that could predispose neurons in the substantia nigra to injury following infection with H1N1. We found that pre-treatment with Mn prior to infection with H1N1 increased the number of unique transcripts and significantly altered the global transcriptional profile in the SNpc compared to H1N1 treatment alone (Fig 4). Amongst unique transcripts in the dual treatment group, STRING analysis revealed an overrepresentation of Interferon-inducible GTPase1, dynein heavy chain, homeobox protein, mitochondrial hydrodoxymethylglutaryl-CoA synthase (Fig 4B–4E), as well as DJ-1/Park 7, other interferon regulatory proteins, autophagy-related proteins (Atg), amyloid-beta precursor binding proteins, NFκB related inflammatory proteins and histone acetyltransferases (S2 Table). Park7 (DJ1) is linked to the preservation of mitochondrial function in PD, and the unique activation of DJ-1, Atg and interferon regulatory protein in the dual treatment group likely represents the activation of the antioxidant activation and autophagy stress response pathways [55], possibly in reponse to underrepresented transcripts involved in antioxidant activity, such as glutathione peroxidase (Fig 4A). These data suggest that increased inflammatory activation of glia in the SNpc in the Mn + H1N1 group was associated with stress reponses in this brain region consistent with inflammation and innate immune function, as well as oxidative stress and mitophagy. Moreover, given that viral infections can inhibit mitochondrial bioenergetics by depressing cellular ATP content and induce oxidative stress in neurons and glia, these data suggest that juvenile exposure to Mn could reduce glutathione-dependent anti-oxidant protection through a reduction of glutathione peroxidase that exacerbates oxidative stress during subsequent infection with H1N1 [42].

Analysis of RNA sequencing data also identified unique patterns of gene expression within each treatment group. Annotated of genes to specific molecular function and biological processes (Fig 4) enabled us to generate a network summarizing predicted associations between gene products to attain a better systems-level understanding of cellular processes in the dual treatment group (Fig 5). Dj-1/Park7, Lingo2, Keap1, Eif2c1, and Pak6 were all unique genes identified using overrepresentation analysis that are relevant to etiology and/or progression of PD (Table 1). Interestingly, LINGO2 is a member of LRR gene family that, along with LRRK2, has been linked to Essential tremor (ET) and PD and has even become a promising therapeutic target in multiple sclerosis (MS) and PD [56–58]. Taken together, these findings support the involvement of distinct proteins and pathways in the neurologically deleterious effects caused by successive environmental challenge with Mn and H1N1. Although there was not significant loss of dopaminergic neurons at the relatively early timepoint evaluated, the unique protein-protein interactions in dual treated animals could represent an early pre-symptomatic stage of neuronal dysfunction corresponding to prodromal disease.

The present studies used a two-hit model of successive exposure to Mn and H1N1 influenza virus to identify mechanisms by which multiple environmental insults and microbial infection could act in concert to increase susceptibility to PD and related neurodegenerative diseases. Exposure to Mn from PN 21–51 did not result in significant loss of dopaminergic neurons in the SNpc but markedly increased neuroinflammatory activation of microglia and astrocytes that could represent a form of innate immune memory in the brain that predisposes glia cells to a neurotoxic reactive phenotype during aging. This is consistent with the patterns of gene expression noted in RNA-Seq studies that revealed transciptional signatures consistent with a neurodegenerative phenotype, particularly for genes related to oxidative stress, mitophagy, protein processing, and immune function. The unique patterns of gene expression noted in animals exposed to both Mn and H1N1 suggests that the mechanism by which multiple environmental exposures modulate neurotoxic injury seen in certain cases of sporadic and post-encephalitic PD may involve epigenetic changes that favor expression of genes associated with inflammation and protein misfolding in astrocytes and microglia.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Complete list of expressed transcripts in each experimental group.

The complete list of annotated transcripts from unbiased global analysis of the transcriptional profile of Control, H1N1 and H1N1+Mn treatment groups is presented in S1 Table.

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

(XLSX)

S2 Table. Uniquely expressed transcripts in each experimental group.

Transcripts uniquely expression in each experimental group are presented in S2 Table. Mapping these transcripts to the major biological pathways that were altered in each treatment group was performed using the Gene Ontology (GO) Consortium and PANTHER Classification System pathway and overrepresentation analyses.

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

(XLSX)

Acknowledgments

The authors acknowledge the technical assistance of animal care and use veterinary personnel at Colorado State University and Thomas Jefferson University.

Citation: Bantle CM, French CT, Cummings JE, Sadasivan S, Tran K, Slayden RA, et al. (2021) Manganese exposure in juvenile C57BL/6 mice increases glial inflammatory responses in the substantia nigra following infection with H1N1 influenza virus. PLoS ONE 16(1): e0245171. https://doi.org/10.1371/journal.pone.0245171

References

1. Beal MF. Parkinson’s disease: a model dilemma. Nature. 2010;466(7310).

2. Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: Online Mendelian Inheritance in Man (OMIM(R)), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43:D789–98. pmid:25428349

3. Reid AH, McCall S, Henry JM, Taubenberger JK. Experimenting on the past: the enigma of von Economo’s encephalitis lethargica. J Neuropathol Exp Neurol. 2001;60(7):663–70. pmid:11444794

4. Dourmashkin RR. What caused the 1918–30 epidemic of encephalitis lethargica? J R Soc Med. 1997;90(9):515–20. pmid:9370993

5. Elizan TS, Schwartz J, Yahr MD, Casals J. Antibodies against arboviruses in postencephalitic and idiopathic Parkinson’s disease. Arch Neurol. 1978;35(5):257–60 pmid:646678

6. Clarke P, Leser JS, Quick ED, Dionne KR, Beckham JD, Tyler KL. Death receptor-mediated apoptotic signaling is activated in the brain following infection with West Nile virus in the absence of a peripheral immune response. J Virol. 2014;88(2):1080–9 pmid:24198425

7. Bantle Collin M. P AT, Smeyne Richard J., Rocha Savannah M., Olson Ken E. and Tjalkens Ronald B. Infecton with mosquito-borne alphavirus induces selective loss of dopaminergic neurons, neuroinflamation and widespread protein aggregation. Nature Parkinson’s Disease 2019. Epub In press 08/08/2019.

8. Jang H, Boltz D, McClaren J, Pani AK, Smeyne M, Korff A, et al. Inflammatory effects of highly pathogenic H5N1 influenza virus infection in the CNS of mice. J Neurosci. 2012;32(5):1545–59. Epub 2012/02/04. pmid:22302798.

9. Ludlow M, Kortekaas J, Herden C, Hoffmann B, Tappe D, Trebst C, et al. Neurotropic virus infections as the cause of immediate and delayed neuropathology. Acta Neuropathol. 2016;131(2):159–84. pmid:26659576.

10. Tjalkens RB, Popichak KA, Kirkley KA. Inflammatory Activation of Microglia and Astrocytes in Manganese Neurotoxicity. Adv Neurobiol. 2017;18:159–81. Epub 2017/09/11. pmid:28889267.

11. Caggiu E, Arru G, Hosseini S, Niegowska M, Sechi G, Zarbo IR, et al. Inflammation, Infectious Triggers, and Parkinson’s Disease. Front Neurol. 2019;10:122. Epub 2019/03/07. pmid:30837941.

12. Wendeln AC, Degenhardt K, Kaurani L, Gertig M, Ulas T, Jain G, et al. Innate immune memory in the brain shapes neurological disease hallmarks. Nature. 2018;556(7701):332–8. pmid:29643512.

13. Jang H, Boltz DA, Webster RG, Smeyne RJ. Viral parkinsonism. Biochim Biophys Acta. 2009;1792(7):714–21. pmid:18760350.

14. Sadasivan S, Sharp B, Schultz-Cherry S, Smeyne RJ. Synergistic effects of influenza and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) can be eliminated by the use of influenza therapeutics: experimental evidence for the multi-hit hypothesis. NPJ Parkinsons Dis. 2017;3:18. pmid:28649618.

15. Henry J, Smeyne RJ, Jang H, Miller B, Okun MS. Parkinsonism and neurological manifestations of influenza throughout the 20th and 21st centuries. Parkinsonism Relat Disord. 2010;16(9):566–71. pmid:20650672

16. Sadasivan S, Pond BB, Pani AK, Qu C, Jiao Y, Smeyne RJ. Methylphenidate exposure induces dopamine neuron loss and activation of microglia in the basal ganglia of mice. PloS one. 2012;7(3):e33693. pmid:22470460.

17. Mattson MP. Infectious agents and age-related neurodegenerative disorders. Ageing Res Rev. 2004;3(1):105–20. pmid:15163105

18. Moreno JA, Streifel KM, Sullivan KA, Legare ME, Tjalkens RB. Developmental exposure to manganese increases adult susceptibility to inflammatory activation of glia and neuronal protein nitration. Toxicol Sci. 2009;112(2):405–15. pmid:19812365.

19. Zhang S, Zhou Z, Fu J. Effect of manganese chloride exposure on liver and brain mitochondria function in rats. Environ Res. 2003;93(2):149–57. Epub 2003/09/10. pmid:12963399.

20. Barcia C. Glial-mediated inflammation underlying parkinsonism. Scientifica (Cairo). 2013;2013:357805. Epub 2013/11/28. pmid:24278772.

21. Kim J, Pajarillo E, Rizor A, Son DS, Lee J, Aschner M, et al. LRRK2 kinase plays a critical role in manganese-induced inflammation and apoptosis in microglia. PloS one. 2019;14(1):e0210248. pmid:30645642.

22. Kirkley KS, Popichak KA, Afzali MF, Legare ME, Tjalkens RB. Microglia amplify inflammatory activation of astrocytes in manganese neurotoxicity. 2017:1–18. pmid:28476157

23. Kirkley KS, Popichak KA, Hammond SL, Davies C, Hunt L, Tjalkens RB. Genetic suppression of IKK2/NF-kappaB in astrocytes inhibits neuroinflammation and reduces neuronal loss in the MPTP-Probenecid model of Parkinson’s disease. Neurobiol Dis. 2019;127:193–209. Epub 2019/03/01. pmid:30818064.

24. Moreno JA, Streifel KM, Sullivan KA, Hanneman WH, Tjalkens RB. Manganese-induced NF-kappaB activation and nitrosative stress is decreased by estrogen in juvenile mice. Toxicol Sci. 2011;122(1):121–33. Epub 2011/04/23. pmid:21512103.

25. Moreno JA, Sullivan KA, Carbone DL, Hanneman WH, Tjalkens RB. Manganese potentiates nuclear factor-kappaB-dependent expression of nitric oxide synthase 2 in astrocytes by activating soluble guanylate cyclase and extracellular responsive kinase signaling pathways. J Neurosci Res. 2008;86(9):2028–38. Epub 2008/03/13. pmid:18335517.

26. al He. Juvenile exposure to manganese exacerbates astrocyte-dependent neuroinflammatory injury following adult challenge with MPTP. 2019.

27. Sadasivan S, Zanin M, O’Brien K, Schultz-Cherry S, Smeyne RJ. Induction of microglia activation after infection with the non-neurotropic A/CA/04/2009 H1N1 influenza virus. PloS one. 2015;10(4):e0124047. pmid:25861024.

28. Ramakrishnan MA. Determination of 50% endpoint titer using a simple formula. World J Virol. 2016;5(2):85–6. pmid:27175354.

29. Baquet ZC, Williams D, Brody J, Smeyne RJ. A comparison of model-based (2D) and design-based (3D) stereological methods for estimating cell number in the substantia nigra pars compacta (SNpc) of the C57BL/6J mouse. Neuroscience. 2009;161(4):1082–90. pmid:19376196.

30. Graeber MB, Streit WJ. Microglia: biology and pathology. Acta Neuropathol. 2010;119(1):89–105. pmid:20012873.

31. West MJ, Slomianka L, Gundersen HJ. Unbiased stereological estimation of the total number of neurons in thesubdivisions of the rat hippocampus using the optical fractionator. Anat Rec. 1991;231(4):482–97. pmid:1793176.

32. Cummings JE, Slayden RA. Transient In Vivo Resistance Mechanisms of Burkholderia pseudomallei to Ceftazidime and Molecular Markers for Monitoring Treatment Response. PLoS Negl Trop Dis. 2017;11(1):e0005209. pmid:28081127.

33. Popichak KA, Afzali MF, Kirkley KS, Tjalkens RB. Glial-neuronal signaling mechanisms underlying the neuroinflammatory effects of manganese. J Neuroinflammation. 2018;15(1):324. Epub 2018/11/23. pmid:30463564.

34. Tanuma N, Sakuma H, Sasaki A, Matsumoto Y. Chemokine expression by astrocytes plays a role in microglia/macrophage activation and subsequent neurodegeneration in secondary progressive multiple sclerosis. Acta Neuropathol. 2006;112(2):195–204. Epub 2006/05/31. pmid:16733654.

35. Harrison IF, Smith AD, Dexter DT. Pathological histone acetylation in Parkinson’s disease: Neuroprotection and inhibition of microglial activation through SIRT 2 inhibition. Neurosci Lett. 2018;666:48–57. Epub 2017/12/24. pmid:29273397.

36. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–D13. Epub 2018/11/27. pmid:30476243.

37. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–D8. Epub 2016/12/08. pmid:27924014.

38. Miner JJ, Sene A, Richner JM, Smith AM, Santeford A, Ban N, et al. Zika Virus Infection in Mice Causes Panuveitis with Shedding of Virus in Tears. Cell Rep. 2016;16(12):3208–18. pmid:27612415.

39. Hahn CS, Lustig S, Strauss EG, Strauss JH. Western equine encephalitis virus is a recombinant virus. Proc Natl Acad Sci U S A. 1988;85(16):5997–6001. pmid:3413072.

40. Seth P, Husain MM, Gupta P, Schoneboom A, Grieder BF, Mani H, et al. Early onset of virus infection and up-regulation of cytokines in mice treated with cadmium and manganese. Biometals: an international journal on the role of metal ions in biology, biochemistry, and medicine. 2003;16(2):359–68. pmid:12572694.

41. Wang C, Guan Y, Lv M, Zhang R, Guo Z, Wei X, et al. Manganese Increases the Sensitivity of the cGAS-STING Pathway for Double-Stranded DNA and Is Required for the Host Defense against DNA Viruses. Immunity. 2018;48(4):675–87.e7. Epub 2018/04/15. pmid:29653696.

42. Silva da Costa L, Pereira da Silva AP, Da Poian AT, El-Bacha T. Mitochondrial bioenergetic alterations in mouse neuroblastoma cells infected with Sindbis virus: implications to viral replication and neuronal death. PloS one. 2012;7(4):2. pmid:22485151

43. Ronca SE, Smith J, Koma T, Miller MM, Yun N, Dineley KT, et al. Mouse Model of Neurological Complications Resulting from Encephalitic Alphavirus Infection. Front Microbiol. 2017;8:188. pmid:28223982.

44. Elizan TS, Casals J. Astrogliosis in von Economo’s and postencephalitic Parkinson’s diseases supports probable viral etiology. Journal of the neurological sciences. 1991;105(2):131–4. Epub 1991/10/01. pmid:1757788.

45. Hirsch EC, Hunot S. Neuroinflammation in Parkinson’s disease: a target for neuroprotection? Lancet Neurol. 2009;8(4):382–97. pmid:19296921

46. Clarke LE, Liddelow SA, Chakraborty C, Munch AE, Heiman M, Barres BA. Normal aging induces A1-like astrocyte reactivity. Proc Natl Acad Sci U S A. 2018;115(8):E1896–E905. pmid:29437957

47. Liddelow SA, Guttenplan KA, Clarke LE, Bennett FC, Bohlen CJ, Schirmer L, et al. Neurotoxic reactive astrocytes are induced by activated microglia. Nature Publishing Group. 2017;541(7638):481–7. pmid:28099414

48. Popichak KA, Hammond SL, Moreno JA, Afzali MF, Backos DS, Slayden RD, et al. Compensatory expression of Nur77 and Nurr1 regulates NF-kappaB-dependent inflammatory signaling in astrocytes. Mol Pharmacol. 2018. 94(4):1174–1186 pmid:30111648

49. Streifel KM, Moreno JA, Hanneman WH, Legare ME, Tjalkens RB. Gene deletion of nos2 protects against manganese-induced neurological dysfunction in juvenile mice. Toxicol Sci. 2012;126(1):183–92. pmid:22174044

50. McCarthy MM, Nugent BM, Lenz KM. Neuroimmunology and neuroepigenetics in the establishment of sex differences in the brain. Nat Rev Neurosci. 2017;18(8):471–84. pmid:28638119

51. Song C, Kanthasamy A, Anantharam V, Sun F, Kanthasamy AG. Environmental neurotoxic pesticide increases histone acetylation to promote apoptosis in dopaminergic neuronal cells: relevance to epigenetic mechanisms of neurodegeneration. Mol Pharmacol. 2010;77(4):621–32. pmid:20097775

52. Song C, Kanthasamy A, Jin H, Anantharam V, Kanthasamy AG. Paraquat induces epigenetic changes by promoting histone acetylation in cell culture models of dopaminergic degeneration. Neurotoxicology. 2011;32(5):586–95. pmid:21777615

53. Guo Z, Zhang Z, Wang Q, Zhang J, Wang L, Zhang Q, et al. Manganese chloride induces histone acetylation changes in neuronal cells: Its role in manganese-induced damage. Neurotoxicology. 2018;65:255–63. pmid:29155171

54. Kostuk EW, Cai J, Iacovitti L. Subregional differences in astrocytes underlie selective neurodegeneration or protection in Parkinson’s disease models in culture. Glia. 2019. Epub 2019/04/27. pmid:31025779

55. Antipova D, Bandopadhyay R. Expression of DJ-1 in Neurodegenerative Disorders. Adv Exp Med Biol. 2017;1037:25–43. pmid:29147901

56. Wu YW, Prakash KM, Rong TY, Li HH, Xiao Q, Tan LC, et al. Lingo2 variants associated with essential tremor and Parkinson’s disease. Hum Genet. 2011;129(6):611–5. pmid:21287203

57. Vilarino-Guell C, Wider C, Ross OA, Jasinska-Myga B, Kachergus J, Cobb SA, et al. LINGO1 and LINGO2 variants are associated with essential tremor and Parkinson disease. Neurogenetics. 2010;11(4):401–8. pmid:20369371

58. Ranger A, Ray S, Szak S, Dearth A, Allaire N, Murray R, et al. Anti-LINGO-1 has no detectable immunomodulatory effects in preclinical and phase 1 studies. Neurol Neuroimmunol Neuroinflamm. 2018;5(1):e417. pmid:29259995

59. Takahashi H, Wakabayashi K. The cellular pathology of Parkinson’s disease. Neuropathology. 2001;21(4):315–22. pmid:11837539

60. Shinbo Y, Niki T, Taira T, Ooe H, Takahashi-Niki K, Maita C, et al. Proper SUMO-1 conjugation is essential to DJ-1 to exert its full activities. Cell Death Differ. 2006;13(1):96–108. pmid:15976810

61. Fan J, Ren H, Fei E, Jia N, Ying Z, Jiang P, et al. Sumoylation is critical for DJ-1 to repress p53 transcriptional activity. FEBS Lett. 2008;582(7):1151–6. pmid:18339323

62. Fan J, Ren H, Jia N, Fei E, Zhou T, Jiang P, et al. DJ-1 decreases Bax expression through repressing p53 transcriptional activity. J Biol Chem. 2008;283(7):4022–30. pmid:18042550

63. Delay C, Tremblay C, Brochu E, Paris-Robidas S, Emond V, Rajput AH, et al. Increased LINGO1 in the cerebellum of essential tremor patients. Mov Disord. 2014;29(13):1637–47. pmid:24531928

64. Civiero L, Cogo S, Kiekens A, Morganti C, Tessari I, Lobbestael E, et al. PAK6 Phosphorylates 14-3-3gamma to Regulate Steady State Phosphorylation of LRRK2. Front Mol Neurosci. 2017;10:417. pmid:29311810

65. Pandey A, Dan I, Kristiansen TZ, Watanabe NM, Voldby J, Kajikawa E, et al. Cloning and characterization of PAK5, a novel member of mammalian p21-activated kinase-II subfamily that is predominantly expressed in brain. Oncogene. 2002;21(24):3939–48. pmid:12032833

66. Schrantz N, da Silva Correia J, Fowler B, Ge Q, Sun Z, Bokoch GM. Mechanism of p21-activated kinase 6-mediated inhibition of androgen receptor signaling. J Biol Chem. 2004;279(3):1922–31. pmid:14573606

67. Robb GB, Brown KM, Khurana J, Rana TM. Specific and potent RNAi in the nucleus of human cells. Nat Struct Mol Biol. 2005;12(2):133–7. pmid:15643423

68. Weinmann L, Hock J, Ivacevic T, Ohrt T, Mutze J, Schwille P, et al. Importin 8 is a gene silencing factor that targets argonaute proteins to distinct mRNAs. Cell. 2009;136(3):496–507. pmid:19167051

69. Doyle M, Badertscher L, Jaskiewicz L, Guttinger S, Jurado S, Hugenschmidt T, et al. The double-stranded RNA binding domain of human Dicer functions as a nuclear localization signal. Rna. 2013;19(9):1238–52. pmid:23882114

70. Nishi K, Nishi A, Nagasawa T, Ui-Tei K. Human TNRC6A is an Argonaute-navigator protein for microRNA-mediated gene silencing in the nucleus. Rna. 2013;19(1):17–35. pmid:23150874

71. Gagnon KT, Li L, Chu Y, Janowski BA, Corey DR. RNAi factors are present and active in human cell nuclei. Cell Rep. 2014;6(1):211–21. pmid:24388755

72. Schraivogel D, Meister G. Import routes and nuclear functions of Argonaute and other small RNA-silencing proteins. Trends Biochem Sci. 2014;39(9):420–31. pmid:25131816

73. De A. Wnt/Ca2+ signaling pathway: a brief overview. Acta Biochim Biophys Sin (Shanghai). 2011;43(10):745–56. pmid:21903638

74. Kohn AD, Moon RT. Wnt and calcium signaling: beta-catenin-independent pathways. Cell Calcium. 2005;38(3–4):439–46. pmid:16099039

75. Sugimura R, Li L. Noncanonical Wnt signaling in vertebrate development, stem cells, and diseases. Birth Defects Res C Embryo Today. 2010;90(4):243–56. pmid:21181886

76. Stambolic V, MacPherson D, Sas D, Lin Y, Snow B, Jang Y, et al. Regulation of PTEN transcription by p53. Mol Cell. 2001;8(2):317–25. pmid:11545734

77. Ye Y, Rape M. Building ubiquitin chains: E2 enzymes at work. Nat Rev Mol Cell Biol. 2009;10(11):755–64. pmid:19851334

78. Reyes-Turcu FE, Ventii KH, Wilkinson KD. Regulation and cellular roles of ubiquitin-specific deubiquitinating enzymes. Annu Rev Biochem. 2009;78:363–97. pmid:19489724

79. Kloetzel PM, Ossendorp F. Proteasome and peptidase function in MHC-class-I-mediated antigen presentation. Curr Opin Immunol. 2004;16(1):76–81. pmid:14734113

80. Princiotta MF, Finzi D, Qian SB, Gibbs J, Schuchmann S, Buttgereit F, et al. Quantitating protein synthesis, degradation, and endogenous antigen processing. Immunity. 2003;18(3):343–54. pmid:12648452

81. Kormish JD, Sinner D, Zorn AM. Interactions between SOX factors and Wnt/beta-catenin signaling in development and disease. Dev Dyn. 2010;239(1):56–68. pmid:19655378

82. Topol L, Chen W, Song H, Day TF, Yang Y. Sox9 inhibits Wnt signaling by promoting beta-catenin phosphorylation in the nucleus. J Biol Chem. 2009;284(5):3323–33. pmid:19047045

83. Chen X, Yang J, Evans PM, Liu C. Wnt signaling: the good and the bad. Acta Biochim Biophys Sin (Shanghai). 2008;40(7):577–94. pmid:18604449

84. Sinner D, Kordich JJ, Spence JR, Opoka R, Rankin S, Lin SC, et al. Sox17 and Sox4 differentially regulate beta-catenin/T-cell factor activity and proliferation of colon carcinoma cells. Mol Cell Biol. 2007;27(22):7802–15. pmid:17875931

85. Bramham CR, Worley PF, Moore MJ, Guzowski JF. The immediate early gene arc/arg3.1: regulation, mechanisms, and function. J Neurosci. 2008;28(46):11760–7. pmid:19005037

86. Chowdhury S, Shepherd JD, Okuno H, Lyford G, Petralia RS, Plath N, et al. Arc/Arg3.1 interacts with the endocytic machinery to regulate AMPA receptor trafficking. Neuron. 2006;52(3):445–59. pmid:17088211

87. Wu J, Petralia RS, Kurushima H, Patel H, Jung MY, Volk L, et al. Arc/Arg3.1 regulates an endosomal pathway essential for activity-dependent beta-amyloid generation. Cell. 2011;147(3):615–28. pmid:22036569

88. Luo X, Huang S, Luo S, Liao H, Wang Y, Deng X, et al. Identification of genes underlying the enhancement of immunity by a formula of lentinan, pachymaran and tremelia polysaccharides in immunosuppressive mice. Sci Rep. 2018;8(1):10082. pmid:29973708

89. Lee HY, Itahana Y, Schuechner S, Fukuda M, Je HS, Ogris E, et al. Ca(2+)-dependent demethylation of phosphatase PP2Ac promotes glucose deprivation-induced cell death independently of inhibiting glycolysis. Sci Signal. 2018;11(512).

90. Santiago ADS, Counago RM, Ramos PZ, Godoi PHC, Massirer KB, Gileadi O, et al. Structural Analysis of Inhibitor Binding to CAMKK1 Identifies Features Necessary for Design of Specific Inhibitors. Sci Rep. 2018;8(1):14800. pmid:30287839

91. Wallace DC. The epigenome and the mitochondrion: bioenergetics and the environment [corrected]. Genes Dev. 2010;24(15):1571–3. pmid:20679390

92. Kukic I, Lee JK, Coblentz J, Kelleher SL, Kiselyov K. Zinc-dependent lysosomal enlargement in TRPML1-deficient cells involves MTF-1 transcription factor and ZnT4 (Slc30a4) transporter. Biochem J. 2013;451(2):155–63. pmid:23368743

93. Subramanian Vignesh K, Landero Figueroa JA, Porollo A, Caruso JA, Deepe GS Jr. Granulocyte macrophage-colony stimulating factor induced Zn sequestration enhances macrophage superoxide and limits intracellular pathogen survival. Immunity. 2013;39(4):697–710. Epub 2013/10/22. pmid:24138881

94. Subramanian Vignesh K, Landero Figueroa JA, Porollo A, Divanovic S, Caruso JA, Deepe GS Jr. IL-4 Induces Metallothionein 3- and SLC30A4-Dependent Increase in Intracellular Zn(2+) that Promotes Pathogen Persistence in Macrophages. Cell Rep. 2016;16(12):3232–46. pmid:27653687

95. Ravikumar B, Sarkar S, Davies JE, Futter M, Garcia-Arencibia M, Green-Thompson ZW, et al. Regulation of mammalian autophagy in physiology and pathophysiology. Physiol Rev. 2010;90(4):1383–435. pmid:20959619

96. Kochl R, Hu XW, Chan EY, Tooze SA. Microtubules facilitate autophagosome formation and fusion of autophagosomes with endosomes. Traffic. 2006;7(2):129–45. pmid:16420522

97. Kimura S, Noda T, Yoshimori T. Dynein-dependent movement of autophagosomes mediates efficient encounters with lysosomes. Cell Struct Funct. 2008;33(1):109–22. pmid:18388399

98. Mizushima N, Yoshimori T, Ohsumi Y. The role of Atg proteins in autophagosome formation. Annu Rev Cell Dev Biol. 2011;27:107–32. pmid:21801009

99. Musiek ES, Lim MM, Yang G, Bauer AQ, Qi L, Lee Y, et al. Circadian clock proteins regulate neuronal redox homeostasis and neurodegeneration. J Clin Invest. 2013;123(12):5389–400. pmid:24270424

100. Naruse Y, Oh-hashi K, Iijima N, Naruse M, Yoshioka H, Tanaka M. Circadian and light-induced transcription of clock gene Per1 depends on histone acetylation and deacetylation. Mol Cell Biol. 2004;24(14):6278–87. Epub 2004/07/01. pmid:15226430

101. Rosensweig C, Reynolds KA, Gao P, Laothamatas I, Shan Y, Ranganathan R, et al. An evolutionary hotspot defines functional differences between CRYPTOCHROMES. Nat Commun. 2018;9(1):1138. pmid:29556064

102. Yoshida K, Hashiramoto A, Okano T, Yamane T, Shibanuma N, Shiozawa S. TNF-alpha modulates expression of the circadian clock gene Per2 in rheumatoid synovial cells. Scand J Rheumatol. 2013;42(4):276–80. pmid:23496259

103. Zhang EE, Liu Y, Dentin R, Pongsawakul PY, Liu AC, Hirota T, et al. Cryptochrome mediates circadian regulation of cAMP signaling and hepatic gluconeogenesis. Nat Med. 2010;16(10):1152–6. pmid:20852621

104. Bonin RP, Zurek AA, Yu J, Bayliss DA, Orser BA. Hyperpolarization-activated current (In) is reduced in hippocampal neurons from Gabra5-/- mice. PLoS One. 2013;8(3):e58679. Epub 2013/03/22. pmid:23516534

105. Butler KM, Moody OA, Schuler E, Coryell J, Alexander JJ, Jenkins A, et al. De novo variants in GABRA2 and GABRA5 alter receptor function and contribute to early-onset epilepsy. Brain. 2018;141(8):2392–405. pmid:29961870

106. Hernandez CC, XiangWei W, Hu N, Shen D, Shen W, Lagrange AH, et al. Altered inhibitory synapses in de novo GABRA5 and GABRA1 mutations associated with early onset epileptic encephalopathies. Brain. 2019. pmid:31056671

107. Auburger G, Gispert S, Brehm N. Methyl-Arginine Profile of Brain from Aged PINK1-KO+A53T-SNCA Mice Suggests Altered Mitochondrial Biogenesis. Parkinsons Dis. 2016;2016:4686185. Epub 2016/04/02. pmid:27034888

108. Hecht I, Toporik A, Podojil JR, Vaknin I, Cojocaru G, Oren A, et al. ILDR2 Is a Novel B7-like Protein That Negatively Regulates T Cell Responses. J Immunol. 2018;200(6):2025–37. Epub 2018/02/13. pmid:29431694.

109. Podojil JR, Hecht I, Chiang MY, Vaknin I, Barbiro I, Novik A, et al. ILDR2-Fc Is a Novel Regulator of Immune Homeostasis and Inducer of Antigen-Specific Immune Tolerance. J Immunol. 2018;200(6):2013–24. pmid:29431690

110. Watanabe K, Watson E, Cremona ML, Millings EJ, Lefkowitch JH, Fischer SG, et al. ILDR2: an endoplasmic reticulum resident molecule mediating hepatic lipid homeostasis. PLoS One. 2013;8(6):e67234. pmid:23826244

111. Forstner P, Rehman R, Anastasiadou S, Haffner-Luntzer M, Sinske D, Ignatius A, et al. Neuroinflammation after Traumatic Brain Injury Is Enhanced in Activating Transcription Factor 3 Mutant Mice. J Neurotrauma. 2018;35(19):2317–29. pmid:29463176

112. Nicke B, Bastien J, Khanna SJ, Warne PH, Cowling V, Cook SJ, et al. Involvement of MINK, a Ste20 family kinase, in Ras oncogene-induced growth arrest in human ovarian surface epithelial cells. Mol Cell. 2005;20(5):673–85. pmid:16337592

113. Murphy LO, Smith S, Chen RH, Fingar DC, Blenis J. Molecular interpretation of ERK signal duration by immediate early gene products. Nat Cell Biol. 2002;4(8):556–64. pmid:12134156

114. Bracken AP, Kleine-Kohlbrecher D, Dietrich N, Pasini D, Gargiulo G, Beekman C, et al. The Polycomb group proteins bind throughout the INK4A-ARF locus and are disassociated in senescent cells. Genes Dev. 2007;21(5):525–30. pmid:17344414

115. Kotake Y, Cao R, Viatour P, Sage J, Zhang Y, Xiong Y. pRB family proteins are required for H3K27 trimethylation and Polycomb repression complexes binding to and silencing p16INK4alpha tumor suppressor gene. Genes Dev. 2007;21(1):49–54. pmid:17210787

116. Kotake Y, Sekiya Y, Okuda K, Ohta S. Cytotoxicity of 17 tetrahydroisoquinoline derivatives in SH-SY5Y human neuroblastoma cells is related to mitochondrial NADH-ubiquinone oxidoreductase inhibition. Neurotoxicology. 2007;28(1):27–32. pmid:16860391

117. Voncken JW, Niessen H, Neufeld B, Rennefahrt U, Dahlmans V, Kubben N, et al. MAPKAP kinase 3pK phosphorylates and regulates chromatin association of the polycomb group protein Bmi1. J Biol Chem. 2005;280(7):5178–87. pmid:15563468

118. Agger K, Cloos PA, Rudkjaer L, Williams K, Andersen G, Christensen J, et al. The H3K27me3 demethylase JMJD3 contributes to the activation of the INK4A-ARF locus in response to oncogene- and stress-induced senescence. Genes Dev. 2009;23(10):1171–6. pmid:19451217

119. Barradas M, Anderton E, Acosta JC, Li S, Banito A, Rodriguez-Niedenfuhr M, et al. Histone demethylase JMJD3 contributes to epigenetic control of INK4a/ARF by oncogenic RAS. Genes Dev. 2009;23(10):1177–82. pmid:19451218

120. Lin SY, Li TY, Liu Q, Zhang C, Li X, Chen Y, et al. GSK3-TIP60-ULK1 signaling pathway links growth factor deprivation to autophagy. Science. 2012;336(6080):477–81. pmid:22539723

121. Sun YW, Herzog CR, Krzeminski J, Amin S, Perdew G, El-Bayoumy K. Effects of the environmental mammary carcinogen 6-nitrochrysene on p53 and p21(Cip1) protein expression and cell cycle regulation in MCF-7 and MCF-10A cells. Chem Biol Interact. 2007;170(1):31–9. pmid:17678638