- Volume 7, Issue 4, 2021
Volume 7, Issue 4, 2021
- Research Articles
-
- Genomic Methodologies
-
-
Application of a strain-level shotgun metagenomics approach on food samples: resolution of the source of a Salmonella food-borne outbreak
Food-borne outbreak investigation currently relies on the time-consuming and challenging bacterial isolation from food, to be able to link food-derived strains to more easily obtained isolates from infected people. When no food isolate can be obtained, the source of the outbreak cannot be unambiguously determined. Shotgun metagenomics approaches applied to the food samples could circumvent this need for isolation from the suspected source, but require downstream strain-level data analysis to be able to accurately link to the human isolate. Until now, this approach has not yet been applied outside research settings to analyse real food-borne outbreak samples. In September 2019, a Salmonella outbreak occurred in a hotel school in Bruges, Belgium, affecting over 200 students and teachers. Following standard procedures, the Belgian National Reference Center for human salmonellosis and the National Reference Laboratory for Salmonella in food and feed used conventional analysis based on isolation, serotyping and MLVA (multilocus variable number tandem repeat analysis) comparison, followed by whole-genome sequencing, to confirm the source of the contamination over 2 weeks after receipt of the sample, which was freshly prepared tartar sauce in a meal cooked at the school. Our team used this outbreak as a case study to deliver a proof of concept for a short-read strain-level shotgun metagenomics approach for source tracking. We received two suspect food samples: the full meal and some freshly made tartar sauce served with this meal, requiring the use of raw eggs. After analysis, we could prove, without isolation, that Salmonella was present in both samples, and we obtained an inferred genome of a Salmonella enterica subsp. enterica serovar Enteritidis that could be linked back to the human isolates of the outbreak in a phylogenetic tree. These metagenomics-derived outbreak strains were separated from sporadic cases as well as from another outbreak circulating in Europe at the same time period. This is, to our knowledge, the first Salmonella food-borne outbreak investigation uniquely linking the food source using a metagenomics approach and this in a fast time frame.
-
-
-
Rapid nanopore-based DNA sequencing protocol of antibiotic-resistant bacteria for use in surveillance and outbreak investigation
Outbreak investigations are essential to control and prevent the dissemination of pathogens. This study developed and validated a complete analysis protocol for faster and more accurate surveillance and outbreak investigations of antibiotic-resistant microbes based on Oxford Nanopore Technologies (ONT) DNA whole-genome sequencing. The protocol was developed using 42 methicillin-resistant Staphylococcus aureus (MRSA) isolates identified from former well-characterized outbreaks. The validation of the protocol was performed using Illumina technology (MiSeq, Illumina). Additionally, a real-time outbreak investigation of six clinical S. aureus isolates was conducted to test the ONT-based protocol. The suggested protocol includes: (1) a 20 h sequencing run; (2) identification of the sequence type (ST); (3) de novo genome assembly; (4) polishing of the draft genomes; and (5) phylogenetic analysis based on SNPs. After the sequencing run, it was possible to identify the ST in 2 h (20 min per isolate). Assemblies were achieved after 4 h (40 min per isolate) while the polishing was carried out in 7 min per isolate (42 min in total). The phylogenetic analysis took 0.6 h to confirm an outbreak. Overall, the developed protocol was able to at least discard an outbreak in 27 h (mean) after the bacterial identification and less than 33 h to confirm it. All these estimated times were calculated considering the average time for six MRSA isolates per sequencing run. During the real-time S. aureus outbreak investigation, the protocol was able to identify two outbreaks in less than 31 h. The suggested protocol enables identification of outbreaks in early stages using a portable and low-cost device along with a streamlined downstream analysis, therefore having the potential to be incorporated in routine surveillance analysis workflows. In addition, further analysis may include identification of virulence and antibiotic resistance genes for improved pathogen characterization.
-
- Functional Genomics and Microbe–Niche Interactions
-
-
Niche-specific adaptation of Lactobacillus helveticus strains isolated from malt whisky and dairy fermentations
More LessLactobacillus helveticus is a well characterized lactobacillus for dairy fermentations that is also found in malt whisky fermentations. The two environments contain considerable differences related to microbial growth, including the presence of different growth inhibitors and nutrients. The present study characterized L. helveticus strains originating from dairy fermentations (called milk strains hereafter) and malt whisky fermentations (called whisky strains hereafter) by in vitro phenotypic tests and comparative genomics. The whisky strains can tolerate ethanol more than the milk strains, whereas the milk strains can tolerate lysozyme and lactoferrin more than the whisky strains. Several plant-origin carbohydrates, including cellobiose, maltose, sucrose, fructooligosaccharide and salicin, were generally metabolized only by the whisky strains, whereas milk-derived carbohydrates, i.e. lactose and galactose, were metabolized only by the milk strains. Milk fermentation properties also distinguished the two groups. The general genomic characteristics, including genomic size, number of coding sequences and average nucleotide identity values, differentiated the two groups. The observed differences in carbohydrate metabolic properties between the two groups correlated with the presence of intact specific enzymes in glycoside hydrolase (GH) families GH1, GH4, GH13, GH32 and GH65. Several GHs in the milk strains were inactive due to the presence of stop codon(s) in genes encoding the GHs, and the inactivation patterns of the genes encoding specific enzymes assigned to GH1 in the milk strains suggested a possible diversification manner of L. helveticus strains. The present study has demonstrated how L. helveticus strains have adapted to their habitats.
-
- Microbial Communities
-
-
Shotgun metagenomics reveals distinct functional diversity and metabolic capabilities between 12 000-year-old permafrost and active layers on Muot da Barba Peider (Swiss Alps)
The warming-induced thawing of permafrost promotes microbial activity, often resulting in enhanced greenhouse gas emissions. The ability of permafrost microorganisms to survive the in situ sub-zero temperatures, their energetic strategies and their metabolic versatility in using soil organic materials determine their growth and functionality upon thawing. Hence, functional characterization of the permafrost microbiome, particularly in the underexplored mid-latitudinal alpine regions, is a crucial first step in predicting its responses to the changing climate, and the consequences for soil–climate feedbacks. In this study, for the first time, the functional potential and metabolic capabilities of a temperate mountain permafrost microbiome from central Europe has been analysed using shotgun metagenomics. Permafrost and active layers from the summit of Muot da Barba Peider (MBP) [Swiss Alps, 2979 m above sea level (a.s.l.)] revealed a strikingly high functional diversity in the permafrost (north-facing soils at a depth of 160 cm). Permafrost metagenomes were enriched in stress-response genes (e.g. cold-shock genes, chaperones), as well as in genes involved in cell defence and competition (e.g. antiviral proteins, antibiotics, motility, nutrient-uptake ABC transporters), compared with active-layer metagenomes. Permafrost also showed a higher potential for the synthesis of carbohydrate-active enzymes, and an overrepresentation of genes involved in fermentation, carbon fixation, denitrification and nitrogen reduction reactions. Collectively, these findings demonstrate the potential capabilities of permafrost microorganisms to thrive in cold and oligotrophic conditions, and highlight their metabolic versatility in carbon and nitrogen cycling. Our study provides a first insight into the high functional gene diversity of the central European mountain permafrost microbiome. Our findings extend our understanding of the microbial ecology of permafrost and represent a baseline for future investigations comparing the functional profiles of permafrost microbial communities at different latitudes.
-
- Pathogens and Epidemiology
-
-
Detection, spread and phylogeny of meticillin-resistant Staphylococcus aureus sequence type 45 in Taiwan
More LessMeticillin-resistant Staphylococcus aureus (MRSA) sequence type (ST) 45 was reported in the literature to have been first identified in 2006 in Taiwan. The present study was carried out to explore and trace the emergence, transmission and evolutional dynamics of MRSA ST45 in Taiwan. We identified MRSA ST45 isolates retrospectively from two collections of MRSA isolates, namely TSAR (Taiwan Surveillance of Antimicrobial Resistance) surveys and the CGMH (Chang Gung Memorial Hospital)-based laboratory collection. Representative ST45 isolates were selected for whole-genome sequencing (WGS) analysis. A total of 9554 MRSA isolates was included in this study. Among the 3766 MRSA isolates biennially collected from TSAR surveys between 1998 and 2014, ST45 accounted for 133 (3.53 %) MRSA isolates, was first identified in 2004, and the prevalence rate peaked in 2010 (up to 10.77 %). Among the 5788 MRSA isolates collected between 1995 and 2017 by the CGMH-based laboratory, 257 isolates (4.44 %) were characterized as ST45, with most identified from nursing homes since 2012. Of the 75 isolates randomly selected for WGS, two clades were identified. The major clade, clade II, comprised 63 isolates and was phylogenetically relatively close to those isolates identified from Singapore. All but one of the isolates in clade I, the minor clade, were identified from non-Taiwanese people, mostly from newly recruited foreign workers in 2017, and were phylogenetically relatively close to one isolate from the USA (CA-347). Conclusively, the emergence of MRSA ST45 strain in Taiwan can be traced back to 2004 and the strain is connected to South-East Asian countries. Since its emergence, transmission and spread of MRSA ST45 has occurred in Taiwan.
-
-
-
Origin and evolution of nonulosonic acid synthases and their relationship with bacterial pathogenicity revealed by a large-scale phylogenetic analysis
More LessNonulosonic acids (NulOs) are a group of nine-carbon monosaccharides with different functions in nature. N-acetylneuraminic acid (Neu5Ac) is the most common NulO. It covers the membrane surface of all human cells and is a central molecule in the process of self-recognition via SIGLECS receptors. Some pathogenic bacteria escape the immune system by copying the sialylation of the host cell membrane. Neu5Ac production in these bacteria is catalysed by the enzyme NeuB. Some bacteria can also produce other NulOs named pseudaminic and legionaminic acids, through the NeuB homologues PseI and LegI, respectively. In Opisthokonta eukaryotes, the biosynthesis of Neu5Ac is catalysed by the enzyme NanS. In this study, we used publicly available data of sequences of NulOs synthases to investigate its distribution within the three domains of life and its relationship with pathogenic bacteria. We mined the KEGG database and found 425 NeuB sequences. Most NeuB sequences (58.74 %) from the KEGG orthology database were classified as from environmental bacteria; however, sequences from pathogenic bacteria showed higher conservation and prevalence of a specific domain named SAF. Using the HMM profile we identified 13 941 NulO synthase sequences in UniProt. Phylogenetic analysis of these sequences showed that the synthases were divided into three main groups that can be related to the lifestyle of these bacteria: (I) predominantly environmental, (II) intermediate and (III) predominantly pathogenic. NeuB was widely distributed in the groups. However, LegI and PseI were more concentrated in groups II and III, respectively. We also found that PseI appeared later in the evolutionary process, derived from NeuB. We use this same methodology to retrieve sialic acid synthase sequences from Archaea and Eukarya. A large-scale phylogenetic analysis showed that while the Archaea sequences are spread across the tree, the eukaryotic NanS sequences were grouped in a specific branch in group II. None of the bacterial NanS sequences grouped with the eukaryotic branch. The analysis of conserved residues showed that the synthases of Archaea and Eukarya present a mutation in one of the three catalytic residues, an E134D change, related to a Neisseria meningitidis reference sequence. We also found that the conservation profile is higher between NeuB of pathogenic bacteria and NanS of eukaryotes than between NeuB of environmental bacteria and NanS of eukaryotes. Our large-scale analysis brings new perspectives on the evolution of NulOs synthases, suggesting their presence in the last common universal ancestor.
-
- Evolution and Responses to Interventions
-
-
Genome-wide analysis in Escherichia coli unravels a high level of genetic homoplasy associated with cefotaxime resistance
Cefotaxime (CTX) is a third-generation cephalosporin (3GC) commonly used to treat infections caused by Escherichia coli . Two genetic mechanisms have been associated with 3GC resistance in E. coli . The first is the conjugative transfer of a plasmid harbouring antibiotic-resistance genes. The second is the introduction of mutations in the promoter region of the ampC β-lactamase gene that cause chromosome-encoded β-lactamase hyperproduction. A wide variety of promoter mutations related to AmpC hyperproduction have been described. However, their link to CTX resistance has not been reported. We recultured 172 cefoxitin-resistant E. coli isolates with known CTX minimum inhibitory concentrations and performed genome-wide analysis of homoplastic mutations associated with CTX resistance by comparing Illumina whole-genome sequencing data of all isolates to a PacBio sequenced reference chromosome. We mapped the mutations on the reference chromosome and determined their occurrence in the phylogeny, revealing extreme homoplasy at the −42 position of the ampC promoter. The 24 occurrences of a T at the −42 position rather than the wild-type C, resulted from 18 independent C>T mutations in five phylogroups. The −42 C>T mutation was only observed in E. coli lacking a plasmid-encoded ampC gene. The association of the −42 C>T mutation with CTX resistance was confirmed to be significant (false discovery rate <0.05). To conclude, genome-wide analysis of homoplasy in combination with CTX resistance identifies the −42 C>T mutation of the ampC promotor as significantly associated with CTX resistance and underlines the role of recurrent mutations in the spread of antibiotic resistance.
-
-
-
Evaluation of bat adenoviruses suggests co-evolution and host roosting behaviour as drivers for diversity
Nkom F. Ntumvi, Joseph Le Doux Diffo, Ubald Tamoufe, Valantine Ngum Ndze, Jean-Michel Takuo, Moctar M. M. Mouiche, Julius Nwobegahay, Matthew LeBreton, Amethyst Gillis, Anne W. Rimoin, Bradley S. Schneider, Corina Monagin, David J. McIver, Damien O. Joly, Nathan D. Wolfe, Edward M. Rubin and Christian E. LangeAdenoviruses (AdVs) are diverse pathogens of humans and animals, with several dozen bat AdVs already identified. Considering that over 100 human AdVs are known, and the huge diversity of bat species, many bat AdVs likely remain undiscovered. To learn more about AdV prevalence, diversity and evolution, we sampled and tested bats in Cameroon using several PCR assays for viral and host DNA. AdV DNA was detected in 14 % of the 671 sampled animals belonging to 37 different bat species. There was a correlation between species roosting in larger groups and AdV DNA detection. The detected AdV DNA belonged to between 28 and 44 different, mostly previously unknown, mastadenovirus species. The novel isolates are phylogenetically diverse and while some cluster with known viruses, others appear to form divergent new clusters. The phylogenetic tree of novel and previously known bat AdVs does not mirror that of the various host species, but does contain structures consistent with a degree of virus–host co-evolution. Given that closely related isolates were found in different host species, it seems likely that at least some bat AdVs have jumped species barriers, probably in the more recent past; however, the tree is also consistent with such events having taken place throughout bat AdV evolution. AdV diversity was highest in bat species roosting in large groups. The study significantly increased the diversity of AdVs known to be harboured by bats, and suggests that host behaviours, such as roosting size, may be what limits some AdVs to one species rather than an inability of AdVs to infect other related hosts.
-
- Short Communications
-
- Pathogens and Epidemiology
-
-
Genomics reveals the novel species placement of industrial contaminant isolates incorrectly identified as Burkholderia lata
The Burkholderia cepacia complex (Bcc) is a closely related group of bacteria, composed of at least 20 different species, the accurate identification of which is essential in the context of infectious diseases. In industry, they can contaminate non-food products, including home and personal care products and cosmetics. The Bcc are problematic contaminants due to their ubiquitous presence and intrinsic antimicrobial resistance, which enables them to occasionally overcome preservation systems in non-sterile products. Burkholderia lata and Burkholderia contaminans are amongst the Bcc bacteria encountered most frequently as industrial contaminants, but their identification is not straightforward. Both species were historically established as a part of a group known collectively as taxon K, based upon analysis of the recA gene and multilocus sequence typing (MLST). Here, we deploy a straightforward genomics-based workflow for accurate Bcc classification using average nucleotide identity (ANI) and core-gene analysis. The workflow was used to examine a panel of 23 Burkholderia taxon K industrial strains, which, based on MLST, comprised 13 B. lata, 4 B. contaminans and 6 unclassified Bcc strains. Our genomic identification showed that the B. contaminans strains retained their classification, whilst the remaining strains were reclassified as Burkholderia aenigmatica sp. nov. Incorrect taxonomic identification of industrial contaminants is a problematic issue. Application and testing of our genomic workflow allowed the correct classification of 23 Bcc industrial strains, and also indicated that B. aenigmatica sp. nov. may have greater importance than B. lata as a contaminant species. Our study illustrates how the non-food manufacturing industry can harness whole-genome sequencing to better understand antimicrobial-resistant bacteria affecting their products.
-
- Methods
-
- Genomic Methodologies
-
-
Detection of plasmid contigs in draft genome assemblies using customized Kraken databases
More LessPlasmids play an important role in bacterial evolution and mediate horizontal transfer of genes including virulence and antimicrobial resistance genes. Although short-read sequencing technologies have enabled large-scale bacterial genomics, the resulting draft genome assemblies are often fragmented into hundreds of discrete contigs. Several tools and approaches have been developed to identify plasmid sequences in such assemblies, but require trade-off between sensitivity and specificity. Here we propose using the Kraken classifier, together with a custom Kraken database comprising known chromosomal and plasmid sequences of Klebsiella pneumoniae species complex (KpSC), to identify plasmid-derived contigs in draft assemblies. We assessed performance using Illumina-based draft genome assemblies for 82 KpSC isolates, for which complete genomes were available to supply ground truth. When benchmarked against five other classifiers (Centrifuge, RFPlasmid, mlplasmids, PlaScope and Platon), Kraken showed balanced performance in terms of overall sensitivity and specificity (90.8 and 99.4 %, respectively, for contig count; 96.5 and >99.9 %, respectively, for cumulative contig length), and the highest accuracy (96.8% vs 91.8-96.6% for contig count; 99.8% vs 99.0-99.7 % for cumulative contig length), and F1-score (94.5 % vs 84.5-94.1 %, for contig count; 98.0 % vs 88.9-96.7 % for cumulative contig length). Kraken also achieved consistent performance across our genome collection. Furthermore, we demonstrate that expanding the Kraken database with additional known chromosomal and plasmid sequences can further improve classification performance. Although we have focused here on the KpSC, this methodology could easily be applied to other species with a sufficient number of completed genomes.
-
- Research Articles
-
- Pathogens and Epidemiology
-
-
Complete, closed and curated genome sequences of Photobacterium damselae subsp. piscicida isolates from Australia indicate mobilome-driven localized evolution and novel pathogenicity determinants
Despite the recent advances in sequencing technologies, the complete assembly of multi-chromosome genomes of the Vibrionaceae , often containing several plasmids, remains challenging. Using a combination of Oxford Nanopore MinION long reads and short Illumina reads, we fully sequenced, closed and curated the genomes of two strains of a primary aquatic pathogen Photobacterium damselae subsp. piscicida isolated in Australia. These are also the first genome sequences of P. damselae subsp. piscicida isolated in Oceania and, to our knowledge, in the Southern hemisphere. We also investigated the phylogenetic relationships between Australian and overseas isolates, revealing that Australian P. damselae subsp. piscicida are more closely related to the Asian and American strains rather than to the European ones. We investigated the mobilome and present new evidence showing that a host specialization process and progressive adaptive evolution to fish are ongoing in P. damselae subsp. piscicida , and are largely mediated by transposable elements, predominantly in chromosome 2, and by plasmids. Finally, we identified two novel potential virulence determinants in P. damselae subsp. piscicida – a chorismate mutase gene, which is ubiquitously retained and co-localized with the AIP56 apoptogenic toxin-encoding gene on the pPHDP10 plasmid, and transfer-messenger RNA gene ssrA located on the main chromosome, homologous to a critical-to-virulence determinant in Yersinia pseudotuberculosis . Our study describes, to our knowledge, the only fully closed and manually curated genomes of P. damselae subsp. piscicida available to date, offering new insights into this important fish pathogen and its evolution.
-