Abstract Parkinson disease (PD) is characterized by the preferential, but poorly understood, vulnerability to degeneration of midbrain dopaminergic (mDA) neurons in the ventral substantia nigra compacta (vSNc). These sensitive mDA neurons express Pitx3, a transcription factor that is critical for their survival during development. We used this dependence to identify, by flow cytometry and expression profiling, the negative regulator of G-protein signaling Rgs6 for its restricted expression in these neurons. In contrast to Pitx3−/− mDA neurons that die during fetal (vSNc) or post-natal (VTA) period, the vSNc mDA neurons of Rgs6−/− mutant mice begin to exhibit unilateral signs of degeneration at around 6 months of age, and by one year cell loss is observed in a fraction of mice. Unilateral cell loss is accompanied by contralateral degenerating neurons that exhibit smaller cell size, altered morphology and reduced dendritic network. The degenerating neurons have low levels of tyrosine hydroxylase (TH) and decreased nuclear Pitx3; accordingly, expression of many Pitx3 target gene products is altered, including Vmat2, Bdnf, Aldh1a1 (Adh2) and Fgf10. These low TH neurons also express markers of increased dopamine signaling, namely increased DAT and phospho-Erk1/2 expression. The late onset degeneration may reflect the protective action of Rgs6 against excessive DA signaling throughout life. Rgs6-dependent protection is thus critical for adult survival and maintenance of the vSNc mDA neurons that are most affected in PD.

Author Summary The locomotor deficits associated with Parkinson disease result from the death of a specific subset of dopamine neurons in the ventral part of the midbrain. The reason for the greater sensitivity to degeneration of those, relative to other, neurons is not clear. Prior work showed that the Pitx3 transcription factor is specifically expressed in these neurons where it has a survival role during development. The present work identified a cell signaling component, Rgs6, that is also restricted to the sensitive neurons in the midbrain and that exerts a protective function, particularly late in life. While the loss of Rgs6 function may predispose or contribute to Parkinson disease, its stimulation may provide a novel therapeutic avenue to treat Parkinson disease.

Citation: Bifsha P, Yang J, Fisher RA, Drouin J (2014) Rgs6 is Required for Adult Maintenance of Dopaminergic Neurons in the Ventral Substantia Nigra. PLoS Genet 10(12): e1004863. https://doi.org/10.1371/journal.pgen.1004863 Editor: Gregory P. Copenhaver, The University of North Carolina at Chapel Hill, United States of America Received: July 31, 2014; Accepted: October 29, 2014; Published: December 11, 2014 Copyright: © 2014 Bifsha et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files. Funding: Funding for this project was provided by NIH CA161882 to RAF and by the Parkinson Society of Canada and Canadian Institutes of Health Research MOP-123213 to JD and FRSQ scholarship to PB. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction Parkinson disease (PD) is characterised by the progressive loss of midbrain dopaminergic (mDA) neurons [1]. Although the clinical manifestations of PD can be variable, the appearance of motor deficits is the hallmark of this neurodegenerative disease. Similarly, the etiology of PD appears to be multifactorial but one consistent feature of this disease is the greater sensitivity of ventral substantia nigra compacta (vSNc) mDA neurons to degenerate [2] as opposed to mDA neurons of the dorsal SNc (dSNc) and ventral tegmental area (VTA). The molecular basis for this preferential sensitivity remains poorly understood although work in animal models has been useful [3]. Vertebrate animal models based on genetic causes of PD, which constitute 10–20% of PD cases, have not been overly successful in reproducing the selective neurodegeneration patterns of the mDA system [4]. For example, mice with transgenic expression of human autosomal dominant mutants of α-synuclein (SNCA) or leucine-rich repeat kinase 2 (LRRK2) rarely produce mDA neurodegeneration [5], [6]. Murine loss-of-function mutations in autosomal recessive gene products for PTEN induced putative kinase 1 (Pink) and Parkinson protein 2 (Park2), have not been enlightening either, with the recent exception of Parkinson protein 7 (Park7 or DJ-1). Indeed, DJ-1−/− mice show progressive adult degeneration of SNc mDA neurons upon backcrossing to an appropriate genetic background [7], indicating that many factors are necessary in order to model polygenic diseases such as PD. However, it is noteworthy that rat models of Pink1 and DJ-1 loss-of-function showed progressive loss of mDA neurons [8]. Both early-onset and late-onset forms of PD bear a major histopathological hallmark, the presence of Lewy bodies, which are α-synuclein-rich protein inclusions that are also found in non-dopaminergic brain regions depending on the stage of disease progression. Many familial PD genes have widespread brain expression, without any preferential expression in mDA subpopulations [5]. In general, they all participate in similar inter-related cellular processes such as in mitochondrial function (PINK1, DJ-1, SNCA), the secretory pathway (PARK2, LRRK2) and the ubiquitin-proteasome degradation pathway (PINK1, SNCA, PARK2). Rodent animal models based on genes that participate in development and survival of mDA neurons (Pitx3, Nurr1, Girk2, En1/2, Otx2, etc) have proved extremely useful, especially in defining the candidate cellular pathways underlying the respective differential vulnerability of SNc versus VTA mDA neuron subpopulations to toxin-induced neurodegeneration and in human pathology [9]–[11]. Notably, mouse mutants for the homeobox transcription factor Pitx3 (Entrez gene ID: 5309) are unique in that they display a selective and stereotypic pattern of mDA cell loss that resembles typical PD [12]–[14]. In particular, Pitx3-deficient mice exhibit developmental loss of Pitx3-positive Calbindin (Calb)-negative mDA neurons of the vSNc (Pitx3-dependent) while Pitx3-negative Calb-positive mDA neurons of dSNc and VTA (Pitx3-independent) remain essentially unaffected by Pitx3 deficiency [15]. This cell loss is associated with a loss of spontaneous movement that can be partially rescued by L-dopa treatment [16], [17]. Human PITX3 polymorphisms are associated with sporadic PD [18]. Pitx3 mediates its effects by regulating the expression of many genes (Aldh1a1, DAT, Drd2, TH, Bdnf) in a subset specific fashion [19]. A well-studied Pitx3 target gene in mDA neurons of the SNc is Aldh1a1, which is important for retinoic acid production and subsequent neuronal maturation and protection through regulation of TH expression [19]–[21]. Other Pitx3 target genes, such as the classical dopaminergic markers DAT, Vmat2, Drd2 are important for neurotransmitter identity and are subject to the cooperative action of Pitx3 with Nurr1 [22]. On the other hand, Pitx3 expression has been reported to be itself regulated by GDNF, especially during development [23]. GDNF is the only neurotrophic factor for which conditional inactivation in the adult mouse has provided strong evidence of its absolute requirement for the cell-autonomous survival of brain cathecholaminergic neurons, including mDA neurons [24]. Although, many studies described the implication of Pitx3 in post-natal maturation and developmental survival of mDA neurons, Pitx3 has not been conclusively linked to mechanisms of survival and maintenance of mDA neurons in the adult, especially as it pertains to degenerative processes. We hypothesized that Pitx3-controlled genes and pathways, in addition to their known role in development, may also be implicated in the neuroprotective pathway required to maintain the integrity of specific subset of mDA neurons throughout adulthood. In a screen of expression profiling data comparing Pitx3-dependent and -independent FACS-purified mDA neurons of SNc and VTA, the Regulator of G-protein Signaling 6 (Rgs6) (Entrez Gene ID: 9628) was identified as a putative survival factor that is preferentially expressed in vSNc mDA neurons and whose expression is positively regulated by Pitx3. Rgs6 belongs to the R7 subfamily of Rgs and functions as a GTPase activating protein to terminate signaling downstream of ligand-bound G-protein coupled receptors (GPCR). It does so by accelerating the conversion from the active Gα-GTP bound state (dissociated from Gβγ subunit) to the inactive Gα-GDP bound state (associated to Gβγ subunit). The R7 subfamily of Rgs regulators, including Rgs6, are known to have preference for catalysis of pertussis toxin-sensitive Gi/o heterotrimeric G-proteins through recognition of their Gαi by a C-terminal Rgs protein domain [25]. Activated Gαi subunits inhibit adenylate cyclase such that cAMP production from ATP is halted and PKA/cAMP-dependent protein kinase pathways are inhibited. Conversely, the activated Gβγ subunit opens Girk channels to allow efflux of potassium ions outside the cell resulting in hyperpolarization [26]. The neuronal GPCRs previously associated with Gi/o proteins include dopamine receptors (Drd2, Drd3), acetylcholine receptors (m2, m4), GABA B receptor, metabotropic glutamate receptors (mGluR2, 3, 4, 6, 7, 8). Thus in cerebellum and heart, phenotypes resulting from inactivation of Rgs6 are consistent with over-activation of signaling downstream of GABA B , serotonin 5-HT 1A , M2 acetylcholine receptors, respectively. [26]–[29]. In the present work, we identified Rgs6 and investigated its role in the midbrain dopaminergic system. Rgs6 is shown to be preferentially expressed in vSNc mDA neurons and its knockout in mice results in progressive loss and alterations of Pitx3-positive mDA neurons specifically within the vSNc of aged animals. This late-onset degeneration is associated with markers of increased Drd2 signaling, down-regulation of Pitx3 expression and deregulated expression of its target genes, Aldh1a1, Bdnf, Vmat2, TH and Fgf10. Further, the pattern of mDA degeneration observed in Rgs6−/− mice is a close phenocopy of DJ-1−/− mice suggesting that these two genes may act through similar pathways.

Discussion A few gene expression profiling studies have compared gene expression in SNc versus VTA [46]–[48]. These studies identified large numbers of SN or VTA restricted genes but did not include criteria to relate these specificities to function. Our reliance on Pitx3 gene dependence to identify genes with preferential expression in Pitx3-independent dSNC versus Pitx3-dependent vSNc neurons allowed the prioritization of candidate genes based on the pro-survival Pitx3 gene. This approach also allowed definition of the unique expression profiles for dorsal compared to ventral SNc mDA neurons. We focused on vSNc-enriched and Pitx3-dependent genes in order to identify candidates for role(s) in vSNc mDA neurons survival. The list of 10 Pitx3-activated and 6 Pitx3-repressed genes in this subset includes only one known gene encoding a regulator of a signaling pathway, Rgs6. Further, Rgs6 is the most dependent on Pitx3 for expression in VTA. The list also includes a transcriptional co-regulator Lmo3 that may contribute to Pitx3-dependent survival pathways but it's putative involvement in survival appeared less likely than Rgs6 to contribute to an age-dependent phenotype. Thus, the present study focused on characterization of neuronal loss and on molecular features of degenerating vSNc neurons as a result of Rgs6 inactivation. We observed late-onset degeneration of vSNc mDA neurons in Rgs6−/− midbrains (Fig. 9B). This phenotype is detected at about 6 months of age and becomes more important in 1y-old mice. At that age, clear mDA cell loss is observed in a subset of mice (Table S2). It is likely that the degenerating THlow mDA phenotype (Fig. 5) and the mDA cell loss (Fig. 4) represent progressive steps of the same defects; accordingly, THlow neurons are also observed in the vSNc of midbrains with cell loss. This late-onset degeneration is similar to another mouse model of monogenic PD, the DJ-1 (Parkin 7) mutant, that also initially exhibits unilateral defects [7]. Indeed, DJ-1−/− mice present unilateral loss of mDA cell bodies as early as 2 months after birth, with a transition to bilateral cell loss occurring at 1 year of age. Our study of Rgs6−/− mice showed unilateral cell loss, as evidence by decreased number of TH+ and NissL+ neurons (Fig. 4A–C, Table S2) at 12 months after birth, while evidence of THlow degenerating neurons is readily apparent earlier at 6 months of age (Fig. 5A, Table S2). The comparison of phenotypes for DJ-1 and Rgs6 knockout mice indicates that: 1) they both have selective degeneration of SNc, but not VTA, mDA neurons, 2) they both have progressive degeneration and loss of mDA neurons, 3) degeneration begins unilaterally, 4) degeneration eventually becomes bilateral with earlier transition in DJ-1−/− than Rgs6−/− mice. A distinguishing feature of the Rgs6−/− model is the bias towards degeneration of Calb-negative vSNc mDA neurons compared to Calb-positive dSNc mDA neurons that remain largely unaffected (Fig. 5A, Figure S4), as is usually observed in PD. The slow degeneration of THlow vSNc neurons provided an opportunity to define the molecular features that accompany the dysmorphology. These neurons display increased expression of markers previously associated with pathological changes such as FluroJade C, LC3B and phospho-p27Kip1. Moreover, degenerating THlow vSNc neurons show decreased DJ-1 and elevated Pink1 and Lrrk2 protein expression, suggesting a relationship between Rgs6 signaling and pathways implicated in PD pathology (Fig. 6). Future studies should address the relationships between DJ-1, Pink1 and Lrrk2 in degeneration pathways of Rgs6−/− mice in terms of their known roles in mitochondria dynamics, calcium, balance, redox state and cell signaling. Collectively, the data show that Rgs6 signaling is necessary for maintenance of vSNc mDA neurons in the aging animal and that its downstream action may be mediated, at least in part, by Pitx3-dependent mechanisms (Fig. 9A). Indeed, the THlow vSNc neurons exhibit low levels of Pitx3 and of its target gene products TH, Aldh1a1, Bdnf, Vmat2, together with enhanced expression of Pitx3-repressed Fgf10 (Fig. 7). The THlow neurons also exhibit cytoplasmic Pitx3 staining suggesting that there may be regulation of nuclear-cytoplasmic localization: this effect could be mediated through phosphorylation of Pitx3 as it was suggested that phosphorylated Pitx1 has greater affinity for nuclear DNA binding than its de-phosphorylated form [49]. One of the documented Pitx3-activated factors, Bdnf, is an important mediator of the neuroprotective action of Pitx3 during development and could contribute to trophic impairment in adult degenerating neurons. The degenerating neurons also exhibit a loss of determinants of postmitotic dopaminergic identity (TH, Aldh1a1, Vmat2) and this likely affects their neuronal activity. Aberrant neuronal activity is a hallmark [45] of the pre-symptomatic stage of PD and this likely transitions to major cellular disruptions (proteosome dysfunction, mitochondrial integrity, calcium permeability…) associated with degeneration and cell death. What could be the target of Rgs6 action? One likely possibility is the dopamine receptor D2 (Drd2). Indeed, Rgs6 is a negative modulator of GPCR activity, including the dopamine Drd2 receptor [25]. Drd2 expression itself was not affected in degenerating Rgs6−/− mDA neurons (Fig. 8A). However, the vSNc THlow neurons showed evidence of increased DA signaling, namely accumulation of phospho-Erk1/2 (Fig. 8B) and enhanced glycosylated dopamine transporter (Slc6a3/DAT) expression in the mutant (Fig. 8C), consistent with a putative loss of negative Rgs6 input on dopamine signaling [50]. Since the Pitx3−/− midbrain exhibits decreased DAT and Drd2 [19], the observed increase in DAT together with phospho-Erk1/2 are consistent with a primary action of Rgs6 inactivation on DA signaling. Rgs6 may thus contribute to the auto-regulatory negative feedback of the dopaminergic system and its absence may lead to dopamine-dependent oxidative stress and neuronal loss [2], [50]. The enhanced phospho-p27Kip1 and FluoroJade staining support the interpretation that these cells are under stress. Alternatively, Rgs6 may have GPCR-independent actions: those could involve the GDNF pathway that is essential for catecholaminergic neuron survival [24] or involve direct action on apoptotic pathways [39]. An important aspect of the expression changes discussed above is that they only occur in THlow dysmorphic neurons that normally express Rgs6 and not in other mDA neurons of VTA and dSNc that are negative for Rgs6 (Fig. 3). Therefore, the concordance between Rgs6 midbrain expression and observed cell degeneration patterns suggests that the changes are cell-autonomous and directly related to Rgs6-dependent signaling operating in Rgs6+Pitx3+ vSNc neurons. We cannot however rule out the contribution of other brain systems affected by Rgs6 deficiency. Collectively, the data show that Rgs6 signaling is necessary for maintenance of vSNc mDA neurons in the aging animal and that its downstream action may be mediated, at least in part, by Pitx3-dependent mechanisms (Fig. 9). The present work identified a critical signalling pathway that controls survival of the mDA neuron subset that preferentially degenerates in PD. Further dissection of this pathway may lead to therapeutically useful insights on the unique properties of this group of mDA neurons.

Materials and Methods Ethics statement All experimental procedures with laboratory animals were approved by the IRCM Animal Protection Committee and followed guidelines and regulations of the Canadian Council of Animal Care. Animal models All mice were maintained as heterozygous carriers in the C57Bl/6J background and maintained on a 12 h light-dark cycle with food and water ad libitum. Rgs6-null [26] and TH-EGFP [51] mice were described previously. Pitx3-null mice were generated in this laboratory [52]. Dissections and flow cytometry Ventral midbrain dissections were performed on WT and Pitx3−/− newborn mice (P1–P4) crossed onto TH-EGFP heterozygous background. Mouse brains were quickly washed in ice-cold PBS and then placed into cold Hibernate-A/1%B27 solution (Gibco) to dissect EGFP+ ventral midbrain (vMB) tissue under the fluorescence stereoscope (Leica DFC300 FX). Tissue blocks of vMB were then further microdissected so as to separate SNc (lateral) from VTA (medial). VTA and SNc tissue blocks were digested using the papain dissociation kit (Worthington). Dissociated cells were then resuspended in warm Hibernate-A/1%B27 solution containing propidium iodide (PI, 1 µg/ml), passed through 100 µm mesh and sorted by flow cytometry using the MOFLO™ instrument (Beckman Coulter). PI−/EGFP+ live sorted cells were deposited in 30 µl of RNAlater solution (Ambion) (max. of 3000 cells per 30 µl of RNAlater) to preserve RNA integrity. RNA extraction, quantitative real-time PCR (qRT-PCR) and microarray analyses Sorted EGFP+ cells in RNAlater were processed in batches of approximately 5000 cells for purification of total RNA using RNeasy Micro kit (Qiagen). Briefly, 350 µl of RLT lysis buffer was added per 30 µl RNAlater-suspended cells. After vortexing for 1 min, 1 volume of 70% ethanol was added and the content loaded into single pre-equilibrated RNeasy MinElute column. This was done in duplicate for each of the four different preparations of EGFP+ cells (Pitx3+/+ SNc, Pitx3+/+ VTA, Pitx3−/− SNc, Pitx3−/− VTA). Column-bound RNA was washed as recommended and eluted with 14 µl of RNAse-free water. Quality of total RNA was verified with the Agilent RNA 6000 Nano kit adapted for Agilent 2100 Bioanalyzer. For RT-qPCR, first-strand cDNA was synthesized using Superscript III RT enzyme and accompanying kit (Invitrogen). Primers for PCR amplification are displayed in Table S3. qPCR was performed using Perfecta reagents (Quanta) on a MX-3005 device (Stratagene), and results were analyzed using the accompanying software. All quantifications were relative to Gapdh mRNA. For microarray hybridization, total RNA with minimal degradation was used. Prior to hybridization onto Affymetrix Gene1.0ST expression arrays, which was done at the Genome Quebec/McGill Innovation Centre, total RNA was linearly amplified using WT-Ovation Pico RNA Amplification kit (NuGen). Gene expression summary values were computed by RMA Express (Bolstad et al, 2003) and raw data was normalized with the LPE algorithm embedded in the FlexArray suite of programs (Genome Quebec). Differentially expressed genes were chosen on the basis of their p-value≥0.05, fold-enrichment ≥1.5 and raw array signal ≥60. Hierarchical clustering and heat map display was done using Genesis (Institute for Genomics and Bioinformatics, Austria) Immunohistochemistry Mice were anesthetized and perfused intra-cardially with fresh 4% paraformaldehyde/PBS buffer. Brains were collected and post-fixed for 24 h at 4°C. After inclusion in paraffin, brains were cut in 5–6 µm coronal sections using microtome and mounted on Superfrost Plus (Fisher Scientific) slides. Immunonohistochemistry was performed after paraffin removal and hydration through xylene and graded alcohol series. Antigen retrieval was performed in 10 mM sodium citrate (pH 8.5) at 80°C water bath for 30 min. Sections in citrate solution were left to cool to room-temperature (RT) after which a step of endogenous biotin block was performed (Streptavidin/biotin kit, Vector Labs). Blocking with 5% normal serum for 1 h at RT preceded primary antibody incubation (overnight at 4°C). Primary antibodies used were against Pitx3 (rabbit home-made, 1∶400), Rgs6 (rabbit home-made 1∶25), Th (Millipore MAB318, 1∶1000), Th (Millipore AB152, 1∶500), Calb1 (R&D Systems AF3320, 1∶40), Fgf10 (Millipore ABN44, 1∶1000), Slc6a3 (Santa Cruz sc-32258, 1∶250), Drd2 (Santa Cruz sc5303, 1∶100), Aldh1a1 (Abcam ab24343, 1∶400), BDNF (Abcam ab108319, 1∶25), phospho-p27 (Abcam ab32096, 1∶50), phospho-Erk1/2 (Cell Signaling 4376, 1∶25), DJ-1 (Abcam ab18257, 1∶500) Lrrk2 (Epitomics 3514-1, 1∶50), Pink1 (Novus Biologicals BC100-494, 1∶50), LC3B (Cell Signaling 3868, 1∶50), and Vmat2 (Millipore AB1598P, 1∶200). Secondary antibodies were either biotinylated (Vector Labs, 1∶250) or directly coupled to fluorochromes such as AlexaFluor 488/546 (Invitrogen, 1∶250). For ABC method of amplification, AlexaFluor350/480/546- or HRP-conjugated (PerkinElmer NEL750, 1∶1000) streptavidin compounds were used. For immunoperoxidase staining, 1% hydrogen peroxide treatment was done just after antigen retrieval step. Mounting of sections was in Mowiol (Sigma) (fluorescence) or in Permount (Fisher Scientific) (DAB chromogen reaction). Immunofluorescence sections were observed on Leica DM6000B light microscope and Carl Zeiss LSM700 confocal microscope. Fluoro-Jade (FJC) and Nissl staining FJC (Molecular Probes, 0.01% stock in water) staining was performed just after the immunohistochemical procedure [33], [34]. Briefly, sections were washed in water, dipped for 5 min in 0.06% of potassium permanganate, rinsed again in water and placed in 0.0001% FJC/0.1% acetic/0.0001% DAPI solution for 10 min. Finally, slides were rinsed in water and mounted in 0.1% acetic acid/80% glycerol (v/v). The Fluorescein/FITC filter system was used to visualize FJC. Nissl staining was performed on deparaffinized sections by immersion into warm 0.1% cresyl violet solution for 10 min, rinsing three times in distilled water and differentiating in 95% ethanol. Slides were then dehydrated in 100% alcohol, cleared in xylene and mounted with Permount. Cell counts and statistics Cell counts were performed using ImageJ (National Institutes of Health). For cell counts of degenerating neurons, TH-stained or TH/FJC/Dapi-stained coronal sections were loaded on ImageJ; the sections spanned regular intervals (30 or 100 µm) across rostro-caudal extent of midbrain of WT (n = 2) and 1-yo Rgs6−/− mice (n = 2). For each section, total numbers of TH+, TH+/FJC+ and TH−/FJC+ cells were separately counted for SNc and VTA in both hemispheres. The percentage of TH+ degenerating neurons for each anatomical region reflects the ratio between the total number of TH+/FJC+ events and the total number of TH+ events for all rostro-caudal series. Values are reported as means +/−S.D. Statistical significance was calculated using Student T-test.

Acknowledgments We thank Abdul-Rasul Tabasum for secretarial assistance. We thank Dominique Lauzier for histology processing, Eric Massicotte for flow cytometry, Konstantin Khetchoumian for brain dissections and Ted Fon for antibodies. TH-EGFP mice were kindly provided by Drs L.E Trudeau (Montréal) and H. Okano (Tokyo).

Author Contributions Conceived and designed the experiments: PB JY RAF JD. Performed the experiments: PB JY. Analyzed the data: PB JY RAF JD. Contributed reagents/materials/analysis tools: PB JY RAF JD. Wrote the paper: PB JY RAF JD.