Identification of essential genes in C. jejuni genome highlights hyper-variable plasticity regions

Functional genomics of Campylobacter -host interactions in an intestinal tissue model reveals a small lipoprotein essential for flagellar assembly

Campylobacter jejuni is currently the most common cause of bacterial gastroenteritis worldwide. However, its genome provides few clues about how it interacts with the host. Moreover, infection screens have often been limited to classical cell culture or animal models. To identify C. jejuni genes involved in host cell interactions, we applied transposon sequencing in a humanized 3D intestinal infection model based on tissue engineering. This revealed key proteins required for host cell adherence and/or internalization, including an Rrf2 family transcriptional regulator as well as three so far uncharacterized genes ( pflC / Cj1643 , pflD / Cj0892c , pflE / Cj0978c ), which we demonstrate to encode factors essential for motility. Deletion mutants of pflC / D / E are non-motile but retain intact, paralysed flagella filaments. We demonstrate that two of these newly identified motility proteins, PflC and PflD, are components of the C. jejuni 's periplasmic disk structures of the high torque motor. The third gene, pflE , encodes a small protein of only 57 aa. Using CryoET imaging we uncovered that the small protein has a striking effect on motor biogenesis, leading to a complete loss of the flagellar disk and motor structures upon its deletion. While PflE does not appear to be a structural component of the motor itself, our data suggests that it is a lipoprotein and supports localization of the main basal disk protein FlgP, which is the first assembly step of the flagellar disk structure. Despite being annotated as a lipoprotein, we find that C. jejuni FlgP instead relies on PflE for its association with the outer membrane. Overall, our genome-wide screen revealed novel C. jejuni host interaction factors including a transcriptional regulator as well as two structural components and a small protein crucial for biogenesis of the C. jejuni high torque flagella motor. Since the flagella machinery is a critical virulence determining factor for C. jejuni , our work demonstrates how such a small protein can, quite literally, bring a bacterial pathogen to a halt.

Monitoring AMR in Campylobacter jejuni from Italy in the last 10 years (2011-2021): Microbiological and WGS data risk assessment

Campylobacter jejuni is considered as the main pathogen in human food-borne outbreaks worldwide. Over the past years, several studies have reported antimicrobial resistance (AMR) in C. jejuni strains. In Europe, the official monitoring of AMR comprises the testing of Campylobacter spp. from food-producing animals because this microorganism is responsible for human infections and usually predominant in poultry. Food-producing animals are considered to be a major source of campylobacteriosis through contamination of food products. Concerns are growing due to the current classification of C. jejuni by the WHO as a 'high priority pathogen' due to the emergence of resistance to multiple drugs such as those belonging to the fluoroquinolones, macrolides and other classes, which limits the treatment alternatives. Knowledge about the contributions of different food sources to gastrointestinal disease is fundamental to prioritise food safety interventions and to establish proper control strategies. Assessing the genetic diversity among Campylobacter species is essential to the understanding of their epidemiology and population structure. Using a population genetic approach and grouping the isolates into sequence types within different clonal complexes, it is possible to investigate the source of the human cases. The work programme was aimed for the fellow to assess the AMR of C. jejuni isolated from humans, poultry and birds from wild and urban Italian habitats. Given the public health concern represented by resistant pathogens in food-producing animals and the paucity of data about this topic in Italy, the aim was to identify correlations between phenotypic and genotypic AMR and comparing the origin of the isolates. The work programme allowed the fellow to acquire knowledge, skills and competencies on the web-based tools used by IZSAM to process the NGS data and perform bioinformatics analyses for the identification of epidemiological clusters, the study of AMR patterns in C. jejuni isolates, and the assessment of the human exposure to such AMR pathogens. Furthermore, the fellow became able to transfer the acquired knowledge through innovative web-based didactical tools applied to WGS and clustering of specific food-borne pathogens, with particular reference to C. jejuni. To achieve this objective, 2,734 C. jejuni strains isolated from domestic and wild animals and humans, during the period 2011-2021 were analysed. The resistance phenotypes of the isolates were determined using the microdilution method with EUCAST breakpoints, for the following antibiotics: nalidixic acid, ciprofloxacin, chloramphenicol, erythromycin, gentamicin, streptomycin, tetracycline. The data were complemented by WGS data for each strain, uploaded in the Italian information system for the collection and analysis of complete genome sequence of pathogens isolated from animal, food and environment (GENPAT) developed and maintained at IZSAM; information like clonal complex and sequence type to understand the phylogenetical distance between strains according to their origins were also considered. This work underlines that a better knowledge of the resistance levels of C. jejuni is necessary, and mandatory monitoring of Campylobacter species in the different animal productions is strongly suggested.

Generating Transposon Insertion Libraries in Gram-Negative Bacteria for High-Throughput Sequencing

Transposon sequencing (Tn-seq) is a powerful method that combines transposon mutagenesis and massive parallel sequencing to identify genes and pathways that contribute to bacterial fitness under a wide range of environmental conditions. Tn-seq applications are extensive and have not only enabled examination of genotype-phenotype relationships at an organism level but also at the population, community and systems levels. Gram-negative bacteria are highly associated with antimicrobial resistance phenotypes, which has increased incidents of antibiotic treatment failure. Antimicrobial resistance is defined as bacterial growth in the presence of otherwise lethal antibiotics. The "last-line" antimicrobial colistin is used to treat Gram-negative bacterial infections. However, several Gram-negative pathogens, including Acinetobacter baumannii can develop colistin resistance through a range of molecular mechanisms, some of which were characterized using Tn-seq. Furthermore, signal transduction pathways that regulate colistin resistance vary within Gram-negative bacteria. Here we propose an efficient method of transposon mutagenesis in A. baumannii that streamlines generation of a saturating transposon insertion library and amplicon library construction by eliminating the need for restriction enzymes, adapter ligation, and gel purification. The methods described herein will enable in-depth analysis of molecular determinants that contribute to A. baumannii fitness when challenged with colistin. The protocol is also applicable to other Gram-negative ESKAPE pathogens, which are primarily associated with drug resistant hospital-acquired infections.

RiboFACSeq: A new method for investigating metabolic and transport pathways in bacterial cells by combining a riboswitch-based sensor, fluorescence-activated cell sorting and next-generation sequencing

The elucidation of the cellular processes involved in vitamin and cofactor biosynthesis is a challenging task. The conventional approaches to these investigations rely on the discovery and purification of the products (i.e proteins and metabolites) of a particular transport or biosynthetic pathway, prior to their subsequent analysis. However, the purification of low-abundance proteins or metabolites is a formidable undertaking that presents considerable technical challenges. As a solution, we present an alternative approach to such studies that circumvents the purification step. The proposed approach takes advantage of: (1) the molecular detection capabilities of a riboswitch-based sensor to detect the cellular levels of its cognate molecule, as a means to probe the integrity of the transport and biosynthetic pathways of the target molecule in cells, (2) the high-throughput screening ability of fluorescence-activated cell sorters to isolate cells in which only these specific pathways are disrupted, and (3) the ability of next-generation sequencing to quickly identify the genes of the FACS-sorted populations. This approach was named “RiboFACSeq”. Following their identification by RiboFACSeq, the role of these genes in the presumed pathway needs to be verified through appropriate functional assays. To demonstrate the utility of our approach, an adenosylcobalamin (AdoCbl)-responsive riboswitch-based sensor was used in this study to demonstrate that RiboFACSeq can be used to track and sort cells carrying genetic mutations in known AdoCbl transport and biosynthesis genes with desirable sensitivity and specificity. This method could potentially be used to elucidate any pathway of interest, as long as a suitable riboswitch-based sensor can be created. We believe that RiboFACSeq would be especially useful for the elucidation of biological pathways in which the proteins and/or their metabolites are present at very low physiological concentrations in cells, as is the case with vitamin and cofactor biosynthesis.

Double-Serine Fluoroquinolone Resistance Mutations Advance Major International Clones and Lineages of Various Multi-Drug Resistant Bacteria

The major international sequence types/lineages of methicillin-resistant Staphylococcus aureus (MRSA), extended-spectrum β-lactamase (ESBL)-producing Klebsiella pneumoniae and ESBL-producing E. coli were demonstrated to have been advanced by favorable fitness balance associated with high-level resistance to fluoroquinolones. The paper shows that favorable fitness in the major STs/lineages of these pathogens was principally attained by the capacity of evolving mutations in the fluoroquinolone-binding serine residues of both the DNA gyrase and topoisomerase IV enzymes. The available information on fitness balance incurred by individual and various combinations of mutations in the enzymes is reviewed in multiple species. Moreover, strong circumstantial evidence is presented that major STs/lineages of other multi-drug resistant bacteria, primarily vancomycin-resistant Enterococcus faecium (VRE), emerged by a similar mechanism. The reason(s) why the major ST/lineage strains of various pathogens proved more adept at evolving favorable mutations than most isolates of the same species remains to be elucidated.

Essential genome of Campylobacter jejuni

BACKGROUND

Campylobacter species are a leading cause of bacterial foodborne illness worldwide. Despite the global efforts to curb them, Campylobacter infections have increased continuously in both developed and developing countries. The development of effective strategies to control the infection by this pathogen is warranted. The essential genes of bacteria are the most prominent targets for this purpose. In this study, we used transposon sequencing (Tn-seq) of a genome-saturating library of Tn5 insertion mutants to define the essential genome of C. jejuni at a high resolution.

RESULT

We constructed a Tn5 mutant library of unprecedented complexity in C. jejuni NCTC 11168 with 95,929 unique insertions throughout the genome and used the genomic DNA of the library for the reconstruction of Tn5 libraries in the same (C. jejuni NCTC 11168) and different strain background (C. jejuni 81-176) through natural transformation. We identified 166 essential protein-coding genes and 20 essential transfer RNAs (tRNA) in C. jejuni NCTC 11168 which were intolerant to Tn5 insertions during in vitro growth. The reconstructed C. jejuni 81-176 library had 384 protein coding genes with no Tn5 insertions. Essential genes in both strain backgrounds were highly enriched in the cluster of orthologous group (COG) categories of 'Translation, ribosomal structure and biogenesis (J)', 'Energy production and conversion (C)', and 'Coenzyme transport and metabolism (H)'.

CONCLUSION

Comparative analysis among this and previous studies identified 50 core essential genes of C. jejuni, which can be further investigated for the development of novel strategies to control the spread of this notorious foodborne bacterial pathogen.

Profiling of Campylobacter jejuni Proteome in Exponential and Stationary Phase of Growth

Campylobacter jejuni has been reported as a major cause of bacterial food-borne enteritides in developed countries during the last decade. Despite its fastidious growth requirements, including low level of oxygen and high level of CO2, this pathogen is able to persist in the environment without permanent loss of its viability and virulence. As C. jejuni is not able to multiply outside a host, the cells spend significant amount of time in stationary phase of growth. The entry into the stationary phase is often correlated to resistance to various stresses in bacteria. The switching between exponential and stationary phases is frequently mediated by the regulator sigma S (RpoS). However, this factor is absent in C. jejuni and molecular mechanisms responsible for transition of cells to the stationary phase remain elusive. In this work, proteomic profiles of cells from exponential and stationary phases were compared using 2-D electrophoresis (2DE) fingerprinting combined with mass spectrometry analysis and qRT-PCR. The identified proteins, whose expression differed between the two phases, are mostly involved in protein biosynthesis, carbon metabolism, stress response and motility. Altered expression was observed also in the pleiotropic regulator CosR that was over-expressed during stationary phase. A shift between transcript and protein level evolution of CosR throughout the growth of C. jejuni was observed using qRT-PCR and (2DE). From these data, we hypothesized that CosR could undergo a negative autoregulation in stationary phase. A consensus sequence resulting from promoter sequence alignment of genes potentially regulated by CosR, including its own upstream region, among C. jejuni strains is proposed. To verify experimentally the potential autoregulation of CosR at the DNA level, electrophoretic mobility shift assay was performed with DNA fragments of CosR promoter region and rCosR. Different migration pattern of the promoter fragments indicates the binding capacity of CosR, suggesting its auto-regulation potential.

Genome-wide fitness analyses of the foodborne pathogen Campylobacter jejuni in in vitro and in vivo models

Campylobacter is the most common cause of foodborne bacterial illness worldwide. Faecal contamination of meat, especially chicken, during processing represents a key route of transmission to humans. There is a lack of insight into the mechanisms driving C. jejuni growth and survival within hosts and the environment. Here, we report a detailed analysis of C. jejuni fitness across models reflecting stages in its life cycle. Transposon (Tn) gene-inactivation libraries were generated in three C. jejuni strains and the impact on fitness during chicken colonisation, survival in houseflies and under nutrient-rich and -poor conditions at 4 °C and infection of human gut epithelial cells was assessed by Tn-insertion site sequencing (Tn-seq). A total of 331 homologous gene clusters were essential for fitness during in vitro growth in three C. jejuni strains, revealing that a large part of its genome is dedicated to growth. We report novel C. jejuni factors essential throughout its life cycle. Importantly, we identified genes that fulfil important roles across multiple conditions. Our comprehensive screens showed which flagella elements are essential for growth and which are vital to the interaction with host organisms. Future efforts should focus on how to exploit this knowledge to effectively control infections caused by C. jejuni.

Dissimilar Fitness Associated with Resistance to Fluoroquinolones Influences Clonal Dynamics of Various Multiresistant Bacteria

Fitness cost associated with resistance to fluoroquinolones was recently shown to vary across clones of methicillin-resistant Staphylococcus aureus and extended-spectrum β-lactamase-producing Klebsiella pneumoniae. The resulting dissimilar fitness should have influenced the clonal dynamics and thereby the rates of resistance for these pathogens. Moreover, a similar mechanism was recently proposed for the emergence of the H30 and H30R lineages of ESBL-producing E. coli and the major international clone (ribotype 027) of Clostridium difficile. Furthermore, several additional international clones of various multiresistant bacteria are suspect to have been selected by an analogous process. An ability to develop favorable mutations in the gyrase and topoisomerase IV genes seems to be a prerequisite for pathogens to retain fitness while showing high-level resistance to fluoroquinolones. Since, the consumption of other "non-fluoroquinolone" groups of antibiotics have also contributed to the rise in resistance rates a more judicious use of antibiotics in general and of fluoroquinolones in particular could ameliorate the international resistance situation.

Biofilm spatial organization by the emerging pathogen Campylobacter jejuni: comparison between NCTC 11168 and 81-176 strains under microaerobic and oxygen-enriched conditions

During the last years, Campylobacter has emerged as the leading cause of bacterial foodborne infections in developed countries. Described as an obligate microaerophile, Campylobacter has puzzled scientists by surviving a wide range of environmental oxidative stresses on foods farm to retail, and thereafter intestinal transit and oxidative damage from macrophages to cause human infection. In this study, confocal laser scanning microscopy (CLSM) was used to explore the biofilm development of two well-described Campylobacter jejuni strains (NCTC 11168 and 81-176) prior to or during cultivation under oxygen-enriched conditions. Quantitative and qualitative appraisal indicated that C. jejuni formed finger-like biofilm structures with an open ultrastructure for 81-176 and a multilayer-like structure for NCTC 11168 under microaerobic conditions (MAC). The presence of motile cells within the biofilm confirmed the maturation of the C. jejuni 81-176 biofilm. Acclimation of cells to oxygen-enriched conditions led to significant enhancement of biofilm formation during the early stages of the process. Exposure to these conditions during biofilm cultivation induced an even greater biofilm development for both strains, indicating that oxygen demand for biofilm formation is higher than for planktonic growth counterparts. Overexpression of cosR in the poorer biofilm-forming strain, NCTC 11168, enhanced biofilm development dramatically by promoting an open ultrastructure similar to that observed for 81-176. Consequently, the regulator CosR is likely to be a key protein in the maturation of C. jejuni biofilm, although it is not linked to oxygen stimulation. These unexpected data advocate challenging studies by reconsidering the paradigm of fastidious requirements for C. jejuni growth when various subpopulations (from quiescent to motile cells) coexist in biofilms. These findings constitute a clear example of a survival strategy used by this emerging human pathogen.

The Campylobacter jejuni CprRS two-component regulatory system regulates aspects of the cell envelope

Campylobacter jejuni is a leading cause of food-borne gastroenteritis in humans. It lives commensally in the gastrointestinal tract of animals, and tolerates variable conditions during transit/colonization of susceptible hosts. The C. jejuni CprRS two-component system contains an essential response regulator (CprR), and deletion of the cprS sensor kinase enhances biofilms. We sought to identify CprRS-regulated genes and better understand how the system affects survival. Expression from the cprR promoter was highest during logarithmic growth and dependent on CprS. CprR(D52A) did not support viability, indicating that CprR phosphorylation is essential despite the dispensability of CprS. We identified a GTAAAC consensus bound by the CprR C-terminus; the Asp52 residue of full-length CprR was required for binding, suggesting phosphorylation is required. Transcripts differing in expression in ΔcprS compared with wildtype (WT) contained a putative CprR binding site upstream of their promoter region and encoded htrA (periplasmic protease upstream of cprRS) and peb4 (SurA-like chaperone). Consistent with direct regulation, the CprR consensus in the htrA promoter was bound by CprR(CTD). Finally, ΔhtrA formed enhanced biofilms, and ΔcprS biofilms were suppressed by Mg(2+). CprRS is the first C. jejuni regulatory system shown to control genes related to the cell envelope, the first line of interaction between pathogen and changing environments.

Microarray transposon tracking for the mapping of conditionally essential genes in Campylobacter jejuni

Although whole genome approaches to the study of bacteria have grown substantially in the past decade, there is still a need for quick and easy methods for the determination of which genes are essential for the growth of these bacteria under specific growth conditions. There are numerous methods to accomplish this depending on the resources and equipment available, each with their own advantages and disadvantages. Here we describe a method we successfully employed to map the essential genes of Campylobacter jejuni using a microarray transposon tracking approach where we constructed a saturated transposon mutant library in the C. jejuni strain NCTC11168 and used a genomic microarray approach to identify genes lacking transposon insertions under standard laboratory growth conditions. With a fully saturated library, the absence of transposon insertions can be used as an indicator of a gene essential for the survival and growth for the conditions used for the mutant library.

A statistical framework for improving genomic annotations of transposon mutagenesis (TM) assigned essential genes

Whole-genome transposon mutagenesis (TM) experiment followed by sequence-based identification of insertion sites is the most popular genome-wise experiment to identify essential genes in Prokaryota. However, due to the limitation of high-throughput technique, this approach yields substantial systematic biases resulting in the incorrect assignments of many essential genes. To obtain unbiased and accurate annotations of essential genes from TM experiments, we developed a novel Poisson model based statistical framework to refine these TM assignments. In the model, first we identified and incorporated several potential factors such as gene length and TM insertion information which may cause the TM assignment biases into the basic Poisson model. Then we calculated the conditional probability of an essential gene given the observed TM insertion number. By factorizing this probability through introducing a latent variable the real insertion number, we formalized the statistical framework. Through iteratively updating and optimizing model parameters to maximize the goodness-of-fit of the model to the observed TM insertion data, we finalized the model. Using this model, we are able to assign the probability score of essentiality to each individual gene given its TM assignment, which subsequently correct the experimental biases. To enable our model widely useable, we established a user-friendly Web-server that is accessible to the public: http://research.cchmc.org/essentialgene/.

Defining the metabolic requirements for the growth and colonization capacity of Campylobacter jejuni

During the last decade Campylobacter jejuni has been recognized as the leading cause of bacterial gastroenteritis worldwide. This facultative intracellular pathogen is a member of the Epsilonproteobacteria and requires microaerobic atmosphere and nutrient rich media for efficient proliferation in vitro. Its catabolic capacity is highly restricted in contrast to Salmonella Typhimurium and other enteropathogenic bacteria because several common pathways for carbohydrate utilization are either missing or incomplete. Despite these metabolic limitations, C. jejuni efficiently colonizes various animal hosts as a commensal intestinal inhabitant. Moreover, C. jejuni is tremendously successful in competing with the human intestinal microbiota; an infectious dose of few hundreds bacteria is sufficient to overcome the colonization resistance of humans and can lead to campylobacteriosis. Besides the importance and clear clinical manifestation of this disease, the pathogenesis mechanisms of C. jejuni infections are still poorly understood. In recent years comparative genome sequence, transcriptome and metabolome analyses as well as mutagenesis studies combined with animal infection models have provided a new understanding of how the specific metabolic capacity of C. jejuni drives its persistence in the intestinal habitat of various hosts. Furthermore, new insights into the metabolic requirements that support the intracellular survival of C. jejuni were obtained. Because C. jejuni harbors distinct properties in establishing an infection in comparison to pathogenic Enterobacteriaceae, it represents an excellent organism for elucidating new aspects of the dynamic interaction and metabolic cross talk between a bacterial pathogen, the microbiota and the host.

Novel components of the flagellar system in epsilonproteobacteria

Motility is essential for the pathogenesis of many bacterial species. Most bacteria move using flagella, which are multiprotein filaments that rotate propelled by a cell wall-anchored motor using chemical energy. Although some components of the flagellar apparatus are common to many bacterial species, recent studies have shown significant differences in the flagellar structures of different bacterial species. The molecular bases for these differences, however, are not understood. The flagella from epsilonproteobacteria, which include the bacterial pathogens Campylobacter jejuni and Helicobacter pylori, are among the most divergent. Using next-generation sequencing combined with transposon mutagenesis, we have conducted a comprehensive high-throughput genetic screen in Campylobacter jejuni, which identified several novel components of its flagellar system. Biochemical analyses detected interactions between the identified proteins and known components of the flagellar machinery, and in vivo imaging located them to the bacterial poles, where flagella assemble. Most of the identified new components are conserved within but restricted to epsilonproteobacteria. These studies provide insight into the divergent flagella of this group of bacteria and highlight the complexity of this remarkable structure, which has adapted to carry out its conserved functions in the context of widely diverse bacterial species.

IMPORTANCE

Motility is essential for the normal physiology and pathogenesis of many bacterial species. Most bacteria move using flagella, which are multiprotein filaments that rotate propelled by a motor that uses chemical energy as fuel. Although some components of the flagellar apparatus are common to many bacterial species, recent studies have shown significant divergence in the flagellar structures across bacterial species. However, the molecular bases for these differences are not understood. The flagella from epsilonproteobacteria, which include the bacterial pathogens Campylobacter jejuni and Helicobacter pylori, are among the most divergent. We conducted a comprehensive genetic screen in Campylobacter jejuni and identified several novel components of the flagellar system. These studies provide important information to understand how flagella have adapted to function in the context of widely diverse sets of bacterial species and bring unique insight into the evolution and function of this remarkable bacterial organelle.

Phenotypic screening of a targeted mutant library reveals Campylobacter jejuni defenses against oxidative stress

During host colonization, Campylobacter jejuni is exposed to harmful reactive oxygen species (ROS) produced from the host immune system and from the gut microbiota. Consequently, identification and characterization of oxidative stress defenses are important for understanding how C. jejuni survives ROS stress during colonization of the gastrointestinal tract. Previous transcriptomic studies have defined the genes belonging to oxidant stimulons within C. jejuni. We have constructed isogenic deletion mutants of these identified genes to assess their role in oxidative stress survival. Phenotypic screening of 109 isogenic deletion mutants identified 22 genes which were either hypersensitive or hyposensitive to oxidants, demonstrating important roles for these genes in oxidant defense. The significance of these genes in host colonization was also assessed in an in vivo chick model of C. jejuni colonization. Overall, our findings identify an indirect role for motility in resistance to oxidative stress. We found that a nonmotile flagellum mutant, the ΔmotAB mutant, displayed increased sensitivity to oxidants. Restoration of sensitivity to superoxide in the ΔmotAB mutant was achieved by fumarate supplementation or tandem deletion of motAB with ccoQ, suggesting that disruption of the proton gradient across the inner membrane resulted in increased superoxide production in this strain. Furthermore, we have identified genes involved in cation transport and binding, detoxification, and energy metabolism that are also important factors in oxidant defense. This report describes the first isogenic deletion mutant library construction for screening of relevant oxidative stress defense genes within C. jejuni, thus providing a comprehensive analysis of the total set of oxidative stress defenses.

Training set selection for the prediction of essential genes

Various computational models have been developed to transfer annotations of gene essentiality between organisms. However, despite the increasing number of microorganisms with well-characterized sets of essential genes, selection of appropriate training sets for predicting the essential genes of poorly-studied or newly sequenced organisms remains challenging. In this study, a machine learning approach was applied reciprocally to predict the essential genes in 21 microorganisms. Results showed that training set selection greatly influenced predictive accuracy. We determined four criteria for training set selection: (1) essential genes in the selected training set should be reliable; (2) the growth conditions in which essential genes are defined should be consistent in training and prediction sets; (3) species used as training set should be closely related to the target organism; and (4) organisms used as training and prediction sets should exhibit similar phenotypes or lifestyles. We then analyzed the performance of an incomplete training set and an integrated training set with multiple organisms. We found that the size of the training set should be at least 10% of the total genes to yield accurate predictions. Additionally, the integrated training sets exhibited remarkable increase in stability and accuracy compared with single sets. Finally, we compared the performance of the integrated training sets with the four criteria and with random selection. The results revealed that a rational selection of training sets based on our criteria yields better performance than random selection. Thus, our results provide empirical guidance on training set selection for the identification of essential genes on a genome-wide scale.

Campylobacter jejuni pdxA affects flagellum-mediated motility to alter host colonization

Vitamin B6 (pyridoxal-5'-phosphate, PLP) is linked to a variety of biological functions in prokaryotes. Here, we report that the pdxA (putative 4-hydroxy-L-threonine phosphate dehydrogenase) gene plays a pivotal role in the PLP-dependent regulation of flagellar motility, thereby altering host colonization in a leading foodborne pathogen, Campylobacter jejuni. A C. jejuni pdxA mutant failed to produce PLP and exhibited a coincident loss of flagellar motility. Mass spectrometric analyses showed a 3-fold reduction in the main flagellar glycan pseudaminic acid (Pse) associated with the disruption of pdxA. The pdxA mutant also exhibited reduced growth rates compared with the WT strain. Comparative metabolomic analyses revealed differences in respiratory/energy metabolism between WT C. jejuni and the pdxA mutant, providing a possible explanation for the differential growth fitness between the two strains. Consistent with the lack of flagellar motility, the pdxA mutant showed impaired motility-mediated responses (bacterial adhesion, ERK1/2 activation, and IL-8 production) in INT407 cells and reduced colonization of chickens compared with the WT strain. Overall, this study demonstrated that the pdxA gene affects the PLP-mediated flagellar motility function, mainly through alteration of Pse modification, and the disruption of this gene also alters the respiratory/energy metabolisms to potentially affect host colonization. Our data therefore present novel implications regarding the utility of PLP and its dependent enzymes as potent target(s) for the control of this pathogen in the poultry host.

A statistical framework for improving genomic annotations of prokaryotic essential genes

Large-scale systematic analysis of gene essentiality is an important step closer toward unraveling the complex relationship between genotypes and phenotypes. Such analysis cannot be accomplished without unbiased and accurate annotations of essential genes. In current genomic databases, most of the essential gene annotations are derived from whole-genome transposon mutagenesis (TM), the most frequently used experimental approach for determining essential genes in microorganisms under defined conditions. However, there are substantial systematic biases associated with TM experiments. In this study, we developed a novel Poisson model–based statistical framework to simulate the TM insertion process and subsequently correct the experimental biases. We first quantitatively assessed the effects of major factors that potentially influence the accuracy of TM and subsequently incorporated relevant factors into the framework. Through iteratively optimizing parameters, we inferred the actual insertion events occurred and described each gene’s essentiality on probability measure. Evaluated by the definite mapping of essential gene profile in Escherichia coli, our model significantly improved the accuracy of original TM datasets, resulting in more accurate annotations of essential genes. Our method also showed encouraging results in improving subsaturation level TM datasets. To test our model’s broad applicability to other bacteria, we applied it to Pseudomonas aeruginosa PAO1 and Francisella tularensis novicida TM datasets. We validated our predictions by literature as well as allelic exchange experiments in PAO1. Our model was correct on six of the seven tested genes. Remarkably, among all three cases that our predictions contradicted the TM assignments, experimental validations supported our predictions. In summary, our method will be a promising tool in improving genomic annotations of essential genes and enabling large-scale explorations of gene essentiality. Our contribution is timely considering the rapidly increasing essential gene sets. A Webserver has been set up to provide convenient access to this tool. All results and source codes are available for download upon publication at http://research.cchmc.org/essentialgene/.

Contribution of amino acid catabolism to the tissue specific persistence of Campylobacter jejuni in a murine colonization model

Campylobacter jejuni is a major cause of food-borne disease in industrialized countries. Carbohydrate utilization by C. jejuni is severely restricted, and knowledge about which substrates fuel C. jejuni infection and growth is limited. Some amino acids have been shown to serve as carbon sources both in vitro and in vivo. In the present study we investigated the contribution of serine and proline catabolism to the invitro and invivo growth of C. jejuni 81-176. We confirmed that the serine transporter SdaC and the serine ammonia-lyase SdaA are required for serine utilization, and demonstrated that a predicted proline permease PutP and a bifunctional proline/delta-1-pyrroline-5-carboxylate dehydrogenase PutA are required for proline utilization by C. jejuni 81-176. C. jejuni 81-176 mutants unable to utilize serine were shown to be severely defective for colonization of the intestine and systemic tissues in a mouse model of infection. In contrast, C. jejuni 81-176 mutants unable to utilize proline were only defective for intestinal colonization. These results further emphasize the importance of amino acid utilization in C. jejuni colonization of various tissues.

Identification of essential genes of the periodontal pathogen Porphyromonas gingivalis

Background

Porphyromonas gingivalis is a Gram-negative anaerobic bacterium associated with periodontal disease onset and progression. Genetic tools for the manipulation of bacterial genomes allow for in-depth mechanistic studies of metabolism, physiology, interspecies and host-pathogen interactions. Analysis of the essential genes, protein-coding sequences necessary for survival of P. gingivalis by transposon mutagenesis has not previously been attempted due to the limitations of available transposon systems for the organism. We adapted a Mariner transposon system for mutagenesis of P. gingivalis and created an insertion mutant library. By analyzing the location of insertions using massively-parallel sequencing technology we used this mutant library to define genes essential for P. gingivalis survival under in vitro conditions.

Results

In mutagenesis experiments we identified 463 genes in P. gingivalis strain ATCC 33277 that are putatively essential for viability in vitro. Comparing the 463 P. gingivalis essential genes with previous essential gene studies, 364 of the 463 are homologues to essential genes in other species; 339 are shared with more than one other species. Twenty-five genes are known to be essential in P. gingivalis and B. thetaiotaomicron only. Significant enrichment of essential genes within Cluster of Orthologous Groups ‘D’ (cell division), ‘I’ (lipid transport and metabolism) and ‘J’ (translation/ribosome) were identified. Previously, the P. gingivalis core genome was shown to encode 1,476 proteins out of a possible 1,909; 434 of 463 essential genes are contained within the core genome. Thus, for the species P. gingivalis twenty-two, seventy-seven and twenty-three percent of the genome respectively are devoted to essential, core and accessory functions.

Conclusions

A Mariner transposon system can be adapted to create mutant libraries in P. gingivalis amenable to analysis by next-generation sequencing technologies. In silico analysis of genes essential for in vitro growth demonstrates that although the majority are homologous across bacterial species as a whole, species and strain-specific subsets are apparent. Understanding the putative essential genes of P. gingivalis will provide insights into metabolic pathways and niche adaptations as well as clinical therapeutic strategies.

Transcriptional regulation of the CmeABC multidrug efflux pump and the KatA catalase by CosR in Campylobacter jejuni

CosR is an essential response regulator in Campylobacter jejuni, a major food-borne pathogen causing enteritis worldwide. A transcriptomic analysis performed in this study discovered 93 genes whose transcriptional levels were changed >2-fold due to the repression of CosR expression by antisense peptide nucleic acid. The identified CosR-regulated genes are involved in various cellular functions, such as energy production, protein synthesis and folding, flagellum biogenesis, and lipid metabolism. Interestingly, 17 of the 93 CosR-regulated genes (18.3%) are predicted essential genes, indicating that CosR may participate in the regulation of vital biological processes in C. jejuni. In particular, CosR knockdown increased the transcriptional levels of cmeA, cmeB, and cmeC genes, whose protein product (CmeABC) is an important determinant conferring multidrug resistance in Campylobacter. Negative regulation of cmeABC by CosR was verified by quantitative real-time PCR (qRT-PCR) and P(cmeABC)::lacZ assay. The results of electrophoretic mobility shift assays (EMSAs) and DNase I footprinting assays demonstrated that CosR directly binds to the cmeABC promoter. Another notable finding is that CosR regulates the transcription of katA, the sole catalase gene in C. jejuni. Further characterization with qRT-PCR, the catalase enzyme assay, EMSA, and DNase I footprinting assays successfully demonstrated that CosR affects the katA transcription and the catalase activity by direct interactions with the katA promoter. The findings in this study clearly demonstrated that CosR regulates resistance mechanisms in C. jejuni by controlling the expression of genes involved in oxidative stress defense and extrusion of toxic compounds out of the cell.

Prediction and identification of sequences coding for orphan enzymes using genomic and metagenomic neighbours

Many characterized metabolic enzymes currently lack associated gene and protein sequences. Here, pathway and genomic neighbour data are used to assign genes to these ‘orphan enzymes,' and the predictions are validated with experimental assays and genome-scale metabolic modelling.

A computational method is developed for assigning candidate sequences to orphan enzymes. The method uses metabolic pathway, genomic neighbourhood, genomic co-occurrence, and protein domain information to predict genes that are likely to perform a particular enzymatic function.

Benchmarking of the scoring scheme based on the 4 features above revealed that some combinations of parameters yielded greater than 70% accuracy, and that high-confidence predictions could be generated for 131 orphan enzymes.

Enzyme assay experiments confirmed the predicted enzymatic activity for two of the high-confidence candidate sequences.

Predicted functions can improve the annotation of genomic and metagenomic data, and can reveal putative genes for enzymes with potential biotechnological applications.

Incorporating the predicted enzymatic reactions into genome-scale metabolic models changed the flux connectivity and improved their ability to correctly predict gene essentiality, supporting the biological relevance of these predictions.

Despite the current wealth of sequencing data, one-third of all biochemically characterized metabolic enzymes lack a corresponding gene or protein sequence, and as such can be considered orphan enzymes. They represent a major gap between our molecular and biochemical knowledge, and consequently are not amenable to modern systemic analyses. As 555 of these orphan enzymes have metabolic pathway neighbours, we developed a global framework that utilizes the pathway and (meta)genomic neighbour information to assign candidate sequences to orphan enzymes. For 131 orphan enzymes (37% of those for which (meta)genomic neighbours are available), we associate sequences to them using scoring parameters with an estimated accuracy of 70%, implying functional annotation of 16 345 gene sequences in numerous (meta)genomes. As a case in point, two of these candidate sequences were experimentally validated to encode the predicted activity. In addition, we augmented the currently available genome-scale metabolic models with these new sequence–function associations and were able to expand the models by on average 8%, with a considerable change in the flux connectivity patterns and improved essentiality prediction.

How a sugary bug gets through the day: recent developments in understanding fundamental processes impacting Campylobacter jejuni pathogenesis

Campylobacter jejuni is a highly prevalent yet fastidious bacterial pathogen that poses a significant health burden worldwide. Lacking many hallmark virulence factors, it is becoming increasingly clear that C. jejuni pathogenesis involves different strategies compared with other well-characterized enteric organisms. This includes the involvement of basic biological processes and cell envelope glycans in a number of aspects related to pathogenesis. The past few years have seen significant progress in the understanding of these pathways and how they relate to C. jejuni fundamental biology, stress survival, colonization, and virulence attributes. This review focuses on recent studies in three general areas where "pathogenesis" and "basic biology" overlap: physiology, stress responses and glycobiology.

In vivo and in silico determination of essential genes of Campylobacter jejuni

BACKGROUND

In the United Kingdom, the thermophilic Campylobacter species C. jejuni and C. coli are the most frequent causes of food-borne gastroenteritis in humans. While campylobacteriosis is usually a relatively mild infection, it has a significant public health and economic impact, and possible complications include reactive arthritis and the autoimmune diseases Guillain-Barré syndrome. The rapid developments in "omics" technologies have resulted in the availability of diverse datasets allowing predictions of metabolism and physiology of pathogenic micro-organisms. When combined, these datasets may allow for the identification of potential weaknesses that can be used for development of new antimicrobials to reduce or eliminate C. jejuni and C. coli from the food chain.

RESULTS

A metabolic model of C. jejuni was constructed using the annotation of the NCTC 11168 genome sequence, a published model of the related bacterium Helicobacter pylori, and extensive literature mining. Using this model, we have used in silico Flux Balance Analysis (FBA) to determine key metabolic routes that are essential for generating energy and biomass, thus creating a list of genes potentially essential for growth under laboratory conditions. To complement this in silico approach, candidate essential genes have been determined using a whole genome transposon mutagenesis method. FBA and transposon mutagenesis (both this study and a published study) predict a similar number of essential genes (around 200). The analysis of the intersection between the three approaches highlights the shikimate pathway where genes are predicted to be essential by one or more method, and tend to be network hubs, based on a previously published Campylobacter protein-protein interaction network, and could therefore be targets for novel antimicrobial therapy.

CONCLUSIONS

We have constructed the first curated metabolic model for the food-borne pathogen Campylobacter jejuni and have presented the resulting metabolic insights. We have shown that the combination of in silico and in vivo approaches could point to non-redundant, indispensable genes associated with the well characterised shikimate pathway, and also genes of unknown function specific to C. jejuni, which are all potential novel Campylobacter intervention targets.