Molecular signatures (unique proteins and conserved indels) that are specific for the epsilon proteobacteria (Campylobacterales)

Application of TraDIS to define the core essential genome of Campylobacter jejuni and Campylobacter coli

Campylobacter species are the major cause of bacterial gastroenteritis. As there is no effective vaccine, combined with the rapid increase in antimicrobial resistant strains, there is a need to identify new targets for intervention. Essential genes are those that are necessary for growth and/or survival, making these attractive targets. In this study, comprehensive transposon mutant libraries were created in six C. jejuni strains, four C. coli strains and one C. lari and C. hyointestinalis strain, allowing for those genes that cannot tolerate a transposon insertion being called as essential. Comparison of essential gene lists using core genome analysis can highlight those genes which are common across multiple strains and/or species. Comparison of C. jejuni and C. coli, the two species that cause the most disease, identified 316 essential genes. Genes of interest highlighted members of the purine pathway being essential for C. jejuni whilst also finding that a functional potassium uptake system is essential. Protein-protein interaction networks using these essential gene lists also highlighted proteins in the purine pathway being major 'hub' proteins which have a large number of interactors across the network. When adding in two more species (C. lari and C. hyointestinalis) the essential gene list reduces to 261. Within these 261 essential genes, there are many genes that have been found to be essential in other bacteria. These include htrB and PEB4, which have previously been found as core virulence genes across Campylobacter species in other studies. There were 21 genes which have no known function with eight of these being associated with the membrane. These surface-associated essential genes may provide attractive targets. The essential gene lists presented will help to prioritise targets for the development of novel therapeutic and preventative interventions.

Phylogenomic and comparative genomic studies robustly demarcate two distinct clades of Pseudomonas aeruginosa strains: proposal to transfer the strains from an outlier clade to a novel species Pseudomonas paraeruginosa sp. nov

The strains of Pseudomonas aeruginosa exhibit considerable differences in their genotypic and pathogenic properties. To clarify their evolutionary/taxonomic relationships, comprehensive phylogenomic and comparative genomic studies were conducted on the genome sequences of 212 P . aeruginosa strains covering their genetic diversity. In a phylogenomic tree based on 118 conserved proteins, the analysed strains formed two distinct clades. One of these clades, Clade-1, encompassing >70 % of the strains including the type strain DSM 50071T, represents the species P. aeruginosa sensu stricto. Clade-2, referred to in earlier work as the outlier group, with NCTC 13628T as its type strain, constitutes a novel species level lineage. The average nucleotide identity, average amino acid identity and digital DNA–DNA hybridization values between the strains from Clade-1 and Clade-2 are in the range of 93.4–93.7, 95.1–95.3 and 52–53 %, respectively. The 16S rRNA gene of P. aeruginosa DSM 50071T also shows 98.3 % similarity to that of NCTC 13628T. These values are lower than the suggested cut-off values for species distinction, indicating that the Clade-2 strains (NCTC 13628T) constitute a new species. We also report the identification of 12 conserved signature indels in different proteins and 24 conserved signature proteins that are exclusively found in either Clade-1 or Clade-2, providing a reliable means for distinguishing these clades. Additionally, in contrast to swimming motility, twitching motility is only present in Clade-1 strains. Based on earlier work, the strains from these two clades also differ in their pathogenic mechanisms (presence/absence of Type III secretion system), production of biosurfactants, phenazines and siderophores, and several other genomic characteristics. Based on the evidence from different studies, we propose that the Clade-2 strains constitute a novel species for which the name Pseudomonas paraeruginosa is proposed. The type strain is NCTC 13628T (=PA7T=ATCC 9027T). The description of Pseudomonas aeruginosa is also emended to include information for different molecular markers specific for this species.

The Bacterial Symbionts of Closely Related Hydrothermal Vent Snails With Distinct Geochemical Habitats Show Broad Similarity in Chemoautotrophic Gene Content

Symbiosis has evolved between a diversity of invertebrate taxa and chemosynthetic bacterial lineages. At the broadest level, these symbioses share primary function: the bacterial symbionts use the energy harnessed from the oxidation of reduced chemicals to power the fixation of inorganic carbon and/or other nutrients, providing the bulk of host nutrition. However, it is unclear to what extent the ecological niche of the host species is influenced by differences in symbiont traits, particularly those involved in chemoautotrophic function and interaction with the geochemical environment. Hydrothermal vents in the Lau Basin (Tonga) are home to four morphologically and physiologically similar snail species from the sister genera Alviniconcha and Ifremeria. Here, we assembled nearly complete genomes from their symbionts to determine whether differences in chemoautotrophic capacity exist among these symbionts that could explain the observed distribution of these snail species into distinct geochemical habitats. Phylogenomic analyses confirmed that the symbionts have evolved from four distinct lineages in the classes γ-proteobacteria or Campylobacteria. The genomes differed with respect to genes related to motility, adhesion, secretion, and amino acid uptake or excretion, though were quite similar in chemoautotrophic function, with all four containing genes for carbon fixation, sulfur and hydrogen oxidation, and oxygen and nitrate respiration. This indicates that differences in the presence or absence of symbiont chemoautotrophic functions does not likely explain the observed geochemical habitat partitioning. Rather, differences in gene expression and regulation, biochemical differences among these chemoautotrophic pathways, and/or differences in host physiology could all influence the observed patterns of habitat partitioning.

Diversity of transducer-like proteins (Tlps) in Campylobacter

Campylobacter transducer-like proteins (Tlps), also known as methyl-accepting chemotaxis proteins (MCPs), are associated with virulence as well as niche and host adaptation. While functional attributes of these proteins are being elucidated, little has been published regarding their sequence diversity or chromosomal locations and context, although they appear to define invertible regions within Campylobacter jejuni genomes. Genome assemblies for several species of Campylobacter were obtained from the publicly available NCBI repositories. Genomes from all isolates were obtained from GenBank and assessed for Tlp content, while data from isolates with complete, finished genomes were used to determine the identity of Tlps as well as the gene content of putative invertible elements (IEs) in C. jejuni (Cj) and C. coli (Cc). Tlps from several Campylobacter species were organized into a nomenclature system and novel Tlps were defined and named for Cj and Cc. The content of Tlps appears to be species-specific, though diverse within species. Cj and Cc carried overlapping, related Tlp content, as did the three C. fetus subspecies. Tlp1 was detected in 88% of Cj isolates and approximately 43% of Cc, and was found in a different conserved chromosomal location and genetic context in each species. Tlp1 and Tlp 3 predominated in genomes from Cj whereas other Tlps were detected less frequently. Tlp13 and Tlp20 predominated in genomes from Cc while some Cj/Cc Tlps were not detected at all. Tlps 2–4 and 11–20 were less frequently detected and many showed sequence heterogeneity that could affect substrate binding, signal transduction, or both. Tlps other than Tlp1, 7, and 10 had substantial sequence identity in the C-terminal half of the protein, creating chromosomal repeats potentially capable of mediating the inversion of large chromosomal DNA. Cj and Cc Tlps were both found in association with only 14 different genes, indicating a limited genomic context. In Cj these Tlps defined IEs that were for the most part found at a single chromosomal location and comprised of a conserved set of genes. Cc IEs were situated at very different chromosomal locations, had different structures than Cj IEs, and were occasionally incomplete, therefore not capable of inversion. Tlps may have a role in Campylobacter genome structure and dynamics as well as acting as chemoreceptors mediating chemotactic responses.

Function of Serine Protease HtrA in the Lifecycle of the Foodborne Pathogen Campylobacter jejuni

Campylobacter jejuni is a major food-borne zoonotic pathogen, responsible for a large proportion of bacterial gastroenteritis cases, as well as Guillian-Barré and Miller-Fisher syndromes. During infection, tissue damage is mainly caused by bacteria invading epithelial cells and traversing the intestinal barrier. C. jejuni is able to enter the lamina propria and the bloodstream and may move into other organs, such as spleen, liver, or mesenteric lymph nodes. However, the involved molecular mechanisms are not fully understood. C. jejuni can transmigrate effectively across polarized intestinal epithelial cells mainly by the paracellular route using the serine protease high-temperature requirement A (HtrA). However, it appears that HtrA has a dual function, as it also acts as a chaperone, interacting with denatured or misfolded periplasmic proteins under stress conditions. Here, we review recent progress on the role of HtrA in C. jejuni pathogenesis. HtrA can be transported into the extracellular space and cleaves cell-to-cell junction factors, such as E-cadherin and probably others, disrupting the epithelial barrier and enabling paracellular transmigration of the bacteria. The secretion of HtrA is a newly discovered strategy also utilized by other pathogens. Thus, secreted HtrA proteases represent highly attractive targets for anti-bacterial treatment and may provide a suitable candidate for vaccine development.

Structure and Function of the Campylobacter jejuni Chromosome Replication Origin

Campylobacter jejuni is the leading bacterial cause of foodborne infections worldwide. However, our understanding of its cell cycle is poor. We identified the probable C. jejuni origin of replication (oriC) - a key element for initiation of chromosome replication, which is also important for chromosome structure, maintenance and dynamics. The herein characterized C. jejuni oriC is monopartite and contains (i) the DnaA box cluster, (ii) the DnaA-dependent DNA unwinding element (DUE) and (iii) binding sites for regulatory proteins. The cluster of five DnaA boxes and the DUE were found in the dnaA-dnaN intergenic region. Binding of DnaA to this cluster of DnaA-boxes enabled unwinding of the DUE in vitro. However, it was not sufficient to sustain replication of minichromosomes, unless the cluster was extended by additional DnaA boxes located in the 3' end of dnaA. This suggests, that C. jejuni oriC requires these boxes to initiate or to regulate replication of its chromosome. However, further detailed mutagenesis is required to confirm the role of these two boxes in initiation of C. jejuni chromosome replication and thus to confirm partial localization of C. jejuni oriC within a coding region, which has not been reported thus far for any bacterial oriC. In vitro DUE unwinding by DnaA was inhibited by Cj1509, an orphan response regulator and a homolog of HP1021, that has been previously shown to inhibit replication in Helicobacter pylori. Thus, Cj1509 might play a similar role as a regulator of C. jejuni chromosome replication. This is the first systematic analysis of chromosome replication initiation in C. jejuni, and we expect that these studies will provide a basis for future research examining the structure and dynamics of the C. jejuni chromosome, which will be crucial for understanding the pathogens' life cycle and virulence.

Acid-regulated gene expression of Helicobacter pylori: Insight into acid protection and gastric colonization

BACKGROUND

The pathogen Helicobacter pylori encounters many stressors as it transits to and infects the gastric epithelium. Gastric acidity is the predominate stressor encountered by the bacterium during initial infection and establishment of persistent infection. H. pylori initiates a rapid response to acid to maintain intracellular pH and proton motive force appropriate for a neutralophile. However, acid sensing by H. pylori may also serve as a transcriptional trigger to increase the levels of other pathogenic factors needed to subvert host defenses such as acid acclimation, antioxidants, flagellar synthesis and assembly, and CagA secretion.

MATERIALS AND METHODS

Helicobacter pylori were acid challenged at pH 3.0, 4.5, 6.0 vs nonacidic pH for 4 hours in the presence of urea, followed by RNA-seq analysis and qPCR. Cytoplasmic pH was monitored under the same conditions.

RESULTS

About 250 genes were induced, and an equal number were repressed at acidic pHs. Genes encoding for antioxidant proteins, flagellar structural proteins, particularly class 2 genes, T4SS/Cag-PAI, F_o F₁ -ATPase, and proteins involved in acid acclimation were highly expressed at acidic pH. Cytoplasmic pH decreased from 7.8 at pH_out of 8.0 to 6.0 at pH_out of 3.0.

CONCLUSIONS

These results suggest that increasing extracellular or intracellular acidity or both are detected by the bacterium and serve as a signal to initiate increased production of protective and pathogenic factors needed to counter host defenses for persistent infection. These changes are dependent on degree of acidity and time of acid exposure, triggering a coordinated response to the environment required for colonization.

Changes in the gut microbiome of the Chinese mitten crab (Eriocheir sinensis) in response to White spot syndrome virus (WSSV) infection

Intestinal microorganisms play important roles in maintaining host health, but their functions in aquatic animal hosts have yet to be fully elucidated. The Chinese mitten crab, Eriocheir sinensis, is one such example. We attempted to identify the shift of gut microbiota that occurred in response to infection of white spot syndrome virus (WSSV), an emerging viral pathogen in the crab aquaculture industry. The microbiota may exert some control over aspects of the viral pathogenesis. We investigated the changes in composition and structure of the crab gut microbiome during various WSSV infection stages of 6 h post-infection (hpi) and 48 hpi, using a 16S rRNA approach on the MiSeq Illumina sequencing platform. Four phyla (Firmicutes, Proteobacteria, Tenericutes and Bacteroidetes) were most dominant in the gut of E. sinensis regardless of the WSSV infection stages. However, further analysis revealed that over 12 bacterial phyla, 44 orders and 68 families were significantly different in abundance at various states of WSSV infection. Several intriguing aspects of E. sinensis gut bacteria that had not been previously reported were also uncovered, such as class Mollicutes was dominant here, but absent in crabs from Yangtze River estuary and Chongming Islands. Overall, this study provided the first evidence that changes in gut microbiome were closely associated with the severity of WSSV infection and that indicator taxa could be used to evaluate the crab health status.

Unique and Universal Features of Epsilonproteobacterial Origins of Chromosome Replication and DnaA-DnaA Box Interactions

In bacteria, chromosome replication is initiated by the interaction of the initiator protein DnaA with a defined region of a chromosome at which DNA replication starts (oriC). While DnaA proteins share significant homology regardless of phylogeny, oriC regions exhibit more variable structures. The general architecture of oriCs is universal, i.e., they are composed of a cluster of DnaA binding sites, a DNA-unwinding element, and sequences that bind regulatory proteins. However, detailed structures of oriCs are shared by related species while being significantly different in unrelated bacteria. In this work, we characterized Epsilonproteobacterial oriC regions. Helicobacter pylori was the only species of the class for which oriC was characterized. A few unique features were found such as bipartite oriC structure, not encountered in any other Gram-negative species, and topology-sensitive DnaA-DNA interactions, which have not been found in any other bacterium. These unusual H. pylori oriC features raised questions of whether oriC structure and DnaA-DNA interactions are unique to this bacterium or whether they are common to related species. By in silico and in vitro analyses we identified putative oriCs in three Epsilonproteobacterial species: pathogenic Arcobacter butzleri, symbiotic Wolinella succinogenes, and free-living Sulfurimonas denitrificans. We propose that oriCs typically co-localize with ruvC-dnaA-dnaN in Epsilonproteobacteria, with the exception of Helicobacteriaceae species. The clusters of DnaA boxes localize upstream (oriC1) and downstream (oriC2) of dnaA, and they likely constitute bipartite origins. In all cases, DNA unwinding was shown to occur in oriC2. Unlike the DnaA box pattern, which is not conserved in Epsilonproteobacterial oriCs, the consensus DnaA box sequences and the mode of DnaA-DnaA box interactions are common to the class. We propose that the typical Epsilonproteobacterial DnaA box consists of the core nucleotide sequence 5′-TTCAC-3′ (4–8 nt), which, together with the significant changes in the DNA-binding motif of corresponding DnaAs, determines the unique molecular mechanism of DnaA-DNA interaction. Our results will facilitate identification of oriCs and subsequent identification of factors which regulate chromosome replication in other Epsilonproteobacteria. Since replication is controlled at the initiation step, it will help to better characterize life cycles of these species, many of which are considered as emerging pathogens.

Impact of genomics on the understanding of microbial evolution and classification: the importance of Darwin's views on classification

Analyses of genome sequences, by some approaches, suggest that the widespread occurrence of horizontal gene transfers (HGTs) in prokaryotes disguises their evolutionary relationships and have led to questioning of the Darwinian model of evolution for prokaryotes. These inferences are critically examined in the light of comparative genome analysis, characteristic synapomorphies, phylogenetic trees and Darwin's views on examining evolutionary relationships. Genome sequences are enabling discovery of numerous molecular markers (synapomorphies) such as conserved signature indels (CSIs) and conserved signature proteins (CSPs), which are distinctive characteristics of different prokaryotic taxa. Based on these molecular markers, exhibiting high degree of specificity and predictive ability, numerous prokaryotic taxa of different ranks, currently identified based on the 16S rRNA gene trees, can now be reliably demarcated in molecular terms. Within all studied groups, multiple CSIs and CSPs have been identified for successive nested clades providing reliable information regarding their hierarchical relationships and these inferences are not affected by HGTs. These results strongly support Darwin's views on evolution and classification and supplement the current phylogenetic framework based on 16S rRNA in important respects. The identified molecular markers provide important means for developing novel diagnostics, therapeutics and for functional studies providing important insights regarding prokaryotic taxa.

Cloning, Expression, Purification, Regulation, and Subcellular Localization of a Mini-protein from Campylobacter jejuni

The Cj1169c-encoded putative protein of Campylobacter was expressed and purified from E. coli after sequence optimization. The purified protein allowed the production of a specific rabbit antiserum that was used to study the protein expression in vitro and its subcellular localization in the bacterial cell and putative interactions with other proteins. This protein is produced in Campylobacter and it clearly localizes into the periplasmic space. The level of protein production depends on factors, including pH, temperature, osmolarity, and growth phase suggesting a role in the Campylobacter environmental adaptation. The cysteine residues present in the sequence are probably involved in disulfide bridges, which may promote covalent interactions with other proteins of the Campylobacter envelope.

The Biological Fight Against Pathogenic Bacteria and Protozoa

The animal gastrointestinal tract is a tube with two open ends; hence, from the microbial point of view it constitutes an open system, as opposed to the circulatory system that must be a tightly closed microbial-free environment. In particular, the human intestine spans ca. 200 m² and represents a massive absorptive surface composed of a layer of epithelial cells as well as a paracellular barrier. The permeability of this paracellular barrier is regulated by transmembrane proteins known as claudins that play a critical role in tight junctions.

An abundance of Epsilonproteobacteria revealed in the gut microbiome of the laboratory cultured sea urchin, Lytechinus variegatus

In this study, we have examined the bacterial community composition of the laboratory cultured sea urchin Lytechinus variegatus gut microbiome and its culture environment using NextGen amplicon sequencing of the V4 segment of the 16S rRNA gene, and downstream bioinformatics tools. Overall, the gut and tank water was dominated by Proteobacteria, whereas the feed consisted of a co-occurrence of Proteobacteria and Firmicutes at a high abundance. The gut tissue represented Epsilonproteobacteria as dominant, with order Campylobacterales at the highest relative abundance (>95%). However, the pharynx tissue was dominated by class Alphaproteobacteria. The gut digesta and egested fecal pellets had a high abundance of class Gammaproteobacteria, from which Vibrio was found to be the primary genus, and Epsilonproteobacteria, with genus Arcobacter occurring at a moderate level. At the class level, the tank water was dominated by Gammaproteobacteria, and the feed by Alphaproteobacteria. Multi-Dimensional Scaling analysis showed that the microbial community of the gut tissue clustered together, as did the pharynx tissue to the feed. The gut digesta and egested fecal pellets showed a similarity relationship to the tank water. Further analysis of Campylobacterales at a lower taxonomic level using the oligotyping method revealed 37 unique types across the 10 samples, where Oligotype 1 was primarily represented in the gut tissue. BLAST analysis identified Oligotype 1 to be Arcobacter sp., Sulfuricurvum sp., and Arcobacter bivalviorum at an identity level >90%. This study showed that although distinct microbial communities are evident across multiple components of the sea urchin gut ecosystem, there is a noticeable correlation between the overall microbial communities of the gut with the sea urchin L. variegatus culture environment.

Bacterial community shift in the coastal Gulf of Mexico salt-marsh sediment microcosm in vitro following exposure to the Mississippi Canyon Block 252 oil (MC252)

In this study, we examined the responses by the indigenous bacterial communities in salt-marsh sediment microcosms in vitro following treatment with Mississippi Canyon Block 252 oil (MC252). Microcosms were constructed of sediment and seawater collected from Bayou La Batre located in coastal Alabama on the Gulf of Mexico. We used an amplicon pyrosequencing approach on microcosm sediment metagenome targeting the V3–V5 region of the 16S rRNA gene. Overall, we identified a shift in the bacterial community in three distinct groups. The first group was the early responders (orders Pseudomonadales and Oceanospirillales within class Gammaproteobacteria), which increased their relative abundance within 2 weeks and were maintained 3 weeks after oil treatment. The second group was identified as early, but transient responders (order Rhodobacterales within class Alphaproteobacteria; class Epsilonproteobacteria), which increased their population by 2 weeks, but returned to the basal level 3 weeks after oil treatment. The third group was the late responders (order Clostridiales within phylum Firmicutes; order Methylococcales within class Gammaproteobacteria; and phylum Tenericutes), which only increased 3 weeks after oil treatment. Furthermore, we identified oil-sensitive bacterial taxa (order Chromatiales within class Gammaproteobacteria; order Syntrophobacterales within class Deltaproteobacteria), which decreased in their population after 2 weeks of oil treatment. Detection of alkane (alkB), catechol (C2,3DO) and biphenyl (bph) biodegradation genes by PCR, particularly in oil-treated sediment metacommunity DNA, delineates proliferation of  the hydrocarbon degrading bacterial community. Overall, the indigenous bacterial communities in our salt-marsh sediment in vitro microcosm study responded rapidly and shifted towards members of the taxonomic groups that are capable of surviving in an MC252 oil-contaminated environment.

Diversity of the Epsilonproteobacteria Dsb (disulfide bond) systems

The bacterial proteins of the Dsb family-important components of the post-translational protein modification system-catalyze the formation of disulfide bridges, a process that is crucial for protein structure stabilization and activity. Dsb systems play an essential role in the assembly of many virulence factors. Recent rapid advances in global analysis of bacteria have thrown light on the enormous diversity among bacterial Dsb systems. While the Escherichia coli disulfide bond-forming system is quite well understood, the mechanisms of action of Dsb systems in other bacteria, including members of class Epsilonproteobacteria that contain pathogenic and non-pathogenic bacteria colonizing extremely diverse ecological niches, are poorly characterized. Here we present a review of current knowledge on Epsilonproteobacteria Dsb systems. We have focused on the Dsb systems of Campylobacter spp. and Helicobacter spp. because our knowledge about Dsb proteins of Wolinella and Arcobacter spp. is still scarce and comes mainly from bioinformatic studies. Helicobacter pylori is a common human pathogen that colonizes the gastric epithelium of humans with severe consequences. Campylobacter spp. is a leading cause of zoonotic enteric bacterial infections in most developed and developing nations. We focus on various aspects of the diversity of the Dsb systems and their influence on pathogenicity, particularly because Dsb proteins are considered as potential targets for a new class of anti-virulence drugs to treat human infections by Campylobacter or Helicobacter spp.

Horizontal transfers and gene losses in the phospholipid pathway of bartonella reveal clues about early ecological niches

Bartonellae are mammalian pathogens vectored by blood-feeding arthropods. Although of increasing medical importance, little is known about their ecological past, and host associations are underexplored. Previous studies suggest an influence of horizontal gene transfers in ecological niche colonization by acquisition of host pathogenicity genes. We here expand these analyses to metabolic pathways of 28 Bartonella genomes, and experimentally explore the distribution of bartonellae in 21 species of blood-feeding arthropods. Across genomes, repeated gene losses and horizontal gains in the phospholipid pathway were found. The evolutionary timing of these patterns suggests functional consequences likely leading to an early intracellular lifestyle for stem bartonellae. Comparative phylogenomic analyses discover three independent lineage-specific reacquisitions of a core metabolic gene—NAD(P)H-dependent glycerol-3-phosphate dehydrogenase (gpsA)—from Gammaproteobacteria and Epsilonproteobacteria. Transferred genes are significantly closely related to invertebrate Arsenophonus-, and Serratia-like endosymbionts, and mammalian Helicobacter-like pathogens, supporting a cellular association with arthropods and mammals at the base of extant Bartonella spp. Our studies suggest that the horizontal reacquisitions had a key impact on bartonellae lineage specific ecological and functional evolution.

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.

Protein based molecular markers provide reliable means to understand prokaryotic phylogeny and support Darwinian mode of evolution

The analyses of genome sequences have led to the proposal that lateral gene transfers (LGTs) among prokaryotes are so widespread that they disguise the interrelationships among these organisms. This has led to questioning of whether the Darwinian model of evolution is applicable to prokaryotic organisms. In this review, we discuss the usefulness of taxon-specific molecular markers such as conserved signature indels (CSIs) and conserved signature proteins (CSPs) for understanding the evolutionary relationships among prokaryotes and to assess the influence of LGTs on prokaryotic evolution. The analyses of genomic sequences have identified large numbers of CSIs and CSPs that are unique properties of different groups of prokaryotes ranging from phylum to genus levels. The species distribution patterns of these molecular signatures strongly support a tree-like vertical inheritance of the genes containing these molecular signatures that is consistent with phylogenetic trees. Recent detailed studies in this regard on the Thermotogae and Archaea, which are reviewed here, have identified large numbers of CSIs and CSPs that are specific for the species from these two taxa and a number of their major clades. The genetic changes responsible for these CSIs (and CSPs) initially likely occurred in the common ancestors of these taxa and then vertically transferred to various descendants. Although some CSIs and CSPs in unrelated groups of prokaryotes were identified, their small numbers and random occurrence has no apparent influence on the consistent tree-like branching pattern emerging from other markers. These results provide evidence that although LGT is an important evolutionary force, it does not mask the tree-like branching pattern of prokaryotes or understanding of their evolutionary relationships. The identified CSIs and CSPs also provide novel and highly specific means for identification of different groups of microbes and for taxonomical and biochemical studies.

Directed evolution and in silico analysis of reaction centre proteins reveal molecular signatures of photosynthesis adaptation to radiation pressure

Evolutionary mechanisms adopted by the photosynthetic apparatus to modifications in the Earth's atmosphere on a geological time-scale remain a focus of intense research. The photosynthetic machinery has had to cope with continuously changing environmental conditions and particularly with the complex ionizing radiation emitted by solar flares. The photosynthetic D1 protein, being the site of electron tunneling-mediated charge separation and solar energy transduction, is a hot spot for the generation of radiation-induced radical injuries. We explored the possibility to produce D1 variants tolerant to ionizing radiation in Chlamydomonas reinhardtii and clarified the effect of radiation-induced oxidative damage on the photosynthetic proteins evolution. In vitro directed evolution strategies targeted at the D1 protein were adopted to create libraries of chlamydomonas random mutants, subsequently selected by exposures to radical-generating proton or neutron sources. The common trend observed in the D1 aminoacidic substitutions was the replacement of less polar by more polar amino acids. The applied selection pressure forced replacement of residues more sensitive to oxidative damage with less sensitive ones, suggesting that ionizing radiation may have been one of the driving forces in the evolution of the eukaryotic photosynthetic apparatus. A set of the identified aminoacidic substitutions, close to the secondary plastoquinone binding niche and oxygen evolving complex, were introduced by site-directed mutagenesis in un-transformed strains, and their sensitivity to free radicals attack analyzed. Mutants displayed reduced electron transport efficiency in physiological conditions, and increased photosynthetic performance stability and oxygen evolution capacity in stressful high-light conditions. Finally, comparative in silico analyses of D1 aminoacidic sequences of organisms differently located in the evolution chain, revealed a higher ratio of residues more sensitive to oxidative damage in the eukaryotic/cyanobacterial proteins compared to their bacterial orthologs. These results led us to hypothesize an archaean atmosphere less challenging in terms of ionizing radiation than the present one.

Networks of gene sharing among 329 proteobacterial genomes reveal differences in lateral gene transfer frequency at different phylogenetic depths

Lateral gene transfer (LGT) is an important mechanism of natural variation among prokaryotes. Over the full course of evolution, most or all of the genes resident in a given prokaryotic genome have been affected by LGT, yet the frequency of LGT can vary greatly across genes and across prokaryotic groups. The proteobacteria are among the most diverse of prokaryotic taxa. The prevalence of LGT in their genome evolution calls for the application of network-based methods instead of tree-based methods to investigate the relationships among these species. Here, we report networks that capture both vertical and horizontal components of evolutionary history among 1,207,272 proteins distributed across 329 sequenced proteobacterial genomes. The network of shared proteins reveals modularity structure that does not correspond to current classification schemes. On the basis of shared protein-coding genes, the five classes of proteobacteria fall into two main modules, one including the alpha-, delta-, and epsilonproteobacteria and the other including beta- and gammaproteobacteria. The first module is stable over different protein identity thresholds. The second shows more plasticity with regard to the sequence conservation of proteins sampled, with the gammaproteobacteria showing the most chameleon-like evolutionary characteristics within the present sample. Using a minimal lateral network approach, we compared LGT rates at different phylogenetic depths. In general, gene evolution by LGT within proteobacteria is very common. At least one LGT event was inferred to have occurred in at least 75% of the protein families. The average LGT rate at the species and class depth is about one LGT event per protein family, the rate doubling at the phylum level to an average of two LGT events per protein family. Hence, our results indicate that the rate of gene acquisition per protein family is similar at the level of species (by recombination) and at the level of classes (by LGT). The frequency of LGT per genome strongly depends on the species lifestyle, with endosymbionts showing far lower LGT frequencies than free-living species. Moreover, the nature of the transferred genes suggests that gene transfer in proteobacteria is frequently mediated by conjugation.

Molecular signatures for the main phyla of photosynthetic bacteria and their subgroups

The bacterial groups corresponding to different photosynthetic prokaryotes are presently identified mainly on the basis of their branching in phylogenetic trees. The availability of genome sequences is enabling identification of many molecular signatures that are specific for different groups of photosynthetic bacteria. Our recent work has identified large numbers of signatures consisting of conserved inserts or deletions (indels) in widely distributed proteins, as well as whole proteins that are specific for various sequenced species/strains from Cyanobacteria, Chlorobi, and Proteobacteria phyla. Based upon these signatures, it is now possible to identify/distinguish bacteria from these phyla of photosynthetic bacteria as well as their major subclades in clear molecular terms. The use of these signatures in conjunction with phylogenomic analyses, summarized here, is leading to a holistic picture concerning the branching order and evolutionary relationships among the above groups of photosynthetic bacteria. Although detailed studies in this regard have not yet been carried on Chloroflexi and Heliobacteriaceae, we have identified some conserved indels that are specific for these groups. Some of the conserved indels for the photosynthetic bacteria are present in photosynthesis-related proteins. These include a 4 aa insert in the pyruvate flavodoxin/ferridoxin oxidoreductase that is specific for the genus Chloroflexus, a 2 aa insert in magnesium chelatase that is uniquely shared by all Cyanobacteria except the deepest branching Clade A (Gloebacterales), a 6 aa insert in an A-type flavoprotein that is specific for various marine unicellular Cyanobacteria, a 2 aa insert in heme oxygenase that is specific for various Prochlorococcus strains/isolates, and 1 aa deletion in the protein protochlorophyllide oxidoreductase that is commonly shared by various Prochlorococcus strains except the deepest branching isolates MIT 9303 and MIT 9313. The identified CSIs are located in the structures of these proteins in surface loops indicating that they may be important in mediating protein-protein interactions. The cellular functions of these conserved indels, or most of the signature proteins are presently unknown, but they provide valuable means for discovering novel properties that are unique to different groups of photosynthetic bacteria.

Signature proteins for the major clades of Cyanobacteria

BACKGROUND

The phylogeny and taxonomy of cyanobacteria is currently poorly understood due to paucity of reliable markers for identification and circumscription of its major clades.

RESULTS

A combination of phylogenomic and protein signature based approaches was used to characterize the major clades of cyanobacteria. Phylogenetic trees were constructed for 44 cyanobacteria based on 44 conserved proteins. In parallel, Blastp searches were carried out on each ORF in the genomes of Synechococcus WH8102, Synechocystis PCC6803, Nostoc PCC7120, Synechococcus JA-3-3Ab, Prochlorococcus MIT9215 and Prochlor. marinus subsp. marinus CCMP1375 to identify proteins that are specific for various main clades of cyanobacteria. These studies have identified 39 proteins that are specific for all (or most) cyanobacteria and large numbers of proteins for other cyanobacterial clades. The identified signature proteins include: (i) 14 proteins for a deep branching clade (Clade A) of Gloebacter violaceus and two diazotrophic Synechococcus strains (JA-3-3Ab and JA2-3-B'a); (ii) 5 proteins that are present in all other cyanobacteria except those from Clade A; (iii) 60 proteins that are specific for a clade (Clade C) consisting of various marine unicellular cyanobacteria (viz. Synechococcus and Prochlorococcus); (iv) 14 and 19 signature proteins that are specific for the Clade C Synechococcus and Prochlorococcus strains, respectively; (v) 67 proteins that are specific for the Low B/A ecotype Prochlorococcus strains, containing lower ratio of chl b/a2 and adapted to growth at high light intensities; (vi) 65 and 8 proteins that are specific for the Nostocales and Chroococcales orders, respectively; and (vii) 22 and 9 proteins that are uniquely shared by various Nostocales and Oscillatoriales orders, or by these two orders and the Chroococcales, respectively. We also describe 3 conserved indels in flavoprotein, heme oxygenase and protochlorophyllide oxidoreductase proteins that are specific for either Clade C cyanobacteria or for various subclades of Prochlorococcus. Many other conserved indels for cyanobacterial clades have been described recently.

CONCLUSIONS

These signature proteins and indels provide novel means for circumscription of various cyanobacterial clades in clear molecular terms. Their functional studies should lead to discovery of novel properties that are unique to these groups of cyanobacteria.

Large-scale parsimony analysis of metazoan indels in protein-coding genes

Insertions and deletions (indels) are considered to be rare evolutionary events, the analysis of which may resolve controversial phylogenetic relationships. Indeed, indel characters are often assumed to be less homoplastic than amino acid and nucleotide substitutions and, consequently, more reliable markers for phylogenetic reconstruction. In this study, we analyzed indels from over 1,000 metazoan orthologous genes. We studied the impact of different species sampling, ortholog data sets, lengths of included indels, and indel-coding methods on the resulting metazoan tree. Our results show that, similar to sequence substitutions, indels are homoplastic characters, and their analysis is sensitive to the long-branch attraction artifact. Furthermore, improving the taxon sampling and choosing a closely related outgroup greatly impact the phylogenetic inference. Our indel-based inferences support the Ecdysozoa hypothesis over the Coelomata hypothesis and suggest that sponges are a sister clade to other animals.

Evolutionary placement of Xanthomonadales based on conserved protein signature sequences

Xanthomonadales comprises one of the largest phytopathogenic bacterial groups, and is currently classified within the gamma-proteobacteria. However, the phylogenetic placement of this group is not clearly resolved, and the results of different studies contradict one another. In this work, the evolutionary position of Xanthomonadales was determined by analyzing the presence of shared insertions and deletions (INDELs) in highly conserved proteins. Several distinctive insertions found in most of the members of the gamma-proteobacteria are absent in Xanthomonadales and groups such as Legionelalles, Chromatiales, Methylococcales, Thiotrichales and Cardiobacteriales. These INDELs were most likely introduced after the branching of Xanthomonadales from most of the gamma-proteobacteria and provide evidence for the phylogenetic placement of the early gamma-proteobacteria. Moreover, other proteins contain insertions exclusive to the Xanthomonadales order, confirming that this is a monophyletic group and provide important specific genetic markers. Thus, the data presented clearly support the Xanthomonadales group as an independent subdivision, and constitute one of the deepest branching lineage within the gamma-proteobacteria clade.

Phylogenomic analyses of clostridia and identification of novel protein signatures that are specific to the genus Clostridium sensu stricto (cluster I)

The species of Clostridium comprise a very heterogeneous assemblage of bacteria that do not form a phylogenetically coherent group. It has been proposed previously that only a subset of the species of Clostridium that form a distinct cluster in the 16S rRNA tree (cluster I) should be regarded as the true representatives of the genus Clostridium (i.e. Clostridium sensu stricto). However, this cluster is presently defined only in phylogenetic terms, and no biochemical, molecular or phenotypic characteristic is known that is unique to species from this cluster. We report here phylogenomic and comparative analyses based on sequenced clostridial genomes in an attempt to bridge this gap and to clarify the evolutionary relationships among species of clostridia. In phylogenetic trees for species of clostridia based on concatenated sequences for 37 highly conserved proteins, the species of Clostridium cluster I formed a strongly supported clade that was separated from all other clostridia by a long branch. Several other Clostridium species that are not part of this cluster grouped reliably with other species of clostridia in a number of well-resolved clades. Our comparative genomic analyses have identified three conserved indels in three highly conserved proteins (a 4 aa insert in DNA gyrase A, a 1 aa deletion in ATP synthase beta subunit and a 1 aa insert in ribosomal protein S2) that are unique to the species of Clostridium cluster I and are not found in any other bacteria. blastp searches on various proteins in the genomes of Clostridium tetani E88 and Clostridium perfringens SM101 have also identified more than 10 proteins that are found uniquely in the cluster I species. These results provide evidence that the species of Clostridium cluster I not only are phylogenetically distinct but also share many unique molecular characteristics. These newly identified molecular markers provide useful tools to define and circumscribe the genus Clostridium sensu stricto in more definitive terms. We have also identified a 7-9 aa conserved insert in the enzyme phosphoglycerate dehydrogenase that is uniquely found in the Clostridium thermocellum, Thermoanaerobacter pseudethanolicus, Thermoanaerobacter tengcogensis and Caldicellulosiruptor saccharolyticus homologues, and is absent from all other bacteria. These species form a well-defined clade in the phylogenetic trees and this indel provides a potential molecular marker for this clostridial cluster.

Phylogenomics and signature proteins for the alpha proteobacteria and its main groups

Background

Alpha proteobacteria are one of the largest and most extensively studied groups within bacteria. However, for these bacteria as a whole and for all of its major subgroups (viz. Rhizobiales, Rhodobacterales, Rhodospirillales, Rickettsiales, Sphingomonadales and Caulobacterales), very few or no distinctive molecular or biochemical characteristics are known.

Results

We have carried out comprehensive phylogenomic analyses by means of Blastp and PSI-Blast searches on the open reading frames in the genomes of several α-proteobacteria (viz. Bradyrhizobium japonicum, Brucella suis, Caulobacter crescentus, Gluconobacter oxydans, Mesorhizobium loti, Nitrobacter winogradskyi, Novosphingobium aromaticivorans, Rhodobacter sphaeroides 2.4.1, Silicibacter sp. TM1040, Rhodospirillum rubrum and Wolbachia (Drosophila) endosymbiont). These studies have identified several proteins that are distinctive characteristics of all α-proteobacteria, as well as numerous proteins that are unique repertoires of all of its main orders (viz. Rhizobiales, Rhodobacterales, Rhodospirillales, Rickettsiales, Sphingomonadales and Caulobacterales) and many families (viz. Rickettsiaceae, Anaplasmataceae, Rhodospirillaceae, Acetobacteraceae, Bradyrhiozobiaceae, Brucellaceae and Bartonellaceae). Many other proteins that are present at different phylogenetic depths in α-proteobacteria provide important information regarding their evolution. The evolutionary relationships among α-proteobacteria as deduced from these studies are in excellent agreement with their branching pattern in the phylogenetic trees and character compatibility cliques based on concatenated sequences for many conserved proteins. These studies provide evidence that the major groups within α-proteobacteria have diverged in the following order: (Rickettsiales(Rhodospirillales (Sphingomonadales (Rhodobacterales (Caulobacterales-Parvularculales (Rhizobiales)))))). We also describe two conserved inserts in DNA Gyrase B and RNA polymerase beta subunit that are distinctive characteristics of the Sphingomonadales and Rhodosprilllales species, respectively. The results presented here also provide support for the grouping of Hyphomonadaceae and Parvularcula species with the Caulobacterales and the placement of Stappia aggregata with the Rhizobiaceae group.

Conclusion

The α-proteobacteria-specific proteins and indels described here provide novel and powerful means for the taxonomic, biochemical and molecular biological studies on these bacteria. Their functional studies should prove helpful in identifying novel biochemical and physiological characteristics that are unique to these bacteria.

Ten years after the first Helicobacter pylori genome: comparative and functional genomics provide new insights in the variability and adaptability of a persistent pathogen

In this review, we summarize how genomic approaches contributed to the understanding of the biology of the recently discovered pathogen Helicobacter pylori. Comparative genomics provided new insights into H. pylori's spectacular genetic diversity and generated exiting hypotheses on its evolutionary history. Transcriptomic studies provided original information on the mechanisms of H. pylori gastric adaptation that are central to its virulence.

Phylogeny and molecular signatures (conserved proteins and indels) that are specific for the Bacteroidetes and Chlorobi species

BACKGROUND

The Bacteroidetes and Chlorobi species constitute two main groups of the Bacteria that are closely related in phylogenetic trees. The Bacteroidetes species are widely distributed and include many important periodontal pathogens. In contrast, all Chlorobi are anoxygenic obligate photoautotrophs. Very few (or no) biochemical or molecular characteristics are known that are distinctive characteristics of these bacteria, or are commonly shared by them.

RESULTS

Systematic blast searches were performed on each open reading frame in the genomes of Porphyromonas gingivalis W83, Bacteroides fragilis YCH46, B. thetaiotaomicron VPI-5482, Gramella forsetii KT0803, Chlorobium luteolum (formerly Pelodictyon luteolum) DSM 273 and Chlorobaculum tepidum (formerly Chlorobium tepidum) TLS to search for proteins that are uniquely present in either all or certain subgroups of Bacteroidetes and Chlorobi. These studies have identified > 600 proteins for which homologues are not found in other organisms. This includes 27 and 51 proteins that are specific for most of the sequenced Bacteroidetes and Chlorobi genomes, respectively; 52 and 38 proteins that are limited to species from the Bacteroidales and Flavobacteriales orders, respectively, and 5 proteins that are common to species from these two orders; 185 proteins that are specific for the Bacteroides genus. Additionally, 6 proteins that are uniquely shared by species from the Bacteroidetes and Chlorobi phyla (one of them also present in the Fibrobacteres) have also been identified. This work also describes two large conserved inserts in DNA polymerase III (DnaE) and alanyl-tRNA synthetase that are distinctive characteristics of the Chlorobi species and a 3 aa deletion in ClpB chaperone that is mainly found in various Bacteroidales, Flavobacteriales and Flexebacteraceae, but generally not found in the homologs from other organisms. Phylogenetic analyses of the Bacteroidetes and Chlorobi species is also reported based on concatenated sequences for 12 conserved proteins by different methods including the character compatibility (or clique) approach. The placement of Salinibacter ruber with other Bacteroidetes species was not resolved by other phylogenetic methods, but this affiliation was strongly supported by the character compatibility approach.

CONCLUSION

The molecular signatures described here provide novel tools for identifying and circumscribing species from the Bacteroidetes and Chlorobi phyla as well as some of their main groups in clear terms. These results also provide strong evidence that species from these two phyla (and also possibly Fibrobacteres) are specifically related to each other and they form a single superphylum. Functional studies on these proteins and indels should aid in the discovery of novel biochemical and physiological characteristics that are unique to these groups of bacteria.

Phylogenomic analysis of proteins that are distinctive of Archaea and its main subgroups and the origin of methanogenesis

BACKGROUND

The Archaea are highly diverse in terms of their physiology, metabolism and ecology. Presently, very few molecular characteristics are known that are uniquely shared by either all archaea or the different main groups within archaea. The evolutionary relationships among different groups within the Euryarchaeota branch are also not clearly understood.

RESULTS

We have carried out comprehensive analyses on each open reading frame (ORFs) in the genomes of 11 archaea (3 Crenarchaeota--Aeropyrum pernix, Pyrobaculum aerophilum and Sulfolobus acidocaldarius; 8 Euryarchaeota--Pyrococcus abyssi, Methanococcus maripaludis, Methanopyrus kandleri, Methanococcoides burtonii, Halobacterium sp. NCR-1, Haloquadratum walsbyi, Thermoplasma acidophilum and Picrophilus torridus) to search for proteins that are unique to either all Archaea or for its main subgroups. These studies have identified 1448 proteins or ORFs that are distinctive characteristics of Archaea and its various subgroups and whose homologues are not found in other organisms. Six of these proteins are unique to all Archaea, 10 others are only missing in Nanoarchaeum equitans and a large number of other proteins are specific for various main groups within the Archaea (e.g. Crenarchaeota, Euryarchaeota, Sulfolobales and Desulfurococcales, Halobacteriales, Thermococci, Thermoplasmata, all methanogenic archaea or particular groups of methanogens). Of particular importance is the observation that 31 proteins are uniquely present in virtually all methanogens (including M. kandleri) and 10 additional proteins are only found in different methanogens as well as A. fulgidus. In contrast, no protein was exclusively shared by various methanogen and any of the Halobacteriales or Thermoplasmatales. These results strongly indicate that all methanogenic archaea form a monophyletic group exclusive of other archaea and that this lineage likely evolved from Archaeoglobus. In addition, 15 proteins that are uniquely shared by M. kandleri and Methanobacteriales suggest a close evolutionary relationship between them. In contrast to the phylogenomics studies, a monophyletic grouping of archaea is not supported by phylogenetic analyses based on protein sequences.

CONCLUSION

The identified archaea-specific proteins provide novel molecular markers or signature proteins that are distinctive characteristics of Archaea and all of its major subgroups. The species distributions of these proteins provide novel insights into the evolutionary relationships among different groups within Archaea, particularly regarding the origin of methanogenesis. Most of these proteins are of unknown function and further studies should lead to discovery of novel biochemical and physiological characteristics that are unique to either all archaea or its different subgroups.

Application of the character compatibility approach to generalized molecular sequence data: branching order of the proteobacterial subdivisions

The character compatibility approach, which removes all homoplasic characters and involves finding the largest clique of compatible characters in a dataset, in principle, provides a powerful means for obtaining correct topology in difficult to resolve cases. However, the usefulness of this approach to generalized molecular sequence data for phylogeny determination has not been studied in the past. We have used this approach to determine the topology of 23 proteobacterial species (6 each of alpha-, beta- and gamma-, 3 delta-, and 2 epsilon-proteobacteria) using sequence data for 10 conserved proteins (Hsp60, Hsp70, EF-Tu, EF-G, alanyl-tRNA synthetase, RecA, GyrA, GyrB, RpoB and RpoC). All sites in the sequence alignments of these proteins where only two amino acids were found, with each amino acid present in at least two species, were selected. Mutual compatibility determination on these binary state sites was carried out by two means. In one case, all of these sites were combined into a large dataset (Set A; 957 characters) prior to compatibility analysis. In the second case, compatibility analysis was carried out on characters from individual proteins and all compatible sites were combined into a large dataset (Set B; 398 characters) for further studies. Upon compatibility analyses, the largest cliques that were obtained from Sets A and B consisted of 337 and 323 compatible characters, respectively. In these cliques, all proteobacterial subgroups were clearly distinguished and branching orders of most of the species were also resolved. The epsilon-proteobacteria exhibited the earliest branching, whereas the beta- and gamma-subgroups were found to have emerged last. The relative placement of the alpha- and delta-subgroups, however, was not resolved. The topology of these species was also determined based on 16S rRNA sequences and a concatenated dataset of sequences for all 10 proteins by means of neighbor-joining, maximum likelihood, and maximum parsimony methods. In the protein trees, all proteobacterial groups were reliably resolved and they branched in the following order: (epsilon(delta(alpha(beta,gamma)))). However, in the rRNA trees, the gamma- and beta-subgroups exhibited polyphyletic branching and many internal nodes were not resolved. These results indicate that the character compatibility analysis using generalized molecular sequence data provides a powerful means for evolutionary studies. Based on molecular sequences, it should be possible to obtain very large datasets of compatible characters that should prove very helpful in clarifying difficult to resolve phylogenetic relationships.