Gene decay in Shigella as an incipient stage of host-adaptation

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Pathogenicity and virulence of Shigella sonnei: A highly drug-resistant pathogen of increasing prevalence

Shigella spp. are the causative agent of shigellosis (or bacillary dysentery), a diarrhoeal disease characterized for the bacterial invasion of gut epithelial cells. Among the 4 species included in the genus, Shigella flexneri is principally responsible for the disease in the developing world while Shigella sonnei is the main causative agent in high-income countries. Remarkably, as more countries improve their socioeconomic conditions, we observe an increase in the relative prevalence of S. sonnei. To date, the reasons behind this change in aetiology depending on economic growth are not understood. S. flexneri has been widely used as a model to study the pathogenesis of the genus, but as more research data are collected, important discrepancies with S. sonnei have come to light. In comparison to S. flexneri, S. sonnei can be differentiated in numerous aspects; it presents a characteristic O-antigen identical to that of one serogroup of the environmental bacterium Plesiomonas shigelloides, a group 4 capsule, antibacterial mechanisms to outcompete and displace gut commensal bacteria, and a poorer adaptation to an intracellular lifestyle. In addition, the World Health Organization (WHO) have recognized the significant threat posed by antibiotic-resistant strains of S. sonnei, demanding new approaches. This review gathers knowledge on what is known about S. sonnei within the context of other Shigella spp. and aims to open the door for future research on understanding the increasing spread of this pathogen.

Role of the MDR Efflux Pump AcrAB in Epithelial Cell Invasion by Shigella flexneri

The tripartite complex AcrAB-TolC is the major RND pump in Escherichia coli and other Enterobacteriaceae, including Shigella, the etiological agent of bacillary dysentery. In addition to conferring resistance to many classes of antibiotics, AcrAB plays a role in the pathogenesis and virulence of several bacterial pathogens. Here, we report data demonstrating that AcrAB specifically contributes to Shigella flexneri invasion of epithelial cells. We found that deletion of both acrA and acrB genes causes reduced survival of S. flexneri M90T strain within Caco-2 epithelial cells and prevents cell-to-cell spread of the bacteria. Infections with single deletion mutant strains indicate that both AcrA and AcrB favor the viability of the intracellular bacteria. Finally, we were able to further confirm the requirement of the AcrB transporter activity for intraepithelial survival by using a specific EP inhibitor. Overall, the data from the present study expand the role of the AcrAB pump to an important human pathogen, such as Shigella, and add insights into the mechanism governing the Shigella infection process.

Characterization of Shigella flexneri Serotype 6 Strains Isolated from Bangladesh and Identification of a New Phylogenetic Cluster

A significant cause of shigellosis in Bangladesh and other developing countries is Shigella flexneri serotype 6. This serotype has been subtyped, on the basis of the absence or presence of a group-specific antigen, E1037, into S. flexneri 6a and 6b, respectively. Here, we provided rationales for the subclassification, using several phenotypic and molecular tools. A set of S. flexneri 6a and 6b strains isolated between 1997 and 2015 were characterized by analyzing their biochemical properties, plasmid profiles, virulence markers, pulsed-field gel electrophoresis (PFGE) results, and ribotype. Additionally, the genomic relatedness of these subserotypes was investigated with global isolates of serotype 6 using publicly available genomes. Both subserotypes of S. flexneri 6 agglutinated with monoclonal antiserum against S. flexneri (MASF) B and type VI-specific antiserum (MASF VI) and were PCR positive for O-antigen flippase-specific genes and virulence markers (ipaH, ial, sen, and sigA). Unlike S. flexneri 6a strains, S. flexneri 6b strains seroagglutinated with anti-E1037 antibodies, MASF IV-I. Notably, these two antigenically distinct subserotypes were clonally diverse, showing two distinct PFGE patterns following the digestion of chromosomal DNA with either XbaI or IceuI. In addition, hybridization of a 16S rRNA gene probe with HindIII-digested genomic DNA yielded two distinguishing ribotypes. Genomic comparison of S. flexneri subserotype 6a and 6b strains from Bangladesh indicated that, although these strains were in genomic synteny, the majority of them formed a unique phylogroup (PG-4) that was missing for the global isolates. This study supports the subserotyping and emphasizes the need for global monitoring of the S. flexneri subserotypes 6a and 6b. IMPORTANCE Shigella flexneri serotype 6 is one of the predominant serotypes among shigellosis cases in Bangladesh. Characterization of a novel subserotype of S. flexneri 6 (VI:E1037), agglutinated with type 6-specific antibody and anti-E1037, indicates a unique evolutionary ancestry. PFGE genotyping supports the finding that these two antigenically distinct subserotypes are clonally diverse. A phylogenetic study based on single-nucleotide polymorphism (SNP) data revealed that these two subserotypes were in genomic synteny, although their genomes were reduced. Interestingly, a majority of the S. flexneri 6 strains isolated from Bangladesh form a novel phylogenetic cluster. Therefore, this report underpins the global monitoring and tracking of the novel subserotype.

Chromosome-encoded IpaH ubiquitin ligases indicate non-human enteroinvasive Escherichia

Until recently, Shigella and enteroinvasive Escherichia coli were thought to be primate-restricted pathogens. The base of their pathogenicity is the type 3 secretion system (T3SS) encoded by the pINV virulence plasmid, which facilitates host cell invasion and subsequent proliferation. A large family of T3SS effectors, E3 ubiquitin-ligases encoded by the ipaH genes, have a key role in the Shigella pathogenicity through the modulation of cellular ubiquitination that degrades host proteins. However, recent genomic studies identified ipaH genes in the genomes of Escherichia marmotae, a potential marmot pathogen, and an E. coli extracted from fecal samples of bovine calves, suggesting that non-human hosts may also be infected by these strains, potentially pathogenic to humans. We performed a comparative genomic study of the functional repertoires in the ipaH gene family in Shigella and enteroinvasive Escherichia from human and predicted non-human hosts. We found that fewer than half of Shigella genomes had a complete set of ipaH genes, with frequent gene losses and duplications that were not consistent with the species tree and nomenclature. Non-human host IpaH proteins had a diverse set of substrate-binding domains and, in contrast to the Shigella proteins, two variants of the NEL C-terminal domain. Inconsistencies between strains phylogeny and composition of effectors indicate horizontal gene transfer between E. coli adapted to different hosts. These results provide a framework for understanding of ipaH-mediated host-pathogens interactions and suggest a need for a genomic study of fecal samples from diseased animals.

Phylogenetic Relatedness and Genome Structure of Yersinia ruckeri Revealed by Whole Genome Sequencing and a Comparative Analysis

Yersinia ruckeri is the causative agent of enteric redmouth disease (ERM), a serious infection that affects global aquaculture with high economic impact. The present study used whole genome sequences to perform a comparative analysis on 10 Y. ruckeri strains and to explore their genetic relatedness to other members of the genus. Y. ruckeri, Yersinia entomophaga, and Yersinia nurmii formed a species complex that constitutes the most basal lineage of the genus. The results showed that the taxonomy of Y. ruckeri strains is better defined by using a core genome alignment and phylogenetic analysis. The distribution of accessory genes in all Yersinia species revealed the presence of 303 distinctive genes in Y. ruckeri. Of these, 169 genes were distributed in 17 genomic islands potentially involved in the pathogenesis of ERM via (1) encoding virulence factors such as Afp18, Yrp1, phage proteins and (2) improving the metabolic capabilities by enhancing utilization and metabolism of iron, amino acids (specifically, arginine and histidine), and carbohydrates. The genome of Y. ruckeri is highly conserved regarding gene structure, gene layout and functional categorization of genes. It contains various components of mobile genetic elements but lacks the CRISPR-Cas system and possesses a stable set of virulence genes possibly playing a critical role in pathogenicity. Distinct virulence plasmids were exclusively restricted to a specific clonal group of Y. ruckeri (CG4), possibly indicating a selective advantage. Phylogenetic analysis of Y. ruckeri genomes revealed the co-presence of multiple genetically distant lineages of Y. ruckeri strains circulating in Germany. Our results also suggest a possible dissemination of a specific group of strains in the United States, Peru, Germany, and Denmark. In conclusion, this study provides new insights into the taxonomy and evolution of Y. ruckeri and contributes to a better understanding of the pathogenicity of ERM in aquaculture. The genomic analysis presented here offers a framework for the development of more efficient control strategies for this pathogen.

Radical genome remodelling accompanied the emergence of a novel host-restricted bacterial pathogen

The emergence of new pathogens is a major threat to public and veterinary health. Changes in bacterial habitat such as a switch in host or disease tropism are typically accompanied by genetic diversification. Staphylococcus aureus is a multi-host bacterial species associated with human and livestock infections. A microaerophilic subspecies, Staphylococcus aureus subsp. anaerobius, is responsible for Morel’s disease, a lymphadenitis restricted to sheep and goats. However, the evolutionary history of S. aureus subsp. anaerobius and its relatedness to S. aureus are unknown. Population genomic analyses of clinical S. aureus subsp. anaerobius isolates revealed a highly conserved clone that descended from a S. aureus progenitor about 1000 years ago before differentiating into distinct lineages that contain African and European isolates. S. aureus subsp. anaerobius has undergone limited clonal expansion, with a restricted population size, and an evolutionary rate 10-fold slower than S. aureus. The transition to its current restricted ecological niche involved acquisition of a pathogenicity island encoding a ruminant host-specific effector of abscess formation, large chromosomal re-arrangements, and the accumulation of at least 205 pseudogenes, resulting in a highly fastidious metabolism. Importantly, expansion of ~87 insertion sequences (IS) located largely in intergenic regions provided distinct mechanisms for the control of expression of flanking genes, including a novel mechanism associated with IS-mediated anti-anti-sense decoupling of ancestral gene repression. Our findings reveal the remarkable evolutionary trajectory of a host-restricted bacterial pathogen that resulted from extensive remodelling of the S. aureus genome through an array of diverse mechanisms in parallel.

The emergence of new pathogens is a major threat to public and veterinary health. Some bacteria such as Staphylococcus aureus, have the capacity to infect many different host species including humans and livestock while others such as the closely-related S. aureus subsp. anaerobius, associated with a single type of pathology called Morel’s disease in small ruminants, are highly niche-restricted. However, our understanding of the genetic basis for such differences in bacterial host-tropism is very limited. Here, we discovered that S. aureus subsp. anaerobius evolved from an S. aureus ancestor and underwent an array of extensive changes to its genome that accompanied the transition to its current restricted lifestyle. We observed genome decay involving loss of function of hundreds of genes, large intra-chromosomal rearrangements affecting most of the genome, acquisition of a pathogenicity island, and expansion of large numbers of insertion sequences that are inserted at intergenic sites around the genome. Importantly, we found that IS elements affect the expression of neighbouring genes in different ways including a novel mechanism of IS-enabled disruption of ancestral gene repression. Taken together, we provide a remarkable example of radical genomic changes associated with evolutionary transition from a multi-host to highly restricted host ecology.

High Rates of Genome Rearrangements and Pathogenicity of Shigella spp

Shigella are pathogens originating within the Escherichia lineage but frequently classified as a separate genus. Shigella genomes contain numerous insertion sequences (ISs) that lead to pseudogenisation of affected genes and an increase of non-homologous recombination. Here, we study 414 genomes of E. coli and Shigella strains to assess the contribution of genomic rearrangements to Shigella evolution. We found that Shigella experienced exceptionally high rates of intragenomic rearrangements and had a decreased rate of homologous recombination compared to pathogenic and non-pathogenic E. coli. The high rearrangement rate resulted in independent disruption of syntenic regions and parallel rearrangements in different Shigella lineages. Specifically, we identified two types of chromosomally encoded E3 ubiquitin-protein ligases acquired independently by all Shigella strains that also showed a high level of sequence conservation in the promoter and further in the 5′-intergenic region. In the only available enteroinvasive E. coli (EIEC) strain, which is a pathogenic E. coli with a phenotype intermediate between Shigella and non-pathogenic E. coli, we found a rate of genome rearrangements comparable to those in other E. coli and no functional copies of the two Shigella-specific E3 ubiquitin ligases. These data indicate that the accumulation of ISs influenced many aspects of genome evolution and played an important role in the evolution of intracellular pathogens. Our research demonstrates the power of comparative genomics-based on synteny block composition and an important role of non-coding regions in the evolution of genomic islands.

Comparative genomics of a novel clade shed light on the evolution of the genus Erysipelothrix and characterise an emerging species

Erysipelothrix sp. isolates obtained from a deadly outbreak in farmed turkeys were sequenced and compared to representatives of the genus. Phylogenetic trees—supported by digital DNA:DNA hybridization and Average Nucleotide Identity—revealed a novel monophyletic clade comprising isolates from pigs, turkeys, and fish, including isolates previously described as E. sp. Strain 2. Genes coding for the SpaC protein, typically found in E. sp. Strain 2, were detected in all isolates of the clade. Therefore, we confirm E. sp. Strain 2 represents a unique species, that despite its official name “Erysipelothrix piscisicarius” (meaning a killer of fish), may be isolated from a broad host range. Core genome analysis showed that the pathogenic species of this genus, E. rhusiopathiae and the clade E. sp. Strain 2, are enriched in core functionalities related to nutrient uptake and transport, but not necessarily homologous pathways. For instance, whereas the aerobic DctA transporter may uptake C₄-dicarboxylates in both species, the anaerobic DcuC transporter is exclusive of the E. sp. Strain 2. Remarkably, the pan-genome analysis uncovered that genes related to transport and metabolism, recombination and repair, translation and transcription in the fish isolate, within the novel clade, have undergone a genomic reduction through pseudogenization. This reflects distinct selective pressures shaping the genome of species and strains within the genus Erysipelothrix while adapting to their respective niches.

A One Health Perspective for Defining and Deciphering Escherichia coli Pathogenic Potential in Multiple Hosts

E. coli is one of the most common species of bacteria colonizing humans and animals. The singularity of E. coli 's genus and species underestimates its multifaceted nature, which is represented by different strains, each with different combinations of distinct virulence factors. In fact, several E. coli pathotypes, or hybrid strains, may be associated with both subclinical infection and a range of clinical conditions, including enteric, urinary, and systemic infections. E. coli may also express DNA-damaging toxins that could impact cancer development. This review summarizes the different E. coli pathotypes in the context of their history, hosts, clinical signs, epidemiology, and control. The pathotypic characterization of E. coli in the context of disease in different animals, including humans, provides comparative and One Health perspectives that will guide future clinical and research investigations of E. coli infections.

Campylobacter jejuni and Campylobacter coli autotransporter genes exhibit lineage-associated distribution and decay

Background

Campylobacter jejuni and Campylobacter coli are major global causes of bacterial gastroenteritis. Whilst several individual colonisation and virulence factors have been identified, our understanding of their role in the transmission, pathogenesis and ecology of Campylobacter has been hampered by the genotypic and phenotypic diversity within C. jejuni and C. coli. Autotransporter proteins are a family of outer membrane or secreted proteins in Gram-negative bacteria such as Campylobacter, which are associated with virulence functions. In this study we have examined the distribution and predicted functionality of the previously described capC and the newly identified, related capD autotransporter gene families in Campylobacter.

Results

Two capC-like autotransporter families, designated capC and capD, were identified by homology searches of genomes of the genus Campylobacter. Each family contained four distinct orthologs of CapC and CapD. The distribution of these autotransporter genes was determined in 5829 C. jejuni and 1347 C. coli genomes. Autotransporter genes were found as intact, complete copies and inactive formats due to premature stop codons and frameshift mutations. Presence of inactive and intact autotransporter genes was associated with C. jejuni and C. coli multi-locus sequence types, but for capC, inactivation was independent from the length of homopolymeric tracts in the region upstream of the capC gene. Inactivation of capC or capD genes appears to represent lineage-specific gene decay of autotransporter genes. Intact capC genes were predominantly associated with the C. jejuni ST-45 and C. coli ST-828 generalist lineages. The capD3 gene was only found in the environmental C. coli Clade 3 lineage. These combined data support a scenario of inter-lineage and interspecies exchange of capC and subsets of capD autotransporters.

Conclusions

In this study we have identified two novel, related autotransporter gene families in the genus Campylobacter, which are not uniformly present and exhibit lineage-specific associations and gene decay. The distribution and decay of the capC and capD genes exemplifies the erosion of species barriers between certain lineages of C. jejuni and C. coli, probably arising through co-habitation. This may have implications for the phenotypic variability of these two pathogens and provide opportunity for new, hybrid genotypes to emerge.

Complete genome sequence and annotation of the laboratory reference strain Shigella flexneri serotype 5a M90T and genome-wide transcriptional start site determination

Background

Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches.

Results

We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome.

Conclusions

We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

The MFS efflux pump EmrKY contributes to the survival of Shigella within macrophages

Efflux pumps are membrane protein complexes conserved in all living organisms. Beyond being involved in antibiotic extrusion in several bacteria, efflux pumps are emerging as relevant players in pathogen-host interactions. We have investigated on the possible role of the efflux pump network in Shigella flexneri, the etiological agent of bacillary dysentery. We have found that S. flexneri has retained 14 of the 20 pumps characterized in Escherichia coli and that their expression is differentially modulated during the intracellular life of Shigella. In particular, the emrKY operon, encoding an efflux pump of the Major Facilitator Superfamily, is specifically and highly induced in Shigella-infected U937 macrophage-like cells and is activated in response to a combination of high K⁺ and acidic pH, which are sensed by the EvgS/EvgA two-component system. Notably, we show that following S. flexneri infection, macrophage cytosol undergoes a mild reduction of intracellular pH, permitting EvgA to trigger the emrKY activation. Finally, we present data suggesting that EmrKY is required for the survival of Shigella in the harsh macrophage environment, highlighting for the first time the key role of an efflux pump during the Shigella invasive process.

Polymorphisms in the CHIT1 gene: Associations with colorectal cancer

Colorectal cancer (CRC) is one of the most common solid tumors worldwide, often associated with inflammation. The microbes in the human intestine have a key role in inflammations and CRC. Chitotriose renders growth advantage to some bacteria, especially some pathogens, and thus has a role in inflammations. The enzyme chitotriosidase, encoded by the CHIT1 gene of the host, may degrade chitotriose with different efficiencies depending on the alleles. We sequenced the CHIT1 gene for 320 Chinese Han CRC patients and 404 normal controls, and focused on variations rs61745299 and rs35920428 within the CHIT1 gene for their possible roles in CRC. Statistical analyses were conducted using Chi-Square Tests as implemented in SPSS (version 19.0). Multiple sequence alignment was conducted using the Vector NTI, and protein expression levels were analyzed by western blotting. The two variations, rs61745299 and rs35920428 within the CDS region of CHIT1 gene, were associated with the risk of CRC (both with P values < 0.001). Western blotting analysis showed that the variations increased the expression levels of the CHIT1 and C-reaction protein genes in the cancer tissue. We conclude that the two variations of CHIT1, rs61745299 and rs35920428, increase expression of the gene and are associated with CRC in Chinese Han populations.

Advances in Understanding Bacterial Pathogenesis Gained from Whole-Genome Sequencing and Phylogenetics

The development of next-generation sequencing as a cost-effective technology has facilitated the analysis of bacterial population structure at a whole-genome level and at scale. From these data, phylogenic trees have been constructed that define population structures at a local, national, and global level, providing a framework for genetic analysis. Although still at an early stage, these approaches have yielded progress in several areas, including pathogen transmission mapping, the genetics of niche colonization and host adaptation, as well as gene-to-phenotype association studies. Antibiotic resistance has proven to be a major challenge in the early 21(st) century, and phylogenetic analyses have uncovered the dramatic effect that the use of antibiotics has had on shaping bacterial population structures. An update on insights into bacterial evolution from comparative genomics is provided in this review.

Genomic analysis of the marine fish pathogen Photobacterium damselae subsp. piscicida: Insertion sequences proliferation is associated with chromosomal reorganisations and rampant gene decay

Photobacterium damselae subsp. piscicida (Pdp) is an intracellular fish pathogen that causes photobacteriosis, a disease proven deadly in farmed fish worldwide. This work focuses on the analysis of genome sequences, chromosomes structure and gene contents of two strains from Sparus aurata (DI21) and Solea senegalensis (L091106-03H), isolated on the Spanish Atlantic coast. The comparative genomic analysis revealed that DI21 and L091106-03H share 98% of their genomes, including two virulence plasmids: pPHDP70 encoding siderophore piscibactin synthesis and pPHDP10 encoding the apoptotic toxin AIP56. Both genomes harbour a surprisingly large number of IS elements accounting for 12-17% of the total genome, representing an IS density of 0.15 elements per kb, one of the highest IS density values in a bacterial pathogen. This massive proliferation of ISs is responsible for the generation of a high number of pseudogenes that caused extensive loss of biological functions. Pseudogene formation is one of the main features of Pdp genome that explains most of the ecological and phenotypic differences with respect to its sibling subspecies P. damselae subsp. damselae and to other Vibrionaceae. Evidence was also found proving the existence of two chromosomal configurations depending on the origin of the strains: an European and an Asian/American types of genome organisation, reinforcing the idea of the existence of two geographically-linked clonal lineages in Pdp. In short, our study suggests that the host-dependent lifestyle of Pdp allowed massive IS proliferation and gene decay processes, which are major evolutionary forces in the shaping of the Pdp genome.

Are Escherichia coli Pathotypes Still Relevant in the Era of Whole-Genome Sequencing?

The empirical and pragmatic nature of diagnostic microbiology has given rise to several different schemes to subtype E.coli, including biotyping, serotyping, and pathotyping. These schemes have proved invaluable in identifying and tracking outbreaks, and for prognostication in individual cases of infection, but they are imprecise and potentially misleading due to the malleability and continuous evolution of E. coli. Whole genome sequencing can be used to accurately determine E. coli subtypes that are based on allelic variation or differences in gene content, such as serotyping and pathotyping. Whole genome sequencing also provides information about single nucleotide polymorphisms in the core genome of E. coli, which form the basis of sequence typing, and is more reliable than other systems for tracking the evolution and spread of individual strains. A typing scheme for E. coli based on genome sequences that includes elements of both the core and accessory genomes, should reduce typing anomalies and promote understanding of how different varieties of E. coli spread and cause disease. Such a scheme could also define pathotypes more precisely than current methods.

A profile-based method for identifying functional divergence of orthologous genes in bacterial genomes

Motivation: Next generation sequencing technologies have provided us with a wealth of information on genetic variation, but predicting the functional significance of this variation is a difficult task. While many comparative genomics studies have focused on gene flux and large scale changes, relatively little attention has been paid to quantifying the effects of single nucleotide polymorphisms and indels on protein function, particularly in bacterial genomics.

Results: We present a hidden Markov model based approach we call delta-bitscore (DBS) for identifying orthologous proteins that have diverged at the amino acid sequence level in a way that is likely to impact biological function. We benchmark this approach with several widely used datasets and apply it to a proof-of-concept study of orthologous proteomes in an investigation of host adaptation in Salmonella enterica. We highlight the value of the method in identifying functional divergence of genes, and suggest that this tool may be a better approach than the commonly used dN/dS metric for identifying functionally significant genetic changes occurring in recently diverged organisms.

Availability and Implementation: A program implementing DBS for pairwise genome comparisons is freely available at: https://github.com/UCanCompBio/deltaBS.

Contact:nicole.wheeler@pg.canterbury.ac.nz or lars.barquist@uni-wuerzburg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Mutational Studies of Putative Biosynthetic Genes for the Cyanobacterial Sunscreen Scytonemin in Nostoc punctiforme ATCC 29133

The heterocyclic indole-alkaloid scytonemin is a sunscreen found exclusively among cyanobacteria. An 18-gene cluster is responsible for scytonemin production in Nostoc punctiforme ATCC 29133. The upstream genes scyABCDEF in the cluster are proposed to be responsible for scytonemin biosynthesis from aromatic amino acid substrates. In vitro studies of ScyA, ScyB, and ScyC proved that these enzymes indeed catalyze initial pathway reactions. Here we characterize the role of ScyD, ScyE, and ScyF, which were logically predicted to be responsible for late biosynthetic steps, in the biological context of N. punctiforme. In-frame deletion mutants of each were constructed (ΔscyD, ΔscyE, and ΔscyF) and their phenotypes studied. Expectedly, ΔscyE presents a scytoneminless phenotype, but no accumulation of the predicted intermediaries. Surprisingly, ΔscyD retains scytonemin production, implying that it is not required for biosynthesis. Indeed, scyD presents an interesting evolutionary paradox: it likely originated in a duplication event from scyE, and unlike other genes in the operon, it has not been subjected to purifying selection. This would suggest that it is a pseudogene, and yet scyD is highly conserved in the scytonemin operon of cyanobacteria. ΔscyF also retains scytonemin production, albeit exhibiting a reduction of the production yield compared with the wild-type. This indicates that ScyF is not essential but may play an adjuvant role for scytonemin synthesis. Altogether, our findings suggest that these downstream genes are not responsible, as expected, for the late steps of scytonemin synthesis and we must look for those functions elsewhere. These findings are particularly important for biotechnological production of this sunscreen through heterologous expression of its genes in more tractable organisms.

Natural Selection in Virulence Genes of Francisella tularensis

A fundamental tenet of evolution is that alleles that are under negative selection are often deleterious and confer no evolutionary advantage. Negatively selected alleles are removed from the gene pool and are eventually extinguished from the population. Conversely, alleles under positive selection do confer an evolutionary advantage and lead to an increase in the overall fitness of the organism. These alleles increase in frequency until they eventually become fixed in the population. Francisella tularensis is a zoonotic pathogen and a potential biothreat agent. The most virulent type of F. tularensis, Type A, is distributed across North America with Type A.I occurring mainly in the east and Type A.II appearing mainly in the west. F. tularensis is thought to be a genome in decay (losing genes) because of the relatively large number of pseudogenes present in its genome. We hypothesized that the observed frequency of gene loss/pseudogenes may be an artifact of evolution in response to a changing environment, and that genes involved in virulence should be under strong positive selection. To test this hypothesis, we sequenced and compared whole genomes of Type A.I and A.II isolates. We analyzed a subset of virulence and housekeeping genes from several F. tularensis subspecies genomes to ascertain the presence and extent of positive selection. Eleven previously identified virulence genes were screened for positive selection along with 10 housekeeping genes. Analyses of selection yielded one housekeeping gene and 7 virulence genes which showed significant evidence of positive selection at loci implicated in cell surface structures and membrane proteins, metabolism and biosynthesis, transcription, translation and cell separation, and substrate binding and transport. Our results suggest that while the loss of functional genes through disuse could be accelerated by negative selection, the genome decay in Francisella could also be the byproduct of adaptive evolution driven by complex interactions between host, pathogen, and thier environment, as evidenced by several of its virulence genes which are undergoing strong, positive selection.

The genomic signatures of Shigella evolution, adaptation and geographical spread

Shigella spp. are some of the key pathogens responsible for the global burden of diarrhoeal disease. These facultative intracellular bacteria belong to the family Enterobacteriaceae, together with other intestinal pathogens, such as Escherichia coli and Salmonella spp. The genus Shigella comprises four different species, each consisting of several serogroups, all of which show phenotypic similarity, including invasive pathogenicity. DNA sequencing suggests that this similarity results from the convergent evolution of different Shigella spp. founders. Here, we review the evolutionary relationships between Shigella spp. and E . coli, and we highlight how the genomic plasticity of these bacteria and their acquisition of a distinctive virulence plasmid have enabled the development of such highly specialized pathogens. Furthermore, we discuss the insights that genotyping and whole-genome sequencing have provided into the phylogenetics and intercontinental spread of Shigella spp.

Trade-Offs of Escherichia coli Adaptation to an Intracellular Lifestyle in Macrophages

The bacterium Escherichia coli exhibits remarkable genomic and phenotypic variation, with some pathogenic strains having evolved to survive and even replicate in the harsh intra-macrophage environment. The rate and effects of mutations that can cause pathoadaptation are key determinants of the pace at which E. coli can colonize such niches and become pathogenic. We used experimental evolution to determine the speed and evolutionary paths undertaken by a commensal strain of E. coli when adapting to intracellular life. We estimated the acquisition of pathoadaptive mutations at a rate of 10⁻⁶ per genome per generation, resulting in the fixation of more virulent strains in less than a hundred generations. Whole genome sequencing of independently evolved clones showed that the main targets of intracellular adaptation involved loss of function mutations in genes implicated in the assembly of the lipopolysaccharide core, iron metabolism and di- and tri-peptide transport, namely rfaI, fhuA and tppB, respectively. We found a substantial amount of antagonistic pleiotropy in evolved populations, as well as metabolic trade-offs, commonly found in intracellular bacteria with reduced genome sizes. Overall, the low levels of clonal interference detected indicate that the first steps of the transition of a commensal E. coli into intracellular pathogens are dominated by a few pathoadaptive mutations with very strong effects.

Characterization of an emergent clone of enteroinvasive Escherichia coli circulating in Europe

Enteroinvasive Escherichia coli (EIEC) cause intestinal illness indistinguishable from that caused by Shigella, mainly in developing countries. Recently an upsurge of cases of EIEC infections has been observed in Europe, with two large outbreaks occurring in Italy and in the United Kingdom. We have characterized phenotypically and genotypically the strains responsible for these epidemics together with an additional isolate from a sporadic case isolated in Spain. The three isolates belonged to the same rare serotype O96:H19 and were of sequence type ST-99, never reported before in EIEC or Shigella. The EIEC strains investigated possessed all the virulence genes harboured on the large plasmid conferring the invasive phenotype to EIEC and Shigella while showing only some of the known chromosomal virulence genes and none of the described pathoadaptative mutations. At the same time, they displayed motility abilities and biochemical requirements resembling more closely those of the non-pathogenic E. coli rather than the EIEC and Shigella strains used as reference. Our observations suggested that the O96:H19 strains belong to an emerging EIEC clone, which could be the result of a recent event of acquisition of the invasion plasmid by commensal E. coli.

Multifactor Regulation of the MdtJI Polyamine Transporter in Shigella

The polyamine profile of Shigella, the etiological agent of bacillary dysentery in humans, differs markedly from that of E. coli, its innocuous commensal ancestor. Pathoadaptive mutations such as the loss of cadaverine and the increase of spermidine favour the full expression of the virulent phenotype of Shigella. Spermidine levels affect the expression of the MdtJI complex, a recently identified efflux pump belonging to the small multi-drug resistance family of transporters. In the present study, we have addressed the regulation of the mdtJI operon in Shigella by asking which factors influence its expression as compared to E. coli. In particular, after identifying the mdtJI promoter by primer extension analysis, in vivo transcription assays and gel-retardation experiments were carried out to get insight on the silencing of mdtJI in E. coli. The results indicate that H-NS, a major nucleoid protein, plays a key role in repressing the mdtJI operon by direct binding to the regulatory region. In the Shigella background mdtJI expression is increased by the high levels of spermidine typically found in this microorganism and by VirF, the plasmid-encoded regulator of the Shigella virulence regulatory cascade. We also show that the expression of mdtJI is stimulated by bile components. Functional analyses reveal that MdtJI is able to promote the excretion of putrescine, the spermidine precursor. This leads us to consider the MdtJI complex as a possible safety valve allowing Shigella to maintain spermidine to a level optimally suited to survival within infected macrophages and, at the same time, prevent toxicity due to spermidine over-accumulation.

Taking the pseudo out of pseudogenes

Pseudogenes are defined as fragments of once-functional genes that have been silenced by one or more nonsense, frameshift or missense mutations. Despite continuing increases in the speed of sequencing and annotating bacterial genomes, the identification and categorisation of pseudogenes remains problematic. Even when identified, pseudogenes are considered to be rare and tend to be ignored. On the contrary, pseudogenes are surprisingly prevalent and can persist for long evolutionary time periods, representing a record of once-functional genetic characteristics. Most importantly, pseudogenes provide an insight into prokaryotic evolutionary history as a record of phenotypic traits that have been lost. Focusing on the intracellular and symbiotic bacteria in which pseudogenes predominate, this review discusses the importance of identifying pseudogenes to fully understand the abilities of bacteria, and to understand prokaryotes within their evolutionary context.

Molecular and functional profiling of the polyamine content in enteroinvasive E. coli : looking into the gap between commensal E. coli and harmful Shigella

Polyamines are small molecules associated with a wide variety of physiological functions. Bacterial pathogens have developed subtle strategies to exploit polyamines or manipulate polyamine-related processes to optimize fitness within the host. During the transition from its innocuous E. coli ancestor, Shigella, the aetiological agent of bacillary dysentery, has undergone drastic genomic rearrangements affecting the polyamine profile. A pathoadaptation process involving the speG gene and the cad operon has led to spermidine accumulation and loss of cadaverine. While a higher spermidine content promotes the survival of Shigella within infected macrophages, the lack of cadaverine boosts the pathogenic potential of the bacterium in host tissues. Enteroinvasive E. coli (EIEC) display the same pathogenicity process as Shigella, but have a higher infectious dose and a higher metabolic activity. Pathoadaption events affecting the cad locus have occurred also in EIEC, silencing cadaverine production. Since EIEC are commonly regarded as evolutionary intermediates between E. coli and Shigella, we investigated on their polyamine profile in order to better understand which changes have occurred along the path to pathogenicity. By functional and molecular analyses carried out in EIEC strains belonging to different serotypes, we show that speG has been silenced in one strain only, favouring resistance to oxidative stress conditions and survival within macrophages. At the same time, we observe that the content of spermidine and putrescine, a relevant intermediate in the synthesis of spermidine, is higher in all strains as compared to E. coli. This may represent an evolutionary response to the lack of cadaverine. Indeed, restoring cadaverine synthesis decreases the expression of the speC gene, whose product affects putrescine production. In the light of these results, we discuss the possible impact of pathoadaptation events on the evolutionary emergence of a polyamine profile favouring to the pathogenic lifestyle of Shigella and EIEC.

Copy-number changes in evolution: rates, fitness effects and adaptive significance

Gene copy-number differences due to gene duplications and deletions are rampant in natural populations and play a crucial role in the evolution of genome complexity. Per-locus analyses of gene duplication rates in the pre-genomic era revealed that gene duplication rates are much higher than the per nucleotide substitution rate. Analyses of gene duplication and deletion rates in mutation accumulation lines of model organisms have revealed that these high rates of copy-number mutations occur at a genome-wide scale. Furthermore, comparisons of the spontaneous duplication and deletion rates to copy-number polymorphism data and bioinformatic-based estimates of duplication rates from sequenced genomes suggest that the vast majority of gene duplications are detrimental and removed by natural selection. The rate at which new gene copies appear in populations greatly influences their evolutionary dynamics and standing gene copy-number variation in populations. The opportunity for mutations that result in the maintenance of duplicate copies, either through neofunctionalization or subfunctionalization, also depends on the equilibrium frequency of additional gene copies in the population, and hence on the spontaneous gene duplication (and loss) rate. The duplication rate may therefore have profound effects on the role of adaptation in the evolution of duplicated genes as well as important consequences for the evolutionary potential of organisms. We further discuss the broad ramifications of this standing gene copy-number variation on fitness and adaptive potential from a population-genetic and genome-wide perspective.

Predicting genetic and ecological characteristics of bacterial species by comparing the trajectories of dN/dS and dI/dS in bacterial genomes

Indel (insertion/deletion) causes gene disruption and is considered to be deleterious like non-synonymous mutation during the evolution of bacterial genomes. The trajectory of dN/dS (the ratio of non-synonymous to synonymous mutation) has been found to decrease exponentially over time, but the trajectory of dI/dS (the ratio of indel to synonymous mutation) has not been thoroughly explored yet. Here we compared the patterns of dN/dS and dI/dS for several bacterial species. The majority of them showed a much steeper dI/dS trajectory than the dN/dS trajectory, suggesting that indel was more deleterious than non-synonymous mutation and therefore was more difficult to fix in genomes. However, the naturally competent bacteria, or those with a much lower genetic barrier for DNA exchange, presented a reverse relationship between dI/dS and dN/dS, indicative of an exceptional ability to tolerate horizontal genetic transfer. The result suggests that plotting of dN/dS and dI/dS trajectories can help to predict the genetic and ecological characteristics of bacterial species.

Recent advances in understanding enteric pathogenic Escherichia coli

Although Escherichia coli can be an innocuous resident of the gastrointestinal tract, it also has the pathogenic capacity to cause significant diarrheal and extraintestinal diseases. Pathogenic variants of E. coli (pathovars or pathotypes) cause much morbidity and mortality worldwide. Consequently, pathogenic E. coli is widely studied in humans, animals, food, and the environment. While there are many common features that these pathotypes employ to colonize the intestinal mucosa and cause disease, the course, onset, and complications vary significantly. Outbreaks are common in developed and developing countries, and they sometimes have fatal consequences. Many of these pathotypes are a major public health concern as they have low infectious doses and are transmitted through ubiquitous mediums, including food and water. The seriousness of pathogenic E. coli is exemplified by dedicated national and international surveillance programs that monitor and track outbreaks; unfortunately, this surveillance is often lacking in developing countries. While not all pathotypes carry the same public health profile, they all carry an enormous potential to cause disease and continue to present challenges to human health. This comprehensive review highlights recent advances in our understanding of the intestinal pathotypes of E. coli.

Partial disruption of translational and posttranslational machinery reshapes growth rates of Bartonella birtlesii

Specialization of bacteria in a new niche is associated with genome repertoire changes, and speciation in bacterial specialists is associated with genome reduction. Here, we tested a signature-tagged mutant library of 3,456 Bartonella birtlesii clones to detect mutants that could grow rapidly in vitro. Overall, we found 124 mutants that grew faster than the parental wild-type strain in vitro. We sequenced the genomes of the four mutants with the most rapid growth (formed visible colonies in only 1 to 2 days compared with 5 days for the wild type) and compared them to the parental isolate genome. We found that the number of disrupted genes associated with translation in the 124 rapid-growth clones was significantly higher than the number of genes involved in translation in the full genome (P < 10⁻⁶). Analysis of transposon integration in the genome of the four most rapidly growing clones revealed that one clone lacked one of the two wild-type RNA ribosomal operons. Finally, one of the four clones did not induce bacteremia in our mouse model, whereas infection with the other three resulted in a significantly lower bacterial count in blood than that with the wild-type strain.

Here, we show that specialization in a specific niche could be caused by the disruption of critical genes. Most of these genes were involved in translation, and we show that evolution of obligate parasitism bacteria was specifically associated with disruption of translation system-encoding genes.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

The chbG gene of the chitobiose (chb) operon of Escherichia coli encodes a chitooligosaccharide deacetylase

The chb operon of Escherichia coli is involved in the utilization of the β-glucosides chitobiose and cellobiose. The function of chbG (ydjC), the sixth open reading frame of the operon that codes for an evolutionarily conserved protein is unknown. We show that chbG encodes a monodeacetylase that is essential for growth on the acetylated chitooligosaccharides chitobiose and chitotriose but is dispensable for growth on cellobiose and chitosan dimer, the deacetylated form of chitobiose. The predicted active site of the enzyme was validated by demonstrating loss of function upon substitution of its putative metal-binding residues that are conserved across the YdjC family of proteins. We show that activation of the chb promoter by the regulatory protein ChbR is dependent on ChbG, suggesting that deacetylation of chitobiose-6-P and chitotriose-6-P is necessary for their recognition by ChbR as inducers. Strains carrying mutations in chbR conferring the ability to grow on both cellobiose and chitobiose are independent of chbG function for induction, suggesting that gain of function mutations in ChbR allow it to recognize the acetylated form of the oligosaccharides. ChbR-independent expression of the permease and phospho-β-glucosidase from a heterologous promoter did not support growth on both chitobiose and chitotriose in the absence of chbG, suggesting an additional role of chbG in the hydrolysis of chitooligosaccharides. The homologs of chbG in metazoans have been implicated in development and inflammatory diseases of the intestine, indicating that understanding the function of E. coli chbG has a broader significance.

Evolution of variation in presence and absence of genes in bacterial pathways

Background

Bacterial genomes exhibit a remarkable degree of variation in the presence and absence of genes, which probably extends to the level of individual pathways. This variation may be a consequence of the significant evolutionary role played by horizontal gene transfer, but might also be explained by the loss of genes through mutation. A challenge is to understand why there would be variation in gene presence within pathways if they confer a benefit only when complete.

Results

Here, we develop a mathematical model to study how variation in pathway content is produced by horizontal transfer, gene loss and partial exposure of a population to a novel environment.

Conclusions

We discuss the possibility that variation in gene presence acts as cryptic genetic variation on which selection acts when the appropriate environment occurs. We find that a high level of variation in gene presence can be readily explained by decay of the pathway through mutation when there is no longer exposure to the selective environment, or when selection becomes too weak to maintain the genes. In the context of pathway variation the role of horizontal gene transfer is probably the initial introduction of a complete novel pathway rather than in building up the variation in a genome without the pathway.