Synaptic deregulation of cholinergic projection neurons causes olfactory dysfunction across five fly Parkinsonism models

  1. Ulrike Pech
  2. Jasper Janssens
  3. Nils Schoovaerts
  4. Sabine Kuenen
  5. Carles Calatayud Aristoy
  6. Sandra F Gallego
  7. Samira Makhzami
  8. Gert J Hulselmans
  9. Suresh Poovathingal
  10. Kristofer Davie
  11. Adekunle T Bademosi
  12. Jef Swerts
  13. Sven Vilain
  14. Stein Aerts
  15. Patrik Verstreken  Is a corresponding author
  1. VIB-KU Leuven Center for Brain and Disease Research, Belgium
  2. KU Leuven, Department of Neurosciences, Leuven Brain Institute, Belgium
  3. KU Leuven, Department of Human Genetics, Leuven Brain Institute, Belgium
  4. VIB-KU Leuven Center for AI and Computational Biology (VIB.AI), Belgium
  5. VIB-KU Leuven Center for Brain and Disease Research Technologies, Single Cell, Microfluidics and Bioinformatics Expertise Units, Belgium
7 figures, 1 table and 8 additional files

Figures

Figure 1 with 1 supplement
A collection of familial Parkinson’s disease (PD) knock-in models.

(a) Scheme of the knock-in strategy where the first common exon in all Drosophila transcripts (left) was replaced by an attP-flanked mini-white gene using CRISPR/Cas9 mediated homologous recombination, creating a null mutant, and then replaced by wild-type or pathogenic mutant human or fly coding DNA sequence (CDS) using PhiC31 mediated integration (right). The chromosomal positions of the genes are indicated. All knock-ins are in the same isogenic genetic background (cantonized w1118). (b) Phenotypic analysis of PD mutants: (1) Mitochondrial membrane potential measured by ratiometric tetramethyl rhodamine ethyl ester (TMRE) fluorescence at neuromuscular junction (NMJ) boutons in third-instar larvae. n≥4 animals per genotype and 10 boutons from ≥3 NMJs per animal. (2) Depolarization amplitude quantified from electroretinograms (ERGs), recorded after 7 days of light exposure. n≥20 animals per genotype. (3) Quantification of startle-induced negative geotaxis (SING) at 5±1 days after eclosion (young) and 25±1 days after eclosion (old). 95% of pink1 mutants died <25 days and were tested at 15 days. Values are normalized to control (see Materials and methods). Variance of control measurements are in gray. Bars are mean ± s.e.m.; *, p<0.05 ANOVA/Dunnett’s test.

Figure 1—figure supplement 1
Characterization of Parkinson’s disease (PD) knock-in models.

(a) RNA expression levels, measured by RT-PCR with primers to the indicated coding DNA sequence (CDS) (above the graphs) in heads of controls and animals with wild-type or pathogenic mutant CDS knocked-in. Data are expressed as percent of endogenous Drosophila gene expression. Gray bars indicate the variance in controls. Bars are mean ± s.e.m. (b) Example images of mitochondrial tetramethyl rhodamine ethyl ester (TMRE) fluorescence at larval neuromuscular junctions (NMJs). Insets: magnified synaptic boutons (see also Figure 1). Scale bar: 10 μm. (c) Example electroretinogram (ERG) recordings (1 s light stimulus); see also Figure 1. Animals were kept for 7 days in constant light prior to the experiment. Red line: the average of the individual raw data traces that are shown.

Figure 2 with 1 supplement
Single-cell RNA sequencing reveals that olfactory projection neurons (OPNs) are impaired across young Parkinson’s disease (PD) fly models.

(a) tSNE of 226k cells characterized by single-cell RNA sequencing of entire young (5-day-old) fly brains (controls and PD knock-in mutants). Colors indicate the cell types; 186 of which were identified (see Materials and methods). Key cell types are encircled, including dopaminergic neurons (DAN, blue), mushroom body Kenyon cells (orange), and OPN (green). (b–f) tSNE of the cells from five selective PD knock-in mutants (5 days of age). Black cells are those with a significant transcriptomic change. Key cell types labeled in (a) are indicated. Note that OPN are consistently affected across mutants, while DAN are not at this early stage.

Figure 2—figure supplement 1
Validation and analysis of single-cell RNA sequencing data.

(a) Frequency of individual cell types (olfactory projection neurons [OPNs], dopaminergic neurons [DAN], visual projection neurons [T1]), and the sex ratios of the sequenced cells detected across the Parkinson’s disease (PD) knock-in models and controls. Note, we only sequenced males from hPINK1P399L mutants. Error bars represent mean±95% CI. (b) Alignment of published transcriptomic profiles defined by bulk sequencing of FAC sorted cell types (y-axis, see also Davie et al., 2018) and the cell types identified in this study (x-axis), validating the identification of cell types in the tSNE clusters. Green indicates central brain cell types; blue is optic lobe; yellow: glia; gray: others. (c) Expression level of knocked-in coding DNA sequence (CDS) across the cell types, indicating the knocked-in genes are broadly expressed. DAN - dopaminergic neurons (PAM cluster), OPNs - olfactory projection neurons, T1 - T1 visual interneurons. (d) Comparison of two different algorithms to detect cell-type-specific differentially expressed gene (DEG) (see Materials and methods) showing a similar number of DEG and high correlation of signed p-values. (e) Graphs showing the number of DEG versus the number of cells in a cluster (each dot is a cell type) for each PD fly model and the model of cell number dependency of p-values (dotted line, see Materials and methods). Confidence intervals in red. Cell types with residuals above confidence intervals are considered affected.

Figure 3 with 1 supplement
Common differentially expressed genes (DEGs) in cholinergic neuron-rich brain regions of human Parkinson’s disease (PD) patients and PD fly models.

(a–a”) Schematic of the sunburst plot indicating Gene Ontology (GO) terms for each sector (a) and the mapping of the DEGs in nucleus basalis of Meynert (NBM), nucleus accumbens, and putamen brain samples idiopathic PD patients (with LRRK2 risk mutations) and controls (a’) and mapping of the DEGs found commonly in fly and human samples (a”). Inner rings represent the different GO categories (indicated in a), with their subcategories in the outer rings, rings (in a’–a”) are color-coded according to enrichment Q-value. (b) GO analysis of DEG in cholinergic neurons of young PD fly models and NBM neurons of PD patients. Redundant terms were removed. Color: adjusted p-value. (c) Schematic of the DEGs found commonly in fly PD models (blue) and human PD samples (black) manually sorted according to their previously described synaptic functions. Supplementary file 7 contains the summarized results of the DEG analysis of fly brains and postmortem human brain samples and the SynGO analysis.

Figure 3—figure supplement 1
Single-cell RNA sequencing of cholinergic regions in postmortem human brain samples.

(a) tSNE of cells characterized in three cholinergic brain regions of 10 brains, annotated are the major cell types (left) and 36 individually identified cell types (numbered, right). N - neurons, OG - oligodendrocytes, ENT - endothelial cells, EPY - ependymal cells, LM - lymphocyte, MG - microglia, AST - astrocytes, OPC - oligodendrocyte precursor cells, U - unknown. (b) Frequency of the 36 cell types in pooled Parkinson’s disease (PD) versus non-PD samples, in individual brain samples, and in the three brain regions. (c) tSNE, as shown in (a), with different neuronal marker gene expression indicated (counts of unique molecular identifier [UMI]). (d) Frequency of cells expressing cholinergic marker genes ACHE and CHAT. NBM - nucleus basalis Meynert, Nac - nucleus accumbens, Put - putamen.

Synaptic defects in olfactory projection neuron (OPN) of Parkinson’s disease (PD) mutants.

(a) Schematic of a head-fixed awake fly for live Ca2+-imaging through a window in the head capsule. (b) Confocal image of GCaMP3-fluorescence expressed in OPN using GH146-Gal4, indicating the locations of the cell bodies, antennal lobe, and the calyx (scale bar: 50 µm). (c–d) Confocal image of the stimulus-induced change of fluorescence (peak amplitude) in the synaptic region of the calyx of a control and an hLRRK2G2019S knock-in animal (c) and quantification of fluorescence change (± SEM) over time (d). Arrowhead: time of stimulus application (10 mM nicotine, see Materials and methods). (e, f) Images of GFP fluorescence marking the synaptic area of OPN in the calyx of control (e) and hLRRK2G2019S knock-ins (f). Scale bar is 20 µm. (g) Quantification of GCaMP3 peak amplitude at OPN synapses in the calyx following stimulation (10 mM nicotine) in controls, wild-type coding DNA sequence (CDS) knock-ins, and in the PD knock-in mutants where the wild-type CDS is not (-) or is (+) expressed in OPN using GH146-Gal4 (OPN>wt CDS). Note that the hPINK1 control could not be determined as the combination of nSyb-Gal4>UAS-hPINK1 (expression in all neurons) in the Pink1 knock-out background interferes with the OPN-specific expression of Gal4 to drive UAS-GCaMP3 expression (g) (left bar; ‘nd’). In contrast, this issue is not present in hPINK1P399L mutant knock-in flies (nor the other flies used in the study) that could be rescued by OPN-specific expression of hPINK1 (g) (right bar); (‘OPN>wt CDS +’). (h) Quantification of the GFP fluorescence area of OPN synapses in the calyx (based on GCaMP3 signal) in controls, wild-type CDS knock-ins, and in the PD knock-in mutants where the wild-type CDS is not (-) or is (+) expressed in OPN (OPN>wt CDS). For (g, h) n≥5 animals per genotype. Bars are mean ± s.e.m. *, p<0.05 in ANOVA/Dunnett.

Olfactory projection neurons (OPN) dysfunction causes hyposmia in Parkinson’s disease (PD) mutants.

(a, b) Olfactory performance of PD knock-in flies and controls when given the choice between a blend of motor oil and banana odors (a) or a blend of beer and wine odors (b). (c, d) Olfactory performance of PD knock-in flies and LRRKKO when the relevant wild-type coding DNA sequence (CDS) is not (-) or is expressed selectively in OPN (‘OPN>wt CDS +’, green label) or T1 cholinergic interneurons (‘T1>hRab39B +’, red label), as a negative control. For (a–d) bars are mean ± s.e.m. n≥5 assays, ≥200 flies per genotype. *, p<0.05 in ANOVA/Dunnett.

Figure 6 with 1 supplement
Cholinergic neuron activity rescues dopaminergic defects that occur at later stages in the life of Parkinson’s disease (PD) mutants.

(a, left) Synaptic area of dopaminergic neuron (DAN) innervating the mushroom body in aged PD models (AuxR927G, synjR258Q, or LRRK2G2019S) and with or without GH146-Gal4-driven expression of the wild-type PD gene (aux or synj) or EndoAS75D, respectively. Bars: mean ± SEM. n≥5, *p<0.05 in ANOVA, Dunnett’s test. (a, right) Startle-induced negative geotaxis (SING) of the PD models with or without GH146-Gal4-driven expression of wild-type gene or endoAS75D. Points: mean ± SEM. n≥5, *p<0.05 in two-way ANOVA. Gray zone: variance of controls. (b) Odor choice performance, stimulus-induced changes in synaptic Ca2+ signal and olfactory projection neuron (OPN) synapse area of young controls and hLRRK2G2019S flies with or without chronic nicotine (Nic) feeding (up to 1 day before testing). Bars: mean ± SEM. n≥5 assays, *p<0.05 in ANOVA, Dunnett’s test. (c) SING, stimulus-induced changes in synaptic Ca2+ and DAN synapse area of aged controls and hLRRKG2019S flies with or without chronic application of nicotine. Bars: mean ± SEM. n≥5 assays, *p<0.05 in ANOVA, Dunnett’s test. (d) Confocal images of differentiated (60 days) wild-type and LRRK2G2019S ventral midbrain DAN labeled with the ventral midbrain marker FOXA2, dopaminergic marker TH, and neuronal marker MAP2. Scale bar: 20 µm. (e) Scheme of the treatment protocol and spontaneous Ca2+ activity (e’) and amplitude (e’’) of human induced DAN, 2 days after 20 days of no treatment (Ctrl), nicotine (Nic) treatment or nicotine+mecamylamine (Nic+Meca) treatment. Bars: mean ± SEM. n≥60 DAN from three independent differentiations, *p<0.05 in ANOVA, Dunnett’s test.

Figure 6—figure supplement 1
Generation and validation of human LRRK2G2019S ventral midbrain dopaminergic neuron (DAN).

(a) Schematic representation of gene editing strategy to knock-in the G2019S mutation in the LRRK2 locus. (b) Sanger sequencing of a single gene edited clone showing successful homozygous editing of the indicated nucleotide. (c) KOLF2-1J wild-type (LRRK2WT) control iPSCs and KOLF2-1J LRRK2G2019S/G2019S (LRRK2G2019S) iPSCs show normal expression of pluripotency markers OCT4 (yellow), SOX2 (purple), NANOG (purple), and TRA-1–81 (green). Nuclei are counterstained with DAPI (blue). Scale bar: 100 μm. (d) KOLF2-1J wild-type (LRRK2WT) control and KOLF2-1J LRRK2G2019S/G2019S (LRRK2G2019S) mutant ventral midbrain neural progenitor cells show normal expression of FOXA2 (magenta), LMX1A (green), and OTX2 (green), confirming that the neural progenitor cells are ventral midbrain-specific and capable of differentiating into DAN. Nuclei are counterstained with DAPI (blue). Scale bar: 100 μm. (e) Quantification of the frequency of TH+/MAP2+DAN. Statistical significance calculated with an ordinary t-test: ns, not significant. Error bars represent mean ± SEM. Data were collected in three independent vmDAN differentiations. (f) Differentiation protocol used.

Author response image 1
Upregulation of protein homeostasis rescues hyposmia across familial models of PD.

Results of a behavioral screen for cell-specific rescue of olfactory preference defects of young PD fly models using up and downregulation of deregulated genes in affected cell types. Genes implicated in the indicated pathways are over expressed or knocked down using GH146-Gal4 (OPN>) and UAS-constructs (over expression or RNAi) . UAS-only (-) and OPN>UAS (+) were scored in parallel and are compared to each other. n.d. not determined; Bars represent mean ± s.e.m.; gray zone indicates the variance of controls; n≥5 independent experiments per genotype, with ~50 flies each; red bars: p<0.05 in ANOVA and Bonferroni-corrected comparison to UAS-only control.

Tables

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (D. melanogaster)LrrkFlyBaseFBgn0038816
Gene (D. melanogaster)Pink1FlyBaseFBgn0029891
Gene (D. melanogaster)Rab39FlyBaseFBgn0029959
Gene (D. melanogaster)AuxFlyBaseFBgn0037218
Gene (D. melanogaster)SynjFlyBaseFBgn0034691
Strain, strain background (Escherichia coli)One Shot top10 cellsThermo Fisher ScientificC404010Competent cells
Genetic reagent (D. melanogaster)GH146-Gal4 (OPN-Gal4)Stocker et al., 1997; Bloomington Drosophila Stock CenterBDSC_30026; FlyBase: FBti0016783
Genetic reagent (D. melanogaster)GMR31F10-Gal4 (T1-Gal4)Jenett et al., 2012; Bloomington Drosophila Stock CenterBDSC_49685; FlyBase: FBti0134817
Genetic reagent (D. melanogaster)GMR58E02-Gal4 (PAM-Gal4)Jenett et al., 2012; Bloomington Drosophila Stock CenterBDSC_41347; FlyBase: FBti0137105
Genetic reagent (D. melanogaster)UAS-GCaMP3Tian et al., 2009; Bloomington Drosophila Stock CenterBDSC_32236; BDSC_32116
Genetic reagent (D. melanogaster)UAS-hLRRK2Venderova et al., 2009
Genetic reagent (D. melanogaster)UAS-hPINK1Yang et al., 2006
Genetic reagent (D. melanogaster)UAS-YFP-Rab39Zhang et al., 2007; Bloomington Drosophila Stock CenterBDSC_9825; FlyBase: FBti0100878
Genetic reagent (D. melanogaster)UAS-SynjVanhauwaert et al., 2017
Genetic reagent (D. melanogaster)UAS-AuxJacquemyn et al., 2023
Genetic reagent (D. melanogaster)LRRKk.o.This study (see Figure 1a and the Materials and methods section)Knock-out fly line carrying an attP-flanked mini-white gene replacing the first common exon in all transcripts of lrrk
Genetic reagent (D. melanogaster)hLRRK2This study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the wild-type human LRRK2 CDS at the lrrk locus
Genetic reagent (D. melanogaster)hLRRK2G2019SThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic human LRRK2G2019S CDS at the lrrk locus
Genetic reagent (D. melanogaster)hLRRK2Y1699CThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic human LRRK2Y1699C CDS at the lrrk locus
Genetic reagent (D. melanogaster)PINK1k.o.This study (see Figure 1a and the Materials and methods section)Knock-out fly line carrying an attP-flanked mini-white gene replacing the first common exon in all transcripts of pink1
Genetic reagent (D. melanogaster)hPINK1P399LThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic human PINK1P399L CDS at the pink1 locus
Genetic reagent (D. melanogaster)hPINK1L347PThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic human PINK1L347P CDS at the pink1 locus
Genetic reagent (D. melanogaster)Rab39k.o.This study (see Figure 1a and the Materials and methods section)Knock-out fly line carrying an attP-flanked mini-white gene replacing the first common exon in all transcripts of rab39
Genetic reagent (D. melanogaster)hRab39BThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the wild-type human Rab39B CDS at the rab39 locus
Genetic reagent (D. melanogaster)hRab39BG192RThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic human Rab39BG192R CDS at the rab39 locus
Genetic reagent (D. melanogaster)hRab39BT168KThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic human Rab39BT168K CDS at the rab39 locus
Genetic reagent (D. melanogaster)Auxk.o.This study (see Figure 1a and the Materials and methods section)Knock-out fly line carrying an attP-flanked mini-white gene replacing the first common exon in all transcripts of aux
Genetic reagent (D. melanogaster)AuxThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the wild-type Drosophila Aux CDS at the aux locus
Genetic reagent (D. melanogaster)AuxR927GThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic Drosophila AuxR927G CDS at the aux locus
Genetic reagent (D. melanogaster)Synjk.o.This study (see Figure 1a and the Materials and methods section)Knock-out fly line carrying an attP-flanked mini-white gene replacing the first common exon in all transcripts of synj
Genetic reagent (D. melanogaster)SynjThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the wild-type Drosophila Synj CDS at the synj locus
Genetic reagent (D. melanogaster)SynjR258QThis study (see Figure 1a and the Materials and methods section)Knock-in fly line carrying the pathogenic Drosophila SynjR258Q CDS at the synj locus
Cell line (H. sapiens)KOLF2-1J wild-type (LRRK2WT) control iPSCsPantazis et al., 2022; from the Jackson Laboratory under the iPSC Neurodegenerative Disease InitiativeProduct code: JIPSC001000
Cell line (H. sapiens)KOLF2-1J LRRK2G2019S/G2019S (LRRK2G2019S) iPSCsThis study (see Figure 6—figure supplement 1 and the Materials and methods section)Knock-in KOLF2-1J iPSC line carrying the pathogenic p.G2019S mutation in the LRRK2 locus
Transfected construct (H. sapiens)pLenti-hSyn1-mRuby2-T2A-GCaMP6fThis study (see the Materials and methods section)RRID:Addgene_197595Lentiviral transfer vector for the neuronal expression of mRuby2-T2A-GCaMP6f (fluorescent reporter for calcium imaging)
Biological sample (H. sapiens)Postmortem brain tissueParkinson's UK Brain Bank at the Imperial College LondonEC reference NH019 2019-02-01
AntibodyMouse monoclonal IgG1 anti-SOX2Santa Cruzsc-365823IF (1:200)
AntibodyRabbit polyclonal anti-OCT4Abcamab19857IF (1:50)
AntibodyMouse monoclonal IgG1 anti-NANOGSanta Cruzsc-293121IF (1:50)
AntibodyMouse monoclonal IgM anti-TRA-1–81SigmaMAB4381IF (1:100)
AntibodyRabbit polyclonal anti-LMX1A/BMilliporeAB10533IF (1:1000)
AntibodyMouse monoclonal IgG2a anti-FOXA2/HNF-3βSanta Cruzsc-101060IF (1:250)
AntibodyMouse monoclonal IgG1k anti-OTX2Santa Cruzsc-514195IF (1:100)
AntibodyRabbit polyclonal anti-THSigmaAB152IF (1:500)
AntibodyMouse monoclonal IgG1 anti-MAP2SigmaM1406IF (1:1000)
AntibodyGoat polyclonal Alexa Fluor-488 conjugated secondary antibodyInvitrogenAnti-rabbit: # A-11034 for OCT4/TH Anti-mouse IgG1: # A-21121 for Nanog/OTX2IF (1:500)
AntibodyGoat polyclonal Alexa Fluor-555 conjugated secondary antibodyInvitrogenAnti-mouse IgG1: # A-21127 for SOX2 Anti-mouse IgG2a: # A-21137 for FOXA2IF (1:500)
AntibodyDonkey/Goat polyclonal Alexa Fluor-647 conjugated secondary antibodyInvitrogenAnti-rabbit: # A-31573 for LMX1A Anti-mouse IgM: # A-21238 for TRA-1–81 Anti-mouse IgG1: # A-21240 for MAP2IF (1:500)
Recombinant DNA reagentpCFD4: U6:1-gRNA U6:3-gRNAPort et al., 2014RRID:Addgene_49411
Recombinant DNA reagentpUC19Norrander et al., 1983RRID:Addgene_50005
Recombinant DNA reagentpUAST.attBBischof et al., 2007GenBank:EF362409.1
Recombinant DNA reagentpFL44S{w+}-attBKhuong et al., 2013
Recombinant DNA reagentpWhite-STARChoi et al., 2009
Recombinant DNA reagent2XMyc-LRRK2-WTGreggio et al., 2008RRID:Addgene_25361
Recombinant DNA reagentpLenti6-DEST PINK1-V5 WTBeilina et al., 2005RRID:Addgene_13320
Sequence-based reagentPrimers, gRNAs, oligos, gBlocks, ssODNIntegrated DNA Technologies (IDT)See Supplementary file 1
Sequence-based reagentiPSC gRNASynthegoSee Supplementary file 1
Peptide, recombinant proteinCas9 HiFi proteinIntegrated DNA Technologies (IDT)Alt-R S.p. HiFi Cas9 Nuclease V3, 500 µg; #1081061
Peptide, recombinant proteinDispaseSigma-AldrichD4818
Peptide, recombinant proteinCollagenaseInvitrogen, CA, USA17100-017
Chemical compound, drugRinger's solutionEstes et al., 1996
Chemical compound, drugTMREAbcam, Cambridge, UKab113852
Chemical compound, drugTrypsin-EDTAInvitrogen, CA, USA25300054
Chemical compound, drugNicotineSigma-AldrichN3876
Commercial assay, kitNEBuilder HiFi DNA AssemblyNew England Biolabs, Inc (NEB)NEB #E2621
Commercial assay, kitMaxwell RSC instrument and kitPromegaAS8500, AS1340
Commercial assay, kitChromium Single Cell 3′ Library & Gel Bead Kit v210x GenomicsPN-120237
Commercial assay, kitChromium Single Cell A Chip Kit10x GenomicsPN-1000009
Commercial assay, kitChromium i7 Multiplex Kit10x GenomicsPN-120262
Commercial assay, kitChromium Next GEM Single Cell 3' Kit v310x GenomicsPN-1000075
Commercial assay, kitChromium Next GEM Chip B Single Cell Kit10x GenomicsPN-1000073
Commercial assay, kitSingle Index Kit T Set A10x GenomicsPN-1000213
Software, algorithmFijiSchindelin et al., 2012RRID:SCR_002285https://imagej.net/Fiji
Software, algorithmCell Ranger (v 3.0.2)10x GenomicsRRID:SCR_017344https://www.10xgenomics.com/support/software/cell-ranger/latest
Software, algorithmGOrillaEden et al., 2009RRID:SCR_006848http://cbl-gorilla.cs.technion.ac.il/
Software, algorithmGraphPad PrismGraphPad software, CA, USARRID:SCR_002798https://www.graphpad.com/scientific-software/prism/
Software, algorithmVib-singlecell-nf/vsn-pipelines (v 0.27.0)Flerin et al., 2021https://github.com/vib-singlecell-nf/vsn-pipelines
Software, algorithmPython (v 3.7.3)PythonRRID:SCR_008394http://www.python.org/
Software, algorithmpandas (v 0.25.1)GitHubRRID:SCR_018214https://github.com/pandas-dev/pandas
Software, algorithmnumpy (v 1.19.1)Harris et al., 2020RRID:SCR_008633http://www.numpy.org/
Software, algorithmmatplotlib (v 3.3.0)GitHubRRID:SCR_008624https://github.com/matplotlib/matplotlib
Software, algorithmseaborn (v 0.9.0)GitHubRRID:SCR_018132https://github.com/mwaskom/seaborn
Software, algorithmstatsmodels (v 0.10.1)Seabold and Perktold, 2010RRID:SCR_016074http://www.statsmodels.org/
Software, algorithmscipy (v 1.4.1)Virtanen et al., 2020RRID:SCR_008058http://www.scipy.org/
Software, algorithmadjustText (v 0.7.3)GitHubRRID:SCR_022260https://github.com/Phlya/adjustText
Software, algorithmloompy (v 3.0.6)GitHubRRID:SCR_016666https://github.com/linnarsson-lab/loompy
Software, algorithmsklearn (v 0.20.1)Pedregosa et al., 2011RRID:SCR_019053https://scikit-learn.org/stable/modules/generated/sklearn.decomposition.NMF.html
Software, algorithmpatsy (v 0.5.1)Smith, 2018https://github.com/pydata/patsy
Software, algorithmvenn (v 0.1.3)Pythonhttps://pypi.org/project/venn/
Software, algorithmR (v 3.6.1)R FoundationRRID:SCR_001905http://www.r-project.org/
Software, algorithmMAST (v 1.12.0)GitHubRRID:SCR_016340https://github.com/RGLab/MAST/
Software, algorithmSCopeLoomR (v 0.5.1)De Waegeneer et al., 2022https://github.com/aertslab/SCopeLoomR
OtherSCope Resource WebsiteDavie et al., 2018Visualization tool for large-scale sc-RNA-seq datasets; http://scope.aertslab.org

Additional files

Supplementary file 1

List of gRNAs, gBlocks, and oligos used to generate and validate the fly lines and hiPSC line generated in this study.

List of full fly genotypes and hiPSC lines used in this study.

https://cdn.elifesciences.org/articles/98348/elife-98348-supp1-v1.xlsx
Supplementary file 2

Full results of the Cell Ranger summaries (10x Genomics).

https://cdn.elifesciences.org/articles/98348/elife-98348-supp2-v1.xlsx
Supplementary file 3

Summary of the number of cells per cell cluster for each mutant and controls, including the raw number of cells, the percentage of cells for each cluster per mutant, and the percentage of cells for each mutant and controls per cluster.

https://cdn.elifesciences.org/articles/98348/elife-98348-supp3-v1.xlsx
Supplementary file 4

Summary of the parameters used to model differentially expressed gene (DEG)-cell number correlation.

https://cdn.elifesciences.org/articles/98348/elife-98348-supp4-v1.xlsx
Supplementary file 5

List of cell types and respective z-residuals (deviation from model) to determine affected cell types and those that are outside of the 95% confidence interval.

https://cdn.elifesciences.org/articles/98348/elife-98348-supp5-v1.xlsx
Supplementary file 6

List of the clusters with significant changes (based on cell number-adjusted differentially expressed genes [DEGs]) for each model, black cell clusters (cells with significant number of DEG) are listed as TRUE.

https://cdn.elifesciences.org/articles/98348/elife-98348-supp6-v1.xlsx
Supplementary file 7

Summarized results of the differentially expressed gene (DEG) analysis in fly model brains and postmortem human samples and the summary of the SynGO analysis.

https://cdn.elifesciences.org/articles/98348/elife-98348-supp7-v1.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/98348/elife-98348-mdarchecklist1-v1.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ulrike Pech
  2. Jasper Janssens
  3. Nils Schoovaerts
  4. Sabine Kuenen
  5. Carles Calatayud Aristoy
  6. Sandra F Gallego
  7. Samira Makhzami
  8. Gert J Hulselmans
  9. Suresh Poovathingal
  10. Kristofer Davie
  11. Adekunle T Bademosi
  12. Jef Swerts
  13. Sven Vilain
  14. Stein Aerts
  15. Patrik Verstreken
(2025)
Synaptic deregulation of cholinergic projection neurons causes olfactory dysfunction across five fly Parkinsonism models
eLife 13:RP98348.
https://doi.org/10.7554/eLife.98348.3