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.

Nitric oxide (NO) contributes to protection from tuberculosis (TB). It is generally assumed that this protection is due to direct inhibition of Mycobacterium tuberculosis (Mtb) growth, which prevents subsequent pathological inflammation. In contrast, we report NO primarily protects mice by repressing an interleukin-1 and 12/15-lipoxygenase dependent neutrophil recruitment cascade that promotes bacterial replication. Using Mtb mutants as indicators of the pathogen's environment, we inferred that granulocytic inflammation generates a nutrient-replete niche that supports Mtb growth. Parallel clinical studies indicate that a similar inflammatory pathway promotes TB in patients. The human 12/15 lipoxygenase ortholog, ALOX12, is expressed in cavitary TB lesions, the abundance of its products correlate with the number of airway neutrophils and bacterial burden, and a genetic polymorphism that increases ALOX12 expression is associated with TB risk. These data suggest that Mtb exploits neutrophilic inflammation to preferentially replicate at sites of tissue damage that promote contagion.

Identifying genomic elements required for viability is central to our understanding of the basic physiology of bacterial pathogens. Recently, the combination of high-density mutagenesis and deep sequencing has allowed for the identification of required and conditionally required genes in many bacteria. Genes, however, make up only a part of the complex genomes of important bacterial pathogens. Here, we use an unbiased analysis to comprehensively identify genomic regions, including genes, domains, and intergenic elements, required for the optimal growth of Mycobacterium tuberculosis, a major global health pathogen. We found that several proteins jointly contain both domains required for optimal growth and domains that are dispensable. In addition, many non-coding regions, including regulatory elements and non-coding RNAs, are critical for mycobacterial growth. Our analysis shows that the genetic requirements for growth are more complex than can be appreciated using gene-centric analysis.

SUMMARY M. tuberculosis (Mtb) survives a hostile environment within the host that is shaped in part by oxidative stress. The mechanisms used by Mtb to resist these stresses remain ill-defined because the complex combination of oxidants generated by host immunity is difficult to accurately recapitulate in vitro. We performed a genome-wide genetic interaction screen to comprehensively delineate oxidative stress resistance pathways necessary for Mtb to resist oxidation during infection. Our analysis predicted functional relationships between the superoxide-detoxifying enzyme (SodA), an integral membrane protein (DoxX), and a predicted thiol-oxidoreductase (SseA). Consistent with that, SodA, DoxX and SseA form a membrane-associated oxidoreductase complex (MRC) that physically links radical detoxification with cytosolic thiol homeostasis. Loss of any MRC component correlated with defective recycling of mycothiol and accumulation of cellular oxidative damage. This previously uncharacterized coordination between oxygen radical detoxification and thiol homeostasis is required to overcome the oxidative environment Mtb encounters in the host.

Transposon sequencing (TnSeq) is a next-generation deep sequencing-based method to quantitatively assess the composition of complex mutant transposon libraries after pressure from selection. Although this method can be used for any organism in which transposon mutagenesis is possible, this chapter describes its use in Mycobacterium tuberculosis. More specifically, the methods for generating complex libraries through transposon mutagenesis, design of selective pressure, extraction of genomic DNA, amplification and quantification of transposon insertions through next-generation deep sequencing are covered. Determining gene essentiality and statistical analysis on data collected are also discussed.