Abstract
Down syndrome (DS), the main genetic cause of intellectual disability, is associated with an imbalance of excitatory/inhibitory neurotransmitter systems. The phenotypic assessment and pharmacotherapy interventions in DS murine models strongly pointed out glutamatergic neurotransmission alterations (specially affecting ionotropic glutamate receptors [iGluRs]) that might contribute to DS pathophysiology, which is in agreement with DS condition. iGluRs play a critical role in fast-mediated excitatory transmission, a process underlying synaptic plasticity. Neuronal plasticity is biochemically modulated by post-translational modifications, allowing rapid and reversible adaptation of synaptic strength. Among these modifications, phosphorylation/dephosphorylation processes strongly dictate iGluR protein–protein interactions, cell surface trafficking, and subsynaptic mobility. Hence, we hypothesized that dysregulation of phosphorylation/dephosphorylation balance might affect neuronal function, which in turn could contribute to the glutamatergic neurotransmitter alterations observed in DS. To address this point, we biochemically purified subsynaptic hippocampal fractions from adult Ts65Dn mice, a trisomic mouse model recapitulating DS phenotypic alterations. Proteomic analysis showed significant alterations of the molecular composition of subsynaptic compartments of hippocampal trisomic neurons. Further, we characterized iGluR phosphopattern in the hippocampal glutamatergic synapse of trisomic mice. Phosphoenrichment-coupled mass spectrometry analysis revealed specific subsynaptic- and trisomy-associated iGluR phosphorylation signature, concomitant with differential subsynaptic kinase and phosphatase composition of Ts65Dn hippocampal subsynaptic compartments. Furthermore, biochemical data were used to build up a genotype-kinome-iGluR phosphopattern matrix in the different subsynaptic compartments. Overall, our results provide a precise profile of iGluR phosphopattern alterations in the glutamatergic synapse of the Ts65Dn mouse model and support their contribution to DS-associated synaptopathy. The alteration of iGluR phosphoresidues in Ts65Dn hippocampi, together with the kinase/phosphatase signature, identifies potential novel therapeutic targets for the treatment of glutamatergic dysfunctions in DS.
Introduction
Down syndrome (DS) is the most common genetic cause of intellectual disability (OMIM #190685). Individuals with DS present a trisomy of the human chromosome 21 (Hsa21), and pivotal genes located in Hsa21 have been proposed as candidate genes for DS pathophysiology (Antonarakis et al., ). To model DS pathophysiology, different animal models have been used, either overexpressing single Hsa21-homologous genes or encompassing a partial trisomy of Hsa21 syntenic region (Dierssen, ; Rueda et al., 2012). Among the latter, Ts65Dn mouse is the most well-established DS murine model and genetically consists of the triplication of the Hsa21 syntenic mouse chromosome Mmu16. Ts65Dn mouse recapitulates cognitive and motor alterations observed in individuals with DS (Costa et al., ; Galdzicki and Siarey, ; Altafaj et al., ). These alterations are neurochemically associated with an imbalance of excitatory/inhibitory neurotransmitter systems (Belichenko et al., ), similar to those described in individuals with DS (Risser et al., 1997; Seidl et al., 2001). This neurotransmitter imbalance resulting from the alteration of GABAergic inhibitory neurotransmission (Fernandez et al., ; Best et al., ; Contestabile et al., ), together with excitatory glutamatergic neurotransmission defects (Siddiqui et al., 2008; Kleschevnikov, ), might disturb synaptic plasticity and impair learning and memory abilities in Ts65Dn mice. Consequently, several attempts have been made toward the restoration of homeostatic excitatory/inhibitory balance. In particular, Ts65Dn treatment with memantine—a reversible N-methyl-d-aspartate glutamate receptor (NMDAR) antagonist—has been shown to partially rescue critical electrophysiological alterations (Scott-McKean and Costa, 2011) and cognitive phenotypic abnormalities of trisomic mice (Costa et al., ; Lockrow et al., ). These findings have been significantly translated to the clinical practice, with an improvement of the cognitive skills of individuals with DS (Boada et al., ). These data strongly support the hypothesis that Hsa21 gene overexpression (e.g., DYRK1A, APP, TIAM1, and ITSN1 gene products) might dysregulate NMDAR expression and/or function, contributing to DS-like pathogenesis in trisomic mice (Siddiqui et al., 2008). NMDARs belong to the ionotropic glutamate receptor (iGluR) family, together with α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid receptors (AMPARs) and kainate receptors (KARs). iGluRs importantly play a critical role in fast-mediated transmission, a process underlying synaptic plasticity and neuronal survival (Paoletti et al., 2013). The precise control of iGluR activity is tuned by numerous post-translational mechanisms (Wang et al., 2006; Traynelis et al., 2010). Among them, (de)phosphorylation events strongly influence iGluR biogenesis (assembly, trafficking, docking, surface diffusion, and internalization) and receptor channel gating. These processes modulate synaptic strength and thus modulate neuronal plasticity (Lee, ; Lau and Zukin, ; Goebel-Goody et al., ; Pahl et al., ). The phosphorylation pattern of iGluR results from the activity of several canonical protein kinases (PKA, PKC, CaMKII) (reviewed by Lussier et al., ) and protein phosphatases (STEP, PP2A, PP1, calcineurin) (Westphal, 1999; Chan and Sucher, ; Braithwaite et al., ; Sanderson et al., 2016) that collectively modulate the gating and subcellular location of iGluR. It has been proposed that the specific subsynaptic location of NMDARs (synaptic vs. extrasynaptic) plays a dichotomic role related with synaptic plasticity and neuronal survival (Hardingham and Bading, ). We have recently showed that the gene product of DS candidate gene DYRK1A (dual-specificity tyrosine phosphorylation-regulated kinase 1) also influences NMDAR surface expression and function similar to NMDAR regulation by Hsa21 gene products (Grau et al., ). To sum up, these data suggest that the overexpression of Hsa21 genes impacts glutamatergic transmission in DS. This potential functional disturbance might molecularly be associated with and/or result from a complex reorganization of the synaptic proteome composition and protein phosphopattern, with protein-specific dysregulation (subsets of proteins under-represented, over-represented, with conserved expression and with protein and site-specific phosphorylation changes). To identify these potential molecular changes, in this study, we have optimized a biochemical subfractioning method allowing the purification of subsynaptic compartments (pre-, extra-, and postsynaptic fractions) for further proteomic and phosphoproteomic analysis of Ts65Dn adult mice hippocampi. Phosphoenrichment-coupled mass spectrometry analysis revealed a subsynaptic-specific and disease-associated iGluR phosphopattern signature. This phosphopattern is in agreement with the concomitant subsynaptic-specific increase of kinase expression and phosphatase downregulation in Ts65Dn mice hippocampi. Overall, our results demonstrate an altered phosphopattern in the glutamatergic synapse, together with the identification of a protein kinase/phosphatase biochemical signature in the Ts65Dn murine model, which may represent novel therapeutic targets for DS synaptopathy.
Methods
Mice
Ts65Dn mouse colony female B6EiC3Sn a/A-Ts(1716)65Dn (Ts65Dn) and male B6C3F1/J mice were purchased from the Jackson Laboratory (Bar Harbor, ME). The mouse colony was housed and bred in the Animal Facilities of the Barcelona Biomedical Research Park (PRBB, Barcelona, Spain, EU). All animal procedures met the guidelines of the European Community Directive 86/609/EEC and were approved by the Local Ethics Committee. The mice were housed under a 12:12 h light–dark schedule (lights on at 8:00 a.m.) in controlled environmental conditions of humidity (60%) and temperature (22 ± 2°C) with food and water ad libitum. Both Ts65Dn and euploid mice were genotyped by qPCR, following the Jackson laboratories protocol (https://www.jax.org/research-and-faculty/tools/cytogenetic-and-down-syndrome-models-resource/protocols/cytogenic-qpcr-protocol).
Antibodies
Antibodies were obtained from commercial suppliers as it follows: anti-PSD95 (Neuromab, #P78352), anti-synaptophysin (Sigma-Aldrich, S5768), anti-β-actin (Sigma-Aldrich, A5316), anti-GluN2A (Sigma-Aldrich, M264), anti-GluN1 (Millipore, #AB9864R), anti-CaMKIIα (CST, #50049), anti-Fyn (CST, #4023), anti-pS890 GluN1 (CST, #3381), anti-pS1284 GluN2B (CST, #5355), anti-rabbit IgG-horseradish peroxidase (HRP; Dako, #P0448), and anti-mouse IgG-HRP (Dako, # P0447).
Subsynaptic fractionation and western blot
Subsynaptic fractionation protocol was adapted from the protocol developed by Philips and collaborators (Phillips et al., 2001). This protocol allows the obtention of the enriched extrasynaptic (Extra), postsynaptic (Post), and presynaptic (Pre) fractions from brain regions. In brief, the frozen tissue was thawed in cold buffer A (50 mM Tris-HCl pH 7.4, 0.32 M sucrose, 5 mM EDTA, 1 mM EGTA, 1 μg/ml aprotinin, 1 μg/ml leupeptin, 1/2500 PMSF, 20 μM ZnCl2, 50 mM NaF, 1 mM sodium orthovanadate, 2.5 mM sodium pyrophosphate, and a cocktail of protease inhibitors [Complete, Roche]) and mechanically homogenized using a potter. The homogenate was centrifuged for 10 min at 1,400 g, and the supernatant was collected. To increase the yield of protein recovery, this step was repeated thrice. The collected supernatants were pooled and centrifuged for 10 min at 700 g. The resulting supernatant was again centrifuged for 15 min at 21,000 g. The resultant pellet containing the crude membrane fraction was kept for further subfractionation, as follows. Crude membrane fraction was solubilized in buffer B (0.32 M sucrose, 50 mM Tris-HCl, pH 7.4) and loaded on a discontinuous sucrose step gradient (0.85 M/1.0 M/1.2 M). Further after centrifugation at 82,500 g for 2 h, the synaptosomes (Syns) were collected from the 1.0 M/1.2 M interface and diluted in 50 mM Tris-HCl pH 7.4 and centrifuged once more for additional 30 min at 21,000 g. Then this pellet was collected and resuspended in synaptosomal resuspension buffer, consisting of 36 mM sucrose 0.1 mM CaCl2 + 10 mM NaF + 1 mM sodium orthovanadate, 2 mM EDTA, protease, and phosphatase inhibitor cocktail (HaltTM, ThermoFisher), leading to the purified synaptosomal fraction (P2) for further subsynaptic fractionation, as follows. In brief, synaptosomes were solubilized by resuspending in 40 mM Tris-HCl pH 6.0 containing 1% triton X-100 and incubated at 4°C for 30 min under agitation. Following centrifugation at 40,000 g for 30 min, the supernatant fraction (corresponding to the extrasynaptic fraction) was acetone precipitated. The pellet, containing the synaptosomal junction, was resuspended in 20 mM Tris-HCl buffer pH 8.0 supplemented with 1% triton X-100. After 30-min incubation at 4°C, the fractions were separated by centrifugation at 40,000 g for 30 min. The resulting supernatant, corresponding to the enriched presynaptic fraction, was acetone precipitated. The pellet, containing the enriched postsynaptic fraction, was resuspended in 50 mM Tris buffer—2% SDS, as well as the resulting acetone-precipitated subsynaptic fractions.
For western blot analysis, equal protein amounts from each sample were loaded and separated using 8 and 10% SDS-PAGE. The proteins were transferred onto nitrocellulose membranes (Invitrogen) using iBlot semi-dry blotting system (Life Technologies). After transient staining of transferred proteins by Ponceau S staining, the membranes were blocked with 10% skim milk in TBST (10 mM Tris-HCl, 100 mM NaCl, 0.1% Tween-20) for 1 h. The membranes were incubated overnight at 4°C with the corresponding primary antibodies (diluted in TBST/+5% skimmed milk) and followed by incubation with HRP-conjugated anti-mouse, anti-guinea pig IgG-HRP, or anti-rabbit IgG secondary antibodies (Dako) for 1 h at room temperature. The detection was performed with chemiluminescence using Amersham ECL Prime Western blotting detection reagent (GE, Amersham) according to the manufacturers' instructions. The immunoreactive bands were visualized using a ChemiDoc MP (BioRad), and the immunoreactive signals were analyzed using the Image Lab Biorad software.
Nano-UPLC mass spectrometry system and phosphopeptide enrichment
Mass spectrometry experiments were performed from two independent biological replicates. Each replicate consisted of a pool of hippocampi (N = 4 mice/genotype) that were subsequently fractionated to obtain subsynaptic fractions, either from euploid or from trisomic mice. Following resuspension and protein quantification, subsynaptic fractions (70–200 μg) were diluted in 1% SDS and digested with the single sequence-specific protease trypsin, following the previously described FASP protocol (Wiśniewski et al., 2009). Digested peptides were subjected to phosphopeptide enrichment using the High-Select™ TiO2 Phosphopeptide Enrichment Kit (Thermo Fisher Scientific). About 45% of each enriched sample was analyzed using an Orbitrap Fusion Lumos with an EASY-Spray nanosource coupled to a nano-UPLC system (EASY-nanoLC 1000 liquid chromatograph) equipped with a 50-cm C18 column (EASY-Spray; 75 μm id, PepMap RSLC C18, 2-μm particles, 45°C). Chromatographic gradients started at 5% buffer B with a flow rate of 300 nl/min and gradually increased to 22% buffer B in 79 min and to 32% in 11 min. After each analysis, the column was washed for 10 min with 95% buffer B (buffer A: 0.1% formic acid in water and buffer B: 0.1% formic acid in acetonitrile). The mass spectrometer was operated in data-dependent acquisition mode, with full MS scans over a mass range of m/z 350–1500 with detection in the Orbitrap (120-K resolution) and with auto gain control (AGC) set to 100,000. In each cycle of data-dependent acquisition analysis, following each survey scan, the most intense ions above a threshold ion count of 10,000 were selected for fragmentation with HCD at normalized collision energy of 28%. The number of selected precursor ions for fragmentation was determined by the “Top Speed” acquisition algorithm (maximum cycle time of 3 s), and a dynamic exclusion of 60 s was set. Fragment ion spectra were acquired in the ion trap with an AGC of 10,000 and a maximum injection time of 200 ms.
Raw data processing and data analysis
Acquired data were analyzed using the Proteome Discoverer software suite (v2.0, Thermo Fisher Scientific), and the Mascot search engine (v2.5.1, Matrix Science) was used for peptide identification. Data were searched against a Mus musculus protein database derived from SwissProt and included the most common contaminants. A precursor ion mass tolerance of 7 ppm at the MS1 level was used, and up to three missed cleavages for trypsin were allowed. The fragment ion mass tolerance was set to 0.5 Da oxidation of methionine; N-terminal protein acetylation and phosphorylation in serine, threonine, and tyrosine were defined as variable modification; and carbamidomethylation of cysteines was set as fixed modification. The identified peptides were filtered by 5%FDR, and only those proteins identified in both biological replicates (for each condition) were retained for further analysis (Supplementary Table 1).
Bioinformatic analysis of proteomic data
The proteomic signature of the experimental groups was represented using Venn diagrams. These diagrams were created using VENNY 2.1 freeware (http://bioinfogp.cnb.csic.es/tools/venny/; Oliveros et al., 2007–2015). This graphical interface allowed the representation of proteins overlapped among groups, together with the visualization of differentially expressed proteins in the different groups/fractions.
To functionally classify the differentially expressed proteins from the different groups (subsynaptic fraction and genotype), the data were submitted to the Database for Annotation, Visualization, and Integrated Discovery (DAVID: http://david.abcc.ncifcrf.gov/). This software, based on the use of gene ontology (GO), was used for grouping the differentially expressed proteins into biological processes, molecular functions, and cellular components.
In Silico analysis of phosphoproteomic data
An updated database of murine kinases and phosphatases was generated before the analysis of the phosphoproteomic data. To this end, Uniprot database was used as the reference library, proteins were sorted (e.g., “kinase” and “phosphatase” criteria), and the protein list was further manually curated. This database was used to analyze potentially differentially expressed kinases and phosphatases between genotypes (Supplementary Table 2).
The analysis of the proteome and phosphoproteome signatures was performed using, as inclusion criteria, the identification of the proteins and phosphosites in both biological replicates. Furthermore, the conserved phosphosites were classified based on their qualitative expression in the different subsynaptic compartments, as well as in their relative expression levels. Regarding the latter, proteins and phosphosites detected in both fractions were classified based on the relative TS:EU ratio, resulting from the mean MS peak area (TS) (replicates 1 and 2)/mean MS peak area (EU) (replicates 1 and 2). Ratios < 0.75 or >1.5 were considered as downregulated and upregulated protein levels, respectively (see Supplementary Table 3).
Following phosphoproteomic profile identification of trisomic subsynaptic fractions, a systems biology approach was conducted to infer the potential kinases associated with the phosphoproteomic profile. To this end, the Human protein Reference Database (Hprd.org) was used, and the putative kinases involved in the phosphoprofile were identified. Finally, an interaction network analysis was performed using the Search Tool for the Retrieval of Interacting Genes/Proteins database (STRING; http://www.string-db.org/). This bioinformatic tool allowed the generation of protein–protein interaction networks, integrating the kinases/phosphatases and iGluR phosphoprofile of Ts65Dn subsynaptic fractions.
Statistical analysis
Comparison of immunolabeling signal intensities between euploid and trisomic hippocampal fractions was evaluated using the GraphPad software. All datasets passed the Kolmogorov–Smirnov normality test and were analyzed applying a two-tailed unpaired t-test with Welch's correction. Graph bars represent the mean ± SEM for each group.
Results
Subsynaptic proteomic signature of the Ts65Dn mice hippocampi
Neuronal communication mainly resides in the synapse, a highly specialized subcellular compartment juxtaposing crosstalking neurons. Within synapses, there is an additional morphological compartmentalization that allows the precise neurotransmitter and structural machinery positioning, underlying neuronal function. Consequently, potential specific changes on protein composition of the subsynaptic compartments can impact synaptic function and ultimately alter Ts65Dn glutamatergic neurons. To address this issue, we performed a subsynaptic fractionation protocol toward the obtention of extrasynaptic-, presynaptic-, and postsynaptic-enriched hippocampal fractions, both from euploid and from trisomic adult mice (Figure 1A).
Figure 1
Hippocampal subsynaptic fractionation protocol was validated by western blot, using specific subsynaptic compartment markers. As expected, extrasynaptic fraction showed an elevated signal of synaptophysin, whereas the postsynaptic fraction was enriched in scaffolding protein PSD95 (Figure 1A). Moreover, we observed a significant enrichment of the NMDA receptor subunits GluN1 and GluN2A in the postsynaptic fraction (Figure 1A). Overall, these biochemical studies confirmed the obtention of enriched subsynaptic compartments from murine brain tissue, and the same protocol was used to obtain enriched subsynaptic hippocampal fractions for mass spectrometry analysis.
Subsynaptic fractions (postsynaptic and extrasynaptic) were subjected to trypsin digestion, and mass spectrometry studies were performed. Data analysis was based on the comparison of the genotype effect on the subsynaptic proteomic/phosphoproteomic signature (defined as the overlapping profile of the biological replicates).
Proteomic analysis allowed the identification of 1922 proteins (M. musculus protein database, SwissProt) consistently present in both biological replicates (Figure 1B). Subsynaptic distribution analysis of these proteins showed that 938/1043 proteins are present in the postsynaptic fraction of euploid/trisomic mice, respectively, whereas 863/1059 proteins are present in the extrasynaptic fractions of euploid/trisomic mice, respectively (Figure 1C and Supplementary Table 1). Despite the proteomic profile showed the presence of conserved proteins in trisomic hippocampal synapses (60.7 and 57.8% protein coincidence for postsynaptic and extrasynaptic fractions, respectively; Figures 2A, 3A), significant protein composition changes were detected in trisomic hippocampi. Further in silico analyses were performed to analyze the potential functional outcomes of these proteomic changes.
Figure 2
Figure 3
Gene ontology enrichment analysis of subsynaptic Ts65Dn proteomic signature
The genotype effect on protein composition of the postsynaptic and extrasynaptic fractions was performed following the GO-term-based protein enrichment analysis. For functional annotation studies, only those proteins showing differential subsynaptic fraction expression were considered (Figures 2A, 3A). Statistical enrichment analysis of “cellular compartment,” “biological processes,” and “molecular functions” revealed the presence of significant genotype-dependent differences, both in the postsynaptic and in the extrasynaptic fractions (Figures 2B–E, 3B–E).
GO-term analysis of the euploid-specific postsynaptic proteome indicated the presence of an enriched representation of dendrite-associated proteins, thereby reflecting a dendritic protein under-representation in the trisomic postsynaptic fraction (Figures 2B–D). Moreover, GO analysis of the trisomic postsynaptic fraction showed significant changes in the following cellular functions: membranes, cytoplasm, transport, and protein binding (Figures 2C–E).
Regarding the extrasynaptic fractions, GO-term analysis of the euploid-specific extrasynaptic fraction showed the presence of mitochondria and protein phosphatase-binding proteins, indicating their relative absent detection in the trisomic extrasynaptic fraction (Figures 3B–D). In contrast, GO analysis of trisomic extrasynaptic fraction indicated changes in membrane, cytoplasm, synapse transport, phosphorylation, protein kinase activity, and protein binding (Figure 3C–E). Overall, GO-term analysis suggested that the subsynaptic proteome of the trisomic hippocampus is altered, in terms of both cellular compartments and functional outcomes. In particular, the analysis strongly suggested that phosphorylation-associated processes could be affected, as described below.
Altered expression of kinases and phosphatases in the Ts65Dn mouse hippocampal subsynaptic fractions
Based on the functional annotation analyses showing alterations in phosphorylation, protein kinase activity, and phosphatase binding in the Ts65Dn synapses, we aimed to elucidate potential expression changes in kinases and phosphatases. To this end, an updated mouse kinome and phosphatome database was developed, using Uniprot as the reference database and filtering the M. musculus protein entries for the extraction of mouse kinases and phosphatases (Supplementary Table 2). Afterward, the genotype effect was compared either for postynaptic or for extrasynaptic fractions. With respect to euploid mice, the postsynaptic fraction of trisomic mice showed the presence of a subset of specific/enriched kinases and phosphatases (50 and 10 proteins, respectively), together with nondetectable/reduced expression of kinases and phosphatases (21 and 7 proteins, respectively) (Figure 4 and Supplementary Figure 1A). The comparative study of the extrasynaptic fractions indicated the presence, in trisomic samples, of a subset of specific/enriched kinases and phosphatases (37 and 7 proteins, respectively), together with non-detectable/reduced expression of kinases and phosphatases that are present in euploid samples (five and eight proteins, respectively) (Figure 5 and Supplementary Figure 1B). We performed western blot analysis of MS-identified kinases with putative differential expression to validate these unbiased findings. In agreement with proteomic data, western blot showed a significant increase of CamKII expression levels in the postsynaptic fraction of trisomic mice (p = 0.033, two-tailed Student's t-test), together with a significant decrease of Fyn kinase levels in the postsynaptic and a concomitant increase in the extrasynaptic fraction of the Ts65Dn mice (p = 0.015 and 0.024, respectively; two-tailed Student's t-test; Supplementary Figure 2A).
Figure 4
Figure 5
Further, we aimed to predict the potential consequences of the dysregulated trisomic synaptic kinome. We focused in silico analysis on the ionotropic family of glutamate receptors similar to the glutamatergic alterations reported in the Ts65Dn mouse model. The potential protein–protein interaction networks between iGluRs and protein kinases upregulated in the trisomic fractions (either postsynaptic or extrasynaptic) were analyzed using STRING software. Interactome analysis of the trisomic postsynaptic fraction unraveled significant kinase-iGluR interaction, as shown by the elevated number of edges (predicted number of edges: 123 vs. expected number of edges: 18) and the PPI enrichment (p-value = 0), indicating the significant interactions within the PPI network. In particular, STRING analysis identified protein kinases from the CaMKII family, Dlgap1, Agap2, MRCKβ distinctively interacting with the detected iGluR subunits (GluN1, GluN2A, GluN2B, GluN2D, GluA1, and GluA2; Figure 4). Interactome analysis of the trisomic extrasynaptic fraction unraveled significant kinase-iGluR interaction, as shown by the significant increase of the network edges (predicted number of edges: 51 vs. expected number of edges: 21) and the PPI enrichment (p-value = 2.49.10−8), indicating the significant interactions of the PPI network. Moreover, several kinases (CaMKII kinase family members, cAMP-dependent protein kinase, PKC family, and Agap2) distinctively interacting with GriA1 and GriA2 subunits of the AMPARs were also significantly detected (Figure 5).
Altered iGluR subsynaptic phosphopattern in the Ts65Dn mouse hippocampus
Hippocampal subsynaptic fractions were analyzed by phosphoenrichment-coupled mass spectrometry to evaluate the predictive kinome-iGluR interactome model. The phosphoproteome of the different fractions was analyzed, and a total of 59 iGluR phosphosites were identified in the different biological replicates (Supplementary Table 3). From these, only those iGluR phosphosites found in both replicates were retained for further genotype effect analysis (Figure 6 and Supplementary Table 3). This technique allowed the identification of several iGluR phosphosites, mostly corresponding to previously reported phosphoresidues (15 individual phosphosites); two novel phosphosites were identified. Importantly, the subsynaptic iGluR phosphopattern was characterized in the hippocampus of Ts65Dn mice. Using this unbiased method, the identified upregulated phosphosites (displaying higher MS peaks in trisomic samples) showed an elevated representation of NMDAR subunits (eight phosphosites for GluN2A, seven for GluN2B, and one for GluN1), whereas only one phosphosite was detected in GluA1 subunit (AMPAR subunit) and no phosphosites were identified in KAR subunits. We performed western blot analysis on selected MS-identified differentially expressed phosphosites to further validate our results. This analysis showed an increase of both GluN1(pS890) and GluN2B(pS1284) levels in the postsynaptic fraction of trisomic mice, although not significant, in agreement with the changes observed along unbiased phosphoproteomic analysis (Supplementary Figure 2B). Finally, the kinases putatively phosphorylating the consensus iGluR phosphoresidues were in silico predicted using “Phosphomotif Finder” software. This analysis predicted a potential enzymatic activity of CaMKII, PKC, CK1 on the identified iGluR phosphoresidues, which would be in agreement with their above mentioned upregulated expression.
Figure 6
Discussion
In the present study, we have performed a phosphoenrichment-coupled mass spectrometry analysis that showed altered functional changes in the hippocampal extrasynaptic and postsynaptic fractions of the Ts65Dn mice, a trisomic mouse model recapitulating DS phenotypic alterations. The majority of kinases interacting with iGluRs remarkably show an overall increase in trisomic samples, whereas the majority of iGluR-related phosphatases show reduction, with some exceptions in both protein categories. These results are concomitant with the increased number of iGluR phosphopeptides identified by the phosphoproteomic analysis of Ts65Dn mice hippocampal postsynaptic fractions. Importantly, in silico kinome-iGluR phosphopattern analysis supported the role of specific kinases and their potential iGluR phosphosites in the neuronal pathophysiology of Ts65Dn mouse.
The initial proteome analysis showed a relative conservation of the number of proteins (roughly 1,000 proteins identified), regardless of the subsynaptic fraction and genotype. The comparative study of the genotype effect on the subsynaptic fraction composition showed a protein coincidence of 60.7 and 57.8% for the postsynaptic and extrasynaptic fractions, respectively. This quantitative analysis indicates a different molecular composition of the Ts65Dn hippocampal synapses, which might result from the overexpression of triplicated genes. Hence, several proteins are differentially expressed in the subsynaptic compartments of Ts65Dn hippocampal neurons, potentially reflecting the molecular substrates associated with Ts65Dn synaptic abnormalities. Interestingly, further GO-term-based protein enrichment analysis supported functional outcomes of the proteomic changes detected in the postsynaptic and extrasynaptic fractions of trisomic mice. In particular, the postsynaptic fraction of trisomic mice showed a relative decrease of proteins related with “dendritic function” GO term. This finding is in agreement with previous studies showing a reduction of dendritic branching and spine density in Ts65Dn mice (Belichenko et al.,
Besides the differential enrichment of structure- and energy-associated proteins, a significant alteration of phosphorylation-associated proteins was detected in the trisomic subsynaptic compartments. Indeed, the extrasynaptic trisomic fraction showed an enrichment of “protein phosphorylation” and “protein kinase activity” GO terms, while “protein phosphatase binding” (enriched in euploid extrasynaptic fraction) was significantly reduced. Although these evidences were suggestive of an altered enzymatic activity of kinases/phosphatases in Ts65Dn hippocampi, a more detailed analysis of kinase/phosphatase subsynaptic expression levels was performed. This analysis roughly showed an increase of kinase expression, together with a reduction of phosphatase levels in the Ts65Dn mouse model. In our study, further detailed analysis of the trisomic kinome showed increased levels of CaMKII subunits (alpha, beta, delta, and gamma subunits) in hippocampal postsynaptic and extrasynaptic fractions of Ts65Dn mice, by proteomics and WB techniques. This increase could result in an enhanced CamKII activity, similar to previous studies showing increased phosphorylated CaMKII in Ts65Dn hippocampi (Siarey et al., 2006). In agreement with this kinase increase, previous work of Fernandez and collaborators showed an increased amount of phosphopeptides in crude synaptosomes from Ts65Dn mice (Fernandez et al.,
Interactome analysis revealed significant interactions between protein kinases enriched in trisomic subsynaptic fractions and iGluR subunits. In the postsynaptic fraction, these kinases interact with both NMDAR and AMPAR subunits. In contrast, the interactome analysis only revealed potential protein–protein interactions of these kinases with GluA1 and GluA2 subunits of the AMPAR, which are proportionally more largely represented than NMDAR subunits in the extrasynaptic fraction. Following this predictive analysis, we performed an accurate characterization of the iGluR phosphopattern in hippocampal subsynaptic fractions of Ts65Dn mice. Further bioinformatic analyses were focused on the phosphorylation of these receptors to unveil potential changes on iGluR phosphopattern. In sample preparation, the fractions followed tryptic digestion with a single protease (trypsin) for further unbiased MS method. Technically, the use of a single protease limited the extent of detectable phosphopeptides (Dephoure et al.,
Furthermore, some of the phosphosites found to be over-represented in our study had previously been described to be involved in synaptic plasticity. For example, GluN1 S890 phosphorylation disperses the surface clusters of GluN1 (Tingley et al., 1997). GluN2A S929 phosphorylation modulates desensitization of GluN2A/NMDA receptors (Maki et al.,
In summary, our study revealed an altered subsynaptic composition of the adult Ts65Dn mouse hippocampus. These molecular subsynaptic alterations might functionally be underlying the Ts65Dn-associated dendritic deterioration, altered synapse organization, and mitochondrial dysfunction. Furthermore, this phosphoproteomic study provides an extensive repertoire of the iGluR phosphosites, together with the respective predicted kinases. Importantly, although the mechanistic studies remain elusive, the characterization of the iGluR phosphopattern and the predicted associated kinases in the Ts65Dn hippocampi may help to identify novel therapeutic targets for DS synaptopathy.
Statements
Author contributions
XA and MG contributed to the conception and design of the study. MG performed the experimental and bioinformatic analyses, with the technical help from CG. MG, FC, and XA wrote the manuscript. All authors contributed to the revision of the manuscript revision and also read and approved the submitted version.
Funding
This work was supported by the Grants La Marató (Project N. 20140210), PI16/00851 (ISCIII), PCIN-2014-105 (MINECO) and Miguel Servet Program (CPII16/00021, ISCIII) to XA. MG and CG contract were funded by Fundació La Marató (Project N. 20140210).
Acknowledgments
The authors would like to thank G. Gou (Hospital de la Santa Creu i de Sant Pau-UAB, Barcelona, Spain) and X. Morató and M. López (Unitat Neurofarmacologia, Universitat Barcelona, Spain) for their technical support in the synaptosomes preparation. We also acknowledge Dr. M. López Heredia (Unitat Genètica, IDIBELL), G. Gou, and Dr. À. Bayés (Hospital de la Santa Creu i de Sant Pau-UAB, Barcelona, Spain) for guidance in bioinformatic analysis. We thank Dr. Eduard Sabidó and Cristina Chiva (Proteomics Unit, Center for Genomic Regulation, Barcelona, Spain) for technical help in MS studies and data managing.
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2018.00226/full#supplementary-material
References
1
AltafajX.MartínE. D.Ortiz-AbaliaJ.ValderramaA.Lao-PeregrínC.DierssenM.et al. (2013). Normalization of Dyrk1A expression by AAV2/1-shDyrk1A attenuates hippocampal-dependent defects in the Ts65Dn mouse model of Down syndrome. Neurobiol. Dis.52, 117–127. 10.1016/j.nbd.2012.11.017
2
AltafajX.Ortiz-AbaliaJ.FernándezM.PotierM. C.LaffaireJ.AndreuN.et al. (2008). Increased NR2A expression and prolonged decay of NMDA-induced calcium transient in cerebellum of TgDyrk1A mice, a mouse model of Down syndrome. Neurobiol. Dis.32, 377–384. 10.1016/j.nbd.2008.07.024
3
AntonarakisS. E.LyleR.DermitzakisE. T.ReymondA.DeutschS. (2004). Chromosome 21 and down syndrome: from genomics to pathophysiology. Nat. Rev. Genet.5, 725–738. 10.1038/nrg1448
4
BelichenkoP. V. (2004). Synaptic structural abnormalities in the Ts65Dn mouse model of down syndrome. J. Comp. Neurol.480, 281–298. 10.1002/cne.20337
5
BelichenkoP. V.KleschevnikovA. M.MasliahE.WuC.Takimoto-KimuraR.SalehiA.et al. (2009). Excitatory-inhibitory relationship in the fascia dentata in the Ts65Dn mouse model of down syndrome. J. Comp. Neurol.512, 453–466. 10.1002/cne.21895
6
BelichenkoP. V.KleschevnikovA. M.SalehiA.EpsteinC. J.MobleyW. C. (2007). Synaptic and cognitive abnormalities in mouse models of down syndrome: exploring genotype-phenotype relationships. J. Comp. Neurol.504, 329–345. 10.1002/cne.21433.
7
Benavides-PiccioneR.Ballesteros-YáñezI.de LagránM. MElstonG.EstivillX.FillatC.et al. (2004). On dendrites in Down syndrome and DS murine models: a spiny way to learn. Prog. Neurobiol. 74, 111–126. 10.1016/j.pneurobio.2004.08.001
8
BestT. K.CramerN. P.ChakrabartiL.HaydarT. F.GaldzickiZ. (2012). Dysfunctional hippocampal inhibition in the Ts65Dn mouse model of Down syndrome. Exp. Neurol.233, 749–757. 10.1016/j.expneurol.2011.11.033
9
BoadaR.Hutaff-LeeC.SchraderA.WeitzenkampD.BenkeT. A.GoldsonE. J.et al. (2012). Antagonism of NMDA receptors as a potential treatment for Down syndrome: a pilot randomized controlledTrial. 2, e141–e111. 10.1038/tp.2012.66
10
BraithwaiteS. P.AdkissonM.LeungJ.NavaA.MastersonB.UrferR.et al. (2006). Regulation of NMDA receptor trafficking and function by striatal-enriched tyrosine phosphatase (STEP). Eur. J. Neurosci.23, 2847–2856. 10.1111/j.1460-9568.2006.04837.x
11
ChanS. F.SucherN. J. (2001). An NMDA receptor signaling complex with protein phosphatase 2A. J. Neurosci.21, 7985–7992. 10.1523/JNEUROSCI.21-20-07985.2001
12
CherguiK. (2005). Physiological role for Casein Kinase 1 in glutamatergic synaptic transmission. J. Neurosci.25, 6601–6609. 10.1523/JNEUROSCI.1082-05.2005
13
ContestabileA.MagaraS.CanceddaL. (2017). The GABAergic hypothesis for cognitive disabilities in down syndrome. Front. Cell. Neurosci.11:54. 10.3389/fncel.2017.00054
14
CostaA. C.Scott-McKeanJ. J.StaskoM. R. (2007). Acute injections of the NMDA receptor antagonist memantine rescue performance deficits of the Ts65Dn mouse model of down syndrome on a fear conditioning test. Neuropsychopharmacology33, 1624–1632. 10.1038/sj.npp.1301535
15
CostaA. C.WalshK.DavissonM. T. (1999). Motor dysfunction in a mouse model for Down syndrome. Physiol. Behav.68, 211–220.
16
CramerN.GaldzickiZ. (2012). From abnormal hippocampal synaptic plasticity in down syndrome mouse models to cognitive disability in down syndrome. Neural Plast.2012, 1–12. 10.1155/2012/101542
17
DephoureN.GouldK. L.GygiS. P.KelloggD. R. (2013). Mapping and analysis of phosphorylation sites: a quick guide for cell biologists. Mol. Biol. Cell24, 535–542. 10.1091/mbc.E12-09-0677
18
DierssenM. (2012). Down syndrome: the brain in trisomic mode. Nat. Rev. Neurosci.13, 844–858. 10.1038/nrn3314
19
FernandezF.TrinidadJ. C.BlankM.FengD. D.BurlingameA. L.GarnerC. C. (2009). Normal protein composition of synapses in Ts65Dn mice: a mouse model of Down syndrome. J. Neurochem.110, 157–169. 10.1111/j.1471-4159.2009.06110.x
20
GaldzickiZ.SiareyR. J. (2003). Understanding mental retardation in Down's syndrome using trisomy 16 mouse models. Genes Brain Behav.2, 167–178. 10.1034/j.1601-183X.2003.00024.x
21
Goebel-GoodyS. M.DaviesK. D.LingerR. M.FreundR. K.BrowningM. D. (2009). Phospho-regulation of synaptic and extrasynaptic N-methyl-d-aspartate receptors in adult hippocampal slices. NSC158, 1446–1459. 10.1016/j.neuroscience.2008.11.006
22
GrauC.AratóK.Fernández-FernándezJ. M.ValderramaA.SindreuC.FillatC.et al. (2014). DYRK1A-mediated phosphorylation of GluN2A at Ser(1048) regulates the surface expression and channel activity of GluN1/GluN2A receptors. Front. Cell. Neurosci.8:331. 10.3389/fncel.2014.00331
23
HansonJ. E.WeberM.MeilandtW. J.WuT.LuuT.DengL.et al. (2013). GluN2B antagonism affects interneurons and leads to immediate and persistent changes in synaptic plasticity, oscillations, and behavior. Neuropsychopharmacology38, 1221–1233. 10.1038/npp.2013.19
24
HardinghamG. E.BadingH. (2010). Synaptic versus extrasynaptic NMDA receptor signalling: implications for neurodegenerative disorders. Nat. Rev. Neurosci.11, 682–696. 10.1038/nrn2911
25
HayashiM. K.TangC.VerpelliC.NarayananR.StearnsM. H.XuR. M.et al. (2009). The Postsynaptic density proteins homer and shank form a polymeric network structure. Cell137, 159–171. 10.1016/j.cell.2009.01.050
26
IzzoA.MancoR.BonfiglioF.CalìG.De CristofaroT.PatergnaniS.et al. (2014). NRIP1/RIP140 siRNA-mediated attenuation counteracts mitochondrial dysfunction in Down syndrome. Hum. Mol. Genet.23, 4406–4419. 10.1093/hmg/ddu157
27
KleschevnikovA. M. (2012). Deficits in Cognition and synaptic plasticity in a mouse model of down syndrome ameliorated by GABAB receptor antagonists. J. Neurosci.32, 9217–9227. 10.1523/JNEUROSCI.1673-12.2012
28
LanJ. Y.SkeberdisV. A.JoverT.GroomsS. Y.LinY.AranedaR. C.et al. (2001). Protein kinase C modulates NMDA receptor trafficking and gating. Nat. Publ. Group4, 382–390. 10.1038/86028
29
LauC. G.ZukinR. S. (2007). NMDA receptor trafficking in synaptic plasticity and neuropsychiatric disorders. Nat. Rev. Neurosci.8, 413–426. 10.1038/nrn2153
30
LeeH.-K. (2006). Synaptic plasticity and phosphorylation. Pharmacol. Ther.112, 810–832. 10.1016/j.pharmthera.2006.06.003
31
LiJ.WilkinsonB.ClementelV. A.HouJ.O'DellT. J.CobaM. P. (2016). Long-term potentiation modulates synaptic phosphorylation networks and reshapes the structure of the postsynaptic interactome. Sci. Signal.9:rs8. 10.1126/scisignal.aaf6716
32
LiaoG. Y.WagnerD. A.HsuM. H.LeonardJ. P. (2001). Evidence for direct protein kinase-C mediated modulation of N-methyl-D-aspartate receptor current. Mol. Pharmacol.59, 960–964. 10.1124/mol.59.5.960
33
LockrowJ.BogerH.Bimonte-NelsonH.GranholmA. C. (2011). Effects of long-term memantine on memory and neuropathology in Ts65Dn mice, a model for Down syndrome. Behav. Brain Res.221, 610–622. 10.1016/j.bbr.2010.03.036
34
LundbyA.SecherA.LageK.NordsborgN. B.DmytriyevA.LundbyC.et al. (2012). Quantitative maps of protein phosphorylation sites across 14 different rat organs and tissues. Nat. Commun.3:876. 10.1038/ncomms1871
35
LussierM. P.Sanz-ClementeA.RocheK. W. (2015). Dynamic Regulation of N-Methyl-d-aspartate (NMDA) and α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic Acid (AMPA) Receptors by Posttranslational Modifications. J. Biol. Chem.290, 28596–28603. 10.1074/jbc.R115.652750
36
MakiB. A.ColeR.PopescuG. K. (2013). Two serine residues on GluN2A C-terminal tails control NMDA receptor current decay times. Channels7, 126–132. 10.4161/chan.23968
37
McQueenJ.RyanT. J.McKayS.MarwickK.BaxterP.CarpaniniS. M.et al. (2017). Pro-death NMDA receptor signaling is promoted by the GluN2B C-terminus independently of Dapk1. Elife6:e17161. 10.7554/eLife.17161
38
MuntonR. P.Tweedie-CullenR.Livingstone-ZatchejM.WeinandyF.WaidelichM.LongoD.et al. (2007). Qualitative and quantitative analyses of protein phosphorylation in naive and stimulated mouse synaptosomal preparations. Mol. Cell. Proteomics6, 283–293. 10.1074/mcp.M600046-MCP200
39
MurphyJ. A.SteinI. S.LauC. G.PeixotoR. T.AmanT. K.KanekoN.ZukinR. S.et al. (2014). Phosphorylation of Ser1166 on GluN2B by PKA is critical to synaptic NMDA receptor function and Ca2+ signaling in spines. J. Neurosci.34, 869–879. 10.1523/JNEUROSCI.4538-13.2014
40
OliverosJ. C. (2007–2015). Venny. An Interactive Tool for Comparing Lists with Venn's Diagrams. Available online at: http://bioinfogp.cnb.csic.es/tools/venny/index.html
41
PahlS.TapkenD.HaeringS. C.HollmannM. (2014). Trafficking of kainate receptors. Membranes4, 565–595. 10.3390/membranes4030565
42
PaolettiP.BelloneC.ZhouQ. (2013). NMDA receptor subunit diversity: impact on receptor properties, synaptic plasticity and disease. Nat. Rev. Neurosci. 14, 383–400. 10.1038/nrn3504
43
ParkerC. E.MocanuV.MocanuM.DichevaN.WarrenM. (2010). Mass spectrometry For post-translational modifications, in Neuroproteomics, ed AlzateO. (Boca Raton, FL: CRC Press/Taylor & Francis).
44
PhillipsG. R.HuangJ. K.WangY.TanakaH.ShapiroL.ZhangW.et al. (2001). The presynaptic particle web: ultrastructure, composition, dissolution, and reconstitution. Neuron32, 63–77. 10.1016/S0896-6273(01)00450-0
45
RisserD.LubecG.CairnsN.Herrera-MarschitzM. (1997). Excitatory amino acids and monoamines in parahippocampal gyrus and frontal cortical pole of adults with Down syndrome. Life Sci.60, 1231–1237.
46
RuedaN.FlórezJ.Martínez-CuéC. (2012). Mouse models of down syndrome as a tool to unravel the causes of mental disabilities. Neural Plast.2012, 1–26. 10.1155/2012/584071
47
SandersonJ. L.GorskiJ. A.Dell'AcquaM. L. (2016). NMDA receptor-dependent LTD requires transient synaptic incorporation of Ca 2+ -permeable AMPARs mediated by AKAP150-Anchored PKA and calcineurin. Neuron89, 1000–1015. 10.1016/j.neuron.2016.01.043
48
Sanz-ClementeA.GrayJ. A.OgilvieK. A.NicollR. A.RocheK. W. (2013). Activated CaMKII couples GluN2B and casein kinase 2 to control synaptic NMDA receptors. Cell Rep.3, 607–614. 10.1016/j.celrep.2013.02.011
49
Scott-McKeanJ. J.CostaA. C. (2011). Exaggerated NMDA mediated LTD in a mouse model of Down syndrome and pharmacological rescuing by memantine. Learn. Mem.18, 774–778. 10.1101/lm.024182.111
50
SeidlR.CairnsN.SingewaldN.KaehlerS. T.LubecG. (2001). Differences between GABA levels in Alzheimer's disease and Down syndrome with Alzheimer-like neuropathology. Naunyn Schmiedebergs Arch. Pharmacol.363, 139–145. 10.1007/s002100000346
51
SiareyR. J.Kline-BurgessA.ChoM.BalboA.BestT. K.HarashimaC.et al. (2006). Altered signaling pathways underlying abnormal hippocampal synaptic plasticity in the Ts65Dn mouse model of Down syndrome. J. Neurochem.98, 1266–1277. 10.1111/j.1471-4159.2006.03971.x
52
SiddiquiA.LacroixT.StaskoM. R.Scott-McKeanJ. J.CostaA. C.GardinerK. J. (2008). Molecular responses of the Ts65Dn and Ts1Cje mouse models of Down syndrome to MK-801. Genes Brain Behav.7, 810–820. 10.1111/j.1601-183X.2008.00428.x
53
TingleyW. G.EhlersM. D.KameyamaK.DohertyC.PtakJ. B.RileyC. T.et al. (1997). Characterization of protein kinase A and protein kinase C phosphorylation of the N-methyl-D-aspartate receptor NR1 subunit using phosphorylation site-specific antibodies. J. Biol. Chem. 272, 5157–5166.
54
TraynelisS. F.WollmuthL. P.McBainC. J.MennitiF. S.VanceK. M.OgdenK. K.et al. (2010). Glutamate receptor ion channels: structure, regulation, and function. Pharmacol. Rev.62, 405–496. 10.1124/pr.109.002451
55
TrinidadJ. C.BarkanD. T.GulledgeB. F.ThalhammerA.SaliA.SchoepferR.et al. (2012). Global identification and characterization of both O-GlcNAcylation and phosphorylation at the murine synapse. Mol. Cell. Proteomics11, 215–229. 10.1074/mcp.O112.018366
56
WangJ. Q.LiuX.ZhangG.ParelkarN. K.AroraA.HainesM.et al. (2006). Phosphorylation of glutamate receptors: a potential mechanism for the regulation of receptor function and psychostimulant action. J. Neurosci. Res.84, 1621–1629. 10.1002/jnr.21050
57
WeitzdoerferR.DierssenM.FountoulakisM.LubecG. (2001). Fetal life in Down syndrome starts with normal neuronal density but impaired dendritic spines and synaptosomal structure. J. Neural Transm.61, 59–70. 10.1007/978-3-7091-6262-0_5
58
WestphalR. S. (1999). Regulation of NMDA receptors by an associated phosphatase-kinase signaling complex. Science285, 93–96. 10.1126/science.285.5424.93
59
WiśniewskiJ. R.NagarajN.ZougmanA.GnadF.MannM. (2010). Brain phosphoproteome obtained by a FASP-based method reveals plasma membrane protein topology. J. Proteome Res.9, 3280–3289. 10.1021/pr1002214
60
WiśniewskiJ. R.ZougmanA.NagarajN.MannM. (2009). Universal sample preparation method for proteome analysis. Nat. Methods6, 359–362. 10.1038/nmeth.1322
61
ZhangF.GuoA.LiuC.CombM.HuB. (2013). Phosphorylation and assembly of glutamate receptors after brain ischemia. Stroke44, 170–176. 10.1161/STROKEAHA.112.667253
Summary
Keywords
Down syndrome, ionotropic glutamate receptors, NMDA receptor, kinases, proteomics, phosphoproteomics
Citation
Gómez de Salazar M, Grau C, Ciruela F and Altafaj X (2018) Phosphoproteomic Alterations of Ionotropic Glutamate Receptors in the Hippocampus of the Ts65Dn Mouse Model of Down Syndrome. Front. Mol. Neurosci. 11:226. doi: 10.3389/fnmol.2018.00226
Received
01 March 2018
Accepted
11 June 2018
Published
25 July 2018
Volume
11 - 2018
Edited by
Inmaculada Maria Gonzalez-Gonzalez, University of Central Lancashire, United Kingdom
Reviewed by
Antonio Sanz-Clemente, Feinberg School of Medicine, Northwestern University, United States; Mariana Vargas-Caballero, University of Southampton, United Kingdom
Updates

Check for updates
Copyright
© 2018 Gómez de Salazar, Grau, Ciruela and Altafaj.
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) and the copyright owner(s) 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: Xavier Altafaj xaltafaj@idibell.cat
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.