A randomized multiplex CRISPRi-Seq approach for the identification of critical combinations of genes
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, United States
- Great Lakes Bioenergy Research Center, University of Wisconsin-Madison, United States
- Wisconsin Energy Institute, University of Wisconsin-Madison, United States
- Department of Bacteriology, University of Wisconsin-Madison, United States
Figures
Figure 1
Intracellular growth of Lp02(dcas9) upon CRISPRi gene silencing with predetermined 10-plex arrays.
(A) Gene silencing efficiency assay. Each of the ten genes targeted by one of the six indicated MC arrays was assayed for gene repression by quantitative polymerase chain reaction (qPCR). RNA levels in the Lp02(dcas9) strain bearing the CRISPR array were compared to that bearing an empty vector. Data are a summary of two replicates shown side-by-side. Horizontal bars indicate mean fold repression of all targets of the array. (B) L. pneumophila intracellular growth assay. Host cells were challenged at a multiplicity of infection (MOI) of 0.05 (U937 macrophages) or 0.03 (Acanthamoeba castellanii) with Lp02(dcas9) bearing the indicated MC arrays. Colony forming units (CFUs) of samples taken 2 hr post infection (hpi) and either 72 (U937 macrophages) or 48 (A. castellanii) hpi were determined, and fold growth of Lp02(dcas9) bearing the indicated MC arrays relative to that of Lp02(dcas9) bearing the empty vector was plotted. Numbers below 1 indicate a growth defect upon gene silencing. BD, below detection (arbitrarily set to –100). Horizontal bars indicate mean fold growth vs. vector using data from three replicates.
Figure 2 with 1 supplement
Cloning strategy for de novo self-assembly of CRISPR array libraries.
(A) Nucleotide sequence of R-S-R building blocks. Each R-S-R element is composed of a top and bottom oligonucleotide. The 24 bp spacer sequence (black font; shown for crRNAmavN as an example) is flanked by sequences of the upstream (red) and downstream (blue) repeat element, with each end containing sticky overhangs (green). (B) Nucleotide sequence of attB4r-R and R-attB3r ‘dead ends’ with only one sticky overhang each (green). (C) Array self-assembly, size selection, and cloning. R-S-Rs and ‘dead end’ elements were allowed to self-assemble in a single tube and then ligated together and incorporated into an interim cloning vector. Arrays excised from an interim cloning vector were subjected to gel electrophoresis to separate them according to size. RE, restriction enzyme. After gel extraction, arrays of 550–650 bps and 650–800 bps in size were cloned into a donor plasmid such that three additional elements, namely the tracrRNA, a tet promoter, and a rrnB T1 terminator, could assemble with the arrays into the destination vector by Invitrogen Multisite Gateway Pro cloning. Final constructs were introduced into Lp02(dcas9) by electroporation.
Figure 2—figure supplement 1
Model of the Legionella-containing vacuole bearing transmembrane effectors.
L. pneumophila translocates >300 effector proteins via a T4SS into a host cell to re-wire cellular pathways to evade the lysosome and establish a replication vacuole known as the Legionella-containing vacuole (LCV). It is hypothesized that transmembrane effectors are incorporated into the LCV for nutrient acquisition, detoxification, and membrane fusion. Diagram made using BioRender (https://biorender.com).
Figure 3 with 1 supplement
Self-assembled CRISPR arrays are diverse in length and spacer composition.
(A) Two CRISPR array libraries were built from 44 R-S-R building blocks. Vectors from each library were harvested four times and sequenced. Venn diagrams show the overlap of unique spacer combinations, with the requirement of five or more raw (pre-sequence-depth normalized) read counts per array in each round of sequencing. (B) The distribution of plex-ness (spacer count) for arrays found in all four rounds of sequencing. (C) Chord diagrams link transmembrane effector (TME) genes, listed around the outside of the circle, each time spacers targeting each is present in an array. Link line width is weighted according to the number of times the combination of spacers was observed. Link line color is unique for each spacer and is constant between the two diagrams.
Figure 3—figure supplement 1
Spacer abundance post CRISPRi induction.
The occurrence of each spacer within the plasmid libraries was quantified 24 hr post induction of expression. The fact that each spacer was well represented within the library indicated that the crRNA products of these spacers were not toxic to axenic L. pneumophila growth.
Figure 4
Schematic overview of the experimental and bioinformatics pipeline of MuRCiS.
(A) During MuRCiS, a pooled population of L. pneumophila (dcas9) bearing the multiplex random CRISPR arrays were subjected to a selective pressure (intracellularity). Vectors were purified from both input and output bacteria populations, linearized, and submitted for long-read PacBio Sequel sequencing. For simplification, an array of mixed spacer population is shown as a single stretch of color, each color representing a different combination of spacers. In this example, the yellow and red arrays are lost in the output suggesting they silence critical combinations of genes. (B) Overview of the custom bioinformatics pipeline used to identify unique spacer combinations causing growth attenuation as defined by a fold reduction of five or more (https://github.com/GLBRC/MuRCiS_pipeline,copy archived at Myers and Donohue, 2023).
Figure 5
Discovery of L. pneumophila gene combinations critical for growth in U937 macrophages.
(A) Correlation grid plotting all pairwise spacer combinations. The number of times two unique spacers were present in the sequenced array library is indicated by color-coding, ranging from light (rare) to dark (abundant). Black boxes, same spacer. (B) Total number of reads bearing the spacer which encodes crRNAmavN. Counts in technical output replicates were summed. (C) Flowchart going from the number of total spacer combinations to the number of critical combinations identified. (D) Histogram indicating the plex-ness of all critical spacer combinations identified. (E) A list of virulence-critical combinations of genes and the corresponding fold reduction observed for each in each experiment. Library 1 results are shown above Library 2 results. (F) Intracellular growth assays of strains with deletions in the indicated genes. Results are given as fold growth (colony forming units [CFUs] harvested 72 hr post infection (hpi) vs. 2 hpi) compared to that of Lp02 (WT). Horizontal bars indicate mean fold growth of the deletion strain vs. WT strain for two or more experimental replicates. NG, no growth (arbitrarily set to –1000). (G) Intracellular growth assays of deletion strains that do not have a synergistic phenotype.
Figure 6
Discovery of L. pneumophila gene combinations critical for growth in A. castellanii.
(A) Correlation grid plotting all pairwise spacer combinations. The number of times two unique spacers were present in the sequenced array library is indicated by color-coding, ranging from light (rare) to dark (abundant). Black boxes, same spacer. (B) Flowchart going from the number of total spacer combinations to the number of critical combinations identified. (C) Histogram indicating the plex-ness of all critical spacer combinations identified. (D) A list of virulence-critical combinations of genes and the corresponding fold reduction observed for each in each experiment. Library 1 results are shown above Library 2 results. (E) Intracellular growth assays of strains+pMME2400 with deletions in the indicated genes. Results are given as fold growth (colony forming units [CFUs] harvested 48 hr post infection [hpi] vs. 2 hpi) compared to that of Lp02+pMME2400 (WT). Horizontal bars indicate mean fold growth of the deletion strain vs. WT strain for two or more experimental replicates. BD, below detection (arbitrarily set to –100,000). (F) Intracellular growth assays of strains with deletions in the indicated genes having been identified in single-round CRISPRi experiments.
Additional files
-
Supplementary file 1
Supplementary tables for construct design and gene combinations discovered by CRISPRi.
(a) 10-plex construct descriptions and list of spacer sequences used in this study. (b) List of transmembrane effector (TME) genes and corresponding spacer sequences targeting them. (c) PacBio long-read sequencing metrics. (d) Read counts and calculations for all constructs tested in the experiments. (e) Subset data for all hits identified in U937 macrophage infection experiments. (f) Subset data for all critical gene combinations identified in U937 macrophage infection experiments. (g) Subset data for all hits identified in A. castellanii infection experiments. (h) Subset data for all critical gene combinations identified in A. castellanii infection experiments. (i) Hits found in single-round CRISPRi experiments. (j) Strains and plasmids used in this study.
- https://cdn.elifesciences.org/articles/86903/elife-86903-supp1-v1.xlsx
-
MDAR checklist
- https://cdn.elifesciences.org/articles/86903/elife-86903-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)
A randomized multiplex CRISPRi-Seq approach for the identification of critical combinations of genes
eLife 12:RP86903.
https://doi.org/10.7554/eLife.86903.3
