- Research article
- Open access
- Published:
Trem2 Y38C mutation and loss of Trem2 impairs neuronal synapses in adult mice
Molecular Neurodegeneration volume 15, Article number: 62 (2020)
Abstract
Background
Triggering receptor expressed on myeloid cells 2 (TREM2) is expressed in the brain exclusively on microglia and genetic variants are linked to neurodegenerative diseases including Alzheimer’s disease (AD), frontotemporal dementia (FTD) and Nasu Hakola Disease (NHD). The Trem2 variant R47H, confers substantially elevated risk of developing late onset Alzheimer’s disease, while NHD-linked Trem2 variants like Y38C, are associated with development of early onset dementia with white matter pathology. However, it is not known how these Trem2 species, predisposes individuals to presenile dementia.
Methods
To investigate if Trem2 Y38C or loss of Trem2 alters neuronal function we generated a novel mouse model to introduce the NHD Trem2 Y38C variant in murine Trem2 using CRISPR/Cas9 technology. Trem2Y38C/Y38C and Trem2−/− mice were assessed for Trem2 expression, differentially expressed genes, synaptic protein levels and synaptic plasticity using biochemical, electrophysiological and transcriptomic approaches.
Results
While mice harboring the Trem2 Y38C exhibited normal expression levels of TREM2, the pathological outcomes phenocopied Trem2−/− mice at 6 months. Transcriptomic analysis revealed altered expression of neuronal and oligodendrocytes/myelin genes. We observed regional decreases in synaptic protein levels, with the most affected synapses in the hippocampus. These alterations were associated with reduced synaptic plasticity.
Conclusion
Our findings provide in vivo evidence that Trem2 Y38C disrupts normal TREM2 functions. Trem2Y38C/Y38C and Trem2−/− mice demonstrated altered gene expression, changes in microglia morphology, loss of synaptic proteins and reduced hippocampal synaptic plasticity at 6 months in absence of any pathological triggers like amyloid. This suggests TREM2 impacts neuronal functions providing molecular insights on the predisposition of individuals with TREM2 variants resulting in presenile dementia.
Background
Microglia are brain-resident macrophages that are appreciated for their dynamic nature during inflammation, aging and disease. Recent genome-wide association studies have identified variants of microglial genes that are associated with high risk of developing neurodegenerative diseases, including variants of triggering receptor expressed on myeloid cells 2 (TREM2) [1,2,3,4]. TREM2, a single-pass trans-membrane glycoprotein with an extracellular V-type Ig domain and a connecting stalk followed by a transmembrane region and a C-terminal tail. In the brain, TREM2 expression is restricted to microglia [5] and acts to modulate a spectrum of innate immune functions. Upon ligand binding it associates with DNAX-activation protein 12 (DAP12) to initiate downstream signaling [2].
Homozygous loss of function mutations in TREM2 or in its cognate signaling element DAP12, has been reported to cause Nasu Hakola disease (NHD), an autosomal recessive disorder, characterized by early onset dementia and bone cysts [6]. Individuals with the rare R47H variant of TREM2 have substantial risk of developing late onset Alzheimer’s disease (AD) [7, 8]. This mutation alters the ligand binding region of TREM2, conferring loss of function [9]. Other TREM2 variants such as Q33X, T66M and Y38C are reported to confer loss of function and have been linked to the development of a frontotemporal dementia (FTD)-like syndrome with early-onset dementia [10]. In vivo, T66M knock-in mice exhibit reduced microglial phagocytosis, cerebral blood flow and brain glucose metabolism [11, 12]. In vitro studies have shown that Trem2 Y38C mutation have altered protein conformation, which results in retention of TREM2 within the endoplasmic reticulum and less translocation to the plasma membrane [9, 13, 14]. However, it’s in vivo expression pattern remains unknown. A compound heterozygous mutation p.[(Y38C)];[(D86V)] has also been reported in a patient who presented with an FTD-like syndrome [15]. Trem2 Y38C mutation has also been shown to have altered glycosylation, diminished binding to TREM2 ligands, and impaired phagocytosis in vitro [12, 16,17,18]. However, the effects of Y38C, in vivo, on brain homeostasis in the presence and absence of pathology have yet to be determined.
Beyond the association with neurodegenerative diseases, TREM2 has also been shown to be crucial in brain development. Recent work indicates that Trem2 deficient mice have increased synapse density at P20 suggesting that TREM2 controls supernumerary synapse pruning by microglia [19]. By 3 months of age, absence of TREM2 leads to altered neuronal connections and an autism spectrum disease-like behavioral phenotype [19]; however, TREM2 mediated mechanisms underlying this neuronal pathology remains unclear. In support of a developmental role for TREM2, the absence of TREM2 in Velocigene Trem2 knock-out mouse model decreased apposed synapse numbers in the motor cortex at 1 month of age, presumably due to altered communication between astrocytes and Trem2-deficient microglia [20]. However, this deficit was normalized by 4 months of age. In the current study, we examine the role of TREM2 during both development and adult stages using two novel mouse models to elucidate how TREM2 influences neuronal functions and can predispose neurons to synaptic alterations. To address this, in the absence of pathological insults like amyloid, we used mice expressing Trem2 Y38C (Trem2Y38C/Y38C) along with those in which Trem2 gene has been knocked out on a C57BL/6 background (Trem2−/−), allowing determination of whether the Y38C variant exhibited altered function. Using these mouse models, transcriptomics data analysis at 6 months indicated alterations in genes and pathways associated with neuronal functions and oligodendrocyte/myelin, thus presenting mice models that mimic physiological conditions of patients with TREM2 variants with FTD-like dementia [10]. We observed expression of TREM2 Y38C and Trem2 deficiency leads to loss of synaptic proteins post developmentally. Vulnerability of pre-synaptic and post-synaptic elements showed region specific differences, with the hippocampus being affected the most. Furthermore, the loss of synaptic protein impacts hippocampal synaptic plasticity indicating effects of Trem2Y38C/Y38C and Trem2−/− on neuronal functions. Overall, this study demonstrates loss of functional TREM2 shows age dependent, brain region specific synaptic alterations in absence of pathological triggers like amyloid. These results provide insights on how NHD Trem2 variant predisposes individual to early onset dementia.
Methods
Experimental model
To generate the Trem2Y38C/Y38C mouse model, CRISPR/Cas9-mediated insertion of the SNP encoding the Trem2 Y38C variant into the mouse Trem2 gene was performed as previously reported [21]. Briefly, embryos were injected with Cas9, short-guide RNA (sgRNA) and replacement oligo. The relevant sequences are as follows: Trem2 targeted region 3′-TACTGCGGAACTTCGTGACCCCC, sgRNA (antisense) 5′- ATGACGCCTTGAAGCACTGGGGG and replacement oligo 3′- TGTGACGCCTTGAAGCACTGGGGA. The first codon in the replacement oligo corresponds to the SNP encoding the Y38C variant, while the third codon corresponds to a silent mutation that ablates the protospacer adjacent motif (PAM) site, necessary for initial binding of CRISPR/Cas9. To identify mice containing the Trem2 Y38C mutation, Sanger sequencing of founder lines was used to identify mice carrying heterozygous or homozygous Trem2 Y38C SNP. Five founder lines were selected to generate subsequent crosses. SNP-based genotyping (Thermo Fisher) of offspring was performed using the following primers: forward primer: 5′-GCCGGCCAGTCCTTGAG, reverse primer: 5′-CACCAGGCCTTGCGTCT, SNP reporter 1: 5′-CAAGGCGTCATAAGTACA, SNP reporter 2: 5′-AGGCGTCACAAGTACA. Trem2 Y38C variant carrying offspring of the founder lines were maintained on a C57BL6/J (B6) background.
Trem2−/− (JAX stock #027197) and B6 wildtype (WT) mice were obtained from the Jackson Laboratories (JAX stock #000664). Both male and female mice at ages postnatal day 20 (P20) and 6 months were utilized in this study.
To analyze CRISPR/Cas9 editing efficiency, potential off target mutations were predicted using prediction tool Synthego (https://www.synthego.com/). The top 4 predicted mutations were considered for analysis. Genomic DNA was extracted using 50 mM NaOH from the tail snips of 3 F9 (mice used for the experiments) Trem2Y38C/Y38C and one WT. Genomic regions with the predicted off-target sites were amplified (Additional file 1: Table S1) and Sanger’s DNA sequencing was performed (GENWIZ) (Additional file 2: Table S2). Sequences were then analyzed for the predicted mutation. No predicted off-target mutations were detected in the Trem2Y38C/Y38C samples.
Mice were housed in the Indiana University School of Medicine (IUSM) animal facilities, which are accredited by the Association and Accreditation of Laboratory Animal Care. Animals were maintained according to USDA standards and the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Experiments were approved by the IUSM Institutional Animal Care and Use Committee.
Brain extraction and tissue processing
Mice were deeply anesthetized with ketamine/xylazine and perfused with sterile ice-cold PBS, and their brains were removed. For immunohistochemistry, one hemisphere was drop-fixed in 4% PFA in PBS for 24–48 h, transferred to ice-cold 30% sucrose and stored at 4 °C in in PBS. After embedding in OCT Compound (VWR), 30 μm thick sections were obtained on a Leica CM 1950 cryostat. Cortex and hippocampus was microdissected from the other half of the brain, snap-frozen and stored at − 80 °C until use. Frozen brain tissue (cortex or hippocampus) was homogenized in ice-cold T-PER homogenization buffer (ThermoFisher, 78,510) supplemented with protein phosphatase (1:100, P5726, Sigma-Aldrich) and protease inhibitor cocktail (1:100, P8340, Sigma-Aldrich). Homogenates were aliquoted and reserved for protein extraction or RNA extraction.
RNA extraction
Homogenized cortical lysates that were prepared in ice-cold T-per homogenization buffer, as described above, were added to an equal volume of QIAzol lysis reagent (QIAGEN RNeasy kits). RNA was isolated using QIAGEN RNeasy kit as per the manufacturer’s protocol and quantified using NanoDrop (Thermo Scientific).
RNA-seq library preparation
Total RNA from cortices of WT, Trem2Y38C/Y38C and Trem2−/− mice were first evaluated for its quantity, and quality, using Agilent Bioanalyzer 2100. One hundred nanograms of total RNA were used. Ribosomal RNA was removed from total RNA using QIAseq FastSelect rRNA Removal HMR Kit (Qiagen, Catalog# 334387 Human/Mouse/Rat). After the depletion of rRNA, cDNA library preparation was carried out including RNA fragmentation, cDNA synthesis, ligation of index adaptors, and amplification, following the KAPA RNA Hyper Prep Kit Technical Data Sheet, KR0961 – v3.15 (Roche Corporate, Catalog #KK8541). Each resulting indexed library was quantified and its quality accessed by Qubit and Agilent Bioanalyzer, and multiple libraries were pooled in equal molarity. The library pool was then sequenced in 100b paired-end read format on NovaSeq 6000 (Illumina, Inc.). More than 30 million reads per sample were generated and 91% of the sequencing reads reached Q30 (99.9% base call accuracy). A Phred quality score (Q score) was used to measure the quality of sequencing.
Mapping QC and data analysis
The sequencing data were first assessed using FastQC (Babraham Bioinformatics, Cambridge, UK) for quality control. All sequenced libraries were mapped to the mouse genome (mm10) (or to the human genome (hg38)) using STAR RNA-seq aligner (v.2.5) [22] with the following parameter: “--outSAMmapqUnique 60”. The reads distribution across the genome was assessed using bamutils (from ngsutils v.0.5.9) [23]. Uniquely mapped sequencing reads were assigned to mm10 refGene genes (or hg38 refGene genes) using featureCounts (from subread v.1.5.1) [24] with the following parameters: “-s 2 –p –Q 10”. Differential expression and visualization were performed using R Bioconductor package DESeq2 by pair-wise comparison of Trem2Y38C/Y38C versus WT and Trem2−/− versus WT. Counts were normalized and examined for outliers in DESeq2. The list was further filtered based on FDR (P-adj) < 0.05 and Log fold changes (logFC), |logFC| < 0.58. All the filtered differentially expressed genes were annotated using Ingenuity Pathway Analysis (IPA) software. Enrichment for Tissue type, Molecular Functions, Cellular Components, and Pathways on the filtered differentially expressed genes were performed using Enrichr web-server tool [25, 26]. The smallest P-value indicates the highest degree of enrichment. Volcano plots and heatmaps were plotted using volcanoPlot and pheatmap functions in R.
Western blotting
For protein extraction, homogenates of microdissected cortices and hippocampi were sonicated and centrifuged at 50,000×g. Protein concentrations were measured using a Bicinchoninic Acid (BCA) assay (Thermo Fisher, 23,225), according to the manufacturer’s instructions. Proteins were denatured at 95 °C for 10 min in 3X denaturing buffer containing LDS sample buffer and 50 mM DTT. 15–20 μg of protein for each sample were loaded onto Novex 4–12% Bis-Tris gels (Invitrogen) and run at 150 V for 1 h in MES running buffer (Thermo Fischer) and transferred onto PVDF membranes (EMD Millipore) in Tris-Glycine transfer buffer containing methanol at 100 V for 2 h on ice. Membranes were blocked in TBST (Tris-buffered saline with 0.1% Tween 20) containing 5% BSA for 1 h at room temperature and incubated in the indicated primary antibodies in blocking buffer overnight at 4 °C with shaking: PSD-95 (NeuroMab, 75–028, 1:20,000), Synaptophysin (Cell Signaling technology, 5461, 1:20,000), Synapsin 2 (Abcam, ab13258, 1:1000), Homer 1 (GeneTex, 103,278, 1:5000), CNPase (abcam, ab6319,1:500), MBP (abcam, ab7349, 1:2000), β actin (Sigma, A1978, 1:200,000). Membranes were then washed in TBST and incubated for 1 h with appropriate secondary antibody diluted in 5% milk in TBST at room temperature. Membranes were imaged using Chemiluminescence HRP substrate (Millipore). Densitometry was performed using Image J.
Immunohistochemistry
30 μm free-floating tissue sections were subjected to antigen retrieval in sodium citrate buffer (10 mM, pH 6) for 10 min at 95 °C. After cooling, sections were blocked in blocking buffer (5% Normal Donkey Serum, 0.3% TritonX-100 in PBS) and incubated with the appropriate primary antibody overnight at 4 °C: Iba-1 (Millipore, MABN92 or Wako, 019–19,741; both at 1:500). After washing, sections were incubated in corresponding Alexa Fluor-conjugated secondary antibodies (1:500). Sections were mounted with ProLong Gold AntiFade with DAPI (Fisher Scientific). For Trem2 immunofluorescent staining, Tyramide Signal Amplification (TSA) Biotin System kit (PerkinElmer, NEL700A001KT) was used. Briefly, following antigen retrieval (as described above), sections were blocked in TNB blocking buffer at room temperature for 1 h. Sections were then incubated in primary antibody for TREM2 (R&D systems, AF1729, 1:300) overnight at 4 °C. After washing in TNT wash buffer, sections were incubated with anti-sheep biotinylated secondary antibody (Vector laboratories, BA-6000, 1:850) mkj for 1 h at room temperature. Sections were washed and incubated in Streptavidin-HRP (1:1000). Sections were washed and then incubated in Tyramide solution (1:250) for 4 min. After washing, the sections were incubated in Streptavidin Alexa Fluor 488 (Invitrogen, S32354, 1:1000) at room temperature for 1 h. After washing, sections were blocked and stained for Iba-1 as described above.
Morphological analysis of Iba-1 positive cells was performed on a total of 4 sections (2 medial tissue sections and 2 lateral sections), from 6 mice (equal number of males and females) per genotype. The number of Iba-1 positive cells and percent area of Iba-1 staining were quantified on 10X images of the entire cortex and hippocampus of a section. After thresholding, soma number and total cell area were quantified in a blinded manner using ImageJ. 31–33 microglia per genotype were selected randomly from the cortical images which were binarized and skeletonized using ImageJ Skeletonize plugin, as described previously [27]. Using Analyze Skeleton feature, total number of branches and junctions per cell were quantified for images from cortex. Immunofluorescent images to visualize TREM2 were acquired with the Nikon AR1 Confocal microscope. Maximum intensity of Z-stacks were obtained using 60X Nyquist view images of microglia from similar cortical regions for all the genotypes to confirm the expression of TREM2 in Iba-1 positive cells.
Enzyme linked immunosorbent assay (ELISA)
To quantify total TREM2 protein levels, cortices from WT, Trem2Y38C/Y38C and Trem2−/− were homogenized in lysis buffer (25 mM Tris pH 7.4, 150 mM NaCl, 1 mM EDTA, 5% glycerol, 1% NP-40) and protein concentration was determined by Bicinchoninic Acid (BCA) assay (Thermo Fisher, 23,225). F8 Maxisorp Nunc-Immuno Module (Thermo Fisher, 468,667) wells were coated with 2 μg/ml of the TREM2 capture antibody (R&D Systems, MAB1729) in 0.05 M carbonate/bicarbonate buffer (pH 9.6), overnight at 4 °C, and blocked with 3% BSA, 0.05% Tween in PBS for 1 h at RT. TREM2 standards were prepared using recombinant mouse TREM2 protein (R&D system, 9228-T2). Standards, and 500 μg of cortical lysates were incubated for 2 h at room temperature. Wells were washed 4 times with 0.05% Tween in PBS and incubated with 0.25μg/ml of the TREM2 biotinylated detection antibody (R&D Systems, BAF1729) for 1 h at RT. After washing, samples were incubated with HRP-conjugated streptavidin (PerkinElmer, NEL750001EA, 1;10,000). The samples were washed and incubated with the chromogenic substrate TMB (3,3′,5,5′-tetramethylbenzidine) (Pierce TMB Substrate kit, ThermoFisher, 34,021). Upon optimal color development, reactions were stopped using 1 N HCL and wells were read at 450 nM using the Epoch2 microplate reader (BioTek). For sTREM2 quantification, frozen mouse cortices were subjected to cryogenic grinding in a mortar and pestle. Cortices were further homogenized in ice-cold Tris-buffered saline (TBS; 50 mM Tris pH 7.4, 150 mM NaCl; 10 w/v) using a 26-G needle and cleared by ultracentrifugation (186,000×g,1 h, 4 °C). Protein concentration was determined using BCA assay and 140 μg of protein was added to the wells. sTREM2 levels were quantified using the ELISA protocol described above.
Quantitative RT-PCR (qRT-PCR)
cDNA was prepared from 700 ng of RNA using the High Capacity RNA-to-cDNA kit (Applied Biosystems, 4,387,406) and quantitative PCR (qPCR) was performed using the QuantStudio 6 Flex (Applied Biosystems). Murine Trem2 mRNA quantification was performed using following primers: Forward 5′- CCTCTCCACCAGTTTCTCCT − 3′ and Reverse 5′- CAGTGCTTCAAGGCGTCATAAG − 3′. Taqman assays used: Mbp (Mm01266402_m1), Olig2 (Mm01210556_m1), Mobp (Mm02745649_m1). Relative gene expression was determined using the ΔΔCT method and graphed as fold change relative to murine GAPDH.
Brain slice preparations
Mice (6-month-old) were deeply anesthetized using isoflurane and sacrificed by rapid decapitation. Brains were removed and placed in a 95% O2 and 5% CO2-saturated, ice-cold tissue cutting solution (194 mM sucrose, 30 mM NaCl, 4.5 mM KCl, 1 mM MgCl2, 26 mM NaHCO3, 1.2 mM NaH2PO4, 10 mM glucose). Sagittal slices containing hippocampus were cut to a thickness of 350 μm using a vibratome (Leica VT1200S) and transferred to an artificial cerebral spinal fluid (aCSF) (124 mM NaCl, 4.5 mM KCl, 1 mM MgCl2, 26 mM NaHCO3, 1.2 mM NaH2PO4, 10 mM glucose, 2 mM CaCl2) saturated with 95% O2 and 5% CO2 at 30 °C for 1 h. Slices were then kept at room temperature until recording.
Field potential recordings
Field excitatory postsynaptic potential (fEPSP) recordings were performed on brain slices containing the stratum radiatum of hippocampal CA1 region in a chamber that was continuously perfused with aCSF at a rate of 1–2 mL/min in 5% CO2 at 30–32 °C. Slices were visualized using an Olympus BX51WI microscope (Olympus Corporation of America). Extracellular recordings were conducted using a Multiclamp 700B amplifier (Molecular Devices). Schaffer collaterals in the hippocampus were stimulated with tungsten stereotrodes (MicroProbes for Life Science). To record responses from hippocampal CA1 neurons, filament-containing borosilicate micropipettes (World Precision Instruments) were prepared using a P-1000 micropipette puller (Sutter Instruments) and filled with 1 M NaCl. A Constant current Isolated Stimulator (Digitimer) was used to produce electrical stimulation.
Input-output (I/O) curves were recorded to determine the maximum stimulus-response, which was determined by calculating the slope of the response (mV/ms). Using the stimulation strength that produced 50% of the maximum intensity, a stable baseline was observed for 10 min before recording the pair pulse ratio (PPR) and field excitatory postsynaptic potential (fEPSPs). The PPRs were obtained every 20 s for a total of 3 min. The fEPSPs were acquired across a 70-min time window: 10 min of pre-stimulation baseline and 60 min post-stimulation. Long term potentiation (LTP) was induced using a high-frequency protocol: 4 trains of 10 pulses, every 20 s at 100 Hz. The investigators were blinded to genotype.
Statistical analysis
False discovery rates (P-adjusted) and P value were calculated using the Benjamini & Hochberg method [28] and Wald test, respectively, within DESeq2. Statistical analyses were performed using Prism (GraphPad) for rest of the data. Statistical significance was determined using a one-way or two-way ANOVA with Bonferroni or Tukey’s post hoc analysis, with P values less than 0.05 considered as significant. Fisher’s exact test were used for contingency tables. Each ‘N’ represents a single biological replicate and details for each experiment can be found in the figure legends. Data shown are representative of at least two independent experiments and are represented as the mean and error bars show the standard error of the mean (SEM) unless otherwise noted.
Results
Trem2Y38C/Y38C mouse model generation
To examine the consequences of a putative loss of function Trem2 variant, Trem2Y38C/Y38C mice containing the NHD Trem2 Y38C variant were generated using CRISPR/Cas9 to introduce a point mutation in the endogenous Trem2 gene, resulting in a substitution of tyrosine to cystine at position 38 (Fig. 1a). To investigate the effects of TREM2 deficiency, we used the CRISPR/Cas9 Trem2−/− mouse model. Trem2 expression in both models was evaluated using PCR primers flanking the inserted stop codon in Trem2−/− mice. Trem2 mRNA levels were unaffected in Trem2Y38C/Y38C mice, while no Trem2 mRNA expression was observed in Trem2−/− mice (Fig. 1b). Since an artifactual alternate splice variant has been reported in the Trem2 R47H variant mouse model that resulted in reduced Trem2 mRNA and protein expression [29], we evaluated the presence of alternate splice variants in the Trem2Y38C/Y38C mouse model but none were found (data not shown). Additionally, because the Velocigene Trem2−/− mouse model shows off-target upregulation of the adjacent Treml1 gene [30], we examined the expression of genes flanking the Trem2 locus, 100 kb upstream and downstream, on chromosome 17 in Trem2Y38C/Y38C and Trem2−/−. We did not detect significant aberrant expression of genes within this region (Additional file 3: Table S3). At the protein level, TREM2 was localized to Iba-1+ microglia in the brains of Trem2Y38C/Y38C mice and was absent in Trem2−/− mice (Fig. 1c). TREM2 protein levels in Trem2Y38C/Y38C cortical lysates, when measured by ELISA, were slightly higher than WT and were negligible in Trem2−/−, indicating that TREM2 protein is expressed and detected in Trem2Y38C/Y38C mice (Fig. 1d). sTREM2 levels in TBS soluble cortical lysates of Trem2Y38C/Y38C were similar to WT and were scant in Trem2−/− mice (Fig. 1e). This suggests, Trem2-Y38C mutation did not alter the expression of TREM2 protein in Trem2Y38C/Y38C mice.
Trem2Y38C/Y38C and Trem2−/− alters microglial morphology
To determine whether the Trem2- Y38C variant and deficiency of Trem2 affects microglia phenotype and number, we examined the total number of Iba-1+ cells and microglial morphology in the cortices of 6 months old WT, Trem2Y38C/Y38C and Trem2−/− mice using immunohistochemistry. The Iba-1+ cell number in both Trem2Y38C/Y38C or Trem2−/− mice were similar to WT (Fig. 1f). Trem2Y38C/Y38C and Trem2−/− microglia both showed increase in the percent area of Iba-1 staining indicating individual cells covered a larger area (Fig. 1g). This could be due to increased microglial branching in Trem2Y38C/Y38C and Trem2−/− mice as represented by the quantification of skeletonized microglia (Fig. 1h-j).
Trem2Y38C/Y38C and Trem2−/− mice show alterations in genes related to oligodendrocyte/myelin and neuronal functions
In order to better understand the effects of loss of Trem2 expression and the Trem2 Y38C variant on various physiological processes in the brain in an unbiased manner, we performed RNA sequencing to explore differences in gene expression in microdissected cortices of Trem2Y38C/Y38C, Trem2−/−, and WT mice at 6 months of age (Fig. 2a, Additional file 4: Figure S1A,B). The list of all differentially expressed genes (DEG) are in Additional files 5 and 6: Table S4 and S5. Compared to WT control mice, we observed 1076 DEG in Trem2Y38C/Y38C and 665 DEG in Trem2−/−. Cortices of Trem2Y38C/Y38C and Trem2−/− mice shared, 331 DEG in common (Additional file 4: Figure S1C). We performed gene ontology enrichment analysis on commonly altered genes using Enrichr. Enrichment for cell type indicated the downregulated genes associated with cortical layers 2/3 (Additional file 4: Figure S1D). We found 177 DEGs in Trem2Y38C/Y38C vs Trem2−/− cortices (Additional file 7: Table S6, Additional file 8: Figure S2A). To understand the pathways affected, genes (FDR < 0.05) were analyzed for pathway enrichment using Enrichr. Down regulated genes in Trem2Y38C/Y38C samples enriched in circadian rhythm pathways whereas upregulated genes were enriched in pathways associated with protein processing in endoplasmic reticulum as compared to the Trem2−/− samples (Additional file 8: Figure S2B-C).
Enrichment analysis of downregulated genes in Trem2Y38C/Y38C and Trem2−/− samples for cell compartments revealed myelin and axons as top enriched terms as compared to WT samples (Fig. 2c, Additional file 4: Figure S1E). Overall, this analysis highlighted alterations in oligodendrocytes/myelin and neuronal cells. To further characterize the effect of Trem2Y38C/Y38C and Trem2−/− on oligodendrocyte/myelin gene expression, we analyzed oligodendrocyte/myelin molecular signatures [31,32,33] to determine what effect Trem2 variant or loss of function would have on known oligodendrocyte/myelin gene expression patterns. Among these markers, we observed significant reduction in the expression of 17 transcripts in Trem2Y38C/Y38C and 6 transcripts in Trem2−/− (Additional file 9.10: Table S7 and S4). White-matter alterations have been reported in NHD patients, including patients with Trem2-Y38C homozygous mutations [10].
To explore the pathways that could elucidate the microglial and myelin alterations in Trem2Y38C/Y38C and Trem2−/−, we performed pathway analysis using Ingenuity Pathway Analysis (IPA) and Enrichr. As expected, microglial pathways associated with Trem2 including phagocytosis and Tyrob network were among the most affected pathways (Additional file 4: Figure S1F). RhoA signaling was commonly altered in Trem2Y38C/Y38C and Trem2−/− mice (Fig. 2d). Rho kinases are involved in pathways regulating cell shape, motility, cell survival and apoptosis. RhoA inhibits oligodendrocyte process extension thus regulating myelination [34]. RhoA/Rho (ROCK) signaling has also been implicated in dendritic spine morphology and axonal growth [35]. Abnormal upregulation of RhoA in Trem2Y38C/Y38C and Trem2−/− mice, thus, suggests impairments in myelination and synapses.
To further understand the gene expression changes, we performed IPA upstream regulator analysis. Based on the differentially expressed gene, this analysis predicts Upstream regulators that could be activated or inhibited resulting in the observed changes in gene expression. Our analysis indicated a predicted activation of TGFβ-1 upstream regulator in Trem2−/− versus WT mice (Fig. 2e). TGFβ-1 signaling has shown to be neuroprotective and increases after injury or stroke [36,37,38]. Additionally, we found MEF2c (Fig. 2e), a transcription factor that regulates synaptic remodeling in brains of adult mice [39,40,41], was predicted to be activated in both Trem2Y38C/Y38C and Trem2−/− mice. Together, these data along with the enrichment analysis provided support that Trem2 Y38C variant and Trem2 deficiency induce neuronal dysfunction at 6 months of age.
Trem2Y38C/Y38C and Trem2−/− adult mice display reduced myelination
To further study the oligodendrocyte/myelin alterations observed in the transcriptomics analysis we examined three oligodendrocyte/myelin associated genes, that were downregulated in both Trem2Y38C/Y38C and Trem2−/−, by qPCR. Myelin basic protein (Mbp) transcripts were significantly reduction in both Trem2Y38C/Y38C and Trem2−/− as compared to WT mice (Fig. 3a). Myelin Associated Oligodendrocyte Basic Protein (Mobp) and Oligodendrocyte Transcription Factor 2 (Olig2) transcripts decreased significantly in Trem2−/− mice as compared to the WT mice (Fig. 3b-c). We further studied oligodendrocyte/myelin pathology by analyzing protein expression of CNPase (2′,3′-Cyclic-nucleotide 3′-phosphodiesterase), an enzyme expressed in myelin associated oligodendrocytes, and MBP in cortical lysates of 6 months old WT, Trem2Y38C/Y38C and Trem2−/− mice. CNPase levels were reduced significantly in Trem2−/− as compared to the WT (Fig. 3d-e). MBP protein levels were lowered in both Trem2Y38C/Y38C and Trem2−/− as compared to WT mice (Fig. 3d-e). Overall, this suggests the presence of altered myelination in Trem2Y38C/Y38C and Trem2−/− brains.
Synaptic elements and synaptic plasticity are impaired in Trem2Y38C/Y38C and Trem2−/− young adults
Transcriptomics analysis also suggested alteration of genes associated with neuronal function. To validate and extend these finding in 6 months old Trem2Y38C/Y38C and Trem2−/− mice, we first evaluated synaptic protein levels in cortex and hippocampus using western blot. In the cortex, no changes were observed in postsynaptic protein PSD-95 levels but presynaptic protein synaptophysin levels were significantly reduced in both Trem2Y38C/Y38C and Trem2−/− mice (Fig. 4a,b). To confirm that the differences in presynaptic and postsynaptic proteins were not specific to these markers, we examined additional synaptic proteins. We observed no changes in Homer1 levels but noted significant decreases in synapsin 2 levels in cortices of Trem2−/− mice (Additional file 11: Figure S3A-B). The hippocampal protein lysates showed decreased pre- and postsynaptic elements with reduction in both PSD-95 and synaptophysin levels in Trem2Y38C/Y38C and Trem2−/− mice (Fig. 4c,d). Thus, region-specific changes in synaptic proteins were observed at 6 months with more alterations in hippocampal synapses.
To investigate if these alterations were initiated during development or were induced post-development, we evaluated synaptic elements in cortical and hippocampal lysates of WT, Trem2Y38C/Y38C, and Trem2−/− mice at P20. Neither group demonstrated alterations in synaptic protein levels at P20 (Additional file 11: Figure S3C-D). Additionally, we compared pre- and postsynaptic proteins in the hippocampi of P20 mice to that of 6 months old mice. PSD95 levels increased significantly with age in WT mice; while Trem2Y38C/Y38C and Trem2−/− mice didn’t display similar increase over the time, which resulted in reduced PSD95 levels in Trem2Y38C/Y38C and Trem2−/− as compared to WT mice at 6 months of age (Fig. 4e). This suggests that the decrease in PSD95 could be due to altered development of post-synaptic elements in adult Trem2Y38C/Y38C and Trem2−/− mice. Synaptophysin levels were maintained in WT hippocampi from P20 to 6 months of age whereas, 6 months old Trem2−/− mice displayed synaptophysin levels lower than their respective P20 levels. This suggests that absence of Trem2 might induce active loss of pre-synaptic elements in adult mice. Thus, this decreased abundance of synaptic proteins in Trem2Y38C/Y38C and Trem2−/− mice at 6 months may stem from a combination of altered synaptic development and loss of synaptic elements in adult mice. Next, we determined if the severe hippocampal alterations in synaptic elements correlated with altered microglial abundance and Iba1+ hippocampal area. Iba-1+ cell counts did not change significantly (Fig Additional file 11: Figure S3E), however, percentage of area occupied by Iba-1+ cells significantly increased in both Trem2Y38C/Y38C and Trem2−/− hippocampus (Fig Additional file 11: Figure S3F).
Lastly, to evaluate the impact of alterations of hippocampal synaptic protein on synaptic functions, we assessed induction of long-term potentiation (LTP) in hippocampal area CA1 by stimulating Schaffer collaterals in brain slices of 6 months old WT, Trem2Y38C/Y38C and Trem2−/− mice. No differences were detected in input/output curves or paired-pulse facilitation (Fig. 5a, Additional file 8: Figure S4). However, a statistically significant reduction in the magnitude of LTP was observed in Trem2Y38C/Y38C and Trem2−/− LTP responses as compared to WT (Fig. 5b-d). These results indicate severe impairment in hippocampal neuronal plasticity in Trem2Y38C/Y38C and Trem2−/− mice at 6 months.
Discussion
A common clinical manifestation in NHD patients and some of the genetic variants of TREM2 or TYROBP is early-onset dementia. To understand how loss of functional TREM2 confers predisposition to early-onset dementia, we used Trem2Y38C/Y38C and Trem2−/− mouse models to study neuronal pathology at a developmental and 6 months age. Importantly, as observed in the brains of humans harboring TREM2 variants [42], the endogenous expression levels of TREM2 Y38C transcript levels were comparable to the WT. This suggests there were no aberrant alterations in Trem2 mRNA, which was reported in Trem2 R47H mouse model [21, 29]. Similarly, no alterations were observed in protein expression of TREM2 and sTREM2 in Trem2Y38C/Y38C mice. Further studies are required to determine the proportion of TREM2 receptor expression at the plasma membrane retention to its retention in endoplasmic reticulum in vivo.
Since TREM2 is expressed exclusively on microglia in the brain, we analyzed microglia number and morphology in the cortices of Trem2Y38C/Y38C and Trem2−/− mice and found no significant changes in the number of Iba-1 positive microglia. However, we found increased area covered by Iba-1 positive microglia with increased number of branches which is in support of the previous findings [19, 20]. The increased microglial ramification of microglial processes suggests in Trem2Y38C/Y38C and Trem2−/− mice as compared to WT microglia is reminiscent of the “homeostatically locked” state that has been described for Trem2 deficient microglia [43]. Studying Trem2Y38C/Y38C and Trem2−/− microglia in presence of a pathological trigger such as amyloid can provide further information on activation capacity of microglia with mutant TREM2 or lack of TREM2.
Another manifestation in NHD patients with TREM2 dysfunction is development of FTD-like dementia with distinct white matter loss [10, 44]. Weighted gene co-expression network analysis of Velocigene Trem2−/− mice also showed oligodendrocyte/myelin related modules as one of the most disrupted cell type modules [45]. Likewise, Trem2Y38C/Y38C and Trem2−/− mice revealed alterations in oligodendrocytes/myelination based on the transcriptomics analysis. Since this suggested white matter dysfunction, we expanded our analysis of oligodendrocyte specific genes in Trem2Y38C/Y38C and Trem2−/− mice brains. Among the many downregulated oligodendrocyte/myelin genes, Mbp (myelin basic protein) and Mobp (myelin associated oligodendrocyte basic protein) were downregulated in both Trem2Y38C/Y38C and Trem2−/−. Expression of myelin proteins such as MBP and CNPase were reduced in Trem2Y38C/Y38C and Trem2−/−. IPA analysis revealed ILK signaling, which is involved in oligodendrocyte development [46], was impacted in both Trem2Y38C/Y38C and Trem2−/− mice. We observed significant reduction in MBP protein expression in the cortical lysates of Trem2Y38C/Y38C and Trem2−/− mice, further supporting our transcriptomics findings. Our observed oligodendrocytes/myelin alterations in cortices of the Trem2Y38C/Y38C and Trem2−/− young adult mice suggest nascent white-matter alterations in these mice. This is concordant with the reports highlighting the role of microglia and TREM2 in white matter physiology. Further detailed studies for analyzing white-matter tracts, either by brain scans or electron microscopy, in older mice are necessary to determine the severity and progression of the pathology.
Microglia support myelination by taking up myelin debris and stimulating oligodendrocyte maturation after demyelinating events [47, 48]. TREM2 contributes to this process by sensing lipids and phagocytosing myelin debris. In the absence of TREM2 in mouse models of demyelination, microglia take up myelin debris but fail to degrade it, resulting in the accumulation of lipids and persistent demyelination [49, 50]. With aging, Trem2 T66M mice exhibit reduced microglial activity and microglial clusters in white matter tracts [11]. Similarly, our findings highlight the role of microglial TREM2 in oligodendrocytes/myelin maintenance. Recently, using single-nuclei transcriptomics of human samples, AD-reactive oligodendrocytes were shown to increase significantly in TREM2-R62H carries as compared to the controls but seem to be lower than TREM2-common variant samples [51]. The study also reported downregulation of genes in Oligo1 cluster (includes genes such as Mbp, Mobp and Olig2) in early-onset AD which corresponds to our findings that indicate reduction in oligodendrocyte/myelin markers in Trem2Y38C/Y38C and Trem2−/− young adult mice.
We found perturbations in gene expression associated with neuronal functions in the cortex of adult Trem2Y38C/Y38C and Trem2−/− mice. IPA upstream regulator analysis was suggestive of synaptic remodeling as it predicted activation of upstream transcription factor Mef2c (Myocyte-specific enhancer factor 2C). Mef2c is regulated by synaptic activity and important in neuronal migration, synapse development and synaptic transmission [52]. Presynaptic proteins such as synaptophysin and synapsin, were reduced in the cortex and hippocampus, while levels of postsynaptic proteins, PSD95 and Homer1, were reduced only in the hippocampus. This suggests enhanced synapse-associated pathology in hippocampus of Trem2Y38C/Y38C and Trem2−/− mice at 6 months. These differences in regional specificity could be due to regional diversity observed in glial cells [53, 54]. To understand if the reduced synaptic protein levels were reflective of a developmental defect or a post-developmental phenomenon, we analyzed the synaptic markers in cortex and hippocampus at P20 and 6 months age. Previous reports have observed increased synaptic proteins in the hippocampus of Trem2−/− mice at P20 [19] or no changes were reported in 1 month old Trem2−/− mice [20]. Our data demonstrate no changes in synaptic protein expression at P20 in Trem2Y38C/Y38C and Trem2−/− mice. These results indicate that synaptic protein levels are not affected during development but decrease in adult mice. Analysis of synaptic proteins across P20 and 6 months ages suggested decreased abundance of synaptic proteins in Trem2Y38C/Y38C and Trem2−/− mice might be a combination of altered synaptic development and active loss of synaptic elements in adult mice. These synaptic alterations further influenced deficits in hippocampal synaptic plasticity as observed by significant reduction in hippocampal LTP. LTP at the Schaffer collateral-CA1 synapse is thought to be crucial for the cognitive process involved in learning and memory. Spatial memory deficits have also been reported in 12 months old Trem2−/− mice [55]. Trem2Y38C/Y38C and Trem2−/− adult mice exhibited abnormal behaviors such as excess self-grooming or barbering leading to hairless patches (data not shown and as reported previously [19]). However, more extensive battery of behavioral tests is needed to characterize impairments in behavior.
Pathway enrichment analysis revealed alterations in RhoA signaling in Trem2Y38C/Y38C and Trem2−/− mice. Rho signaling is involved in cytoskeleton remodeling and inhibition of RhoA triggers axonal growth cone [56]. RhoA inhibitors not only increase myelination in injury models but also restore synaptic plasticity in an anxiety mouse model [57, 58]. Altered Rho signaling and microglial activation could potentially provide molecular insights on underlying process the oligodendrocyte/myelin and synaptic dysfunction in Trem2Y38C/Y38C and Trem2−/− mice. However, further investigation is needed to understand if this pathway is altered in a specific or multiple cell type.
As previously reported [19, 20], our data indicate TREM2 is required to maintain proper synaptic balance. We expand these findings to demonstrate that loss of functional TREM2 negatively impacts neuronal synapses and oligodendrocytes/myelin post-developmentally. Interestingly, mice deficient in DAP12, a signaling adapter of various immunoreceptors including TREM2, also show hypomyelination and synaptic degeneration in thalamus as early as 3 months [59]. As dysmyelination has shown to alter synaptic transmission and in turn impact synapse integrity [60,61,62], we suspect that the alteration in synaptic integrity and transmission in absence of TREM2 could be a secondary effect of oligodendrocyte/myelin impairment due to reduced microglial activity. RhoA could potentially be the pathway orchestrating this pathology. Trem2Y38C/Y38C mice demonstrated pathological changes that were similar to Trem2−/− mice including several overlapping differentially expressed genes and significant changes in microglial morphology and synaptic protein levels. Functional outcomes were analogous since the reduction in LTP responses were comparable Trem2Y38C/Y38C and Trem2−/− mice. Recently, TREM2-R47H mutation has also shown to reduce LTP in young rats [63]. In normal human brains, Trem2 expression is higher in white matter and hippocampus as compared to other brain regions [64]. Interestingly, our findings display Trem2Y38C/Y38C and Trem2−/− mice have reduced myelination and synaptic alterations are most prominent in the hippocampus. Overall, Trem2Y38C/Y38C and Trem2−/− young adult mice present phenotypic characteristics reminiscent of the clinical manifestations in patients with NHD. The synaptic dysfunction along with oligodendrocyte/myelin impairment observed at 6 months in Trem2Y38C/Y38C and Trem2−/− mice could help explain why Trem2 variants are predisposed to early onset dementia. Collectively, the continuation of studies that dissect the mechanisms by which loss of functional TREM2 in microglia directly and indirectly, through microglial modulation of other cell types, impact synaptic functions and white matter could unveil mechanistic role of TREM2 in brain homeostasis and disease.
Conclusion
In summary, our findings contribute in vivo evidence that TREM2 Y38C disrupts normal TREM2 functions. Trem2Y38C/Y38C and Trem2−/− mice demonstrated altered gene expression, changes in microglia morphology, loss of synaptic proteins, reduced oligodendrocyte/myelin specific transcripts and impaired hippocampal synaptic plasticity. These results explain some of the early events leading to presenile dementia. Loss of TREM2 signaling negatively impacts neuronal function post-developmentally and offers insight to how TREM2 alterations can confer synaptic impairment. Future studies are needed to determine if these events precede neuronal susceptibility to pathological triggers such as amyloid.
Availability of data and materials
The datasets used and/or analyzed during the current study are provided as additional files and are available from the corresponding author on reasonable request.
Change history
20 November 2020
The original article has been revised to correct the presentation of Figure 4d.
Abbreviations
- NHD:
-
Nasu Hakola Disease
- FTD:
-
Fronto-Temporal Dementia
- AD:
-
Alzheimer’s disease
- CRISPR:
-
Clustered regularly interspaced short palindromic repeats
- PAM:
-
Protospacer adjacent motif
- Trem2 :
-
Triggering Receptor expressed on myeloid cells 2
- Iba-1:
-
Ionized calcium-binding adapter molecule 1
- RhoA:
-
Ras homolog family member A
- P20:
-
Postnatal day 20
- Tyrobp:
-
TYRO protein tyrosine kinase binding protein
- DAP12:
-
DNAX-Activation Protein 12
- LTP:
-
Long term potentiation
- PSD95:
-
Postsynaptic density protein 95
- Mef2c:
-
Myocyte Enhancer Factor 2C
- IPA:
-
Ingenuity pathway analysis
References
Hollingworth P, Harold D, Sims R, Gerrish A, Lambert J-C, Carrasquillo MM, et al. Common variants in ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer’s disease. Nat Genet. 2011;43:429–35.
Lambert J-C, Heath S, Even G, Campion D, Sleegers K, Hiltunen M, et al. Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer’s disease. Nat Genet. 2009;41:1094–9.
Naj AC, Jun G, Beecham GW, Wang L-S, Vardarajan BN, Buros J, et al. Common variants in MS4A4/MS4A6E, CD2uAP, CD33, and EPHA1 are associated with late-onset Alzheimer’s disease. Nat Genet. 2011;43:436–41.
Karch CM, Goate AM. Alzheimer’s disease risk genes and mechanisms of disease pathogenesis. Biol Psychiatry. 2015;77:43–51.
Distribution and signaling of TREM2/DAP12, the receptor system mutated in human polycystic lipomembraneous osteodysplasia with sclerosing leukoence... - PubMed - NCBI. Available from: https://www.ncbi.nlm.nih.gov/pubmed? Db=pubmed&Cmd=ShowDetailView&TermToSearch=15548205. [cited 2020 Feb 19].
Paloneva J, Manninen T, Christman G, Hovanes K, Mandelin J, Adolfsson R, et al. Mutations in two genes encoding different subunits of a receptor signaling complex result in an identical disease phenotype. Am J Hum Genet. 2002;71:656–62.
Jonsson T, Stefansson H, Steinberg S, Jonsdottir I, Jonsson PV, Snaedal J, et al. Variant of TREM2 associated with the risk of Alzheimer’s disease. N Engl J Med. 2013;368:107–16.
Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, et al. TREM2 variants in Alzheimer’s disease. N Engl J Med. 2013;368:117–27.
Kober DL, Alexander-Brett JM, Karch CM, Cruchaga C, Colonna M, Holtzman MJ, et al. Neurodegenerative disease mutations in TREM2 reveal a functional surface and distinct loss-of-function mechanisms. Elife. 2016;5:e20391.
Guerreiro RJ, Lohmann E, Brás JM, Gibbs JR, Rohrer JD, Gurunlian N, et al. Using exome sequencing to reveal mutations in TREM2 presenting as a Frontotemporal dementia–like syndrome without bone involvement. JAMA Neurol. 2013;70:78–84.
Kleinberger G, Brendel M, Mracsko E, Wefers B, Groeneweg L, Xiang X, et al. The FTD-like syndrome causing TREM2 T66M mutation impairs microglia function, brain perfusion, and glucose metabolism. EMBO J. 2017;36:1837–53.
Yao H, Coppola K, Schweig JE, Crawford F, Mullan M, Paris D. Distinct Signaling Pathways Regulate TREM2 Phagocytic and NFκB Antagonistic Activities. Front Cell Neurosci. 2019; 13. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6795686/. [cited 2020 Feb 21].
Kleinberger G, Yamanishi Y, Suárez-Calvet M, Czirr E, Lohmann E, Cuyvers E, et al. TREM2 mutations implicated in neurodegeneration impair cell surface transport and phagocytosis. Sci Transl Med. 2014;6:243ra86.
Sirkis DW, Aparicio RE, Schekman R. Neurodegeneration-associated mutant TREM2 proteins abortively cycle between the ER and ER-Golgi intermediate compartment. Mol Biol Cell. 2017;28:2723–33.
Guerreiro R, Bilgic B, Guven G, Brás J, Rohrer J, Lohmann E, et al. A novel compound heterozygous mutation in TREM2 found in a Turkish frontotemporal dementia-like family. Neurobiol Aging. 2013;34:2890.e1–5.
Park J-S, Ji IJ, An HJ, Kang M-J, Kang S-W, Kim D-H, et al. Disease-associated mutations of TREM2 Alter the processing of N-linked oligosaccharides in the Golgi apparatus. Traffic. 2015;16:510–8.
Schlepckow K, Kleinberger G, Fukumori A, Feederle R, Lichtenthaler SF, Steiner H, et al. An Alzheimer-associated TREM2 variant occurs at the ADAM cleavage site and affects shedding and phagocytic function. EMBO Mol Med. 2017;9:1356–65.
Yeh FL, Wang Y, Tom I, Gonzalez LC, Sheng M. TREM2 binds to Apolipoproteins, including APOE and CLU/APOJ, and thereby facilitates uptake of amyloid-Beta by microglia. Neuron. 2016;91:328–40.
Filipello F, Morini R, Corradini I, Zerbi V, Canzi A, Michalski B, et al. The Microglial Innate Immune Receptor TREM2 Is Required for Synapse Elimination and Normal Brain Connectivity. Immunity. 2018;48:979–991.e8.
Jay TR, von Saucken VE, Muñoz B, Codocedo JF, Atwood BK, Lamb BT, et al. TREM2 is required for microglial instruction of astrocytic synaptic engulfment in neurodevelopment. Glia. 2019;67:1873–92.
Cheng-Hathaway PJ, Reed-Geaghan EG, Jay TR, Casali BT, Bemiller SM, Puntambekar SS, et al. The Trem2 R47H variant confers loss-of-function-like phenotypes in Alzheimer’s disease. Mol Neurodegener. 2018;13:29.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Breese MR, Liu Y. NGSUtils: a software suite for analyzing and manipulating next-generation sequencing datasets. Bioinformatics. 2013;29:494–6.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44:W90–7.
Bemiller SM, McCray TJ, Allan K, Formica SV, Xu G, Wilson G, et al. TREM2 deficiency exacerbates tau pathology through dysregulated kinase signaling in a mouse model of tauopathy. Mol Neurodegener. 2017;12:74.
Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing - Benjamini - 1995 - Journal of the Royal Statistical Society: Series B (Methodological) - Wiley Online Library. Available from: https://rss.onlinelibrary.wiley.com/doi/abs/10.1111/j.2517-6161.1995.tb02031.x. [cited 2020 May 26].
Xiang X, Piers TM, Wefers B, Zhu K, Mallach A, Brunner B, et al. The Trem2 R47H Alzheimer’s risk variant impairs splicing and reduces Trem2 mRNA and protein in mice but not in humans. Mol Neurodegener. 2018;13:49.
Kang SS, Kurti A, Baker KE, Liu C-C, Colonna M, Ulrich JD, et al. Behavioral and transcriptomic analysis of Trem2-null mice: not all knockout mice are created equal. Hum Mol Genet. 2018;27:211–23.
Zhang Y, Sloan SA, Clarke LE, Caneda C, Plaza CA, Blumenthal PD, et al. Purification and characterization of progenitor and mature human astrocytes reveals transcriptional and functional differences with mouse. Neuron. 2016;89:37–53.
Avey D, Sankararaman S, Yim AKY, Barve R, Milbrandt J, Mitra RD. Single-Cell RNA-Seq Uncovers a Robust Transcriptional Response to Morphine by Glia. Cell Rep. 2018;24:3619–3629.e4.
Zhang Y, Chen K, Sloan SA, Bennett ML, Scholze AR, O’Keeffe S, et al. An RNA-sequencing Transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci Society Neurosci. 2014;34:11929–47.
Feltri ML, Suter U, Relvas JB. The function of RhoGTPases in axon ensheathment and myelination. Glia. 2008;56:1508–17.
Gross RE, Mei Q, Gutekunst C-A, Torre E. The pivotal role of RhoA GTPase in the molecular signaling of axon growth inhibition after CNS injury and targeted therapeutic strategies. Cell Transplant. 2007;16:245–62.
Buckwalter MS, Wyss-Coray T. Modelling neuroinflammatory phenotypes in vivo. J Neuroinflammation. 2004;1:10.
Brionne TC, Tesseur I, Masliah E, Wyss-Coray T. Loss of TGF-β1 leads to increased neuronal cell death and microgliosis in mouse brain. Neuron. 2003;40:1133–45.
Doyle KP, Cekanaviciute E, Mamer LE, Buckwalter MS. TGFβ signaling in the brain increases with aging and signals to astrocytes and innate immune cells in the weeks after stroke. J Neuroinflammation. 2010;7:62.
Harrington AJ, Raissi A, Rajkovich K, Berto S, Kumar J, Molinaro G, et al. MEF2C regulates cortical inhibitory and excitatory synapses and behaviors relevant to neurodevelopmental disorders. eLife. 5. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5094851/. [cited 2020 May 9].
Akhtar MW, Kim M-S, Adachi M, Morris MJ, Qi X, Richardson JA, et al. In vivo analysis of MEF2 transcription factors in synapse regulation and neuronal survival. PLoS One. 2012;7:e34863.
Tu S, Akhtar MW, Escorihuela RM, Amador-Arjona A, Swarup V, Parker J, et al. NitroSynapsin therapy for a mouse MEF2C haploinsufficiency model of human autism. Nat Commun. 2017;8:1488.
Del-Aguila JL, Benitez BA, Li Z, Dube U, Mihindukulasuriya KA, Budde JP, et al. TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers. Mol Neurodegener. 2019;14:18.
Mazaheri F, Snaidero N, Kleinberger G, Madore C, Daria A, Werner G, et al. TREM2 deficiency impairs chemotaxis and microglial responses to neuronal injury. EMBO Rep. 2017;18:1186–98.
Klünemann HH, Ridha BH, Magy L, Wherrett JR, Hemelsoet DM, Keen RW, et al. The genetic causes of basal ganglia calcification, dementia, and bone cysts: DAP12 and TREM2. Neurology. 2005;64:1502–7.
Carbajosa G, Malki K, Lawless N, Wang H, Ryder JW, Wozniak E, et al. Loss of Trem2 in microglia leads to widespread disruption of cell coexpression networks in mouse brain. Neurobiol Aging. 2018;69:151–66.
Hussain R, Macklin WB. Integrin-linked kinase (ILK) deletion disrupts Oligodendrocyte development by altering cell cycle. J Neurosci. 2017;37:397–412.
Erblich B, Zhu L, Etgen AM, Dobrenis K, Pollard JW. Absence of Colony Stimulation Factor-1 Receptor Results in Loss of Microglia, Disrupted Brain Development and Olfactory Deficits. Plos One. 2011;6:e26317.
Neumann H, Kotter MR, Franklin RJM. Debris clearance by microglia: an essential link between degeneration and regeneration. Brain. 2009;132:288–95.
Nugent AA, Lin K, B van L, Lianoglou S, Przybyla L, Davis SS, et al. TREM2 Regulates Microglial Cholesterol Metabolism upon Chronic Phagocytic Challenge. Neuron. 2020;105:837–854.e9.
Poliani PL, Wang Y, Fontana E, Robinette ML, Yamanishi Y, Gilfillan S, et al. TREM2 sustains microglial expansion during aging and response to demyelination. J Clin Invest. 2015;125:2161–70.
Zhou Y, Song WM, Andhey PS, Swain A, Levy T, Miller KR, et al. Human and mouse single-nucleus transcriptomics reveal TREM2-dependent and TREM2-independent cellular responses in Alzheimer’s disease. Nat Med. 2020;26:131–42.
Assali A, Harrington AJ, Cowan CW. Emerging roles for MEF2 in brain development and mental disorders. Curr Opin Neurobiol. 2019;59:49–58.
Tan Y-L, Yuan Y, Tian L. Microglial regional heterogeneity and its role in the brain. Mol Psychiatry. 2020;25:351–67.
Xin W, Bonci A. Functional astrocyte heterogeneity and implications for their role in shaping neurotransmission. Front Cell Neurosci. 2018;12:141.
Kim S-M, Mun B-R, Lee S-J, Joh Y, Lee H-Y, Ji K-Y, et al. TREM2 promotes Aβ phagocytosis by upregulating C/EBPα-dependent CD36 expression in microglia. Sci Rep. 2017;7 Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5593901/. [cited 2020 May 23].
Hall A, Lalli G. Rho and Ras GTPases in Axon Growth, Guidance, and Branching. Cold Spring Harb Perspect Biol. 2010;2 Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2828272/. [cited 2020 May 20].
Koch JC, Tatenhorst L, Roser A-E, Saal K-A, Tönges L, Lingor P. ROCK inhibition in models of neurodegeneration and its potential for clinical translation. Pharmacol Ther. 2018;189:1–21.
Deguchi Y, Harada M, Shinohara R, Lazarus M, Cherasse Y, Urade Y, et al. mDia and ROCK mediate actin-dependent presynaptic remodeling regulating synaptic efficacy and anxiety. Cell Rep. 2016;17:2405–17.
Kaifu T, Nakahara J, Inui M, Mishima K, Momiyama T, Kaji M, et al. Osteopetrosis and thalamic hypomyelinosis with synaptic degeneration in DAP12-deficient mice. J Clin Invest. 2003;111:323–32.
Wang F, Yang Y-J, Yang N, Chen X-J, Huang N-X, Zhang J, et al. Enhancing Oligodendrocyte Myelination Rescues Synaptic Loss and Improves Functional Recovery after Chronic Hypoxia. Neuron. 2018;99:689–701.e5.
Barron T, Saifetiarova J, Bhat MA, Kim JH. Myelination of Purkinje axons is critical for resilient synaptic transmission in the deep cerebellar nucleus. Sci Rep. 2018;8:1022.
Papuć E, Rejdak K. The role of myelin damage in Alzheimer’s disease pathology. Arch Med Sci. 2018;16:345–51.
Ren S, Yao W, Tambini MD, Yin T, Norris KA, D’Adamio L. In: MM MC, Rath S, Levy E, Gan L, editors. Microglia TREM2R47H Alzheimer-linked variant enhances excitatory transmission and reduces LTP via increased TNF-α levels, vol. 9: eLife Sciences Publications, Ltd; 2020. p. e57513.
Forabosco P, Ramasamy A, Trabzuni D, Walker R, Smith C, Bras J, et al. Insights into TREM2 biology by network analysis of human brain gene expression data. Neurobiol Aging. 2013;34:2699–714.
Acknowledgements
This work was supported by grants from the Alzheimer’s Association (BTL and GEL); BrightFocus Foundation A2015296S (BTL); NIH RF1 AG051495 (BTL and GEL); and NIH U54 AG054345 (BTL). In addition, this study was supported by generous donations from the Jane and Lee Seidman Fund, Chet and Jane Scholtz, and Dave and Susan Roberts. We would like to acknowledge Center for Medical Genomics at Indiana University School of Medicine and Indiana Precision Health Initiative program for conducting sequencing experiments and analysis. We thank McCray TJ and Reddy SV for technical assistance; Codocedo JF for helpful advice; Cheng-Hathaway PJ for contributing in Trem2 Y38C model generation; Dr. Meyer J and Huang KC for their guidance on assessing CRIPSR/Cas9 off-target effects; laboratories of Drs. Kim J and Lasagna-Reeves C for sharing reagents; Karandikar RH for his assistance with coding.
Author information
Authors and Affiliations
Contributions
Conceptualization, VSJ, BTL, GEL, ALO; Methodology, VSJ, BTL, GEL, ALO, SJB, SSP and MM; Formal Analysis, VSJ, BTL, GEL, ALO, SJB, SSP; RNA sequencing analysis, JZ, AJ, VSJ; Electrophysiology experiments, BKA, TP, GVDP; Investigation, VSJ, BTL, GEL, PBCL, GX; Writing – Original Draft, VSJ, SJB, ALO; Writing – Review & Editing, VJ, BTL, GEL, SJB, ALO, SSP, MM, BKA, JZ; Funding acquisition, BTL, GEL; Resources, ALO and GX. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Mice were housed in the Indiana University School of Medicine (IUSM) animal facilities, which are accredited by the Association and Accreditation of Laboratory Animal Care. Animals were maintained according to USDA standards and the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Experiments were approved by the IUSM Institutional Animal Care and Use Committee.
Consent for publication
Consent for publication on animals used in this study are detailed in the Material Transfer Agreement (MTA) between organizations.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Table S1.
Primers for amplification of genomic regions predicted to have off-target mutations from CRISPR/Cas9.
Additional file 2: Table S2.
Primers for Sanger’s sequencing of predicted off-target mutations from CRISPR/Cas9.
Additional file 3: Table S3.
Trem2Y38C/Y38C mice do not display aberrant expression of genes upstream or downstream of Trem2. Log fold change (logFC) with their FDRs are shown for genes within 100 kb upstream and downstream of Trem2. Results are normalized to gene expression in WT mice. N/A = not available.
Additional file 4: Figure S1.
(A,B) Heatmap representations for differentially expressed genes in (A) Trem2Y38C/Y38C versus WT mice and (B) Trem2−/− versus WT mice. (C) Venn diagram showing the overlap of significantly upregulated (blue) and downregulated (gray) differentially expressed genes between Trem2Y38C/Y38C and Trem2−/− mice compared to wildtype mice . (D-G) Using Enrichr web- tools, downregulated genes shared between Trem2Y38C/Y38C and Trem2−/− mice were analyzed for tissue enrichment (D) enrichment of genes highlighting association of downregulated genes with cortical layers. Cellular compartment enrichment includes (E) myelin sheath, (F) Pathway enrichment showing alterations in microglial phagocytosis-associated pathway. (D-F) P values are indicated at the end of each colored bar. Shaded gets lighter with smaller P value. Sample sizes: WT mice, N = 8 (4 males, 4 females); Trem2Y38C/Y38C mice, N = 8 (4 males, 4 females); Trem2−/− mice, N = 5 (1 male, 4 females).
Additional file 5: Table S4.
CSV file. List of differentially expressed genes in cortices of Trem2Y38C/Y38C versus WT mice. N=8 (equal number of males and females).
Additional file 6: Table S5.
CSV file. List of differentially expressed genes in cortices of Trem2-/- versus WT mice mice. N=5 (1 male, 4 females).
Additional file 7: Table S6.
CSV file. List of differentially expressed genes in cortices of Trem2Y38C/Y38C versus Trem2-/- mice. N=5 (1 male, 4 females).
Additional file 8: Figure S2.
(A) Volcano plot illustrating differentially expressed genes (P<0.05) in Trem2Y38C/Y38C versus Trem2-/- cortices. (B-C) Pathway enrichment of downregulated and upregulated genes (FDR<0.05) by Enrichr. (B) Pathway enrichment for downregulated genes indicates alteration in circadian rhythm associated pathways in Trem2Y38C/Y38C as compared to Trem2-/- mice (C) Upregulated genes show enrichment in pathways associated with protein processing in endoplasmic reticulum in Trem2Y38C/Y38C versus Trem2-/-. (B-C) P values are indicated at the end of each colored bar. Shaded gets lighter with smaller P value. Red bars represent enrichment for downregulated genes and green bars represent enrichment for upregulated genes. Sample sizes: WT mice, N = 8 (4 males, 4 females); Trem2Y38C/Y38C mice, N = 8 (4 males, 4 females); Trem2-/- mice, N = 5 (1 male, 4 females).
Additional file 9: Table S7.
Downregulated genes associated with oligodendrocyte/myelin in Trem2Y38C/Y38C versus WT mice. logFC = Log fold change.
Additional file 10: Table S8.
Downregulated genes associated with oligodendrocyte/myelin genes Trem2-/- versus WT mice. logFC = Log fold change.
Additional file 11: Figure S3.
(A,B) Western blot analysis was performed on cortical lysates of 6 months old WT, Trem2Y38C/Y38C and Trem2-/- mice. (B) Quantification of western blot bands shows decreased synapsin 2 levels in Trem2-/- mice with no changes in Homer1 levels. (C,D) Western blots and quantification of PSD95 and synaptophysin in (C) cortical and (D) hippocampal lysates of P20 WT, Trem2Y38C/Y38C and Trem2-/- mice. (E) Iba-1+ cell count in the hippocampus of 6 months WT, Trem2Y38C/Y38C and Trem2-/- mice showed no significant change (F) Iba-1 percent area in the hippocampus increased significantly in both Trem2Y38C/Y38C and Trem2-/- as compared to the WT. (B-D) Data represented as mean ± SEM. Sample size: N = 8 (4 males and 4 females). One-way ANOVA with Tukey’s post hoc test. *** p< 0.001. (E-F).
Additional file 12: Figure S4.
WT, Trem2Y38C/Y38C and Trem2−/− mice exhibit similar paired-pulse ratio (PPR). Paired-pulse ratio is the amplitude of excitatory postsynaptic current (EPSC) 2/EPSC1. ISI = interstimulus interval. Data are shown as mean ± SEM. Sample size: WT mice, N = 15 (7 females and 8 males), n = 33 recordings; Trem2Y38C/Y38C mice, N = 7 (4 males and 3 females), n = 20 recordings; Trem2−/− mice, N = 8 (4 males and 4 females), n = 15 recordings.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jadhav, V.S., Lin, P.B.C., Pennington, T. et al. Trem2 Y38C mutation and loss of Trem2 impairs neuronal synapses in adult mice. Mol Neurodegeneration 15, 62 (2020). https://doi.org/10.1186/s13024-020-00409-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13024-020-00409-0