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

Genomic census of invasive nontyphoidal Salmonella infections reveals global and local human-to-human transmission

Extraintestinal infections caused by Enterobacteriaceae represent a global concern, further exacerbated by the growing prevalence of antimicrobial resistance (AMR). Among these, invasive nontyphoidal Salmonella (iNTS) infections have become increasingly challenging to manage, and their global spread remains poorly understood. Here we compiled 1,115 patient records and generated a comprehensive genomic dataset on iNTS. Age and sex emerged as significant risk factors, with Salmonella Enteritidis identified as a major cause. We observed serovar-specific AMR patterns, with notable resistance to fluoroquinolones and third-generation cephalosporins. A global phylogenomic analysis of Enteritidis revealed three distinct clades, highlighting the accumulation of AMR determinants during its international spread. Importantly, our genomic and transmission analyses suggest that iNTS infections may involve human-to-human transmission, with diarrheal patients acting as potential intermediaries, deviating from typical zoonotic pathways. Collectively, our newly generated cohort and iNTS genomic dataset provide a framework for precise local iNTS burden and underscore emerging transmission trends.

Genomic epidemiology and public health implications of zoonotic monophasic Salmonella Typhimurium ST34

Background

Monophasic Salmonella Typhimurium sequence type 34 (mSTM ST34) has emerged as a significant global health threat, but our understanding of its genomic epidemiology and potential public health implications in international and regional contexts remains limited. This study aims to fill this crucial gap by assessing the genomic epidemiology of multidrug resistance (MDR) mSTM ST34, as well as its clinical characteristics and virulence.

Methods

To achieve the objectives of this study, we conducted a comprehensive genomic analysis of mSTM ST34 isolates. We obtained a global dataset comprising 13,844 strains from public databases, along with 339 strains from a regional surveillance collection in Zhejiang Province, China. This dataset aims to provide in-depth insights into antimicrobial resistance, mobile genetic elements, and pathogenicity. Additionally, we meticulously assessed the association between phenotypic profiles and clinical presentations.

Results

Our findings revealed that the prevalence of mSTM ST34 has surpassed that of the previously dominant ST19. In addition, we observed an increase in the detection of the IncQ1 plasmid, which is responsible for disseminating MDR. The prevalence of mSTM ST34 carriage was exceptionally high among children (≤12 years old) and elderly individuals (≥65 years old), with 92.6% of the isolates exhibiting MDR, including resistance to frontline antimicrobials such as third-generation cephalosporins and ciprofloxacin. Additionally, the human mSTM ST34 strain demonstrates a remarkable capacity for biofilm formation, which increases its virulence in animal models and complicates therapeutic interventions.

Conclusions

mSTM ST34 has surpassed the previously dominant ST19, and its ability to transmit across multi-species increases its potential for further human transmission. This study addresses critical gaps in our understanding of mSTM ST34 prevalence, highlighting the importance of whole genome sequencing in surveilling zoonotic pathogens.

A genomic and phenotypic investigation of pigeon-adaptive Salmonella

Salmonella, a significant threat to public safety, inflicts substantial economic losses on the poultry industry. The unique "parental feeding" breeding model of pigeon farms, against the "all-in & all-out" biosecurity strategy, makes them susceptible to Salmonella infections and subsequent outbreaks of pigeon paratyphoid. This study initially studied three pigeon paratyphoid outbreak incidents in Henan, China, in which 53 strains of pigeon-origin Salmonella Typhimurium (STM) were identified. Whole-genome sequencing (WGS) and antimicrobial-resistant profile analysis revealed that the three outbreaks were caused by distinct STM clones (ST128-DT2, ST19-DT99). Global phylogenetic analysis suggested that the United States is a possible origin, indicating a risk of intercontinental transmission via pigeon eggs. Further bacterial virulence and invasion assays, including in vitro and in vivo assays, revealed that pigeon-host-adaptive STM, compared to broad-host-range STM, carried fewer resistance genes, exhibited higher invasion indices and pseudogene levels, displayed a non-rdar (red dry and rough) phenotype, and had strong biofilm formation capability. Additionally, they showed reduced virulence and invasiveness in mice but a pigeon-adaptive feature in cogent models. The collective results support the host adaptation for pigeons among DT2 and DT99 phage-type isolates.

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.

Using GWAS and Machine Learning to Identify and Predict Genetic Variants Associated with Foodborne Bacteria Phenotypic Traits

One of the main challenges in food microbiology is to prevent the risk of outbreaks by avoiding the distribution of food contaminated by bacteria. This requires constant monitoring of the circulating strains throughout the food production chain. Bacterial genomes contain signatures of natural evolution and adaptive markers that can be exploited to better understand the behavior of pathogen in the food industry. The monitoring of foodborne strains can therefore be facilitated by the use of these genomic markers capable of rapidly providing essential information on isolated strains, such as the source of contamination, risk of illness, potential for biofilm formation, and tolerance or resistance to biocides. The increasing availability of large genome datasets is enhancing the understanding of the genetic basis of complex traits such as host adaptation, virulence, and persistence. Genome-wide association studies have shown very promising results in the discovery of genomic markers that can be integrated into rapid detection tools. In addition, machine learning has successfully predicted phenotypes and classified important traits. Genome-wide association and machine learning tools have therefore the potential to support decision-making circuits intending at reducing the burden of foodborne diseases. The aim of this chapter review is to provide knowledge on the use of these two methods in food microbiology and to recommend their use in the field.

Genomic and Transcriptomic Comparison Between Invasive Non-Typhoidal Salmonella and Non-Invasive Isolates

Invasive non-typhoidal Salmonella (iNTS) poses a significant threat to global public health. Salmonella enterica Enteritidis and Typhimurium are the primary serovars responsible for both invasive diseases and gastroenteritis. This study aims to investigate the genomic and transcriptomic differences between isolates associated with these contrasting clinical presentations. We retrieved genomes of Salmonella Enteritidis and Salmonella Typhimurium from Enterobase, utilizing blood and stool isolates as representatives for iNTS and non-iNTS, respectively. An indistinguishable phylogenetic relationship was revealed between the blood and stool isolates for both serovars. Few genes were specifically identified in iNTS. Random forest and principal coordinates analysis permitted moderate discrimination between the two sources of isolates based on overall genome content. Notably, the blood isolates of Salmonella Typhimurium displayed an elevated level of antimicrobial resistance and genome degradation compared to stool isolates. Meanwhile, transcriptome sequencing identified few genes that were differentially expressed between blood and stool isolates. Hierarchical clustering and principal component analysis did not effectively differentiate the expression profile of iNTS from non-iNTS. In summary, few genes could serve as reliable biomarkers to distinguish iNTS and non-iNTS at either the genomic or transcriptomic level. Nevertheless, iNTS has indeed accumulated subtle genomic differences from non-iNTS, which might contribute to invasiveness.

Use of whole genome sequencing for surveillance and control of foodborne diseases: status quo and quo vadis

Improvements in sequencing quality, availability, speed and costs results in an increased presence of genomics in infectious disease applications. Nevertheless, there are still hurdles in regard to the optimal use of WGS for public health purposes. Here, we discuss the current state ("status quo") and future directions ("quo vadis") based on literature regarding the use of genomics in surveillance, hazard characterization and source attribution of foodborne pathogens. The future directions include the application of new techniques, such as machine learning and network approaches that may overcome the current shortcomings. These include the use of fixed genomic distances in cluster delineation, disentangling similarity or lack thereof in source attribution, and difficulties ascertaining function in hazard characterization. Although, the aforementioned methods can relatively easily be applied technically, an overarching challenge is the inference and biological/epidemiological interpretation of these large amounts of high-resolution data. Understanding the context in terms of bacterial isolate and host diversity allows to assess the level of representativeness in regard to sources and isolates in the dataset, which in turn defines the level of certainty associated with defining clusters, sources and risks. This also marks the importance of metadata (clinical, epidemiological, and biological) when using genomics for public health purposes.

Characterization of the prevalence of Salmonella in different retail chicken supply modes using genome-wide and machine-learning analyses

Salmonella is a foodborne pathogen that causes salmonellosis, of which retail chicken meat is a major source. However, the prevalence of Salmonella in different retail chicken supply modes and the threat posed to consumers remains unclear. The prevalence, serotype distribution, antibiotic resistance, and genomic characteristics of Salmonella in three supply modes of retail chicken (live poultry, frozen, and chilled) were investigated using whole-genome sequencing (WGS) and machine learning (ML). In this study, 480 retail chicken samples from live poultry, frozen, and chilled supply modes in Guangzhou from 2020 to 2021, as well as 253 Salmonella isolates (total isolation rate = 53.1 %), were collected. The prevalence of isolates in the live poultry mode (67.5 %, 81/120) was statistically higher than in the frozen (50.0 %, 120/240) and chilled (43.3 %, 52/120) (P < 0.05) modes. Serotype identification showed significant differences in the serotype distribution of Salmonella in different supply modes. S. Enteritis (46.7 %) and S. Indiana (14.2 %) were predominant in the frozen mode. S. Agona (23.5 %) and S. Saintpaul (13.6 %) were predominant in live poultry, while S. Enteritis (40.4 %) and S. Kentucky (17.3 %) were predominant in chilled mode. Antibiotic testing showed that frozen mode isolates were more resistant; the multidrug-resistant (MDR) rate of isolates in the frozen mode reached 91.8 %, significantly higher than in the chilled (86.5 %) and live (74.1 %) (P < 0.05) modes. WGS was performed on 155 top serotypes (S. Enteritidis, S. Kentucky, S. Indiana, and S. Agona). The antibiotic resistance gene analysis showed that the abundance and carrying rate of antibiotic resistance genes of Salmonella in the frozen mode (54 types, 16.1 %) were significantly higher than in other modes (live poultry: 36 types, 9.4 %, P < 0.05; chilled: 31 types, 11.6 %). The blaNDM-1 and blaNDM-9 genes encoding carbapenem resistance were found in frozen mode isolates on a complex transposon consisting of TnAS3-IS26. Virulence factors and plasmid replicons were abundant in the studied frozen mode isolates. In addition, single nucleotide polymorphism (SNP) phylogenetic tree results showed that in the frozen supply mode, the S. Enteritidis clonal clade continued to contaminate retail chicken meat and was homologous to S. Enteritidis strains found in farm chicken embryos, slaughterhouse chicken carcasses, and patients from hospitals in China (SNP 0 = 10). Notably, the pan-genome-based ML model showed that characteristic genes in frozen and live poultry isolates differed. The narZ gene was a key characteristic gene in frozen isolates, encoding nitrate reductase, relating to anaerobic bacterial growth. The ydgJ gene is a key characteristic gene in the live mode and encodes an oxidoreductase related to oxidative function in bacteria. The high prevalence of live poultry mode Salmonella and the transmission of frozen mode MDR Salmonella in this study pose serious risks to food safety and public health, emphasizing the importance of improving disinfection and cold storage measures to reduce Salmonella contamination and transmission. In conclusion, the continued surveillance of Salmonella across different supply models and the development of an epidemiological surveillance system based on WGS is necessary.

Generalizability of machine learning in predicting antimicrobial resistance in E. coli: a multi-country case study in Africa

BACKGROUND

Antimicrobial resistance (AMR) remains a significant global health threat particularly impacting low- and middle-income countries (LMICs). These regions often grapple with limited healthcare resources and access to advanced diagnostic tools. Consequently, there is a pressing need for innovative approaches that can enhance AMR surveillance and management. Machine learning (ML) though underutilized in these settings, presents a promising avenue. This study leverages ML models trained on whole-genome sequencing data from England, where such data is more readily available, to predict AMR in E. coli, targeting key antibiotics such as ciprofloxacin, ampicillin, and cefotaxime. A crucial part of our work involved the validation of these models using an independent dataset from Africa, specifically from Uganda, Nigeria, and Tanzania, to ascertain their applicability and effectiveness in LMICs.

RESULTS

Model performance varied across antibiotics. The Support Vector Machine excelled in predicting ciprofloxacin resistance (87% accuracy, F1 Score: 0.57), Light Gradient Boosting Machine for cefotaxime (92% accuracy, F1 Score: 0.42), and Gradient Boosting for ampicillin (58% accuracy, F1 Score: 0.66). In validation with data from Africa, Logistic Regression showed high accuracy for ampicillin (94%, F1 Score: 0.97), while Random Forest and Light Gradient Boosting Machine were effective for ciprofloxacin (50% accuracy, F1 Score: 0.56) and cefotaxime (45% accuracy, F1 Score:0.54), respectively. Key mutations associated with AMR were identified for these antibiotics.

CONCLUSION

As the threat of AMR continues to rise, the successful application of these models, particularly on genomic datasets from LMICs, signals a promising avenue for improving AMR prediction to support large AMR surveillance programs. This work thus not only expands our current understanding of the genetic underpinnings of AMR but also provides a robust methodological framework that can guide future research and applications in the fight against AMR.

Holistic One Health Surveillance Framework: Synergizing Environmental, Animal, and Human Determinants for Enhanced Infectious Disease Management

Recent pandemics, including the COVID-19 outbreak, have brought up growing concerns about transmission of zoonotic diseases from animals to humans. This highlights the requirement for a novel approach to discern and address the escalating health threats. The One Health paradigm has been developed as a responsive strategy to confront forthcoming outbreaks through early warning, highlighting the interconnectedness of humans, animals, and their environment. The system employs several innovative methods such as the use of advanced technology, global collaboration, and data-driven decision-making to come up with an extraordinary solution for improving worldwide disease responses. This Review deliberates environmental, animal, and human factors that influence disease risk, analyzes the challenges and advantages inherent in using the One Health surveillance system, and demonstrates how these can be empowered by Big Data and Artificial Intelligence. The Holistic One Health Surveillance Framework presented herein holds the potential to revolutionize our capacity to monitor, understand, and mitigate the impact of infectious diseases on global populations.

Genomic surveillance for antimicrobial resistance - a One Health perspective

Antimicrobial resistance (AMR) - the ability of microorganisms to adapt and survive under diverse chemical selection pressures - is influenced by complex interactions between humans, companion and food-producing animals, wildlife, insects and the environment. To understand and manage the threat posed to health (human, animal, plant and environmental) and security (food and water security and biosecurity), a multifaceted 'One Health' approach to AMR surveillance is required. Genomic technologies have enabled monitoring of the mobilization, persistence and abundance of AMR genes and mutations within and between microbial populations. Their adoption has also allowed source-tracing of AMR pathogens and modelling of AMR evolution and transmission. Here, we highlight recent advances in genomic AMR surveillance and the relative strengths of different technologies for AMR surveillance and research. We showcase recent insights derived from One Health genomic surveillance and consider the challenges to broader adoption both in developed and in lower- and middle-income countries.

The leading global health challenges in the artificial intelligence era

Millions of people's health is at risk because of several factors and multiple overlapping crises, all of which hit the vulnerable the most. These challenges are dynamic and evolve in response to emerging health challenges and concerns, which need effective collaboration among countries working toward achieving Sustainable Development Goals (SDGs) and securing global health. Mental Health, the Impact of climate change, cardiovascular diseases (CVDs), diabetes, Infectious diseases, health system, and population aging are examples of challenges known to pose a vast burden worldwide. We are at a point known as the "digital revolution," characterized by the expansion of artificial intelligence (AI) and a fusion of technology types. AI has emerged as a powerful tool for addressing various health challenges, and the last ten years have been influential due to the rapid expansion in the production and accessibility of health-related data. The computational models and algorithms can understand complicated health and medical data to perform various functions and deep-learning strategies. This narrative mini-review summarizes the most current AI applications to address the leading global health challenges. Harnessing its capabilities can ultimately mitigate the Impact of these challenges and revolutionize the field. It has the ability to strengthen global health through personalized health care and improved preparedness and response to future challenges. However, ethical and legal concerns about individual or community privacy and autonomy must be addressed for effective implementation.

The advantage of intergenic regions as genomic features for machine-learning-based host attribution of Salmonella Typhimurium from the USA

Salmonella enterica is a taxonomically diverse pathogen with over 2600 serovars associated with a wide variety of animal hosts including humans, other mammals, birds and reptiles. Some serovars are host-specific or host-restricted and cause disease in distinct host species, while others, such as serovar S. Typhimurium (STm), are generalists and have the potential to colonize a wide variety of species. However, even within generalist serovars such as STm it is becoming clear that pathovariants exist that differ in tropism and virulence. Identifying the genetic factors underlying host specificity is complex, but the availability of thousands of genome sequences and advances in machine learning have made it possible to build specific host prediction models to aid outbreak control and predict the human pathogenic potential of isolates from animals and other reservoirs. We have advanced this area by building host-association prediction models trained on a wide range of genomic features and compared them with predictions based on nearest-neighbour phylogeny. SNPs, protein variants (PVs), antimicrobial resistance (AMR) profiles and intergenic regions (IGRs) were extracted from 3883 high-quality STm assemblies collected from humans, swine, bovine and poultry in the USA, and used to construct Random Forest (RF) machine learning models. An additional 244 recent STm assemblies from farm animals were used as a test set for further validation. The models based on PVs and IGRs had the best performance in terms of predicting the host of origin of isolates and outperformed nearest-neighbour phylogenetic host prediction as well as models based on SNPs or AMR data. However, the models did not yield reliable predictions when tested with isolates that were phylogenetically distinct from the training set. The IGR and PV models were often able to differentiate human isolates in clusters where the majority of isolates were from a single animal source. Notably, IGRs were the feature with the best performance across multiple models which may be due to IGRs acting as both a representation of their flanking genes, equivalent to PVs, while also capturing genomic regulatory variation, such as altered promoter regions. The IGR and PV models predict that ~45 % of the human infections with STm in the USA originate from bovine, ~40 % from poultry and ~14.5 % from swine, although sequences of isolates from other sources were not used for training. In summary, the research demonstrates a significant gain in accuracy for models with IGRs and PVs as features compared to SNP-based and core genome phylogeny predictions when applied within the existing population structure. This article contains data hosted by Microreact.

Predictive modeling of antibiotic eradication therapy success for new-onset Pseudomonas aeruginosa pulmonary infections in children with cystic fibrosis

Chronic Pseudomonas aeruginosa (Pa) lung infections are the leading cause of mortality among cystic fibrosis (CF) patients; therefore, the eradication of new-onset Pa lung infections is an important therapeutic goal that can have long-term health benefits. The use of early antibiotic eradication therapy (AET) has been shown to clear the majority of new-onset Pa infections, and it is hoped that identifying the underlying basis for AET failure will further improve treatment outcomes. Here we generated machine learning models to predict AET outcomes based on pathogen genomic data. We used a nested cross validation design, population structure control, and recursive feature selection to improve model performance and showed that incorporating population structure control was crucial for improving model interpretation and generalizability. Our best model, controlling for population structure and using only 30 recursively selected features, had an area under the curve of 0.87 for a holdout test dataset. The top-ranked features were generally associated with motility, adhesion, and biofilm formation.

Screening under infection-relevant conditions reveals chemical sensitivity in multidrug resistant invasive non-typhoidal Salmonella (iNTS)

Bloodstream infections caused by invasive, non-typhoidal Salmonella (iNTS) are a major global health concern, particularly in Africa where the pathogenic variant of Salmonella Typhimurium sequence type (ST) 313 is dominant. Unlike S. Typhimurium strains that cause gastroenteritis, iNTS strains cause bloodstream infections and are resistant to multiple first-line antibiotics, thus limiting current treatment options. Here, we developed and implemented multiple small molecule screens under physiological, infection-relevant conditions to reveal chemical sensitivities in ST313 and to identify host-directed therapeutics as entry points to drug discovery to combat the clinical burden of iNTS. Screening ST313 iNTS under host-mimicking growth conditions identified 92 compounds with antimicrobial activity despite inherent multidrug resistance. We characterized the antimicrobial activity of the nucleoside analog 3'-azido-3'-deoxythymidine as an exemplary compound from this screen, which depended on bacterial thymidine kinase activity for antimicrobial activity. In a companion macrophage-based screening platform designed to enrich for host-directed therapeutics, we identified three compounds (amodiaquine, berbamine, and indatraline) as actives that required the presence of host cells for antibacterial activity. These three compounds had antimicrobial activity only in the presence of host cells that significantly inhibited intracellular ST313 iNTS replication in macrophages. This work provides evidence that despite high invasiveness and multidrug resistance, ST313 iNTS remains susceptible to unconventional drug discovery approaches.

Leveraging artificial intelligence in the fight against infectious diseases

Despite advances in molecular biology, genetics, computation, and medicinal chemistry, infectious disease remains an ominous threat to public health. Addressing the challenges posed by pathogen outbreaks, pandemics, and antimicrobial resistance will require concerted interdisciplinary efforts. In conjunction with systems and synthetic biology, artificial intelligence (AI) is now leading to rapid progress, expanding anti-infective drug discovery, enhancing our understanding of infection biology, and accelerating the development of diagnostics. In this Review, we discuss approaches for detecting, treating, and understanding infectious diseases, underscoring the progress supported by AI in each case. We suggest future applications of AI and how it might be harnessed to help control infectious disease outbreaks and pandemics.

Plasmid permissiveness of wastewater microbiomes can be predicted from 16S rRNA sequences by machine learning

MOTIVATION

Wastewater treatment plants (WWTPs) harbor a dense and diverse microbial community. They constantly receive antimicrobial residues and resistant strains, and therefore provide conditions for horizontal gene transfer (HGT) of antimicrobial resistance (AMR) determinants. This facilitates the transmission of clinically important genes between, e.g. enteric and environmental bacteria, and vice versa. Despite the clinical importance, tools for predicting HGT remain underdeveloped.

RESULTS

In this study, we examined to which extent water cycle microbial community composition, as inferred by partial 16S rRNA gene sequences, can predict plasmid permissiveness, i.e. the ability of cells to receive a plasmid through conjugation, based on data from standardized filter mating assays using fluorescent bio-reporter plasmids. We leveraged a range of machine learning models for predicting the permissiveness for each taxon in the community, representing the range of hosts a plasmid is able to transfer to, for three broad host-range resistance IncP plasmids (pKJK5, pB10, and RP4). Our results indicate that the predicted permissiveness from the best performing model (random forest) showed a moderate-to-strong average correlation of 0.49 for pB10 [95% confidence interval (CI): 0.44-0.55], 0.43 for pKJK5 (0.95% CI: 0.41-0.49), and 0.53 for RP4 (0.95% CI: 0.48-0.57) with the experimental permissiveness in the unseen test dataset. Predictive phylogenetic signals occurred despite the broad host-range nature of these plasmids. Our results provide a framework that contributes to the assessment of the risk of AMR pollution in wastewater systems.

AVAILABILITY AND IMPLEMENTATION

The predictive tool is available as an application at https://github.com/DaneshMoradigaravand/PlasmidPerm.

Rapid geographical source attribution of Salmonella enterica serovar Enteritidis genomes using hierarchical machine learning

Salmonella enterica serovar Enteritidis is one of the most frequent causes of Salmonellosis globally and is commonly transmitted from animals to humans by the consumption of contaminated foodstuffs. In the UK and many other countries in the Global North, a significant proportion of cases are caused by the consumption of imported food products or contracted during foreign travel, therefore, making the rapid identification of the geographical source of new infections a requirement for robust public health outbreak investigations. Herein, we detail the development and application of a hierarchical machine learning model to rapidly identify and trace the geographical source of S. Enteritidis infections from whole genome sequencing data. 2313 S. Enteritidis genomes, collected by the UKHSA between 2014-2019, were used to train a 'local classifier per node' hierarchical classifier to attribute isolates to four continents, 11 sub-regions, and 38 countries (53 classes). The highest classification accuracy was achieved at the continental level followed by the sub-regional and country levels (macro F1: 0.954, 0.718, 0.661, respectively). A number of countries commonly visited by UK travelers were predicted with high accuracy (hF1: >0.9). Longitudinal analysis and validation with publicly accessible international samples indicated that predictions were robust to prospective external datasets. The hierarchical machine learning framework provided granular geographical source prediction directly from sequencing reads in <4 min per sample, facilitating rapid outbreak resolution and real-time genomic epidemiology. The results suggest additional application to a broader range of pathogens and other geographically structured problems, such as antimicrobial resistance prediction, is warranted.

Development of a novel machine learning-based weighted modeling approach to incorporate Salmonella enterica heterogeneity on a genetic scale in a dose-response modeling framework

Estimating microbial dose-response is an important aspect of a food safety risk assessment. In recent years, there has been considerable interest to advance these models with potential incorporation of gene expression data. The aim of this study was to develop a novel machine learning model that considers the weights of expression of Salmonella genes that could be associated with illness, given exposure, in hosts. Here, an elastic net-based weighted Poisson regression method was proposed to identify Salmonella enterica genes that could be significantly associated with the illness response, irrespective of serovar. The best-fit elastic net model was obtained by 10-fold cross-validation. The best-fit elastic net model identified 33 gene expression-dose interaction terms that added to the predictability of the model. Of these, nine genes associated with Salmonella metabolism and virulence were found to be significant by the best-fit Poisson regression model (p < 0.05). This method could improve or redefine dose-response relationships for illness from relative proportions of significant genes from a microbial genetic dataset, which would help in refining endpoint and risk estimations.

Computational biology: Role and scope in taming antimicrobial resistance

BACKGROUND

Infectious diseases pose many challenges due to increasing threat of antimicrobial resistance, which necessitates continuous research to develop novel strategies for development of new molecules with antibacterial activity. In the era of computational biology there are tools and techniques available to address and solve the disease management issues in the field of clinical microbiology. The sequencing techniques, structural biology and machine learning can be implemented collectively to tackle infectious diseases e.g. for the diagnosis, epidemiological typing, pathotyping, antimicrobial resistance detection as well as the discovery of novel drugs and vaccine biomarkers.

OBJECTIVES

The present review is a narrative review representing a comprehensive literature-based assessment regarding the use of whole genome sequencing, structural biology and machine learning for the diagnosis, molecular typing and antibacterial drug discovery.

CONTENT

Here, we seek to present an overview of molecular and structural basis of resistance to antibiotics, while focusing on the recent bioinformatics approaches in whole genome sequencing and structural biology. The application of next generation sequencing in management of bacterial infections focusing on investigation of microbial population diversity, genotypic resistance testing and scope for the identification of targets for novel drug and vaccine candidates, has been addressed along with the use of structural biophysics and artificial intelligence.

Paving the way for precise diagnostics of antimicrobial resistant bacteria

The antimicrobial resistance (AMR) crisis from bacterial pathogens is frequently emerging and rapidly disseminated during the sustained antimicrobial exposure in human-dominated communities, posing a compelling threat as one of the biggest challenges in humans. The frequent incidences of some common but untreatable infections unfold the public health catastrophe that antimicrobial-resistant pathogens have outpaced the available countermeasures, now explicitly amplified during the COVID-19 pandemic. Nowadays, biotechnology and machine learning advancements help create more fundamental knowledge of distinct spatiotemporal dynamics in AMR bacterial adaptation and evolutionary processes. Integrated with reliable diagnostic tools and powerful analytic approaches, a collaborative and systematic surveillance platform with high accuracy and predictability should be established and implemented, which is not just for an effective controlling strategy on AMR but also for protecting the longevity of valuable antimicrobials currently and in the future.

KLFDAPC: a supervised machine learning approach for spatial genetic structure analysis

Geographic patterns of human genetic variation provide important insights into human evolution and disease. A commonly used tool to detect and describe them is principal component analysis (PCA) or the supervised linear discriminant analysis of principal components (DAPC). However, genetic features produced from both approaches could fail to correctly characterize population structure for complex scenarios involving admixture. In this study, we introduce Kernel Local Fisher Discriminant Analysis of Principal Components (KLFDAPC), a supervised non-linear approach for inferring individual geographic genetic structure that could rectify the limitations of these approaches by preserving the multimodal space of samples. We tested the power of KLFDAPC to infer population structure and to predict individual geographic origin using neural networks. Simulation results showed that KLFDAPC has higher discriminatory power than PCA and DAPC. The application of our method to empirical European and East Asian genome-wide genetic datasets indicated that the first two reduced features of KLFDAPC correctly recapitulated the geography of individuals and significantly improved the accuracy of predicting individual geographic origin when compared to PCA and DAPC. Therefore, KLFDAPC can be useful for geographic ancestry inference, design of genome scans and correction for spatial stratification in GWAS that link genes to adaptation or disease susceptibility.

A Machine Learning Model for Food Source Attribution of Listeria monocytogenes

Despite its low morbidity, listeriosis has a high mortality rate due to the severity of its clinical manifestations. The source of human listeriosis is often unclear. In this study, we investigate the ability of machine learning to predict the food source from which clinical Listeria monocytogenes isolates originated. Four machine learning classification algorithms were trained on core genome multilocus sequence typing data of 1212 L. monocytogenes isolates from various food sources. The average accuracies of random forest, support vector machine radial kernel, stochastic gradient boosting, and logit boost were found to be 0.72, 0.61, 0.7, and 0.73, respectively. Logit boost showed the best performance and was used in model testing on 154 L. monocytogenes clinical isolates. The model attributed 17.5 % of human clinical cases to dairy, 32.5% to fruits, 14.3% to leafy greens, 9.7% to meat, 4.6% to poultry, and 18.8% to vegetables. The final model also provided us with genetic features that were predictive of specific sources. Thus, this combination of genomic data and machine learning-based models can greatly enhance our ability to track L. monocytogenes from different food sources.

Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections

During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.

Salmonella enterica serovar Typhimurium from Wild Birds in the United States Represent Distinct Lineages Defined by Bird Type

Salmonella enterica serovar Typhimurium is typically considered a host generalist; however, certain isolates are associated with specific hosts and show genetic features of host adaptation. Here, we sequenced 131 S. Typhimurium isolates from wild birds collected in 30 U.S. states during 1978–2019. We found that isolates from broad taxonomic host groups including passerine birds, water birds (Aequornithes), and larids (gulls and terns) represented three distinct lineages and certain S. Typhimurium CRISPR types presented in individual lineages. We also showed that lineages formed by wild bird isolates differed from most isolates originating from domestic animal sources, and that genomes from these lineages substantially improved source attribution of Typhimurium genomes to wild birds by a machine learning classifier. Furthermore, virulence gene signatures that differentiated S. Typhimurium from passerines, water birds, and larids were detected. Passerine isolates tended to lack S. Typhimurium-specific virulence plasmids. Isolates from the passerine, water bird, and larid lineages had close genetic relatedness with human clinical isolates, including those from a 2021 U.S. outbreak linked to passerine birds. These observations indicate that S. Typhimurium from wild birds in the United States are likely host-adapted, and the representative genomic data set examined in this study can improve source prediction and facilitate outbreak investigation.

IMPORTANCE Within-host evolution of S. Typhimurium may lead to pathovars adapted to specific hosts. Here, we report the emergence of disparate avian S. Typhimurium lineages with distinct virulence gene signatures. The findings highlight the importance of wild birds as a reservoir for S. Typhimurium and contribute to our understanding of the genetic diversity of S. Typhimurium from wild birds. Our study indicates that S. Typhimurium may have undergone adaptive evolution within wild birds in the United States. The representative S. Typhimurium genomes from wild birds, together with the virulence gene signatures identified in these bird isolates, are valuable for S. Typhimurium source attribution and epidemiological surveillance.

Exploring the predictive capability of advanced machine learning in identifying severe disease phenotype in Salmonella enterica

The past few years have seen a significant increase in availability of whole genome sequencing information, allowing for its incorporation in predictive modeling for foodborne pathogens to account for inter- and intra-species differences in their virulence. However, this is hindered by the inability of traditional statistical methods to analyze such large amounts of data compared to the number of observations/isolates. In this study, we have explored the applicability of machine learning (ML) models to predict the disease outcome, while identifying features that exert a significant effect on the prediction. This study was conducted on Salmonella enterica, a major foodborne pathogen with considerable inter- and intra-serovar variation. WGS of isolates obtained from various sources (i.e., human, chicken, and swine) were used as input in four machine learning models (logistic regression with ridge, random forest, support vector machine, and AdaBoost) to classify isolates based on disease severity (extraintestinal vs. gastrointestinal) in the host. The predictive performances of all models were tested with and without Elastic Net regularization to combat dimensionality issues. Elastic Net-regularized logistic regression model showed the best area under the receiver operating characteristic curve (AUC-ROC; 0.86) and outcome prediction accuracy (0.76). Additionally, genes coding for transcriptional regulation, acidic, oxidative, and anaerobic stress response, and antibiotic resistance were found to be significant predictors of disease severity. These genes, which were significantly associated with each outcome, could possibly be input in amended, gene-expression-specific predictive models to estimate virulence pattern-specific effect of Salmonella and other foodborne pathogens on human health.

A genomic data resource for predicting antimicrobial resistance from laboratory-derived antimicrobial susceptibility phenotypes

Antimicrobial resistance (AMR) is a major global health threat that affects millions of people each year. Funding agencies worldwide and the global research community have expended considerable capital and effort tracking the evolution and spread of AMR by isolating and sequencing bacterial strains and performing antimicrobial susceptibility testing (AST). For the last several years, we have been capturing these efforts by curating data from the literature and data resources and building a set of assembled bacterial genome sequences that are paired with laboratory-derived AST data. This collection currently contains AST data for over 67 000 genomes encompassing approximately 40 genera and over 100 species. In this paper, we describe the characteristics of this collection, highlighting areas where sampling is comparatively deep or shallow, and showing areas where attention is needed from the research community to improve sampling and tracking efforts. In addition to using the data to track the evolution and spread of AMR, it also serves as a useful starting point for building machine learning models for predicting AMR phenotypes. We demonstrate this by describing two machine learning models that are built from the entire dataset to show where the predictive power is comparatively high or low. This AMR metadata collection is freely available and maintained on the Bacterial and Viral Bioinformatics Center (BV-BRC) FTP site ftp://ftp.bvbrc.org/RELEASE_NOTES/PATRIC_genomes_AMR.txt.

Genomic signatures of host adaptation in group B Salmonella enterica ST416/ST417 from harbour porpoises

A type of monophasic group B Salmonella enterica with the antigenic formula 4,12:a:- (“Fulica-like”) has been described as associated with harbour porpoises (Phocoena phocoena), most frequently recovered from lung samples. In the present study, lung tissue samples from 47 porpoises found along the Swedish coast or as bycatch in fishing nets were analysed, two of which were positive for S. enterica. Pneumonia due to the infection was considered the likely cause of death for one of the two animals. The recovered isolates were whole genome sequenced and found to belong to sequence type (ST) 416 and to be closely related to ST416/ST417 porpoise isolates from UK waters as determined by core-genome MLST. Serovars Bispebjerg, Fulica and Abortusequi were identified as distantly related to the porpoise isolates, but no close relatives from other host species were found. All ST416/417 isolates had extensive loss of function mutations in key Salmonella pathogenicity islands, but carried accessory genetic elements associated with extraintestinal infection such as iron uptake systems. Gene ontology and pathway analysis revealed reduced secondary metabolic capabilities and loss of function in terms of signalling and response to environmental cues, consistent with adaptation for the extraintestinal niche. A classification system based on machine learning identified ST416/417 as more invasive than classical gastrointestinal serovars. Genome analysis results are thus consistent with ST416/417 as a host-adapted and extraintestinal clonal population of S. enterica, which while found in porpoises without associated pathology can also cause severe opportunistic infections.

Sugar-Phosphate Toxicities

Accumulation of phosphorylated intermediates during cellular metabolism can have wide-ranging toxic effects on many organisms, including humans and the pathogens that infect them. These toxicities can be induced by feeding an upstream metabolite (a sugar, for instance) while simultaneously blocking the appropriate metabolic pathway with either a mutation or an enzyme inhibitor. Here, we survey the toxicities that can arise in the metabolism of glucose, galactose, fructose, fructose-asparagine, glycerol, trehalose, maltose, mannose, mannitol, arabinose, and rhamnose. Select enzymes in these metabolic pathways may serve as novel therapeutic targets. Some are conserved broadly among prokaryotes and eukaryotes (e.g., glucose and galactose) and are therefore unlikely to be viable drug targets. However, others are found only in bacteria (e.g., fructose-asparagine, rhamnose, and arabinose), and one is found in fungi but not in humans (trehalose). We discuss what is known about the mechanisms of toxicity and how resistance is achieved in order to identify the prospects and challenges associated with targeted exploitation of these pervasive metabolic vulnerabilities.

Machine Learning Prediction of Resistance to Subinhibitory Antimicrobial Concentrations from Escherichia coli Genomes

Escherichia coli is an important cause of bacterial infections worldwide, with multidrug-resistant strains incurring substantial costs on human lives. Besides therapeutic concentrations of antimicrobials in health care settings, the presence of subinhibitory antimicrobial residues in the environment and in clinics selects for antimicrobial resistance (AMR), but the underlying genetic repertoire is less well understood. Here, we used machine learning to predict the population doubling time and cell growth yield of 1,407 genetically diverse E. coli strains expanding under exposure to three subinhibitory concentrations of six classes of antimicrobials from single-nucleotide genetic variants, accessory gene variation, and the presence of known AMR genes. We predicted cell growth yields in the held-out test data with an average correlation (Spearman's ρ) of 0.63 (0.36 to 0.81 across concentrations) and cell doubling times with an average correlation of 0.59 (0.32 to 0.92 across concentrations), with moderate increases in sample size unlikely to improve predictions further. This finding points to the remaining missing heritability of growth under antimicrobial exposure being explained by effects that are too rare or weak to be captured unless sample size is dramatically increased, or by effects other than those conferred by the presence of individual single-nucleotide polymorphisms (SNPs) and genes. Predictions based on whole-genome information were generally superior to those based only on known AMR genes and were accurate for AMR resistance at therapeutic concentrations. We pinpointed genes and SNPs determining the predicted growth and thereby recapitulated many known AMR determinants. Finally, we estimated the effect sizes of resistance genes across the entire collection of strains, disclosing the growth effects for known resistance genes in each individual strain. Our results underscore the potential of predictive modeling of growth patterns from genomic data under subinhibitory concentrations of antimicrobials, although the remaining missing heritability poses a challenge for achieving the accuracy and precision required for clinical use. IMPORTANCE Predicting bacterial growth from genome sequences is important for a rapid characterization of strains in clinical diagnostics and to disclose candidate novel targets for anti-infective drugs. Previous studies have dissected the relationship between bacterial growth and genotype in mutant libraries for laboratory strains, yet no study so far has examined the predictive power of genome sequence in natural strains. In this study, we used a high-throughput phenotypic assay to measure the growth of a systematic collection of natural Escherichia coli strains and then employed machine learning models to predict bacterial growth from genomic data under nontherapeutic subinhibitory concentrations of antimicrobials that are common in nonclinical settings. We found a moderate to strong correlation between predicted and actual values for the different collected data sets. Moreover, we observed that the known resistance genes are still effective at sublethal concentrations, pointing to clinical implications of these concentrations.

Campylobacter jejuni genotypes are associated with post-infection irritable bowel syndrome in humans

Campylobacter enterocolitis may lead to post-infection irritable bowel syndrome (PI-IBS) and while some C. jejuni strains are more likely than others to cause human disease, genomic and virulence characteristics promoting PI-IBS development remain uncharacterized. We combined pangenome-wide association studies and phenotypic assays to compare C. jejuni isolates from patients who developed PI-IBS with those who did not. We show that variation in bacterial stress response (Cj0145_phoX), adhesion protein (Cj0628_CapA), and core biosynthetic pathway genes (biotin: Cj0308_bioD; purine: Cj0514_purQ; isoprenoid: Cj0894c_ispH) were associated with PI-IBS development. In vitro assays demonstrated greater adhesion, invasion, IL-8 and TNFα secretion on colonocytes with PI-IBS compared to PI-no-IBS strains. A risk-score for PI-IBS development was generated using 22 genomic markers, four of which were from Cj1631c, a putative heme oxidase gene linked to virulence. Our finding that specific Campylobacter genotypes confer greater in vitro virulence and increased risk of PI-IBS has potential to improve understanding of the complex host-pathogen interactions underlying this condition.

Stephanie Peters, Ben Pascoe, et al. use whole-genome sequencing and phenotypic analysis of clinical strains from patients to identify potential genetic factors involved in irritable bowel syndrome resulting from Campylobacter jejuni infection. Their data suggest that genes involved in the bacterial stress response and biosynthetic pathways may contribute toward irritable bowel syndrome, providing further insight into links between Campylobacter genotypes and risk of disease.

The Avian Pathogenic Escherichia coli (APEC) pathotype is comprised of multiple distinct, independent genotypes

Avian Pathogenic E. coli (APEC) is the causative agent of avian colibacillosis, resulting in economic losses to the poultry industry through morbidity, mortality and carcass condemnation, and impacts the welfare of poultry. Colibacillosis remains a complex disease to manage, hampered by diagnostic and classification strategies for E. coli that are inadequate for defining APEC. However, increased accessibility of whole genome sequencing (WGS) technology has enabled phylogenetic approaches to be applied to the classification of E. coli and genomic characterization of the most common APEC serotypes associated with colibacillosis O1, O2 and O78. These approaches have demonstrated that the O78 serotype is representative of two distinct APEC lineages, ST-23 in phylogroup C and ST-117 in phylogroup G. The O1 and O2 serotypes belong to a third lineage comprised of three sub-populations in phylogroup B2; ST-95, ST-140 and ST-428/ST-429. The frequency with which these genotypes are associated with colibacillosis implicates them as the predominant APEC populations and distinct from those causing incidental or opportunistic infections. The fact that these are disparate clusters from multiple phylogroups suggests that these lineages may have become adapted to the poultry niche independently. WGS studies have highlighted the limitations of traditional APEC classification and can now provide a path towards a robust and more meaningful definition of the APEC pathotype. Future studies should focus on characterizing individual APEC populations in detail and using this information to develop improved diagnostics and interventions.

Salmonella pathogenesis and host-adaptation in farmed animals

Salmonella is an animal and zoonotic pathogen of global importance. Depending on pathogen and host factors, infections can be asymptomatic or involve acute gastroenteritis or invasive disease. Genomic signatures associated with host-range, tissue tropism or differential virulence of Salmonella enterica serovars, and their variants, have emerged. In turn, it is becoming feasible to predict invasive potential, host-adaptation and zoonotic risk of Salmonella from sequence data to improve outbreak investigation, risk assessment and control strategies. Functional annotation of Salmonella genomes has accelerated with the screening of high-density mutant libraries, revealing host-specific, niche-specific and serovar-specific virulence factors. As natural hosts and reservoirs, farmed animals provide powerful insights into host-adaptation and pathogenesis of Salmonella not always evident from surrogate rodent or cell-based models.

ProkEvo: an automated, reproducible, and scalable framework for high-throughput bacterial population genomics analyses

Whole Genome Sequence (WGS) data from bacterial species is used for a variety of applications ranging from basic microbiological research, diagnostics, and epidemiological surveillance. The availability of WGS data from hundreds of thousands of individual isolates of individual microbial species poses a tremendous opportunity for discovery and hypothesis-generating research into ecology and evolution of these microorganisms. Flexibility, scalability, and user-friendliness of existing pipelines for population-scale inquiry, however, limit applications of systematic, population-scale approaches. Here, we present ProkEvo, an automated, scalable, reproducible, and open-source framework for bacterial population genomics analyses using WGS data. ProkEvo was specifically developed to achieve the following goals: (1) Automation and scaling of complex combinations of computational analyses for many thousands of bacterial genomes from inputs of raw Illumina paired-end sequence reads; (2) Use of workflow management systems (WMS) such as Pegasus WMS to ensure reproducibility, scalability, modularity, fault-tolerance, and robust file management throughout the process; (3) Use of high-performance and high-throughput computational platforms; (4) Generation of hierarchical-based population structure analysis based on combinations of multi-locus and Bayesian statistical approaches for classification for ecological and epidemiological inquiries; (5) Association of antimicrobial resistance (AMR) genes, putative virulence factors, and plasmids from curated databases with the hierarchically-related genotypic classifications; and (6) Production of pan-genome annotations and data compilation that can be utilized for downstream analysis such as identification of population-specific genomic signatures. The scalability of ProkEvo was measured with two datasets comprising significantly different numbers of input genomes (one with ~2,400 genomes, and the second with ~23,000 genomes). Depending on the dataset and the computational platform used, the running time of ProkEvo varied from ~3-26 days. ProkEvo can be used with virtually any bacterial species, and the Pegasus WMS uniquely facilitates addition or removal of programs from the workflow or modification of options within them. To demonstrate versatility of the ProkEvo platform, we performed a hierarchical-based population structure analyses from available genomes of three distinct pathogenic bacterial species as individual case studies. The specific case studies illustrate how hierarchical analyses of population structures, genotype frequencies, and distribution of specific gene functions can be integrated into an analysis. Collectively, our study shows that ProkEvo presents a practical viable option for scalable, automated analyses of bacterial populations with direct applications for basic microbiology research, clinical microbiological diagnostics, and epidemiological surveillance.

Machine learning approach to support taxonomic species discrimination based on helminth collections data

Background

There are more than 300 species of capillariids that parasitize various vertebrate groups worldwide. Species identification is hindered because of the few taxonomically informative structures available, making the task laborious and genus definition controversial. Thus, its taxonomy is one of the most complex among Nematoda. Eggs are the parasitic structures most viewed in coprological analysis in both modern and ancient samples; consequently, their presence is indicative of positive diagnosis for infection. The structure of the egg could play a role in genera or species discrimination. Institutional biological collections are taxonomic repositories of specimens described and strictly identified by systematics specialists.

Methods

The present work aims to characterize eggs of capillariid species deposited in institutional helminth collections and to process the morphological, morphometric and ecological data using machine learning (ML) as a new approach for taxonomic identification. Specimens of 28 species and 8 genera deposited at Coleção Helmintológica do Instituto Oswaldo Cruz (CHIOC, IOC/FIOCRUZ/Brazil) and Collection de Nématodes Zooparasites du Muséum National d’Histoire Naturelle de Paris (MNHN/France) were examined under light microscopy. In the morphological and morphometric analyses (MM), the total length and width of eggs as well as plugs and shell thickness were considered. In addition, eggshell ornamentations and ecological parameters of the geographical location (GL) and host (H) were included.

Results

The performance of the logistic model tree (LMT) algorithm showed the highest values in all metrics compared with the other algorithms. Algorithm J48 produced the most reliable decision tree for species identification alongside REPTree. The Majority Voting algorithm showed high metric values, but the combined classifiers did not attenuate the errors revealed in each algorithm alone. The statistical evaluation of the dataset indicated a significant difference between trees, with GL + H + MM and MM only with the best scores.

Conclusions

The present research proposed a novel procedure for taxonomic species identification, integrating data from centenary biological collections and the logic of artificial intelligence techniques. This study will support future research on taxonomic identification and diagnosis of both modern and archaeological capillariids.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04721-6.

Stepwise evolution of Salmonella Typhimurium ST313 causing bloodstream infection in Africa

Bloodstream infections caused by nontyphoidal Salmonella are a major public health concern in Africa, causing ~49,600 deaths every year. The most common Salmonella enterica pathovariant associated with invasive nontyphoidal Salmonella disease is Salmonella Typhimurium sequence type (ST)313. It has been proposed that antimicrobial resistance and genome degradation has contributed to the success of ST313 lineages in Africa, but the evolutionary trajectory of such changes was unclear. Here, to define the evolutionary dynamics of ST313, we sub-sampled from two comprehensive collections of Salmonella isolates from African patients with bloodstream infections, spanning 1966 to 2018. The resulting 680 genome sequences led to the discovery of a pan-susceptible ST313 lineage (ST313 L3), which emerged in Malawi in 2016 and is closely related to ST313 variants that cause gastrointestinal disease in the United Kingdom and Brazil. Genomic analysis revealed degradation events in important virulence genes in ST313 L3, which had not occurred in other ST313 lineages. Despite arising only recently in the clinic, ST313 L3 is a phylogenetic intermediate between ST313 L1 and L2, with a characteristic accessory genome. Our in-depth genotypic and phenotypic characterization identifies the crucial loss-of-function genetic events that occurred during the stepwise evolution of invasive S. Typhimurium across Africa.

Comparative genomic and phenotypic characterization of invasive non-typhoidal Salmonella isolates from Siaya, Kenya

Non-typhoidal Salmonella (NTS) is a major global health concern that often causes bloodstream infections in areas of the world affected by malnutrition and comorbidities such as HIV and malaria. Developing a strategy to control the emergence and spread of highly invasive and antimicrobial resistant NTS isolates requires a comprehensive analysis of epidemiological factors and molecular pathogenesis. Here, we characterize 11 NTS isolates that caused bloodstream infections in pediatric patients in Siaya, Kenya from 2003-2010. Nine isolates were identified as S. Typhimurium sequence type 313 while the other two were S. Enteritidis. Comprehensive genotypic and phenotypic analyses were performed to compare these isolates to those previously identified in sub-Saharan Africa. We identified a S. Typhimurium isolate referred to as UGA14 that displayed novel plasmid, pseudogene and resistance features as compared to other isolates reported from Africa. Notably, UGA14 is able to ferment both lactose and sucrose due to the acquisition of insertion elements on the pKST313 plasmid. These findings show for the first time the co-evolution of plasmid-mediated lactose and sucrose metabolism along with cephalosporin resistance in NTS further elucidating the evolutionary mechanisms of invasive NTS phenotypes. These results further support the use of combined genomic and phenotypic approaches to detect and characterize atypical NTS isolates in order to advance biosurveillance efforts that inform countermeasures aimed at controlling invasive and antimicrobial resistant NTS.

Emerging Applications of Machine Learning in Food Safety

Food safety continues to threaten public health. Machine learning holds potential in leveraging large, emerging data sets to improve the safety of the food supply and mitigate the impact of food safety incidents. Foodborne pathogen genomes and novel data streams, including text, transactional, and trade data, have seen emerging applications enabled by a machine learning approach, such as prediction of antibiotic resistance, source attribution of pathogens, and foodborne outbreak detection and risk assessment. In this article, we provide a gentle introduction to machine learning in the context of food safety and an overview of recent developments and applications. With many of these applications still in their nascence, general and domain-specific pitfalls and challenges associated with machine learning have begun to be recognized and addressed, which are critical to prospective use and future deployment of large data sets and their associated machine learning models for food safety applications.

Machine learning and statistics to qualify environments through multi-traits in Coffea arabica

Several factors such as genotype, environment, and post-harvest processing can affect the responses of important traits in the coffee production chain. Determining the influence of these factors is of great relevance, as they can be indicators of the characteristics of the coffee produced. The most efficient models choice to be applied should take into account the variety of information and the particularities of each biological material. This study was developed to evaluate statistical and machine learning models that would better discriminate environments through multi-traits of coffee genotypes and identify the main agronomic and beverage quality traits responsible for the variation of the environments. For that, 31 morpho-agronomic and post-harvest traits were evaluated, from field experiments installed in three municipalities in the Matas de Minas region, in the State of Minas Gerais, Brazil. Two types of post-harvest processing were evaluated: natural and pulped. The apparent error rate was estimated for each method. The Multilayer Perceptron and Radial Basis Function networks were able to discriminate the coffee samples in multi-environment more efficiently than the other methods, identifying differences in multi-traits responses according to the production sites and type of post-harvest processing. The local factors did not present specific traits that favored the severity of diseases and differentiated vegetative vigor. Sensory traits acidity and fragrance/aroma score also made little contribution to the discrimination process, indicating that acidity and fragrance/aroma are characteristic of coffee produced and all coffee samples evaluated are of the special type in the Mata of Minas region. The main traits responsible for the differentiation of production sites are plant height, fruit size, and bean production. The sensory trait "Body" is the main one to discriminate the form of post-harvest processing.

Salmonella versus the Microbiome

A balanced gut microbiota contributes to health, but the mechanisms maintaining homeostasis remain elusive. Microbiota assembly during infancy is governed by competition between species and by environmental factors, termed habitat filters, that determine the range of successful traits within the microbial community. These habitat filters include the diet, host-derived resources, and microbiota-derived metabolites, such as short-chain fatty acids. Once the microbiota has matured, competition and habitat filtering prevent engraftment of new microbes, thereby providing protection against opportunistic infections. Competition with endogenous Enterobacterales, habitat filtering by short-chain fatty acids, and a host-derived habitat filter, epithelial hypoxia, also contribute to colonization resistance against Salmonella serovars. However, at a high challenge dose, these frank pathogens can overcome colonization resistance by using their virulence factors to trigger intestinal inflammation. In turn, inflammation increases the luminal availability of host-derived resources, such as oxygen, nitrate, tetrathionate, and lactate, thereby creating a state of abnormal habitat filtering that enables the pathogen to overcome growth inhibition by short-chain fatty acids. Thus, studying the process of ecosystem invasion by Salmonella serovars clarifies that colonization resistance can become weakened by disrupting host-mediated habitat filtering. This insight is relevant for understanding how inflammation triggers dysbiosis linked to noncommunicable diseases, conditions in which endogenous Enterobacterales expand in the fecal microbiota using some of the same growth-limiting resources required by Salmonella serovars for ecosystem invasion. In essence, ecosystem invasion by Salmonella serovars suggests that homeostasis and dysbiosis simply represent states where competition and habitat filtering are normal or abnormal, respectively.

Identification of Nitrogen Fixation Genes in Lactococcus Isolated from Maize Using Population Genomics and Machine Learning

Sierra Mixe maize is a landrace variety from Oaxaca, Mexico, that utilizes nitrogen derived from the atmosphere via an undefined nitrogen fixation mechanism. The diazotrophic microbiota associated with the plant’s mucilaginous aerial root exudate composed of complex carbohydrates was previously identified and characterized by our group where we found 23 lactococci capable of biological nitrogen fixation (BNF) without containing any of the proposed essential genes for this trait (nifHDKENB). To determine the genes in Lactococcus associated with this phenotype, we selected 70 lactococci from the dairy industry that are not known to be diazotrophic to conduct a comparative population genomic analysis. This showed that the diazotrophic lactococcal genomes were distinctly different from the dairy isolates. Examining the pangenome followed by genome-wide association study and machine learning identified genes with the functions needed for BNF in the maize isolates that were absent from the dairy isolates. Many of the putative genes received an ‘unknown’ annotation, which led to the domain analysis of the 135 homologs. This revealed genes with molecular functions needed for BNF, including mucilage carbohydrate catabolism, glycan-mediated host adhesion, iron/siderophore utilization, and oxidation/reduction control. This is the first report of this pathway in this organism to underpin BNF. Consequently, we proposed a model needed for BNF in lactococci that plausibly accounts for BNF in the absence of the nif operon in this organism.

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.

A Genome-Based Model to Predict the Virulence of Pseudomonas aeruginosa Isolates

Variation in the genome of Pseudomonas aeruginosa, an important pathogen, can have dramatic impacts on the bacterium's ability to cause disease. We therefore asked whether it was possible to predict the virulence of P. aeruginosa isolates based on their genomic content. We applied a machine learning approach to a genetically and phenotypically diverse collection of 115 clinical P. aeruginosa isolates using genomic information and corresponding virulence phenotypes in a mouse model of bacteremia. We defined the accessory genome of these isolates through the presence or absence of accessory genomic elements (AGEs), sequences present in some strains but not others. Machine learning models trained using AGEs were predictive of virulence, with a mean nested cross-validation accuracy of 75% using the random forest algorithm. However, individual AGEs did not have a large influence on the algorithm's performance, suggesting instead that virulence predictions are derived from a diffuse genomic signature. These results were validated with an independent test set of 25 P. aeruginosa isolates whose virulence was predicted with 72% accuracy. Machine learning models trained using core genome single-nucleotide variants and whole-genome k-mers also predicted virulence. Our findings are a proof of concept for the use of bacterial genomes to predict pathogenicity in P. aeruginosa and highlight the potential of this approach for predicting patient outcomes.IMPORTANCE Pseudomonas aeruginosa is a clinically important Gram-negative opportunistic pathogen. P. aeruginosa shows a large degree of genomic heterogeneity both through variation in sequences found throughout the species (core genome) and through the presence or absence of sequences in different isolates (accessory genome). P. aeruginosa isolates also differ markedly in their ability to cause disease. In this study, we used machine learning to predict the virulence level of P. aeruginosa isolates in a mouse bacteremia model based on genomic content. We show that both the accessory and core genomes are predictive of virulence. This study provides a machine learning framework to investigate relationships between bacterial genomes and complex phenotypes such as virulence.

Assessment of vector-host-pathogen relationships using data mining and machine learning

Infectious diseases, including vector-borne diseases transmitted by arthropods, are a leading cause of morbidity and mortality worldwide. In the era of big data, addressing broad-scale, fundamental questions regarding the complex dynamics of these diseases will increasingly require the integration of diverse datasets to produce new biological knowledge. This review provides a current snapshot of the systematic assessment of the relationships between microbial pathogens, arthropod vectors and mammalian hosts using data mining and machine learning. We employ PRISMA to identify 32 key papers relevant to this topic. Our analysis shows an increasing use of data mining and machine learning tasks and techniques, including prediction, classification, clustering, association rules mining, and deep learning, over the last decade. However, it also reveals a number of critical challenges in applying these to the study of vector-host-pathogen interactions at various systems biology levels. Here, relevant studies, current limitations and future directions are discussed. Furthermore, the quality of data in relevant papers was assessed using the FAIR (Findable, Accessible, Interoperable, Reusable) compliance criteria to evaluate and encourage reproducibility and shareability of research outcomes. Although shortcomings in their application remain, data mining and machine learning have significant potential to break new ground in understanding fundamental aspects of vector-host-pathogen relationships and their application in this field should be encouraged. In particular, while predictive modeling, feature engineering and supervised machine learning are already being used in the field, other data mining and machine learning methods such as deep learning and association rules analysis lag behind and should be implemented in combination with established methods to accelerate hypothesis and knowledge generation in the domain.

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).

A Novel Approach for Biofilm Detection Based on a Convolutional Neural Network

Rhinology studies anatomy, physiology and diseases affecting the nasal region: one of the most modern techniques to diagnose these diseases is nasal cytology or rhinocytology, which involves analyzing the cells contained in the nasal mucosa under a microscope and researching of other elements such as bacteria, to suspect a pathology. During the microscopic observation, bacteria can be detected in the form of biofilm, that is, a bacterial colony surrounded by an organic extracellular matrix, with a protective function, made of polysaccharides. In the field of nasal cytology, the presence of biofilm in microscopic samples denotes the presence of an infection. In this paper, we describe the design and testing of interesting diagnostic support, for the automatic detection of biofilm, based on a convolutional neural network (CNN). To demonstrate the reliability of the system, alternative solutions based on isolation forest and deep random forest techniques were also tested. Texture analysis is used, with Haralick feature extraction and dominant color. The CNN-based biofilm detection system shows an accuracy of about 98%, an average accuracy of about 100% on the test set and about 99% on the validation set. The CNN-based system designed in this study is confirmed as the most reliable among the best automatic image recognition technologies, in the specific context of this study. The developed system allows the specialist to obtain a rapid and accurate identification of the biofilm in the slide images.

Genetic Markers in S. Paratyphi C Reveal Primary Adaptation to Pigs

Salmonella enterica with the identical antigenic formula 6,7:c:1,5 can be differentiated biochemically and by disease syndrome. One grouping, Salmonella Paratyphi C, is currently considered a typhoidal serovar, responsible for enteric fever in humans. The human-restricted typhoidal serovars (S. Typhi and Paratyphi A, B and C) typically display high levels of genome degradation and are cited as an example of convergent evolution for host adaptation in humans. However, S. Paratyphi C presents a different clinical picture to S. Typhi/Paratyphi A, in a patient group with predisposition, raising the possibility that its natural history is different, and that infection is invasive salmonellosis rather than enteric fever. Using whole genome sequencing and metabolic pathway analysis, we compared the genomes of 17 S. Paratyphi C strains to other members of the 6,7:c:1,5 group and to two typhoidal serovars: S. Typhi and Paratyphi A. The genome degradation observed in S. Paratyphi C was much lower than S. Typhi/Paratyphi A, but similar to the other 6,7:c:1,5 strains. Genomic and metabolic comparisons revealed little to no overlap between S. Paratyphi C and the other typhoidal serovars, arguing against convergent evolution and instead providing evidence of a primary adaptation to pigs in accordance with the 6,7:c:1.5 strains.