Approaches to querying bacterial genomes with transposon-insertion sequencing

Clinical and laboratory-induced colistin-resistance mechanisms in Acinetobacter baumannii

The increasing incidence and emergence of multi-drug resistant (MDR) Acinetobacter baumannii has become a major global health concern. Colistin is a historic antimicrobial that has become commonly used as a treatment for MDR A. baumannii infections. The increase in colistin usage has been mirrored by an increase in colistin resistance. We aimed to identify the mechanisms associated with colistin resistance in A. baumannii using multiple high-throughput-sequencing technologies, including transposon-directed insertion site sequencing (TraDIS), RNA sequencing (RNAseq) and whole-genome sequencing (WGS) to investigate the genotypic changes of colistin resistance in A. baumannii. Using TraDIS, we found that genes involved in drug efflux (adeIJK), and phospholipid (mlaC, mlaF and mlaD) and lipooligosaccharide synthesis (lpxC and lpsO) were required for survival in sub-inhibitory concentrations of colistin. Transcriptomic (RNAseq) analysis revealed that expression of genes encoding efflux proteins (adeI, adeC, emrB, mexB and macAB) was enhanced in in vitro generated colistin-resistant strains. WGS of these organisms identified disruptions in genes involved in lipid A (lpxC) and phospholipid synthesis (mlaA), and in the baeS/R two-component system (TCS). We additionally found that mutations in the pmrB TCS genes were the primary colistin-resistance-associated mechanisms in three Vietnamese clinical colistin-resistant A. baumannii strains. Our results outline the entire range of mechanisms employed in A. baumannii for resistance against colistin, including drug extrusion and the loss of lipid A moieties by gene disruption or modification.

Basal transcription profiles of the rhamnose-inducible promoter PLRA3 and the development of efficient PLRA3-based systems for markerless gene deletion and a mutant library in Pichia pastoris

An ideal inducible promoter presents inducibility with an inducer and no basal transcription without inducer. Previous studies have shown that P_LRA3 in Pichia pastoris is a strong rhamnose-inducible promoter for driving the industrial production of recombinant proteins. However, another important profile of P_LRA3, the basal transcription, was not investigated yet. In this study, the basal transcription of P_LRA3 was assessed according to the profiles of two test strains grown in media lacking rhamnose: (1) the production of secretory β-galactosidase in P. pastoris GS115/P_LRA3-LacB, in which lacB expression was regulated by P_LRA3, and (2) growth in P. pastoris GS115/P_LRA3-MazF, in which the expression of mazF, which encodes an intracellular toxic protein, was controlled by P_LRA3. Analyses revealed low β-galactosidase production and non-obviously inhibited growth of the test strains, which suggests that there was a low basal transcription level of P_LRA3 when rhamnose was absent. Thus, P_LRA3 was an excellent candidate for genetic manipulation in P. pastoris due to its strict regulation, a strong and a low transcriptional activity with and without rhamnose, respectively. Subsequently, two systems were developed based on P_LRA3 in P. pastoris: (1) an efficient markerless gene deletion system for single or multiple genes and (2) a high efficient piggyBac transposase-mediated mutation system for investigating the functions of unknown genes, as well as for the screening of expected mutants.

Application of Whole-Genome Sequencing to Transposon Mutants of Salmonella Heidelberg

Transposons are elements widely dispersed among organisms which are able to move and replicate fragments of genomes. The extensive variability in transposons present in most organisms requires extensive identification and interpretation of the resulting transposon mutants after transposon mutagenesis. However, much of this is reliant on utilizing randomness characteristics of transposon to identify essential genes for the organism of interest. Integration of the transposon mutant approach with commercialized next-generation sequencing (NGS) technology has helped to advance transposon identification by sequencing millions of reads generated from a single run on an NGS platform. Transposon sequencing is defined as a gene sequencing methodology that allows for the identification of nonessential genes and the determination of gene function using a random transposon insertional mutagenesis followed by massively parallel sequencing. The detailed protocol will be outlined in this chapter. The genomic DNA integrated with the transposons is sequenced using an NGS platform in order to determine the location of each mutation.

Querying Legionella Genomes Using Transposition-Sequencing

Transposition-sequencing (Tn-seq) has recently emerged as a powerful technique to query bacterial genomes. Tn-seq can be used to query the bacterial genome with unprecedented resolution, allowing the identification of small genes (e.g., noncoding RNA) that may be missed in conventional screening approaches. Tn-seq can be used to predict genes essential for in vitro growth and to directly identify genetic requirements for survival under multiple conditions. For instance, Tn-seq can be applied to determine the genes, and cellular processes, required to resist an antibacterial treatment or to acquire new resistance genes, to adapt to intracellular life or to compete with other bacteria. Virtually any assay that involves a selection pressure can be used to identify the associated genetic determinants. So far, genome-wide Tn-seq has not been applied to Legionella species. Here, we provide a protocol covering all the different steps to conduct a Tn-seq analysis in L. pneumophila. This includes generating a high-density library of insertional mutants, setting up a selection screen, sequencing the libraries, mapping the insertion sites, and analyzing the data to obtain the list of genes involved in surviving the applied selection.

Implementation of Transposon Mutagenesis in Bifidobacterium

Random transposon mutagenesis allows for relatively rapid, genome-wide surveys to detect genes involved in functional traits, by performing screens of mutant libraries. This approach has been widely applied to identify genes responsible for activities of interest in multiple eukaryote and prokaryote organisms, although most studies on microorganisms have focused on pathogenic and clinically relevant bacteria. In this chapter we describe the implementation of an in vitro Tn5-based transposome strategy to generate a large collection of random mutants in the gut commensal Bifidobacterium breve UCC2003, and discuss considerations when applying this mutagenesis system to other Bifidobacterium species or strains of interest.

A global genomic approach uncovers novel components for twitching motility-mediated biofilm expansion in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an extremely successful pathogen able to cause both acute and chronic infections in a range of hosts, utilizing a diverse arsenal of cell-associated and secreted virulence factors. A major cell-associated virulence factor, the Type IV pilus (T4P), is required for epithelial cell adherence and mediates a form of surface translocation termed twitching motility, which is necessary to establish a mature biofilm and actively expand these biofilms. P. aeruginosa twitching motility-mediated biofilm expansion is a coordinated, multicellular behaviour, allowing cells to rapidly colonize surfaces, including implanted medical devices. Although at least 44 proteins are known to be involved in the biogenesis, assembly and regulation of the T4P, with additional regulatory components and pathways implicated, it is unclear how these components and pathways interact to control these processes. In the current study, we used a global genomics-based random-mutagenesis technique, transposon directed insertion-site sequencing (TraDIS), coupled with a physical segregation approach, to identify all genes implicated in twitching motility-mediated biofilm expansion in P. aeruginosa. Our approach allowed identification of both known and novel genes, providing new insight into the complex molecular network that regulates this process in P. aeruginosa. Additionally, our data suggest that the flagellum-associated gene products have a differential effect on twitching motility, based on whether components are intra- or extracellular. Overall the success of our TraDIS approach supports the use of this global genomic technique for investigating virulence genes in bacterial pathogens.

Standardization of G. mellonella Larvae to Provide Reliable and Reproducible Results in the Study of Fungal Pathogens

In the past decade, Galleria mellonella (wax moth) larvae have become widely used as a non-mammalian infection model. However, the full potential of this infection model has yet to be realised, limited by the variable quality of larvae used and the lack of standardised procedures. Here, we review larvae suitable for research, protocols for dosing larvae, and methods for scoring illness in larvae infected with fungal pathogens. The development of standardised protocols for carrying out our experimental work will allow high throughput screens to be developed, changing the way in which we evaluate panels of mutants and strains. It will also enable the in vivo screening of potential antimicrobials at an earlier stage in the research and development cycle.

Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157

Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.

Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157

Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.

PiggyBac transposon-mediated mutagenesis and application in yeast Komagataella phaffii

OBJECTIVE

Around one-fourth of the Komagataella phaffii genes encode hypothetical proteins with unknown functions. However, lack of powerful tools for genetic screening in K. phaffii significantly limits the functional analysis of these unknown genes. Transposon mutagenesis has been utilized as an insertional mutagenesis tool in many other organisms and would be extremely valuable if it could be applied in K. phaffii.

RESULTS

In this study, we investigated in K. phaffii the transposition activity and efficiency of piggyBac (PB) transposon, a DNA transposon from the cabbage looper moth Trichoplusia ni through the integrated-plasmid system. We also designed a binary-plasmid system which could generate stable mutants. Finally we evaluated the quality of this mutagenesis system by a simple screening for functional genes involved in K. phaffii carbon catabolite repression.

CONCLUSIONS

Our results demonstrate that PB-mediated mutagenesis could be a feasible and useful tool for functional gene screening in K. phaffii.

Regulatory RNAs in Virulence and Host-Microbe Interactions

Bacterial regulatory RNAs are key players in adaptation to changing environmental conditions and response to diverse cellular stresses. However, while regulatory RNAs of bacterial pathogens have been intensely studied under defined conditions in vitro, characterization of their role during the infection of eukaryotic host organisms is lagging behind. This review summarizes our current understanding of the contribution of the different classes of regulatory RNAs and RNA-binding proteins to bacterial virulence and illustrates their role in infection by reviewing the mechanisms of some prominent representatives of each class. Emerging technologies are described that bear great potential for global, unbiased studies of virulence-related RNAs in bacterial model and nonmodel pathogens in the future. The review concludes by deducing common principles of RNA-mediated gene expression control of virulence programs in different pathogens, and by defining important open questions for upcoming research in the field.

Transposon Insertion Site Sequencing for Synthetic Lethal Screening

Transposon insertion site sequencing (TIS) permits genome-wide, quantitative fitness assessment of individual genomic loci. In addition to the identification of essential genes in given growth conditions, TIS enables the elucidation of genetic networks such as synthetic lethal or suppressor gene combinations. Therefore, TIS becomes an exceptionally powerful tool for the high-throughput determination of genotype-phenotype relationships in bacteria. Here, we describe a protocol for the generation of high-density transposon insertion libraries and subsequent preparation of DNA samples for Illumina sequencing using the Gram-negative bacterium Vibrio cholerae as an example.

Magic Pools: Parallel Assessment of Transposon Delivery Vectors in Bacteria

Molecular genetics is indispensable for interrogating the physiology of bacteria. However, the development of a functional genetic system for any given bacterium can be time-consuming. Here, we present a streamlined approach for identifying an effective transposon mutagenesis system for a new bacterium. Our strategy first involves the construction of hundreds of different transposon vector variants, which we term a “magic pool.” The efficacy of each vector in a magic pool is monitored in parallel using a unique DNA barcode that is introduced into each vector design. Using archived DNA “parts,” we next reassemble an effective vector for making a whole-genome transposon mutant library that is suitable for large-scale interrogation of gene function using competitive growth assays. Here, we demonstrate the utility of the magic pool system to make mutant libraries in five genera of bacteria.

Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach for discovering the functions of bacterial genes. However, the development of a suitable TnSeq strategy for a given bacterium can be costly and time-consuming. To meet this challenge, we describe a part-based strategy for constructing libraries of hundreds of transposon delivery vectors, which we term “magic pools.” Within a magic pool, each transposon vector has a different combination of upstream sequences (promoters and ribosome binding sites) and antibiotic resistance markers as well as a random DNA barcode sequence, which allows the tracking of each vector during mutagenesis experiments. To identify an efficient vector for a given bacterium, we mutagenize it with a magic pool and sequence the resulting insertions; we then use this efficient vector to generate a large mutant library. We used the magic pool strategy to construct transposon mutant libraries in five genera of bacteria, including three genera of the phylum Bacteroidetes.

IMPORTANCE Molecular genetics is indispensable for interrogating the physiology of bacteria. However, the development of a functional genetic system for any given bacterium can be time-consuming. Here, we present a streamlined approach for identifying an effective transposon mutagenesis system for a new bacterium. Our strategy first involves the construction of hundreds of different transposon vector variants, which we term a “magic pool.” The efficacy of each vector in a magic pool is monitored in parallel using a unique DNA barcode that is introduced into each vector design. Using archived DNA “parts,” we next reassemble an effective vector for making a whole-genome transposon mutant library that is suitable for large-scale interrogation of gene function using competitive growth assays. Here, we demonstrate the utility of the magic pool system to make mutant libraries in five genera of bacteria.

Molecular phenotyping of infection-associated small non-coding RNAs

Infection is a complicated balance, with both pathogen and host struggling to tilt the result in their favour. Bacterial infection biology has relied on forward genetics for many of its advances, defining phenotype in terms of replication in model systems. However, many known virulence factors fail to produce robust phenotypes, particularly in the systems most amenable to genetic manipulation, such as cell-culture models. This has particularly been limiting for the study of the bacterial regulatory small RNAs (sRNAs) in infection. We argue that new sequencing-based technologies can work around this problem by providing a 'molecular phenotype', defined in terms of the specific transcriptional dysregulation in the infection system induced by gene deletion. We illustrate this using the example of our recent study of the PinT sRNA using dual RNA-seq, that is, simultaneous RNA sequencing of host and pathogen during infection. We additionally discuss how other high-throughput technologies, in particular genetic interaction mapping using transposon insertion sequencing, may be used to further dissect molecular phenotypes. We propose a strategy for how high-throughput technologies can be integrated in the study of non-coding regulators as well as bacterial virulence factors, enhancing our ability to rapidly generate hypotheses with regards to their function.This article is part of the themed issue 'The new bacteriology'.

Genome-wide identification of tolerance mechanisms toward p-coumaric acid in Pseudomonas putida

The soil bacterium Pseudomonas putida KT2440 has gained increasing biotechnological interest due to its ability to tolerate different types of stress. Here, the tolerance of P. putida KT2440 toward eleven toxic chemical compounds was investigated. P. putida was found to be significantly more tolerant toward three of the eleven compounds when compared to Escherichia coli. Increased tolerance was for example found toward p‐coumaric acid, an interesting precursor for polymerization with a significant industrial relevance. The tolerance mechanism was therefore investigated using the genome‐wide approach, Tn‐seq. Libraries containing a large number of miniTn5‐Km transposon insertion mutants were grown in the presence and absence of p‐coumaric acid, and the enrichment or depletion of mutants was quantified by high‐throughput sequencing. Several genes, including the ABC transporter Ttg2ABC and the cytochrome c maturation system (ccm), were identified to play an important role in the tolerance toward p‐coumaric acid of this bacterium. Most of the identified genes were involved in membrane stability, suggesting that tolerance toward p‐coumaric acid is related to transport and membrane integrity.

RNA-seq and Tn-seq reveal fitness determinants of vancomycin-resistant Enterococcus faecium during growth in human serum

Background

The Gram-positive bacterium Enterococcus faecium is a commensal of the human gastrointestinal tract and a frequent cause of bloodstream infections in hospitalized patients. The mechanisms by which E. faecium can survive and grow in blood during an infection have not yet been characterized. Here, we identify genes that contribute to growth of E. faecium in human serum through transcriptome profiling (RNA-seq) and a high-throughput transposon mutant library sequencing approach (Tn-seq).

Results

We first sequenced the genome of E. faecium E745, a vancomycin-resistant clinical isolate, using a combination of short- and long read sequencing, revealing a 2,765,010 nt chromosome and 6 plasmids, with sizes ranging between 9.3 kbp and 223.7 kbp. We then compared the transcriptome of E. faecium E745 during exponential growth in rich medium and in human serum by RNA-seq. This analysis revealed that 27.8% of genes on the E. faecium E745 genome were differentially expressed in these two conditions. A gene cluster with a role in purine biosynthesis was among the most upregulated genes in E. faecium E745 upon growth in serum. The E. faecium E745 transposon mutant library was then used to identify genes that were specifically required for growth of E. faecium in serum. Genes involved in de novo nucleotide biosynthesis (including pyrK_2, pyrF, purD, purH) and a gene encoding a phosphotransferase system subunit (manY_2) were thus identified to be contributing to E. faecium growth in human serum. Transposon mutants in pyrK_2, pyrF, purD, purH and manY_2 were isolated from the library and their impaired growth in human serum was confirmed. In addition, the pyrK_2 and manY_2 mutants were tested for their virulence in an intravenous zebrafish infection model and exhibited significantly attenuated virulence compared to E. faecium E745.

Conclusions

Genes involved in carbohydrate metabolism and nucleotide biosynthesis of E. faecium are essential for growth in human serum and contribute to the pathogenesis of this organism. These genes may serve as targets for the development of novel anti-infectives for the treatment of E. faecium bloodstream infections.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4299-9) contains supplementary material, which is available to authorized users.

Advances and perspectives in computational prediction of microbial gene essentiality

The minimal subset of genes required for cellular growth, survival and viability of an organism are classified as essential genes. Knowledge of essential genes gives insight into the core structure and functioning of a cell. This might lead to more efficient antimicrobial drug discovery, to elucidation of the correlations between genotype and phenotype, and a better understanding of the minimal requirements for a (synthetic) cell. Traditionally, constructing a catalog of essential genes for a given microbe involved costly and time-consuming laboratory experiments. While experimental methods have produced abundant gene essentiality data for model organisms like Escherichia coli and Bacillus subtilis, the knowledge generated cannot automatically be extrapolated to predict essential genes in all bacteria. In addition, essential genes identified in the laboratory are by definition 'conditionally essential', as they are essential under the specified experimental conditions: these might not resemble conditions in the microorganisms' natural habitat(s). Also, large-scale experimental assaying for essential genes is not always feasible because of the time investment required to setup these assays. The ability to rapidly and precisely identify essential genes in silico is therefore important and has great potential for applications in medicine, biotechnology and basic biological research. Here, we review the advances made in the use of computational methods to predict microbial gene essentiality, perspectives for the future of these techniques and the possible practical applications of essential genes.

Time-Resolved Transposon Insertion Sequencing Reveals Genome-Wide Fitness Dynamics during Infection

Transposon insertion sequencing (TIS) is a powerful high-throughput genetic technique that is transforming functional genomics in prokaryotes, because it enables genome-wide mapping of the determinants of fitness. However, current approaches for analyzing TIS data assume that selective pressures are constant over time and thus do not yield information regarding changes in the genetic requirements for growth in dynamic environments (e.g., during infection). Here, we describe structured analysis of TIS data collected as a time series, termed pattern analysis of conditional essentiality (PACE). From a temporal series of TIS data, PACE derives a quantitative assessment of each mutant’s fitness over the course of an experiment and identifies mutants with related fitness profiles. In so doing, PACE circumvents major limitations of existing methodologies, specifically the need for artificial effect size thresholds and enumeration of bacterial population expansion. We used PACE to analyze TIS samples of Edwardsiella piscicida (a fish pathogen) collected over a 2-week infection period from a natural host (the flatfish turbot). PACE uncovered more genes that affect E. piscicida’s fitness in vivo than were detected using a cutoff at a terminal sampling point, and it identified subpopulations of mutants with distinct fitness profiles, one of which informed the design of new live vaccine candidates. Overall, PACE enables efficient mining of time series TIS data and enhances the power and sensitivity of TIS-based analyses.

Transposon insertion sequencing (TIS) enables genome-wide mapping of the genetic determinants of fitness, typically based on observations at a single sampling point. Here, we move beyond analysis of endpoint TIS data to create a framework for analysis of time series TIS data, termed pattern analysis of conditional essentiality (PACE). We applied PACE to identify genes that contribute to colonization of a natural host by the fish pathogen Edwardsiella piscicida. PACE uncovered more genes that affect E. piscicida’s fitness in vivo than were detected using a terminal sampling point, and its clustering of mutants with related fitness profiles informed design of new live vaccine candidates. PACE yields insights into patterns of fitness dynamics and circumvents major limitations of existing methodologies. Finally, the PACE method should be applicable to additional “omic” time series data, including screens based on clustered regularly interspaced short palindromic repeats with Cas9 (CRISPR/Cas9).

Genome Scale Mutational Analysis of Geobacter sulfurreducens Reveals Distinct Molecular Mechanisms for Respiration and Sensing of Poised Electrodes versus Fe(III) Oxides

Geobacter sulfurreducens generates electrical current by coupling intracellular oxidation of organic acids to the reduction of proteins on the cell surface that are able to interface with electrodes. This ability is attributed to the bacterium's capacity to respire other extracellular electron acceptors that require contact, such as insoluble metal oxides. To directly investigate the genetic basis of electrode-based respiration, we constructed Geobacter sulfurreducens transposon-insertion sequencing (Tn-Seq) libraries for growth, with soluble fumarate or an electrode as the electron acceptor. Libraries with >33,000 unique insertions and an average of 9 insertions/kb allowed an assessment of each gene's fitness in a single experiment. Mutations in 1,214 different genomic features impaired growth with fumarate, and the significance of 270 genes unresolved by annotation due to the presence of one or more functional homologs was determined. Tn-Seq analysis of −0.1 V versus standard hydrogen electrode (SHE) electrode-grown cells identified mutations in a subset of genes encoding cytochromes, processing systems for proline-rich proteins, sensory networks, extracellular structures, polysaccharides, and metabolic enzymes that caused at least a 50% reduction in apparent growth rate. Scarless deletion mutants of select genes identified via Tn-Seq revealed a new putative porin-cytochrome conduit complex (extABCD) crucial for growth with electrodes, which was not required for Fe(III) oxide reduction. In addition, four mutants lacking components of a putative methyl-accepting chemotaxis–cyclic dinucleotide sensing network (esnABCD) were defective in electrode colonization but grew normally with Fe(III) oxides. These results suggest that G. sulfurreducens possesses distinct mechanisms for recognition, colonization, and reduction of electrodes compared to Fe(III) oxides.

IMPORTANCE Since metal oxide electron acceptors are insoluble, one hypothesis is that cells sense and reduce metals using the same molecular mechanisms used to form biofilms on electrodes and produce electricity. However, by simultaneously comparing thousands of Geobacter sulfurreducens transposon mutants undergoing electrode-dependent respiration, we discovered new cytochromes and chemosensory proteins supporting growth with electrodes that are not required for metal respiration. This supports an emerging model where G. sulfurreducens recognizes surfaces and forms conductive biofilms using mechanisms distinct from those used for growth with metal oxides. These findings provide a possible explanation for studies that correlate electricity generation with syntrophic interspecies electron transfer by Geobacter and reveal many previously unrecognized targets for engineering this useful capability in other organisms.

Genome-wide discovery of novel M1T1 group A streptococcal determinants important for fitness and virulence during soft-tissue infection

The Group A Streptococcus remains a significant human pathogen causing a wide array of disease ranging from self-limiting to life-threatening invasive infections. Epithelium (skin or throat) colonization with progression to the subepithelial tissues is the common step in all GAS infections. Here, we used transposon-sequencing (Tn-seq) to define the GAS 5448 genetic requirements for in vivo fitness in subepithelial tissue. A near-saturation transposon library of the M1T1 GAS 5448 strain was injected subcutaneously into mice, producing suppurative inflammation at 24 h that progressed to prominent abscesses with tissue necrosis at 48 h. The library composition was monitored en masse by Tn-seq and ratios of mutant abundance comparing the output (12, 24 and 48 h) versus input (T₀) mutant pools were calculated for each gene. We identified a total of 273 subcutaneous fitness (scf) genes with 147 genes (55 of unknown function) critical for the M1T1 GAS 5448 fitness in vivo; and 126 genes (53 of unknown function) potentially linked to in vivo fitness advantage. Selected scf genes were validated in competitive subcutaneous infection with parental 5448. Two uncharacterized genes, scfA and scfB, encoding putative membrane-associated proteins and conserved among Gram-positive pathogens, were further characterized. Defined scfAB mutants in GAS were outcompeted by wild type 5448 in vivo, attenuated for lesion formation in the soft tissue infection model and dissemination to the bloodstream. We hypothesize that scfAB play an integral role in enhancing adaptation and fitness of GAS during localized skin infection, and potentially in propagation to other deeper host environments.

The WHO ranks the Group A Streptococcus (GAS) in the top 10 leading causes of morbidity and mortality from infectious diseases worldwide. GAS is a strict human pathogen causing both benign superficial infections as well as life-threatening invasive diseases. All GAS infections begin by colonization of an epithelium (throat or skin) followed by propagation into subepithelial tissues. The genetic requirements for M1T1 GAS 5448 within this niche were interrogated by in vivo transposon sequencing (Tn-seq), identifying 273 subcutaneous fitness (scf) genes with 108 of those previously of “unknown function”. Two yet uncharacterized genes, scfA and scfB, were shown to be critical during GAS 5448 soft tissue infection and dissemination into the bloodstream. Thus, this study improves the functional annotation of the GAS genome, providing new insights into GAS pathophysiology and enhancing the development of novel GAS therapeutics.

Role of sapA and yfgA in Susceptibility to Antibody-Mediated Complement-Dependent Killing and Virulence of Salmonella enterica Serovar Typhimurium

The ST313 pathovar of Salmonella enterica serovar Typhimurium contributes to a high burden of invasive disease among African infants and HIV-infected adults. It is characterized by genome degradation (loss of coding capacity) and has increased resistance to antibody-dependent complement-mediated killing compared with enterocolitis-causing strains of S. Typhimurium. Vaccination is an attractive disease-prevention strategy, and leading candidates focus on the induction of bactericidal antibodies. Antibody-resistant strains arising through further gene deletion could compromise such a strategy. Exposing a saturating transposon insertion mutant library of S. Typhimurium to immune serum identified a repertoire of S. Typhimurium genes that, when interrupted, result in increased resistance to serum killing. These genes included several involved in bacterial envelope biogenesis, protein translocation, and metabolism. We generated defined mutant derivatives using S. Typhimurium SL1344 as the host. Based on their initial levels of enhanced resistance to killing, yfgA and sapA mutants were selected for further characterization. The S. Typhimurium yfgA mutant lost the characteristic Salmonella rod-shaped appearance, exhibited increased sensitivity to osmotic and detergent stress, lacked very long lipopolysaccharide, was unable to invade enterocytes, and demonstrated decreased ability to infect mice. In contrast, the S. Typhimurium sapA mutants had similar sensitivity to osmotic and detergent stress and lipopolysaccharide profile and an increased ability to infect enterocytes compared with the wild type, but it had no increased ability to cause in vivo infection. These findings indicate that increased resistance to antibody-dependent complement-mediated killing secondary to genetic deletion is not necessarily accompanied by increased virulence and suggest the presence of different mechanisms of antibody resistance.

Small RNAs in Bacterial Virulence and Communication

Bacterial pathogens must endure or adapt to different environments and stresses during transmission and infection. Posttranscriptional gene expression control by regulatory RNAs, such as small RNAs and riboswitches, is now considered central to adaptation in many bacteria, including pathogens. The study of RNA-based regulation (riboregulation) in pathogenic species has provided novel insight into how these bacteria regulate virulence gene expression. It has also uncovered diverse mechanisms by which bacterial small RNAs, in general, globally control gene expression. Riboregulators as well as their targets may also prove to be alternative targets or provide new strategies for antimicrobials. In this article, we present an overview of the general mechanisms that bacteria use to regulate with RNA, focusing on examples from pathogens. In addition, we also briefly review how deep sequencing approaches have aided in opening new perspectives in small RNA identification and the study of their functions. Finally, we discuss examples of riboregulators in two model pathogens that control virulence factor expression or survival-associated phenotypes, such as stress tolerance, biofilm formation, or cell-cell communication, to illustrate how riboregulation factors into regulatory networks in bacterial pathogens.

Natural resource landscapes of a marine bacterium reveal distinct fitness-determining genes across the genome

Heterotrophic bacteria exploit diverse microhabitats in the ocean, from particles to transient gradients. Yet the degree to which genes and pathways can contribute to an organism's fitness on such complex and variable natural resource landscapes remains poorly understood. Here, we determine the gene-by-gene fitness of a generalist saprophytic marine bacterium (Vibrio sp. F13 9CS106) on complex resources derived from its natural habitats - copepods (Apocyclops royi) and brown algae (Fucus vesiculosus) - and as reference substrates, glucose and the polysaccharide alginate, derived from brown algal cell walls. We find that resource complexity strongly buffers fitness costs of mutations, and that anabolic rather than catabolic pathways are more stringently required, likely due to functional redundancy in the latter. Moreover, while carbohydrate-rich algae requires several synthesis pathways, protein-rich Apocyclops does not, suggesting this ancestral habitat for Vibrios is a replete medium with metabolically redundant substrates. We also identify a candidate fitness trade-off for algal colonization: deletion of mshA increases mutant fitness. Our results demonstrate that gene fitness depends on habitat composition, and suggest that this generalist uses distinct resources in different natural habitats. The results further indicate that substrate replete conditions may lead to relatively relaxed selection on catabolic genes.

Defining the ABC of gene essentiality in streptococci

BACKGROUND

Utilising next generation sequencing to interrogate saturated bacterial mutant libraries provides unprecedented information for the assignment of genome-wide gene essentiality. Exposure of saturated mutant libraries to specific conditions and subsequent sequencing can be exploited to uncover gene essentiality relevant to the condition. Here we present a barcoded transposon directed insertion-site sequencing (TraDIS) system to define an essential gene list for Streptococcus equi subsp. equi, the causative agent of strangles in horses, for the first time. The gene essentiality data for this group C Streptococcus was compared to that of group A and B streptococci.

RESULTS

Six barcoded variants of pGh9:ISS1 were designed and used to generate mutant libraries containing between 33,000-66,000 unique mutants. TraDIS was performed on DNA extracted from each library and data were analysed separately and as a combined master pool. Gene essentiality determined that 19.5% of the S. equi genome was essential. Gene essentialities were compared to those of group A and group B streptococci, identifying concordances of 90.2% and 89.4%, respectively and an overall concordance of 83.7% between the three species.

CONCLUSIONS

The use of barcoded pGh9:ISS1 to generate mutant libraries provides a highly useful tool for the assignment of gene function in S. equi and other streptococci. The shared essential gene set of group A, B and C streptococci provides further evidence of the close genetic relationships between these important pathogenic bacteria. Therefore, the ABC of gene essentiality reported here provides a solid foundation towards reporting the functional genome of streptococci.

Physical enrichment of transposon mutants from saturation mutant libraries using the TraDISort approach

Transposon-insertion sequencing methods are finding their way into the molecular toolbox of many fields of microbiology. These methods can identify the genomic locations and density of transposon insertions in saturated transposon mutant libraries and can be used to make inferences on gene function. For example, where no insertions or very few insertions are identified within a gene in a mutant library grown under permissive conditions, the gene may be essential. Furthermore, where mutations are enriched or lost in a gene after passaging the library through a selective process, the gene is likely to be involved in phenotypes linked to the process. Typically, a fitness based selection such as a stress condition is used in these experiments and the processed sequencing data are used to identify genes required for fitness under the selection. Our research team recently expanded the utility of the transposon directed insertion sequencing (TraDIS) method by applying a physical separation of a transposon mutant library mediated by fluorescence activated cell sorting, rather than a fitness-based selection. This approach, which we have named "TraDISort" is significant because it allows the study of phenotypes that are not linked to cell survival. The TraDISort approach has a broad range of future applications, in drug development, metabolic engineering and in studies of basic bacterial cell physiology.

The secondary resistome of multidrug-resistant Klebsiella pneumoniae

Klebsiella pneumoniae causes severe lung and bloodstream infections that are difficult to treat due to multidrug resistance. We hypothesized that antimicrobial resistance can be reversed by targeting chromosomal non-essential genes that are not responsible for acquired resistance but essential for resistant bacteria under therapeutic concentrations of antimicrobials. Conditional essentiality of individual genes to antimicrobial resistance was evaluated in an epidemic multidrug-resistant clone of K. pneumoniae (ST258). We constructed a high-density transposon mutant library of >430,000 unique Tn5 insertions and measured mutant depletion upon exposure to three clinically relevant antimicrobials (colistin, imipenem or ciprofloxacin) by Transposon Directed Insertion-site Sequencing (TraDIS). Using this high-throughput approach, we defined three sets of chromosomal non-essential genes essential for growth during exposure to colistin (n = 35), imipenem (n = 1) or ciprofloxacin (n = 1) in addition to known resistance determinants, collectively termed the "secondary resistome". As proof of principle, we demonstrated that inactivation of a non-essential gene not previously found linked to colistin resistance (dedA) restored colistin susceptibility by reducing the minimum inhibitory concentration from 8 to 0.5 μg/ml, 4-fold below the susceptibility breakpoint (S ≤ 2 μg/ml). This finding suggests that the secondary resistome is a potential target for developing antimicrobial "helper" drugs that restore the efficacy of existing antimicrobials.

A Noise Trimming and Positional Significance of Transposon Insertion System to Identify Essential Genes in Yersinia pestis

Massively parallel sequencing technology coupled with saturation mutagenesis has provided new and global insights into gene functions and roles. At a simplistic level, the frequency of mutations within genes can indicate the degree of essentiality. However, this approach neglects to take account of the positional significance of mutations - the function of a gene is less likely to be disrupted by a mutation close to the distal ends. Therefore, a systematic bioinformatics approach to improve the reliability of essential gene identification is desirable. We report here a parametric model which introduces a novel mutation feature together with a noise trimming approach to predict the biological significance of Tn5 mutations. We show improved performance of essential gene prediction in the bacterium Yersinia pestis, the causative agent of plague. This method would have broad applicability to other organisms and to the identification of genes which are essential for competitiveness or survival under a broad range of stresses.

Genetic Regulation of Virulence and Antibiotic Resistance in Acinetobacter baumannii

Multidrug resistant microorganisms are forecast to become the single biggest challenge to medical care in the 21st century. Over the last decades, members of the genus Acinetobacter have emerged as bacterial opportunistic pathogens, in particular as challenging nosocomial pathogens because of the rapid evolution of antimicrobial resistances. Although we lack fundamental biological insight into virulence mechanisms, an increasing number of researchers are working to identify virulence factors and to study antibiotic resistance. Here, we review current knowledge regarding the regulation of virulence genes and antibiotic resistance in Acinetobacter baumannii. A survey of the two-component systems AdeRS, BaeSR, GacSA and PmrAB explains how each contributes to antibiotic resistance and virulence gene expression, while BfmRS regulates cell envelope structures important for pathogen persistence. A. baumannii uses the transcription factors Fur and Zur to sense iron or zinc depletion and upregulate genes for metal scavenging as a critical survival tool in an animal host. Quorum sensing, nucleoid-associated proteins, and non-classical transcription factors such as AtfA and small regulatory RNAs are discussed in the context of virulence and antibiotic resistance.

The Use of Transposon Insertion Sequencing to Interrogate the Core Functional Genome of the Legume Symbiont Rhizobium leguminosarum

The free-living legume symbiont Rhizobium leguminosarum is of significant economic value because of its ability to provide fixed nitrogen to globally important leguminous food crops, such as peas and lentils. Discovery based research into the genetics and physiology of R. leguminosarum provides the foundational knowledge necessary for understanding the bacterium's complex lifestyle, necessary for augmenting its use in an agricultural setting. Transposon insertion sequencing (INSeq) facilitates high-throughput forward genetic screening at a genomic scale to identify individual genes required for growth in a specific environment. In this study we applied INSeq to screen the genome of R. leguminosarum bv. viciae strain 3841 (RLV3841) for genes required for growth on minimal mannitol containing medium. Results from this study were contrasted with a prior INSeq experiment screened on peptide rich media to identify a common set of functional genes necessary for basic physiology. Contrasting the two growth conditions indicated that approximately 10% of the chromosome was required for growth, under both growth conditions. Specific genes that were essential to singular growth conditions were also identified. Data from INSeq screening on mannitol as a sole carbon source were used to reconstruct a metabolic map summarizing growth impaired phenotypes observed in the Embden-Meyerhof-Parnas pathway, Entner-Doudoroff pathway, pentose phosphate pathway, and tricarboxylic acid cycle. This revealed the presence of mannitol dependent and independent metabolic pathways required for growth, along with identifying metabolic steps with isozymes or possible carbon flux by-passes. Additionally, genes were identified on plasmids pRL11 and pRL12 that are likely to encode functional activities important to the central physiology of RLV3841.

Essential Genes for In Vitro Growth of the Endophyte Herbaspirillum seropedicae SmR1 as Revealed by Transposon Insertion Site Sequencing

The interior of plants contains microorganisms (referred to as endophytes) that are distinct from those present at the root surface or in the surrounding soil. Herbaspirillum seropedicae strain SmR1, belonging to the betaproteobacteria, is an endophyte that colonizes crops, including rice, maize, sugarcane, and sorghum. Different approaches have revealed genes and pathways regulated during the interactions of H. seropedicae with its plant hosts. However, functional genomic analysis of transposon (Tn) mutants has been hampered by the lack of genetic tools. Here we successfully employed a combination of in vivo high-density mariner Tn mutagenesis and targeted Tn insertion site sequencing (Tn-seq) in H. seropedicae SmR1. The analysis of multiple gene-saturating Tn libraries revealed that 395 genes are essential for the growth of H. seropedicae SmR1 in tryptone-yeast extract medium. A comparative analysis with the Database of Essential Genes (DEG) showed that 25 genes are uniquely essential in H. seropedicae SmR1. The Tn mutagenesis protocol developed and the gene-saturating Tn libraries generated will facilitate elucidation of the genetic mechanisms of the H. seropedicae endophytic lifestyle.

IMPORTANCE

A focal point in the study of endophytes is the development of effective biofertilizers that could help to reduce the input of agrochemicals in croplands. Besides the ability to promote plant growth, a good biofertilizer should be successful in colonizing its host and competing against the native microbiota. By using a systematic Tn-based gene-inactivation strategy and massively parallel sequencing of Tn insertion sites (Tn-seq), it is possible to study the fitness of thousands of Tn mutants in a single experiment. We have applied the combination of these techniques to the plant-growth-promoting endophyte Herbaspirillum seropedicae SmR1. The Tn mutant libraries generated will enable studies into the genetic mechanisms of H. seropedicae-plant interactions. The approach that we have taken is applicable to other plant-interacting bacteria.

DNA Transposition at Work

DNA transposons are defined segments of DNA that are able to move from one genomic location to another. Movement is facilitated by one or more proteins, called the transposase, typically encoded by the mobile element itself. Here, we first provide an overview of the classification of such mobile elements in a variety of organisms. From a mechanistic perspective, we have focused on one particular group of DNA transposons that encode a transposase with a DD(E/D) catalytic domain that is topologically similar to RNase H. For these, a number of three-dimensional structures of transpososomes (transposase-nucleic acid complexes) are available, and we use these to describe the basics of their mechanisms. The DD(E/D) group, in addition to being the largest and most common among all DNA transposases, is the one whose members have been used for a wide variety of genomic applications. Therefore, a second focus of the article is to provide a nonexhaustive overview of transposon applications. Although several non-transposon-based approaches to site-directed genome modifications have emerged in the past decade, transposon-based applications are highly relevant when integration specificity is not sought. In fact, for many applications, the almost-perfect randomness and high frequency of integration make transposon-based approaches indispensable.

No-Go Zones for Mariner Transposition

The property of transposons to randomly insert into target DNA has long been exploited for generalized mutagenesis and forward genetic screens. Newer applications that monitor the relative abundance of each transposon insertion in large libraries of mutants can be used to evaluate the roles in cellular fitness of all genes of an organism, provided that transposition is in fact random across all genes. In a recent article, Kimura and colleagues identified an important exception to the latter assumption [S. Kimura, T. P. Hubbard, B. M. Davis, M. K. Waldor, mBio 7(4):e01351-16, 2016, doi:10.1128/mBio.01351-16]. They provide evidence that the Mariner transposon exhibits locus-specific site preferences in the presence of the histone-like nucleoid structuring protein H-NS. This effect was shown to bias results for important virulence loci in Vibrio cholerae and to result in misidentification of genes involved in growth in vitro. Fortunately, the bulk of this bacterium’s genome was unaffected by this bias, and recognizing the H-NS effect allows filtering to improve the accuracy of the results.

The Nucleoid Binding Protein H-NS Biases Genome-Wide Transposon Insertion Landscapes

Transposon insertion sequencing (TIS; also known as TnSeq) is a potent approach commonly used to comprehensively define the genetic loci that contribute to bacterial fitness in diverse environments. A key presumption underlying analyses of TIS datasets is that loci with a low frequency of transposon insertions contribute to fitness. However, it is not known whether factors such as nucleoid binding proteins can alter the frequency of transposon insertion and thus whether TIS output may systematically reflect factors that are independent of the role of the loci in fitness. Here, we investigated whether the histone-like nucleoid structuring (H-NS) protein, which preferentially associates with AT-rich sequences, modulates the frequency of Mariner transposon insertion in the Vibrio cholerae genome, using comparative analysis of TIS results from wild-type (wt) and Δhns V. cholerae strains. These analyses were overlaid on gene classification based on GC content as well as on extant genome-wide identification of H-NS binding loci. Our analyses revealed a significant dearth of insertions within AT-rich loci in wt V. cholerae that was not apparent in the Δhns insertion library. Additionally, we observed a striking correlation between genetic loci that are overrepresented in the Δhns insertion library relative to their insertion frequency in wt V. cholerae and loci previously found to physically interact with H-NS. Collectively, our findings reveal that factors other than genetic fitness can systematically modulate the frequency of transposon insertions in TIS studies and add a cautionary note to interpretation of TIS data, particularly for AT-rich sequences.

IMPORTANCE

Transposon insertion sequencing (TIS) is often used to assess the relative frequency with which genetic loci can be disrupted, which is taken as an indicator of their importance for bacterial fitness. Here, we report that biological factors other than the relative levels of fitness of insertion mutants can influence TIS output. We found that the presence of the DNA binding protein H-NS, which preferentially recognizes AT-rich sequences, is linked to significant underrepresentation of mutations within AT-rich loci in transposon insertion libraries. Furthermore, there is a marked correspondence between loci bound by H-NS and loci with an increased frequency of disruption in a Δhns insertion library relative to a wt library. Our data suggest that factors other than genetic fitness (e.g., DNA binding proteins such as H-NS) can systematically modulate the frequency of transposon insertions in TIS studies and add a note of caution for interpretation of TIS data.

Lessons from Digestive-Tract Symbioses Between Bacteria and Invertebrates

In most animals, digestive tracts harbor the greatest number of bacteria in the animal that contribute to its health: by aiding in the digestion of nutrients, provisioning essential nutrients and protecting against colonization by pathogens. Invertebrates have been used to enhance our understanding of metabolic processes and microbe-host interactions owing to experimental advantages. This review describes how advances in DNA sequencing technologies have dramatically altered how researchers investigate microbe-host interactions, including 16S rRNA gene surveys, metagenome experiments, and metatranscriptome studies. Advantages and challenges of each of these approaches are described herein. Hypotheses generated through omics studies can be directly tested using site-directed mutagenesis, and findings from transposon studies and site-directed experiments are presented. Finally, unique structural aspects of invertebrate digestive tracts that contribute to symbiont specificity are presented. The combination of omics approaches with genetics and microscopy allows researchers to move beyond correlations to identify conserved mechanisms of microbe-host interactions.

Right on Q: genetics begin to unravel Coxiella burnetii host cell interactions

Invasion of macrophages and replication within an acidic and degradative phagolysosome-like vacuole are essential for disease pathogenesis by Coxiella burnetii, the bacterial agent of human Q fever. Previous experimental constraints imposed by the obligate intracellular nature of Coxiella limited knowledge of pathogen strategies that promote infection. Fortunately, new genetic tools facilitated by axenic culture now allow allelic exchange and transposon mutagenesis approaches for virulence gene discovery. Phenotypic screens have illuminated the critical importance of Coxiella's type 4B secretion system in host cell subversion and discovered genes encoding translocated effector proteins that manipulate critical infection events. Here, we highlight the cellular microbiology and genetics of Coxiella and how recent technical advances now make Coxiella a model organism to study macrophage parasitism.

The design and analysis of transposon insertion sequencing experiments

Transposon insertion sequencing (TIS) is a powerful approach that can be extensively applied to the genome-wide definition of loci that are required for bacterial growth under diverse conditions. However, experimental design choices and stochastic biological processes can heavily influence the results of TIS experiments and affect downstream statistical analysis. In this Opinion article, we discuss TIS experimental parameters and how these factors relate to the benefits and limitations of the various statistical frameworks that can be applied to the computational analysis of TIS data.

Genetic analysis of Vibrio parahaemolyticus intestinal colonization

Vibrio parahaemolyticus is the most common cause of seafood-borne gastroenteritis worldwide and a blight on global aquaculture. This organism requires a horizontally acquired type III secretion system (T3SS2) to infect the small intestine, but knowledge of additional factors that underlie V. parahaemolyticus pathogenicity is limited. We used transposon-insertion sequencing to screen for genes that contribute to viability of V. parahaemolyticus in vitro and in the mammalian intestine. Our analysis enumerated and controlled for the host infection bottleneck, enabling robust assessment of genetic contributions to in vivo fitness. We identified genes that contribute to V. parahaemolyticus colonization of the intestine independent of known virulence mechanisms in addition to uncharacterized components of T3SS2. Our study revealed that toxR, an ancestral locus in Vibrio species, is required for V. parahaemolyticus fitness in vivo and for induction of T3SS2 gene expression. The regulatory mechanism by which V. parahaemolyticus ToxR activates expression of T3SS2 resembles Vibrio cholerae ToxR regulation of distinct virulence elements acquired via lateral gene transfer. Thus, disparate horizontally acquired virulence systems have been placed under the control of this ancestral transcription factor across independently evolved human pathogens.

The TraDIS toolkit: sequencing and analysis for dense transposon mutant libraries

Summary: Transposon insertion sequencing is a high-throughput technique for assaying large libraries of otherwise isogenic transposon mutants providing insight into gene essentiality, gene function and genetic interactions. We previously developed the Transposon Directed Insertion Sequencing (TraDIS) protocol for this purpose, which utilizes shearing of genomic DNA followed by specific PCR amplification of transposon-containing fragments and Illumina sequencing. Here we describe an optimized high-yield library preparation and sequencing protocol for TraDIS experiments and a novel software pipeline for analysis of the resulting data. The Bio-Tradis analysis pipeline is implemented as an extensible Perl library which can either be used as is, or as a basis for the development of more advanced analysis tools. This article can serve as a general reference for the application of the TraDIS methodology.

Availability and implementation: The optimized sequencing protocol is included as supplementary information. The Bio-Tradis analysis pipeline is available under a GPL license at https://github.com/sanger-pathogens/Bio-Tradis

Contact:parkhill@sanger.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Bacterial Hypoxic Responses Revealed as Critical Determinants of the Host-Pathogen Outcome by TnSeq Analysis of Staphylococcus aureus Invasive Infection

Staphylococcus aureus is capable of infecting nearly every organ in the human body. In order to infiltrate and thrive in such diverse host tissues, staphylococci must possess remarkable flexibility in both metabolic and virulence programs. To investigate the genetic requirements for bacterial survival during invasive infection, we performed a transposon sequencing (TnSeq) analysis of S. aureus during experimental osteomyelitis. TnSeq identified 65 genes essential for staphylococcal survival in infected bone and an additional 148 mutants with compromised fitness in vivo. Among the loci essential for in vivo survival was SrrAB, a staphylococcal two-component system previously reported to coordinate hypoxic and nitrosative stress responses in vitro. Healthy bone is intrinsically hypoxic, and intravital oxygen monitoring revealed further decreases in skeletal oxygen concentrations upon S. aureus infection. The fitness of an srrAB mutant during osteomyelitis was significantly increased by depletion of neutrophils, suggesting that neutrophils impose hypoxic and/or nitrosative stresses on invading bacteria. To more globally evaluate staphylococcal responses to changing oxygenation, we examined quorum sensing and virulence factor production in staphylococci grown under aerobic or hypoxic conditions. Hypoxic growth resulted in a profound increase in quorum sensing-dependent toxin production, and a concomitant increase in cytotoxicity toward mammalian cells. Moreover, aerobic growth limited quorum sensing and cytotoxicity in an SrrAB-dependent manner, suggesting a mechanism by which S. aureus modulates quorum sensing and toxin production in response to environmental oxygenation. Collectively, our results demonstrate that bacterial hypoxic responses are key determinants of the staphylococcal-host interaction.

The evolution of early cellular systems viewed through the lens of biological interactions

The minimal cell concept represents a pragmatic approach to the question of how few genes are required to run a cell. This is a helpful way to build a parts-list, and has been more successful than attempts to deduce a minimal gene set for life by inferring the gene repertoire of the last universal common ancestor, as few genes trace back to this hypothetical ancestral state. However, the study of minimal cellular systems is the study of biological outliers where, by practical necessity, coevolutionary interactions are minimized or ignored. In this paper, we consider the biological context from which minimal genomes have been removed. For instance, some of the most reduced genomes are from endosymbionts and are the result of coevolutionary interactions with a host; few such organisms are "free-living." As few, if any, biological systems exist in complete isolation, we expect that, as with modern life, early biological systems were part of an ecosystem, replete with organismal interactions. We favor refocusing discussions of the evolution of cellular systems on processes rather than gene counts. We therefore draw a distinction between a pragmatic minimal cell (an interesting engineering problem), a distributed genome (a system resulting from an evolutionary transition involving more than one cell) and the looser coevolutionary interactions that are ubiquitous in ecosystems. Finally, we consider the distributed genome and coevolutionary interactions between genomic entities in the context of early evolution.

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

TRANSIT--A Software Tool for Himar1 TnSeq Analysis

TnSeq has become a popular technique for determining the essentiality of genomic regions in bacterial organisms. Several methods have been developed to analyze the wealth of data that has been obtained through TnSeq experiments. We developed a tool for analyzing Himar1 TnSeq data called TRANSIT. TRANSIT provides a graphical interface to three different statistical methods for analyzing TnSeq data. These methods cover a variety of approaches capable of identifying essential genes in individual datasets as well as comparative analysis between conditions. We demonstrate the utility of this software by analyzing TnSeq datasets of M. tuberculosis grown on glycerol and cholesterol. We show that TRANSIT can be used to discover genes which have been previously implicated for growth on these carbon sources. TRANSIT is written in Python, and thus can be run on Windows, OSX and Linux platforms. The source code is distributed under the GNU GPL v3 license and can be obtained from the following GitHub repository: https://github.com/mad-lab/transit.

High-throughput bacterial functional genomics in the sequencing era

High-throughput functional genomic technologies are accelerating progress in understanding the diversity of bacterial life and in developing a systems-level understanding of model bacterial organisms. Here we highlight progress in deep-sequencing-based functional genomics, show how whole genome sequencing is enabling phenotyping in organisms recalcitrant to genetic approaches, recount the rapid proliferation of functional genomic approaches to non-growth phenotypes, and discuss how advances are enabling genome-scale resource libraries for many different bacteria.

Rapid mapping of insertional mutations to probe cell wall regulation in Cryptococcus neoformans

Random insertional mutagenesis screens are important tools in microbial genetics studies. Investigators in fungal systems have used the plant pathogen Agrobacterium tumefaciens to create tagged, random mutations for genetic screens in their fungal species of interest through a unique process of trans-kingdom cellular transconjugation. However, identifying the locations of insertion has traditionally required tedious PCR-based methods, limiting the effective throughput of this system. We have developed an efficient genomic sequencing and analysis method (AIM-Seq) to facilitate identification of randomly generated genomic insertions in microorganisms. AIM-Seq combines batch sampling, whole genome sequencing, and a novel bioinformatics pipeline, AIM-HII, to rapidly identify sites of genomic insertion. We have specifically applied this technique to Agrobacterium-mediated transconjugation in the human fungal pathogen Cryptococcus neoformans. With this approach, we have screened a library of C. neoformans cell wall mutants, selecting twenty-seven mutants of interest for analysis by AIM-Seq. We identified thirty-five putative genomic insertions in known and previously unknown regulators of cell wall processes in this pathogenic fungus. We confirmed the relevance of a subset of these by creating independent mutant strains and analyzing resulting cell wall phenotypes. Through our sequence-based analysis of these mutations, we observed "typical" insertions of the Agrobacterium transfer DNA as well as atypical insertion events, including large deletions and chromosomal rearrangements. Initially applied to C. neoformans, this mutant analysis tool can be applied to a wide range of experimental systems and methods of mutagenesis, facilitating future microbial genetic screens.

General and condition-specific essential functions of Pseudomonas aeruginosa

The essential functions of a bacterial pathogen reflect the most basic processes required for its viability and growth, and represent potential therapeutic targets. Most screens for essential genes have assayed a single condition--growth in a rich undefined medium--and thus have not distinguished genes that are generally essential from those that are specific to this particular condition. To help define these classes for Pseudomonas aeruginosa, we identified genes required for growth on six different media, including a medium made from cystic fibrosis patient sputum. The analysis used the Tn-seq circle method to achieve high genome coverage and analyzed more than 1,000,000 unique insertion positions (an average of one insertion every 6.0 bp). We identified 352 general and 199 condition-specific essential genes. A subset of assignments was verified in individual strains with regulated expression alleles. The profile of essential genes revealed that, compared with Escherichia coli, P. aeruginosa is highly vulnerable to mutations disrupting central carbon-energy metabolism and reactive oxygen defenses. These vulnerabilities may arise from the stripped-down architecture of the organism's carbohydrate utilization pathways and its reliance on respiration for energy generation. The essential function profile thus provides fundamental insights into P. aeruginosa physiology as well as identifying candidate targets for new antibacterial agents.