The development of new drug regimens that allow rapid, sterilizing treatment of tuberculosis has been limited by the complexity and time required for genetic manipulations in Mycobacterium tuberculosis. CRISPR interference (CRISPRi) promises to be a robust, easily engineered and scalable platform for regulated gene silencing. However, in M. tuberculosis, the existing Streptococcus pyogenes Cas9-based CRISPRi system is of limited utility because of relatively poor knockdown efficiency and proteotoxicity. To address these limitations, we screened eleven diverse Cas9 orthologues and identified four that are broadly functional for targeted gene knockdown in mycobacteria. The most efficacious of these proteins, the CRISPR1 Cas9 from Streptococcus thermophilus (dCas9Sth1), typically achieves 20- to 100-fold knockdown of endogenous gene expression with minimal proteotoxicity. In contrast to other CRISPRi systems, dCas9Sth1-mediated gene knockdown is robust when targeted far from the transcriptional start site, thereby allowing high-resolution dissection of gene function in the context of bacterial operons. We demonstrate the utility of this system by addressing persistent controversies regarding drug synergies in the mycobacterial folate biosynthesis pathway. We anticipate that the dCas9Sth1 CRISPRi system will have broad utility for functional genomics, genetic interaction mapping and drug-target profiling in M. tuberculosis.

For decades, identifying the regions of a bacterial chromosome that are necessary for viability has relied on mapping integration sites in libraries of random transposon mutants to find loci that are unable to sustain insertion. To date, these studies have analyzed subsaturated libraries, necessitating the application of statistical methods to estimate the likelihood that a gap in transposon coverage is the result of biological selection and not the stochasticity of insertion. As a result, the essentiality of many genomic features, particularly small ones, could not be reliably assessed. We sought to overcome this limitation by creating a completely saturated transposon library in Mycobacterium tuberculosis. In assessing the composition of this highly saturated library by deep sequencing, we discovered that a previously unknown sequence bias of the Himar1 element rendered approximately 9% of potential TA dinucleotide insertion sites less permissible for insertion. We used a hidden Markov model of essentiality that accounted for this unanticipated bias, allowing us to confidently evaluate the essentiality of features that contained as few as 2 TA sites, including open reading frames (ORF), experimentally identified noncoding RNAs, methylation sites, and promoters. In addition, several essential regions that did not correspond to known features were identified, suggesting uncharacterized functions that are necessary for growth. This work provides an authoritative catalog of essential regions of the M. tuberculosis genome and a statistical framework for applying saturating mutagenesis to other bacteria. IMPORTANCE: Sequencing of transposon-insertion mutant libraries has become a widely used tool for probing the functions of genes under various conditions. The Himar1 transposon is generally believed to insert with equal probabilities at all TA dinucleotides, and therefore its absence in a mutant library is taken to indicate biological selection against the corresponding mutant. Through sequencing of a saturated Himar1 library, we found evidence that TA dinucleotides are not equally permissive for insertion. The insertion bias was observed in multiple prokaryotes and influences the statistical interpretation of transposon insertion (TnSeq) data and characterization of essential genomic regions. Using these insights, we analyzed a fully saturated TnSeq library for M. tuberculosis, enabling us to generate a comprehensive catalog of in vitro essentiality, including ORFs smaller than those found in any previous study, small (noncoding) RNAs (sRNAs), promoters, and other genomic features.