Single-cell transcriptomics reveals a multiphasic Wolbachia host infection trajectory
Data files
Jul 17, 2026 version files 14.12 GB
-
JW18DOX_EGFP.tiff
1.65 MB
-
JW18DOX_ER-tracker.tiff
2.28 MB
-
JW18DOX_HOECHST_NucBlue.tiff
1.57 MB
-
JW18DOX_Lysotracker.png
2.58 MB
-
JW18DOX_MitoTracker_Alexa647.tiff
1.47 MB
-
JW18DOX_mitotracker.png
2.75 MB
-
JW18wMel_EGFP.tiff
1.50 MB
-
JW18wMel_ER-tracker.tiff
2.17 MB
-
JW18wMel_HOECHST_NucBlue.tiff
1.61 MB
-
JW18wMel_Lysotracker.png
2.86 MB
-
JW18wMel_MitoTracker_Alexa647.tiff
1.46 MB
-
JW18wMel_mitotracker.png
2.55 MB
-
Lysotracker.zip
5.76 GB
-
Mitotracker.zip
8.33 GB
-
README.md
4.52 KB
Abstract
To establish persistent infections, intracellular bacterial symbionts must navigate complex host cellular environments, co-opt host biology, and evade immune clearance, yet the molecular mechanisms underlying infection remain poorly characterized. The alphaproteobacterial endosymbiont Wolbachia pipientis, prevalent in diverse arthropods and nematodes, exemplifies this challenge. The wMel strain, natively associated with Drosophila melanogaster, has proven particularly effective in biocontrol applications for blocking viral transmission and suppressing reproduction when transinfected into mosquito populations. Despite these successful applications, the molecular mechanisms underlying initial infection establishment remain poorly understood in wMel’s native host. Importantly, Wolbachia demonstrates host cell-type preferences and induces tissue-specific phenotypes, yet bulk RNA-sequencing studies have yielded inconsistent signals of differential expression, likely due to signal averaging across cells with heterogeneous infection titers, states, and developmental stages. Here, we employed single-cell RNA sequencing to investigate how wMel colonization influences the host transcriptome during the establishment of stable infection of D. melanogaster JW18 somatic cell lines. First, we use 10X Genomics Chromium 3’ single cell RNA-seq to validate the lower-cost Illumina 3’ single cell RNA-seq (PIPseq) platform and show that mis-priming of symbiont and host ribosomal RNAs can be leveraged as a proxy to measure intracellular bacterial titer.Using the PIPseq platform to characterize six timepoints across the three months required to establish stable infection, we demonstrate that nascent wMel infections induce distinct transcriptional changes that create novel cellular states diverged from uninfected controls. Through unsupervised clustering, pseudotime analysis, and marker gene identification, we reveal how Wolbachia reprograms host somatic cells to evade immune responses, establish infection, and facilitate nutrient acquisition. These findings uncover the hidden role Wolbachia plays in eukaryotic cellular biology and provide mechanistic insights from the native host cell system that will enable future efforts to develop advanced control measures leveraging cell-type specific phenotypes of Wolbachia strains.
Dataset DOI: 10.5061/dryad.612jm64md
Description of the data and file structure
High-content imaging data for Wolbachia-infected and uninfected JW18 cells
Overview
This dataset contains high-content fluorescence microscopy images and derived per-cell quantification data from Wolbachia wMel-infected and uninfected Drosophila melanogaster JW18 somatic cell lines. Cells were stained for lysosomal and mitochondrial markers and imaged on a PerkinElmer/Revvity Opera Phenix High-Content Screening System to assess endosomal/organellar signal in relation to Wolbachia infection status.
Associated publication
Jacobs J, Lum A, Nykamp J, Lagousis CRM, Russell SL. Single-cell transcriptomics reveals a multiphasic Wolbachia host infection trajectory.
This imaging dataset corresponds to Figure 5 and the associated Methods subsection describing lysosomal/mitochondrial imaging.
Experimental design
wMel-infected and uninfected JW18 cells were seeded into 96-well imaging plates (100 µL per well) in Shields and Sang M3 (or Schneider's) medium supplemented with 10% FBS. Lysosomes and mitochondria were labeled in separate parallel wells (i.e., each well received only one organelle-specific tracker dye, alongside the Hoechst nuclear/Wolbachia-chromosome counterstain).
Staining protocol
- Lysosomes: LysoTracker Red DND-99 (Invitrogen, L7528; ex/em 577/590 nm), diluted from a 1 mM stock to a 50 nM working concentration in room-temperature medium, applied for 1.5 h.
- Mitochondria: MitoTracker Deep Red FM (Invitrogen, M22426; ex/em 644/665 nm), applied at 250 nM for 20 min.
- Nuclear/Wolbachia chromosome counterstain (all wells): NucBlue Live ReadyProbes Reagent (Hoechst 33342; Invitrogen, R37605), 2 drops per mL, added in the same medium as the tracker dye.
- The lowest effective dye concentration was used for each tracker to minimize labeling artifacts.
- Culture medium was replaced with dye-containing medium for the incubation period, then removed and replaced with fresh (dye-free) medium prior to imaging.
Imaging parameters
- Instrument: Opera Phenix High-Content Screening System (PerkinElmer/Revvity)
- Objectives: 20X or 63X (see file naming/metadata for objective used per image set)
- Exposure: 120 ms
- Excitation power: 50 %
- Channels acquired: brightfield, Hoechst (nuclear/DNA), and the relevant tracker channel (LysoTracker Red or MitoTracker Deep Red) per well
Image analysis
Images were analyzed in Harmony software (PerkinElmer/Revvity):
- Cells were segmented from the brightfield channel outline.
- Nuclei were segmented from the Hoechst channel.
- Per-cell lysosomal or mitochondrial signal intensity was quantified from the respective tracker channel.
Statistical comparisons between cell types (infected vs. uninfected) and infection states were performed using Welch's two-tailed t-test (scipy.stats.ttest_ind).
Contact
Jodie Jacobs, jacobs.jodiem@gmail.com, Russell Lab, Department of Biomolecular Engineering, University of California, Santa Cruz
Shelbi Russell, shelbilrussell@gmail.com, Russell Lab, Department of Biomolecular Engineering, University of California, Santa Cruz
Files and variables
File: JW18DOX_EGFP.tiff
Description: Representative Image
File: JW18DOX_HOECHST_NucBlue.tiff
Description: Representative Image
File: JW18DOX_MitoTracker_Alexa647.tiff
Description: Representative Image
File: JW18DOX_Lysotracker.png
Description: Representative Image
File: JW18DOX_ER-tracker.tiff
Description: Representative Image
File: JW18DOX_mitotracker.png
Description: Representative Image
File: JW18wMel_EGFP.tiff
Description: Representative Image
File: JW18wMel_ER-tracker.tiff
Description: Representative Image
File: JW18wMel_Lysotracker.png
Description: Representative Image
File: JW18wMel_HOECHST_NucBlue.tiff
Description: Representative Image
File: JW18wMel_MitoTracker_Alexa647.tiff
Description: Representative Image
File: JW18wMel_mitotracker.png
Description: Representative Image
File: Lysotracker.zip
Description: Analysis and associated data
File: Mitotracker.zip
Description: Analysis and associated data
We seeded wMel-infected and uninfected JW18 cells into 96-well imaging plates (100 µL per well) in Shields and Sang M3 (or Schneider's) medium supplemented with 10% FBS. We labeled lysosomes and mitochondria in separate parallel wells. We diluted LysoTracker Red DND-99 (Invitrogen, L7528; ex/em 577/590 nm) from a 1 mM stock to a 50 nM working concentration in room-temperature medium and applied it for 1.5 h; we applied MitoTracker Deep Red FM (Invitrogen, M22426; ex/em 644/665 nm) at 250 nM for 20 min. We used the lowest effective dye concentration to minimize labeling artifacts. We counter-stained Wolbachia chromosomes and Drosophila nuclei in the same medium with NucBlue Live ReadyProbes Reagent (Hoechst 33342; Invitrogen, R37605) at 2 drops per mL. We replaced culture medium with dye-containing medium for the incubation, then removed it and replaced it with fresh medium prior to imaging. We imaged plates on an Opera Phenix High-Content Screening System (PerkinElmer/Revvity) using a 20X or 63X objective (120 ms exposure, 50 % excitation power), acquiring brightfield, Hoechst, and the relevant LysoTracker or MitoTracker channel. We analyzed images in Harmony software: we segmented cells from the brightfield outline and nuclei from the Hoechst channel, and quantified per-cell lysosomal or mitochondrial signal from the respective tracker channel. We tested for statistical differences in these measurements between cell types and infection states with Welch's two-tailed t-test from stats.ttest_ind.
