ORIGINAL RESEARCH article
Sec. Neuroplasticity and Development
Volume 10 - 2017 | https://doi.org/10.3389/fnmol.2017.00278
Glutaminase and MMP-9 Downregulation in Cortex and Hippocampus of LPA1 Receptor Null Mice Correlate with Altered Dendritic Spine Plasticity
- 1Canceromics Laboratory, Departamento de Biología Molecular y Bioquímica, Facultad de Ciencias, Universidad de Málaga, Campus de Teatinos, Málaga, Spain
- 2Unidad de Gestión Clínica de Salud Mental, Instituto de Investigación Biomédica de Málaga (IBIMA), Hospital Regional Universitario de Málaga, Málaga, Spain
- 3Unidad de Gestión Clínica de Neurociencias, Instituto de Investigación Biomédica de Málaga (IBIMA), Hospital Regional Universitario de Málaga, Málaga, Spain
- 4Departamento de Biología Celular, Genética y Fisiología, Facultad de Ciencias, Instituto de Investigación Biomédica de Málaga (IBIMA), Centro de Investigación Biomédica en Red sobre Enfermedades Neurodegenerativas (CIBERNED), Universidad de Málaga, Campus de Teatinos, Málaga, Spain
- 5Departamento de Farmacología, Facultad de Medicina, Universidad Complutense, Instituto de Investigación Sanitaria Hospital 12 de Octubre, Madrid, Spain
Lysophosphatidic acid (LPA) is an extracellular lipid mediator that regulates nervous system development and functions acting through G protein-coupled receptors (GPCRs). Here we explore the crosstalk between LPA1 receptor and glutamatergic transmission by examining expression of glutaminase (GA) isoforms in different brain areas isolated from wild-type (WT) and KOLPA1 mice. Silencing of LPA1 receptor induced a severe down-regulation of Gls-encoded long glutaminase protein variant (KGA) (glutaminase gene encoding the kidney-type isoforms, GLS) protein expression in several brain regions, particularly in brain cortex and hippocampus. Immunohistochemical assessment of protein levels for the second type of glutaminase (GA) isoform, glutaminase gene encoding the liver-type isoforms (GLS2), did not detect substantial differences with regard to WT animals. The regional mRNA levels of GLS were determined by real time RT-PCR and did not show significant variations, except for prefrontal and motor cortex values which clearly diminished in KO mice. Total GA activity was also significantly reduced in prefrontal and motor cortex, but remained essentially unchanged in the hippocampus and rest of brain regions examined, suggesting activation of genetic compensatory mechanisms and/or post-translational modifications to compensate for KGA protein deficit. Remarkably, Golgi staining of hippocampal regions showed an altered morphology of glutamatergic pyramidal cells dendritic spines towards a less mature filopodia-like phenotype, as compared with WT littermates. This structural change correlated with a strong decrease of active matrix-metalloproteinase (MMP) 9 in cerebral cortex and hippocampus of KOLPA1 mice. Taken together, these results demonstrate that LPA signaling through LPA1 influence expression of the main isoenzyme of glutamate biosynthesis with strong repercussions on dendritic spines maturation, which may partially explain the cognitive and learning defects previously reported for this colony of KOLPA1 mice.
Lysophosphatidic acid (LPA; 1 or 2-acyl-sn-glycero-3-phosphate) is a naturally occurring lysophospholipid derived from cell membranes that belongs to the growing list of bioactive lipids. It may act as an intercellular signaling molecule after recognition by G protein-coupled receptor(s) (GPCRs). In mammals, LPA is enriched in brain tissue, produced by platelets and cells from the nervous system and is also found in serum (Choi and Chun, 2013; Yung et al., 2015). Six types of LPA receptors (LPARs) have been identified so far named LPA1-LPA6, with different expression patterns and involved in different functions (Chun et al., 2010; Mutoh et al., 2012; Yung et al., 2014).
Multiple cellular responses have been reported for this bioactive natural lipid in the central nervous system (CNS) as well as in neural cell lines. LPA elicited sharp changes in cell’s shape and morphology, growth cone collapse and neurite retraction in many neuronal cell lines through activation of the small GTPase Rho which in turn contracts cytoskeleton (Ye et al., 2002; Choi et al., 2010; Fukushima et al., 2011). Furthermore, LPA also affect cell motility in glioma cell lines through Rho/Rho-associated kinases (ROKC) signaling cascades, which may contribute to their invasive phenotype in vivo (Manning et al., 2000). In cultured astrocytes, LPA induces a plethora of responses including proliferation and inhibition of glutamate uptake (Steiner et al., 2002; Shano et al., 2008). Of note, LPA-primed astrocytes promote neuronal differentiation of cerebral cortical progenitors and developing cortical neurons, which showed increases in arborization and neurite outgrowth (Spohr et al., 2008, 2011).
Although many types of LPARs have been detected in brain, LPA1 is abundantly expressed and considered the most prevalent receptor type in both embryonic and adult brains of humans and mice (Hecht et al., 1996; Choi et al., 2010; Mutoh et al., 2012). Genetic silencing of LPA1 in mice causes a reduced ventricular zone along with loss of cortical layer cellularity (Estivill-Torrús et al., 2008), reduces neurogenesis in the dentate gyrus (DG; Matas-Rico et al., 2008), and shows altered neurotransmitter homeostasis (Musazzi et al., 2011; Blanco et al., 2012a) which were previously related to psychiatric diseases (Harrison et al., 2003; Roberts et al., 2005). Behavioral studies with LPA1 deficient mice reported deficiencies in spatial memory retention and abnormal use of searching orientation strategies (Santin et al., 2009), defective working and reference memory independently of exploratory and emotional impairments attributed to hippocampal malfunction (Castilla-Ortega et al., 2010). Lastly, mice lacking LPA1 receptor exhibit an endophenotype for alcohol preference associated with glutamate receptor alterations in the prefrontal cortex (PFC), including a decrease in their relative mRNA levels of Gls-encoded long glutaminase protein variant (KGA; Castilla-Ortega et al., 2016).
Studies dealing with KO models of LPARs are giving new insights into the role of LPA in developmental and differentiation processes of neurons and glial cells and their cross-talk in the tripartite synapsis; however, little is known about how LPA signaling regulates synaptic function and neurotransmission (Ye et al., 2002; García-Morales et al., 2015). LPA may modify N-methyl-D-aspartate (NMDA) receptor functions in hippocampal neurons (Lu et al., 1999) as well as calcium intracellular levels (Holtsberg et al., 1997). Furthermore, LPA inhibited Na+, K+-ATPase activity in rat cortical synaptosomes (Nishikawa et al., 1989) and was related with synaptic vesicle formation (Schmidt et al., 1999). In addition, LPA through LPA1 receptor has been proposed as a potential candidate to control short-term synaptic plasticity both in excitatory and inhibitory synapses, but using different mechanisms: decreasing the synaptic vesicle pool and reducing the number of postsynaptic receptors, respectively (García-Morales et al., 2015). The aim of this study was to elucidate further relationships between LPA signaling and synaptic plasticity of glutamatergic excitatory transmission using KOLPA1 mice as a model. We found marked regional down-regulations of the main glutaminase (GA; EC 22.214.171.124) isoform (KGA) involved in the synthesis of neurotransmitter glutamate, altered dendritic spines morphology at the glutamatergic pyramidal cells of the hippocampus and a decrease in MMP-9 activity and expression, an enzyme crucially required for dendritic spine maturation and remodeling. Altogether, the results point to important implications of LPA in synaptic excitatory plasticity which may contribute to the cognitive and memory deficits shown by LPA1-deficient mice.
Materials and Methods
The generation and characterization of maLPA1-null mice has been described previously (Estivill-Torrús et al., 2008). An LPA1-null mouse colony, termed maLPA1 after the Málaga variant of the LPA1 knockout, was spontaneously derived during the original colony (Contos et al., 2000) expansion by crossing heterozygous (HT) foundation parents (maintained in the original hybrid C57BL/6J6x1291/SvJ background). Intercrosses were performed with these mice and were subsequently backcrossed for 20 generations with mice generated within this mixed background. MaLPA1-null mice carrying the LPA1 deletion were born at the expected Mendelian ratio and survived to adulthood. Targeted disruption of the Lpa1 gene was confirmed by genotyping (Contos et al., 2000), and immunochemistry confirmed the absence of LPA1 protein expression. Three month-old male littermates from wild-type (WT), heterozygous and malLPA1-null homozygous (KOLPA1) genotypes were employed for immunohistochemistry. In all other experiments, WT and KOLPA1 genotypes were employed. Mice were singly housed on a 12-h light/dark cycle (lights on at 07:00 a.m.). Water and food were provided ad libitum. All animal experiments were performed in accordance with the animal research regulations (RD53/2013 and 2010/63/UE) from Spain and European Union, and with the approval of the Committees of Animal Research from the University of Malaga (Spain).
Immunohistochemistry of Glutaminase Isoforms
For immunohistochemical studies animals (n = 8 mice per genotype) were deeply anesthetized with sodium pentobarbital (60 mg/kg) and perfused transcardially with 0.1 M phosphate buffer (PB) saline (PBS), pH 7.4, followed by fixative solution (4% paraformaldehyde, 75 mM lysine and 10 mM sodium metaperiodate in 0.1 M PB, pH 7.4). After overnight postfixation in the same fixative at 4°C, brains were cryoprotected (30% sucrose), sectioned (40 μm) on a freezing microtome and serially collected (six series for each brain, with a separation of 240 μm between consecutive sections) in wells containing cold PBS and 0.02% sodim azide. The rabbit polyclonal isoform-specific antibodies anti-KGA and anti-GAB (glutaminase gene encoding the liver-type isoforms (Gls2)-encoded long GA isoform) were generated and characterized in our laboratory (Olalla et al., 2002, 2008; Campos et al., 2003; Cardona et al., 2015). Briefly, free-floating sections, after blocking the endogenous peroxidase activity (3% H2O2/10% methanol in 0.1 M PBS pH 7.4 for 20 min) and avidin, biotin, and biotin-binding proteins for 30 min with avidin-biotin blocking Kit (Vector Labs, United Kingdom), were incubated with the affinity-purified rabbit polyclonal anti-KGA antibody (1:1000) or rabbit polyclonal anti-GAB antibody (1:500 dilution) for 48 h at room temperature. The tissue-bound primary antibody was then detected by incubating with biotinylated anti-rabbit IgG (Vector Labs, United Kingdom) for 1 h at room temperature (1:500 dilution) and then with ExtrAvidin-peroxidase conjugate (Sigma-Aldrich, Spain) for 1 h at room temperature (1:2000 dilution). Immunoreaction product was visualized with 0.05% 3-3′-diaminobenzidine tetrahydrochloride (DAB; Sigma-Aldrich, Spain), 0.03% nickel ammonium sulfate, and 0.01% H2O2 in PBS. Sections were then mounted onto gelatin-coated slides, dehydrated in graded ethanol, cleared in xylene and coverslipped with DPX mounting medium (VWR BDH Prolabo, VWR International). Immunostaining was observed under a Nikon Eclipse 80i microscope and images acquired with a Nikon DS-5M high-resolution digital camera. All immunohistochemical studies by light microscopy were tested using negative controls (omitting the primary antisera) and no immunoreaction products were detected in any case.
Golgi staining with FD Rapid GolgiStain™ kit (Golgi-Cox technique; FD Neurotechnologies) was performed to detect morphological alterations in neuronal dendrites and dendritic spines. The kit was used following the provider instructions. Briefly, adult male mice from both genotypes, WT and KOLPA1 (n = 5 per group), were deeply anesthetized with sodium pentobarbital (60 mg/kg) and brains were immediately removed. Tissues were rinsed in distilled water and immersed in the fixative solution (A + B) during 2 weeks in the dark. Afterwards, brains were transferred to solution C for 4 days. Before sectioning on a freezing microtome (Leica, CM 1325), tissue was rapidly frozen in isopentane precooled to −70°C. Sections of 150 μm were mounted with solution C on gelatine-coated slides and air-dried at RT. Then, the kit staining procedure was followed, and slides were dehydrated in graded ethanol, cleared in xylene and coverslipped with DPX (BDH) mounting medium.
Microphotographs of basal pyramidal dendrites in CA1 stratum oriens from hippocampus were acquired using a 63× oil-immersion objective of a SP5 II microscope (Leica) with a 405 Diode, UV Laser. Series of Z-stacks of dendrites were taken with 0.3 μm spacing. LAS AF Lite (Leica) software was used for quantification and morphometric analysis. For quantitative analysis, we assessed nonprimary dendritic segments, using the following criteria: (1) the quality of Golgi staining; (2) the relative isolation of the dendrites; (3) we selected only the spines that protruded in the transverse direction and were clearly distinguished. Spine density (number of spines per 10 μm), as well as dendrite diameter at several points of each counted segment, were measured. Other parameters such as spine morphology and length were also evaluated. Protrusions from dendrites were classified based on predefined parameters as mushroom spines, stubby spines, thin spines, or filopodia (Ghani et al., 2017). A total of 527 dendritic spines (283 spines for WT and 244 for KOLPA1 mice) were analyzed. Finally, mushroom head size was assessed, defining the head width (μm) as the diameter of the largest spine section that was perpendicular to a virtual axis of the spine trajectory (Ruszczycki et al., 2012).
Real-Time Quantitative RT-PCR
For qRT-PCR analyses (and enzymatic assays), animals were sacrificed by decapitation. Brains were quickly removed, frozen (−80°C) and dissected in coronal brain slices (2 mm thickness) with razor blades in a mouse brain slicer matrix (Zivic Instruments). The discrete brain regions (striatum, hippocampus, PFC and cerebellum) were picked up by free hand dissection using a scalpel. These brain regions were identified according to Mouse Brain in Stereotaxic Coordinates (Paxinos and Franklin, 2001) as described previously (Blanco et al., 2012a): striatum from the Bregma 1.54 mm to −0.46 mm, hippocampus from the Bregma −1.22 mm to −3.52 mm, PFC from the Bregma 2.46 mm to 1.34 mm and cerebellum from the Bregma −5.52 mm to −7.80 mm.
Total RNA from brain regions was isolated according to manufacturer’s specifications (AllPrep DNA/RNA/Protein Mini Kit, Qiagen). Total RNA (1 μg) were reverse transcribed using the Quantitec Reverse Transcription kit (Qiagen) in a 20 μl reaction, according to manufacturer’s instructions. RNA samples were tested for genomic DNA contamination by including reverse transcriptase (RT) omitted controls, and for reagents and aerosol contamination by including two non-template controls. Real time PCR was carried out with the CFX thermocycler (Bio-Rad) as described elsewhere (Martín-Rufián et al., 2012). GLS transcripts were amplified using mouse GLS-specific primers: GLS forward (5′-GCGAGGGCAAGGAGATGGTG-3′) and GLS reverse (5′-CTCTTTCAACCTGGGATCAGA TGTTC-3′), as described previously (Martín-Rufián et al., 2012). Results were evaluated using CFX Bio-Rad Analysis Software. Amplicon size and the absence of non-specific products were confirmed by agarose gel electrophoresis. For absolute quantification of GA transcripts, PCR transcript-specific amplicons were prepared as standards. PCR products were analyzed by 2% agarose gel electrophoresis and purified (Illustra™ GFXTMPCR DNA and Gel Band Purification Kit, GE Healthcare). The standard curve was created by plotting the threshold cycle (Ct) corresponding to each standard, vs. the value of their corresponding log of starting quantity (DNA copy number). GLS transcripts were expressed as numbers of molecules (± SEM) of mRNA per ng of total RNA. Five animals of each genotype (WT and KOLPA1) were employed and independent experiments were performed in triplicate.
Glutaminase Activity Assay
Brain regions were isolated as described above. The discrete brain regions from dissections were resuspended in TES buffer (25 mM Tris–HCl, 0.2 mM EDTA, 0.33 M sucrose, pH 8.0) containing the complete protease inhibitor cocktail (Roche), homogenized with a Pestle motor mixer (VWR) and solubilized with TX-100 at a final concentration of 1% (v/v) for 1 h at 4°C. After centrifugation at 100,000× g for 30 min at 4°C, the supernatants were divided into aliquots and kept at −80°C until analysis. The protein content in each sample was determined by the Bradford method. GA activity was assayed by measuring the ammonia produced in the catalytic reaction as described by Heini et al. (1987). Briefly, samples of 25 μL were added to 35 μL of a mixture of 100 mM potassium phosphate, 171 mM L-Gln and 1.5 mM NH4Cl, pH 8.0 and incubated for 1 h at 37°C. The reaction was terminated by adding 10 μL of trichloroacetic acid (TCA) 10% (w/v), kept on ice for 15 min and centrifuged at 12,000× g. Aliquots of 5 μL were then mixed with 150 μL of o-phthalaldehyde/mercaptoethanol reagent (10 mL 0.2 M potassium phosphate, pH 7.4, 0.56 mL of 72 mM mercaptoethanol in ethanol and 0.56 mL of 186 mM o-phthalaldehyde in ethanol). The samples were kept at room temperature in the dark and their absorbance measured at 410 nm after 45 min together with an NH4Cl standard. For blanks, the samples and substrate solution were incubated separately and mixed after the addition of TCA.
Gel Zymography of MMP-9 and MMP-2 Activity
Brain regions were homogenized in nonyl phenoxypolyethoxylethanol (NP-40, Thermo Fisher) lysis buffer (150 mM NaCl, 50 mM Tris-HCl, 1% NP-40, pH 7.4) supplemented by 5% protease inhibitor cocktail (Sigma-Aldrich, Spain) and 1% phosphatase inhibitor cocktail-2 (Sigma-Aldrich, Spain). Protein concentrations were measured by DC Protein Assay kit (Bio-Rad, Spain). Samples (20 μg of total protein) were mixed with 4× non reducing loading buffer, heated 15 min at 37°C and subjected to 9% SDS-PAGE gel with 0.1% gelatin (w/v) (Sigma, Spain). Gels were washed twice for 30 min at room temperature with slight shaking in enzymatic activation buffer (50 mM Tris-HCl, 6 mM CaCl2, 1.5 μM ZnCl2, pH 7.4) with 2.5% Tween-20 to remove SDS and restore gelatinase activity. Gels were then incubated at 37°C for 48 h in enzymatic activation buffer without Tween-20. Gels were stained in Coomasie brilliant blue (Bio-Rad, Spain) until gels were dark blue colored, followed by immersion in discoloring solution (40% methanol, 10% acetic acid) until bands were clear. Gels were digitized and bands were quantified by ImageJ Software version 1.42q (NIH, Bethesda, MD, USA).
Gel Zymography of Tissue Plasminogen Activator (tPA) Activity
Tissue plasminogen activator (tPA) was measured because of its role on MPP-9 proteolytic activation. tPA activity was assessed in the same samples mentioned above. Samples were mixed with 4× non reducing loading buffer, heated 15 min at 37°C and subjected to 9% SDS-PAGE gel with 0.1% plasminogen (Millipore, Spain) and 0.1% casein (Sigma-Aldrich, Spain). Gels were washed twice for 30 min at room temperature with slight shaking in modified enzymatic activation buffer (0.1 M Tris-HCl, 0.1 M glycine, 0.01 M EDTA, pH 8.0) with 2.5% Tween-20 to remove SDS and restore gelatinase activity. Gels were then incubated at 37°C for 48 h in modified enzymatic activation buffer without Tween-20. tPA activity bands were developed and quantified in the same way that we mentioned above.
To measure levels of GA and MMP-9 proteins in each brain region, Western blotting was performed with purified specific polyclonal antibodies against KGA and GAB, as well as commercial antibodies against MMP-9 (Millipore, Spain AB19016, 1/1000). Protein samples (40 μg for KGA and GAB or 20 μg for MMP-9 detection) were separated on SDS-PAGE gels and transferred to nitrocellulose or PVDF membranes. After blocking at room temperature for 1 h with 5% bovine serum albumin (BSA) in TBST buffer (0.1% Tween 20 in Tris-buffered saline, TBS, 100 mM Tris, 150 mM NaCl, pH 7.5), membranes were incubated with primary antibodies in blocking buffer with 5% BSA, overnight at 4°C. After incubation with secondary antibody, the blots were developed by enhanced chemiluminescense technique as recommended by the supplier (Pierce) and bands quantified using the Chemi Doc System (BioRad). β-actin was quantified and used as a loading control.
All statistical analyses were conducted using two-tailed unpaired t-tests (WT vs. KOLPA1) except when the assumptions of normality and/or homoscedasticicity were violated (Golgi staining). In the latter case, non-parametric Mann-Whitney tests were carried out to compare WT and KOLPA1 genotypes. The significance level was set at 5%. GraphPad Prism 7 software (La Jolla, CA, USA) was used to perform all statistical analyses.
Silencing of LPA1 Sharply Decreases KGA, but Not GLS2, Protein Expression
Four different GA isoenzymes are expressed in mammalian brain (Márquez et al., 2016): KGA and Gls-encoded short glutaminase protein variant (GAC) proteins are encoded by the Gls gene (GLS isoforms), whilst GAB and Gls2-encoded short glutaminase protein variant (LGA) arise from the Gls2 gene (GLS2 isoforms; Aledo et al., 2000; Matés et al., 2013). The influence of silencing LPA1 receptor on brain GA protein expression was assessed by immunohistochemistry of the KGA and GLS2 (GAB/LGA) isoforms. The isoform-specific anti-KGA antibodies employed were raised for a unique C-terminal region of this isoform, which is absent in the GAC isoform, and hence they cannot recognize this second GLS isoform (GAC); however, GAB antibodies would recognize both GLS2 isoforms (GAB and LGA). The immunostaining obtained for KGA (Figure 1) revealed a cytoplasmic punctate pattern typical for a mitochondrial localization, in accordance with previous reports on KGA localization in mammalian brain (Olalla et al., 2002; Cardona et al., 2015). However, a drastic decrease in the intensity of the staining was clearly seen in several cerebral regions from HT and KO mice, like cortex (Figures 1B,C, details in Figures 1E,F, respectively) and hippocampus (Figures 1H,I, respectively) in comparison to the same areas in WT mice (Figures 1A,D,G). Specifically, in the hippocampus, CA3 pyramidal layer neurons (Figures 1J3–L3) and the neuropile of both CA1 stratum lacunosum-moleculare and DG molecular layer (Figures 1J1–L1) were among the most affected. In addition, neuronal KGA-immunolabeling in certain cortical regions, like deep layers of agranular insular cortex (Figures 1D–F), was also notably reduced, whereas in others, as the cerebellum, non-specific alteration was observed (Figures 1P–R). Interestingly, these reductions in the KGA labeling were generally gradual from control to HT and KO animals, which seem to reflect a proportional decline on KGA protein expression when LPA1 receptor levels were diminished. Also progressive was the appearance of intense perivascular immunoreactivity in animals with decreased LPA1 expression, especially in the KO mice (Figures 1C,F,M–O).
Figure 1. Comparative analysis of immunostaining for Gls-encoded long glutaminase protein variant (KGA) in wild-type (WT), heterozygous (HT) and KOLPA1 mice brain. Specific immunolabeling for KGA was detected within the somata/neuropile of several brain regions in WT brain (left column). In general, there was a strong reduction in the staining together with a progressive switch to a perivascular location from control to HT and KO mice (A–C), panoramic views of sections containing cerebral cortex and striatum, black arrows indicate KGA-positive blood vessels). (D–F) Progressive change of KGA-staining in deep layers of agranular insular cortex (square; detail of cortical somatic staining, inset in (D); detail of blood vessel, inset in E). (G–I) KGA-positive staining in the CA1-CA3 hippocampal subfields and DG was dramatically decreased in both HT and KO genotypes in comparison to WT. (J–L) DG and slm neuropile stainings were found to be very affected in HT (K1) and KOLPA1 (L1) mice in comparison to WT (J1), as well as CA3 principal neurons (J2–L2), and to a lesser extent equivalent CA1 neurons (J3–L3). (M–O) Ventrally, in the hypothalamic levels, we found good examples of this gradual change to perivascular staining around the 3V. No differences were detected among the different genotypes in cerebellum (P–R). Ctx, cortex; St, striatum; bv, blood vessel; CA1, CA3, hippocampal subfields; DG, dentate gyrus; Th, thalamus; so, stratum oriens; sp, stratum pyramidale; sr, stratum radiatum; slm, stratum lacunosum-moleculare; m, molecular layer; mo, outer molecular layer; mi, inner molecular layer; h, hilus; g, granular layer; Pe, periventricular hypothalamic nucleus; 3V, third ventricle; G, cerebelar granular layer; M, cerebelar molecular layer. Scale bar: A–C, 1 mm; D–F, 100 μm; G–I, 500 μm; J1–L1, 50 μm; J2–L2 and J3–L3, 20 μm; M–R, 200 μm; inset: D, 25 μm; inset: E, 40 μm.
In contrast with KGA, the immunolabeling patterns against GLS2 isoforms were roughly similar for the three types of mouse strains studied (Figure 2). No apparent change on GLS2 expression was found in the different brain regions analyzed by immunohistochemistry such as brain cortex, striatum (Figures 2A–C), hippocampus (Figures 2D–F) or cerebellum (Figures 2G–I). As expected, the appearance of the immunolabeling revealed a distinct pattern compared with KGA, showing a broad nuclear staining characteristic of GLS2 in mammalian brain (insets in Figure 2), in agreement with previously reported data (Olalla et al., 2002; Cardona et al., 2015). Finally, supporting experimental evidences were obtained by immunoblot analysis of total protein extracts isolated from whole cerebral cortex and hippocampus and revealed with isoform-specific KGA antibodies (Figure 3). In agreement with immunohistochemical data revealing a decreased KGA protein expression, Western blot results also showed a significant reduction of KGA protein in both regions (Figure 3). Meanwhile, GLS2 protein levels were not significantly altered in cerebral cortex and hippocampus as assessed by immunoblot assays (results not shown).
Figure 2. Immunohistochemal staining for glutaminase gene encoding the liver-type isoforms (GLS2) in WT, HT and KOLPA1 mice brain. Immunolabeling for GLS2 isoforms yielded a nuclear staining in all the brain regions analyzed by this technique. (A–C) Cerebral cortex and insets showing details from layer III of the motor cortex (white squares). (D–F) Hippocampus and details from CA1 stratum pyramidale in insets. (G–I) Details from cerebellum are depicted. No differences were detected among the different genotypes in any of these cerebral areas. Ctx, cortex; St, striatum; CA1, CA3, hippocampal subfields; sp, stratum pyramidale; DG, dentate gyrus; Th, thalamus; G, cerebelar granular layer; P, Purkinje layer; M, cerebelar molecular layer. Scale bar: A–C: 500 μm, D–F: 250 μm, G–I: 100 μm; Insets: 20 μm.
Figure 3. Immunoblot analysis of KGA protein expression in cerebral cortex and hippocampus. Total protein extracts (40 μg) were isolated from cerebral cortex and hippocampus of WT and KOLPA1 mice and analyzed by Western blotting using isoform-specific KGA antibodies. A representative Western blot is shown in the top panel, using β-actin as a loading control. Densitometry analysis of the blots was done using the Chemi Doc System (Bio-Rad). KGA expression values normalized to β-actin are shown in the bottom panel. Values are means ± SEM of nine animals of each phenotype (WT and KOLPA1). ***p < 0.001; **p < 0.01.
Enzymatic Activity and GLS mRNA Levels in KOLPA1 Mice
The results of real-time PCR indicate that mRNA levels of GLS (KGA + GAC) remained essentially unchanged in KOLPA1 mice, with the exception of prefrontal and motor cortex where a significant decrease was observed (Figure 4; see also Castilla-Ortega et al., 2016). With regard to GA activity, we assayed enzyme catalytic activity in total protein extracts isolated from individual brain areas. The GA activity measured in whole regions from WT and KOLPA1 mice did not show significant variations, except for the motor cortex and PFC of LPA1-deficient mice which showed a significant bias to lower activity levels (Figure 5). Altogether, the results suggest a lack of correlation between GLS mRNA levels and protein/activity values in most cerebral regions studied, except for two cerebral cortex subregions: motor and PFC, where diminished mRNA levels correlated with reduced KGA protein and GA activity values.
Figure 4. Quantitative measurement of GLS mRNA levels by real-time RT-PCR. Total RNA was isolated from the different brain regions and absolute quantifications of GLS mRNA levels were performed as described in the Methods section. Five animals of each genotype were analyzed. The GLS transcripts are expressed as numbers of molecules of mRNA per ng of total RNA ± SEM (five independent experiments done in triplicate). *p < 0.05. PFC, prefrontal cortex; MC, motor cortex; STR, striatum; HIP, hippocampus; DG, dentate gyrus; CB, cerebellum.
Figure 5. Determination of total glutaminase (GA) activity in mouse brain regions. Significant differences in GA specific activity were found in motor cortex (WT n = 9; KO n = 5) and PFC (WT n = 12; KO n = 8), but not in hippocampus (WT n = 16; KO n = 12) and striatum (WT n = 15; KO n = 8). ***p < 0.001; *p < 0.05.
Anomalous Dendritic Spine Morphology in KOLPA1 in Mice
The potential involvement of LPA signaling in developmental process and neurotransmitter handling has been previously reported. Glutamate has been also involved in trophic functions and neuronal development. Therefore, we wanted to check whether the absence of LPA1 receptor and/or the strong deficit observed in GLS protein expression would affect the dendritic spines. Since hippocampus exhibited a pronounced KGA deficit and the behavioral phenotype of the LPA1 deficient mice are mainly attributed to hippocampal malfunction (Santin et al., 2009; Castilla-Ortega et al., 2010; Blanco et al., 2012a), we decided to assess dendritic spine density and morphology using the Golgi-Cox staining in the CA1 region. Representative images of basal pyramidal dendrites in CA1 are shown in Figure 6A (WT) and Figure 6B (KOLPA1). First, we quantified the spine density (number of spines per 10-μm long dendritic segments) in nonprimary dendrites. Although we did not find statistically significant differences (not shown) neither in spine density nor dendrite diameter (Mann-Whitney, p = 0.12 and p = 0.86, respectively), between KOLPA1 and WT mice, microscopic analysis revealed morphological alterations in the KO dendritic spines. Subsequently, we performed a detailed morphometric analysis of CA1 pyramidal dendrites, and first observed that dendritic spines from this region were longer, tortuous and thinner in KO (Figures 6F–H) than in WT mice (Figures 6C–E), exhibiting the typical characteristics of immature states (Ebrahimi and Okabe, 2014). Overall spine length values evidenced significant differences between genotypes in this parameter: a clear trend was found towards longer spines in the KOLPA1 mice (Figure 6I). Then, comparing lengths between the subtypes of spines (Table 1 below), we also found significant differences in thin spines between genotypes. Moreover, quantitative analysis of spine types (mushroom, thin, stubby or filopodia) indicated drastic changes in spine populations: filopodia was significantly more frequent in KO mice in detriment of the mushroom and stubby categories, which showed marked decreases in the null LPA1 mice (Figure 6J). In addition, the spine head-width of mushroom subtype was significant smaller in KOLPA1 spines (Figure 6L) than in WT (Figure 6K; t-test, ****p < 0.0001; Figure 6M).
Figure 6. Golgi-Cox impregnated basal tufts of CA1 pyramidal neurons from WT (A) and KOLPA1 (B) mice. Segments from basal dendrites (stratum oriens, white rectangles in (A,B) highlight representative dendritic segments) demonstrating the morphological differences between KOLPA1 (F–H) and WT (C–E) spines. Black arrows point to dendritic filopodia. (I) Overall spine length values showed significant differences between genotypes (WT 1.4 ± 0.05 μm, KOLPA1 2.12 ± 0.09 μm; Mann-Whitney, ****p < 0.0001). (J) Significant differences were found in the distribution of dendritic spine morphologies within stratum oriens from KOLPA1 model compared to the same region of WT mice. Filopodia type was significantly more frequent in KOLPA1 mice in detriment of both mushroom and stubby categories (t-test, Filopodia ****p < 0.0001, Stubby **p = 0.0073 and Thin, n.s.; Mann-Whitney, Mushroom ****p < 0.0001). (K,L) Spine head of mushroom subtype was smaller in KOLPA1 (L; L2,L3 depict details from L1 and L4, respectively) than in WT mice (K; K2,K3 depict details from K1 and K4, respectively). (M) Quantitative analysis demonstrated a significant decrease in the head-width of mushroom spines from KOLPA1 mice (0.52 ± 0.02 μm) compared to WT (0.91 ± 0.02 μm); t-test, ****p < 0.0001. Scale bar: A,B (10 μm); C–H (2 μm); K1,K4,L1 and L4 (2 μm); K2,K3,L2 and L3 (1 μm).
Table 1. Length (mean ± SEM; μm) of dendritic spine subtypes analyzed from wild-type (WT) and KOLPA1 mice.
Therefore, these results strongly suggest that the formation/maturation of dendritic spines could be significantly impaired in the hippocampal pyramidal neurons of KOLPA1 mice.
Silencing of LPA1 Reduces MMP-9, but Not MMP-2 and tPA Activities
To identify molecular candidates which entrain this structural change, we investigated alterations in protein expression profiles of dendritic markers known to be involved in the formation and modulation of the dendritic spine synapse. Immunoblot analysis of post-synaptic density protein 95 (PSD95) and synaptophysin revealed no significant changes in their expression levels between WT and mutant KOLPA1 mice (results not shown). Matrix metalloproteinases (MMP) have attracted attention recently as important modulators of spine morphology, along with other constitutive proteins well known to be involved in synapsis formation. Several of these proteins were analyzed in order to search for altered expression patterns which might underlie some of the changes observed in spine structure. Gelatin zymography reveals a band at approximately 82 kDa corresponding to the active form of MMP-9 (Figure 7). Quantitative image analysis revealed a significant decrease on MMP-9 activity in cortex (Figure 7A; p = 0.0001) and hippocampus (Figure 7B; p < 0.0001) of KOLPA1 mice. The overall decrease observed was around a 50% in the LPA1-deficient mice. There was no change in MMP-2 activity (cortex, p = 0.16; hippocampus, p = 0.26).
Figure 7. Zymograms of cortical (A) and hippocampal (B) gelatinase activity from WT and KOLPA1 mice. Matrix metalloproteinase (MMP)-9 and MMP-2 activities were determined in hippocampus and cerebral cortex (whole hemisphere) from both genotypes (WT n = 8; KO n = 11). Significant differences were found in MMP-9 proteolytic activity of both brain areas, showing a decreased activity in KO mice (****p < 0.0001), while MMP-2 remained unchanged.
To study the possibility that reduced MMP-9 activity was due to diminished enzyme expression, Western blot analysis was performed. Pro-active and active MMP-9 bands were detected at 92 and 82 kDa, respectively (Figure 8). Quantitative image analysis of each of the two bands revealed a reduction in active MMP-9 (Figure 8, cortex, p = 0.03; hippocampus, p = 0.0025) of KOLPA1 mice. No changes were detected in pro-active MMP-9 brain expression (Figure 8, cortex, p = 0.22; hippocampus, p = 0.07).
Figure 8. Immunoblot analysis of MMP-9 protein expression in cerebral cortex (A) and hippocampus (B) from WT and KOLPA1 mice. Quantitive analysis of proactive (92 kDa) and active (82 kDa) MMP-9 were determined by densitometry and relative to β-actin expression (WT n = 4; KO n = 5). A reduction in active MMP-9 expression were detected in KOLPA1 compared to WT mice (*p < 0.05, **p < 0.01), while no changes were detected in pro-active MMP-9 expression.
Although tPA/plasmin system has been involved in numerous aspects of brain function, including dendritic remodeling in mouse hippocampus (Pawlak et al., 2005) and visual cortex (Mataga et al., 2004), as well as inducing MMP-9 proteolytic activation since pro-active form (Wang et al., 2003), we did not detect any significant changes in tPA activity in cortex (Figure 9A, p = 0.60) and hippocampus (Figure 9B, p = 0.16) of KOLPA1 compared with WT mice.
Figure 9. Zymograms of tissue plasminogen activator (tPA) activity in the cortex (A) and hippocampus (B) of WT and KOLPA1 mice. Quantitative analysis of tPA activity was determined by densitometry of the 68 kDa band of substrate degradation (WT n = 4; KO n = 6). No changes in tPA activity were detected between both type of mice.
LPA1 Deficit Decreases KGA Protein Expression in Key Glutamatergic Regions
Using affinity-purified isoform-specific antibodies and immunohistochemistry we revealed a marked deficit of the most abundant cerebral GA: the KGA isoform encoded by the Gls gene (Blanco et al., 2012b; Martín-Rufián et al., 2012). This drastic decrease of KGA protein expression, particularly relevant in cerebral cortex and hippocampus, was completely unexpected and cannot be explained by previous experimental data reporting alterations in the brain of this KO mouse model. Although a decrease of adult cortical neurons (due to an increased cortical cell death) was previously described in KOLPA1 mice (Estivill-Torrús et al., 2008), the reported lower number of cortical neurons only occurs in layers II/III, V and VI and seems insufficient to completely justify the strong reduction found for KGA protein in cerebral cortex. Similarly, the reduced neurogenesis in the mouse DG previously noted for KOLPA1 mice (Matas-Rico et al., 2008) may contribute to part of the deficit shown by KGA in this subregion, but it does not explain the dramatic silencing of KGA expression in other hippocampal regions. Despite the overall fall in KGA protein levels, one striking feature observed was the intense perivascular mark, which positively correlated with LPA1 silencing. This gradual shift of KGA staining toward regions surrounding blood vessels may suggest a mechanism to control vascular permeability in LPA1-deficient mice. It is well known that LPA is a critical factor for regulating angiogenesis in vertebrates (Yukiura et al., 2011), as well as tight junction permeability in cultured endothelial cells (Schulze et al., 1997). Therefore, HT and KO animals may have defects on brain vascular development and/or blood-brain barrier permeability. Interestingly, glutamate is a vasoactive compound and glutamate receptors are expressed in perivascular glia and vascular endothelial cells; hence, the intense vascular KGA reactivity would be a functional response aimed at regulating the vascular tone under pathophysiological conditions (reduced or null LPA1 expression). Of note, a role of GA in cerebrovascular regulation was first postulated by our group after demonstrating concrete perivascular locations for brain GA isoforms (Olalla et al., 2008; Márquez et al., 2009).
Correlation between GLS mRNA abundance, KGA protein levels and GA activity was only apparent in two brain regions: prefrontal and motor cortex, which showed significant decreases at the three levels examined. Mice lacking LPA1 exhibit motor alterations; thus, reductions in locomotor activity, hypolocomotion during first exposure to the open field and reduced exploratory reactivity to a novel environment have been previously reported (Harrison et al., 2003; Santin et al., 2009). In contrast, the GLS mRNA levels in other brain areas, like hippocampus, remained essentially unaltered; therefore, transcript abundance did not parallel the behavior shown by KGA protein levels. This dissociation between transcript and protein levels for cerebral GA isoforms was previously found in studies dealing with their involvement in the addictive behavior of cocaine abuse (Blanco et al., 2012b) or in neuronal differentiation (Agostini et al., 2016). Hence, the mechanism(s) underlying the huge decline of KGA protein levels in brain hippocampus and cerebral cortex regions of KOLPA1 mice, other than motor and PFC, should be mostly posttranscriptional.
The maintenance of GA activity in key glutamatergic regions of KOLPA1 mice, despite substantial decrements in KGA protein, immediately suggests the existence of compensatory mechanisms. Accordingly, the glutamate content in PFC and hippocampus of KOLPA1 mice did not change as compared with WT mice (Harrison et al., 2003), reinforcing the view that other GA isozymes may be activated to keep constant glutamate levels. Although the contribution of each GA isoform to the pool of neurotransmitter Glu is currently unknown, KGA is the most abundant GA in mammalian brain (Olalla et al., 2002; Martín-Rufián et al., 2012) and the principal source of neurotransmitter glutamate (Márquez et al., 2016 and references therein). Thereby, a plausible explanation is that neurons attempt to compensate for the loss of KGA protein by increasing the process of glutamate generation through other GA isoforms expressed in brain (e.g., enhanced local activity of GLS2 isozymes). Indeed, two distinct LPA-related mechanisms are known to enhance the activity of GLS isoforms and might be operative here via other LPARs distinct from LPA1. The alternative spliced isoform of KGA, named GAC, is activated by Rho GTPases (Wang et al., 2010); interestingly, LPA stimulates activation of Rho in several neuronal cell lines, cortical progenitors and cortical neurons (Ye et al., 2002). Also, the Raf-Mek-Erk signaling pathway has been shown to increase KGA activity by phosphorylation after stimulation by several mitogens, including LPA (Thangavelu et al., 2012); therefore, post-translational modifications of the remaining KGA protein pool would cooperate in keeping basal GA activity in mutant mice. On the other hand, the existence of compensatory processes affecting glutamate signaling enzymes and receptors in KO mouse models has been previously reported. For example, null mice for the dopamine D4 receptor displayed enhanced expression of dopamine D1 and glutamate NMDA receptors (Gan et al., 2004), while mice lacking the mGlu1 receptor showed a compensatory overexpression of functional mGlu5 (Rossi et al., 2013). Finally, it is worth mentioning that null mutants mice lacking Gls gene –and therefore without KGA protein- died in the first postnatal day (Masson et al., 2006). Nevertheless, the synaptic terminals of Gls−/− neurons showed apparently normal glutamate content on individual presynaptic terminals; this unexpected result was explained by compensatory mechanisms of glutamate synthesis (Masson et al., 2006).
Role of Glutaminases in Neuronal Development and Differentiation
Axonal growth, neurite outgrowth and dendrite remodeling have been described in studies dealing with the role of LPA in developmental processes of the nervous system (Fujiwara et al., 2003). Although the post-transcriptional mechanisms underlying down-regulation of brain KGA after LPA1 silencing are currently unknown, the shortage in KGA protein may also give rise to developmental defects in cortical and hippocampal neurons which, in turn, may contribute to strong deficits in neurite outgrowth, axon growth and maturation of dendritic spines. Accordingly, in the KOGLS mouse model, which shows a slight effect on glutamatergic synaptic transmission, the authors concluded that KGA deficit should principally impact the implementation of more active neural network circuits, stressing the importance of KGA in neuronal maturation, development and differentiation (Masson et al., 2006). In fact, KGA has been associated with neuronal differentiation and proliferation of neural cells. Thus, both GLS isoforms (KGA and GAC) were upregulated during neurogenesis of human neural progenitor cells (NPCs), and their expression pattern positively correlated with the neuronal marker microtubule associated protein 2 (MAP-2; Wang et al., 2014). Most important, studies of cultured human NPCs after siRNA silencing of GLS suggest a critical role of GLS isoforms for proliferation and survival of NPCs (Wang et al., 2014). Furthermore, inhibition of GA activity in cultures of mouse embryonary cortical neurons resulted in impaired neuritogenesis (Velletri et al., 2013), while both mRNA and protein KGA levels strongly increased when cerebellar granule cells differentiate in culture (Thomas et al., 1989). This last study also demonstrated KGA upregulation in parallel with Ca2+-dependent glutamate release and formation of neurites, synaptic vesicles and synapses (Thomas et al., 1989).
In recent years, the involvement of KGA in the metabolic reprogramming underlying neuronal differentiation has been revealed. Thus, during ex vivo differentiation of primary cortical neurons, the glutamate-glutamine pathway was notably increased in parallel with upregulation of the GLS and GLS2 mRNA levels (Agostini et al., 2016). Strikingly, novel neuronal locations, distinct from mitochondria, were recently discovered for KGA through interactions with new scaffolding proteins. Thus, KGA relocalized from neuronal cell body mitochondria to neurite terminals after interaction with the protein BNIP-H (Caytaxin), which is highly expressed at the hippocampus and cerebellum; this specific KGA targeting occurs in PC12 cells when they differentiated to a neuronal phenotype, supporting a role in neuronal cell growth (Buschdorf et al., 2006). A second KGA interacting-partner was later discovered in brain: Bmcc1s, a protein also involved in the regulation of cell morphology (Boulay et al., 2013). In mouse cortical neurons, Bmcc1s interacts with the precursor KGA protein avoiding its proper targeting to the mitochondria and leading to its accumulation within the cytoplasm. Thus, KGA may traffic in neurons to distant sites and is not exclusively confined to the mitochondria; instead of it, KGA seems to be needed for neurite outgrowth, dendrites formation/maturation and establishment of novel synapses. Therefore, it is tempting to speculate that an impairment of KGA targeting to synaptic regions, as expected from the drastic downregulation of KGA shown by KOLPA1 mice, would have a negative impact in these mutant mice during brain development.
Regulation of Dendritic Structure by LPA
It is well documented the role of LPA in synapse formation through modulation of dendritic spine dynamics. For example, in cultured hippocampal neurons, the LPA1 receptor was localized in a punctate manner in dendritic spines and its overexpression correlated with increases in the density and size of spines (Pilpel and Segal, 2006). Therefore, we hypothesized that KOLPA1 mice would display changes in the morphology of their excitatory synaptic spines. Golgi staining showed hippocampal dendritic spines with a clear immature phenotype in KOLPA1 mice: the predominance of filopodia, in detriment of the mushroom and stubby spines, and the marked reduction in mushroom spine head-width are relevant structural changes with potential repercussion in synaptic neurotransmission. Of note, during the first few weeks of postnatal life, immature hippocampal dendrites show numerous filopodia while, in contrast, they are rarely observed in the mature hippocampus, although blocking synaptic transmission triggers again filopodia formation (for a review see Bourne and Harris, 2008).
On the other hand, a considerable body of evidence has recently shown a key role for MMP-9, and other extracellular proteases such as tPA, in synaptic remodeling, synaptic plasticity and cognitive processes (Mataga et al., 2004; Rivera et al., 2010; Michaluk et al., 2011). In neurons, MMP-9 is located at the postsynaptic domains of excitatory synapses and its inhibition has been linked to memory deficits in behavioral learning paradigms (Nagy et al., 2007). Of note, MMP-9 is a key regulator of dendritic spine morphology (Michaluk et al., 2011; Dziembowska and Wlodarczyk, 2012). We therefore examined the effects of LPA1 depletion on the activity of three extracellular matrix proteases in hippocampus and cortex: MMP-9, MMP-2 and tPA. Zymography revealed a marked reduction in MMP-9 activity in hippocampus and cortex of KOLPA1 mice compared with WT counterparts, without significant changes in MMP-2 nor tPA activities.
MMP-9 activation is strictly regulated during spine maturation (Tian et al., 2007). This inducible MMP is extracellularly secreted in the pro-active form in which catalytic site is blocked by a signal peptide. Several factors induce its activation from pro-active form through the cleavage of the propeptide. Western blot analysis indicates that brain active MMP-9 protein expression is reduced in KOLPA1 respect to WT mice, although there is no significant differences in the expression of brain pro-active MMP-9. These results indicate that MMP-9 activity decrement arises due to a decrease in active MMP-9 processing and release by post-translational mechanisms. In sharp contrast, no relevant changes were found for MMP-2 protein in our LPA1-deficient animal models. Remarkably, in animal models of drug addiction an activation of MMP-2 occurs after extended withdrawal from drug use: increases in MMP-2 activity, along with enhanced spine head diameter and spine number, are persistent neuroadaptations that occur after chronic drug exposure and underpin the enduring vulnerability to drug relapse (Mulholland et al., 2016).
The marked reduction of MMP-9 activity in KOLPA1 mice would prevent transient synaptic plasticity changes needed for synaptic strength, due to the lack of adequate re-shaping of spines in hippocampal neurons, thereby having important functional correlates in the neurophysiological phenotype of these animals. Interestingly, MMP-9 activity promoted maturation of dendritic spines and drive spine enlargement and synaptic potentiation (Wang et al., 2008). Moreover, spines undergo actin-dependent increases in spine-head volume during hippocampal long-term potentiation (LTP), whilst the levels of both the precursor and the active form of MMP-9 are notably elevated in the CA1 region, suggesting that upregulation of MMP-9 is related mechanistically to LTP (Huntley, 2012 and references therein). In addition, MMP-9 needs to be expressed at very high levels in the postnatal hippocampus and regulation of its expression control neuron numbers (Murase and McKay, 2012). Thus, the downregulation of MMP-9 activity shown by LPA1 null mice may play an important role in their cognitive performance: the lack of synaptic plasticity induced by MMP-9-dependent hippocampal spine remodeling would partially explain the cognitive and memory deficits reported in these animals. The reduced MMP-9 activity would be a downstream effect of the absence of LPA-induced Rho/Rho kinase (ROCK) signaling pathway, which is known to produce proteolytic enzymes like MMP-9 (Yu et al., 2014). Nevertheless, other authors found that spine morphology is affected by LPA1 independently of its three main signaling cascades; therefore, they suggest this effect was mostly mediated by interaction with other synaptic proteins (LPA1 has a PDZ-binding module at its C-terminal region; Pilpel and Segal, 2006).
Our results suggest that LPA1 exerts an important impact on KGA expression mainly through a posttranscriptional mechanism, stressing a key relationship between neuroactive lipids as LPA and glutamatergic transmission. Furthermore, the absence of LPA1 signaling downregulates expression of active MMP-9 and provokes drastic changes in hippocampal dendritic spines toward an immature phenotype. Taken together, these changes may provide a molecular basis for the role of LPA in influencing synaptic plasticity associated to cognitive and memory processes. Further studies are needed to address the concrete mechanism linking LPA1 signaling and GA expression in key glutamatergic regions and its relevance for cognitive function.
AG, MIC, FRF and JM: conceived and designed the experiments. JAC-S, AP, RS-V and CC: immunohistochemistry. JAC-S, LC and AP: western blots and activity measurements. MM-R and FJA: real-time RT-PCR. GE-T, EB and FRF: maintenance, selection and supply of KO mice and brain tissues. AP and EB: quantification and statistical analysis. AP, JAC-S and RS-V: golgi staining. MP-H and MIC: metalloproteinases and tPA assays. AG, FRF and JM: analyzed the data. JM: wrote the article.
This work was financially supported by Red de Trastornos Adictivos, Instituto de Salud Carlos III (ISC-III) and European Regional Development Funds-European Union (ERDF-EU) Grants RD16/0017/001 (JM and FRF), RD12/0028/0013 (JM), RD12/0028/0001 (FRF), and RD12/0028/0002 (MIC). Additional funds come from Ministerio de Economía y Competitividad, Proyectos de investigación en Salud (PI16/01698 to FRF, PI16/01953) Grant FIS PI12/01431 and FIS PI15/00796 (to AG), Grant SAF2015-64501-R (to JM) and SAF2013-40592-R (to MIC) and Excellence Grant-Regional Andalusian government CVI-6656 (to JM) and CTS-8221 (to FRF).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
DAB, diaminobenzidine tetrahydrochloride; GA, glutaminase; GAB, Gls2-encoded long glutaminase protein variant; GAC, Gls-encoded short glutaminase protein variant; Gls, glutaminase gene encoding the kidney-type isoforms; Gls2, glutaminase gene encoding the liver-type isoforms; KGA, Gls-encoded long glutaminase protein variant; LGA, Gls2-encoded short glutaminase protein variant; LPA, lysophosphatidic acid; LPARs, lysophosphatidic acid receptors; LTP, long-term potentiation; MMP, matrix metalloproteinase; NMDA, N-methyl-D-aspartate; NPCs, neural progenitor cells; NP-40, Nonidet-P40; PBS, phosphate-buffered saline; PFC, prefrontal cortex; PSD95, post-synaptic density protein 95; ROCK, Rho kinase; TBS, Tris-buffered saline; TCA, trichloroacetic acid; tPA, tissue plasminogen activator.
Agostini, M., Romeo, F., Inoue, S., Niklison-Chirou, M. V., Elia, A. J., Dinsdale, D., et al. (2016). Metabolic reprogramming during neuronal differentiation. Cell Death Differ. 23, 1502–1514. doi: 10.1038/cdd.2016.36
Aledo, J. C., Gómez-Fabre, P. M., Olalla, L., and Márquez, J. (2000). Identification of two human glutaminase loci and tissue-specific expression of the two related genes. Mamm. Genome 11, 1107–1110. doi: 10.1007/s003350010190
Blanco, E., Bilbao, A., Luque-Rojas, M. J., Palomino, A., Bermúdez-Silva, F. J., Suárez, J., et al. (2012a). Attenuation of cocaine-induced conditioned locomotion is associated with altered expression of hippocampal glutamate receptors in mice lacking LPA1 receptors. Psychopharmacology 220, 27–42. doi: 10.1007/s00213-011-2446-6
Blanco, E., Campos-Sandoval, J. A., Palomino, A., Luque-Rojas, M. J., Bilbao, A., Suárez, J., et al. (2012b). Cocaine modulates both glutaminase gene expression and glutaminase activity in the brain of cocaine-sensitized mice. Psychopharmacology 219, 933–944. doi: 10.1007/s00213-011-2418-x
Boulay, A.-C., Burbassi, S., Lorenzo, H.-K., Loew, D., Ezan, P., Giaume, C., et al. (2013). Bmcc1s interacts with the phosphate-activated glutaminase in the brain. Biochimie 95, 799–807. doi: 10.1016/j.biochi.2012.11.016
Bourne, J. N., and Harris, K. M. (2008). Balancing structure and function at hippocampal dendritic spines. Annu. Rev. Neurosci. 31, 47–67. doi: 10.1146/annurev.neuro.31.060407.125646
Buschdorf, J. P., Chew, L. L., Zhang, B., Cao, Q., Liang, F. Y., Liou, Y. C., et al. (2006). Brain-specific BNIP-2-homology protein Caytaxin relocalises glutaminase to neurite terminals and reduces glutamate levels. J. Cell Sci. 119, 3337–3350. doi: 10.1242/jcs.03061
Campos, J. A., Aledo, J. C., Segura, J. A., Alonso, F. J., Gómez-Fabre, P. M., Núñez de Castro, I., et al. (2003). Expression of recombinant human L-glutaminase in Escherichia coli: polyclonal antibodies production and immunological analysis of mouse tissues. Biochim. Biophys. Acta 1648, 17–23. doi: 10.1016/s1570-9639(03)00026-8
Cardona, C., Sánchez-Mejías, E., Dávila, J. C., Martín-Rufián, M., Campos-Sandoval, J. A., Vitorica, J., et al. (2015). Expression of Gls and Gls2 glutaminase isoforms in astrocytes. Glia 63, 365–382. doi: 10.1002/glia.22758
Castilla-Ortega, E., Pavón, F. J., Sánchez-Marín, L., Estivill-Torrús, G., Pedraza, C., Blanco, E., et al. (2016). Both genetic deletion and pharmacological blockade of lysophosphatidic acid LPA1 receptor results in increased alcohol consumption. Neuropharmacology 103, 92–103. doi: 10.1016/j.neuropharm.2015.12.010
Castilla-Ortega, E., Sánchez-López, J., Hoyo-Becerra, C., Matas-Rico, E., Zambrana-Infantes, E., Chun, J., et al. (2010). Exploratory, anxiety and spatial memory impairments are dissociated in mice lacking the LPA1 receptor. Neurobiol. Learn. Mem. 94, 73–82. doi: 10.1016/j.nlm.2010.04.003
Choi, J. W., and Chun, J. (2013). Lysophospholipids and their receptors in the central nervous system. Biochim. Biophys. Acta 1831, 20–32. doi: 10.1016/j.bbalip.2012.07.015
Choi, J. W., Herr, D. R., Noguchi, K., Yung, Y. C., Lee, C. W., Mutoh, T., et al. (2010). LPA receptors: subtypes and biological actions. Annu. Rev. Pharmacol. Toxicol. 50, 157–186. doi: 10.1146/annurev.pharmtox.010909.105753
Chun, J., Hla, T., Lynch, K. R., Spiegel, S., and Moolenaar, W. H. (2010). International union of basic and clinical pharmacology. LXXVIII. Lysophospholipid receptor nomenclature. Pharmacol. Rev. 62, 579–587. doi: 10.1124/pr.110.003111
Contos, J. J., Fukushima, N., Weiner, J. A., Kaushal, D., and Chun, J. (2000). Requirement for the lpA1 lysophosphatidic acid receptor gene in normal suckling behavior. Proc. Natl. Acad. Sci. U S A 97, 13384–13389. doi: 10.1073/pnas.97.24.13384
Dziembowska, M., and Wlodarczyk, J. (2012). MMP9: a novel function in synaptic plasticity. Int. J. Biochem. Cell Biol. 44, 709–713. doi: 10.1016/j.biocel.2012.01.023
Ebrahimi, S., and Okabe, S. (2014). Structural dynamics of dendritic spines: molecular composition, geometry and functional regulation. Biochim. Biophys. Acta 1838, 2391–2398. doi: 10.1016/j.bbamem.2014.06.002
Estivill-Torrús, G., Llebrez-Zayas, P., Matas-Rico, E., Santín, L., Pedraza, C., De Diego, I., et al. (2008). Absence of LPA1 signaling results in defective cortical development. Cereb. Cortex 18, 938–950. doi: 10.1093/cercor/bhm132
Fujiwara, Y., Sebök, A., Meakin, S., Kobayashi, T., Murakami-Murofushi, K., and Tigyi, G. (2003). Cyclic phosphatidic acid elicits neurotrophin-like actions in embryonic hippocampal neurons. J. Neurochem. 87, 1272–1283. doi: 10.1046/j.1471-4159.2003.02106.x
Fukushima, N., Furuta, D., and Tsujiuchi, T. (2011). Coordinated interactions between actin and microtubules through crosslinkers in neurite retraction induced by lysophosphatidic acid. Neurochem. Int. 59, 109–113. doi: 10.1016/j.neuint.2011.04.020
Gan, L., Falzone, T. L., Zhang, K., Rubinstein, M., Baldessarini, R. J., and Tarazi, F. I. (2004). Enhanced expression of dopamine D1 and glutamate NMDA receptors in dopamine D4 receptor knockout mice. J. Mol. Neurosci. 22, 167–178. doi: 10.1385/jmn:22:3:167
García-Morales, V., Montero, F., González-Forero, D., Rodríguez-Bey, G., Gómez-Pérez, L., Medialdea-Wandossell, M. J., et al. (2015). Membrane-derived phospholipids control synaptic neurotransmission and plasticity. PLoS Biol. 13:e1002153. doi: 10.1371/journal.pbio.1002153
Ghani, M. U., Mesadi, F., Kanik, S. D., Argunsah, A. Ö., Hobbiss, A. F., Israely, I., et al. (2017). Dendritic spine classification using shape and appearance features based on two-photon microscopy. J. Neurosci. Methods 279, 13–21. doi: 10.1016/j.jneumeth.2016.12.006
Harrison, S. M., Reavill, C., Brown, G., Brown, J. T., Cluderay, J. E., Crook, B., et al. (2003). LPA1 receptor-deficient mice have phenotypic changes observed in psychiatric disease. Mol. Cell. Neurosci. 24, 1170–1179. doi: 10.1016/j.men.2003.09.001
Hecht, J. H., Weiner, J. A., Post, S. R., and Chun, J. (1996). Ventricular zone gene-1 (vzg-1) encodes a lysophosphatidic acid receptor expressed in neurogenic regions of the developing cerebral cortex. J. Cell Biol. 135, 1071–1083. doi: 10.1083/jcb.135.4.1071
Heini, H. G., Gebhardt, R., Brecht, A., and Mecke, D. (1987). Purification and characterization of rat liver glutaminase. Eur. J. Biochem. 162, 541–546. doi: 10.1111/j.1432-1033.1987.tb10673.x
Holtsberg, F. W., Steiner, M. R., Furukawa, K., Keller, J. N., Mattson, M. P., and Steiner, S. M. (1997). Lysophosphatidic acid induces a sustained elevation of neuronal intracellular calcium. J. Neurochem. 69, 68–75. doi: 10.1046/j.1471-4159.1997.69010068.x
Huntley, G. W. (2012). Synaptic circuit remodeling by matrix metalloproteinases in health and disease. Nat. Rev. Neurosci. 13, 743–757. doi: 10.1038/nrn3320
Lu, W. Y., Xiong, Z. G., Lei, S., Orser, B. A., Dudek, E., Browning, M. D., et al. (1999). G-protein-coupled receptors act via protein kinase C and Src to regulate NMDA receptors. Nat. Neurosci. 2, 331–338. doi: 10.1038/7243
Manning, T. J. Jr., Parker, J. C., and Sontheimer, H. (2000). Role of lysophosphatidic acid and rho in glioma cell motility. Cell Motil. Cytoskeleton 45, 185–199. doi: 10.1002/(SICI)1097-0169(200003)45:3<185::AID-CM2>3.0.CO;2-G
Márquez, J., Matés, J. M., and Campos-Sandoval, J. A. (2016). “Glutaminases,” in Advances in Neurobiology, 13. The Glutamate/GABA/Glutamine Cycle: Amino Acid Neurotransmitter Homeostasis, eds U. Sonnewald and A. Schousboe (Vienna: Springer Verlag), 133–172.
Márquez, J., Tosina, M., de la Rosa, V., Segura, J. A., Alonso, F. J., Matés, J. M., et al. (2009). New insights into brain glutaminases: beyond their role on glutamatergic transmission. Neurochem. Int. 55, 64–70. doi: 10.1016/j.neuint.2009.02.022
Martín-Rufián, M., Tosina, M., Campos-Sandoval, J. A., Manzanares, E., Lobo, C., Segura, J. A., et al. (2012). Mammalian glutaminase Gls2 gene encodes two functional alternative transcripts by a surrogate promoter usage mechanism. PLoS One 7:e38380. doi: 10.1371/journal.pone.0038380
Masson, J., Darmon, M., Conjard, A., Chuhma, N., Ropert, N., Thoby-Brisson, M., et al. (2006). Mice lacking brain/kidney phosphate-activated glutaminase have impaired glutamatergic synaptic transmission, altered breathing, disorganized goal-directed behavior and die shortly after birth. J. Neurosci. 26, 4660–4671. doi: 10.1523/JNEUROSCI.4241-05.2006
Mataga, N., Mizuguchi, Y., and Hensch, T. K. (2004). Experience-dependent pruning of dendritic spines in visual cortex by tissue plasminogen activator. Neuron 44, 1031–1041. doi: 10.1016/j.neuron.2004.11.028
Matas-Rico, E., García-Diaz, B., Llebrez-Zayas, P., López-Barroso, D., Santín, L., Pedraza, C., et al. (2008). Deletion of lysophosphatidic acid receptor LPA1 reduces neurogenesis in the mouse dentate gyrus. Mol. Cell. Neurosci. 39, 342–355. doi: 10.1016/j.mcn.2008.07.014
Matés, J. M., Segura, J. A., Martín-Rufián, M., Campos-Sandoval, J. A., Alonso, F. J., and Márquez, J. (2013). Glutaminase isoenzymes as key regulators in metabolic and oxidative stress against cancer. Curr. Mol. Med. 13, 514–534. doi: 10.2174/1566524011313040005
Michaluk, P., Wawrzyniak, M., Alot, P., Szczot, M., Wyrembek, P., Mercik, K., et al. (2011). Influence of matrix metalloproteinase MMP-9 on dendritic spine morphology. J. Cell Sci. 124, 3369–3380. doi: 10.1242/jcs.090852
Mulholland, P. J., Chandler, L. J., and Kalivas, P. W. (2016). Signals from the fourth dimension regulate drug relapse. Trends Neurosci. 39, 472–485. doi: 10.1016/j.tins.2016.04.007
Murase, S., and McKay, R. D. (2012). Matrix metalloproteinase-9 regulates survival of neurons in newborn hippocampus. J. Biol. Chem. 287, 12184–12194. doi: 10.1074/jbc.M111.297671
Musazzi, L., Di Daniel, E., Maycox, P., Racagni, G., and Popoli, M. (2011). Abnormalities in α/β-CaMKII and related mechanisms suggest synaptic dysfunction in hippocampus of LPA1 receptor knockout mice. Int. J. Neuropsychopharmacol. 14, 941–953. doi: 10.1017/S1461145710001240
Mutoh, T., Rivera, R., and Chun, J. (2012). Insights into the pharmacological relevance of lysophospholipid receptors. Br. J. Pharmacol. 165, 829–844. doi: 10.1111/j.1476-5381.2011.01622.x
Nagy, V., Bozdagi, O., and Huntley, G. W. (2007). The extracellular protease matrix metalloproteinase-9 is activated by inhibitory avoidance learning and required for long-term memory. Learn. Mem. 14, 655–664. doi: 10.1101/lm.678307
Nishikawa, T., Tomori, Y., Yamashita, S., and Shimizu, S. (1989). Inhibition of Na+, K+-ATPase activity by phospholipase A2 and several lysophospholipids: possible role of phospholipase A2 in noradrenaline release from cerebral cortical synaptosomes. J. Pharm. Pharmacol. 41, 450–458. doi: 10.1111/j.2042-7158.1989.tb06499.x
Olalla, L., Gutiérrez, A., Campos, J. A., Khan, Z. U., Alonso, F., Segura, J. A., et al. (2002). Nuclear localization of L-glutaminase in mammalian brain. J. Biol. Chem. 277, 38939–38944. doi: 10.1074/jbc.c200373200
Olalla, L., Gutiérrez, A., Jiménez, A. J., López-Téllez, J. F., Khan, Z. U., Pérez, J., et al. (2008). Expression of scaffolding PDZ protein GIP (Glutaminase-Interacting-Protein) in mammalian brain. J. Neurosci. Res. 86, 281–292. doi: 10.1002/jnr.21505
Pawlak, R., Rao, B. S. S., Melchor, J. P., Chattarji, S., McEwen, B., and Strickland, S. (2005). Tissue plasminogen activator and plasminogen mediate stress-induced decline of neuronal and cognitive functions in the mouse hippocampus. Proc. Natl. Acad. Sci. U S A 102, 18201–18206. doi: 10.1073/pnas.0509232102
Paxinos, G., and Franklin, K. B. J. (2001). The Mouse Brain in Stereotaxic Coordinates. 2nd Edn. San Diego, CA: Academic Press.
Pilpel, Y., and Segal, M. (2006). The role of LPA1 in formation of synapses among cultured hippocampal neurons. J. Neurochem. 97, 1379–1392. doi: 10.1111/j.1471-4159.2006.03825.x
Rivera, S., Khrestchatisky, M., Kaczmarek, L., Rosenberg, G. A., and Jaworski, D. M. (2010). Metzincin proteases and their inhibitors: foes or friends in nervous system physiology? J. Neurosci. 30, 15337–15357. doi: 10.1523/JNEUROSCI.3467-10.2010
Roberts, C., Winter, P., Shilliam, C. S., Hughes, Z. A., Langmead, C., Maycox, P. R., et al. (2005). Neurochemical changes in LPA1 receptor deficient mice—a putative model of schizophrenia. Neurochem. Res. 30, 371–377. doi: 10.1007/s11064-005-2611-6
Rossi, P. I., Musante, I., Summa, M., Pittaluga, A., Emionite, L., Ikehata, M., et al. (2013). Compensatory molecular and functional mechanisms in nervous system of the Grm1crv4 mouse lacking the mGlu1 receptor: a model for motor coordination deficits. Cereb. Cortex 23, 2179–2189. doi: 10.1093/cercor/bhs200
Ruszczycki, B., Szepesi, Z., Wilczynski, G. M., Bijata, M., Kalita, K., Kaczmarek, L., et al. (2012). Sampling issues in quantitative analysis of dendritic spines morphology. BMC Bioinformatics 13:213. doi: 10.1186/1471-2105-13-213
Santin, L. J., Bilbao, A., Pedraza, C., Matas-Rico, E., López-Barroso, D., Castilla-Ortega, E., et al. (2009). Behavioral phenotype of maLPA1-null mice: increased anxiety-like behavior and spatial memory deficits. Genes Brain Behav. 8, 772–784. doi: 10.1111/j.1601-183x.2009.00524.x
Schmidt, A., Wolde, M., Thiele, C., Fest, W., Kratzin, H., Podtelejnikov, A. V., et al. (1999). Endophilin I mediates synaptic vesicle formation by transfer of arachidonate to lysophosphatidic acid. Nature 401, 133–141. doi: 10.1038/43613
Schulze, C., Smales, C., Rubin, L. L., and Staddon, J. M. (1997). Lysophosphatidic acid increases tight junction permeability in cultured brain endothelial cells. J. Neurochem. 68, 991–1000. doi: 10.1046/j.1471-4159.1997.68030991.x
Shano, S., Moriyama, R., Chun, J., and Fukushima, N. (2008). Lysophosphatidic acid stimulates astrocyte proliferation through LPA1. Neurochem. Int. 52, 216–220. doi: 10.1016/j.neuint.2007.07.004
Spohr, T. C., Choi, J. W., Gardell, S. E., Herr, D. R., Rehen, S. K., Gomes, F. C., et al. (2008). Lysophosphatidic acid receptor-dependent secondary effects via astrocytes promote neuronal differentiation. J. Biol. Chem. 283, 7470–7479. doi: 10.1074/jbc.M707758200
Spohr, T. C., Dezonne, R. S., Rehen, S. K., and Gomes, F. C. (2011). Astrocytes treated by lysophosphatidic acid induce axonal outgrowth of cortical progenitors through extracellular matrix protein and epidermal growth factor signaling pathway. J. Neurochem. 119, 113–123. doi: 10.1111/j.1471-4159.2011.07421.x
Steiner, M. R., Urso, J. R., Klein, J., and Steiner, S. M. (2002). Multiple astrocyte responses to lysophosphatidic acids. Biochim. Biophys. Acta 1582, 154–160. doi: 10.1016/s1388-1981(02)00150-6
Thangavelu, K., Pan, C. Q., Karlberg, T., Balaji, G., Uttamchandani, M., Suresh, V., et al. (2012). Structural basis for the allosteric inhibitory mechanism of human kidney-type glutaminase (KGA) and its regulation by Raf-Mek-Erk signaling in cancer cell metabolism. Proc. Natl. Acad. Sci. U S A 109, 7705–7710. doi: 10.1073/pnas.1116573109
Thomas, J. W., Novelli, A., Tao-Cheng, J. H., Henneberry, R., Smith, H. H., and Banner, C. (1989). Developmental induction of glutaminase in primary cultures of cerebellar granule cells. Mol. Brain Res. 6, 47–54. doi: 10.1016/0169-328x(89)90027-2
Tian, L., Stefanidakis, M., Ning, L., Van Lint, P., Nyman-Huttunen, H., Libert, C., et al. (2007). Activation of NMDA receptors promotes dendritic spine development through MMP-mediated ICAM-5 cleavage. J. Cell. Biol 178, 687–700. doi: 10.1083/jcb.200612097
Velletri, T., Romeo, F., Tucci, P., Peschiaroli, A., Annicchiarico-Petruzzelli, M., Niklison-Chirou, M. V., et al. (2013). GLS2 is transcriptionally regulated by p73 and contributes to neuronal differentiation. Cell Cycle 12, 3564–3573. doi: 10.4161/cc.26771
Wang, X. B., Bozdagi, O., Nikitczuk, J. S., Zhai, Z. W., Zhou, Q., and Huntley, G. W. (2008). Extracellular proteolysis by matrix metalloproteinase-9 drives dendritic spine enlargement and long-term potentiation coordinately. Proc. Natl. Acad. Sci. U S A 105, 19520–19525. doi: 10.1073/pnas.0807248105
Wang, J.-B., Erickson, J. W., Fuji, R., Ramachandran, S., Gao, P., Dinavahi, R., et al. (2010). Targeting mitochondrial glutaminase activity inhibits oncogenic transformation. Cancer Cell 18, 207–219. doi: 10.1016/j.ccr.2010.08.009
Wang, Y., Huang, Y., Zhao, L., Li, Y., and Zheng, J. (2014). Glutaminase 1 is essential for the differentiation, proliferation, and survival of human neural progenitor cells. Stem Cells Dev. 23, 2782–2790. doi: 10.1089/scd.2014.0022
Wang, X., Lee, S.-R., Arai, K., Lee, S.-R., Tsuji, K., Rebeck, G. W., et al. (2003). Lipoprotein receptor-mediated induction of matrix metalloproteinase by tissue plasminogen activator. Nat. Med. 9, 1313–1317. doi: 10.1038/nm926
Ye, X., Fukushima, N., Kingsbury, M. A., and Chun, J. (2002). Lysophosphatidic acid in neural signaling. Neuroreport 13, 2169–2175. doi: 10.1097/00001756-200212030-00002
Yu, Y., Qin, J., Liu, M., Ruan, Q., Li, Y., and Zhang, Z. (2014). Role of Rho kinase in lysophosphatidic acid-induced altering of blood-brain barrier permeability. Int. J. Mol. Med. 33, 661–669. doi: 10.3892/ijmm.2014.1618
Yukiura, H., Hama, K., Nakanaga, K., Tanaka, M., Asaoka, Y., Okudaira, S., et al. (2011). Autotaxin regulates vascular development via multiple lysophosphatidic acid (LPA) receptors in zebrafish. J. Biol. Chem. 286, 43972–43983. doi: 10.1074/jbc.M111.301093
Yung, Y. C., Stoddard, N. C., and Chun, J. (2014). LPA receptor signaling: pharmacology, physiology, and pathophysiology. J. Lipid Res. 55, 1192–1214. doi: 10.1194/jlr.R046458
Yung, Y. C., Stoddard, N. C., Mirendil, H., and Chun, J. (2015). Lysophosphatidic acid (LPA) signalling in the nervous system. Neuron 85, 669–682. doi: 10.1016/j.neuron.2015.01.009
Keywords: lysophosphatidic acid, glutaminase, glutamate, matrix metalloproteinases, rodent knockout model, synaptic plasticity
Citation: Peñalver A, Campos-Sandoval JA, Blanco E, Cardona C, Castilla L, Martín-Rufián M, Estivill-Torrús G, Sánchez-Varo R, Alonso FJ, Pérez-Hernández M, Colado MI, Gutiérrez A, de Fonseca FR and Márquez J (2017) Glutaminase and MMP-9 Downregulation in Cortex and Hippocampus of LPA1 Receptor Null Mice Correlate with Altered Dendritic Spine Plasticity. Front. Mol. Neurosci. 10:278. doi: 10.3389/fnmol.2017.00278
Received: 26 May 2017; Accepted: 17 August 2017;
Published: 05 September 2017.
Edited by:Isabel Varela-Nieto, Consejo Superior de Investigaciones Científicas (CSIC), Spain
Reviewed by:Magdalena Zielińska, Mossakowski Medical Research Centre (PAN), Poland
Chiara Donati, University of Florence, Italy
Copyright © 2017 Peñalver, Campos-Sandoval, Blanco, Cardona, Castilla, Martín-Rufián, Estivill-Torrús, Sánchez-Varo, Alonso, Pérez-Hernández, Colado, Gutiérrez, de Fonseca and Márquez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fernando Rodríguez de Fonseca, email@example.com
Javier Márquez, firstname.lastname@example.org
† These authors have contributed equally to this work.