PGAdb-builder: A web service tool for creating pan-genome allele database for molecular fine typing

Antimicrobial resistance and genetic relatedness among Campylobacter coli and Campylobacter jejuni from humans and retail chicken meat in Taiwan

OBJECTIVES

Campylobacter is a significant zoonotic pathogen primarily transmitted through poultry. Our study aimed to assess antimicrobial resistance and genetic relationships among Campylobacter isolates from retail chicken meat and humans in Taiwan.

METHODS

Campylobacter isolates were analysed using whole-genome sequencing to investigate their antimicrobial resistance, genetic determinants of resistance, and genotypes.

RESULTS

Campylobacter coli and Campylobacter jejuni accounted for 44.9% and 55.1% of chicken meat isolates, and 11.4% and 88.6% of human isolates, respectively. C. coli displayed significantly higher resistance levels. Furthermore, isolates from chicken meat exhibited higher levels of resistance to most tested antimicrobials compared to isolates from humans. Multidrug resistance was observed in 96.3% of C. coli and 43.3% of C. jejuni isolates from chicken meat and 80.6% of C. coli and 15.8% of C. jejuni isolates from humans. Macrolide resistance was observed in 85.5% of C. coli isolates, primarily attributed to the erm(B) rather than the A2075G mutation in 23S rRNA. Among the 511 genomes, we identified 133 conventional MLST sequence types, indicating significant diversity among Campylobacter strains. Notably, hierarchical Core-genome multilocus sequence typing clustering, including HC0, HC5, and HC10, revealed a significant proportion of closely related isolates from chicken meat and humans.

CONCLUSIONS

Our research highlights significant associations in antimicrobial resistance and genetic relatedness between Campylobacter isolates from chicken meat and humans in Taiwan. The genetic analysis data suggest that campylobacteriosis outbreaks may occur more frequently in Taiwan than previously assumed. Our study emphasizes the need for strategies to control multidrug-resistant strains and enhance outbreak prevention.

Pan-genomics: Insight into the Functional Genome, Applications, Advancements, and Challenges

A pan-genome is a compilation of the common and unique genomes found in a given species. It incorporates the genetic information from all of the genomes sampled, producing a big and diverse set of genetic material. Pan-genomic analysis has various advantages over typical genomics research. It creates a vast and varied spectrum of genetic material by combining the genetic data from all the sampled genomes. Comparing pan-genomics analysis to conventional genomic research, there are a number of benefits. Although the most recent era of pan-genomic studies has used cutting-edge sequencing technology to shed fresh light on biological variety and improvement, the potential uses of pan-genomics in improvement have not yet been fully realized. Pan-genome research in various organisms has demonstrated that missing genetic components and the detection of significant Structural Variants (SVs) can be investigated using pan-genomic methods. Many individual-specific sequences have been linked to biological adaptability, phenotypic, and key economic attributes. This study aims to focus on how pangenome analysis uncovers genetic differences in various organisms, including human, and their effects on phenotypes, as well as how this might help us comprehend the diversity of species. The review also concentrated on potential problems and the prospects for future pangenome research.

Insights on the genomic diversity, virulence and resistance profile of a Campylobacter jejuni strain isolated from a hospitalized patient in Brazil

Campylobacteriosis is currently recognized as one of the major causes of foodborne bacterial diseases worldwide. In Brazil, there is insufficient data to estimate the impact of Campylobacter in public health. The aim of this present study was to characterize a C. jejuni CJ-HBSJRP strain isolated from a hospitalized patient in Brazil by its ability to invade human Caco-2 epithelial cells, to survive in U937 human macrophages, and to assess its phenotypic antimicrobial resistance profile. In addition, prophages, virulence and antimicrobial resistance genes were search using whole-genome sequencing data. The genetic relatedness was evaluated by MLST and cgMLST analysis by comparison with 29 other C. jejuni genomes isolated from several countries. The CJ-HBSJRP strain showed an invasion percentage of 50% in Caco-2 polarized cells, 37.5% of survivability in U937 cells and was phenotypically resistant to ampicillin, ciprofloxacin and nalidixic acid. A total of 94 virulence genes related to adherence, biofilm, chemotaxis, immune modulation, invasion process, metabolism, motility and toxin were detected. The resistance genes blaOXA-605 (blaOXA-61), cmeB and mutations in the QRDR region of gyrA were also found and none prophages were detected. The MLST analysis showed 23 different STs among the strains studied. Regarding cgMLST analysis, the CJ-HBSJRP strain was genetically distinct and did not group closely to any other isolate. The results obtained reinforce the pathogenic potential of the CJHBSJRP strain and highlighted the need for more careful attention to Campylobacter spp. infections in Brazil since this pathogen has been the most commonly reported zoonosis in several countries worldwide.

Comparative genomics reveals high genetic similarity among strains of Salmonella enterica serovar Infantis isolated from multiple sources in Brazil

Background

Salmonella enterica serovar Infantis (Salmonella Infantis) is a zoonotic, ubiquitous and foodborne pathogen of worldwide distribution. Despite Brazil's relevance as a major meat exporter, few studies were conducted to characterize strains of this serovar by genomic analyses in this country. Therefore, this study aimed to assess the diversity of 80 Salmonella Infantis strains isolated from veterinary, food and human sources in Brazil between 2013 and 2018 by comparative genomic analyses. Additional genomes of non-Brazilian countries (n = 18) were included for comparison purposes in some analyses.

Methods

Analyses of whole-genome multi-locus sequence typing (wgMLST), using PGAdb-builder, and of fragmented genomes, using Gegenees, were conducted to compare the 80 Brazilian strains to the 18 non-Brazilian genomes. Pangenome analyses and calculations were performed for all Salmonella Infantis genomes analyzed. The presence of prophages was determined using PHASTER for the 80 Brazilian strains. The genome plasticity using BLAST Ring Image Generator (BRIG) and gene synteny using Mauve were evaluated for 20 selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. Unique orthologous protein clusters were searched in ten selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries.

Results

wgMLST and Gegenees showed a high genomic similarity among some Brazilian Salmonella Infantis genomes, and also the correlation of some clusters with non-Brazilian genomes. Gegenees also showed an overall similarity >91% among all Salmonella Infantis genomes. Pangenome calculations revealed an open pangenome for all Salmonella Infantis subsets analyzed and a high gene content in the core genomes. Fifteen types of prophages were detected among 97.5% of the Brazilian strains. BRIG and Mauve demonstrated a high structural similarity among the Brazilian and non-Brazilian isolates. Unique orthologous protein clusters related to biological processes, molecular functions, and cellular components were detected among Brazilian and non-Brazilian genomes.

Conclusion

The results presented using different genomic approaches emphasized the significant genomic similarity among Brazilian Salmonella Infantis genomes analyzed, suggesting wide distribution of closely related genotypes among diverse sources in Brazil. The data generated contributed to novel information regarding the genomic diversity of Brazilian and non-Brazilian Salmonella Infantis in comparison. The different genetically related subtypes of Salmonella Infantis from Brazil can either occur exclusively within the country, or also in other countries, suggesting that some exportation of the Brazilian genotypes may have already occurred.

High-Throughput Short Sequence Typing Schemes for Pseudomonas aeruginosa and Stenotrophomonas maltophilia Pure Culture and Environmental DNA

Molecular typing techniques are utilized to determine genetic similarities between bacterial isolates. However, the use of environmental DNA profiling to assess epidemiologic links between patients and their environment has not been fully explored. This work reports the development and validation of two high-throughput short sequence typing (HiSST) schemes targeting the opportunistic pathogens Pseudomonas aeruginosa and Stenotrophomonas maltophilia, along with a modified SM2I selective medium for the specific isolation of S. maltophilia. These HiSST schemes are based on four discriminative loci for each species and demonstrate high discriminating power, comparable to pairwise whole-genome comparisons. Each scheme includes species-specific PCR primers for precise differentiation from closely related taxa, without the need for upstream culture-dependent methods. For example, the primers targeting the bvgS locus make it possible to distinguish P. aeruginosa from the very closely related Pseudomonas paraeruginosa sp. nov. The selected loci included in the schemes are adapted to massive parallel amplicon sequencing technology. An R-based script implemented in the DADA2 pipeline was assembled to facilitate HiSST analyses for efficient and accurate genotyping of P. aeruginosa and S. maltophilia. We demonstrate the performance of both schemes through in silico validations, assessments against reference culture collections, and a case study involving environmental samples.

ProbioMinServer: an integrated platform for assessing the safety and functional properties of potential probiotic strains

Motivation

ProbioMinServer is a platform designed to help researchers access information on probiotics regarding a wide variety of characteristics, such as safety (e.g. antimicrobial resistance, virulence, pathogenic, plasmid, and prophage genes) and functionality (e.g. functional classes, carbohydrate-active enzyme, and metabolite gene cluster profile). Because probiotics are functional foods, their safety and functionality are a crucial part of health care. Genomics has become a crucial methodology for investigating the safety and functionality of probiotics in food and feed. This shift is primarily attributed to the growing affordability of next-generation sequencing technologies. However, no integrated platform is available for simultaneously evaluating probiotic strain safety, investigating probiotic functionality, and identifying known phylogenetically related strains.

Results

Thus, we constructed a new platform, ProbioMinServer, which incorporates these functions. ProbioMinServer accepts whole-genome sequence files in the FASTA format. If the query genome belongs to the 25 common probiotic species collected in our database, the server performs a database search and analyzes the core-genome multilocus sequence typing. Front-end applications were implemented in JavaScript with a bootstrap framework, and back-end programs were implemented using PHP, Perl, and Python. ProbioMinServer can help researchers quickly and easily retrieve information on the safety and functionality of various probiotics.

Availability and implementation

The platform is available at https://probiomindb.imst.nsysu.edu.tw.

Comparative genomics of Bordetella pertussis and prediction of new vaccines and drug targets

Pertussis is a highly contagious respiratory disease caused by Bordetella pertussis, a Gram-negative bacterium described over a century ago. Despite broad vaccine coverage and treatment options, the disease is remerging as a public health problem especially in infants and older children. Recent data indicate re-emergence of the disease is related to bacterial resistance to immune defences and decreased vaccine effectiveness, which obviously suggests the need of new effective vaccines and drugs. In an attempt to contribute with solutions to this great challenge, bioinformatics tools were used to genetically comprehend the species of these bacteria and predict new vaccines and drug targets. In fact, approaches were used to analysis genomic plasticity, gene synteny and species similarities between the 20 genomes of Bordetella pertussis already available. Furthermore, it was conducted reverse vaccinology and docking analysis to identify proteins with potential to become vaccine and drug targets, respectively. The analyses showed the 20 genomes belongs to a homogeneous group that has preserved most of the genes over time. Besides that, were found genomics islands and good proteins to be candidates for vaccine and drugs. Taken together, these results suggests new possibilities that may be useful to develop new vaccines and drugs that will help the prevention and treatment strategies of pertussis disease caused by these Bordetella strains. Communicated by Ramaswamy H. Sarma.

5NosoAE: a web server for nosocomial bacterial antibiogram investigation and epidemiology survey

5NosoAE is a webserver that can be used for nosocomial bacterial analysis including the identification of similar strains based on antimicrobial resistance profiles (antibiogram) and the spatiotemporal distribution visualization and phylogenetic analysis of identified strains with similar antibiograms. The extensive use of antibiotics has caused many pathogenic bacteria to develop multiple drug resistance, resulting in clinical infection treatment challenges and posing a major threat to global public health. Relevant studies have investigated the key determinants of antimicrobial resistance in the whole-genome sequence of bacteria. However, a web server is currently not available for performing large-scale strain searches according to antimicrobial resistance profiles and visualizing epidemiological information including the spatiotemporal distribution, antibiogram heatmap, and phylogeny of identified strains. Here, we implemented these functions in the new server, referred to as 5NosoAE. This server accepts the genome sequence file in the FASTA format of five nosocomial bacteria, namely Acinetobacter baumannii, Pseudomonas aeruginosa, Klebsiella pneumoniae, Enterococcus faecium and Staphylococcus aureus for query. All visualizations are implemented in JavaScript and PHP. This server will be useful for physicians and epidemiologists involved in research on infectious disease. The 5NosoAE platform is available at https://nosoae.imst.nsysu.edu.tw.

LmTraceMap: A Listeria monocytogenes fast-tracing platform for global surveillance

Listeria monocytogenes can cause listeriosis, and people with hypoimmunity such as pregnant women, infants and fetuses are at high risk of invasive infection. Although the incidence of listeriosis is low, the fatality rate is high. Therefore, continual surveillance and rapid epidemiological investigation are crucial for addressing L. monocytogenes. Because of the popularity of next-generation sequencing, obtaining the whole-genome sequence of a bacterium is easy. Several genome-based typing methods are available, and core-genome multilocus sequence typing (cgMLST) is the most recognized methods. Using cgMLST typing to compare L. monocytogenes whole-genome sequences (WGS) with those obtained across distinct regions is beneficial. However, the concern is how to incorporate the powerful cgMLST method into investigations, such as by using source tracing. Herein, we present an easy-to-use web service called–LmTraceMap (http://lmtracemap.cgu.edu.tw/hua_map/test/upload.php; http://120.126.17.192/hua_map/test/upload.php) that can help public-health professionals rapidly trace closely related isolates worldwide and visually inspect them in search results on a world map with labeled epidemiological data. We expect the proposed service to improve the convenience of public health investigations.

Pangenomics in Microbial and Crop Research: Progress, Applications, and Perspectives

Advances in sequencing technologies and bioinformatics tools have fueled a renewed interest in whole genome sequencing efforts in many organisms. The growing availability of multiple genome sequences has advanced our understanding of the within-species diversity, in the form of a pangenome. Pangenomics has opened new avenues for future research such as allowing dissection of complex molecular mechanisms and increased confidence in genome mapping. To comprehensively capture the genetic diversity for improving plant performance, the pangenome concept is further extended from species to genus level by the inclusion of wild species, constituting a super-pangenome. Characterization of pangenome has implications for both basic and applied research. The concept of pangenome has transformed the way biological questions are addressed. From understanding evolution and adaptation to elucidating host–pathogen interactions, finding novel genes or breeding targets to aid crop improvement to design effective vaccines for human prophylaxis, the increasing availability of the pangenome has revolutionized several aspects of biological research. The future availability of high-resolution pangenomes based on reference-level near-complete genome assemblies would greatly improve our ability to address complex biological problems.

Characterization, Pathogenicity, Phylogeny, and Comparative Genomic Analysis of Pseudomonas tolaasii Strains Isolated from Various Mushrooms in China

Since 2016, devastating bacterial blotch affecting the fruiting bodies of Agaricus bisporus, Cordyceps militaris, Flammulina filiformis, and Pleurotus ostreatus in China has caused severe economic losses. We isolated 102 bacterial strains and characterized them polyphasically. We identified the causal agent as Pseudomonas tolaasii and confirmed the pathogenicity of the strains. A host range test further confirmed the pathogen's ability to infect multiple hosts. This is the first report in China of bacterial blotch in C. militaris caused by P. tolaasii. Whole-genome sequences were generated for three strains: Pt11 (6.48 Mb), Pt51 (6.63 Mb), and Pt53 (6.80 Mb), and pangenome analysis was performed with 13 other publicly accessible P. tolaasii genomes to determine their genetic diversity, virulence, antibiotic resistance, and mobile genetic elements. The pangenome of P. tolaasii is open, and many more gene families are likely to emerge with further genome sequencing. Multilocus sequence analysis using the sequences of four common housekeeping genes (glns, gyrB, rpoB, and rpoD) showed high genetic variability among the P. tolaasii strains, with 115 strains clustered into a monophyletic group. The P. tolaasii strains possess various genes for secretion systems, virulence factors, carbohydrate-active enzymes, toxins, secondary metabolites, and antimicrobial resistance genes that are associated with pathogenesis and adapted to different environments. The myriad of insertion sequences, integrons, prophages, and genome islands encoded in the strains may contribute to genome plasticity, virulence, and antibiotic resistance. These findings advance understanding of the determinants of virulence, which can be targeted for the effective control of bacterial blotch disease.

A High-Throughput Short Sequence Typing Scheme for Serratia marcescens Pure Culture and Environmental DNA

Molecular typing methods are used to characterize the relatedness between bacterial isolates involved in infections. These approaches rely mostly on discrete loci or whole-genome sequencing (WGS) analyses of pure cultures. On the other hand, their application to environmental DNA profiling to evaluate epidemiological relatedness among patients and environments has received less attention. We developed a specific, high-throughput short sequence typing (HiSST) method for the opportunistic human pathogen Serratia marcescens. Genes displaying the highest polymorphism were retrieved from the core genome of 60 S. marcescens strains. Bioinformatics analyses showed that use of only three loci (within bssA, gabR, and dhaM) distinguished strains with a high level of efficiency. This HiSST scheme was applied to an epidemiological survey of S. marcescens in a neonatal intensive care unit (NICU). In a first case study, a strain responsible for an outbreak in the NICU was found in a sink drain of this unit, by using HiSST scheme and confirmed by WGS. The HiSST scheme was also applied to environmental DNA extracted from sink-environment samples. Diversity of S. marcescens was modest, with 11, 6, and 4 different sequence types (ST) of gabR, bssA, and dhaM loci among 19 sink drains, respectively. Epidemiological relationships among sinks were inferred on the basis of pairwise comparisons of ST profiles. Further research aimed at relating ST distribution patterns to environmental features encompassing sink location, utilization, and microbial diversity is needed to improve the surveillance and management of opportunistic pathogens. IMPORTANCE Serratia marcescens is an important opportunistic human pathogen, often multidrug resistant and involved in outbreaks of nosocomial infections in neonatal intensive care units. Here, we propose a quick and user-friendly method to select the best typing scheme for nosocomial outbreaks in relating environmental and clinical sources. This method, named high-throughput short sequence typing (HiSST), allows to distinguish strains and to explore the diversity profile of nonculturable S. marcescens. The application of HiSST profile analysis for environmental DNA offers new possibilities to track opportunistic pathogens, identify their origin, and relate their distribution pattern with environmental features encompassing sink location, utilization, and microbial diversity. Adaptation of the method to other opportunistic pathogens is expected to improve knowledge regarding their ecology, which is of significant interest for epidemiological risk assessment and elaborate outbreak mitigation strategies.

A machine learning-based typing scheme refinement for Listeria monocytogenes core genome multilocus sequence typing with high discriminatory power for common source outbreak tracking

Background

As whole-genome sequencing for pathogen genomes becomes increasingly popular, the typing methods of gene-by-gene comparison, such as core genome multilocus sequence typing (cgMLST) and whole-genome multilocus sequence typing (wgMLST), are being routinely implemented in molecular epidemiology. However, some intrinsic problems remain. For example, genomic sequences with varying read depths, read lengths, and assemblers influence the genome assemblies, introducing error or missing alleles into the generated allelic profiles. These errors and missing alleles might create “specious discrepancy” among closely related isolates, thus making accurate epidemiological interpretation challenging. In addition, the rapid growth of the cgMLST allelic profile database can cause problems related to storage and maintenance as well as long query search times.

Methods

We attempted to resolve these issues by decreasing the scheme size to reduce the occurrence of error and missing alleles, alleviate the storage burden, and improve the query search time. The challenge in this approach is maintaining the typing resolution when using fewer loci. We achieved this by using a popular artificial intelligence technique, XGBoost, coupled with Shapley additive explanations for feature selection. Finally, 370 loci from the original 1701 cgMLST loci of Listeria monocytogenes were selected.

Results

Although the size of the final scheme (LmScheme_370) was approximately 80% lower than that of the original cgMLST scheme, its discriminatory power, tested for 35 outbreaks, was concordant with that of the original cgMLST scheme. Although we used L. monocytogenes as a demonstration in this study, the approach can be applied to other schemes and pathogens. Our findings might help elucidate gene-by-gene–based epidemiology.

Pan-Genome-Wide Analysis of Pantoea ananatis Identified Genes Linked to Pathogenicity in Onion

Pantoea ananatis, a gram negative and facultative anaerobic bacterium is a member of a Pantoea spp. complex that causes center rot of onion, which significantly affects onion yield and quality. This pathogen does not have typical virulence factors like type II or type III secretion systems but appears to require a biosynthetic gene-cluster, HiVir/PASVIL (located chromosomally comprised of 14 genes), for a phosphonate secondary metabolite, and the 'alt' gene cluster (located in plasmid and comprised of 11 genes) that aids in bacterial colonization in onion bulbs by imparting tolerance to thiosulfinates. We conducted a deep pan-genome-wide association study (pan-GWAS) to predict additional genes associated with pathogenicity in P. ananatis using a panel of diverse strains (n = 81). We utilized a red-onion scale necrosis assay as an indicator of pathogenicity. Based on this assay, we differentiated pathogenic (n = 51)- vs. non-pathogenic (n = 30)-strains phenotypically. Pan-genome analysis revealed a large core genome of 3,153 genes and a flexible accessory genome. Pan-GWAS using the presence and absence variants (PAVs) predicted 42 genes, including 14 from the previously identified HiVir/PASVIL cluster associated with pathogenicity, and 28 novel genes that were not previously associated with pathogenicity in onion. Of the 28 novel genes identified, eight have annotated functions of site-specific tyrosine kinase, N-acetylmuramoyl-L-alanine amidase, conjugal transfer, and HTH-type transcriptional regulator. The remaining 20 genes are currently hypothetical. Further, a core-genome SNPs-based phylogeny and horizontal gene transfer (HGT) studies were also conducted to assess the extent of lateral gene transfer among diverse P. ananatis strains. Phylogenetic analysis based on PAVs and whole genome multi locus sequence typing (wgMLST) rather than core-genome SNPs distinguished red-scale necrosis inducing (pathogenic) strains from non-scale necrosis inducing (non-pathogenic) strains of P. ananatis. A total of 1182 HGT events including the HiVir/PASVIL and alt cluster genes were identified. These events could be regarded as a major contributing factor to the diversification, niche-adaptation and potential acquisition of pathogenicity/virulence genes in P. ananatis.

Pan-Genome of Novel Pantoea stewartii subsp. indologenes Reveals Genes Involved in Onion Pathogenicity and Evidence of Lateral Gene Transfer

Pantoea stewartii subsp. indologenes (Psi) is a causative agent of leafspot on foxtail millet and pearl millet; however, novel strains were recently identified that are pathogenic on onions. Our recent host range evaluation study identified two pathovars; P. stewartii subsp. indologenes pv. cepacicola pv. nov. and P. stewartii subsp. indologenes pv. setariae pv. nov. that are pathogenic on onions and millets or on millets only, respectively. In the current study, we developed a pan-genome using the whole genome sequencing of newly identified/classified Psi strains from both pathovars [pv. cepacicola (n = 4) and pv. setariae (n = 13)]. The full spectrum of the pan-genome contained 7030 genes. Among these, 3546 (present in genomes of all 17 strains) were the core genes that were a subset of 3682 soft-core genes (present in ≥16 strains). The accessory genome included 1308 shell genes and 2040 cloud genes (present in ≤2 strains). The pan-genome showed a clear linear progression with >6000 genes, suggesting that the pan-genome of Psi is open. Comparative phylogenetic analysis showed differences in phylogenetic clustering of Pantoea spp. using PAVs/wgMLST approach in comparison with core genome SNPs-based phylogeny. Further, we conducted a horizontal gene transfer (HGT) study using Psi strains from both pathovars along with strains from other Pantoea species, namely, P. stewartii subsp. stewartii LMG 2715T, P. ananatis LMG 2665T, P. agglomerans LMG L15, and P. allii LMG 24248T. A total of 317 HGT events among four Pantoea species were identified with most gene transfer events occurring between Psi pv. cepacicola and Psi pv. setariae. Pan-GWAS analysis predicted a total of 154 genes, including seven gene-clusters, which were associated with the pathogenicity phenotype (necrosis on seedling) on onions. One of the gene-clusters contained 11 genes with known functions and was found to be chromosomally located.

Phylogenetic relationship and genomic characterization of Salmonella Typhimurium strains isolated from swine in Brazil

Salmonella Typhimurium has been transmitted between humans and animals. Although, Brazil has been one of the largest pork meat exporters worldwide, there are few studies that characterized epidemiologically S. Typhimurium strains from swine. The aims of this work were to study the phylogenetic relationship of S. Typhimurium genomes isolated from swine in Brazil among themselves and with other genomes isolated from several sources and countries using wgMLST and cgMLST and to perform the search of Salmonella pathogenicity islands (SPIs). In addition, for S. Typhimurium strains from swine to compare the virulence and antimicrobial resistance genes by VFDB and ResFinder, genetic content by BLAST Atlas and orthologous proteins clusters by OrthoVenn. The constructed phylogenetic trees by wgMLST and cgMLST grouped the majority (92.3% and 80.7%, respectively) of the strains isolated from swine in Brazil into the same group. All the isolates contained important SPIs (SPI-1, SPI-2, SPI-3, SPI-5 and SPI-9). A total of 100 and 31 virulence and resistance genes were detected in the S. Typhimurium strains isolated from swine, respectively. The BLAST Atlas and orthologous proteins analysis found regions of phages and differences in metabolic, regulatory and cellular processes among S. Typhimurium LT2 and S. Typhimurium isolates from swine. In conclusion, molecular typing based in the wgMLST and cgMLST suggested that the S. Typhimurium isolates from swine studied were genetically related. The pathogenic potential of the strains studied was corroborated by the presence of important SPIs and virulence genes. The high number of antimicrobial resistance genes detected is worrying and reinforced their potential risk in swine in Brazil. The comparison by BLAST Atlas suggested differences in mobile genetic elements among S. Typhimurium LT2 and S. Typhimurium isolates from swine in Brazil. The orthologous proteins analysis revealed unique genes related to important cellular processes in the strains from swine.

Pan-genomic analyses of 47 complete genomes of the Rickettsia genus and prediction of new vaccine targets and virulence factors of the species

The genus Rickettsia belongs to the Proteobacteria phylum and these bacteria infect animals and humans causing a range of diseases worldwide. The genus is divided into 4 groups and despite the public health threat and the knowledge accumulated so far, the mandatory intracellular bacteria behaviour and limitation for in vitro culture makes it difficult to create new vaccines and drug targets to these bacteria. In an attempt to overcome these limitations, pan-genomic approaches has used 47 genomes of the genus Rickettsia, in order to describe species similarities and genomics islands. Moreover, we conducted reverse vaccinology and docking analysis aiming the identification of proteins that have great potential to become vaccine and drug targets. We found out that the bacteria of the four Rickettsia groups have a high similarity with each other, with about 90 to 100% of identity. A pathogenicity island and a resistance island were predicted. In addition, 8 proteins were also predicted as strong candidates for vaccine and 9 as candidates for drug targets. The prediction of the proteins leads us to believe in a possibility of prospecting potential drugs or creating a polyvalent vaccine, which could reach most strains of this large group of bacteria.Communicated by Ramaswamy H. Sarma.

Serratia marcescens outbreak in a neonatal intensive care unit and the potential of whole-genome sequencing

BACKGROUND

Serratia marcescens is notorious for its increasing antimicrobial resistance and potential to cause outbreaks in neonatal intensive care units (NICUs). A promising tool in outbreak investigations is whole-genome sequencing (WGS).

OBJECTIVES

To describe a S. marcescens outbreak (2018-2019) in an NICU and discuss which infection control measures contributed to containment, addressing the potential of WGS.

METHODS

S. marcescens isolates from patients and the environment isolated during the 2018-2019 NICU outbreak were analysed. In comparison, isolates from previous presumed NICU outbreaks and adult blood cultures were included. WGS and whole-genome multi-locus sequence typing analysis were performed.

RESULTS

Sixty-three S. marcescens isolates were analysed. The 2018-2019 outbreak was divided into three clusters, including four environmental strains (drains, N=3; baby scale, N=1). The strains differed significantly from those of an NICU outbreak in 2014 and adult blood cultures. Besides standard infection control measures, the siphons were replaced and weekly decontamination was performed with acetic acid 10%. Seven acquired-resistance genes and 29 virulence-associated genes were detected.

CONCLUSIONS

It was assumed that both neonates and drains were reservoirs of S. marcescens cross-contamination via the hands of healthcare workers and parents. Initially, standard measures, including hand hygiene, were reinforced. However, definitive containment was achieved only after replacement of the siphons and weekly decontamination with acetic acid. WGS enables faster recognition of an outbreak with accurate mapping of the spread, facilitating the implementation of infection control measures. WGS also provides interesting information about the spread of antibiotic resistance and virulence genes.

cgMLST@Taiwan: A web service platform for Vibrio cholerae cgMLST profiling and global strain tracking

BACKGROUND

Cholera, a rapidly dehydrating diarrheal disease caused by toxigenic Vibrio cholerae, is a leading cause of morbidity and mortality in some regions of the world. Core genome multilocus sequence typing (cgMLST) is a promising approach in generating genetic fingerprints from whole-genome sequencing (WGS) data for strain comparison among laboratories.

METHODS

We constructed a V. cholerae core gene allele database using an in-house developed computational pipeline, a database with cgMLST profiles converted from genomic sequences from the National Center for Biotechnology Information, and built a REST-based web accessible via the Internet.

RESULTS

We built a web service platform-cgMLST@Taiwan and installed a V. cholerae allele database, a cgMLST profile database, and computational tools for generating V. cholerae cgMLST profiles (based on 3,017 core genes), performing rapid global strain tracking, and clustering analysis of cgMLST profiles. This web-based platform provides services to researchers, public health microbiologists, and physicians who use WGS data for the investigation of cholera outbreaks and tracking of V. cholerae strain transmission across countries and geographic regions. The cgMLST@Taiwan is accessible at http://rdvd.cdc.gov.tw/cgMLST.

Genomic Characterization of Multidrug-Resistant Escherichia coli BH100 Sub-strains

The rapid emergence of multidrug-resistant (MDR) bacteria is a global health problem. Mobile genetic elements like conjugative plasmids, transposons, and integrons are the major players in spreading resistance genes in uropathogenic Escherichia coli (UPEC) pathotype. The E. coli BH100 strain was isolated from the urinary tract of a Brazilian woman in 1974. This strain presents two plasmids carrying MDR cassettes, pBH100, and pAp, with conjugative and mobilization properties, respectively. However, its transposable elements have not been characterized. In this study, we attempted to unravel the factors involved in the mobilization of virulence and drug-resistance genes by assessing genomic rearrangements in four BH100 sub-strains (BH100 MG2014, BH100 MG2017, BH100L MG2017, and BH100N MG2017). Therefore, the complete genomes of the BH100 sub-strains were achieved through Next Generation Sequencing and submitted to comparative genomic analyses. Our data shows recombination events between the two plasmids in the sub-strain BH100 MG2017 and between pBH100 and the chromosome in BH100L MG2017. In both cases, IS3 and IS21 elements were detected upstream of Tn21 family transposons associated with MDR genes at the recombined region. These results integrated with Genomic island analysis suggest pBH100 might be involved in the spreading of drug resistance through the formation of resistance islands. Regarding pathogenicity, our results reveal that BH100 strain is closely related to UPEC strains and contains many IS3 and IS21-transposase-enriched genomic islands associated with virulence. This study concludes that those IS elements are vital for the evolution and adaptation of BH100 strain.

Comparative Genomic Analysis Provides Insights into the Phylogeny, Resistome, Virulome, and Host Adaptation in the Genus Ewingella

Ewingella americana is a cosmopolitan bacterial pathogen that has been isolated from many hosts. Here, we sequenced a high-quality genome of E. americana B6-1 isolated from Flammulina filiformis, an important cultivated mushroom, performed a comparative genomic analysis with four other E. americana strains from various origins, and tested the susceptibility of B6-1 to antibiotics. The genome size, predicted genes, and GC (guanine-cytosine) content of B6-1 was 4.67 Mb, 4301, and 53.80%, respectively. The origin of the strains did not significantly affect the phylogeny, but mobile genetic elements shaped the evolution of the genus Ewingella. The strains encoded a set of common genes for type secretion, virulence effectors, CAZymes, and toxins required for pathogenicity in all hosts. They also had antibiotic resistance, pigments to suppress or evade host defense responses, as well as genes for adaptation to different environmental conditions, including temperature, oxidation, and nutrients. These findings provide a better understanding of the virulence, antibiotic resistance, and host adaptation strategies of Ewingella, and they also contribute to the development of effective control strategies.

Activity of Imipenem-Relebactam and Meropenem-Vaborbactam against Carbapenem-Resistant, SME-Producing Serratia marcescens

The Serratia marcescens enzyme (SME) is a chromosomally encoded carbapenemase with no known optimal treatment. Various β-lactam/β-lactamase inhibitors and comparators were evaluated against 8 SME producers via broth microdilution. Four isolates were subsequently tested via time-kill analyses. All isolates were resistant to imipenem, imipenem-relebactam, and meropenem but susceptible to ceftazidime, ceftazidime-avibactam, and meropenem-vaborbactam. Ceftazidime, imipenem-relebactam, and meropenem-vaborbactam were bactericidal against 3, 0, and 4 isolates, respectively. Meropenem-vaborbactam may be a potential option for severe SME-producing infections.

Phylogenetic analysis revealed that Salmonella Typhimurium ST313 isolated from humans and food in Brazil presented a high genomic similarity

Salmonella Typhimurium sequence type 313 (S. Typhimurium ST313) has caused invasive disease mainly in sub-Saharan Africa. In Brazil, ST313 strains have been recently described, and there is a lack of studies that assessed by whole genome sequencing (WGS)-the relationship of these strains. The aims of this work were to study the phylogenetic relationship of 70 S. Typhimurium genomes comparing strains of ST313 (n = 9) isolated from humans and food in Brazil among themselves, with other STs isolated in this country (n = 31) and in other parts of the globe (n = 30) by 16S rRNA sequences, the Gegenees software, whole genome multilocus sequence typing (wgMLST), and average nucleotide identity (ANI) for the genomes of ST313. Additionally, pangenome analysis was performed to verify the heterogeneity of these genomes. The phylogenetic analyses showed that the ST313 genomes were very similar among themselves. However, the ST313 genomes were usually clustered more distantly to other STs of strains isolated in Brazil and in other parts of the world. By pangenome calculation, the core genome was 2,880 CDSs and 4,171 CDSs singletons for all the 70 S. Typhimurium genomes studied. Considering the 10 ST313 genomes analyzed the core genome was 4,112 CDSs and 76 CDSs singletons. In conclusion, the ST313 genomes from Brazil showed a high similarity among them which information might eventually help in the development of vaccines and antibiotics. The pangenome analysis showed that the S. Typhimurium genomes studied presented an open pangenome, but specifically tending to become close for the ST313 strains.

Whole-Genome Sequencing of Corynebacterium diphtheriae Isolates Recovered from an Inner-City Population Demonstrates the Predominance of a Single Molecular Strain

In some parts of the world, Corynebacterium diphtheriae has reemerged as a pathogen, especially as a cause of infections among impoverished and marginalized populations. We performed whole-genome sequencing (WGS) on all cutaneous C. diphtheriae isolates (n = 56) from Vancouver's inner-city population over a 3-year time period (2015 to 2018). All isolates with complete genome assembly were toxin negative, contained a common set of 22 virulence factors, and shared a highly conserved accessory genome. One of our isolates harbored a novel plasmid conferring macrolide and lincosamide resistance. Fifty-two out of 56 isolates were multilocus sequence type 76, and single nucleotide variants (SNV) and core-genome multilocus sequence typing (cgMLST) analysis demonstrated tight clustering of our isolates relative to all publicly available C. diphtheriae genomes. All sequence type 76 (ST76) study isolates were within a median of 22 SNVs and 13 cgMLST alleles of each other, while NCBI genomes were within a median of 17,436 SNVs and 1,552 cgMLST alleles of each other (both P < 2.2 × 10^-16). A single strain of C. diphtheriae appears to be causing cutaneous infections in the low-income population of Vancouver. Further research is needed to elucidate transmission networks in our study population and standardize C. diphtheriae epidemiological typing when whole genomes are sequenced.

Pangenome of Serratia marcescens strains from nosocomial and environmental origins reveals different populations and the links between them

Serratia marcescens is a Gram-negative bacterial species that can be found in a wide range of environments like soil, water and plant surfaces, while it is also known as an opportunistic human pathogen in hospitals and as a plant growth promoting bacteria (PGPR) in crops. We have used a pangenome-based approach, based on publicly available genomes, to apply whole genome multilocus sequence type schemes to assess whether there is an association between source and genotype, aiming at differentiating between isolates from nosocomial sources and the environment, and between strains reported as PGPR from other environmental strains. Most genomes from a nosocomial setting and environmental origin could be assigned to the proposed nosocomial or environmental MLSTs, which is indicative of an association between source and genotype. The fact that a few genomes from a nosocomial source showed an environmental MLST suggests that a minority of nosocomial strains have recently derived from the environment. PGPR strains were assigned to different environmental types and clades but only one clade comprised strains accumulating a low number of known virulence and antibiotic resistance determinants and was exclusively from environmental sources. This clade is envisaged as a group of promissory MLSTs for selecting prospective PGPR strains.

Persistent Infection and Long-Term Carriage of Typhoidal and Nontyphoidal Salmonellae

The ability of pathogenic bacteria to affect higher organisms and cause disease is one of the most dramatic properties of microorganisms. Some pathogens can establish transient colonization only, but others are capable of infecting their host for many years or even for a lifetime. Long-term infection is called persistence, and this phenotype is fundamental for the biology of important human pathogens, including Helicobacter pylori, Mycobacterium tuberculosis, and Salmonella enterica Both typhoidal and nontyphoidal serovars of the species Salmonella enterica can cause persistent infection in humans; however, as these two Salmonella groups cause clinically distinct diseases, the characteristics of their persistent infections in humans differ significantly. Here, following a general summary of Salmonella pathogenicity, host specificity, epidemiology, and laboratory diagnosis, I review the current knowledge about Salmonella persistence and discuss the relevant epidemiology of persistence (including carrier rate, duration of shedding, and host and pathogen risk factors), the host response to Salmonella persistence, Salmonella genes involved in this lifestyle, as well as genetic and phenotypic changes acquired during prolonged infection within the host. Additionally, I highlight differences between the persistence of typhoidal and nontyphoidal Salmonella strains in humans and summarize the current gaps and limitations in our understanding, diagnosis, and curing of persistent Salmonella infections.

Pantoea ananatis Genetic Diversity Analysis Reveals Limited Genomic Diversity as Well as Accessory Genes Correlated with Onion Pathogenicity

Pantoea ananatis is a member of the family Enterobacteriaceae and an enigmatic plant pathogen with a broad host range. Although P. ananatis strains can be aggressive on onion causing foliar necrosis and onion center rot, previous genomic analysis has shown that P. ananatis lacks the primary virulence secretion systems associated with other plant pathogens. We assessed a collection of fifty P. ananatis strains collected from Georgia over three decades to determine genetic factors that correlated with onion pathogenic potential. Previous genetic analysis studies have compared strains isolated from different hosts with varying diseases potential and isolation sources. Strains varied greatly in their pathogenic potential and aggressiveness on different cultivated Allium species like onion, leek, shallot, and chive. Using multi-locus sequence analysis (MLSA) and repetitive extragenic palindrome repeat (rep)-PCR techniques, we did not observe any correlation between onion pathogenic potential and genetic diversity among strains. Whole genome sequencing and pan-genomic analysis of a sub-set of 10 strains aided in the identification of a novel series of genetic regions, likely plasmid borne, and correlating with onion pathogenicity observed on single contigs of the genetic assemblies. We named these loci Onion Virulence Regions (OVR) A-D. The OVR loci contain genes involved in redox regulation as well as pectate lyase and rhamnogalacturonase genes. Previous studies have not identified distinct genetic loci or plasmids correlating with onion foliar pathogenicity or pathogenicity on a single host pathosystem. The lack of focus on a single host system for this phytopathgenic disease necessitates the pan-genomic analysis performed in this study.

Identification of a novel botulinum neurotoxin gene cluster in Enterococcus

The deadly neurotoxins of Clostridium botulinum (BoNTs) comprise eight serotypes (A-G; X). The neurotoxin gene cluster encoding BoNT and its accessory proteins includes an operon containing an ntnh gene upstream of the boNT gene. Another operon contains either ha (haemagglutinin) or orfX genes (of unknown function). Here we describe a novel boNT gene cluster from Enterococcus sp. 3G1_DIV0629, with a typical ntnh gene and an uncommon orfX arrangement. The neurotoxin (designated putative eBoNT/J) contains a metallopeptidase zinc-binding site, a translocation domain and a target cell attachment domain. Structural properties of the latter suggest a novel targeting mechanism with consequent implications for application by the pharmaceutical industry. This is the first complete boNT gene cluster identified in a non-clostridial genome.

New Insights into the Diversity of the Genus Faecalibacterium

Faecalibacterium prausnitzii is a commensal bacterium, ubiquitous in the gastrointestinal tracts of animals and humans. This species is a functionally important member of the microbiota and studies suggest it has an impact on the physiology and health of the host. F. prausnitzii is the only identified species in the genus Faecalibacterium, but a recent study clustered strains of this species in two different phylogroups. Here, we propose the existence of distinct species in this genus through the use of comparative genomics. Briefly, we performed analyses of 16S rRNA gene phylogeny, phylogenomics, whole genome Multi-Locus Sequence Typing (wgMLST), Average Nucleotide Identity (ANI), gene synteny, and pangenome to better elucidate the phylogenetic relationships among strains of Faecalibacterium. For this, we used 12 newly sequenced, assembled, and curated genomes of F. prausnitzii, which were isolated from feces of healthy volunteers from France and Australia, and combined these with published data from 5 strains downloaded from public databases. The phylogenetic analysis of the 16S rRNA sequences, together with the wgMLST profiles and a phylogenomic tree based on comparisons of genome similarity, all supported the clustering of Faecalibacterium strains in different genospecies. Additionally, the global analysis of gene synteny among all strains showed a highly fragmented profile, whereas the intra-cluster analyses revealed larger and more conserved collinear blocks. Finally, ANI analysis substantiated the presence of three distinct clusters—A, B, and C—composed of five, four, and four strains, respectively. The pangenome analysis of each cluster corroborated the classification of these clusters into three distinct species, each containing less variability than that found within the global pangenome of all strains. Here, we propose that comparison of pangenome subsets and their associated α values may be used as an alternative approach, together with ANI, in the in silico classification of new species. Altogether, our results provide evidence not only for the reconsideration of the phylogenetic and genomic relatedness among strains currently assigned to F. prausnitzii, but also the need for lineage (strain-based) differentiation of this taxon to better define how specific members might be associated with positive or negative host interactions.

PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Although PacBio third-generation sequencers have improved the read lengths of genome sequencing which facilitates the assembly of complete genomes, no study has reported success in using PacBio data alone to completely sequence a two-chromosome bacterial genome from a single library in a single run. Previous studies using earlier versions of sequencing chemistries have at most been able to finish bacterial genomes containing only one chromosome with de novo assembly. In this study, we compared the robustness of PacBio RS II, using one SMRT cell and the latest P6-C4 chemistry, with Illumina HiSeq 1500 in sequencing the genome of Burkholderia pseudomallei, a bacterium which contains two large circular chromosomes, very high G+C content of 68–69%, highly repetitive regions and substantial genomic diversity, and represents one of the largest and most complex bacterial genomes sequenced, using a reference genome generated by hybrid assembly using PacBio and Illumina datasets with subsequent manual validation. Results showed that PacBio data with de novo assembly, but not Illumina, was able to completely sequence the B. pseudomallei genome without any gaps or mis-assemblies. The two large contigs of the PacBio assembly aligned unambiguously to the reference genome, sharing >99.9% nucleotide identities. Conversely, Illumina data assembled using three different assemblers resulted in fragmented assemblies (201–366 contigs), sharing only 92.2–100% and 92.0–100% nucleotide identities to chromosomes I and II reference sequences, respectively, with no indication that the B. pseudomallei genome consisted of two chromosomes with four copies of ribosomal operons. Among all assemblies, the PacBio assembly recovered the highest number of core and virulence proteins, and housekeeping genes based on whole-genome multilocus sequence typing (wgMLST). Most notably, assembly solely based on PacBio outperformed even hybrid assembly using both PacBio and Illumina datasets. Hybrid approach generated only 74 contigs, while the PacBio data alone with de novo assembly achieved complete closure of the two-chromosome B. pseudomallei genome without additional costly bench work and further sequencing. PacBio RS II using P6-C4 chemistry is highly robust and cost-effective and should be the platform of choice in sequencing bacterial genomes, particularly for those that are well-known to be difficult-to-sequence.

In-silico Taxonomic Classification of 373 Genomes Reveals Species Misidentification and New Genospecies within the Genus Pseudomonas

The genus Pseudomonas has one of the largest diversity of species within the Bacteria kingdom. To date, its taxonomy is still being revised and updated. Due to the non-standardized procedure and ambiguous thresholds at species level, largely based on 16S rRNA gene or conventional biochemical assay, species identification of publicly available Pseudomonas genomes remains questionable. In this study, we performed a large-scale analysis of all Pseudomonas genomes with species designation (excluding the well-defined P. aeruginosa) and re-evaluated their taxonomic assignment via in silico genome-genome hybridization and/or genetic comparison with valid type species. Three-hundred and seventy-three pseudomonad genomes were analyzed and subsequently clustered into 145 distinct genospecies. We detected 207 erroneous labels and corrected 43 to the proper species based on Average Nucleotide Identity Multilocus Sequence Typing (MLST) sequence similarity to the type strain. Surprisingly, more than half of the genomes initially designated as Pseudomonas syringae and Pseudomonas fluorescens should be classified either to a previously described species or to a new genospecies. Notably, high pairwise average nucleotide identity (>95%) indicating species-level similarity was observed between P. synxantha-P. libanensis, P. psychrotolerans–P. oryzihabitans, and P. kilonensis- P. brassicacearum, that were previously differentiated based on conventional biochemical tests and/or genome-genome hybridization techniques.

Construction of a Pan-Genome Allele Database of Salmonella enterica Serovar Enteritidis for Molecular Subtyping and Disease Cluster Identification

We built a pan-genome allele database with 395 genomes of Salmonella enterica serovar Enteritidis and developed computer tools for analysis of whole genome sequencing (WGS) data of bacterial isolates for disease cluster identification. A web server (http://wgmlst.imst.nsysu.edu.tw) was set up with the database and the tools, allowing users to upload WGS data to generate whole genome multilocus sequence typing (wgMLST) profiles and to perform cluster analysis of wgMLST profiles. The usefulness of the database in disease cluster identification was demonstrated by analyzing a panel of genomes from 55 epidemiologically well-defined S. Enteritidis isolates provided by the Minnesota Department of Health. The wgMLST-based cluster analysis revealed distinct clades that were concordant with the epidemiologically defined outbreaks. Thus, using a common pan-genome allele database, wgMLST can be a promising WGS-based subtyping approach for disease surveillance and outbreak investigation across laboratories.