Antigen Discovery for Next-Generation Pertussis Vaccines Using Immunoproteomics and Transposon-Directed Insertion Sequencing

BACKGROUND

Despite high vaccination rates, the United States has experienced a resurgence in reported cases of pertussis after switching to the acellular pertussis vaccine, indicating a need for improved vaccines that enhance infection control.

METHODS

Bordetella pertussis antigens recognized by convalescent-baboon serum and nasopharyngeal wash were identified by immunoproteomics and their subcellular localization predicted. Genes essential or important for persistence in the baboon airway were identified by transposon-directed insertion-site sequencing (TraDIS) analysis.

RESULTS

In total, 314 B. pertussis antigens were identified by convalescent baboon serum and 748 by nasopharyngeal wash. Thirteen antigens were identified as immunogenic in baboons, essential for persistence in the airway by TraDIS, and membrane-localized: BP0840 (OmpP), Pal, OmpA2, BP1485, BamA, Pcp, MlaA, YfgL, BP2197, BP1569, MlaD, ComL, and BP0183.

CONCLUSIONS

The B. pertussis antigens identified as immunogenic, essential for persistence in the airway, and membrane-localized warrant further investigation for inclusion in vaccines designed to reduce or prevent carriage of bacteria in the airway of vaccinated individuals.

A Minimally Morphologically Destructive Approach for DNA Retrieval and Whole-Genome Shotgun Sequencing of Pinned Historic Dipteran Vector Species

Museum collections contain enormous quantities of insect specimens collected over the past century, covering a period of increased and varied insecticide usage. These historic collections are therefore incredibly valuable as genomic snapshots of organisms before, during, and after exposure to novel selective pressures. However, these samples come with their own challenges compared with present-day collections, as they are fragile and retrievable DNA is low yield and fragmented. In this article, we tested several DNA extraction procedures across pinned historic Diptera specimens from four disease vector genera: Anopheles, Aedes, Culex, and Glossina. We identify an approach that minimizes morphological damage while maximizing DNA retrieval for Illumina library preparation and sequencing that can accommodate the fragmented and low yield nature of historic DNA. We identify several key points in retrieving sufficient DNA while keeping morphological damage to a minimum: an initial rehydration step, a short incubation without agitation in a modified low salt Proteinase K buffer (referred to as "lysis buffer C" throughout), and critical point drying of samples post-extraction to prevent tissue collapse caused by air drying. The suggested method presented here provides a solid foundation for exploring the genomes and morphology of historic Diptera collections.

Functional analysis of colonization factor antigen I positive enterotoxigenic Escherichia coli identifies genes implicated in survival in water and host colonization

Enterotoxigenic Escherichia coli (ETEC) expressing the colonization pili CFA/I are common causes of diarrhoeal infections in humans. Here, we use a combination of transposon mutagenesis and transcriptomic analysis to identify genes and pathways that contribute to ETEC persistence in water environments and colonization of a mammalian host. ETEC persisting in water exhibit a distinct RNA expression profile from those growing in richer media. Multiple pathways were identified that contribute to water survival, including lipopolysaccharide biosynthesis and stress response regulons. The analysis also indicated that ETEC growing in vivo in mice encounter a bottleneck driving down the diversity of colonizing ETEC populations.

Fundamental differences in physiology of Bordetella pertussis dependent on the two-component system Bvg revealed by gene essentiality studies

The identification of genes essential for a bacterium's growth reveals much about its basic physiology under different conditions. Bordetella pertussis, the causative agent of whooping cough, adopts both virulent and avirulent states through the activity of the two-component system, Bvg. The genes essential for B. pertussis growth in vitro were defined using transposon sequencing, for different Bvg-determined growth states. In addition, comparison of the insertion indices of each gene between Bvg phases identified those genes whose mutation exerted a significantly different fitness cost between phases. As expected, many of the genes identified as essential for growth in other bacteria were also essential for B. pertussis. However, the essentiality of some genes was dependent on Bvg. In particular, a number of key cell wall biosynthesis genes, including the entire mre/mrd locus, were essential for growth of the avirulent (Bvg minus) phase but not the virulent (Bvg plus) phase. In addition, cell wall biosynthesis was identified as a fundamental process that when disrupted produced greater fitness costs for the Bvg minus phase compared to the Bvg plus phase. Bvg minus phase growth was more susceptible than Bvg plus phase growth to the cell wall-disrupting antibiotic ampicillin, demonstrating the increased susceptibility of the Bvg minus phase to disruption of cell wall synthesis. This Bvg-dependent conditional essentiality was not due to Bvg-regulation of expression of cell wall biosynthesis genes; suggesting that this fundamental process differs between the Bvg phases in B. pertussis and is more susceptible to disruption in the Bvg minus phase. The ability of a bacterium to modify its cell wall synthesis is important when considering the action of antibiotics, particularly if developing novel drugs targeting cell wall synthesis.

Genome-wide Analysis of Salmonella enterica serovar Typhi in Humanized Mice Reveals Key Virulence Features

Salmonella enterica serovar Typhi causes typhoid fever only in humans. Murine infection with S. Typhimurium is used as a typhoid model, but its relevance to human typhoid is limited. Non-obese diabetic-scid IL2rγnull mice engrafted with human hematopoietic stem cells (hu-SRC-SCID) are susceptible to lethal S. Typhi infection. In this study, we use a high-density S. Typhi transposon library in hu-SRC-SCID mice to identify virulence loci using transposon-directed insertion site sequencing (TraDIS). Vi capsule, lipopolysaccharide (LPS), and aromatic amino acid biosynthesis were essential for virulence, along with the siderophore salmochelin. However, in contrast to the murine S. Typhimurium model, neither the PhoPQ two-component system nor the SPI-2 pathogenicity island was required for lethal S. Typhi infection, nor was the CdtB typhoid toxin. These observations highlight major differences in the pathogenesis of typhoid and non-typhoidal Salmonella infections and demonstrate the utility of humanized mice for understanding the pathogenesis of a human-specific pathogen.

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.

Retrospective application of transposon-directed insertion-site sequencing to investigate niche-specific virulence of Salmonella Typhimurium in cattle

Background

Salmonella enterica subspecies enterica is an animal and zoonotic pathogen of global importance. Cattle are a significant reservoir of human non-typhoidal salmonellosis and can suffer enteric and systemic disease owing to the ability of Salmonella to survive within the bovine lymphatic system and intestines. Contamination of food can occur due to the incorporation of contaminated peripheral lymph nodes or by direct contamination of carcasses with gut contents. It is essential to understand the mechanisms used by Salmonella to enter and persist within the bovine lymphatic system and how they differ from those required for intestinal colonization to minimize zoonotic infections.

Results

Transposon-directed insertion site sequencing (TraDIS) was applied to pools of mutants recovered from mesenteric lymph nodes (MLNs) draining the distal ileum of calves after oral inoculation with a library of 8550 random S. Typhimurium mini-Tn5Km2 mutants in pools of 475 mutants per calf. A total of 8315 mutants representing 2852 different genes were detected in MLNs and their in vivo fitness was calculated. Using the same improved algorithm for analysis of transposon-flanking sequences, the identity and phenotype of mutants recovered from the distal ileal mucosa of the same calves was also defined, enabling comparison with previously published data and of mutant phenotypes across the tissues. Phenotypes observed for the majority of mutants were highly significantly correlated in the two tissues. However, 32 genes were identified in which transposon insertions consistently resulted in differential fitness in the ileal wall and MLNs, suggesting niche-specific roles for these genes in pathogenesis. Defined null mutations affecting ptsN and spvC were confirmed to result in tissue-specific phenotypes in calves, thus validating the TraDIS dataset.

Conclusions

This validation of the role of thousands of Salmonella genes and identification of genes with niche-specific roles in a key target species will inform the design of control strategies for bovine salmonellosis and zoonotic infections, for which efficacious and cross-protective vaccines are currently lacking.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5319-0) contains supplementary material, which is available to authorized users.

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.

Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesis

Severe malaria is caused by the apicomplexan parasite Plasmodium falciparum. Despite decades of research, the distinct biology of these parasites has made it challenging to establish high-throughput genetic approaches to identify and prioritize therapeutic targets. Using transposon mutagenesis of P. falciparum in an approach that exploited its AT-rich genome, we generated more than 38,000 mutants, saturating the genome and defining mutability and fitness costs for over 87% of genes. Of 5399 genes, our study defined 2680 genes as essential for optimal growth of asexual blood stages in vitro. These essential genes are associated with drug resistance, represent leading vaccine candidates, and include approximately 1000 Plasmodium-conserved genes of unknown function. We validated this approach by testing proteasome pathways for individual mutants associated with artemisinin sensitivity.

Morphological, genomic and transcriptomic responses of Klebsiella pneumoniae to the last-line antibiotic colistin

Colistin remains one of the few antibiotics effective against multi-drug resistant (MDR) hospital pathogens, such as Klebsiella pneumoniae. Yet resistance to this last-line drug is rapidly increasing. Characterized mechanisms of col^R in K. pneumoniae are largely due to chromosomal mutations in two-component regulators, although a plasmid-mediated col^R mechanism has recently been uncovered. However, the effects of intrinsic colistin resistance are yet to be characterized on a whole-genome level. Here, we used a genomics-based approach to understand the mechanisms of adaptive col^R acquisition in K. pneumoniae. In controlled directed-evolution experiments we observed two distinct paths to colistin resistance acquisition. Whole genome sequencing identified mutations in two colistin resistance genes: in the known col^R regulator phoQ which became fixed in the population and resulted in a single amino acid change, and unstable minority variants in the recently described two-component sensor crrB. Through RNAseq and microscopy, we reveal the broad range of effects that colistin exposure has on the cell. This study is the first to use genomics to identify a population of minority variants with mutations in a col^R gene in K. pneumoniae.

Genome-wide transposon screening and quantitative insertion site sequencing for cancer gene discovery in mice

Transposon-mediated forward genetics screening in mice has emerged as a powerful tool for cancer gene discovery. It pinpoints cancer drivers that are difficult to find with other approaches, thus complementing the sequencing-based census of human cancer genes. We describe here a large series of mouse lines for insertional mutagenesis that are compatible with two transposon systems, PiggyBac and Sleeping Beauty, and give guidance on the use of different engineered transposon variants for constitutive or tissue-specific cancer gene discovery screening. We also describe a method for semiquantitative transposon insertion site sequencing (QiSeq). The QiSeq library preparation protocol exploits acoustic DNA fragmentation to reduce bias inherent to widely used restriction-digestion-based approaches for ligation-mediated insertion site amplification. Extensive multiplexing in combination with next-generation sequencing allows affordable ultra-deep transposon insertion site recovery in high-throughput formats within 1 week. Finally, we describe principles of data analysis and interpretation for obtaining insights into cancer gene function and genetic tumor evolution.

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.

A conditional piggyBac transposition system for genetic screening in mice identifies oncogenic networks in pancreatic cancer

Here we describe a conditional piggyBac transposition system in mice and report the discovery of large sets of new cancer genes through a pancreatic insertional mutagenesis screen. We identify Foxp1 as an oncogenic transcription factor that drives pancreatic cancer invasion and spread in a mouse model and correlates with lymph node metastasis in human patients with pancreatic cancer. The propensity of piggyBac for open chromatin also enabled genome-wide screening for cancer-relevant noncoding DNA, which pinpointed a Cdkn2a cis-regulatory region. Histologically, we observed different tumor subentities and discovered associated genetic events, including Fign insertions in hepatoid pancreatic cancer. Our studies demonstrate the power of genetic screening to discover cancer drivers that are difficult to identify by other approaches to cancer genome analysis, such as downstream targets of commonly mutated human cancer genes. These piggyBac resources are universally applicable in any tissue context and provide unique experimental access to the genetic complexity of cancer.

Generation of a Tn5 transposon library in Haemophilus parasuis and analysis by transposon-directed insertion-site sequencing (TraDIS)

Haemophilus parasuis is an important respiratory tract pathogen of swine and the etiological agent of Glässer's disease. The molecular pathogenesis of H. parasuis is not well studied, mainly due to the lack of efficient tools for genetic manipulation of this bacterium. In this study we describe a Tn5-based random mutagenesis method for use in H. parasuis. A novel chloramphenicol-resistant Tn5 transposome was electroporated into the virulent H. parasuis serovar 5 strain 29755. High transposition efficiency of Tn5, up to 10(4) transformants/μg of transposon DNA, was obtained by modification of the Tn5 DNA in the H. parasuis strain HS071 and establishment of optimal electrotransformation conditions, and a library of approximately 10,500 mutants was constructed. Analysis of the library using transposon-directed insertion-site sequencing (TraDIS) revealed that the insertion of Tn5 was evenly distributed throughout the genome. 10,001 individual mutants were identified, with 1561 genes being disrupted (69.4% of the genome). This newly-developed, efficient mutagenesis approach will be a powerful tool for genetic manipulation of H. parasuis in order to study its physiology and pathogenesis.

Evaluation and optimisation of preparative semi-automated electrophoresis systems for Illumina library preparation

Size selection can be a critical step in preparation of next-generation sequencing libraries. Traditional methods employing gel electrophoresis lack reproducibility, are labour intensive, do not scale well and employ hazardous interchelating dyes. In a high-throughput setting, solid-phase reversible immobilisation beads are commonly used for size-selection, but result in quite a broad fragment size range. We have evaluated and optimised the use of two semi-automated preparative DNA electrophoresis systems, the Caliper Labchip XT and the Sage Science Pippin Prep, for size selection of Illumina sequencing libraries.

Comparative genomics of the emerging human pathogen Photorhabdus asymbiotica with the insect pathogen Photorhabdus luminescens

Background

The Gram-negative bacterium Photorhabdus asymbiotica (Pa) has been recovered from human infections in both North America and Australia. Recently, Pa has been shown to have a nematode vector that can also infect insects, like its sister species the insect pathogen P. luminescens (Pl). To understand the relationship between pathogenicity to insects and humans in Photorhabdus we have sequenced the complete genome of Pa strain ATCC43949 from North America. This strain (formerly referred to as Xenorhabdus luminescens strain 2) was isolated in 1977 from the blood of an 80 year old female patient with endocarditis, in Maryland, USA. Here we compare the complete genome of Pa ATCC43949 with that of the previously sequenced insect pathogen P. luminescens strain TT01 which was isolated from its entomopathogenic nematode vector collected from soil in Trinidad and Tobago.

Results

We found that the human pathogen Pa had a smaller genome (5,064,808 bp) than that of the insect pathogen Pl (5,688,987 bp) but that each pathogen carries approximately one megabase of DNA that is unique to each strain. The reduced size of the Pa genome is associated with a smaller diversity in insecticidal genes such as those encoding the Toxin complexes (Tc's), Makes caterpillars floppy (Mcf) toxins and the Photorhabdus Virulence Cassettes (PVCs). The Pa genome, however, also shows the addition of a plasmid related to pMT1 from Yersinia pestis and several novel pathogenicity islands including a novel Type Three Secretion System (TTSS) encoding island. Together these data suggest that Pa may show virulence against man via the acquisition of the pMT1-like plasmid and specific effectors, such as SopB, that promote its persistence inside human macrophages. Interestingly the loss of insecticidal genes in Pa is not reflected by a loss of pathogenicity towards insects.

Conclusion

Our results suggest that North American isolates of Pa have acquired virulence against man via the acquisition of a plasmid and specific virulence factors with similarity to those shown to play roles in pathogenicity against humans in other bacteria.

Post-transcriptional control of nuclear-encoded cytochrome oxidase subunits in Trypanosoma brucei: evidence for genome-wide conservation of life-cycle stage-specific regulatory elements

Trypanosomes represent an excellent model for the post-transcriptional regulation of gene expression because their genome is organized into polycistronic transcription units. However, few signals governing developmental stage-specific expression have been identified, with there being no compelling evidence for widespread conservation of regulatory motifs. As a tool to search for common regulatory sequences we have used the nuclear-encoded components of the cytochrome oxidase (COX) complex of the trypanosome respiratory chain. Components of this complex represent a form of post-transcriptional operon because trypanosome mitochondrial activity is unusual in being developmentally programmed. By genome analysis we identified the genes for seven components of the COX complex. Each mRNA exhibits bloodstream stage-specific instability, which is not mediated by the RNA silencing pathway but which is alleviated by cycloheximide. Reporter assays have identified regulatory regions within the 3′-untranslated regions of three COX mRNAs operating principally at the translational level, but also via mRNA stability. Interrogation of the mapped regions via oligonucleotide frequency scoring provides evidence for genome-wide conservation of regulatory sequences among a large cohort of procyclic-enriched transcripts. Analysis of the co-regulated subunits of a stage-specific enzyme is therefore a novel approach to uncover cryptic regulatory sequences controlling gene expression at the post-transcriptional level.