The atypical RhoGTPase RhoE/Rnd3 is a key molecule to acquire a neuroprotective phenotype in microglia
Journal of Neuroinflammation volume 15, Article number: 343 (2018)
Over-activated microglia play a central role during neuroinflammation, leading to neuronal cell death and neurodegeneration. Reversion of over-activated to neuroprotective microglia phenotype could regenerate a healthy CNS-supporting microglia environment. Our aim was to identify a dataset of intracellular molecules in primary microglia that play a role in the transition of microglia to a ramified, neuroprotective phenotype.
We exploited the anti-inflammatory and neuroprotective properties of conditioned medium of adipose-derived mesenchymal stem cells (CM) as a tool to generate the neuroprotective phenotype of microglia in vitro, and we set up a microscopy-based siRNA screen to identify its hits by cell morphology.
We initially assayed an array of 157 siRNAs against genes that codify proteins and factors of cytoskeleton and activation/inflammatory pathways in microglia. From them, 45 siRNAs significantly inhibited the CM-induced transition from a neurotoxic to a neuroprotective phenotype of microglia, and 50 siRNAs had the opposite effect. As a proof-of-concept, ten of these targets were validated with individual siRNAs and by downregulation of protein expression. This validation step resulted essential, because three of the potential targets were false positives. The seven validated targets were assayed in a functional screen that revealed that the atypical RhoGTPase RhoE/Rnd3 is necessary for BDNF expression and plays an essential role in controlling microglial migration.
Besides the identification of RhoE/Rnd3 as a novel inducer of a potential neuroprotective phenotype in microglia, we propose a list of potential targets to be further confirmed with selective activators or inhibitors.
Neuroinflammation is a fundamental process contributing to the death of neurons in neurodegerenative diseases, for example, in Parkinson’s Disease (PD) , Alzheimer’s Disease, and Multiple Sclerosis . The underlying neuroinflammation in these neurodegenerative diseases are characterized by microglia activation [2, 3]. During this inflammatory process, microglia acquire an amoeboid cell shape , and amongst other mechanisms, they decrease the release of neurotrophic factors, such as the brain-derived neurotrophic factor (BDNF), and secrete cytotoxic substances, which lead to neuronal death . In contrast, in a healthy brain, microglia are found in a ramified morphology and support a variety of central nervous system (CNS) functions by secreting neurotrophic factors. Interestingly, microglia are not permanent in one or another activity state, but are rather able to switch between different phenotypes and activity states . Therefore, understanding the molecular machinery that reverses the inflammatory activation of microglia is essential to protect from neurodegeneration. Hence, to identify the intracellular players of this reversion in a systematic manner, we combined two experimental approaches in this study.
Firstly, our previous work indicates that adipose-derived mesenchymal stem cells (ASCs) exert important anti-inflammatory actions on microglia. We observed that primary microglia exposed to ASCs or their secreted factors (conditioned medium, CM) underwent a dramatic cell shape change into a highly elongated morphology in vitro , similar to the phenotype of microglia observed in a healthy brain . The elongation induced by CM was associated with an upregulation of neurotrophic factors, such as BDNF, which is indicative of the acquisition of a neuroprotective phenotype . Thus, CM-stimulated microglia represent an ideal tool to study the intracellular events necessary for the transition from inflammatory activated to non-inflammatory neuroprotective microglia. Indeed, we previously identified that the small RhoGTPases Rac1 and Cdc42, which are important regulators of the actin cytoskeleton [9, 10], play major roles in this phenotypic transition .
Secondly, we present here our novel results of a small interference RNA (siRNA) screen targeting a panel of cytoskeletal and inflammatory genes in primary microglia. For our siRNA screen, we chose a microscopy-based screen as readout, since the phenotype transition of microglia is easily detectable by light microscopy and its outcome is directly visible. Although RNAi screens have been applied to create big datasets of image-based analysis of cell lines in high-throughput format [11, 12], in the case of microglia available, cell lines do not always resemble the physiology and morphology of primary microglia  and our own experience. Therefore, despite the technical challenge associated with the use of primary cells, we chose primary murine microglia as our model system. We were able to develop a reliable protocol to transfect primary microglia with a panel of siRNA, and optimized the experimental assay for its reproducible readout. Then, we exploited the anti-inflammatory properties of CM on microglia , as described above, in combination with the siRNA screen, and in this way, we identified a list of molecules that were implicated in the reversion from activated to neuroprotective microglia. From this list, we validated seven candidates, of which three of them also downregulated BDNF expression and other three molecules enhanced microglia migration. We found that the atypical RhoGTPase RhoE/Rnd3 was a common hit in both secondary screens enabling us to assign a novel function as a potential regulator of a neuroprotective microglia phenotype to this molecule.
Animals and ethical statement
Female C57Bl/6 mice (7–8 weeks old) were purchased from Charles River. They were used for obtaining P0 to P2 newborns to isolate primary microglia and for obtaining fat tissue to isolate ASCs. Mice were housed in a controlled-temperature/humidity environment (22 ± 1 °C, 60–70% relative humidity) in individual cages (10 mice per cage, with wood shaving bedding and nesting material), with a 12-h light/dark cycle (lights on at 7:00 a.m.) and fed with rodent chow (GlobalDiet 2018, Harlan) and tap water ad libitum. The experimental protocols of this study conform to EU Directive 2010/63 and followed the ethical guidelines for investigations with experimental animals approved by the Ethics Review Committee for Animal Experimentation of Spanish Council of Scientific Research. Animal studies are reported in compliance with the ARRIVE guidelines .
Cell isolation and cultures
Primary microglia were prepared from P0 to P2 newborn C57Bl/6 mice, as described previously . Briefly, dissected brains were discarded of olfactory bulb, cerebellum, hindbrain, and meninges. Then, they were homogenized in microglia growth medium consisting in DMEM (Invitrogen), supplemented with 10% fetal bovine serum (FBS, Gibco), 10% horse serum, and 1% penicillin/streptomycin (Gibco), using a Pasteur pipette and a 23-G syringe. Brain homogenates were centrifuged. The resulting cells were plated in poly-d-lysine-coated flasks and incubated with microglia growth medium for 10–12 days at 37 °C and 5% CO2. Microglia were harvested on a shaker for 2 h at 200 rpm and plated on poly-d-lysine-coated cover slips at a density of 12,000 cells/cm2 or tissue culture 6-well plates in microglia growth medium at a density of 37,000 cells/cm2. After 24 h, microglia growth medium was replaced by fresh microglia growth medium.
ASCs were isolated from adipose tissue of adult C57Bl/6 mice as previously described . These cells showed a fibroblast-like morphology and differentiation capacity to the adipocytic and osteocytic lineages and expressed the phenotype MHC−II−CD14−CD18−CD31−CD34−CD45−CD80−CD117−CD144−CD13+CD44+CD29+CD54+CD73+CD90+CD105+CD106+CD166+. CM was collected from passage 2 until passage 6 of ASC cultures, which were plated at a cell density of 15,000 cells/cm2, grown for 2 days before collecting their supernatant and then stored at − 20 °C. Before use, CM was quickly thawed and passed through a 0.2-μm filter.
A custom-made siRNA library against 157 target genes (with three individual siRNAs against each gene) was ordered from Life Technologies. Primary microglia were plated for 48 h as described above. Then, the growth medium was replaced and cells grown on cover slips in 48-well plates were transfected with 12.5 pmol siRNA and 0.375 μl of Lipofectamine 3000 (Invitrogen) according to the manufacturer’s instructions. For 24-, 12-, or 6-well plate cultures, reagents were proportionally up-scaled. Four days after the siRNA transfection, cells were treated with CM for 4 h and fixed to determine their cell shape. Alternatively, after 4 days, cells were harvested for Western blot or RT-qPCR analysis. Not transfected cells or cells that were transfected with a control Scrambled siRNA were used as reference in all studies.
Determination of cell shape and circularity
Microglia plated on cover slips were fixed for 4 min in ice-cold methanol and stained with Isolectin B4, labeled with AlexaFluor 568 (Invitrogen) for 30 min at room temperature and washed two times with PBS. Images were acquired with a × 10 objective on an Olympus IX 81 fluorescence microscope. To determine the circularity or “form factor” [7, 16], photos of three fields of view were taken per siRNA and analyzed by a Fiji/ImageJ macro to determine the circularity of each individual cell. Briefly, each image was processed by the median filter at a radius of 8 pixels, then a black and white threshold image was generated, cell surroundings were drawn, and the circularity within shape descriptions was determined as 4π*area/(perimeter)2. Cells touching the borders of the image were excluded from the quantification procedure. In total, the siRNA screen was performed three times so that the effect of each target gene on circularity was calculated from three independent experiments. Subsequently, the circularity values were normalized within each experiment and compared to microglia transfected with Scrambled siRNA. The fold change difference between the circularity value for each gene and the Scrambled siRNA was called the differential circularity, since it indicates the difference between the Scrambled siRNA and gene siRNA, and is given as logarithm (log2).
Migration of microglia was determined using a scratch assay. Cells were plated at 18,000 cells/cm2 and transfected as described above. Growth medium was replaced 3 days after transfection to CM, and the cover slip was scratched with a yellow pipette tip. Bright field images were taken 1 and 24 h after the scratch. The scratch width was analyzed with the Fiji/ImageJ macro Wound healing tool (MRI_Wound_Healing_Tool.ijm). The scratch width at 24 h was expressed as a percentage of the original scratch width at 1 h. For the normalization, the scratch width of the Scrambled siRNA-transfected cells was set to 1.0 in each experiment, and the scratch width of the siRNA-transfected cells were set to values proportional to 1.0.
Western blot analysis
Microglia were transfected as described above and then harvested with cold lysis buffer consisting in 10 mM Tris-HCl pH 8.0, 150 mM NaCl, 1% Nonidet-P40, 1 mM EDTA, 10 mMNaF, 1 mM Na3VO4, and a cocktail of commercially available protease inhibitors (Sigma, containing 104 mM AEBSF, 80 μM Aprotinin, 4 mM Bestatin, 1.4 mM E-64, 2 mM Leupeptin, 1.5 mM Pepstatin A). After centrifugation for 15 min at 14,000 rpm, the protein concentration of the supernatants was determined by Bradford assay (Bio-Rad), and samples were prepared for SDS-PAGE with Laemmli SDS sample buffer. After semi-dry blotting or Western blotting in a wet chamber for Tiam1, PVDF membranes were blocked with 5% milk in TBS-Tween (0.1%) and incubated 8 h at 4 °C with the primary antibodies, rabbit anti-Ahrgef4, (Antibodies-Online), rabbit anti-IκBα (Cell Signaling), mouse anti-Rac1 (BD Transduction laboratories), rabbit anti-GAPDH (Sigma), mouse anti-Tiam1, and mouse anti-Map3k2 (both from Santa Cruz) diluted in 2% BSA/TBS-Tween (0.1%) or incubated 1 h at RT with the primary antibodies mouse anti-Creb1, rabbit anti-RhoE, rabbit anti-GM-CSFR, and rabbit anti-Mapk11 (all from Antibodies-Online) diluted in the blocking solution. Horseradish peroxidase-conjugated secondary antibodies (DakoCytomation) and ECL (GE Biotech) were used for detection. If necessary, membranes were stripped with stripping buffer (100 mM β-mercaptoethanol, 2% SDS, 62.5 mMTris pH 6.8) for 30 min at 55 °C.
RNA extraction and RT-qPCR
Total RNA was extracted using Tripure (Roche) from microglia plated in 6-well plates. After DNase I treatment (Sigma), RNA (1 μg/sample) was reverse transcribed using RevertAid First Strand cDNA Synthesis kit (Fermentas) and random hexamer primers. The cDNA was analyzed by qPCR in triplicates on a Cfx96-Cycler (Bio-Rad) with the SensiFAST™ SYBR® No-ROX Kit (Bioline) and 2.5 pmol of the following primers: mouse Rac1 forward, CCC AAT ACT CCT ATC ATC CTC G; mouse Rac1 reverse, CAG CAG GCA TTT TCT CTT CC; mouse BDNF forward, CCC TCC CCC TTT TAA CTG AA; mouse BDNF reverse, GCC TTC ATG CAA CCG AAG TA with a primer efficiency of 100.7% and GAPDH forward and reverse primers from the RevertAid First Strand cDNA Synthesis kit (Fermentas), with a primer efficiency of 86.3%. After 42 cycles, the Ct values were determined. To normalize the samples, ΔCt between the gene of interest and GAPDH Ct values as reference gene was calculated. The x-fold difference in expression between the different treatments was then determined by subtraction of the ΔCt values and called ΔΔCt. Finally, the total change was calculated as 2−ΔΔct and the relative amount compared to Scrambled siRNA-transfected cells was deducted.
BDNF protein levels in culture supernatants of siRNA-transfected primary microglia were determined using the Quantikine® ELISA Kit (R+D Systems®), a quantitative sandwich enzyme immunoassay, according to the manufacture’s recommendations.
All data are expressed as the mean ± SEM. The data and statistical analysis comply with the recommendations on experimental design and analysis in pharmacology (Curtis et al., 2018). Statistical analysis was carried out with two-way ANOVA followed by Student’s t test. We assumed significance at p < 0.05. For the statistical analysis of the siRNA screen, cell HTS2 package version 2.38.0  from R version 3.3.2 was applied to calculate and normalize all circularity values within each experiment, as it is the case in the normalization by experiments. Finally, these normalized values were compared for each gene to a Scrambled siRNA used as negative control sample using limma software version 3.30.11 . Achieved results are presented as a “fold change” in the circularity, given as logarithm (log2) value. The threshold value between significant and not significant genes was determined by an adjusted P value: the Fold Discovery Rate (FDR). We applied the standardized FDR, as previously described .
Rac1 siRNA efficiently inhibited CM-induced ramification
First of all, we established the experimental conditions to evaluate the efficiency of the siRNA knockdown in primary mouse microglia. Since a dominant-negative mutant of the small RhoGTPase Rac1 inhibited the CM-induced conversion to an elongated neuroprotective microglia , we expected the same phenotype for its siRNA knock down. Therefore, we used Rac1 as test protein. Indeed, Rac1 protein and mRNA expression were efficiently downregulated (66.1 ± 9.0% protein and 59.7 ± 5.1% mRNA at 72 h) after transfection with Silencer-Select Rac1 siRNAs in comparison to scrambled siRNA (Additional file 1). Our initial assumption was confirmed as CM-induced cell elongation was inhibited in Rac1 siRNA-transfected microglia (Fig. 1a, b). The quantification shows that CM-treated microglia that were previously transfected with three individual siRNAs targeting Rac1 or a pool of them have a higher value for circularity, corresponding to a rounder cell shape, compared to control CM-treated microglia that were not transfected or transfected with scrambled siRNAs (Fig. 1c). Thus, these experiments show that our assay is functional and that Rac1 siRNA could be used as specific positive control in this siRNA screen.
siRNA screening in CM-treated microglia identifies various potential targets to regulate microglial activation
Using this experimental strategy, we next evaluated the effect on CM-induced microglia elongation of a custom-made Silencer-Select siRNA library against 157 mouse proteins including 100 cytoskeletal proteins and 57 regulators of microglia activation/inflammatory pathways. We selected 22 small RhoGTPases  and 78 of their activators, Rho guanine nucleotide exchange factors RhoGEFs , because they are involved in regulation of cytoskeletal architecture and changes in the cell morphology are intrinsically related to the cytoskeleton. On the other hand, we selected major regulators of inflammatory pathways, including various members of the mitogen-activated protein kinases (MAPK), phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K) and nuclear factor kappa B subunit 1 (NF-κB) signaling routes, because these proteins are likely to interfere with the CM-induced anti-inflammatory phenotype of microglia and hence have a higher probability to influence the inflammatory state of microglia. The supplementary Tables S1 and S2 in Additional files 2 and 3, respectively, show a list of all genes that were targeted by an initial screening using our siRNA library (a pool of three individual siRNAs) including their effects on the differential circularity. A positive differential circularity (Additional file 2) signifies that the cells were rounder than cells transfected with Scrambled siRNA, and negative values of circularity (Additional file 3) imply that microglia transfected with this gene showed an enhanced ramification pattern. Since we were interested in genes that inhibited the neuroprotective ramified phenotype of microglia, the statistically significant positive hits (genes marked in bold in Additional file 2) were compared to publicly available datasets of PD-associated gene expression  and pathways . We also compared them to differential transcriptome expression data from post mortem tissue of the substantia nigra between healthy and PD patients . Overlapping genes between our siRNA screen and those databases as well as the seven hits with the highest differential circularities were subjected to a literature search to collect information about their potential physiological role in microglia. From there, the selected hits Asef (Arhgef4), Tiam1, GM-CSFR (Csf2rA), IκBα (Nfκbia), RhoE (Rnd3), p38β (Mapk11), Creb1, MEK3 (Map2k3), Map3k2, and JNK2 (Mapk9) were then validated in the following step.
Since siRNA screens in mammalian cells might be influenced by off-target effects , we confirmed selected positive hits by transfecting individually the three single siRNAs of the original pool. To exclude these off-target effects, at least two siRNA oligonucleotides targeted against the same gene should exhibit the same phenotypic effect on CM-treated microglia (assayed by changes in circularity values) and efficiency to knock down its gene target (assayed by Western blot analysis). As shown in Table 1, not all positive hits of the previous siRNA screening were verified in this step. None of the three siRNAs against MEK3, and the two siRNAs against JNK2, were able to individually change circularity values in CM-treated microglia (Table 1), and thus, the genes MEK3 and JNK2 were marked as false positives. In addition, the protein levels of all tested genes, except Arhgef4, were significantly downregulated by at least two single siRNAs (Fig. 2). Since none of the individual siRNAs targeting Arhgef4 were able to downregulate the protein Arhgef4, this hit had to be considered as false positive.
Selected hits differentially regulate BDNF production and migration in microglia
To pinpoint the functional consequences of the ramification impairment, we ran secondary screens with the seven validated hits (Tiam1, GM-CSFR, IκBα, RhoE, p38β, Creb1, and Map3k2) by investigating their effects in two processes that are related to microglia activation or with mechanisms through which microglia might exert their neuroprotective functions, namely cell migration capacity and production of neurotrophic factors. Since CM-induced BDNF production in primary microglia , we transfected primary microglia with the two individual siRNAs of the seven positive hits that efficiently downregulated their target protein (see Fig. 2), and measured the expression of BDNF mRNA 4 days later. siRNAs against Creb1, RhoE, and p38β were able to downregulate basal BDNF mRNA expression by microglia (Fig. 3a), suggesting that they could achieve their neuroprotective function via BDNF production. However, siRNAs against Tiam1, GM-CSFR, IκBα, and Map3k2 failed to regulate BDNF mRNA expression (data not shown). To confirm the downregulation of BDNF mRNA at the protein level, we performed an ELISA to measure the total amount of BDNF protein in the culture supernatants of Creb1, RhoE, and p38β siRNA-transfected cells. Indeed, downregulation of these three proteins also resulted in a significant reduction of BDNF protein (Fig. 3b).
On the other hand, since the amoeboid cell shape in microglia is associated with cell migration [26, 27], we used a scratch assay to test the effect of our seven validated hits on microglia migration. We transfected primary microglia with two individual siRNAs of our positive hits, and 3 days later, we added CM and scratched the cover slip. We observed that microglia transfected with siRNAs against RhoE, Tiam1, or IκBα siRNAs closed their scratch prior to Scrambled siRNA (Fig. 4), indicating that these three genes are negative regulators for microglia migration. In contrast, siRNAs against GM-CSFR, p38β, Creb1, and Map3k2 did not affect their migratory capacity (data not shown). The downregulation of RhoE, Tiam1, or IκBα by siRNA transfection did not affect MTT reduction in microglia cultures (data not shown), suggesting that their effects measured by the scratch assay are not due to an action on cell viability/proliferation. Table 2 summarizes the results of both secondary screens.
Highly elongated and ramified microglia are predominantly found in a healthy brain . But it is commonly accepted that microglia can display an amoeboid morphology during a neuroinflammatory response , for example, as seen in PD . Although microglia ramification is equally important as the development of axons and dendrites for neurons to exert their functions, very little is known about the pathways that induce this ramification. To systematically unravel this process, we performed a microscopy-based siRNA screen in primary microglia and have identified a list of 45 genes that are positively involved in the ramification of microglia. Since we are aware of high numbers of false positive hits in siRNA screens , we have applied a rigorous validation procedure in our siRNA screen in order to make our hits reliable and reproducible results. The execution of this validation step seems to be essential, because we observed that three of ten initially selected positive hits were false positives.
Two of the seven positive hits, the Rho guanine nucleotide exchange factor (RhoGEF) Tiam1 and the small RhoGTPase RhoE (also known as Rnd3) are regulators of the cytoskeleton. The NF-κB inhibitor IκBα, the kinases p38β (also known as Mapk11) and Map3k2, the transcription factor Creb1, and the cytokine receptor for granulocyte macrophage colony stimulating factor (GM-CSFR, also known as Csf2ra) are molecules involved in inflammatory responses. Since their knock down inhibited microglia ramification upon CM treatment, a potent anti-inflammatory neuroprotective stimulus , these seven proteins must play an essential role in this process and their identification sheds the first light on these so far poorly characterized signaling routes.
The fact that the BDNF expression is affected by the downregulation of three out of seven positive hits indicates that this might be a possible mechanism for acquiring a neuroprotective phenotype. Indeed, a decreased release of BDNF by microglia is one of the mechanisms by which microglia can kill neurons . Since BDNF transcription is mediated by Creb1 , our result that BDNF mRNA expression was reduced by Creb1 siRNA was actually not surprising. Interestingly, RhoE has been identified as BDNF-induced Creb1-regulated gene in hippocampal neurons . Now, we report here for the first time that RhoE vice-versa also regulates BDNF expression, establishing a potential positive feedback similar to that previously described for Creb1 activation-BDNF expression . Because this positive feedback-loop is one of the reasons to consider to Creb1 as a key player in neuronal survival , our new findings with RhoE could result of relevance from a physiological point of view.
The MAPK p38 is expressed as two major isoforms in the brain, namely p38α and p38β, which have around 70% homology. Most of the published work so far refers to p38α and only few data is available distinguishing both isoforms, mainly because the main p38 inhibitors target both isoforms . However, a recent report showed that siRNA-mediated downregulation of p38β, but not of p38α, resulted in impaired Creb1 activation in microglia . These results are in good concordance with our results showing that p38α (Mapk14, Additional file 3) did not come up as positive hit in our primary screen for deactivating-microglia drug targets. Thus, these findings suggest a differential role of both isoforms of p38 in microglia activity states, with p38β mainly being implicated in a neuroprotective phenotype.
Besides affecting particular changes in cell morphology and gene expression, positive hits also changed microglia migration, being a more global cellular read-out. Migration is one of the important functions of microglia to assess a lesion site once injury or insult in the brain has occurred and is typically associated with an inflammatory phenotype . Three of our positive hits, namely Tiam1, RhoE, and IκBα, enhanced cell migration when downregulated by siRNA. It was previously described that downregulation of RhoE promoted glioblastoma cell migration and invasion . On the other hand, although the RhoGEF Tiam1 is traditionally believed to enhance cell migration, its role in cell adhesion makes its involvement in cell migration controversial depending on cell type and extracellular matrix used . Under our experimental conditions, Tiam1 would be rather involved in cell adhesion, since its downregulation enhances primary microglia migration. Regarding IκBα, it is likely that it exerts its neuroprotective role in microglia by inhibiting the inflammatory signaling of the transcription factor NF-κB, probably by retaining it in the cytoplasm . The role of IκBα itself in cell migration is less studied, but for example curcumin suppresses inflammatory-induced NF-κB signaling and macrophage migration . In addition, pharmacological inhibition of nuclear translocation of NF-κB blocked endothelial cell migration in a scratch assay . These results correlate with ours, namely that the downregulation of IκBα, which probably increases nuclear NF-κB, enhanced microglia migration.
Since in our secondary screens not all hits overlap, we must assume that our seven validated hits probably do not act on the same pathway, probably because the CM we applied is a mix of several growth factors and cytokines and can stimulate several pathways at once. It is however remarking that RhoE came up as hit in both analyses, placing RhoE in a central role to act on microglia activation states. RhoE is an atypical RhoGTPase, because its GTPase domain is not able to hydrolyze GTP, being its activity regulated by its expression level and phosphorylation . As mentioned above, its gene expression can be regulated by BDNF . RhoE’s main function is to act as endogenous antagonist to RhoA-mediated actin cytoskeleton remodeling, which is especially important in cell migration and polarization. For example, polarization and axonal and dendritic lengths are reduced in RhoE-null neurons , a comparable phenotype to that we observed in RhoE-downregulated microglia, which were less ramified. Taken all these data together, we suggest RhoE is a promising candidate gene, whose action is required for the reversion of microglia into a neuroprotective phenotype, being especially indicative its role in BDNF production. However, further studies are necessary to establish its precise molecular role in this process. Since RhoE, similar to p38β, has not been previously identified publicly available databases related to PD, these findings are completely novel in the context of neurodegenerative diseases.
Finally, our primary siRNA screen also identified 50 genes (Additional file 3) that potentially block the acquisition of a ramified morphology in microglia. This suggests that selective inhibitors against these molecules would favor this neuroprotective phenotype. Although the validation of these negative hits is out of the scope of the present study, it offers this possibility to other researchers.
In summary, with the first siRNA screen performed in primary microglia, we generated a validated database of intracellular molecules that play a role in the transition of microglia from an activated damaging phenotype to a potential neuroprotective one. Hence, our hits might represent possible novel drug targets in neurodegenerative diseases with an underlying neuroinflammatory process. With these important elementary results, we open up new investigation lines to reveal the mechanism involved in this process and generate knowledge in order to discover innovative treatments against neurodegenerative disorders by preventing neuronal death mediated by over-activated microglia.
Adipose-derived mesenchymal stem cells
Brain-derived neurotrophic factor
Conditioned medium of ASCs
Central nervous system
Fold discovery rate
Receptor for granulocyte macrophage colony stimulating factor
Inhibitor of NF-κB
Mitogen-activated protein kinases
Nuclear factor kappa B subunit 1
Rho guanine nucleotide exchange factor
Small interference RNA
T-lymphoma invasion and metastasis-inducing protein 1
Qian L, Flood PM, Hong JS. Neuroinflammation is a key player in Parkinson’s disease and a prime target for therapy. J Neural Transm. 2010;117(8):971–9.
Gonzalez H, Elgueta D, Montoya A, Pacheco R. Neuroimmune regulation of microglial activity involved in neuroinflammation and neurodegenerative diseases. J Neuroimmunol. 2014;274(1–2):1–13.
Wolf SA, Boddeke HW, Kettenmann H. Microglia in physiology and disease. Annu Rev Physiol. 2017;79:619–43.
Doorn KJ, Goudriaan A, Blits-Huizinga C, Bol JG, Rozemuller AJ, Hoogland PV, et al. Increased amoeboid microglial density in the olfactory bulb of Parkinson’s and Alzheimer’s patients. Brain Pathol. 2014;24(2):152–65.
Block ML, Zecca L, Hong JS. Microglia-mediated neurotoxicity: uncovering the molecular mechanisms. Nat Rev Neurosci. 2007;8(1):57–69.
Kettenmann H, Hanisch UK, Noda M, Verkhratsky A. Physiology of microglia. Physiol Rev. 2011;91(2):461–553.
Neubrand VE, Pedreno M, Caro M, Forte-Lago I, Delgado M, Gonzalez-Rey E. Mesenchymal stem cells induce the ramification of microglia via the small RhoGTPases Cdc42 and Rac1. Glia. 2014;62(12):1932–42.
Hanisch UK, Kettenmann H. Microglia: active sensor and versatile effector cells in the normal and pathologic brain. Nat Neurosci. 2007;10(11):1387–94.
Nobes CD, Hall A. Rho, rac, and cdc42 GTPases regulate the assembly of multimolecular focal complexes associated with actin stress fibers, lamellipodia, and filopodia. Cell. 1995;81(1):53–62.
Bishop AL, Hall A. Rho GTPases and their effector proteins. Biochem J. 2000;348(Pt 2):241–55.
Krausz E. High-content siRNA screening. Mol BioSyst. 2007;3(4):232–40.
Demir K, Boutros M. Cell perturbation screens for target identification by RNAi. Methods Mol Biol. 2012;910:1–13.
Horvath RJ, Nutile-McMenemy N, Alkaitis MS, Deleo JA. Differential migration, LPS-induced cytokine, chemokine, and NO expression in immortalized BV-2 and HAPI cell lines and primary microglial cultures. J Neurochem. 2008;107(2):557–69.
Kilkenny C, Browne W, Cuthill IC, Emerson M, Altman DG, Group NCRRGW. Animal research: reporting in vivo experiments: the ARRIVE guidelines. Br J Pharmacol. 2010;160(7):1577–9.
Anderson P, Souza-Moreira L, Morell M, Caro M, O'Valle F, Gonzalez-Rey E, et al. Adipose-derived mesenchymal stromal cells induce immunomodulatory macrophages which protect from experimental colitis and sepsis. Gut. 2013;62(8):1131–41.
Wilms H, Hartmann D, Sievers J. Ramification of microglia, monocytes and macrophages in vitro: influences of various epithelial and mesenchymal cells and their conditioned media. Cell Tissue Res. 1997;287(3):447–58.
Boutros M, Bras LP, Huber W. Analysis of cell-based RNAi screens. Genome Biol. 2006;7(7):R66.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article 3.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57(1):289–300.
Bustelo XR, Sauzeau V, Berenjeno IM. GTP-binding proteins of the Rho/Rac family: regulation, effectors and functions in vivo. BioEssays. 2007;29(4):356–70.
Rossman KL, Der CJ, Sondek J. GEF means go: turning on RHO GTPases with guanine nucleotide-exchange factors. Nat Rev Mol Cell Biol. 2005;6(2):167–80.
Foulger RE, Denny P, Hardy J, Martin MJ, Sawford T, Lovering RC. Using the gene ontology to annotate key players in Parkinson’s disease. Neuroinformatics. 2016;14(3):297–304.
Fujita KA, Ostaszewski M, Matsuoka Y, Ghosh S, Glaab E, Trefois C, et al. Integrating pathways of Parkinson’s disease in a molecular interaction map. Mol Neurobiol. 2014;49(1):88–102.
Glaab E, Schneider R. Comparative pathway and network analysis of brain transcriptome changes during adult aging and in Parkinson’s disease. Neurobiol Dis. 2015;74:1–13.
Singh S, Narang AS, Mahato RI. Subcellular fate and off-target effects of siRNA, shRNA, and miRNA. Pharm Res. 2011;28(12):2996–3015.
Jonas RA, Yuan TF, Liang YX, Jonas JB, Tay DK, Ellis-Behnke RG. The spider effect: morphological and orienting classification of microglia in response to stimuli in vivo. PLoS One. 2012;7(2):e30763.
Brown GC, Vilalta A. How microglia kill neurons. Brain Res. 2015;1628(Pt B):288–97.
Birmingham A, Selfors LM, Forster T, Wrobel D, Kennedy CJ, Shanks E, et al. Statistical methods for analysis of high-throughput RNA interference screens. Nat Methods. 2009;6(8):569–75.
West AE, Griffith EC, Greenberg ME. Regulation of transcription factors by neuronal activity. Nat Rev Neurosci. 2002;3(12):921–31.
Lesiak A, Pelz C, Ando H, Zhu M, Davare M, Lambert TJ, et al. A genome-wide screen of CREB occupancy identifies the RhoA inhibitors Par6C and Rnd3 as regulators of BDNF-induced synaptogenesis. PLoS One. 2013;8(6):e64658.
Walton MR, Dragunow I. Is CREB a key to neuronal survival? Trends Neurosci. 2000;23(2):48–53.
Yasuda S, Sugiura H, Tanaka H, Takigami S, Yamagata K. p38 MAP kinase inhibitors as potential therapeutic drugs for neural diseases. Cent Nerv Syst Agents Med Chem. 2011;11(1):45–59.
Wu HY, Mao XF, Fan H, Wang YX. p38beta mitogen-activated protein kinase signaling mediates exenatide-stimulated microglial beta-endorphin expression. Mol Pharmacol. 2017;91(5):451–63.
Liu B, Dong H, Lin X, Yang X, Yue X, Yang J, et al. RND3 promotes Snail 1 protein degradation and inhibits glioblastoma cell migration and invasion. Oncotarget. 2016;7(50):82411–23.
Boissier P, Huynh-Do U. The guanine nucleotide exchange factor Tiam1: a Janus-faced molecule in cellular signaling. Cell Signal. 2014;26(3):483–91.
Li Q, Verma IM. NF-kappaB regulation in the immune system. Nat Rev Immunol. 2002;2(10):725–34.
Young NA, Bruss MS, Gardner M, Willis WL, Mo X, Valiente GR, et al. Oral administration of nano-emulsion curcumin in mice suppresses inflammatory-induced NFkB signaling and macrophage migration. PLoS One. 2014;9(11):e111559.
Liu C, Tsai AL, Li PC, Huang CW, Wu CC. Endothelial differentiation of bone marrow mesenchyme stem cells applicable to hypoxia and increased migration through Akt and NFkappaB signals. Stem Cell Res Ther. 2017;8(1):29.
Jie W, Andrade KC, Lin X, Yang X, Yue X, Chang J. Pathophysiological functions of Rnd3/RhoE. Compr Physiol. 2015;6(1):169–86.
Peris B, Gonzalez-Granero S, Ballester-Lurbe B, Garcia-Verdugo JM, Perez-Roger I, Guerri C, et al. Neuronal polarization is impaired in mice lacking RhoE expression. J Neurochem. 2012;121(6):903–14.
We thank Dr. Eduardo Andres-Leon and Dr. Jenny Campos-Salinas (IPBLN, CSIC) for statistical analyses and technical help in cell shape determination, respectively. We also thank Dr. Rainer Pepperkok and the Advanced Light Microscopy Facility (ALMF) at EMBL Heidelberg for technical advice.
This work was financed by the Michael J Fox Foundation (MJFF grant ID 9404 to VEN and MD) and an Excellence Grant of the Andalusian Government (to MD). Both funding bodies supported the collection, analysis, and interpretation of data.
Availability of data and materials
The datasets analyzed during the current study are all available on the PD map webpage of the University of Luxembourg, https://pdmap.uni.lu/minerva/. In detail, we refer to the Parkinson’s UK Gene Ontology Project (Foulger et al., 2016), the PD map tool of the University of Luxembourg (Fujita et al., 2014), and the differential transcriptome expression data from post mortem tissue of the substantia nigra (Glaab &Schneider, 2015), as referenced here:
Foulger RE, Denny P, Hardy J, Martin MJ, Sawford T, Lovering RC. Using the Gene Ontology to Annotate Key Players in Parkinson’s Disease. Neuroinformatics. 2016;14(3):297–304. doi:https://doi.org/10.1007/s12021-015-9293-2
Fujita KA, Ostaszewski M, Matsuoka Y, Ghosh S, Glaab E, Trefois C, et al. Integrating pathways of Parkinson’s disease in a molecular interaction map. Molecular Neurobiology. 2014;49(1):88–102. doi:1007/s12035-013-8489-4
Glaab E, Schneider R. Comparative pathway and network analysis of brain transcriptome changes during adult aging and in Parkinson’s disease. Neurobiology of Disease. 2015;74:1–13. doi:https://doi.org/10.1016/j.nbd.2014.11.002
Ethics approval and consent to participate
The experimental protocols of this study conform to EU Directive 2010/63 and followed the ethical guidelines for investigations with experimental animals approved by the Ethics Review Committee for Animal Experimentation of Spanish Council of Scientific Research.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rac1 siRNA transfection efficiently downregulated Rac1 protein and gene in microglia. (A) Primary microglia (140,000 cells/well, 12-well plates) were transfected with 50 pmol Rac1 of Silencer Select Rac1 siRNA or Scrambled siRNA (used as control) and 1.5 μl Lipofectamine 3000. Rac1 protein expression was determined by Western blot analysis 96 h later. (B) and (C) Primary microglia were transfected as in (A) for the indicated time periods, and Rac1 downregulation was evaluated by RT-qPCR (B) and by Western blot (C). Data are mean ± SEM of three independent experiments.**p < 0.01, ***p < 0.001* vs. Scrambled siRNA. (JPG 301 kb)
Table S1. List of siRNA targets resulting in positive differential circularity. Primary microglia were transfected with a pool of siRNA oligonucleotides against each gene and tested for their effect on CM-induced microglia ramification. The result was expressed as differential circularity. siRNAs that caused positive differential circularity values were selected and sorted by their fold discovery rate (FDR), which represents an adjusted P value. Statistically significant genes were marked in bold. Genes that were annotated by the Parkinson’s UK Gene Ontology Project  are marked in the PD gene column. Genes that were listed in the PD map tool of the University of Luxembourg (http://pdmap.uni.lu/minerva/) are marked in the PD map column . Genes overlapping with differential transcriptome expression data from post mortem tissue of the substantia nigra  published on the PD map webpage (http://pdmap.uni.lu/minerva/) are marked in the last column. *Positive hits that were selected for validation. (XLSX 55 kb)
Table S2. List of siRNA targets resulting in negative differential circularity. Primary microglia were transfected with a pool of siRNA oligonucleotides against each gene and tested for their effect on CM-induced microglia ramification, which was expressed as differential circularity. siRNAs that resulted in negative differential circularity values were selected and sorted by their fold discovery rate (FDR). Statistically significant genes were marked in bold. Genes that were annotated by the Parkinson’s UK Gene Ontology Project  are marked in the PD genes column. Genes that were listed in the PD map tool of the University of Luxembourg (http://pdmap.uni.lu/minerva/) are marked in the PD map column . Genes overlapping with differential transcriptome expression data from post mortem tissue of the substantia nigra  published on the PD map webpage (http://pdmap.uni.lu/minerva/) are marked in the last column. (XLSX 56 kb)
About this article
Cite this article
Neubrand, V.E., Forte-Lago, I., Caro, M. et al. The atypical RhoGTPase RhoE/Rnd3 is a key molecule to acquire a neuroprotective phenotype in microglia. J Neuroinflammation 15, 343 (2018). https://doi.org/10.1186/s12974-018-1386-z