Comparative genomics of 274 Vibrio cholerae genomes reveals mobile functions structuring three niche dimensions

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

The Promise and Pitfalls of Prophages

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

Beyond to the Stable: Role of the Insertion Sequences as Epidemiological Descriptors in Corynebacterium striatum

In recent years, epidemiological studies of infectious agents have focused mainly on the pathogen and stable components of its genome. The use of these stable components makes it possible to know the evolutionary or epidemiological relationships of the isolates of a particular pathogen. Under this approach, focused on the pathogen, the identification of resistance genes is a complementary stage of a bacterial characterization process or an appendix of its epidemiological characterization, neglecting its genetic components’ acquisition or dispersal mechanisms. Today we know that a large part of antibiotic resistance is associated with mobile elements. Corynebacterium striatum, a bacterium from the normal skin microbiota, is also an opportunistic pathogen. In recent years, reports of infections and nosocomial outbreaks caused by antimicrobial multidrug-resistant C. striatum strains have been increasing worldwide. Despite the different existing mobile genomic elements, there is evidence that acquired resistance genes are coupled to insertion sequences in C. striatum. This perspective article reviews the insertion sequences linked to resistance genes, their relationship to evolutionary lineages, epidemiological characteristics, and the niches the strains inhabit. Finally, we evaluate the potential of the insertion sequences for their application as a descriptor of epidemiological scenarios, allowing us to anticipate the emergence of multidrug-resistant lineages.

Genomic insights into Vibrio cholerae O1 responsible for cholera epidemics in Tanzania between 1993 and 2017

BACKGROUND

Tanzania is one of seven countries with the highest disease burden caused by cholera in Africa. We studied the evolution of Vibrio cholerae O1 isolated in Tanzania during the past three decades.

METHODOLOGY/PRINCIPAL FINDINGS

Genome-wide analysis was performed to characterize V. cholerae O1 responsible for the Tanzanian 2015-2017 outbreak along with strains causing outbreaks in the country for the past three decades. The genomes were further analyzed in a global context of 590 strains of the seventh cholera pandemic (7PET), as well as environmental isolates from Lake Victoria. All Tanzanian cholera outbreaks were caused by the 7PET lineage. The T5 sub-lineage (ctxB3) dominated outbreaks until 1997, followed by the T10 atypical El Tor (ctxB1) up to 2015, which were replaced by the T13 atypical El Tor of the current third wave (ctxB7) causing most cholera outbreaks until 2017 with T13 being phylogenetically related to strains from East African countries, Yemen and Lake Victoria. The strains were less drug resistant with approximate 10-kb deletions found in the SXT element, which encodes resistance to sulfamethoxazole and trimethoprim. Nucleotide deletions were observed in the CTX prophage of some strains, which warrants further virulence studies. Outbreak strains share 90% of core genes with V. cholerae O1 from Lake Victoria with as low as three SNPs difference and a significantly similar accessory genome, composed of genomic islands namely the CTX prophage, Vibrio Pathogenicity Islands; toxin co-regulated pilus biosynthesis proteins and the SXT-ICE element.

CONCLUSION/SIGNIFICANCE

Characterization of V. cholerae O1 from Tanzania reveals genetic diversity of the 7PET lineage composed of T5, T10 and T13 sub-lineages with introductions of new sequence types from neighboring countries. The presence of these sub-lineages in environmental isolates suggests that the African Great Lakes may serve as aquatic reservoirs for survival of V. cholerae O1 favoring continuous human exposure.

Revisiting the Global Epidemiology of Cholera in Conjuction With the Genomics of Vibrio cholerae

Toxigenic Vibrio cholerae is responsible for 1.4 to 4.3 million cases with about 21,000-143,000 deaths per year. Dominance of O1 and O139 serogroups, classical and El tor biotypes, alterations in CTX phages and the pathogenicity Islands are some of the major features of V. cholerae isolates that are responsible for cholera epidemics. Whole-genome sequencing (WGS) based analyses of single-nucleotide polymorphisms (SNPs) and other infrequent genetic variants provide a robust phylogenetic framework. Recent studies on the global transmission of pandemic V. cholerae O1 strains have shown the existence of eight different phyletic lineages. In these, the classical and El Tor biotype strains were separated as two distinctly evolved lineages. The frequency of SNP accumulation and the temporal and geographical distribution supports the perception that the seventh cholera pandemic (7CP) has spread from the Bay of Bengal region in three independent but overlapping waves. The 2010 Haitian outbreak shared a common ancestor with South-Asian wave-3 strains. In West Africa and East/Southern Africa, cholera epidemics are caused by single expanded lineage, which has been introduced several times since 1970. The Latin American epidemics that occurred in 1991 and 2010 were the result of introductions of two 7CP sublineages. Sublineages representing wave-3 have caused huge outbreaks in Haiti and Yemen. The Ogawa-Inaba serotype switchover in several cholera epidemics are believed to be due to the involvement of certain selection mechanism(s) rather than due to random events. V. cholerae O139 serogroup is phylogenetically related to the 7CP El Tor, and almost all these isolates belonged to the multilocus sequence type-69. Additional phenotypic and genotypic information have been generated to understand the pathogenicity of classical and El Tor vibrios. Presence of integrative conjugative elements (ICE) with antibiotic resistance gene cassettes, clustered regularly interspaced short palindromic repeats-associated protein system and ctxAB promoter based ToxRS expression of cholera toxin (CT) separates classical and El Tor biotypes. With the availability of WGS information, several important applications including, molecular typing, antimicrobial resistance, new diagnostics, and vaccination strategies could be generated.

Antimicrobial resistance of Vibrio cholerae from sub-Saharan Africa: A systematic review

Background

The World Health Assembly adopted the Global Action Plan on Antimicrobial Resistance, which includes improving the knowledge base through surveillance and research. Noteworthily, the World Health Organization has advocated a Global Antimicrobial Resistance Surveillance System to address the plan’s surveillance objective, with most African countries enrolling in or after 2017.

Aim

The aim of this article was to review prior data on antimicrobial resistance of Vibrio cholerae from sub-Saharan Africa with a view for future control and intervention strategies.

Methods

We used the Preferred Reporting Items for Systematic Review and Meta-Analysis (or ‘PRISMA’) guidelines to search the PubMed and African Journals Online databases, as well as additional articles provided by the Nigeria Centre for Disease Control, for articles reporting on the antibiotic susceptibility of V. cholerae between January 2000 and December 2017.

Results

We identified 340 publications, of which only 25 (reporting from 16 countries within the sub-Saharan African region) were eligible. The majority (20; 80.0%) of the cholera toxigenic V. cholerae isolates were of the serogroup O1 of the El Tor biotype with Ogawa and Inaba serotypes predominating. Resistance was predominantly documented to trimethoprim-sulphamethoxazole (50% of the studies), ampicillin (43.3% of the studies), chloramphenicol (43.3% of the studies) and streptomycin (30% of the studies). Resistance mechanisms were reported in 40% of the studies.

Conclusion

Our results demonstrate a documented antimicrobial resistance of V. cholerae to multiple antibiotic classes, including cell wall active agents and antimetabolites with evidence of phenotypic/genotypic resistance to fluoroquinolones.

Comparative genomics for non-O1/O139 Vibrio cholerae isolates recovered from the Yangtze River Estuary versus V. cholerae representative isolates from serogroup O1

Vibriocholerae, which is autochthonous to estuaries worldwide, can cause human cholera that is still pandemic in developing countries. A number of V. cholerae isolates of clinical and environmental origin worldwide have been subjected to genome sequencing to address their phylogenesis and bacterial pathogenesis, however, little genome information is available for V. cholerae isolates derived from estuaries, particularly in China. In this study, we determined the complete genome sequence of V. cholerae CHN108B (non-O1/O139 serogroup) isolated from the Yangtze River Estuary, China and performed comparative genome analysis between CHN108B and other eight representative V. cholerae isolates. The 4,168,545-bp V. cholerae CHN108B genome (47.2% G+C) consists of two circular chromosomes with 3,691 predicted protein-encoding genes. It has 110 strain-specific genes, the highest number among the eight representative V. cholerae whole genomes from serogroup O1: there are seven clinical isolates linked to cholera pandemics (1937-2010) and one environmental isolate from Brazil. Various mobile genetic elements (such as insertion sequences, prophages, integrative and conjugative elements, and super-integrons) were identified in the nine V. cholerae genomes of clinical and environmental origin, indicating that the bacterium undergoes extensive genetic recombination via lateral gene transfer. Comparative genomics also revealed different virulence and antimicrobial resistance gene patterns among the V. cholerae isolates, suggesting some potential virulence factors and the rising development of resistance among pathogenic V. cholerae. Additionally, draft genome sequences of multiple V. cholerae isolates recovered from the Yangtze River Estuary were also determined, and comparative genomics revealed many genes involved in specific metabolism pathways, which are likely shaped by the unique estuary environment. These results provide additional evidence of V. cholerae genome plasticity and will facilitate better understanding of the genome evolution and pathogenesis of this severe water-borne pathogen worldwide.

Analysis of Vibrio seventh pandemic island II and novel genomic islands in relation to attachment sequences among a wide variety of Vibrio cholerae strains

Vibrio cholerae O1 El Tor, the pathogen responsible for the current cholera pandemic, became pathogenic by acquiring virulent factors including Vibrio seventh pandemic islands (VSP)‐I and −II. Diversity of VSP‐II is well recognized; however, studies addressing attachment sequence left (attL) sequences of VSP‐II are few. In this report, a wide variety of V. cholerae strains were analyzed for the structure and distribution of VSP‐II in relation to their attachment sequences. Of 188 V. cholerae strains analyzed, 81% (153/188) strains carried VSP‐II; of these, typical VSP‐II, and a short variant was found in 36% (55/153), and 63% (96/153), respectively. A novel VSP‐II was found in two V. cholerae non‐O1/non‐O139 strains. In addition to the typical 14‐bp attL, six new attL‐like sequences were identified. The 14‐bp attL was associated with VSP‐II in 91% (139/153), whereas the remaining six types were found in 9.2% (14/153) of V. cholerae strains. Of note, six distinct types of the attL‐like sequence were found in the seventh pandemic wave 1 strains; however, only one or two types were found in the wave 2 or 3 strains. Interestingly, 86% (24/28) of V. cholerae seventh pandemic strains harboring a 13‐bp attL‐like sequence were devoid of VSP‐II. Six novel genomic islands using two unique insertion sites to those of VSP‐II were identified in 11 V. cholerae strains in this study. Four of those shared similar gene clusters with VSP‐II, except integrase gene.

Computational pan-genomics: status, promises and challenges

Many disciplines, from human genetics and oncology to plant breeding, microbiology and virology, commonly face the challenge of analyzing rapidly increasing numbers of genomes. In case of Homo sapiens, the number of sequenced genomes will approach hundreds of thousands in the next few years. Simply scaling up established bioinformatics pipelines will not be sufficient for leveraging the full potential of such rich genomic data sets. Instead, novel, qualitatively different computational methods and paradigms are needed. We will witness the rapid extension of computational pan-genomics, a new sub-area of research in computational biology. In this article, we generalize existing definitions and understand a pan-genome as any collection of genomic sequences to be analyzed jointly or to be used as a reference. We examine already available approaches to construct and use pan-genomes, discuss the potential benefits of future technologies and methodologies and review open challenges from the vantage point of the above-mentioned biological disciplines. As a prominent example for a computational paradigm shift, we particularly highlight the transition from the representation of reference genomes as strings to representations as graphs. We outline how this and other challenges from different application domains translate into common computational problems, point out relevant bioinformatics techniques and identify open problems in computer science. With this review, we aim to increase awareness that a joint approach to computational pan-genomics can help address many of the problems currently faced in various domains.

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.

A genomic investigation of ecological differentiation between free-living and Drosophila-associated bacteria

Various bacterial taxa have been identified both in association with animals and in the external environment, but the extent to which related bacteria from the two habitat types are ecologically and evolutionarily distinct is largely unknown. This study investigated the scale and pattern of genetic differentiation between bacteria of the family Acetobacteraceae isolated from the guts of Drosophila fruit flies, plant material and industrial fermentations. Genome-scale analysis of the phylogenetic relationships and predicted functions was conducted on 44 Acetobacteraceae isolates, including newly sequenced genomes from 18 isolates from wild and laboratory Drosophila. Isolates from the external environment and Drosophila could not be assigned to distinct phylogenetic groups, nor are their genomes enriched for any different sets of genes or category of predicted gene functions. In contrast, analysis of bacteria from laboratory Drosophila showed they were genetically distinct in their universal capacity to degrade uric acid (a major nitrogenous waste product of Drosophila) and absence of flagellar motility, while these traits vary among wild Drosophila isolates. Analysis of the competitive fitness of Acetobacter discordant for these traits revealed a significant fitness deficit for bacteria that cannot degrade uric acid in culture with Drosophila. We propose that, for wild populations, frequent cycling of Acetobacter between Drosophila and the external environment prevents genetic differentiation by maintaining selection for traits adaptive in both the gut and external habitats. However, laboratory isolates bear the signs of adaptation to persistent association with the Drosophila host under tightly defined environmental conditions.

Genetic diversity of environmental Vibrio cholerae O1 strains isolated in Northern Vietnam

Cholera epidemics have been recorded periodically in Vietnam during the seventh cholera pandemic. Since cholera is a water-borne disease, systematic monitoring of environmental waters for Vibrio cholerae presence is important for predicting and preventing cholera epidemics. We conducted monitoring, isolation, and genetic characterization of V. cholerae strains in Nam Dinh province of Northern Vietnam from Jul 2013 to Feb 2015. In this study, four V. cholerae O1 strains were detected and isolated from 110 analyzed water samples (3.6%); however, none of them carried the cholera toxin gene, ctxA, in their genomes. Whole genome sequencing and phylogenetic analysis revealed that the four O1 isolates were separated into two independent clusters, and one of them diverged from a common ancestor with pandemic strains. The analysis of pathogenicity islands (CTX prophage, VPI-I, VPI-II, VSP-I, and VSP-II) indicated that one strain (VNND_2014Jun_6SS) harbored an unknown prophage-like sequence with high homology to vibriophage KSF-1 phi and VCY phi, identified from Bangladesh and the USA, respectively, while the other three strains carried tcpA gene with a distinct sequence demonstrating a separate clonal lineage. These results suggest that the aquatic environment can harbor highly divergent V. cholera strains and serve as a reservoir for multiple V. cholerae virulence-associated genes which may be exchanged via mobile genetic elements. Therefore, continuous monitoring and genetic characterization of V. cholerae strains in the environment should contribute to the early detection of the sources of infection and prevention of cholera outbreaks as well as to understanding the natural ecology and evolution of V. cholerae.

Insights of Phage-Host Interaction in Hypersaline Ecosystem through Metagenomics Analyses

Bacteriophages, as the most abundant biological entities on Earth, place significant predation pressure on their hosts. This pressure plays a critical role in the evolution, diversity, and abundance of bacteria. In addition, phages modulate the genetic diversity of prokaryotic communities through the transfer of auxiliary metabolic genes. Various studies have been conducted in diverse ecosystems to understand phage-host interactions and their effects on prokaryote metabolism and community composition. However, hypersaline environments remain among the least studied ecosystems and the interaction between the phages and prokaryotes in these habitats is poorly understood. This study begins to fill this knowledge gap by analyzing bacteriophage-host interactions in the Great Salt Lake, the largest prehistoric hypersaline lake in the Western Hemisphere. Our metagenomics analyses allowed us to comprehensively identify the bacterial and phage communities with Proteobacteria, Firmicutes, and Bacteroidetes as the most dominant bacterial species and Siphoviridae, Myoviridae, and Podoviridae as the most dominant viral families found in the metagenomic sequences. We also characterized interactions between the phage and prokaryotic communities of Great Salt Lake and determined how these interactions possibly influence the community diversity, structure, and biogeochemical cycles. In addition, presence of prophages and their interaction with the prokaryotic host was studied and showed the possibility of prophage induction and subsequent infection of prokaryotic community present in the Great Salt Lake environment under different environmental stress factors. We found that carbon cycle was the most susceptible nutrient cycling pathways to prophage induction in the presence of environmental stresses. This study gives an enhanced snapshot of phage and prokaryote abundance and diversity as well as their interactions in a hypersaline complex ecosystem, which can pave the way for further research studies.

Biotic Interactions Shape the Ecological Distributions of Staphylococcus Species

Many metagenomic sequencing studies have observed the presence of closely related bacterial species or genotypes in the same microbiome. Previous attempts to explain these patterns of microdiversity have focused on the abiotic environment, but few have considered how biotic interactions could drive patterns of microbiome diversity. We dissected the patterns, processes, and mechanisms shaping the ecological distributions of three closely related Staphylococcus species in cheese rind biofilms. Paradoxically, the most abundant species (S. equorum) is the slowest colonizer and weakest competitor based on growth and competition assays in the laboratory. Through in vitro community reconstructions, we determined that biotic interactions with neighboring fungi help resolve this paradox. Species-specific stimulation of the poor competitor by fungi of the genus Scopulariopsis allows S. equorum to dominate communities in vitro as it does in situ. Results of comparative genomic and transcriptomic experiments indicate that iron utilization pathways, including a homolog of the S. aureus staphyloferrin B siderophore operon pathway, are potential molecular mechanisms underlying Staphylococcus-Scopulariopsis interactions. Our integrated approach demonstrates that fungi can structure the ecological distributions of closely related bacterial species, and the data highlight the importance of bacterium-fungus interactions in attempts to design and manipulate microbiomes.

Decades of culture-based studies and more recent metagenomic studies have demonstrated that bacterial species in agriculture, medicine, industry, and nature are unevenly distributed across time and space. The ecological processes and molecular mechanisms that shape these distributions are not well understood because it is challenging to connect in situ patterns of diversity with mechanistic in vitro studies in the laboratory. Using tractable cheese rind biofilms and a focus on coagulase-negative Staphylococcus (CNS) species, we demonstrate that fungi can mediate the ecological distributions of closely related bacterial species. One of the Staphylococcus species studied, S. saprophyticus, is a common cause of urinary tract infections. By identifying processes that control the abundance of undesirable CNS species, cheese producers will have more precise control on the safety and quality of their products. More generally, Staphylococcus species frequently co-occur with fungi in mammalian microbiomes, and similar bacterium-fungus interactions may structure bacterial diversity in these systems.

Next-generation sequencing (NGS) for assessment of microbial water quality: current progress, challenges, and future opportunities

Water quality is an emergent property of a complex system comprised of interacting microbial populations and introduced microbial and chemical contaminants. Studies leveraging next-generation sequencing (NGS) technologies are providing new insights into the ecology of microbially mediated processes that influence fresh water quality such as algal blooms, contaminant biodegradation, and pathogen dissemination. In addition, sequencing methods targeting small subunit (SSU) rRNA hypervariable regions have allowed identification of signature microbial species that serve as bioindicators for sewage contamination in these environments. Beyond amplicon sequencing, metagenomic and metatranscriptomic analyses of microbial communities in fresh water environments reveal the genetic capabilities and interplay of waterborne microorganisms, shedding light on the mechanisms for production and biodegradation of toxins and other contaminants. This review discusses the challenges and benefits of applying NGS-based methods to water quality research and assessment. We will consider the suitability and biases inherent in the application of NGS as a screening tool for assessment of biological risks and discuss the potential and limitations for direct quantitative interpretation of NGS data. Secondly, we will examine case studies from recent literature where NGS based methods have been applied to topics in water quality assessment, including development of bioindicators for sewage pollution and microbial source tracking, characterizing the distribution of toxin and antibiotic resistance genes in water samples, and investigating mechanisms of biodegradation of harmful pollutants that threaten water quality. Finally, we provide a short review of emerging NGS platforms and their potential applications to the next generation of water quality assessment tools.

Microbial taxonomy in the post-genomic era: rebuilding from scratch?

Microbial taxonomy should provide adequate descriptions of bacterial, archaeal, and eukaryotic microbial diversity in ecological, clinical, and industrial environments. Its cornerstone, the prokaryote species has been re-evaluated twice. It is time to revisit polyphasic taxonomy, its principles, and its practice, including its underlying pragmatic species concept. Ultimately, we will be able to realize an old dream of our predecessor taxonomists and build a genomic-based microbial taxonomy, using standardized and automated curation of high-quality complete genome sequences as the new gold standard.