Analysis of Essential Genes in Clostridioides difficile by CRISPRi and Tn-seq

Abstract

Essential genes are interesting in their own right and as potential antibiotic targets. To date, only one report has identified essential genes on a genome-wide scale in Clostridioides difficile, a problematic pathogen for which treatment options are limited. That foundational study used large-scale transposon mutagenesis to identify 404 protein-encoding genes as likely to be essential for vegetative growth of the epidemic strain R20291. Here, we revisit the essential genes of strain R20291 using a combination of CRISPR interference (CRISPRi) and transposon-sequencing (Tn-seq). First, we targeted 181 of the 404 putatively essential genes with CRISPRi. We confirmed essentiality for >90% of the targeted genes and observed morphological defects for >80% of them. Second, we conducted a new Tn-seq analysis, which identified 346 genes as essential, of which 283 are in common with the previous report and might be considered a provisional essential gene set that minimizes false positives. We compare the list of essential genes to those of other bacteria, especially Bacillus subtilis, highlighting some noteworthy differences. Finally, we used fusions to red fluorescent protein (RFP) to identify 18 putative new cell division proteins, three of which are conserved in Bacillota but of largely unknown function. Collectively, our findings provide new tools and insights that advance our understanding of C. difficile.

Keywords: CRISPRi; Tn-seq.

Fig. 1.. Summary of gene essentiality determined by CRISPRi and Tn-seq.

(A) CRISPRi-induced viability defects were determined from spot titer assays on TY-thiamphenicol with 1% xylose. Viability defects were scored as strong (≥ 1000-fold), moderate (≥100-fold), weak (≥10-fold or full viability but colonies were small), or none (full viability, normal colony size). In cases where the two sgRNAs produced different results, the stronger viability defect was used. (B) Comparison of Tn-seq datasets. Of the 346 genes determined to be essential in our study, 283 were essential in the Dembek et al. set, and 10 were called ambiguous. Conversely, of the 404 Dembek et al. essential genes, 283 were essential in our dataset, 74 were non-essential, 12 were not found owing to use of different genome annotations, and 35 were ambiguous (Here, ambiguous combines three categories from Supplemental Table 3: unclear (11 genes), short (6 genes), unclear/non-essential (18 genes)). (C) Viability defects in CRISPRi correlate with likelihood a gene will be scored as essential by Tn-seq. Viability defects are from Table S1. Tn-seq calls come from Table S3.

Fig. 2.. Morphology of CRISPRi strains with sgRNAs targeting genes in select functional pathways.

Left: pathway. Middle: Predicted transcription unit. Targeted genes are boxed and indicated above the operon diagrams. Numbers are R20291 locus tags. Genes are color coded to indicate essentiality based on Tn-seq calls in Table S3. Blue: essential. Light blue: ambiguous. White: non-essential. Operon structure is not to scale. Right: Morphological changes based on phase contrast and fluorescence micrographs of cells scraped from viability plates. Membranes were stained with FM4–64 and DNA was stained with Hoechst 33342. Size bars are 5 μm. The control strain expressed an sgRNA that does not target anywhere in the genome. Micrographs are representative of at least two experiments. Figure S1 shows microscopy of more genes.

Fig. 3.. Essentiality follow-up.

(A) Transposon insertion profile for murJ2. Vertical lines represent mapped insertion sites and are scaled to indicate the number of sequence reads mapping to that site. Although murJ2 sustained numerous insertions, ~80% were in the last 10% of the gene, suggesting the non-essentiality call by TRANSIT2 is incorrect. (B) Transposon insertion profile for polA indicating that only the N-terminal domain is essential. Read frequency was scaled to 5 to highlight the absence of reads in the N-terminal domain. The average number of reads per polA site with at least one read was 173. (C) Spot titer assays of CRISPRi strains targeting genes in the atp operon. Serial dilutions of overnight cultures were spotted on TY-Thi10 plates with 1% xylose. Plates were imaged after incubation at 37˚C for ~18h. Silencing atpB and atpD resulted in small colonies, while growth after silencing atpI and atpF was comparable to the negative control. (D) Pre-depletion of ATP synthase proteins impairs growth. Starter cultures were grown overnight in TY-Thi10 without (left) or with (right) 1% xylose, then subcultured into TY-Thi10 with 1% xylose and growth was followed by measuring optical density at 600 nm. To prolong growth, cultures in the left panel were back-diluted at 7h. (E) Zone of inhibition assays reveal CRISPRi knockdown of the atp operon increases sensitivity to CCCP. Plates were imaged after incubation at 37˚C for ~18h. Novobiocin (Novo) served as a control. Guides in panels C-E were: atpI (5531), atpB (5583), atpF (5581), atpD (5579) or a negative control that does not target anywhere in the gemone.

Fig. 4.

Representative fluorescence micrographs of fixed cells that produced the indicated proteins fused to RFP. Percentages indicate the fraction of cells scored positive for septal localization (n ≥ 202 cells). Space bar = 10 μm.