A two-step enrichment process facilitates robust transcriptional profiling of Bb during the tick bloodmeal.

(A) Schematic of Bb during nymphal I. scapularis feeding. Bb in the nymphal tick midgut respond to the nutrient-rich bloodmeal by multiplying and changing their transcriptional state (Ouyang et al., 2012; de Silva and Fikrig, 1995). At the same time, the tick gut undergoes numerous changes to digest the bloodmeal (Caimano et al., 2015; Sonenshine and Anderson, 2014). After two to three days of feeding, a small number of Bb leave the midgut and enter the salivary glands (blue), while the majority are left behind in the gut after engorgement (Dunham-Ems et al., 2009). (B) Schematic of Bb enrichment process from feeding ticks. Whole ticks are dissociated, αBb antibodies are added to lysates, and antibodies and Bb are captured magnetically. RNA is extracted and RNA-seq libraries are prepared. DASH is then used to remove rRNA before sequencing. This process increases Bb reads in the resulting sequencing data. (C) RT-qPCR results showing the percentage of Bb flaB and I. scapularis gapdh RNA in the enriched versus depleted fractions after the enrichment process. Data come from 4 replicates each from day 2, day 3, and day 4, mean +/− SE. ****p-value<0.0001, paired t test. Nearly all Bb flaB RNA was found in the enriched fraction. (D) The percentage of reads mapping to rRNA before and after DASH. n=4. Data are shown as mean +/− SD. ****p-value<0.0001, paired t test. rRNA reads are drastically reduced after DASH. (E) The percentage of reads in RNA-seq libraries mapping to Bb. Bb mRNA reads make up a larger proportion of libraries than without enrichment. Data are shown as mean +/− SD. (F) The number of reads in millions (M) mapped to Bb for each day. n=4. Data are shown as mean +/− SD. An average of 4.3 million reads per sample mapped to Bb genes, covering 92% of genes with at least 10 reads.

Global ex vivo profiling of Bb reveals extent and kinetics of transcriptional changes.

(A) Principal component analysis of samples from across feeding. PC1 correlates strongly with day of feeding. (B) Schematic depicting how data was analyzed, as pairwise comparisons between the first day after attachment and all other days. (C-E) Volcano plots of differentially expressed genes comparing day 2 versus day 1(C), day 3 versus day 1 (D), and day 4 versus day 1 (E). The total number of upregulated genes is shown in the top right and the number of downregulated genes is shown in the top left. Yellow dots are genes that first change expression between day 1 and day 2, red dots are genes that first change expression between day 1 and day 3, and purple dots are genes that first change expression between day 1 and day 4. Two genes with log2 fold changes >4 are shown at x=4, and five genes with -log10(padj) >60 are shown at y=60. Only genes with p-value < 0.05 from Wald tests and at least a two-fold change are highlighted. n=4. By day 4 of feeding, 153 genes are upregulated and 33 genes are downregulated from day 1 baseline levels.

Bb genes upregulated during feeding are found predominantly on plasmids.

(A) Schematic of the chromosome and plasmids in the Bb B31-S9 genome. Plasmid names denote whether the plasmid is linear (lp) or circular (cp) and the length of plasmids in kilobases (kb). For example, lp17 is a 17 kb linear plasmid. Genome is shown approximately to scale. (B-C) The number of genes from each chromosome or plasmid that increased (B) or decreased expression (C) twofold during feeding. Upregulated genes are distributed across plasmids, while downregulated genes are found on the chromosome and lp54.

Bb genes encoding outer surface proteins are highly prevalent among upregulated genes.

(A) The number of Bb genes that change over the course of tick feeding sorted into functional categories. Genes that first change 2 days after attachment are shown in yellow, 3 days after attachment in red, and 4 days after attachment in purple. A majority of upregulated genes fall into cell envelope and unknown categories. (B) Schematic of the outer membrane of Bb showing outer surface lipoproteins. Lipoproteins can also reside in the periplasmic space. (C) Heat map of expression levels of all genes encoding outer surface lipoproteins as average Transcripts Per Million (TPM) across the 4 days of tick feeding. Gene names highlighted in blue were two-fold upregulated and genes in pink two-fold downregulated over feeding (see Figure 2). A majority of genes encoding outer surface proteins increased in expression throughout feeding, while having different magnitudes of expression.

Identification of candidate tick interaction partners of Bb cells ex vivo

(A) Schematic of experiment to determine candidate tick proteins interacting with Bb over the course of feeding. Ticks were collected 1 day after attachment and 4 days after attachment. Uninfected ticks at the same time points were mixed with cultured Bb as controls. Bb was enriched with αBb antibody as in RNA-seq experiments and then subjected to mass spectrometry to identify tick proteins present in the samples. Venn diagram depicts the proteins enriched in day 1 and day 4 samples over controls in at least two of three replicates. Tick proteins that are enriched with Bb vary greatly over the course of feeding. (B) Tick proteins uniquely identified one day after attachment that are annotated as extracellular matrix (ECM) proteins. (C) Tick proteins uniquely identified four days after attachment that are annotated as low-density lipoprotein receptors. ECM and membrane proteins may be good candidates for Bb-interacting proteins.

αBb antibody recognizes OspA and binds Bb in the tick throughout the bloodmeal

(A) Western blot with αBb on lysate from cultured Bb: wildtype (A3, left), a mutant lacking ospA (ospA1), and the mutant with ospA restored (ospA+B1)(Battisti et al., 2008). αBb recognizes OspA among other proteins. (B) Immunofluorescence microscopy with αBb (green, left) and propidium iodide (PI) (DNA, red, center) on each day of feeding (merge is yellow, right). αBb antibody recognizes Bb in the tick across the bloodmeal.

Enrichment process does not induce large scale gene expression changes in in vitro cultured Bb

Log2 fold changes vs. mean normalized number of counts comparing cultured Bb input and samples after enrichment with αBb. n=3. Red dots, p-value < 0.05, Wald tests. The gene expression changes induced during processing are much smaller than those observed between days of feeding.

Ex vivo RNA-seq corroborates transcriptional programs in the tick

(A) Tukey style boxplot of Transcripts Per Million (TPM) on each day for rpoS. Black dots represent replicates. ****p-value < 0.00001, Wald test. rpoS expression increases over the course of feeding. (B) Volcano plot of DE genes comparing day 4 to day 1, with RpoS-upregulated genes (blue). Genes upregulated by RpoS in ticks increase during feeding. (C) Volcano plot of DE genes comparing day 4 to day 1, with Rrp1-upregulated (blue) and downregulated (pink) genes. Rrp1-regulated genes correlate well with genes up and downregulated during feeding. (D) Volcano plot of DE genes comparing day 4 to day 1, with RelBbu-upregulated (blue) and downregulated (pink) genes. About half of RelBbu genes change in the expected direction over feeding.

Genes changing over tick feeding overlap with genes that change expression in previously probed tick feeding contexts

(A) The overlap of two-fold changed genes with genes that changed expression in Bb cultures grown in conditions mimicking feeding ticks. Genes are grouped based on the first day that they changed two-fold from day 1. “Fed tick” culture conditions were 37°C, pH 6.8 in Revel et al. 2002 and 35°C, pH 7.4 in Ojaimi et al. 2003, and genes elevated in these conditions in one or both studies are highlighted in red. “Unfed tick” culture conditions were 23°C, pH 7.5 in Revel et al. 2002 and 23°C, pH 7.4 in Ojaimi et al. 2003, and genes elevated in these conditions in one or both studies are highlighted in teal. Genes that were not elevated in either condition in those studies are in gray. Particularly for genes that increase on day 2, there is a large overlap with genes elevated in “fed tick” culture conditions in previous studies. (B) The overlap of two-fold changed genes with genes that changed expression between fed nymphs and dialysis membrane chambers (DMCs) mimicking mammalian conditions. Genes are grouped based on the first day that they changed two-fold from day 1. Bb expression in fed nymphs versus in DMCs was compared by Iyer et al. 2015 using bacterial RNA amplification and microarray, while Grassmann et al. 2023 used TBDCapSeq. Genes elevated in fed nymphs in one or both studies are highlighted in red, while genes elevated in DMCs in one or both studies are highlighted in purple. Genes that were elevated in conflicting conditions between the two studies are in dark gray, while genes not elevated in either condition are in light gray. For genes that increase on day 2, there is a large overlap with genes elevated in fed nymphs, while genes that increase first on days 3 and 4 have a larger overlap with genes elevated in DMCs.