Abigail U. Carbonell1†
Carmen Freire-Cobo1†Ilana V. Deyneko1
Saunil Dobariya1
Hediye Erdjument-Bromage2Amy E. Clipperton-Allen3Damon T. Page3
Thomas A. Neubert2
Bryen A. Jordan1,4*- 1Dominick P. Purpura Department of Neuroscience, Albert Einstein College of Medicine, Bronx, NY, United States
- 2Department of Neuroscience and Physiology, New York University Grossman School of Medicine, New York, NY, United States
- 3Department of Neuroscience, The Scripps Research Institute Florida, Jupiter, FL, United States
- 4Department of Psychiatry and Behavioral Sciences, Albert Einstein College of Medicine, Bronx, NY, United States
Specific and effective treatments for autism spectrum disorder (ASD) are lacking due to a poor understanding of disease mechanisms. Here we test the idea that similarities between diverse ASD mouse models are caused by deficits in common molecular pathways at neuronal synapses. To do this, we leverage the availability of multiple genetic models of ASD that exhibit shared synaptic and behavioral deficits and use quantitative mass spectrometry with isobaric tandem mass tagging (TMT) to compare their hippocampal synaptic proteomes. Comparative analyses of mouse models for Fragile X syndrome (Fmr1 knockout), cortical dysplasia focal epilepsy syndrome (Cntnap2 knockout), PTEN hamartoma tumor syndrome (Pten haploinsufficiency), ANKS1B syndrome (Anks1b haploinsufficiency), and idiopathic autism (BTBR+) revealed several common altered cellular and molecular pathways at the synapse, including changes in oxidative phosphorylation, and Rho family small GTPase signaling. Functional validation of one of these aberrant pathways, Rac1 signaling, confirms that the ANKS1B model displays altered Rac1 activity counter to that observed in other models, as predicted by the bioinformatic analyses. Overall similarity analyses reveal clusters of synaptic profiles, which may form the basis for molecular subtypes that explain genetic heterogeneity in ASD despite a common clinical diagnosis. Our results suggest that ASD-linked susceptibility genes ultimately converge on common signaling pathways regulating synaptic function and propose that these points of convergence are key to understanding the pathogenesis of this disorder.
Introduction
Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by social-communication deficits and restrictive and repetitive behaviors. Although the specific cause of ASD is unknown, autism is highly heritable, with a monozygotic twin concordance rate of 40–80% (Gaugler et al., 2014). The genetic architecture of ASD is extraordinarily complex, with common inherited variants and rare de novo mutations working together to confer genetic risk (Weiner et al., 2017). This polygenic etiology presents challenges for elucidating the molecular pathogenesis of autism. However, ASD-related syndromes with defined genetic causes for autistic phenotypes present the best opportunities for elucidating the underlying mechanisms of ASD and identifying possible therapeutic targets (Sztainberg and Zoghbi, 2016). These monogenic syndromes, including Fragile X syndrome (FMR1), Rett syndrome (MECP2), PTEN hamartoma tumor syndrome (PTEN), tuberous sclerosis complex (TSC1/TSC2), Phelan McDermid syndrome (SHANK3), and cortical dysplasia focal epilepsy syndrome (CNTNAP2), play a key role in furthering our understanding of ASD.
Animal models of these syndromes have proven essential for studying the underlying neuropathology of ASD, especially changes in the complex processes of mammalian brain development and function. Studies in rodents have shown that ASD risk genes converge on transcriptional regulation, protein homeostasis, and synaptic structure and function (De Rubeis et al., 2014; Pinto et al., 2014; Ruzzo et al., 2019). Accordingly, mice demonstrating the loss of Fmr1, Mecp2, Pten, Tsc1/Tsc2, Shank genes, or the Nrxn and Nlgn families all demonstrate changes in synaptic excitability or plasticity, and most also show altered dendritic growth or spine dynamics (Hulbert and Jiang, 2016; Parikshak et al., 2016; Varghese et al., 2017; Verma et al., 2019). Non-syndromic ASD models, such as those induced by valproic acid or maternal immune activation, also reveal structural and functional synaptic deficits, showing that environmental factors can lead to similar synaptic phenotypes (Sui and Chen, 2012; Martin and Manzoni, 2014; Cellot et al., 2016; Patrich et al., 2016; Li et al., 2018; Wang et al., 2018; Andoh et al., 2019).
We recently characterized ANKS1B haploinsufficiency, a new genetic syndrome presenting with autism and other neurodevelopmental disorders that is caused by a monogenic deletion in the ANKS1B gene (Carbonell et al., 2019). This gene, which encodes for the protein AIDA-1, was previously identified in ASD risk gene networks (Li et al., 2014). To facilitate our work, we developed a mouse model for ANKS1B haploinsufficiency and found that it recapitulates behavioral correlates of the syndrome (Carbonell et al., 2019). Altogether, we found that AIDA-1 regulates activity-induced protein synthesis (Jordan et al., 2007), hippocampal synaptic plasticity (Tindi et al., 2015), and NMDA receptor subunit composition (Tindi et al., 2015; Carbonell et al., 2019). AIDA-1 is a core protein of the postsynaptic density and interacts with PSD95 in a complex that contains other factors associated with neurodevelopmental disorders, including Grin2b, Syngap1, and Nlgn (Kaizuka and Takumi, 2018; Carbonell et al., 2019). Our research provides an additional example corroborating the idea that molecular mechanisms regulating synaptic function underlie ASD pathobiology (Zoghbi and Bear, 2012; Bourgeron, 2015).
Despite the heterogeneous genetic architecture of autism and a complex etiology with contributions from environmental factors, ASD is diagnosed by distinct clinical criteria. Therefore, convergent cellular processes at the circuit, synaptic, or molecular level could underlie these shared behavioral phenotypes (Sestan and State, 2018). In seeking convergent mechanisms among syndromic and non-syndromic forms of autism, comparative studies often narrowly focus on selected behaviors, specific synaptic phenomenology, or shared responses to preclinical pharmacological interventions (Barnes et al., 2015; Heise et al., 2018; Schoen et al., 2019). While broader comparisons have been made using transcriptomic analyses (Forés-Martos et al., 2019; Quesnel-Vallières et al., 2019; Hernandez et al., 2020), these changes may not reflect mechanisms of disease due to multiple downstream levels of regulation such as protein translation, degradation, and transport (Vogel and Marcotte, 2012). Indeed, widespread discrepancies between transcript and protein abundance are well-known (Wang et al., 2017; Hoogendijk et al., 2019) but typically ignored when drawing conclusions from discovery-based transcriptomic screens. Proteomic approaches can therefore yield dramatically different results from transcriptomic profiles, as shown in a mouse model of Rett syndrome (Pacheco et al., 2017), and are particularly favored for ASD-linked genes that regulate protein translation (Fmr1, Tsc1/Tsc2, Pten) or degradation (Ube3a; Louros and Osterweil, 2016). ASD models in which synaptic scaffolding and membrane localization are altered can also be investigated at the proteomic level, especially when the synaptic compartments and complexes are isolated by fractionation or immunoprecipitation (Lee et al., 2017; Reim et al., 2017; Brown et al., 2018; Murtaza et al., 2020). Indeed, analysis of postmortem brain tissue from ASD patients showed that some transcriptomic changes were not observed while changes in synaptic proteins predicted by ASD risk genes and animal models were confirmed (Abraham et al., 2019).
Here, we compare the postsynaptic proteomes of five mouse models for autism using a particularly rigorous quantitative proteomic method employing isobaric tags. The models selected display synaptic deficits and represent models of ASD-related syndromes Fragile X syndrome (Fmr1 -/Y), PTEN hamartoma tumor syndrome (Pten +/−), cortical dysplasia focal epilepsy syndrome (Cntnap2 −/−), ANKS1B syndrome (Anks1b +/−), and the BTBR+ inbred model of idiopathic autism. These models have demonstrated face validity for autism, displaying hallmark behavioral correlates of ASD including social interaction deficits and restrictive behaviors (Kazdoba et al., 2016; Zhou et al., 2016; Kabitzke et al., 2018; Carbonell et al., 2019). In each model, we find evidence for upstream regulators that alter synaptic composition with predicted functional effects consistent with clinical phenotypes in ASD and individual syndromes. Mouse models can be clustered according to similarities and differences in upstream regulators, functional effects, and canonical pathways predicted by changes in their synaptic proteomes. Notably, we identified groups of models with shared changes in the synaptic proteome, which may indicate the presence molecular subtypes of ASD. Among diverse and shared molecular deficits observed in all ASD models tested, deficits in oxidative phosphorylation, and RhoA family signaling were predicted by bioinformatic analyses, and functional assays for Rac1 activity corroborate model-specific bioinformatic predictions. An appreciation of convergent synaptic changes underlying ASD is critical for defining pathogenic mechanisms and prioritizing treatments that can have broad efficacy (Sestan and State, 2018). Our results identify specific druggable pathways that may lead to the design of effective therapeutic interventions for diverse forms of ASD.
Materials and methods
Animal models and fractionation
Fmr1 (stock #003025, MGI:1857169) and Cntnap2 knockout (stock #017482, MGI: 2677631) mouse models and the BTBR+ (stock #002282), B6129SF2/J (stock #101045), and C57BL/6J (stock #000664) mouse strains were purchased from the Jackson Laboratory. Mice with Pten haploinsufficiency (MGI:2151804) and wild-type mice on the C57BL/6J background were obtained from the Page lab at the Scripps Research Institute Florida. Heterozygotes from the Anks1b conditional knockout line (stock #035048, MGI:5779292; Tindi et al., 2015) and wild-type mice from the Nestin-cre transgenic line (stock #003771, MGI:2176173) were bred in house (Carbonell et al., 2019). Mice from this line expressing the Nestin-cre transgene were genotyped for the Anks1bwt allele (forward: 5′-CACCCACAGCTCCATAGACAG-3′, reverse: 5′-GCACCTATTCCCTTCACCCTG-3′) and Anks1bfl allele (forward: 5′-AGTTGCCAGCCATCTGTTGT-3′, reverse: 5′-GGGTTCCGGATCAGCTTGAT-3′). Postsynaptic density (PSD) enriched fractions were isolated based on the method of Carlin et al. (1980) and Cohen et al. (1977) as we have done before (Jordan et al., 2004, 2007; Zhang et al., 2012; Tindi et al., 2015). Briefly, mice were euthanized in compliance with the Institutional Animal Care and Use Committee of the Albert Einstein College of medicine and both hippocampi were rapidly removed, placed in ice-cold solution A (0.32 M sucrose, 10 mM HEPES pH 7.4, 1 mM MgCl2, 0.5 mM CaCl2, protease and phosphatase inhibitors, and 0.1 mM PMSF), and Dounce homogenized. Samples were then centrifuged at 800g for 10 min to remove nuclei and other cellular debris, and the supernatant containing light membrane fractions was subjected to a second centrifugation at 30,000g for 25 min. The pellet was then resuspended in solution B (0.32 M sucrose, 10 mM HEPES pH 7.4) and layered on top of a 0.85, 1, and 1.2 M sucrose step gradient and centrifuged at 82,500g for 2 h. The synaptosomal fraction was collected at the 1–1.2 M sucrose interface and lysed using an equal volume of 1% Triton X-100 (in 0.32 M sucrose, 12 mM Tris pH 8.1) at 4°C for 15 min. Lysed synaptosomes were then centrifuged at 120,000g for 25 min to collect the PSD-enriched fraction.
Antibodies
SDS-PAGE and Western blot were performed under standard conditions using the LI-COR fluorescence-based system with the following antibodies: rabbit anti-Rac1 1:1,000 (Proteintech #24702-1-AP), rabbit anti-RhoA 1:1,000 (ABclonal), rabbit anti-Cdc42 (ABclonal), mouse anti-PSD95 1:1,000 (NeuroMab), and rat anti-tubulin 1:1,000 (Cell Signaling Tech). Statistical analysis for Western blot was performed in JMP 16 (SAS).
Tandem mass tag labeling of peptides and determination of labeling efficiency
PSD-enriched samples isolated from each mouse model were electrophoresed briefly (dye front about 5 mm beyond the gel well) using a 4–12% SDS-PAGE gel to remove SDS and other mass spectrometry (MS)-incompatible chemicals. Proteins were visualized by staining overnight with GelCode® Coomassie blue reagent (Pierce). A gel fragment containing all stacked protein bands was excised, reduced with DTT, alkylated with iodoacetamide, and digested using 5 ng/μL mass spectrometry-grade trypsin (Trypsin Gold, Promega). The resulting peptides were desalted using a Stage Tip manually packed with Empora C18 High Performance Extraction Disks (3 M) (Rappsilber et al., 2007) and eluted peptide solutions were dried under vacuum. Peptides were then resuspended in 18 μL acetonitrile (ACN), 57 μL of 0.2 M HEPES pH 8.5 was added to each sample and were reacted with unique isobaric labels within a 10-plex tandem mass tag set (TMT10-plex; Thompson et al., 2003). TMT10-plex amine reactive reagents (Thermo Fisher, 5 mg per vial) were re-suspended in 1024 μL anhydrous acetonitrile and 25 μL of reagent was added to each sample (TMT label: peptide [w/w] = 12:1). The mixture was incubated at RT for 1 h, quenched by the addition of 10 μL 5% hydroxylamine for 15 min, and acidified by the addition of 10 μL 10% formic acid. To calculate labeling efficiency, a 5-μL aliquot from each reaction was desalted on a StageTip, analyzed by LC–MS/MS with a Q-Exactive Orbitrap HF (high field), and the resulting spectra searched with MaxQuant using its corresponding TMT label as variable modifications on the N-terminus and lysine residues. Labeling efficiency was 95% or greater for all samples labeled (Supplementary Figure 5). To ensure that equal amounts of labeled peptides from each channel were mixed together, a two-step mixing strategy was employed: in the first step, an identical ~1 μL volume of peptides from each channel was mixed and analyzed, and the value of the median ratio (median of the ratios of all peptide intensities of one channel over their corresponding peptide average intensities of all channels) for each channel was determined as the correction factor. In the second step, the rest of the peptides were mixed by adjusting their volume using the correction factors. In this way, median ratios ranging from 0.97 to 1.02 were achieved as previously reported (Erdjument-Bromage et al., 2018). The final mixture of reaction products from 10 TMT channels were desalted on a Sep-Pak tC18 1 mL Vac Cartridge (Waters, #WAT03820). Eluted peptides were dried by vacuum centrifugation and stored at −20°C.
Hydrophilic interaction liquid chromatography fractionation of peptides
The TMT-labeled peptides from the PSD samples were then fractionated by offline HILIC to increase depth of coverage of synaptic proteins and to decrease ratio compression resulting from co-fragmenting peptides (Huang et al., 2017). To do this, the final TMT mixture was dissolved in 90% acetonitrile with 0.1% TFA and peptide separation was carried out on an Agilent pump equipped with a TSK gel amide-80 column (4.6 mm ID, 25 cm long) from TOSOH Bioscience, LLC, PA, United States. A gradient of 90% acetonitrile with 0.1% TFA was introduced over 65 min, and a fraction was collected every 2 min. Concatenated pools of peptides (10 pools) were finally created by pooling non-adjacent peptide fractions; about 10% of each pool was used for LC–MS/MS analysis.
Liquid chromatography–tandem mass spectrometry
Online chromatography was performed with a Thermo Easy nLC 1000 ultrahigh-pressure UPLC system (Thermo Fisher) coupled online to a Q-Exactive HF with a NanoFlex source (Thermo Fisher). Analytical columns (~23 cm long and 75 μm inner diameter) were packed in-house with ReproSil-Pur C18 AQ 3 μm reversed-phase resin (Dr. Maisch GmbH, Ammerbuch-Entringen). The analytical column was placed in a column heater (Sonation GmbH, Biberach) regulated to a temperature of 45°C. The TMT peptide mixture was loaded onto the analytical column with buffer A (0.1% formic acid) at a maximum backpressure of 300 bar. Peptides were eluted with a 2-step gradient of 3–40% buffer B (100% ACN and 0.1% formic acid) in 180 min and 40–90% B in 20 min, at a flow rate of 250 nL/min over 200 min using a 1D online LC–MS2 data-dependent analysis (DDA) method as follows: MS data were acquired using a data-dependent top-10 method, dynamically choosing the most abundant not-yet-sequenced precursor ions from the survey scans (300–1,750 Th). Peptide fragmentation was performed via higher energy collisional dissociation with a target value of 1 × 105 ions determined with predictive automatic gain control. Isolation of precursors was performed with a window of 1 Th. Survey scans were acquired at a resolution of 120,000 at m/z 200. Resolution for HCD spectra was set to 60,000 at m/z 200 with a maximum ion injection time of 128 ms. The normalized collision energy was 35. The underfill ratio specifying the minimum percentage of the target ion value likely to be reached at the maximum fill time was defined as 0.1%. Precursor ions with single, unassigned, or seven and higher charge states were excluded from fragmentation selection. Dynamic exclusion time was set at 30 s. Each of the TMT 10-plex samples was analyzed in triplicate.
All data were analyzed with the MaxQuant proteomics data analysis workflow (version 1.5.5.7) with the Andromeda search engine (Cox et al., 2011; Tyanova et al., 2016). The type of the group specific analysis was set to Reporter ion MS2 with 10plex TMT as isobaric labels for Q Exactive HF MS2 data. Reporter ion mass tolerance was set to 0.01 Da, with activated Precursor Intensity Fraction (PIF) value set at 0.75. False discovery rate (the rate at which identified positives are null) was set to 1% for protein, peptide spectrum match, and site decoy fraction levels. Peptides were required to have a minimum length of eight amino acids and a maximum mass of 4,600 Da. MaxQuant was used to score fragmentation scans for identification based on a search with an allowed mass deviation of the precursor ion of up to 4.5 ppm after time-dependent mass calibration. The allowed fragment mass deviation was 20 ppm. MS2 spectra were used by Andromeda within MaxQuant to search the UniProt mouse database (01092015, 16,699 entries) combined with 262 common contaminants. Enzyme specificity was set as C-terminal to arginine and lysine, and a maximum of two missed cleavages were allowed. Carbamidomethylation of cysteine was set as a fixed modification and N-terminal protein acetylation, deamidated (N, Q), and oxidation (M) as variable modifications. The reporter ion intensities were defined as intensities multiplied by injection time (to obtain the total signal) for each isobaric labeling channel summed over all MS/MS spectra matching to the protein group as previously validated (Tyanova et al., 2016). Following MaxQuant analysis, the protein and peptide.txt files were imported into Perseus (version 1.5.6.0) software which was used for statistical analysis of all the proteins identified.
Experimental samples and bioinformatic analyses
ASD models (Fmr1, Pten, Cntnap2, Anks1b, BTBR+) and appropriate controls (C57BL/6J, B6129SF2/J X Nestin-cre C57BL/6J) for each model were analyzed in biological triplicates, representing 24 independent samples. Each biological replicate was comprised enriched hippocampal PSD fractions pooled from 2 male mice of each genotype at 6–10 weeks of age (N = 3 biological replicates using a total of 6 mice for each genotype). As we could only label 10 samples at a time using the 10-plex TMT isobaric tags, the samples (Fmr1, Pten, Cntnap2, Anks1b, BTBR+ models, C57BL/6J controls, and Nestin-cre C57BL/6J controls; in triplicate) were analyzed in three independent LC–MS/MS runs, with each run containing the animal model and corresponding control.
Following the identification of proteins and TMT-dependent quantitation of their abundance in each sample, the fold-change was calculated by dividing the relative abundance of a protein in the ASD sample by its abundance in the appropriate control sample. Fold-changes were calculated only between samples within each independent MS analyses. For the Fmr1, Pten, Cntnap2, and BTBR+ models, we used the average of 2 samples of C57BL/6J mice as the control values. For the Anks1b Het mice on the Nestin-cre C57BL/6J background, we used wild-type littermates on the Nestin-cre C57BL/6J background as the control value. The fold-changes were plotted, and values were normalized to a value of 1. To increase the rigor of subsequent analyses, we filtered the protein sets in each experiment for proteins identified and quantified by at least 3 unique peptides, and with a Q-value of 0, representing the highest significance for intensity and score in MaxQuant. We used this list of protein as input in StringDB under the species Mus musculus for gene ontology (GO) enrichments analyses (Szklarczyk et al., 2015).
Functional annotation, activation prediction, and regulatory network construction were performed for each mouse model using Ingenuity Pathway Analysis (IPA, QIAGEN Bioinformatics, Fall Release 2022). Fold-change values for each protein in the PSD proteome were used as input for a Core Analysis in IPA using the following default settings: Expression Analysis; General Settings = Ingenuity Knowledge Base (Genes Only), Direct Relationships; Networks = Interaction networks, include endogenous chemicals; Node Types = All; Data Sources = All; Confidence = Experimentally Observed; Species = All; Tissues & Cell Lines = All; Mutations = All. To score similarity between models overall and in each domain (Upstream Regulators, Downstream Effects, and Canonical Pathways), Analysis Match was used. In Analysis Match, activity signatures for each Core Analysis are generated by taking the top 50 (Upstream Regulators and Downstream Effects) or 10 (Canonical Pathways) entities that are activated (z-score > 2) and inhibited (z-score < −2) in each domain. Similarity z-scores between Core Analyses are defined using the following formula, where N is the total number of overlapping entities in each analysis, N+ is the number of correct matches, and N− the number of incorrect matches:
The raw z-score is divided by a hypothetical perfect match (N+ = N <= 100, N− = 0) and multiplied by 100% to obtain a normalized z-score, where a score of 100 is a perfect match to itself. Raw and normalized z-scores are negative when N− > N+, yielding an opposite activation signature. To generate overall value of p scores, the −log10 of the p-values (maximum value of 50) were calculated for each domain and expressed as a percentage of the maximum possible −log10p (a perfect match to itself). To compare models, the Core Analysis of each model were used as input for a Comparison Analysis in IPA. Hierarchical analysis was used to cluster annotations in Upstream Analysis, Downstream Effects, and Canonical Pathways, and to cluster models for similarity relationships.
G-LISA
G-protein linked immunoassays (G-LISA) were used to quantify active Rac1 in hippocampal tissue lysates. Active Rac1 substrate (the G-protein binding domain of Pak1) fused to GST protein (GST-Pak1) was generated in bacteria and purified using glutathione beads. This substrate was added at a concentration of 1 μg to 50 mM Ammonium Carbonate and absorbed onto high-binding 96-well plates overnight to create binding surfaces for GTP-bound Rac1. Hippocampal tissues from the five different mouse models were prepared using the method described by the Cytoskeleton Inc. BK-128 Kit. (N = 6 each of BTBR+, Fmr1, and Cntnap2 mice, N = 5 Pten mice, N = 9 Anks1b mice, and N = 9 each of the two control Mice, C57BL/6J and B6129SF2/J X Nestin-cre C57BL/6J mice were used, all prepared from 6 to 10 weeks of age). Briefly, tissues were lysed for 1 min in the kit lysis buffer, and lysates were snap-frozen to preserve GTP-bound active Rac1. All steps were conducted on ice. Lysates were then incubated on the GST-Pak1 coated 96-wells for exactly 30 min on ice. Primary and secondary antibodies as well as HRP-substrate were provided by the BK-128 Kit, and active Rac1 was detected by standard ELISA using immunodetection at 490 OD and HRP-substrate colorimetric measurements.
Results
TMT-MS analyses of postsynaptic density fractions from mouse models of autism
To compare the synaptic proteomes from different ASD mouse models, we obtained Fmr1 knockout (Fmr1 KO; Consortium, T.D.-B.F.X, 1994), Pten haploinsufficiency (Pten Het) (Clipperton-Allen and Page, 2014; Clipperton-Allen and Page, 2015), Cntnap2 knockout (Cntnap2 KO; Peñagarikano et al., 2011), Anks1b haploinsufficiency (Anks1b Het; Carbonell et al., 2019), and the BTBR+ inbred mouse strains (McFarlane et al., 2008). All selected models display abnormal behaviors in domains relevant to ASD, including social approach and interaction, stereotyped movements, learning and memory, and sensorimotor function. These models have also been used to illustrate synaptic dysfunction, including altered synapse formation, excitatory/inhibitory balance, synaptic plasticity, and glutamatergic signaling (Hulbert and Jiang, 2016; Varghese et al., 2017; Bagni and Zukin, 2019; Carbonell et al., 2019). We analyzed changes in the synaptic proteomes of hippocampal postsynaptic density (PSD)-enriched fractions isolated from each of the five mouse models of ASD. PSD-enriched fractions were collected as we have done before (Jordan et al., 2004, 2007; Zhang et al., 2012; Tindi et al., 2015; see Methods) and modified to a single Triton X-100 extraction step (Tindi et al., 2015) as opposed to two. This allowed us to collect enough material from mouse hippocampi at the expense of purity. However, Western blots show that isolated fractions were highly enriched for the post synaptic marker PSD95 (Figure 1A). Moreover, the proteins identified in the three independent experiments (described below) show ~50–80% overlap with the Genes to Cognition database of PSD components (G2Cdb; Croning et al., 2009) and ~30–50% overlap with the Synaptic Gene Ontologies curated synaptic database (SynGO) (Koopmans et al., 2019; Supplementary Table 1). Proteins in each sample were then labeled with a unique isobaric tag within a 10-plex Tandem Mass Tag set (TMT). TMTs are reactive isobaric tags with unique permutations of 14C and 15N stable isotopes that imbue peptides with small differences in mass that can then be resolved by mass spectrometry (MS) (Figure 1B). This approach allows for up to 10 samples to be mixed and analyzed simultaneously by MS, which reduces variability, enables robust quantitation of identified proteins, and allows for rigorous quantitative and comparative analyses between samples. Overall, this method allowed us to overcome poor quantitation traditionally associated with MS-based proteomics. The proteins identified in each run, along with number of unique peptides, scores, and other MS-based information are listed in Supplementary Table 1. Using all genes as a search space, gene ontology (GO) analyses using STRING (Szklarczyk et al., 2019) revealed that isolated fractions are overwhelmingly enriched for synaptic terms (Figure 1C). GO analyses also revealed a similar distribution of Cellular Compartments and Reactome Pathways (Figure 1D) across each of the three independent MS experiments, and that each of the specific pathways was enriched to a similar extent. We used SynGO (Koopmans et al., 2019) to restrict the search space to brain-specific genes and found that Synapse, Process in the synapse, and Postsynapse were the most significant GO terms (Figure 1E). Hierarchical plots for GO terms show an overwhelming enrichment (high Q-scores) for synaptic processes for all three experiments (Figure 1F). Overall, these results show that all experimental samples were highly enriched in synaptic components and were similar to each other, which increases the strength of subsequent comparisons.
Figure 1. Methodology for isobaric-tag based quantitative comparison of synaptic fractions of ASD mouse models. (A) Synapse-enriched fractions from sample mouse models of autism spectrum disorder (ASD) demonstrate qualitative enrichment for the synaptic marker PSD95 by Western blot (4 μg each sample) compared to total brain fractions (T). (B) Isobaric labeling and comparative proteomics schematic showing the ASD mouse models used (De Rubeis et al., 2014; Gaugler et al., 2014; Sztainberg and Zoghbi, 2016; Weiner et al., 2017; Ruzzo et al., 2019) and appropriate control samples (Pinto et al., 2014; Hulbert and Jiang, 2016; Parikshak et al., 2016; Varghese et al., 2017; Verma et al., 2019) modified from the schematic from the Thermo Scientific Catalog for the TMT10plex™ Isobaric Kit. Samples were reacted with unique isobaric tags, whose identifier (e.g., 128 N) corresponds to the molecular weight (Daltons) and stable isotope amino acid (N or C) present in that tag. (Right) Idealized MS output showing the same protein in each of the 10 samples identified as 10 contiguous peaks. This permits immediate and rigorous relative quantitation across all samples. m/z represents mass/charge ratio of peptides. N = 6 total animals for each genotype. (C) Bar graph representing the top cellular components sorted by statistical significance (False Discovery Rate), identified by gene ontology (GO) analysis of proteins in a sample MS run (Exp 1) and using all genes as a search space. (D) Same as figure C but for showing top Cellular Compartment and Reactome Pathways for each MS run. (E) Enrichment analyses using SynGO and using only brain-specific genes as a search space. (F) Hierarchical plots showing the top GO terms enriched in each experiment and associated q-values (p-values adjusted for multiple comparisons) for each term.
Broad comparison of synaptic proteomes across ASD mouse models shows evidence of molecular subtypes
To compare findings across ASD models, we first characterized each model by calculating the ratio of synaptic protein abundance in the model to its appropriate control (fold-change). C57BL/6J (duplicates) mice were used as controls for all models except for the Anks1b mice, which used wildtype B6129SF2/J X Nestin-cre C57BL/6J controls from the same colony. To preserve the quantitative rigor afforded by the isobaric tags, fold-changes were only calculated within experiments (between samples in each independent MS analysis). Moreover, we used a highly stringent filter to increase the significance of findings, considering only proteins identified and quantified by at least 3 unique peptides, and whose MaxQuant calculated Q-values were 0, representing the lowest possible false discovery rate. Finally, only proteins that were identified and quantified in at least 2 independent experiments were considered for subsequent analyses (Supplementary Table 2; 97.6 and 91.2% of proteins found and quantified in Experiment 1 were found and quantified in Experiments 2 and 3, respectively). Distribution plots of protein fold-changes revealed that most proteins were unchanged at hippocampal synapses in the ASD mouse models (fold-changes were normally distributed across a value of 1; Figure 2A). This is expected as large changes in the synaptic proteome would likely be lethal. To measure the overall similarities across models, we performed multivariate analyses using PCA (Principal Components Analyses; Figure 2B) and calculated correlations using protein identities and fold-changes (Figure 2C). Results show that Fmr1 and Cntnap2 models are most similar, with Anks1b and Pten models appearing to comprise a second molecular subgroup. The inbred BTBR+ strain was most dissimilar.
Figure 2. Similarity analyses between mouse models of ASD. (A) Distribution histograms for all proteins identified each animal model. Results show that fold-values show distribution with outliers. (B) Principal components analysis for the five ASD models plotted across (left) one and (right) two dimensions reveal broad similarities between Cntnap2 and Fmr1 models, as well as between Anks1b and Pten models. (C) Correlation matrices (left- data points, right- Pearson’s correlation) based on genes and fold-changes. For all panels, proteins identified = ~1,300–1,500 from n = 3 MS runs of n = 6 mice per genotype.
Dysregulated molecular pathways for oxidative metabolism and Rho GTPase signaling are shared across ASD mouse models
To analyze the functional implications of the altered synaptic proteomes observed, we used Ingenuity Pathway Analysis (IPA, QIAGEN Bioinformatics). IPA derives functional insights and pathway nomenclatures from a continuously updated knowledge base drawn from the literature (Krämer et al., 2014). We first performed a Comparison Analyses in IPA. Sorted by overall z-scores (standard deviation away from expected values), we found that several canonical pathways were altered across all models in a highly significant manner (Figure 3A; Supplementary Table 3). Moreover, several pathways were predicted to be regulated in the same way (either activated or inhibited) across all ASD models. We hypothesize that these canonical pathways represent critical or potentially causative molecular mechanisms in autism that could serve as novel targets for therapeutic intervention. All tested models show strongly activated (positive z-scores) oxidative phosphorylation and strongly downregulated Granzyme A, Sirtuin, RhoGDI, and HIPPO signaling pathways. These results corroborate long-standing research linking neurodevelopmental disorders to mitochondrial dysfunction (oxidative phosphorylation; Giulivi et al., 2010; Rossignol and Frye, 2012) immune regulatory systems (Granzyme A, Sirtuin; Careaga et al., 2017; Meltzer and Van de Water, 2017; Hughes et al., 2018), and cell proliferation and apoptosis (Granzyme A, HIPPO, RhoGDI, Sirtuin; Dong et al., 2018; Courchesne et al., 2019). In addition, our analyses identified specific dysregulated molecules that may represent the core mechanisms linked to deficits within these broad cellular pathways (Supplementary Table 4). Sirtuins, like granzyme signaling, have been implicated in neuroprotection, cognition, and metabolism, as well as Alzheimer’s disease and aging (Bonda et al., 2011). Other pathways that were altered in all five models showed multi-directional dysregulation. As expected, based on prevailing literature and samples tested, synaptogenesis signaling and glutamate receptor pathways were also significantly affected in all ASD models, corroborating recent hypotheses classifying ASDs as synaptopathies. We also found a significant number of canonical pathways linked to the activity of the small G-protein family RhoA. Interestingly, these pathways showed a similar pattern where the Anks1b model was downregulated while the other four models were upregulated.
Figure 3. Distinct pattern of similarities and differences in Canonical Pathways among mouse models of autism. (A) Heat maps showing the predicted activation or inhibition (see LUT legend) of the top 30 significant (z-scores) canonical pathways shared in all ASD models. (B) Twenty datasets with randomized and normally distributed fold-ratio values were analyzed by IPA. Results reveal top enriched pathways as expected, but with significantly lower z-scores, and no pathways show changes in the same direction (z-score) for all ASD models. (lower) Plots of z-scores (Y-axis) for real datasets (green diamonds) for the identified canonical pathways (X-axis). Overlaid are the distribution of z-scores for the randomized datasets (pink circles; mean ± SD). A significant number of outliers were identified in the real datasets (green diamonds outside the light blue area). (C) Heat map showing the outlier canonical pathways (>3 standard deviations away from the mean of randomized dataset in at least 2 animal models). (D) Analysis Match in IPA to determine the similarities between models based on shared canonical pathways, upstream and downstream regulators. Numbers represent correlation values with 100% representing perfect overlap and −100% being anti-correlated. Results track the broad similarities seen in Figure 2.
While IPA incorporates p-values, z-scores, and other analyses to address statistical biases based on population subsets (Krämer et al., 2014), the exact methodology is proprietary. To independently test the significance of our findings, we created 20 randomized datasets containing the same proteins as in our real samples, but with randomly assigned fold-ratios. These datasets were generated to reflect the distribution and standard deviations calculated from the real datasets of each animal model. IPA analyses of randomized datasets revealed “significant” canonical pathways, but with z-score ranges that were substantially lower than those calculated for the actual sample datasets (average random z-score range = ± ~ 2.5; average real z-score range = −5.01 to 6.48; Figure 3B). Plotting the z-scores for the actual dataset together with the mean and standard deviation of the 20 randomized datasets reveals outlier canonical pathways that were at least three standard deviations removed from the mean (Figure 3B). The standard deviations and z-scores calculated for both actual and randomized data were almost identical, suggesting a standardized method for calculating expectation values. Pathways that were greater than three standard deviations removed from the mean include Oxidative phosphorylation, Sirtuin signaling, Granzyme A signaling, Synaptogenesis, and Rho signaling pathways, among others (Figure 3C).
We then used Analysis Match to create a similarity matrix for all ASD models tested based on Upstream Regulators, Downstream Effects (Diseases and Functions), Canonical Pathways from each Core Analysis, and overall z-scores (Figure 3D; Supplementary Figures 1-3). The top pathways in each category and for each animal model were based on p-values, z-scores, and activation values (Krämer et al., 2014; Supplementary Table 3). Similarity z-scores reflect the percentage of maximum similarity to a given model (out of 100%, a perfect match to itself). Analysis Match based on functional annotations reveals significant similarities across ASD models, which corroborates results presented in Figure 2 and shows that BTBR+ mice represent a distinct molecular subgroup.
Altered expression and function of Rho GTPases in ASD mouse models validates bioinformatic analyses of synaptic proteomes
To validate our proteomic and bioinformatic analyses, we performed biochemical analyses for one altered cellular pathway among all identified. Across all mouse models tested, RhoA family G-protein signaling pathways are involved in all top canonical pathway hits (Figure 3A) including Synaptogenesis Signaling Pathway, Signaling by Rho family GTPases, Rac Signaling, RhoGDI signaling, and Regulation of Actin-based motility by Rho. Bioinformatic predictions validate previous studies as Rho GTPases have crucial roles in synaptic function and have been investigated extensively in both human-based experiments and mouse models of ASD (Guo et al., 2020). Directionality of changes was largely consistent among mouse models, predicting activated signaling by Rho GTPases and inhibited regulation by RhoGDI in Pten, Cntnap2, BTBR+, and Fmr1 mice, but with opposite changes observed in the Anks1b mice (Figure 3A). Depictions of the Rac1 signaling pathways altered in the ASD models show that receptor tyrosine kinases and other synaptic pathways primarily drive changes in Rac activity (Supplementary Figure 6). Components of the WAVE complex and other regulators of actin nucleation and polymerization feature prominently among altered molecules in this pathway. Western blots for Rho family G-proteins Rac1, RhoA, and Cdc42 in synaptic fractions detected Rac1, but not RhoA or Cdc42 in synaptic fractions (Figure 4A). This corroborates the MS analyses showing that Rac1 was identified with more unique peptides in enriched PSD fractions compared to RhoA and CDC42 (Supplementary Table 1). Rac1 expression was increased in the Fmr1 KO model in synaptic fractions, but largely unchanged across all models (Figures 4A,B). Western blots of whole brain lysates reveal significant changes in several Rho family GTPases across models (Figure 4C). For a more informative test measuring function, we performed G-protein activation ELISAs (G-LISAs) to directly measure Rac1 activity in total hippocampal lysates. As predicted by bioinformatic analyses, Rac1 activity in Anks1b model mice was substantially downregulated compared to controls, as well as significantly reduced compared to Cntnap2, Pten, and Fmr1 models (Figure 4D). While Rac1 activity in Pten, Fmr1, Cntnap2, and BTBR+ mice was not significantly upregulated as predicted, a difference in the variance of the WT and ASD dataset was striking. We used F-test calculations for unequal variances to show the differences were statistically significant (Figure 4E).
Figure 4. Validation of RhoA family GTPases dysfunctions in ASD mouse models (A) (Left) Rac1 was the only Rho GTPase detected in whole lysate (WL) and synaptic (PSD enriched) fractions (5 μg each). (Right) Western blots (5 μg PSD enriched fraction) showing Rac1 and synaptic marker PSD95 in the different models. (B) Quantitation shows increased Rac1 in synaptic fractions from Fmr1 KO, but not from other mouse models, compared to controls. (C) Quantitation of RhoA family GTPases in total lysates shows differences throughout. Expression of Cdc42 in total lysate was reduced for all models except in Anks1b Het, where no significant change was observed (30 μg whole lysate). (D) G-LISA showing difference in fold-change in active Rac1 compared to appropriate controls. Only Anks1b showed significant reduction in active Rac1. (E) F-test calculations for unequal variances show significance. Bar graphs show mean ± SEM, n = 3–10 samples for each mouse model and wild-type control strain, Student’s t-test, *p < 0.05 **p < 0.01 ***p < 0.001. For B,C sample values are compared to controls.
Discussion
ASDs are among the most heritable neuropsychiatric conditions, with early epidemiological studies estimating concordance rates of ~90% in monozygotic twins, compared to ~10% for dizygotic twins (Steffenburg et al., 1989; Bailey et al., 1995; Hallmayer et al., 2011). Hundreds of genes and chromosomal loci for ASD susceptibility have been reported, suggesting a highly heterogeneous genetic architecture (Chen et al., 2015; de la Torre-Ubieta et al., 2016). While these studies highlight the notion that there is no single “autism gene,” functional and bioinformatic analyses of genetic studies have identified convergent cellular pathways, including those regulating transcription (Sanders et al., 2012; de Rubeis et al., 2013), excitatory/inhibitory (E/I) balance (Blatt et al., 2001; Hussman, 2001; Rubenstein and Merzenich, 2003; Gao and Penzes, 2015; Nelson and Valakh, 2015), and especially synaptic function (Bucan et al., 2009; Glessner et al., 2009; Hussman et al., 2011; Buxbaum et al., 2012; Sanders et al., 2012; Yu et al., 2013; Brett et al., 2014; Buxbaum et al., 2014; Cukier et al., 2014; Iossifov et al., 2014; Pinto et al., 2014; Ronemus et al., 2014; Toma et al., 2014; Yuen et al., 2015). Autism gene databases show clear enrichment for synaptic pathways, and syndromic ASD mouse models are primarily associated with mutations in synaptic proteins or their regulators, such as SHANKs, neurexins, neuroligins, SynGAP1, FMRP, and TSC1/2 (Pinto et al., 2014; Yoo, 2015; Ramaswami and Geschwind, 2018). Indeed, most (if not all) identified ASD mouse models such as those for Fragile X (Bhakar et al., 2012), tuberous sclerosis (Auerbach et al., 2011), Angelman (Mabb et al., 2011) Phelan McDermid (Sarowar and Grabrucker, 2016), and ANKS1B syndromes show deficits in diverse forms of synaptic function as well as altered neuronal and synaptic morphology (Tindi et al., 2015; Carbonell et al., 2019). Here, we compared the hippocampal postsynaptic proteomes from five mouse models for autism that display synaptic dysfunction, including altered synapse formation, excitatory/inhibitory balance, synaptic plasticity, and glutamatergic signaling (Hulbert and Jiang, 2016; Varghese et al., 2017; Bagni and Zukin, 2019; Carbonell et al., 2019). The goal of this work was to test the hypothesis that common structural and functional deficits among diverse ASD mouse models are caused by similar underlying synaptic deficits and ASD-linked susceptibility factors that ultimately converge on common signaling pathways. By leveraging widely studied ASD mouse models, we hope to identify high-value therapeutic targets for autism.
Wide-spread genetic- and transcriptomic-based screening and comparative studies of disease models suffer from important caveats that limit their interpretability. Often, little is known about how disease-linked genes or other chromosomal loci ultimately correlate to protein abundance and function. Moreover, multiple studies confirm the weak to nil correlation between transcript and protein abundance (Greenbaum et al., 2003; Maier et al., 2009). Transcriptomes display 100-fold ranges in translation efficiency (Ingolia et al., 2009), proteomes reveal >1,000-fold ranges in half-lives (Doherty et al., 2009), and coupled transcriptomic and proteomic analyses reveal that proteins are ~900 times more abundant than corresponding mRNAs, but with ratios that span over five orders of magnitude (Schwanhäusser et al., 2011). Finally, transcriptomic methods are unlikely to provide meaningful information for ASDs such as Angelman syndrome, Fragile X syndrome, and tuberous sclerosis, where the primary contributors to disease etiology are thought to be altered regulation of protein translation and degradation. These important caveats are often ignored in discovery-based genomic research (Wang et al., 2017; Hoogendijk et al., 2019). To overcome these concerns, we have taken a quantitative proteomic approach based on 10-plex Tandem Mass Tags (TMT; Erdjument-Bromage et al., 2018). By directly measuring synaptic protein abundance, we bypass the concerns about poor transcript and protein correlation, as well as the need to understand the complex functional consequences of altered gene expression.
We have used IPA (Qiagen) to analyze the proteomic data and generate inferences on pathways altered in each disease model. IPA derives functional insights and pathway nomenclatures from a continuously updated knowledge base drawn from the literature (Krämer et al., 2014). While this implies a broad context for functional annotations, it also means that proteins are likely to be associated with scientific trends. Perhaps not surprisingly given the current pandemic, several Diseases and Functions associated with the ASD models are related to viral infections (Supplementary Figure 3). Moreover, there is the possibility for missed hits given the lack of comprehensive annotations for cellular and molecular pathways. One example is the surprising finding that the most downregulated proteins across four of the five ASD models were associated with oligodendrocyte function (Supplementary Table 2). The appearance of proteins related to myelin was surprising in our synaptic fractions but may reflect the presence of neuron-oligodendrocyte synapses that may co-purify using our methods. Despite these limitations, IPA lists the specific proteins and their fold-changes leading to the assigned pathways. For example, the prediction for activated Oxidative phosphorylation across ASD models is primarily driven by increases in the abundance of ATPases and different subunits of the cytochrome c oxidases and the NADH:ubiquinone oxidoreductase. Changes in synaptogenesis signaling pathways were driven by both increases and decreases in the abundance of regulators of the actin cytoskeleton which include several small G-proteins (Supplementary Table 4).
Regulatory networks predicted by synaptic changes in each model reveal both novel results and patterns consistent with the literature (Figures 3, 4; Supplementary Figures 1-4; Supplementary Tables 1-4). In the Fmr1 KO model, fold-changes reveal that most altered proteins were upregulated as seen by the rightward skew of the ratio distributions (Figure 2A), which is consistent with the primary role of FMRP as a translational repressor (Darnell et al., 2011). In the model of constitutive Pten haploinsufficiency, synaptic dysfunction is likely due to neuronal Pten, since altered synaptic composition was also observed in a neuron-specific conditional knockout mouse that displayed changes in activity patterns and repetitive behavior (Lugo et al., 2014). Importantly, we directly compared synaptic proteomes in parallel, a necessary approach for finding shared phenotypes and broadly applicable therapeutic targets (Sestan and State, 2018; Silverman and Ellegood, 2018). Drawn from internally controlled measurements, our results demonstrate quantifiable differences in synaptic composition that support the presence of molecular subtypes in ASD rather than universal changes in the same direction.
In the five autism mouse models chosen for this study, changes in the synaptic proteome highlighted Rho family small GTPase signaling as a commonly altered cellular pathway (Figures 3, 4). These results complement a growing body of evidence showing that altered activation of Rho GTPases is an important mechanism of disease in autism and other neurodevelopmental disorders (Pinto et al., 2010; Zeidán-Chuliá et al., 2013; Guo et al., 2020). Functional validation results showed significantly reduced active Rac1 in Anks1b mice compared to wildtype, and compared to Pten, Fmr1, and Cntnap2 mice, which was consistent with proteomic analysis (Figures 3, 4). These results suggest that downregulation of Rac1 pathways may inhibit the processes that require GTP-bound Rac1, thereby potentially causing some of the manifestations of ANKS1B syndrome, including ASD. Consistent with our findings that synaptic proteins in the Rac pathway (Supplementary Figure 6) and Rac1 itself is differentially regulated in ASD mouse models (Figure 4), further evidence continues to emerge that Rac1 signaling plays an important role in the pathobiology of ASD and other neurodevelopmental disorders. Rac1 mutations associated with intellectual disability impair synaptic plasticity (Tian et al., 2018), and variants in the RhoGEF TRIO can cause ASD, intellectual disability, schizophrenia, or macrocephaly, with severe phenotypes predicted by Rac1 overactivation (Katrancha et al., 2017; Sadybekov et al., 2017; Barbosa et al., 2020). Both up- and down-regulation of Rac1 in mice causes social deficits, an important symptom of ASD (Ma et al., 2022). Downstream of Rac1 in actin polymerization, mutations in PAKs and WASF1 can cause macrocephaly, seizures, intellectual disability, or ASD (Ito et al., 2018; Horn et al., 2019). Overall, Rac1 pathways are commonly dysregulated in ASD and other neurodevelopmental disorders and are therefore attractive targets for pharmacotherapy (Zamboni et al., 2018; Guo et al., 2020).
Caveats to our work include comparisons at a single developmental timepoint and single brain region. Models of ASD show behavioral and biological phenotypes that differ based on the spatiotemporal targeting of the gene (Del Pino et al., 2018). A proteomic study of cortical synapses noted smaller differences in Fmr1 KO mice after 3 weeks of age (Tang et al., 2015), and both Shank3 and Syngap1 show changes in interaction partners throughout development (Li et al., 2016, 2017). Moreover, we focused on hippocampal synapses, but there are important roles for the neocortex, striatum, and cerebellum in the pathophysiology of ASD animal models (Golden et al., 2018). Fractionation of the PSD limits our findings to synaptic profiles in glutamatergic neurons, but ASD risk genes also influence interneuron development (Cntnap2; Peñagarikano et al., 2011) and oligodendrocyte maturation (Pten; Lee et al., 2019). The number of isobaric tags available in the TMT set limited the number of simultaneous comparisons possible by mass spectrometry. We therefore had to analyze the biological replicates as three independent experiments, which increases variability. To minimize this variability, fold-ratios were calculated only between samples in the same MS run. Future experiments using increased isobaric tags such as the new TMT16-plex system could be used to compare proteomes across cell types and developmental stages, and additional ASD models (Stessman et al., 2016; Sestan and State, 2018; Iakoucheva et al., 2019).
Understanding the molecular mechanisms underlying ASDs is critical for the development of novel therapies. Our proteomic approach overcomes critical confounds associated with gene-based studies of ASD etiology that improperly equate transcript levels, epigenetic modifications, single-nucleotide polymorphisms (SNPs), or copy number variations (CNVs) to changes in protein abundance. This work shows that synaptic proteomes, as identifiable and quantifiable phenotypes of diverse ASD models, can lead to the identification of molecular convergence. We propose that this can lead to important insights into ASD etiology and yield high-value targets to pursue for future therapies.
Data availability statement
The raw mass spectrometry data generated during this study are available at MassIVE (UCSD), https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp under the deposition number MSV000091848.
Author contributions
AC, CF-C, ID, and BJ designed and performed postsynaptic fractionation in autism mouse models. HE-B and TN performed tandem-mass-tag mass spectrometry and protein identification. AC and ID performed bioinformatics analysis and Western blots. ID and SD performed the GLISA assays. AC-A and DP developed and characterized Pten haploinsufficiency mouse model. AC, CF-C, ID, SD, and BJ interpreted results and wrote the paper. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by NIH R01AG039521, NIH R01NS118820, and NIH R56MH115201 to BJ, T32GM007288 to AC, NIH S10RR027990 to TN, and NIH R01MH108519 to DP. Significant support for this work came from the Rose F. Kennedy Intellectual and Developmental Disabilities Research Center (IDDRC), which is funded through the center grant NIH U54HD090260.
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.
Publisher’s note
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2023.1152562/full#supplementary-material
References
Abraham, J. R., Szoko, N., Barnard, J., Rubin, R. A., Schlatzer, D., Lundberg, K., et al. (2019). Proteomic investigations of autism brain identify known and novel Pathogenetic processes. Sci. Rep. 9:13118. doi: 10.1038/s41598-019-49533-y
Andoh, M., Shibata, K., Okamoto, K., Onodera, J., Morishita, K., Miura, Y., et al. (2019). Exercise reverses behavioral and synaptic abnormalities after maternal inflammation. Cell Rep. 27, 2817–2825.e5. doi: 10.1016/j.celrep.2019.05.015
Auerbach, B. D., Osterweil, E. K., and Bear, M. F. (2011). Mutations causing syndromic autism define an axis of synaptic pathophysiology. Nature 480, 63–68. doi: 10.1038/nature10658
Bagni, C., and Zukin, R. S. (2019). A synaptic perspective of fragile X syndrome and autism Spectrum disorders. Neuron 101, 1070–1088. doi: 10.1016/j.neuron.2019.02.041
Bailey, A., le Couteur, A., Gottesman, I., Bolton, P., Simonoff, E., Yuzda, E., et al. (1995). Autism as a strongly genetic disorder: evidence from a British twin study. Psychol. Med. 25, 63–77. doi: 10.1017/S0033291700028099
Barbosa, S., Greville-Heygate, S., Bonnet, M., Godwin, A., Fagotto-Kaufmann, C., Kajava, A. V., et al. (2020). Opposite modulation of RAC1 by mutations in TRIO is associated with distinct, domain-specific neurodevelopmental disorders. Am. J. Hum. Genet. 106, 338–355. doi: 10.1016/j.ajhg.2020.01.018
Barnes, S. A., Wijetunge, L. S., Jackson, A. D., Katsanevaki, D., Osterweil, E. K., Komiyama, N. H., et al. (2015). Convergence of hippocampal pathophysiology in Syngap+/− and Fmr1−/y mice. J. Neurosci. 35, 15073–15081. doi: 10.1523/JNEUROSCI.1087-15.2015
Bhakar, A. L., Dolen, G., and Bear, M. F. (2012). The pathophysiology of fragile X (and what it teaches us about synapses). Annu. Rev. Neurosci. 35, 417–443. doi: 10.1146/annurev-neuro-060909-153138
Blatt, G. J., Fitzgerald, C. M., Guptill, J. T., Booker, A. B., Kemper, T. L., and Bauman, M. L. (2001). Density and distribution of hippocampal neurotransmitter receptors in autism: an autoradiographic study. J. Autism Dev. Disord. 31, 537–543. doi: 10.1023/A:1013238809666
Bonda, D. J., Lee, H. G., Camins, A., Pallàs, M., Casadesus, G., Smith, M. A., et al. (2011). The sirtuin pathway in ageing and Alzheimer disease: mechanistic and therapeutic considerations. Lancet Neurol. 10, 275–279. doi: 10.1016/S1474-4422(11)70013-8
Bourgeron, T. (2015). From the genetic architecture to synaptic plasticity in autism spectrum disorder. Nat. Rev. Neurosci. 16, 551–563. doi: 10.1038/nrn3992
Brett, M., McPherson, J., Zang, Z. J., Lai, A., Tan, E. S., Ng, I., et al. (2014). Massively parallel sequencing of patients with intellectual disability, congenital anomalies and/or autism spectrum disorders with a targeted gene panel. PLoS One 9:e93409. doi: 10.1371/journal.pone.0093409
Brown, E. A., Lautz, J. D., Davis, T. R., Gniffke, E. P., VanSchoiack, A. A. W., Neier, S. C., et al. (2018). Clustering the autisms using glutamate synapse protein interaction networks from cortical and hippocampal tissue of seven mouse models. Mol. Autism. 9:48. doi: 10.1186/s13229-018-0229-1
Bucan, M., Abrahams, B. S., Wang, K., Glessner, J. T., Herman, E. I., Sonnenblick, L. I., et al. (2009). Genome-wide analyses of exonic copy number variants in a family-based study point to novel autism susceptibility genes. PLoS Genet. 5:e1000536. doi: 10.1371/journal.pgen.1000536
Buxbaum, J. D., Bolshakova, N., Brownfeld, J. M., Anney, R. J. L., Bender, P., Bernier, R., et al. (2014). The autism simplex collection: an international, expertly phenotyped autism sample for genetic and phenotypic analyses. Mol. Autism. 5:34. doi: 10.1186/2040-2392-5-34
Buxbaum, J. D., Daly, M. J., Devlin, B., Lehner, T., Roeder, K., State, M. W., et al. (2012). The autism sequencing consortium: large-scale, high-throughput sequencing in autism spectrum disorders. Neuron 76, 1052–1056. doi: 10.1016/j.neuron.2012.12.008
Carbonell, A. U., Cho, C. H., Tindi, J. O., Counts, P. A., Bates, J. C., Erdjument-Bromage, H., et al. (2019). Haploinsufficiency in the ANKS1B gene encoding AIDA-1 leads to a neurodevelopmental syndrome. Nat. Commun. 10:3529. doi: 10.1038/s41467-019-11437-w
Careaga, M., Rogers, S., Hansen, R. L., Amaral, D. G., van de Water, J., and Ashwood, P. (2017). Immune Endophenotypes in children with autism Spectrum disorder. Biol. Psychiatry 81, 434–441. doi: 10.1016/j.biopsych.2015.08.036
Carlin, R. K., Grab, D. J., Cohen, R. S., and Siekevitz, P. (1980). Isolation and characterization of postsynaptic densities from various brain regions: enrichment of different types of postsynaptic densities. J. Cell Biol. 86, 831–845. doi: 10.1083/jcb.86.3.831
Cellot, G., Maggi, L., di Castro, M. A., Catalano, M., Migliore, R., Migliore, M., et al. (2016). Premature changes in neuronal excitability account for hippocampal network impairment and autistic-like behavior in neonatal BTBR T+tf/J mice. Sci. Rep. 6:31696. doi: 10.1038/srep31696
Chen, J. A., Peñagarikano, O., Belgard, T. G., Swarup, V., and Geschwind, D. H. (2015). The emerging picture of autism spectrum disorder: genetics and pathology. Annu. Rev. Pathol. 10, 111–144. doi: 10.1146/annurev-pathol-012414-040405
Clipperton-Allen, A. E., and Page, D. T. (2014). Pten haploinsufficient mice show broad brain overgrowth but selective impairments in autism-relevant behavioral tests. Hum. Mol. Genet. 23, 3490–3505. doi: 10.1093/hmg/ddu057
Clipperton-Allen, A. E., and Page, D. T. (2015). Decreased aggression and increased repetitive behavior in Pten haploinsufficient mice. Genes Brain Behav. 14, 145–157. doi: 10.1111/gbb.12192
Cohen, R. S., Blomberg, F., Berzins, K., and Siekevitz, P. (1977). The structure of postsynaptic densities isolated from dog cerebral cortex. I. Overall morphology and protein composition. J. Cell Biol. 74, 181–203. doi: 10.1083/jcb.74.1.181
Consortium, T.D.-B.F.X (1994). Fmr1 knockout mice: a model to study fragile X mental retardation. The Dutch-Belgian fragile X consortium. Cells 78, 23–33.
Courchesne, E., Pramparo, T., Gazestani, V. H., Lombardo, M. V., Pierce, K., and Lewis, N. E. (2019). The ASD living biology: from cell proliferation to clinical phenotype. Mol. Psychiatry 24, 88–107. doi: 10.1038/s41380-018-0056-y
Cox, J., Neuhauser, N., Michalski, A., Scheltema, R. A., Olsen, J. V., and Mann, M. (2011). Andromeda: a peptide search engine integrated into the MaxQuant environment. J. Proteome Res. 10, 1794–1805. doi: 10.1021/pr101065j
Croning, M. D., Marshall, M. C., McLaren, P., Armstrong, J. D., and Grant, S. G. N. (2009). G2Cdb: the genes to cognition database. Nucleic Acids Res. 37, D846–D851. doi: 10.1093/nar/gkn700
Cukier, H. N., Dueker, N. D., Slifer, S. H., Lee, J. M., Whitehead, P. L., Lalanne, E., et al. (2014). Exome sequencing of extended families with autism reveals genes shared across neurodevelopmental and neuropsychiatric disorders. Mol. Autism. 5:1. doi: 10.1186/2040-2392-5-1
Darnell, J. C., van Driesche, S. J., Zhang, C., Hung, K. Y. S., Mele, A., Fraser, C. E., et al. (2011). FMRP stalls ribosomal translocation on mRNAs linked to synaptic function and autism. Cells 146, 247–261. doi: 10.1016/j.cell.2011.06.013
de la Torre-Ubieta, L., Won, H., Stein, J. L., and Geschwind, D. H. (2016). Advancing the understanding of autism disease mechanisms through genetics. Nat. Med. 22, 345–361. doi: 10.1038/nm.4071
De Rubeis, S., He, X., Goldberg, A. P., Poultney, C. S., Samocha, K., Ercument Cicek, A., et al. (2014). Synaptic, transcriptional and chromatin genes disrupted in autism. Nature 515, 209–215. doi: 10.1038/nature13772
de Rubeis, S., Pasciuto, E., Li, K. W., Fernández, E., di Marino, D., Buzzi, A., et al. (2013). CYFIP1 coordinates mRNA translation and cytoskeleton remodeling to ensure proper dendritic spine formation. Neuron 79, 1169–1182. doi: 10.1016/j.neuron.2013.06.039
Del Pino, I., Rico, B., and Marin, O. (2018). Neural circuit dysfunction in mouse models of neurodevelopmental disorders. Curr. Opin. Neurobiol. 48, 174–182. doi: 10.1016/j.conb.2017.12.013
Doherty, M. K., Hammond, D. E., Clague, M. J., Gaskell, S. J., and Beynon, R. J. (2009). Turnover of the human proteome: determination of protein intracellular stability by dynamic SILAC. J. Proteome Res. 8, 104–112. doi: 10.1021/pr800641v
Dong, D., Zielke, H. R., Yeh, D., and Yang, P. (2018). Cellular stress and apoptosis contribute to the pathogenesis of autism spectrum disorder. Autism Res. 11, 1076–1090. doi: 10.1002/aur.1966
Erdjument-Bromage, H., Huang, F. K., and Neubert, T. A. (2018). Sample preparation for relative quantitation of proteins using tandem mass tags (TMT) and mass spectrometry (MS). Methods Mol. Biol. 1741, 135–149. doi: 10.1007/978-1-4939-7659-1_11
Forés-Martos, J., Catalá-López, F., Sánchez-Valle, J., Ibáñez, K., Tejero, H., Palma-Gudiel, H., et al. (2019). Transcriptomic metaanalyses of autistic brains reveals shared gene expression and biological pathway abnormalities with cancer. Mol. Autism. 10:17. doi: 10.1186/s13229-019-0262-8
Gao, R., and Penzes, P. (2015). Common mechanisms of excitatory and inhibitory imbalance in schizophrenia and autism spectrum disorders. Curr. Mol. Med. 15, 146–167. doi: 10.2174/1566524015666150303003028
Gaugler, T., Klei, L., Sanders, S. J., Bodea, C. A., Goldberg, A. P., Lee, A. B., et al. (2014). Most genetic risk for autism resides with common variation. Nat. Genet. 46, 881–885. doi: 10.1038/ng.3039
Giulivi, C., Zhang, Y. F., Omanska-Klusek, A., Ross-Inta, C., Wong, S., Hertz-Picciotto, I., et al. (2010). Mitochondrial dysfunction in autism. JAMA 304, 2389–2396. doi: 10.1001/jama.2010.1706
Glessner, J. T., Wang, K., Cai, G., Korvatska, O., Kim, C. E., Wood, S., et al. (2009). Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature 459, 569–573. doi: 10.1038/nature07953
Golden, C. E., Buxbaum, J. D., and De Rubeis, S. (2018). Disrupted circuits in mouse models of autism spectrum disorder and intellectual disability. Curr. Opin. Neurobiol. 48, 106–112. doi: 10.1016/j.conb.2017.11.006
Greenbaum, D., Colangelo, C., Williams, K., and Gerstein, M. (2003). Comparing protein abundance and mRNA expression levels on a genomic scale. Genome Biol. 4:117. doi: 10.1186/gb-2003-4-9-117
Guo, D., Yang, X., and Shi, L. (2020). Rho GTPase regulators and effectors in autism Spectrum disorders: animal models and insights for therapeutics. Cells 9:835. doi: 10.3390/cells9040835
Hallmayer, J., Cleveland, S., Torres, A., Phillips, J., Cohen, B., Torigoe, T., et al. (2011). Genetic heritability and shared environmental factors among twin pairs with autism. Arch. Gen. Psychiatry 68, 1095–1102. doi: 10.1001/archgenpsychiatry.2011.76
Heise, C., Preuss, J. M., Schroeder, J. C., Battaglia, C. R., Kolibius, J., Schmid, R., et al. (2018). Heterogeneity of cell surface glutamate and GABA receptor expression in shank and CNTN4 autism mouse models. Front. Mol. Neurosci. 11:212. doi: 10.3389/fnmol.2018.00212
Hernandez, L. M., Kim, M., Hoftman, G. D., Haney, J. R., de la Torre-Ubieta, L., Pasaniuc, B., et al. (2020). Transcriptomic insight into the polygenic mechanisms underlying psychiatric disorders. Biol. Psychiatry 89, 54–64. doi: 10.1016/j.biopsych.2020.06.005
Hoogendijk, A. J., Pourfarzad, F., Aarts, C. E. M., Tool, A. T. J., Hiemstra, I. H., Grassi, L., et al. (2019). Dynamic Transcriptome-proteome correlation networks reveal human myeloid differentiation and neutrophil-specific programming. Cell Rep. 29, 2505–2519.e4. doi: 10.1016/j.celrep.2019.10.082
Horn, S., Au, M., Basel-Salmon, L., Bayrak-Toydemir, P., Chapin, A., Cohen, L., et al. (2019). De novo variants in PAK1 lead to intellectual disability with macrocephaly and seizures. Brain 142, 3351–3359. doi: 10.1093/brain/awz264
Huang, F. K., Zhang, G., Lawlor, K., Nazarian, A., Philip, J., Tempst, P., et al. (2017). Deep coverage of global protein expression and phosphorylation in breast tumor cell lines using TMT 10-plex isobaric labeling. J. Proteome Res. 16, 1121–1132. doi: 10.1021/acs.jproteome.6b00374
Hughes, H. K., Mills Ko, E., Rose, D., and Ashwood, P. (2018). Immune dysfunction and autoimmunity as pathological mechanisms in autism Spectrum disorders. Front. Cell. Neurosci. 12:405. doi: 10.3389/fncel.2018.00405
Hulbert, S. W., and Jiang, Y. H. (2016). Monogenic mouse models of autism spectrum disorders: common mechanisms and missing links. Neuroscience 321, 3–23. doi: 10.1016/j.neuroscience.2015.12.040
Hussman, J. P. (2001). Suppressed GABAergic inhibition as a common factor in suspected etiologies of autism. J. Autism Dev. Disord. 31, 247–248. doi: 10.1023/A:1010715619091
Hussman, J. P., Chung, R. H., Griswold, A. J., Jaworski, J. M., Salyakina, D., Ma, D., et al. (2011). A noise-reduction GWAS analysis implicates altered regulation of neurite outgrowth and guidance in autism. Mol. Autism. 2:1. doi: 10.1186/2040-2392-2-1
Iakoucheva, L. M., Muotri, A. R., and Sebat, J. (2019). Getting to the cores of autism. Cells 178, 1287–1298. doi: 10.1016/j.cell.2019.07.037
Ingolia, N. T., Ghaemmaghami, S., Newman, J. R. S., and Weissman, J. S. (2009). Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science 324, 218–223. doi: 10.1126/science.1168978
Iossifov, I., O’Roak, B. J., Sanders, S. J., Ronemus, M., Krumm, N., Levy, D., et al. (2014). The contribution of de novo coding mutations to autism spectrum disorder. Nature 515, 216–221. doi: 10.1038/nature13908
Ito, Y., Carss, K. J., Duarte, S. T., Hartley, T., Keren, B., Kurian, M. A., et al. (2018). De novo truncating mutations in WASF1 cause intellectual disability with seizures. Am. J. Hum. Genet. 103, 144–153. doi: 10.1016/j.ajhg.2018.06.001
Jordan, B. A., Fernholz, B. D., Boussac, M., Xu, C., Grigorean, G., Ziff, E. B., et al. (2004). Identification and verification of novel rodent postsynaptic density proteins. Mol. Cell. Proteomics 3, 857–871. doi: 10.1074/mcp.M400045-MCP200
Jordan, B. A., Fernholz, B. D., Khatri, L., and Ziff, E. B. (2007). Activity-dependent AIDA-1 nuclear signaling regulates nucleolar numbers and protein synthesis in neurons. Nat. Neurosci. 10, 427–435. doi: 10.1038/nn1867
Kabitzke, P. A., Brunner, D., He, D., Fazio, P. A., Cox, K., Sutphen, J., et al. (2018). Comprehensive analysis of two Shank3 and the Cacna1c mouse models of autism spectrum disorder. Genes Brain Behav. 17, 4–22. doi: 10.1111/gbb.12405
Kaizuka, T., and Takumi, T. (2018). Postsynaptic density proteins and their involvement in neurodevelopmental disorders. J. Biochem. 163, 447–455. doi: 10.1093/jb/mvy022
Katrancha, S. M., Wu, Y., Zhu, M., Eipper, B. A., Koleske, A. J., and Mains, R. E. (2017). Neurodevelopmental disease-associated de novo mutations and rare sequence variants affect TRIO GDP/GTP exchange factor activity. Hum. Mol. Genet. 26, 4728–4740. doi: 10.1093/hmg/ddx355
Kazdoba, T. M., Leach, P. T., and Crawley, J. N. (2016). Behavioral phenotypes of genetic mouse models of autism. Genes Brain Behav. 15, 7–26. doi: 10.1111/gbb.12256
Koopmans, F., van Nierop, P., Andres-Alonso, M., Byrnes, A., Cijsouw, T., Coba, M. P., et al. (2019). SynGO: an evidence-based, expert-curated Knowledge Base for the synapse. Neuron 103, 217–234.e4. doi: 10.1016/j.neuron.2019.05.002
Krämer, A., Green, J., Pollard, J. Jr., and Tugendreich, S. (2014). Causal analysis approaches in ingenuity pathway analysis. Bioinformatics 30, 523–530. doi: 10.1093/bioinformatics/btt703
Lee, Y., Kang, H., Lee, B., Zhang, Y., Kim, Y., Kim, S., et al. (2017). Integrative analysis of brain region-specific Shank3 Interactomes for understanding the heterogeneity of neuronal pathophysiology related to SHANK3 mutations. Front. Mol. Neurosci. 10:110. doi: 10.3389/fnmol.2017.00110
Lee, H., Thacker, S., Sarn, N., Dutta, R., and Eng, C. (2019). Constitutional mislocalization of Pten drives precocious maturation in oligodendrocytes and aberrant myelination in model of autism spectrum disorder. Transl. Psychiatry 9:13. doi: 10.1038/s41398-018-0364-7
Li, Y., Missig, G., Finger, B. C., Landino, S. M., Alexander, A. J., Mokler, E. L., et al. (2018). Maternal and early postnatal immune activation produce dissociable effects on neurotransmission in mPFC-amygdala circuits. J. Neurosci. 38, 3358–3372. doi: 10.1523/JNEUROSCI.3642-17.2018
Li, J., Shi, M., Ma, Z., Zhao, S., Euskirchen, G., Ziskin, J., et al. (2014). Integrated systems analysis reveals a molecular network underlying autism spectrum disorders. Mol. Syst. Biol. 10:774. doi: 10.15252/msb.20145487
Li, J., Wilkinson, B., Clementel, V. A., Hou, J., O’Dell, T. J., and Coba, M. P. (2016). Long-term potentiation modulates synaptic phosphorylation networks and reshapes the structure of the postsynaptic interactome. Sci. Signal. 9:rs8. doi: 10.1126/scisignal.aaf6716
Li, J., Zhang, W., Yang, H., Howrigan, D. P., Wilkinson, B., Souaiaia, T., et al. (2017). Spatiotemporal profile of postsynaptic interactomes integrates components of complex brain disorders. Nat. Neurosci. 20, 1150–1161. doi: 10.1038/nn.4594
Louros, S. R., and Osterweil, E. K. (2016). Perturbed proteostasis in autism spectrum disorders. J. Neurochem. 139, 1081–1092. doi: 10.1111/jnc.13723
Lugo, J. N., Smith, G. D., Arbuckle, E. P., White, J., Holley, A. J., Floruta, C. M., et al. (2014). Deletion of PTEN produces autism-like behavioral deficits and alterations in synaptic proteins. Front. Mol. Neurosci. 7:27. doi: 10.3389/fnmol.2014.00027
Ma, B., Shan, X., Yu, J., Zhu, T., Li, R., Lv, H., et al. (2022). Social deficits via dysregulated Rac1-dependent excitability control of prefrontal cortical neurons and increased GABA/glutamate ratios. Cell Rep. 41:111722. doi: 10.1016/j.celrep.2022.111722
Mabb, A. M., Judson, M. C., Zylka, M. J., and Philpot, B. D. (2011). Angelman syndrome: insights into genomic imprinting and neurodevelopmental phenotypes. Trends Neurosci. 34, 293–303. doi: 10.1016/j.tins.2011.04.001
Maier, T., Guell, M., and Serrano, L. (2009). Correlation of mRNA and protein in complex biological samples. FEBS Lett. 583, 3966–3973. doi: 10.1016/j.febslet.2009.10.036
Martin, H. G., and Manzoni, O. J. (2014). Late onset deficits in synaptic plasticity in the valproic acid rat model of autism. Front. Cell. Neurosci. 8:23. doi: 10.3389/fncel.2014.00023
McFarlane, H. G., Kusek, G. K., Yang, M., Phoenix, J. L., Bolivar, V. J., and Crawley, J. N. (2008). Autism-like behavioral phenotypes in BTBR T+tf/J mice. Genes Brain Behav. 7, 152–163. doi: 10.1111/j.1601-183X.2007.00330.x
Meltzer, A., and Van de Water, J. (2017). The role of the immune system in autism Spectrum disorder. Neuropsychopharmacology 42, 284–298. doi: 10.1038/npp.2016.158
Murtaza, N., Uy, J., and Singh, K. K. (2020). Emerging proteomic approaches to identify the underlying pathophysiology of neurodevelopmental and neurodegenerative disorders. Mol. Autism. 11:27. doi: 10.1186/s13229-020-00334-5
Nelson, S. B., and Valakh, V. (2015). Excitatory/inhibitory balance and circuit homeostasis in autism Spectrum disorders. Neuron 87, 684–698. doi: 10.1016/j.neuron.2015.07.033
Pacheco, N. L., Heaven, M. R., Holt, L. M., Crossman, D. K., Boggio, K. J., Shaffer, S. A., et al. (2017). RNA sequencing and proteomics approaches reveal novel deficits in the cortex of Mecp2-deficient mice, a model for Rett syndrome. Mol. Autism. 8:56. doi: 10.1186/s13229-017-0174-4
Parikshak, N. N., Swarup, V., Belgard, T. G., Irimia, M., Ramaswami, G., Gandal, M. J., et al. (2016). Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism. Nature 540, 423–427. doi: 10.1038/nature20612
Patrich, E., Piontkewitz, Y., Peretz, A., Weiner, I., and Attali, B. (2016). Maternal immune activation produces neonatal excitability defects in offspring hippocampal neurons from pregnant rats treated with poly I:C. Sci. Rep. 6:19106. doi: 10.1038/srep19106
Peñagarikano, O., Abrahams, B. S., Herman, E. I., Winden, K. D., Gdalyahu, A., Dong, H., et al. (2011). Absence of CNTNAP2 leads to epilepsy, neuronal migration abnormalities, and core autism-related deficits. Cells 147, 235–246. doi: 10.1016/j.cell.2011.08.040
Pinto, D., Delaby, E., Merico, D., Barbosa, M., Merikangas, A., Klei, L., et al. (2014). Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am. J. Hum. Genet. 94, 677–694. doi: 10.1016/j.ajhg.2014.03.018
Pinto, D., Pagnamenta, A. T., Klei, L., Anney, R., Merico, D., Regan, R., et al. (2010). Functional impact of global rare copy number variation in autism spectrum disorders. Nature 466, 368–372. doi: 10.1038/nature09146
Quesnel-Vallières, M., Weatheritt, R. J., Cordes, S. P., and Blencowe, B. J. (2019). Autism spectrum disorder: insights into convergent mechanisms from transcriptomics. Nat. Rev. Genet. 20, 51–63. doi: 10.1038/s41576-018-0066-2
Ramaswami, G., and Geschwind, D. H. (2018). Genetics of autism spectrum disorder. Handb. Clin. Neurol. 147, 321–329. doi: 10.1016/B978-0-444-63233-3.00021-X
Rappsilber, J., Mann, M., and Ishihama, Y. (2007). Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat. Protoc. 2, 1896–1906. doi: 10.1038/nprot.2007.261
Reim, D., Distler, U., Halbedl, S., Verpelli, C., Sala, C., Bockmann, J., et al. (2017). Proteomic analysis of post-synaptic density fractions from Shank3 mutant mice reveals brain region specific changes relevant to autism Spectrum disorder. Front. Mol. Neurosci. 10:26. doi: 10.3389/fnmol.2017.00026
Ronemus, M., Iossifov, I., Levy, D., and Wigler, M. (2014). The role of de novo mutations in the genetics of autism spectrum disorders. Nat. Rev. Genet. 15, 133–141. doi: 10.1038/nrg3585
Rossignol, D. A., and Frye, R. E. (2012). Mitochondrial dysfunction in autism spectrum disorders: a systematic review and meta-analysis. Mol. Psychiatry 17, 290–314. doi: 10.1038/mp.2010.136
Rubenstein, J. L., and Merzenich, M. M. (2003). Model of autism: increased ratio of excitation/inhibition in key neural systems. Genes Brain Behav. 2, 255–267. doi: 10.1034/j.1601-183X.2003.00037.x
Ruzzo, E. K., Pérez-Cano, L., Jung, J. Y., Wang, L. K., Kashef-Haghighi, D., Hartl, C., et al. (2019). Inherited and De novo genetic risk for autism impacts shared networks. Cells 178, 850–866.e26. doi: 10.1016/j.cell.2019.07.015
Sadybekov, A., Tian, C., Arnesano, C., Katritch, V., and Herring, B. E. (2017). An autism spectrum disorder-related de novo mutation hotspot discovered in the GEF1 domain of Trio. Nat. Commun. 8:601. doi: 10.1038/s41467-017-00472-0
Sanders, S. J., Murtha, M. T., Gupta, A. R., Murdoch, J. D., Raubeson, M. J., Willsey, A. J., et al. (2012). De novo mutations revealed by whole-exome sequencing are strongly associated with autism. Nature 485, 237–241. doi: 10.1038/nature10945
Sarowar, T., and Grabrucker, A. M. (2016). Actin-dependent alterations of dendritic spine morphology in Shankopathies. Neural Plast. 2016:8051861. doi: 10.1155/2016/8051861
Schoen, M., Asoglu, H., Bauer, H. F., Müller, H. P., Abaei, A., Sauer, A. K., et al. (2019). Shank3 transgenic and prenatal zinc-deficient autism mouse models show convergent and individual alterations of brain structures in MRI. Front. Neural. Circuits 13:6. doi: 10.3389/fncir.2019.00006
Schwanhäusser, B., Busse, D., Li, N., Dittmar, G., Schuchhardt, J., Wolf, J., et al. (2011). Global quantification of mammalian gene expression control. Nature 473, 337–342. doi: 10.1038/nature10098
Sestan, N., and State, M. W. (2018). Lost in translation: traversing the complex path from genomics to therapeutics in autism Spectrum disorder. Neuron 100, 406–423. doi: 10.1016/j.neuron.2018.10.015
Silverman, J. L., and Ellegood, J. (2018). Behavioral and neuroanatomical approaches in models of neurodevelopmental disorders: opportunities for translation. Curr. Opin. Neurol. 31, 126–133. doi: 10.1097/WCO.0000000000000537
Steffenburg, S., Gillberg, C., Hellgren, L., Andersson, L., Gillberg, I. C., Jakobsson, G., et al. (1989). A twin study of autism in Denmark, Finland, Iceland, Norway and Sweden. J. Child Psychol. Psychiatry 30, 405–416. doi: 10.1111/j.1469-7610.1989.tb00254.x
Stessman, H. A., Turner, T. N., and Eichler, E. E. (2016). Molecular subtyping and improved treatment of neurodevelopmental disease. Genome Med. 8:22. doi: 10.1186/s13073-016-0278-z
Sui, L., and Chen, M. (2012). Prenatal exposure to valproic acid enhances synaptic plasticity in the medial prefrontal cortex and fear memories. Brain Res. Bull. 87, 556–563. doi: 10.1016/j.brainresbull.2012.01.011
Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447–D452. doi: 10.1093/nar/gku1003
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Sztainberg, Y., and Zoghbi, H. Y. (2016). Lessons learned from studying syndromic autism spectrum disorders. Nat. Neurosci. 19, 1408–1417. doi: 10.1038/nn.4420
Tang, B., Wang, T., Wan, H., Han, L., Qin, X., Zhang, Y., et al. (2015). Fmr1 deficiency promotes age-dependent alterations in the cortical synaptic proteome. Proc. Natl. Acad. Sci. U. S. A. 112, E4697–E4706. doi: 10.1073/pnas.1502258112
Thompson, A., Schäfer, J., Kuhn, K., Kienle, S., Schwarz, J., Schmidt, G., et al. (2003). Tandem mass tags: a novel quantification strategy for comparative analysis of complex protein mixtures by MS/MS. Anal. Chem. 75, 1895–1904. doi: 10.1021/ac0262560
Tian, C., Kay, Y., Sadybekov, A., Rao, S., Katritch, V., and Herring, B. E. (2018). An intellectual disability-related missense mutation in Rac1 prevents LTP induction. Front. Mol. Neurosci. 11:223. doi: 10.3389/fnmol.2018.00223
Tindi, J. O., Chavez, A. E., Cvejic, S., Calvo-Ochoa, E., Castillo, P. E., and Jordan, B. A. (2015). ANKS1B gene product AIDA-1 controls hippocampal synaptic transmission by regulating GluN2B subunit localization. J. Neurosci. 35, 8986–8996. doi: 10.1523/JNEUROSCI.4029-14.2015
Toma, C., Torrico, B., Hervás, A., Valdés-Mas, R., Tristán-Noguero, A., Padillo, V., et al. (2014). Exome sequencing in multiplex autism families suggests a major role for heterozygous truncating mutations. Mol. Psychiatry 19, 784–790. doi: 10.1038/mp.2013.106
Tyanova, S., Temu, T., and Cox, J. (2016). The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 11, 2301–2319. doi: 10.1038/nprot.2016.136
Varghese, M., Keshav, N., Jacot-Descombes, S., Warda, T., Wicinski, B., Dickstein, D. L., et al. (2017). Autism spectrum disorder: neuropathology and animal models. Acta Neuropathol. 134, 537–566. doi: 10.1007/s00401-017-1736-4
Verma, V., Paul, A., Amrapali Vishwanath, A., Vaidya, B., and Clement, J. P. (2019). Understanding intellectual disability and autism spectrum disorders from common mouse models: synapses to behaviour. Open Biol. 9:180265. doi: 10.1098/rsob.180265
Vogel, C., and Marcotte, E. M. (2012). Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat. Rev. Genet. 13, 227–232. doi: 10.1038/nrg3185
Wang, J., Ma, Z., Carr, S. A., Mertins, P., Zhang, H., Zhang, Z., et al. (2017). Proteome profiling outperforms Transcriptome profiling for Coexpression based gene function prediction. Mol. Cell. Proteomics 16, 121–134. doi: 10.1074/mcp.M116.060301
Wang, R., Tan, J., Guo, J., Zheng, Y., Han, Q., So, K. F., et al. (2018). Aberrant development and synaptic transmission of cerebellar cortex in a VPA induced mouse autism model. Front. Cell. Neurosci. 12:500. doi: 10.3389/fncel.2018.00500
Weiner, D. J., Wigdor, E. M., Ripke, S., Walters, R. K., Kosmicki, J. A., Grove, J., et al. (2017). Polygenic transmission disequilibrium confirms that common and rare variation act additively to create risk for autism spectrum disorders. Nat. Genet. 49, 978–985. doi: 10.1038/ng.3863
Yoo, H. (2015). Genetics of autism Spectrum disorder: current status and possible clinical applications. Exp. Neurobiol. 24, 257–272. doi: 10.5607/en.2015.24.4.257
Yu, T. W., Chahrour, M. H., Coulter, M. E., Jiralerspong, S., Okamura-Ikeda, K., Ataman, B., et al. (2013). Using whole-exome sequencing to identify inherited causes of autism. Neuron 77, 259–273. doi: 10.1016/j.neuron.2012.11.002
Yuen, R. K., Thiruvahindrapuram, B., Merico, D., Walker, S., Tammimies, K., Hoang, N., et al. (2015). Whole-genome sequencing of quartet families with autism spectrum disorder. Nat. Med. 21, 185–191. doi: 10.1038/nm.3792
Zamboni, V., Armentano, M., Berto, G., Ciraolo, E., Ghigo, A., Garzotto, D., et al. (2018). Hyperactivity of Rac1-GTPase pathway impairs neuritogenesis of cortical neurons by altering actin dynamics. Sci. Rep. 8:7254. doi: 10.1038/s41598-018-25354-3
Zeidán-Chuliá, F., Rybarczyk-Filho, J. L., Salmina, A. B., de Oliveira, B. H. N., Noda, M., and Moreira, J. C. F. (2013). Exploring the multifactorial nature of autism through computational systems biology: calcium and the rho GTPase RAC1 under the spotlight. NeuroMolecular Med. 15, 364–383. doi: 10.1007/s12017-013-8224-3
Zhang, G., Neubert, T. A., and Jordan, B. A. (2012). RNA binding proteins accumulate at the postsynaptic density with synaptic activity. J. Neurosci. 32, 599–609. doi: 10.1523/JNEUROSCI.2463-11.2012
Zhou, Y., Kaiser, T., Monteiro, P., Zhang, X., van der Goes, M. S., Wang, D., et al. (2016). Mice with Shank3 mutations associated with ASD and schizophrenia display both shared and distinct defects. Neuron 89, 147–162. doi: 10.1016/j.neuron.2015.11.023
Keywords: proteomic convergence, PSD, Rho GTPase, Rac, tandem mass tags, TMT, comparative analyses
Citation: Carbonell AU, Freire-Cobo C, Deyneko IV, Dobariya S, Erdjument-Bromage H, Clipperton-Allen AE, Page DT, Neubert TA and Jordan BA (2023) Comparing synaptic proteomes across five mouse models for autism reveals converging molecular similarities including deficits in oxidative phosphorylation and Rho GTPase signaling. Front. Aging Neurosci. 15:1152562. doi: 10.3389/fnagi.2023.1152562
Edited by:
Dario Maschi, Washington University in St. Louis, United StatesReviewed by:
Sehyoun Yoon, Northwestern Medicine, United StatesTilmann Achsel, Université de Lausanne, Switzerland
Copyright © 2023 Carbonell, Freire-Cobo, Deyneko, Dobariya, Erdjument-Bromage, Clipperton-Allen, Page, Neubert and Jordan. 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: Bryen A. Jordan, YnJ5ZW4uam9yZGFuQGVpbnN0ZWlubWVkLmVkdQ==
†Present addressesAbigail U. Carbonell, Department of Neurology, Columbia University Irving Medical Center, New York, NY, United States
Carmen Freire-Cobo, Nash Family Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, United States