In-situ glial cell-surface proteomics identifies pro-longevity factors in Drosophila
Figures
Glial cell-surface proteomics of the intact brains of young and old flies.
(A) Transgenic expression of horseradish peroxidase (HRP) fused to the transmembrane protein CD2 targets CD2-HRP to the extracellular side of the glial plasma membrane. In the presence of H2O2, CD2-HRP converts BxxP substrate into phenoxyl radicals that biotinylate nearby endogenous proteins. Biotinylated proteins are then enriched from brain lysates using streptavidin beads, followed by on-bead trypsin digestion, tandem mass tag (TMT) labeling, and liquid chromatography-tandem mass spectrometry analysis (LC-MS/MS). (B) To confirm HRP-dependent biotinylation of glial cells, Repo-GAL4 were crossed with either UAS-HRP-CD2 or wild-type (w1118) flies and adult brains were stained with neutravidin. Proteins were biotinylated in the experimental group (UAS-HRP-CD2) but not in the control group (w1118). (C) Biochemical characterization by streptavidin blot showed that the glial cell-surface-targeted HRP biotinylated a wide range of proteins compared to the negative control lacking HRP-CD2 expression. (D) For each time point (5 days and 50 days), two biological replicates and two negative controls (one lacking HRP-CD2 expression and one lacking H2O2 in the reaction) were profiled, for a total of eight samples. Each sample consisted of 300 mixed male and female fly brains, for a total of 2400 brains. In total, we detected 2007 proteins with two or more unique peptides. Potential contaminants were removed by comparing differential enrichment between the experimental and negative control groups (i.e. the TMT ratio for each protein identified; see Methods). In the TMT labeling system, 126, 127 N, 127 C, 128 N, 128 C etc. were used to index and compare different samples. (E) Strong correlations between biological replicates were observed, further validating the quality of the proteome profiles. (F) The receiver operating characteristic (ROC) curve depicts the true-positive rate against the false-positive rate of detected proteins, for each biological replicate (only 126:127 C pictured). (G) Gene Ontology (GO) analysis was performed based on the 872 proteins identified following ratiometric and cutoff (see Methods). Top five GO terms on cellular compartments with the lowest false discovery rates are shown. See Figure 1—figure supplement 1.
Glial cell-surface proteomics.
(A) We searched the 872 proteins included in our proteomic profiles against three databases: FlyBase, UniProt, and the Gene Ontology Data Archive (DOI 10.5281/zenodo.1205166; Version 2023-05-10; Ashburner et al., 2000) to determine which proteins had previous plasma membrane or extracellular matrix annotations. We found that 682 out of the 872 proteins could be validated by at least one database, further validating the quality of our cell-surface proteomic labeling. (B) We plotted the receiver operating characteristic (ROC) curve for each biological replicate and observed high specificity. (C) We identified the top 10 ‘hub’ genes in our protein-protein interaction (PPI) network and found that of the top 10 hub genes identified, seven had reported function in cell adhesion (klg, kst) and synaptic wiring (DIP-β, DIP-δ, DIP-γ, DIP-κ, DIP-ɑ), while the remaining three had reported function in neuron projection (CG33543), axon guidance (CG34353), and dendrite guidance (Toll-6).
Quantitative comparison of young and old proteomes and selection of candidate genes.
(A) To select approximately 50 candidates for further screening, we generated a volcano plot, where the tandem mass tag (TMT) ratio (log2FoldChange) for each of the 872 proteins identified was plotted on the x-axis and the log10(Ratio P-value) was plotted on the y-axis. Genes associated with the largest change in protein levels between time points were identified using an arbitrary cutoff targeting about 50 genes for screening. In total, 48 candidate genes were selected for further screening, most of which have at least one strong human ortholog. Down-regulated candidate genes selected for functional screening are labeled in blue whereas up-regulated candidate genes are labeled in red. (B) A Gene Ontology (GO) analysis investigating proteome biological processes revealed that genes encoding the most up-regulated proteins were primarily associated with cell localization, transport, and maintaining homeostasis, while most down-regulated proteins were primarily associated with neural development, synapse organization, and neuroplasticity. (C) Genes encoding the most up- and down-regulated proteins were clustered by their reported protein-protein interactions (PPIs) and corresponding confidence scores (using a Markov clustering algorithm, inflation value set to 2.5). Each protein in the PPI plot was then assessed individually to determine if it had reported function in synapse organization, synapse regulation, cell adhesion, development and cell migration, or localization and transport.
Functional screening of candidate genes predicted to be involved in normal brain aging.
(A) Confocal images showing that the glial-GeneSwitch (GS) line only induced transgene expression in the presence of the inducible drug Ru486. (B) Forty eight candidate genes were selected for functional screening. For each of the 19 up-regulated genes, two RNAi lines were used to control for possible off-target effects. Of the 29 down-regulated genes selected, we were able to acquire and screen lines for 19 genes. (C) No obvious lifespan effect was observed in control lines: a w1118 wild-type control, a UAS-luciferase RNAi control, a UAS-lacZ control, and a flySAM control. (D) Pie charts summarizing the results of our broad-scale lifespan screen for male and female flies. (E) Of the 19 up-regulated candidate genes screened, RNAi knockdown of four genes led to consistent lifespan decreases for males. Alternatively, RNAi knockdown of three genes led to consistent lifespan increases for females. (F) Of the 19 down-regulated candidate genes screened, overexpression of two genes led to lifespan increases for males, whereas overexpression of three genes led to lifespan decreases. In parallel, overexpression of 10 genes led to lifespan increases for females, whereas overexpression of one gene led to decreased lifespan. (G) Summary of lifespan effect for both up- and down-regulated genes. See Figure 3—figure supplements 1–7.
Functional screening of candidate genes.
(A–B) In our initial functional lifespan screen, glial overexpression of DIP-β was found to significantly extend lifespan for both male and female flies, compared to the negative controls (only the luciferase RNAi and flySAM controls are pictured here).
Functional screening of up-regulated candidate genes.
(A–E) We selected 19 up-regulated candidate genes for functional screening, which we selectively knocked down in glia to examine their influence on lifespan. For each up-regulated gene, we screened two RNAi lines (to control for possible off-target effects). Of the 19 up-regulated candidates screened, we found that RNAi knockdown of four genes led to consistent lifespan decreases for males. Alternatively, for female flies, RNAi knockdown of three genes led to consistent lifespan increases. The result of the first five genes are shown.
Functional screening of up-regulated candidate genes.
(A–E) We selected 19 up-regulated candidate genes for functional screening, which we selectively knocked down in glia to examine their influence on lifespan. For each up-regulated gene, we screened two RNAi lines (to control for possible off-target effects). Of the 19 up-regulated candidates screened, we found that RNAi knockdown of four genes led to consistent lifespan decreases for males. Alternatively, for female flies, RNAi knockdown of three genes led to consistent lifespan increases. The result of the second five genes is shown.
Functional screening of up-regulated candidate genes.
(A–E) We selected 19 up-regulated candidate genes for functional screening, which we selectively knocked down in glia to examine their influence on lifespan. For each up-regulated gene, we screened two RNAi lines (to control for possible off-target effects). Of the 19 up-regulated candidates screened, we found that RNAi knockdown of four genes led to consistent lifespan decreases for males. Alternatively, for female flies, RNAi knockdown of three genes led to consistent lifespan increases. The result of the third five genes is shown.
Functional screening of up-regulated candidate genes.
(A–E) We selected 19 up-regulated candidate genes for functional screening, which we selectively knocked down in glia to examine their influence on lifespan. For each up-regulated gene, we screened two RNAi lines (to control for possible off-target effects). Of the 19 up-regulated candidates screened, we found that RNAi knockdown of four genes led to consistent lifespan decreases for males. Alternatively, for female flies, RNAi knockdown of three genes led to consistent lifespan increases. The result of the last four genes is shown.
Functional screening of down-regulated candidate genes.
(A–J) We selected 29 down-regulated candidate genes for functional screening, which we selectively overexpressed in glia to examine their influence on lifespan. For each of the 29 down-regulated genes, we screened one overexpression line (using either UAS-cDNA lines or flySAM2.0 lines), when available. In total, we were able to acquire and screen lines for 19 genes. Of these candidates, overexpression of two genes led to lifespan increases for males, whereas overexpression of three genes led to lifespan decreases. Alternatively, overexpression of 10 genes led to lifespan increases for females, whereas overexpression of one gene led to decreased lifespan. We found that constitutive overexpression of DIP-β in glia significantly increased lifespan in both males and females compared to controls. The result of the first 10 genes are shown.
Functional screening of down-regulated candidate genes.
(A–I) We selected 29 down-regulated candidate genes for functional screening, which we selectively overexpressed in glia to examine their influence on lifespan. For each of the 29 down-regulated genes, we screened one overexpression line (using either UAS-cDNA lines or flySAM2.0 lines), when available. In total, we were able to acquire and screen lines for 19 genes. Of these candidates, overexpression of two genes led to lifespan increases for males, whereas overexpression of three genes led to lifespan decreases. Alternatively, overexpression of 10 genes led to lifespan increases for females, whereas overexpression of one gene led to decreased lifespan. We found that constitutive overexpression of DIP-β in glia significantly increased lifespan in both males and females compared to controls. The result of the remaining nine genes is shown.
Validating that pan-glial overexpression of DIP-β increases lifespan and serves a protective function during aging.
(A–D) Glial overexpression of DIP-β (using both a flySAM2.0 and a UAS overexpression line) significantly extends lifespan for both male and female flies, compared to the negative controls (only the w1118 and UAS-lacZ controls are shown here). (E) To further validate that glial overexpression of DIP-β is serving a protective function during brain aging, climbing (i.e. negative geotaxis) assays were conducted. Glial overexpression of DIP-β was also associated with improved climbing and motor control ability in aged flies.
Single-nucleus RNA-sequencing - investigating the downstream transcriptomic effects of glial overexpression of DIP-β during aging.
(A) Two whole-head samples (20 heads each) were prepared for snRNA-seq: one where DIP-β had been overexpressed in glia from day 5–50 and one control. Single-nucleus RNA-seq libraries were prepared using the 10 X Genomics Chromium Single Cell 3’ HT v3.1 protocol. (B) In total, approximately 49,000 high-quality cells were obtained after quality control and filtering from our fly head samples, which consist of 8 broad cell types. (C) Broad cell-type cluster compositions were largely consistent between the glial DIP-β overexpression treatment (Ru486+) and the control (Ru486-) groups. (D) A broad analysis of differentially expressed genes (DEGs) revealed that the glial cell-type cluster showed the second highest number of DEGs between the glial DIP-β overexpression treatment (Ru486+) and the negative control (Ru486-) groups. (E) A more detailed analysis investigating DEGs across 90 different cell subtypes revealed that of the glial cell subtypes, cortex glia (i.e. ‘adult brain cell body glia’) and ensheathing glia showed the highest number of DEGs between the DIP-β overexpression treatment (Ru486+) and the control (Ru486-) groups. (F) Of the 252 DEGs identified for glial cells, 186 were unique to cortex glia, 27 were unique to ensheathing glia, and 40 were found to be differentially expressed in both cortex and ensheathing glia. (G) As fat cells exhibited the highest number of DEGs between the DIP-β overexpression and control groups, we conducted a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to evaluate the biological functions of these genes. Most fat cell DEGs were up-regulated, and the top altered metabolic and signaling pathways for these genes were found to be related to homeostasis, signal integration, stress resistance, and longevity. (H) Cell-cell communication (CCC) analysis was performed using FlyPhoneDB2. DIP-β overexpression was generally associated with an increase in cell-cell communication, particularly glia-neuron and glia-fat cell communication. The top altered glial signaling pathways are shown. (I) Our proposed model: Glial overexpression of DIP-β leads to improved glia-neuron and glia-fat cell communication as flies age, which may contribute to the observed lifespan extensions. See Figure 5—figure supplement 1.
Single-nucleus RNA-sequencing - investigating the downstream transcriptomic effects of pan-glial DIP-β overexpression.
(A) In total, we sequenced approximately 49,000 cells from our fly head samples, with approximately 23,000 from the glial DIP-β overexpression (i.e. Ru486+) sample and approximately 26,000 from the negative control (i.e. Ru486-) sample. (B) Detailed cell subtype cluster composition annotations. (C) A detailed analysis investigating differentially expressed genes (DEGs) across cell subtypes.
Tables
| Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
|---|---|---|---|---|
| Strain, strain background (Drosophila) | UAS-HRP-CD2 | Liqun Luo lab | Li et al., 2020 | |
| Strain, strain background (Drosophila) | GSG3285-1 (pan-glial GeneSwitch) | Linda Partridge lab | BDSC 59948 | Nicholson et al., 2008 |
| Strain, strain background (Drosophila) | DIP-β flySAM | Norbert Perrimon lab | BDSC 98929 | Jia et al., 2018 |
| Strain, strain background (Drosophila) | UAS-DIP-β | Chris Doe lab | ||
| Chemical compound, drug | Ru486 drug (Mifepristone) | Cayman Chemical | 10006317 | 4 mM concentration |
| Chemical compound, drug | Neutravidin | Thermo Fisher | 31000 | 1:250 dilution |
| Antibody | chicken polyclonal, anti-GFP | Aves Labs | NC9510598 | 1:1000 dilution |
| Antibody | mouse monoclonal, anti-Repo | DSHB | 8D12 anti-Repo | 1:100 dilution |
| Antibody | rat monoclonal, anti-Ncad | DSHB | DN-Ex #8 | 1:40 dilution |
| Antibody | Rabbit polyclonal, anti-DsRed | Takara Bio | 632496 | 1:200 dilution |
| Antibody | Chicken Alexa Fluor 488, secondary antibody | Jackson ImmunoResearch | 703-545-155 | 1:250 dilution |
| Antibody | Mouse Alexa Fluor 568, secondary antibody | Invitrogen | A10037 | 1:250 dilution |
| Antibody | Rat Alexa Fluor 647, secondary antibody | Jackson ImmunoResearch | 112-605-167 | 1:250 dilution |
| other | Slowfade mount buffer | Invitrogen | S36936 | See Methods/ Immunohistochemistry |
| Commercial assay, kit | Pierce Streptavidin Magnetic Beads | Thermo Fisher | 88817 | NA |
| Commercial assay, kit | TMT10plex Isobaric Label Reagents and Kits | Thermo Scientific | 90110 | NA |
Additional files
-
Supplementary file 1
Original proteomic data, candidate genes and fly lines.
(A) Original proteomic data prior to analyses. (B) Candidate genes selected for further investigation. (a) In total, we selected 29 down-regulated candidate genes for functional screening. Candidate genes were selected for further screening only if they were known (or predicted) to be located either on the cell surface or in the extracellular space. (b) In total, we selected 19 up-regulated candidate genes for functional screening. Candidate genes were selected for further screening only if they were known (or predicted) to be located either on the cell surface or extracellular space. Most candidate genes selected for further screening were found to have at least one strong human ortholog (Chen et al., 2017; Wang et al., 2019a; Wang et al., 2019b). (C) RNAi lines used to screen up-regulated candidate genes and OE lines used to screen down-regulated candidate genes. For each of the 19 up-regulated candidate genes, we screened two RNAi transgenic lines to control for potential off-target effects. For each of the 29 down-regulated genes, we screened one overexpression line (using either UAS-cDNA lines or flySAM2.0 lines), when available. Of the 29 down-regulated genes selected, we were able to acquire and screen lines for 19 genes.
- https://cdn.elifesciences.org/articles/109422/elife-109422-supp1-v1.xlsx
-
MDAR checklist
- https://cdn.elifesciences.org/articles/109422/elife-109422-mdarchecklist1-v1.docx