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

Potential links between human bloodstream infection by Salmonella enterica serovar Typhimurium and international transmission to Colombia

Salmonella enterica serovar Typhimurium is a prevalent food-borne pathogen that is usually associated with gastroenteritis infection. S. Typhimurium is also a major cause of bloodstream infections in sub-Saharan Africa, and is responsible for invasive non-typhoidal Salmonella (iNTS) disease. The pathogen also causes bloodstream infection in Colombia, but there has been a lack of information about the S. Typhimurium isolates that were responsible. Here, we investigated the genomic characteristics of 270 S. Typhimurium isolates from bloodstream infection patients in Colombia, collected between 1997 and 2017. We used whole-genome sequencing to analyse multidrug-resistant (MDR) profiles, plasmid distribution, and to define phylogenetic relationships. The study identified the distinct sequence types and phylogenetic clusters of S. Typhimurium prevalent in Colombia. The majority of isolates (90.8%) were ST19, which is distinct from the iNTS-associated S. Typhimurium in sub-Saharan Africa (ST313). The two prominent clusters of MDR S. Typhimurium were either DT104 or closely related to the LT2 reference strain. We used a phylogenetic approach to associate the Colombian clusters with global S. Typhimurium isolates from public databases. By putting the Colombian S. Typhimurium isolates in the context of the global spread of DT104, ST313 and LT2-related variants, we found that the Colombian clusters were introduced to the country via multiple independent events that were consistent with international transmission. We suggest that the acquisition of quinolone and chloramphenicol resistance by the Colombian S. Typhimurium isolates was driven by horizontal gene transfer. Three ST313 isolates that caused bloodstream infection in Colombia were identified. These ST313 isolates were related to the Malawian ST313 lineage 3 & UK ST313, and shared a similarly high invasiveness index. To our knowledge, this is the first report of ST313 in Colombia.

A genomic appraisal of invasive Salmonella Typhimurium and associated antibiotic resistance in sub-Saharan Africa

Invasive non-typhoidal Salmonella (iNTS) disease manifesting as bloodstream infection with high mortality is responsible for a huge public health burden in sub-Saharan Africa. Salmonella enterica serovar Typhimurium (S. Typhimurium) is the main cause of iNTS disease in Africa. By analysing whole genome sequence data from 1303 S. Typhimurium isolates originating from 19 African countries and isolated between 1979 and 2017, here we show a thorough scaled appraisal of the population structure of iNTS disease caused by S. Typhimurium across many of Africa's most impacted countries. At least six invasive S. Typhimurium clades have already emerged, with ST313 lineage 2 or ST313-L2 driving the current pandemic. ST313-L2 likely emerged in the Democratic Republic of Congo around 1980 and further spread in the mid 1990s. We observed plasmid-borne as well as chromosomally encoded fluoroquinolone resistance underlying emergences of extensive-drug and pan-drug resistance. Our work provides an overview of the evolution of invasive S. Typhimurium disease, and can be exploited to target control measures.

Predicting gene expression divergence between single-copy orthologs in two species

Predicting gene expression divergence is integral to understanding the emergence of new biological functions and associated traits. Whereas several sophisticated methods have been developed for this task, their applications are either limited to duplicate genes or require expression data from more than two species. Thus, here we present PiXi, the first machine learning framework for predicting gene expression divergence between single-copy orthologs in two species. PiXi models gene expression evolution as an Ornstein-Uhlenbeck process, and overlays this model with multi-layer neural network, random forest, and support vector machine architectures for making predictions. It outputs the predicted class "conserved" or "diverged" for each pair of orthologs, as well as their predicted expression optima in the two species. We show that PiXi has high power and accuracy in predicting gene expression divergence between single-copy orthologs, as well as high accuracy and precision in estimating their expression optima in the two species, across a wide range of evolutionary scenarios, with the globally best performance achieved by a multi-layer neural network. Moreover, application of our best performing PiXi predictor to empirical gene expression data from single-copy orthologs residing at different loci in two species of Drosophila reveals that approximately 23% underwent expression divergence after positional relocation. Further analysis shows that several of these "diverged" genes are involved in the electron transport chain of the mitochondrial membrane, suggesting that new chromatin environments may impact energy production in Drosophila. Thus, by providing a toolkit for predicting gene expression divergence between single-copy orthologs in two species, PiXi can shed light on the origins of novel phenotypes across diverse biological processes and study systems.

Species-Scale Genomic Analysis of Staphylococcus aureus Genes Influencing Phage Host Range and Their Relationships to Virulence and Antibiotic Resistance Genes

Phage therapy has been proposed as a possible alternative treatment for infections caused by the ubiquitous bacterial pathogen Staphylococcus aureus. However, successful therapy requires understanding the genetic basis of host range—the subset of strains in a species that could be killed by a particular phage. We searched diverse sets of S. aureus public genome sequences against a database of genes suggested from prior studies to influence host range to look for patterns of variation across the species. We found that genes encoding biosynthesis of molecules that were targets of S. aureus phage adsorption to the outer surface of the cell were the most conserved in the pangenome. Putative phage resistance genes that were core components of the pangenome genes had similar nucleotide diversity, ratio of nonsynonymous to synonymous substitutions, and functionality (measured by delta-bitscore) to other core genes. However, phage resistance genes that were not part of the core genome were significantly less consistent with the core genome phylogeny than all noncore genes in this set, suggesting more frequent movement between strains by horizontal gene transfer. Only superinfection immunity genes encoded by temperate phages inserted in the genome correlated with experimentally determined temperate phage resistance. Taken together, these results suggested that, while phage adsorption genes are heavily conserved in the S. aureus species, HGT may play a significant role in strain-specific evolution of host range patterns.

IMPORTANCEStaphylococcus aureus is a widespread, hospital- and community-acquired pathogen that is commonly antibiotic resistant. It causes diverse diseases affecting both the skin and internal organs. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent, such as phage therapy, in which viruses specific to infecting bacteria clear infection. S. aureus phage host range not only determines whether phage therapy will be successful by killing bacteria but also horizontal gene transfer through transduction of host genetic material by phages. In this work, we comprehensively reviewed existing literature to build a list of S. aureus phage resistance genes and searched our database of almost 43,000 S. aureus genomes for these genes to understand their patterns of evolution, finding that prophages’ superinfection immunity correlates best with phage resistance and HGT. These findings improved our understanding of the relationship between known phage resistance genes and phage host range in the species.

Duo: A Signature Based Method to Batch-Analyze Functional Similarities of Proteins

With the rapid advancement of sequencing technology, handling of large sequencing data to analyze for protein coding capacity and functionality of predicted proteins has become an urgent demand. There is a lack of simple and effective tools to functionally annotate large number of unknown proteins in a personalized and customized workflow. To address this, we developed Duo, which batch-analyze functional similarities of predicted proteins. Duo can screen query proteins with specific characteristics based on highly flexible and customizable reference inputs from the user. In the current study, Duo was applied to screen for virulence associated proteins in the genome-sequence of Salmonella Typhimurium. Based on the analysis, recommendation for choice of Seed_database in order to get a reasonable number of predicted proteins for further analysis, and recommendation for preparing a Reference_proteins set for Duo was given. Delta-bitscore analysis was shown to be useful tool to focus the follow-up on predicted proteins. A successful screen for virulence proteins in the bacterial genome-sequence was further performed in a selection of 32 pathogenic bacteria, documenting the ability of Duo to work on a broad collection of bacteria. We anticipate that Duo will be a useful auxiliary tool for personalized and customized protein function research in the future.

Genetic Variation and Preliminary Indications of Divergent Niche Adaptation in Cryptic Clade II of Escherichia

The evolution, habitat, and lifestyle of the cryptic clade II of Escherichia, which were first recovered at low frequency from non-human hosts and later from external environments, were poorly understood. Here, the genomes of selected strains were analyzed for preliminary indications of ecological differentiation within their population. We adopted the delta bitscore metrics to detect functional divergence of their orthologous genes and trained a random forest classifier to differentiate the genomes according to habitats (gastrointestinal vs external environment). Model was built with inclusion of other Escherichia genomes previously demonstrated to have exhibited genomic traits of adaptation to one of the habitats. Overall, gene degradation was more prominent in the gastrointestinal strains. The trained model correctly classified the genomes, identifying a set of predictor genes that were informative of habitat association. Functional divergence in many of these genes were reflective of ecological divergence. Accuracy of the trained model was confirmed by its correct prediction of the habitats of an independent set of strains with known habitat association. In summary, the cryptic clade II of Escherichia displayed genomic signatures that are consistent with divergent adaptation to gastrointestinal and external environments.

Evolution of Salmonella enterica serotype Typhimurium driven by anthropogenic selection and niche adaptation

Salmonella enterica serotype Typhimurium (S. Typhimurium) is a leading cause of gastroenteritis and bacteraemia worldwide, and a model organism for the study of host-pathogen interactions. Two S. Typhimurium strains (SL1344 and ATCC14028) are widely used to study host-pathogen interactions, yet genotypic variation results in strains with diverse host range, pathogenicity and risk to food safety. The population structure of diverse strains of S. Typhimurium revealed a major phylogroup of predominantly sequence type 19 (ST19) and a minor phylogroup of ST36. The major phylogroup had a population structure with two high order clades (α and β) and multiple subclades on extended internal branches, that exhibited distinct signatures of host adaptation and anthropogenic selection. Clade α contained a number of subclades composed of strains from well characterized epidemics in domesticated animals, while clade β contained multiple subclades associated with wild avian species. The contrasting epidemiology of strains in clade α and β was reflected by the distinct distribution of antimicrobial resistance (AMR) genes, accumulation of hypothetically disrupted coding sequences (HDCS), and signatures of functional diversification. These observations were consistent with elevated anthropogenic selection of clade α lineages from adaptation to circulation in populations of domesticated livestock, and the predisposition of clade β lineages to undergo adaptation to an invasive lifestyle by a process of convergent evolution with of host adapted Salmonella serotypes. Gene flux was predominantly driven by acquisition and recombination of prophage and associated cargo genes, with only occasional loss of these elements. The acquisition of large chromosomally-encoded genetic islands was limited, but notably, a feature of two recent pandemic clones (DT104 and monophasic S. Typhimurium ST34) of clade α (SGI-1 and SGI-4).

Genomic, transcriptomic, and proteomic insights into the symbiosis of deep-sea tubeworm holobionts

Deep-sea hydrothermal vents and methane seeps are often densely populated by animals that host chemosynthetic symbiotic bacteria, but the molecular mechanisms of such host-symbiont relationship remain largely unclear. We characterized the symbiont genome of the seep-living siboglinid Paraescarpia echinospica and compared seven siboglinid-symbiont genomes. Our comparative analyses indicate that seep-living siboglinid endosymbionts have more virulence traits for establishing infections and modulating host-bacterium interaction than the vent-dwelling species, and have a high potential to resist environmental hazards. Metatranscriptome and metaproteome analyses of the Paraescarpia holobiont reveal that the symbiont is highly versatile in its energy use and efficient in carbon fixation. There is close cooperation within the holobiont in production and supply of nutrients, and the symbiont may be able to obtain nutrients from host cells using virulence factors. Moreover, the symbiont is speculated to have evolved strategies to mediate host protective immunity, resulting in weak expression of host innate immunity genes in the trophosome. Overall, our results reveal the interdependence of the tubeworm holobiont through mutual nutrient supply, a pathogen-type regulatory mechanism, and host-symbiont cooperation in energy utilization and nutrient production, which is a key adaptation allowing the tubeworm to thrive in deep-sea chemosynthetic environments.

An African Salmonella Typhimurium ST313 sublineage with extensive drug-resistance and signatures of host adaptation

Bloodstream infections by Salmonella enterica serovar Typhimurium constitute a major health burden in sub-Saharan Africa (SSA). These invasive non-typhoidal (iNTS) infections are dominated by isolates of the antibiotic resistance-associated sequence type (ST) 313. Here, we report emergence of ST313 sublineage II.1 in the Democratic Republic of the Congo. Sublineage II.1 exhibits extensive drug resistance, involving a combination of multidrug resistance, extended spectrum β-lactamase production and azithromycin resistance. ST313 lineage II.1 isolates harbour an IncHI2 plasmid we name pSTm-ST313-II.1, with one isolate also exhibiting decreased ciprofloxacin susceptibility. Whole genome sequencing reveals that ST313 II.1 isolates have accumulated genetic signatures potentially associated with altered pathogenicity and host adaptation, related to changes observed in biofilm formation and metabolic capacity. Sublineage II.1 emerged at the beginning of the 21st century and is involved in on-going outbreaks. Our data provide evidence of further evolution within the ST313 clade associated with iNTS in SSA.

Convergent regulatory evolution and loss of flight in paleognathous birds

A core question in evolutionary biology is whether convergent phenotypic evolution is driven by convergent molecular changes in proteins or regulatory regions. We combined phylogenomic, developmental, and epigenomic analysis of 11 new genomes of paleognathous birds, including an extinct moa, to show that convergent evolution of regulatory regions, more so than protein-coding genes, is prevalent among developmental pathways associated with independent losses of flight. A Bayesian analysis of 284,001 conserved noncoding elements, 60,665 of which are corroborated as enhancers by open chromatin states during development, identified 2355 independent accelerations along lineages of flightless paleognaths, with functional consequences for driving gene expression in the developing forelimb. Our results suggest that the genomic landscape associated with morphological convergence in ratites has a substantial shared regulatory component.

Genome-Based Prediction of Bacterial Antibiotic Resistance

Clinical microbiology has long relied on growing bacteria in culture to determine antimicrobial susceptibility profiles, but the use of whole-genome sequencing for antibiotic susceptibility testing (WGS-AST) is now a powerful alternative. This review discusses the technologies that made this possible and presents results from recent studies to predict resistance based on genome sequences. We examine differences between calling antibiotic resistance profiles by the simple presence or absence of previously known genes and single-nucleotide polymorphisms (SNPs) against approaches that deploy machine learning and statistical models. Often, the limitations to genome-based prediction arise from limitations of accuracy of culture-based AST in addition to an incomplete knowledge of the genetic basis of resistance. However, we need to maintain phenotypic testing even as genome-based prediction becomes more widespread to ensure that the results do not diverge over time. We argue that standardization of WGS-AST by challenge with consistently phenotyped strain sets of defined genetic diversity is necessary to compare the efficacy of methods of prediction of antibiotic resistance based on genome sequences.

Genomic correlates of extraintestinal infection are linked with changes in cell morphology in Campylobacter jejuni

Campylobacter jejuni is the most common cause of bacterial diarrheal disease in the world. Clinical outcomes of infection can range from asymptomatic infection to life-threatening extraintestinal infections. This variability in outcomes for infected patients has raised questions as to whether genetic differences between C. jejuni isolates contribute to their likelihood of causing severe disease. In this study, we compare the genomes of ten C. jejuni isolates that were implicated in extraintestinal infections with reference gastrointestinal isolates, in order to identify unusual patterns of sequence variation associated with infection outcome. We identified a collection of genes that display a higher burden of uncommon mutations in invasive isolates compared with gastrointestinal close relatives, including some that have been previously linked to virulence and invasiveness in C. jejuni. Among the top genes identified were mreB and pgp1, which are both involved in determining cell shape. Electron microscopy confirmed morphological differences in isolates carrying unusual sequence variants of these genes, indicating a possible relationship between extraintestinal infection and changes in cell morphology.

Genome Variation and Molecular Epidemiology of Salmonella enterica Serovar Typhimurium Pathovariants

Salmonella enterica serovar Typhimurium is one of approximately 2,500 distinct serovars of the genus Salmonella but is exceptional in its wide distribution in the environment, livestock, and wild animals. S Typhimurium causes a large proportion of nontyphoidal Salmonella (NTS) infections, accounting for a quarter of infections, second only to S. enterica serovar Enteritidis in incidence. S Typhimurium was once considered the archetypal broad-host-range Salmonella serovar due to its wide distribution in livestock and wild animals, and much of what we know of the interaction of Salmonella with the host comes from research using a small number of laboratory strains of the serovar (LT2, SL1344, and ATCC 14028). But it has become clear that these strains do not reflect the genotypic or phenotypic diversity of S Typhimurium. Here, we review the epidemiological record of S Typhimurium and studies of the host-pathogen interactions of diverse strains of S Typhimurium. We present the concept of distinct pathovariants of S Typhimurium that exhibit diversity of host range, distribution in the environment, pathogenicity, and risk to food safety. We review recent evidence from whole-genome sequencing that has revealed the extent of genomic diversity of S Typhimurium pathovariants, the genomic basis of differences in the level of risk to human and animal health, and the molecular epidemiology of prominent strains. An improved understanding of the impact of genome variation of bacterial pathogens on pathogen-host and pathogen-environment interactions has the potential to improve quantitative risk assessment and reveal how new pathogens evolve.

Evolution of Salmonella within Hosts

Within-host evolution has resulted in thousands of variants of Salmonella that exhibit remarkable diversity in host range and disease outcome, from broad host range to exquisite host restriction, causing gastroenteritis to disseminated disease such as typhoid fever. Within-host evolution is a continuing process driven by genomic variation that occurs during each infection, potentiating adaptation to a new niche resulting from changes in animal husbandry, the use of antimicrobials, and emergence of immune compromised populations. We discuss key advances in our understanding of the evolution of Salmonella within the host, inferred from (i) the process of host adaptation of Salmonella pathovars in the past, and (ii) direct observation of the generation of variation and selection of beneficial traits during single infections.

Salmonella is a bacterial pathogen with remarkable diversity in its host range and pathogenicity due to past within-host evolution in vertebrate species that modified ancestral mechanisms of pathogenesis.

Variation arising during infection includes point mutations, new genes acquired through horizontal gene transfer (HGT), deletions, and genomic rearrangements.

Beneficial mutations increase in frequency within the host and, if they retain the ability to be transmitted to subsequent hosts, may become fixed in the population.

Whole-genome sequencing of sequential isolates from clinical infections reveals within-host HGT and point mutations that impact therapy and clinical management.

HGT is the primary mechanism for evolution in prokaryotes and is synergised by complex networks of transfer involving the microbiome.

Within-host evolution of Salmonella, resulting in new pathovars, can proceed in the absence of HGT.

Machine learning identifies signatures of host adaptation in the bacterial pathogen Salmonella enterica

Emerging pathogens are a major threat to public health, however understanding how pathogens adapt to new niches remains a challenge. New methods are urgently required to provide functional insights into pathogens from the massive genomic data sets now being generated from routine pathogen surveillance for epidemiological purposes. Here, we measure the burden of atypical mutations in protein coding genes across independently evolved Salmonella enterica lineages, and use these as input to train a random forest classifier to identify strains associated with extraintestinal disease. Members of the species fall along a continuum, from pathovars which cause gastrointestinal infection and low mortality, associated with a broad host-range, to those that cause invasive infection and high mortality, associated with a narrowed host range. Our random forest classifier learned to perfectly discriminate long-established gastrointestinal and invasive serovars of Salmonella. Additionally, it was able to discriminate recently emerged Salmonella Enteritidis and Typhimurium lineages associated with invasive disease in immunocompromised populations in sub-Saharan Africa, and within-host adaptation to invasive infection. We dissect the architecture of the model to identify the genes that were most informative of phenotype, revealing a common theme of degradation of metabolic pathways in extraintestinal lineages. This approach accurately identifies patterns of gene degradation and diversifying selection specific to invasive serovars that have been captured by more labour-intensive investigations, but can be readily scaled to larger analyses.

Genomic, Transcriptomic, and Phenotypic Analyses of Neisseria meningitidis Isolates from Disease Patients and Their Household Contacts

Neisseria meningitidis causes meningococcal disease but is frequently carried in the throats of healthy individuals; the factors that determine whether invasive disease develops are not completely understood. We carried out detailed studies of isolates, collected from patients and their household contacts, to identify differences between commensal throat isolates and those that caused invasive disease. Though isolates were identical by laboratory typing methods, we uncovered many differences in their genomes, in gene expression, and in their interactions with host cells. In particular, we found that several carriage isolates had lost their type IV pili, a surprising finding since pili are often described as essential for colonization. However, loss of type IV pili correlated with reduced secretion of a proinflammatory cytokine, TNF-α, when meningococci were cocultured with human bronchial epithelial cells; hence, the loss of pili could provide an advantage to meningococci, by resulting in a dampened localized host immune response.

Neisseria meningitidis (meningococcus) can cause meningococcal disease, a rapidly progressing and often fatal disease that can occur in previously healthy children. Meningococci are found in healthy carriers, where they reside in the nasopharynx as commensals. While carriage is relatively common, invasive disease, associated with hypervirulent strains, is a comparatively rare event. The basis of increased virulence in some strains is not well understood. New Zealand suffered a protracted meningococcal disease epidemic, from 1991 to 2008. During this time, a household carriage study was carried out in Auckland: household contacts of index meningococcal disease patients were swabbed for isolation of carriage strains. In many households, healthy carriers harbored strains identical, as determined by laboratory typing, to the ones infecting the associated patient. We carried out more-detailed analyses of carriage and disease isolates from a select number of households. We found that isolates, although indistinguishable by laboratory typing methods and likely closely related, had many differences. We identified multiple genome variants and transcriptional differences between isolates. These studies enabled the identification of two new phase-variable genes. We also found that several carriage strains had lost their type IV pili and that this loss correlated with reduced tumor necrosis factor alpha (TNF-α) expression when cultured with epithelial cells. While nonpiliated meningococcal isolates have been previously found in carriage strains, this is the first evidence of an association between type IV pili from meningococci and a proinflammatory epithelial response. We also identified potentially important metabolic differences between carriage and disease isolates, including the sulfate assimilation pathway.

IMPORTANCENeisseria meningitidis causes meningococcal disease but is frequently carried in the throats of healthy individuals; the factors that determine whether invasive disease develops are not completely understood. We carried out detailed studies of isolates, collected from patients and their household contacts, to identify differences between commensal throat isolates and those that caused invasive disease. Though isolates were identical by laboratory typing methods, we uncovered many differences in their genomes, in gene expression, and in their interactions with host cells. In particular, we found that several carriage isolates had lost their type IV pili, a surprising finding since pili are often described as essential for colonization. However, loss of type IV pili correlated with reduced secretion of a proinflammatory cytokine, TNF-α, when meningococci were cocultured with human bronchial epithelial cells; hence, the loss of pili could provide an advantage to meningococci, by resulting in a dampened localized host immune response.