Gegenees: fragmented alignment of multiple genomes for determining phylogenomic distances and genetic signatures unique for specified target groups

Characterization of the diversity of type IV secretion system-encoding plasmids in Acinetobacter

The multi-drug resistant pathogen Acinetobacter baumannii has gained global attention as an important clinical challenge. Owing to its ability to survive on surfaces, its capacity for horizontal gene transfer, and its resistance to front-line antibiotics, A. baumannii has established itself as a successful pathogen. Bacterial conjugation is a central mechanism for pathogen evolution. The epidemic multidrug-resistant A. baumannii ACICU harbours a plasmid encoding a Type IV Secretion System (T4SS) with homology to the E. coli F-plasmid, and plasmids with homologous gene clusters have been identified in several A. baumannii sequence types. However the genetic and host strain diversity, global distribution, and functional ability of this group of plasmids is not fully understood. Using systematic analysis, we show that pACICU2 belongs to a group of almost 120 T4SS-encoding plasmids within four different species of Acinetobacter and one strain of Klebsiella pneumoniae from human and environmental origin, and globally distributed across 20 countries spanning 4 continents. Genetic diversity was observed both outside and within the T4SS-encoding cluster, and 47% of plasmids harboured resistance determinants, with two plasmids harbouring eleven. Conjugation studies with an extensively drug-resistant (XDR) strain showed that the XDR plasmid could be successfully transferred to a more divergent A. baumanii, and transconjugants exhibited the resistance phenotype of the plasmid. Collectively, this demonstrates that these T4SS-encoding plasmids are globally distributed and more widespread among Acinetobacter than previously thought, and that they represent an important potential reservoir for future clinical concern.

Insights into the genomic traits of Yersinia frederiksenii, Yersinia intermedia and Yersinia kristensenii isolated from diverse sources in Brazil

Yersinia is an important genus comprising foodborne, zoonotic and pathogenic bacteria. On the other hand, species of the so-called group Yersinia enterocolitica-like are understudied and mostly characterized as non-pathogenic, despite of some reports of human infections. The present study aimed to provide genomic insights of Yersinia frederiksenii (YF), Yersinia intermedia (YI) and Yersinia kristensenii (YK) isolated worldwide. A total of 22 YF, 20 YI and 14 YK genomes were searched for antimicrobial resistance genes, plasmids, prophages, and virulence factors. Their phylogenomic relatedness was analyzed by Gegenees and core-genome multi-locus sequence typing. Beta-lactam resistance gene blaTEM-116 and five plasmids replicons (pYE854, ColRNAI, ColE10, Col(pHAD28) and IncN3) were detected in less than five genomes. A total of 59 prophages, 106 virulence markers of the Yersinia genus, associated to adherence, antiphagocytosis, exoenzymes, invasion, iron uptake, proteases, secretion systems and the O-antigen, and virulence factors associated to other 20 bacterial genera were detected. Phylogenomic analysis revealed high inter-species distinction and four highly diverse YF clusters. In conclusion, the results obtained through the analyses of YF, YI and YK genomes suggest the virulence potential of these strains due to the broad diversity and high frequency of prophages and virulence factors found. Phylogenetic analyses were able to correctly distinguish these closely related species and show the presence of different genetic subgroups. These data contributed for a better understanding of YF, YI and YK virulence-associated features and global genetic diversity, and reinforced the need for better characterization of these Y. enterocolitica-like species considered non-pathogenic.

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.

Molecular characterization and environmental impact of newly isolated lytic phage SLAM_phiST1N3 in the Cornellvirus genus for biocontrol of a multidrug-resistant Salmonella Typhimurium in the swine industry chain

Salmonella Typhimurium is a highly lethal pathogenic bacterium in weaned piglets, causing significant treatment costs and economic losses in the swine industry. Additionally, due to its ability to induce zoonotic diseases, resulting in harm to humans through the transmission of the pathogen from pork, it presents a serious public health issue. Bacteriophages (phages), viruses that infect specific bacterial strains, have been proposed as an alternative to antibiotics for controlling pathogenic bacteria. In this study, we isolated SLAM_phiST1N3, a phage infecting a multidrug-resistant (MDR) S. Typhimurium wild-type strain isolated from diseased pigs. First, comparative genomics and phylogenetic analysis revealed that SLAM_phiST1N3 belongs to the Cornellvirus genus. Moreover, utilizing a novel classification approach introduced in this study, SLAM_phiST1N3 was classified at the species level. Host range experiments demonstrated that SLAM_phiST1N3 did not infect other pathogenic bacteria or probiotics derived from pigs or other livestock. While complete eradication of Salmonella was not achievable in the liquid inhibition assay, surprisingly, we succeeded in largely eliminating Salmonella in the FIMM analysis, a gut simulation system using weaned piglet feces. Furthermore, using the C. elegans model, we showcased the potential of SLAM_phiST1N3 to prevent S. Typhimurium infection in living organisms. In addition, it was confirmed that bacterial control could be achieved when phage was applied to Salmonella-contaminated pork. pH and temperature stability experiments demonstrated that SLAM_phiST1N3 can endure swine industry processes and digestive conditions. In conclusion, SLAM_phiST1N3 demonstrates potential environmental impact as a substance for Salmonella prevention across various aspects of the swine industry chain.

Genomic analysis of seven mycobacteriophages identifies three novel species with differing phenotypic stabilities

Recently, case studies have been published regarding the application of mycobacteriophage (MP) therapy (MPT) in patients with multi-antibiotic-resistant infections. A major limitation in the development of MPT is the paucity of therapeutically useful MP. As there are approximately 10,000 MP that have yet to be sequenced, it is possible that characterization of this cohort would increase the repertoire of useful MP. This study aims to contribute to such a strategy, by characterizing a cohort of 7 mycobacteriophages. Sequencing analyses revealed that the MP have unique sequences, and subsequent gene annotation revealed differences in gene organization. Notably, MP LOCARD has the largest genome and operons encoding for glycosyltransferases. Taxonomic analysis executed with VIRIDIC, Gegenees and VICTOR revealed that LOCARD belongs to a different genus than the other phages and is the foundational member of one of three novel species identified in this study. LOCARD, LOCV2, and LOCV5 were selected as representative members of their species and subjected to phenotypic analyses to compare their stability under biologically and industrially relevant conditions. Again LOCARD stood out, as it was unaffected by the typical temperatures (37 °C) and salinity (0.9%) experienced in mammals, while the viability of LOCV2 and LOCV5 was significantly reduced. LOCARD was also tolerant to pH 10, low levels of antiviral detergent and was the least impacted by a single freeze-thaw cycle. When all these results are considered, it indicates that LOCARD in particular, has potential therapeutic and/or diagnostics applications, given its resilience towards physiological and storage conditions.

Microbial Diversity in Children with Gastroenteritis in the Amazon Region of Brazil: Development and Validation of a Molecular Method for Complete Sequencing of Viral Genomes

INTRODUCTION

Metagenomic sequencing is a powerful tool that is widely used in laboratories worldwide for taxonomic characterization of microorganisms in clinical and environmental samples. In this study, we utilized metagenomics to investigate comprehensively the microbial diversity in fecal samples of children over a four-year period. Our methods were carefully designed to ensure accurate and reliable results.

MATERIAL AND METHODS

Validated and analyzed were metagenomic data obtained from sequencing 27 fecal samples from children under 10 years old with gastroenteritis over a four-year period (2012-2016). The fecal specimens were collected from patients who received care at public health facilities in the northern region of Brazil. Sequencing libraries were prepared from cDNA and sequenced on the Illumina HiSeq. Kraken-2 was utilized to classify bacterial taxonomy based on the 16S rRNA gene, using the Silva rRNA database. Additionally, the Diamond program was used for mapping to the non-redundant protein database (NR database). Phylogenomic analyses were conducted using Geneious R10 and MEGA X software, and Bayesian estimation of phylogeny was performed using the MrBayes program. The results indicate significant heterogeneity among norovirus strains, with evidence of recombination and point mutations. This study presents the first complete genome of parechovirus 8 in the region. Additionally, it describes the bacterial populations and bacteriophages present in feces, with a high abundance of Firmicutes and Proteobacteria, including an increased proportion of the Enterobacteriaceae family. The presented data demonstrate the genetic diversity of microbial populations and provide a comprehensive report on viral molecular characterization. These findings are relevant for genomic studies in gastrointestinal infections. The metagenomic approach is a powerful tool for investigating microbial diversity in children with gastroenteritis. However, further studies are imperative to conduct genomic analysis of identified bacterial strains and thoroughly analyze antimicrobial resistance genes.

Vancomycin-resistant Enterococcus faecium in Japan, 2007-2015: a molecular epidemiology analysis focused on examining strain characteristics over time

IMPORTANCE

Our study emphasizes the efficacy of whole-genome sequencing (WGS) in addressing outbreaks of vancomycin-resistant enterococci. WGS enables the identification and tracking of resistant bacterial strains, early detection and management of novel infectious disease outbreaks, and the appropriate selection and use of antibiotics. Furthermore, our approach deepens our understanding of how resistant bacteria transfer genes and adapt to their environments or hosts. For modern medicine, these insights have significant implications for controlling infections and effectively managing antibiotic use in the current era, where antibiotic resistance is progressing.

Bioinformatic Analysis of Staphylococcus Phages: A Key Step for Safe Cocktail Development

An in-depth analysis of phage genomic sequences is essential for the proposal of a cocktail for therapeutic uses. With the burst of publications on phage isolation and genetic studies during the last decade, several different bioinformatics programs have been used. Here we describe our studies on the genetic organization of phages infecting Staphylococcus aureus, a pathogen of human importance, by using an assembly of tools for gene annotation, identification of expression components, and phylogeny analysis.

Identification of novel genera and subcluster classifications for mycobacteriophages

Aim: To identify novel genera amongst mycobacteriophages (MP) and verify a hypothesised correlation between the taxonomy set by the International Committee on Taxonomy of Viruses (ICTV) and the National Centre for Biotechnology Information (NCBI) with that of the Actinobacteriophage Database, which may help formalise subcluster assignment. Methods: A dataset of 721 MP genomes was analysed using VIRIDIC, a nucleotide alignment-based software that predicts genus assignments. Potentially novel genera were analysed using Gegenees and VICTOR, respectively. These genera were then compared to the subclusters assigned by the Actinobacteriophage Database to verify a hypothesis that one genus can be assigned to one subcluster (i.e., the genus-subcluster hypothesis). Results: Initially, when comparing the current genus classifications of the 721 MP dataset to the Actinobacteriophage database subcluster assignments, 83.3% of subclusters supported the genus-subcluster hypothesis. Following the sequential VIRIDIC, Gegenees and VICTOR analyses, a total of 20 novel genera were identified based on a ≥ 70% and ~ 50% similarity threshold for VIRIDIC and Gegenees, respectively, and a monophyletic nature in the VICTOR output. Interestingly, these criteria also appear to support the creation of 13 novel subclusters, which would increase the support for the genus-subcluster hypothesis to 97.6%. Conclusion: The link between genus and subcluster classifications appears robust, as most subclusters can be assigned a single genus and vice versa. By relating the taxonomic and clustering classification systems, they can be easily kept up to date to best reflect MP diversity, which could aid the rapid selection of related (or diverse) phages for research, therapeutic and diagnostic purposes.

Enrichment Culture but Not Metagenomic Sequencing Identified a Highly Prevalent Phage Infecting Lactiplantibacillus plantarum in Human Feces

Lactiplantibacillus plantarum (previously known as Lactobacillus plantarum) is increasingly used as a probiotic to treat human diseases, but its phages in the human gut remain unexplored. Here, we report its first gut phage, Gut-P1, which we systematically screened using metagenomic sequencing, virus-like particle (VLP) sequencing, and enrichment culture from 35 fecal samples. Gut-P1 is virulent, belongs to the Douglaswolinvirus genus, and is highly prevalent in the gut (~11% prevalence); it has a genome of 79,928 bp consisting of 125 protein coding genes and displaying low sequence similarities to public L. plantarum phages. Physiochemical characterization shows that it has a short latent period and adapts to broad ranges of temperatures and pHs. Furthermore, Gut-P1 strongly inhibits the growth of L. plantarum strains at a multiplicity of infection (MOI) of 1e-6. Together, these results indicate that Gut-P1 can greatly impede the application of L. plantarum in humans. Strikingly, Gut-P1 was identified only in the enrichment culture, not in our metagenomic or VLP sequencing data nor in any public human phage databases, indicating the inefficiency of bulk sequencing in recovering low-abundance but highly prevalent phages and pointing to the unexplored hidden diversity of the human gut virome despite recent large-scale sequencing and bioinformatics efforts. IMPORTANCE As Lactiplantibacillus plantarum (previously known as Lactobacillus plantarum) is increasingly used as a probiotic to treat human gut-related diseases, its bacteriophages may pose a certain threat to their further application and should be identified and characterized more often from the human intestine. Here, we isolated and identified the first gut L. plantarum phage that is prevalent in a Chinese population. This phage, Gut-P1, is virulent and can strongly inhibit the growth of multiple L. plantarum strains at low MOIs. Our results also show that bulk sequencing is inefficient at recovering low-abundance but highly prevalent phages such as Gut-P1, suggesting that the hidden diversity of human enteroviruses has not yet been explored. Our results call for innovative approaches to isolate and identify intestinal phages from the human gut and to rethink our current understanding of the enterovirus, particularly its underestimated diversity and overestimated individual specificity.

Complete genome sequence analysis, morphology and structural protein identification of two Bacillus subtilis phages, BSTP4 and BSTP6, which may form a new species in the genus Salasvirus

In the present study, two new Bacillus subtilis phages, BSTP4 and BSTP6, were isolated and studied further. Morphologically, BSTP4 and BSTP6 are podoviruses with complete genome of 19,145 (39.9% G + C content) and 19,367 bp (39.8% G + C content), respectively, which became among the smallest Bacillus phages. Three most prominent structural proteins were separated and identified as pre-neck appendage, major head, and head fiber proteins using LC-MS/MS. Both phages encode putative terminal proteins (TP) and contain short inverted terminal repeats (ITRs) which may be important for their replication. In addition, non-coding RNA (pRNA) and parS sites were identified which may be required for DNA packaging and their maintenance inside the host, respectively. Furthermore, the phage genome sequences show significant similarity to B. subtilis group species genome sequences. Finally, phylogenomic and phylogenetic analyses suggest that BSTP4 and BSTP6 may form a new species in the genus Salasvirus, subfamily Picovirinae of family Salasmaviridae. Considering the small numbers of ICTV-accepted B. subtilis phages and the importance of the host in the food industry and biotechnology, the current study helps to improve our understanding of the diversity of B. subtilis phages and shed light on the phage-host relationships.

Multicenter Preparedness Exercise Enables Rapid Development of Cluster-Specific PCR-Based Screening Assays from Bacterial Genomic Data

PCR-based screening assays targeting strain-specific genetic markers allow the timely detection and specific differentiation of bacterial strains. Especially in situations where an infection cluster occurs, fast assay development is crucial for supporting targeted control measures. However, the turnaround times (TATs) for assay setup may be high due to insufficient knowledge about screening assay methods, workflows, and software tools. Here, two blind-coded and quality-controlled ring trials were performed in which five German laboratories established PCR-based screening assays from genomic data that specifically target selected bacterial clusters within two bacterial monospecies sample panels. While the first ring trial was conducted without a time limit to train the participants and assess assay feasibility, in the second ring trial, a challenging time limit of 2 weeks was set to force fast assay development as soon as genomic data were available. During both ring trials, we detected high interlaboratory variability regarding the screening assay methods and targets, the TATs for assay setup, and the number of screening assays. The participants designed between one and four assays per cluster that targeted cluster-specific unique genetic sequences, genes, or single nucleotide variants using conventional PCRs, high-resolution melting assays, or TaqMan PCRs. Assays were established within the 2-week time limit, with TATs ranging from 4 to 13 days. TaqMan probe delivery times strongly influenced TATs. In summary, we demonstrate that a specific exercise improved the preparedness to develop functional cluster-specific PCR-based screening assays from bacterial genomic data. Furthermore, the parallel development of several assays enhances assay availability.

Genome, biology and stability of the Thurquoise phage – A new virus from the Bastillevirinae subfamily

Bacteriophages from the Bastillevirinae subfamily (Herelleviridae family) have proven to be effective against bacteria from the Bacillus genus including organisms from the B. cereus group, which cause food poisoning and persistent contamination of industrial installations. However, successful application of these phages in biocontrol depends on understanding of their biology and stability in different environments. In this study, we isolated a novel virus from garden soil in Wrocław (Poland) and named it ‘Thurquoise’. The genome of that phage was sequenced and assembled into a single continuous contig with 226 predicted protein-coding genes and 18 tRNAs. The cryo-electron microscopy revealed that Thurquoise has complex virion structure typical for the Bastillevirinae family. Confirmed hosts include selected bacteria from the Bacillus cereus group–specifically B. thuringiensis (isolation host) and B. mycoides, but susceptible strains display different efficiency of plating (EOP). The eclipse and latent periods of Thurquoise in the isolation host last ~ 50 min and ~ 70 min, respectively. The phage remains viable for more than 8 weeks in variants of the SM buffer with magnesium, calcium, caesium, manganese or potassium and can withstand numerous freeze–thaw cycles if protected by the addition of 15% glycerol or, to a lesser extent, 2% gelatine. Thus, with proper buffer formulation, this virus can be safely stored in common freezers and refrigerators for a considerable time. The Thurquoise phage is the exemplar of a new candidate species within the Caeruleovirus genus in the Bastillevirinae subfamily of the Herelleviridae family with a genome, morphology and biology typical for these taxa.

Special Issue “Bacteriophage Genomics”: Editorial

Virus genomics as a separate branch of biology has emerged relatively recently [...]

Comparative genomic analysis of the Dietzia genus: an insight into genomic diversity, and adaptation

Dietzia strains are widely distributed in the environment, presenting an opportunistic role, and some species have undetermined taxonomic characteristics. Here, we propose the existence of errors in the classification of species in this genus using comparative genomics. We performed ANI, dDDH, pangenome and genomic plasticity analyses better to elucidate the phylogenomic relationships between Dietzia strains. For this, we used 55 genomes of Dietzia downloaded from public databases that were combined with a newly sequenced. Sequence analysis of a phylogenetic tree based on genome similarity comparisons and dDDH, ANI analyses supported grouping different Dietzia species into four distinct groups. The pangenome analysis corroborated the classification of these groups, supporting the idea that some species of Dietzia could be reassigned in a possible classification into three distinct species, each containing less variability than that found within the global pangenome of all strains. Additionally, analysis of genomic plasticity based on groups containing Dietzia strains found differences in the presence and absence of symbiotic Islands and pathogenic islands related to their isolation site. We propose that the comparison of pangenome subsets together with phylogenomic approaches can be used as an alternative for the classification and differentiation of new species of the genus Dietzia.

Clostridioides difficile S-Layer Protein A (SlpA) Serves as a General Phage Receptor

Therapeutic bacteriophages (phages) are being considered as alternatives in the fight against Clostridioides difficile infections. To be efficient, phages should have a wide host range, buthe lack of knowledge about the cell receptor used by C. difficile phages hampers the rational design of phage cocktails. Recent reports suggested that the C. difficile surface layer protein A (SlpA) is an important phage receptor, but available data are still limited. Here, using the epidemic R20291 strain and its FM2.5 mutant derivative lacking a functional S-layer, we show that the absence of SlpA renders cells completely resistant to infection by ϕCD38-2, ϕCD111, and ϕCD146, which normally infect the parental strain. Complementation with 12 different S-layer cassette types (SLCTs) expressed from a plasmid revealed that SLCT-6 also allowed infection by ϕCD111 and SLCT-11 enabled infection by ϕCD38-2 and ϕCD146. Of note, the expression of SLCT-1, -6, -8, -9, -10, or -12 conferred susceptibility to infection by 5 myophages that normally do not infect the R20291 strain. Also, deletion of the D2 domain within the low-molecular-weight fragment of SlpA was found to abolish infection by ϕCD38-2 and ϕCD146 but not ϕCD111. Altogether, our data suggest that many phages use SlpA as their receptor and, most importantly, that both siphophages and myophages target SlpA despite major differences in their tail structures. Our study therefore represents an important step in understanding the interactions between C. difficile and its phages. IMPORTANCE Phage therapy represents an interesting alternative to treat Clostridioides difficile infections because, contrary to antibiotics, most phages are highly species specific, thereby sparing the beneficial gut microbes that protect from infection. However, currently available phages against C. difficile have a narrow host range and target members from only one or a few PCR ribotypes. Without a clear comprehension of the factors that define host specificity, and in particular the host receptor recognized by phages, it is hard to develop therapeutic cocktails in a rational manner. In our study, we provide clear and unambiguous experimental evidence that SlpA is a common receptor used by many siphophages and myophages. Although work is still needed to define how a particular phage receptor-binding protein binds to a specific SLCT, the identification of SlpA as a common receptor is a major keystone that will facilitate the rational design of therapeutic phage cocktails against clinically important strains.

Lactic Bacteria in Artisanal Cheese: Characterization through Metagenomics

Artisanal cheese, produced with raw milk by a predominantly manual approach, has a historical and cultural tradition related to the region of origin. Given its economic and cultural importance, the main objective of this study was to investigate and characterize the diversity of lactic acid bacteria (LAB) of artisanal cheeses produced and traded by family agro-industries in a region of southern Brazil. The LAB composition of artisanal cheese samples, belonging to different municipalities of the Region of Vale do Taquari, were characterized by the next-generation sequencing (NGS) method, amplifying the V3/V4 region of the 16S rRNA gene. A total of 35 LAB species, distributed in seven genera, were identified, and rarefaction analysis suggested that the total diversity assessed by 16S rRNA analysis was high in the analyzed samples. The average Ph ranged from 4.6 to 6.6, and a correlation with the genus Lactococcus (r = 0.62) was the most expressive. The LAB genera identified in the cheese samples were Bavariicococcus, Enterococcus, Lactobacillus, Lactococcus, Leuconostoc, Marinillactibacillus, and Pediococcus. Lactococcus lactis was the most predominant species, present in all samples. Although some species have been identified in the three altitudes studied, the abundance varied according to geographic environments. Enterococcus italicus is more present at high altitudes, unlike Lactococcus plantarum and Lactococcus raffinolactis at low altitudes. Lactococcus lactis was present in the three geographic environments evaluated, but the highest abundance was observed at high altitudes. The identification of LAB present in fermented cheeses is essential to understand the organoleptic quality during the maturation process as well as to establish the shelf life, including the safety and the overall quality of the cheese. This specific microbiota contributes to the flavor and unique characteristics of the regional dairy products, and on the other hand can be a source of specific starter cultures that guarantee the product’s identity.

Temperate bacteriophages infecting the mucin-degrading bacterium Ruminococcus gnavus from the human gut

Ruminococcus gnavus is a prevalent gut microbe reported to occur in higher abundance among individuals with inflammatory bowel disease (IBD). This study reports the isolation and characterization of six bacteriophages (phages) isolated from human fecal material and environmental samples that infect this species. Isolated phages have a siphovirus morphology, with genomes ranging between 36.5 and 37.8 kbp. Genome analysis indicates that the phages have a temperate lifestyle, which was confirmed by their ability to form lysogens on their host bacterial species. In contrast to the finding that phages lyse their host in liquid medium, results from a mouse trial indicate these phages can co-exist with the host bacterium in the gut without causing a significant reduction of R. gnavus. The bacterial counts in the feces of phage-treated mice did not significantly differ in the presence of phage. Furthermore, analysis of publicly available gut virome sequence data indicates a high abundance of these phages among individuals suffering from IBD. This work provides the first insight into how phages interact with R. gnavus in the human gut microbiome.

Single Nucleotide Polymorphism-Based Real-Time PCR Screening Assay for Rapid Tracking of Bacterial Infection Clusters To Complement Whole-Genome Sequencing Efforts during Outbreak Investigations

Infection clusters of multidrug-resistant bacteria increase mortality and entail expensive infection control measures. Whereas whole-genome sequencing (WGS) is the current gold standard to confirm infection clusters, PCR-based assays targeting cluster-specific signatures, such as single nucleotide polymorphisms (SNPs) derived from WGS data, are more suitable to initially screen for cluster isolates within large sample sizes. Here, we evaluated four software tools (SeqSphere+, RUCS, Gegenees, and Find Differential Primers) regarding their efficiency to find SNPs within WGS data sets that were specific for two bacterial monospecies infection clusters but were absent from a WGS reference data set comprising several hundred diverse genotypes of the same bacterial species. Cluster-specific SNPs were subsequently used to establish a probe-based real-time PCR screening assay for in vitro differentiation between cluster and noncluster isolates. SeqSphere+ and RUCS found 2 and 24 SNPs for clusters 1 and 14 and 24 SNPs for cluster 2, respectively. However, some signatures detected by RUCS were not cluster specific. Interestingly, all SNPs identified by SeqSphere+ were also detected by RUCS. In contrast, analyses with the remaining tools either resulted in no SNPs (with Find Differential Primers) or failed (Gegenees). Design of six cluster-specific real-time PCR assays enabled reliable cluster screening in vitro. Our evaluation revealed that SeqSphere+ and RUCS identified cluster-specific SNPs that could be used for large-scale screening in surveillance samples via real-time PCR, thereby complementing WGS efforts. This faster and simplified approach for the surveillance of bacterial clusters will improve infection control measures and will enhance protection of patients and physicians. IMPORTANCE Infection clusters of multidrug-resistant bacteria threaten medical facilities worldwide and cause immense health care costs. In recent years, whole-genome sequencing (WGS) has been increasingly applied to detect and to further control bacterial clusters. However, as WGS is still expensive and time-consuming, its exclusive application for screening and confirmation of bacterial infection clusters contributes to high costs and enhanced turnaround times, which many hospitals cannot afford. Therefore, there is need for alternative methods that can enable further surveillance of bacterial clusters that are initially detected by WGS in a faster and more cost-efficient way. Here, we established a system based on real-time PCR that enables rapid large-scale sample screening for bacterial cluster isolates within 7 days after the initial detection of an infection cluster, thereby complementing WGS efforts. This faster and simplified surveillance of bacterial clusters will improve infection control measures and will enhance protection of patients and physicians.

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.

Genomic diversity and comprehensive taxonomical classification of 61 Bacillus subtilis group member infecting bacteriophages, and the identification of ortholog taxonomic signature genes

BACKGROUND

Despite the applications of Bacillus subtilis group species in various sectors, limited information is available regarding their phages. Here, 61 B. subtilis group species-infecting phages (BSPs) were studied for their taxonomic classification considering the genome-size, genomic diversity, and the host, followed by the identification of orthologs taxonomic signature genes.

RESULTS

BSPs have widely ranging genome sizes that can be bunched into groups to demonstrate correlations to family and subfamily classifications. Comparative analysis re-confirmed the existing, BSPs-containing 14 genera and 21 species and displayed inter-genera similarities within existing subfamilies. Importantly, it also revealed the need for the creation of new taxonomic classifications, including 28 species, nine genera, and two subfamilies (New subfamily1 and New subfamily2) to accommodate inter-genera relatedness. Following pangenome analysis, no ortholog shared by all BSPs was identified, while orthologs, namely, the tail fibers/spike proteins and poly-gamma-glutamate hydrolase, that are shared by more than two-thirds of the BSPs were identified. More importantly, major capsid protein (MCP) type I, MCP type II, MCP type III and peptidoglycan binding proteins that are distinctive orthologs for Herelleviridae, Salasmaviridae, New subfamily1, and New subfamily2, respectively, were identified and analyzed which could serve as signatures to distinguish BSP members of the respective taxon.

CONCLUSIONS

In this study, we show the genomic diversity and propose a comprehensive classification of 61 BSPs, including the proposition for the creation of two new subfamilies, followed by the identification of orthologs taxonomic signature genes, potentially contributing to phage taxonomy.

Genomic and Phenotypic Evaluation of Potential Probiotic Pediococcus Strains with Hypocholesterolemic Effect Isolated from Traditional Fermented Food

The use of probiotic microorganisms in food with the aim to confer health benefits to the host is one of the most critical roles of functional foods. Many pediococci bacteria frequently related to the meat environment, have technological properties, and are therefore commercially used as starter in the production of fermented meat products, such as different types of sausages. In this study, different lactic acid bacteria were isolated, identified to the species level, and then evaluated for their safety and functionality as possible probiotics. Different properties, such as resistance to simulated human gastrointestinal conditions, antimicrobial activity, and cholesterol-lowering effects, have been studied. Finally, the complete genome of one strain, namely P. acidilactici IRZ12B, which showed interesting features as a promising probiotic candidate, was sequenced and further studied. The results revealed that IRZ12B possesses interesting probiotic properties, particularly cholesterol-lowering capability and antimicrobial activity. In silico analysis evidenced the absence of plasmids, transmissible antibiotic resistance genes, and virulence factors. We also detected a bacteriocin encoding gene and a cholesterol assimilation-related protein. The phenotypical and genomic outcomes described in this study make P. acidilactici IRZ12B a very interesting cholesterol-lowering potential probiotic strain to be considered for the development of novel non-dairy-based functional foods.

Comparative polyphasic characterization of Weissella strains isolated from beaked whale and rainbow trout (Oncorhynchus mykiss): confirmation of Weissella ceti sp. nov. and description of the novel Weissella tructae sp. nov. isolated from farmed rainbow trout

The weissellosis agent bacterium (WS08T = CBMAI 2730) was isolated from diseased rainbow trout (Oncorhynchus mykiss) in Brazil. The whole genome sequence of this strain was compared with the Mexican W-1 strain, also isolated from diseased rainbow trout, and with the Weissella ceti type strain CECT 7719 T (= 1119-1A-09 T = CCUG 59653 T), recovered from the beaked whale. Digital DNA-DNA hybridization pairwise analyses scored 98.7% between the Mexican W-1 and Brazilian WS08T but just 24.4% for both fish isolates compared to the W. ceti type strain CECT 7719 T. The 16S rRNA gene sequence comparisons with isolates of W. ceti, available at GenBank, were conducted. All rainbow trout-pathogenic isolates grouped close (97% bootstrap confirmation), but when this group was compared to the W. ceti type strain CECT 7719 T the similarity varied from 78.9 to 79.1%. Phenotypic assays were also conducted, and the W. ceti type strain diverged from WS08T and W-1 in the hydrolysis of aesculin, D-mannose, and potassium gluconate and in the hydrolysis of hippurate. Moreover, WS08T and W-1 showed weak growth at 5 °C whereas no growth was observed for W. ceti CECT 7719 T. The major fatty acids (> 10% total fatty acids) presented by WS08T and W-1 were summed feature 8 (C18:1 ω7c/C18:1 ω6c), summed feature 3 (C16:1 ω6c/C16:1ω7c), and C16:0. The results of phylogenetic and phenotypic analyses clearly differentiated the W. ceti CECT 7719 T type strain from the assessed pathogenic strains obtained from rainbow trout. Therefore, Weissella strains isolated from rainbow trout, here represented by strain WS08T (= CBMAI 2730), should be known as members of a novel species for which the name Weissella tructae sp. nov. is proposed.

Genome Mining Reveals High Biosynthetic Potential of Biocontrol Agent Bacillus velezensis B.BV10

The present study demonstrates the biocontrol potential of a plant growth-promoting bacterial strain using three different approaches: (i) an in vitro evaluation of antagonistic activity against important phytopathogenic fungi; (ii) an evaluation under greenhouse conditions with strawberry plants to assess the control of gray mold; and (iii) an in silico whole genome sequence mining to assign genetic features such as gene clusters or isolated genes to the strain activity. The in vitro assay showed that the B.BV10 strain presented antagonistic activity, inhibiting the mycelial growth in all the phytopathogenic fungi evaluated. The application of the Bacillus velezensis strain B.BV10 under greenhouse conditions reduced the presence of Botrytis cinerea and increased the mean fruit biomass. The genome of B.BV10 was estimated at 3,917,533 bp, with a GC content of 46.6% and 4088 coding DNA sequences, and was identified as B. velezensis. Biosynthetic gene clusters related to the synthesis of the molecules with antifungal activity were found in its genome. Genes related to the regulation/formation of biofilms, motility, and the important properties for the rhizospheric colonization were also found in the genome. The current study offers a comprehensive understanding of the genomic architecture and control activity of phytopathogenic fungi by the B. velezensis strain B.BV10 that may substantiate the industrialization of this strain in the future.

Antifungal activity and genomic characterization of the biocontrol agent Bacillus velezensis CMRP 4489

The development of bio-based products has increased in recent years, and species of the Bacillus genus have been widely used for product development due to their elevated production of antimicrobial molecules and resistance to extreme environmental conditions through endospore formation. In this context, the antifungal potential of Bacillus velezensis CMRP 4489 was investigated using in silico predictions of secondary metabolites in its genome and in vitro tests against the following phytopathogenic fungi: Sclerotinia sclerotiorum, Macrophomina phaseolina, and Botrytis cinerea. The in-silico predictions indicated that CMRP 4489 possesses several Biosynthetic Gene Clusters (BGCs) capable of producing molecules with antifungal properties and other non-identified BGCs. The in vitro assay results evidenced strong antifungal activity, inhibiting more than 60% of the tested fungi, and the isolate's molecules were stable under diverse physicochemical conditions. The in vitro assay evidenced significant antifungal activity, deformation of the hyphal structure in SS, biofilm formation capacity, and swarming motility. In the colonization assay, we observed attachment, colonization, and net-shaped biofilm formation, with the strain transitioning from the seeds to nearby structures. Therefore, CMRP 4489 showed to be a potential biocontrol agent against various diseases with agronomic importance and can be used under adverse environmental conditions.

Bioinformatic Analysis of a Set of 14 Temperate Bacteriophages Isolated from Staphylococcus aureus Strains Highlights Their Massive Genetic Diversity

Epidemiology and virulence studies of Staphylococcus aureus showed that temperate bacteriophages are one of the most powerful drivers for its evolution not only because of their abundance but also because of the richness of their genetic payload. Here, we report the isolation, genome sequencing, and bioinformatic analysis of 14 bacteriophages induced from lysogenic S. aureus strains from human or veterinary (cattle) origin. The bacteriophages belonged to the Siphoviridae family; were of similar genome size (40 to 45 kbp); and fell into clusters B2, B3, B5, and B7 according to a recent clustering proposal. One of the phages, namely, vB_SauS_308, was the most unusual one, belonging to the sparsely populated subcluster B7 but showing differences in protein family contents compared with the rest of the members. This phage contains a type I endolysin (one catalytic domain and noncanonical cell wall domain [CBD]) and a host recognition module lacking receptor binding protein, cell wall hydrolase, and tail fiber proteins. This phage also lacked virulence genes, which is opposite to what has been reported for subcluster B6 and B7 members. None of six phages, taken as representatives of each of the four subclusters, showed activity on coagulase-negative staphylococci (excepted for two Staphylococcus hominis strains in which propagation and a very slow adsorption rate were observed) nor transducing ability. Immunity tests on S. aureus RN4220 lysogens with each of these phages showed no cross immunity. IMPORTANCE To the best of our knowledge, this set of sequenced bacteriophages is the largest one in South America. Our report describes for the first time the utilization of MultiTwin software to analyze the relationship between phage protein families. Notwithstanding the fact that most of the genetic information obtained correlated with recently published information, due to their geographical origin, the reported analysis adds up to and confirms currently available knowledge of Staphylococcus aureus temperate bacteriophages in terms of phylogeny and role in host evolution.

Impact of a phage cocktail targeting Escherichia coli and Enterococcus faecalis as members of a gut bacterial consortium in vitro and in vivo

Escherichia coli and Enterococcus faecalis have been implicated as important players in human gut health that have been associated with the onset of inflammatory bowel disease (IBD). Bacteriophage (phage) therapy has been used for decades to target pathogens as an alternative to antibiotics, but the ability of phage to shape complex bacterial consortia in the lower gastrointestinal tract is not clearly understood. We administered a cocktail of six phages (either viable or heat-inactivated) targeting pro-inflammatory Escherichia coli LF82 and Enterococcus faecalis OG1RF as members of a defined community in both a continuous fermenter and a murine colitis model. The two target strains were members of a six species simplified human microbiome consortium (SIHUMI-6). In a 72-h continuous fermentation, the phage cocktail caused a 1.1 and 1.5 log (log10 genome copies/mL) reduction in E. faecalis and E. coli numbers, respectively. This interaction was accompanied by changes in the numbers of other SIHUMI-6 members, with an increase of Lactiplantibacillus plantarum (1.7 log) and Faecalibacterium prausnitzii (1.8 log). However, in germ-free mice colonized by the same bacterial consortium, the same phage cocktail administered twice a week over nine weeks did not cause a significant reduction of the target strains. Mice treated with active or inactive phage had similar levels of pro-inflammatory cytokines (IFN-y/IL12p40) in unstimulated colorectal colonic strip cultures. However, histology scores of the murine lower GIT (cecum and distal colon) were lower in the viable phage-treated mice, suggesting that the phage cocktail did influence the functionality of the SIHUMI-6 consortium. For this study, we conclude that the observed potential of phages to reduce host populations in in vitro models did not translate to a similar outcome in an in vivo setting, with this effect likely brought about by the reduction of phage numbers during transit of the mouse GIT.

Efficiency of Combining Strains Ag87 (Bacillus megaterium) and Ag94 (Lysinibacillus sp.) as Phosphate Solubilizers and Growth Promoters in Maize

Increasing phosphorus (P) use efficiency in agricultural systems is urgent and essential to significantly reduce the global demand for this nutrient. Applying phosphate-solubilizing and plant growth-promoting bacteria in the rhizosphere represents a strategy worthy of attention. In this context, the present work aimed to select and validate bacterial strains capable of solubilizing phosphorous and promoting maize growth, aiming to develop a microbial inoculant to be used in Brazilian agriculture. Bacterial strains from the maize rhizosphere were evaluated based on their ability to solubilize phosphate and produce indole acetic acid. Based on these characteristics, 24 strains were selected to be further evaluated under laboratory, greenhouse, and field conditions. Among the selected strains, four (I04, I12, I13, and I17) showed a high potential to increase maize root growth and shoot P content. Strains I13 (Ag87) and I17 (Ag94) were identified by genomic sequencing as Bacillus megaterium and Lysinibacillus sp., respectively. These strains presented superior yield increments relative to the control treatment with 30% P. In addition, combining Ag87 and Ag94 resulted in even higher yield gains, indicating a synergistic effect that could be harnessed in a commercial inoculant for Brazilian agriculture.

Multiplex PCR assay for correct identification of the fish pathogenic species of Edwardsiella genus reveals the presence of E. anguillarum in South America in strains previously characterized as E. tarda

AIMS

Develop a species-specific multiplex PCR to correctly identify Edwardsiella species in routine diagnostic for fish bacterial diseases.

METHODS AND RESULTS

The genomes of 62 Edwardsiella spp. isolates available from the National Center for Biotechnology Information (NCBI) database were subjected to taxonomic and pan-genomic analyses to identify unique regions that could be exploited by species-specific PCR. The designed primers were tested against isolated Edwardsiella spp. strains, revealing errors in commercial biochemical tests for bacterial classification regarding Edwardsiella species.

CONCLUSION

Some of the genomes of Edwardsiella spp. in the NCBI platform were incorrectly classified, which can lead to errors in some research. A functional mPCR was developed to differentiate between phenotypically and genetically ambiguous Edwardsiella, with which, we detected the presence of Edwardsiella anguillarum affecting fish in Brazil.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study shows that the misclassification of Edwardsiella spp in Brazil concealed the presence of E. anguillarum in South America. Also, this review of the taxonomic classification of the Edwardsiella genus is a contribution to the field to help researchers with their sequencing and identification of genomes, showing some misclassifications in online databases that must be corrected, as well as developing an easy assay to characterize Edwardsiella species in an end-point mPCR.

Prevalence, Diversity and UV-Light Inducibility Potential of Prophages in Bacillus subtilis and Their Possible Roles in Host Properties

Bacillus subtilis is an important bacterial species due to its various industrial, medicinal, and agricultural applications. Prophages are known to play vital roles in host properties. Nevertheless, studies on the prophages and temperate phages of B. subtilis are relatively limited. In the present study, an in silico analysis was carried out in sequenced B. subtilis strains to investigate their prevalence, diversity, insertion sites, and potential roles. In addition, the potential for UV induction and prevalence was investigated. The in silico prophage analysis of 164 genomes of B. subtilis strains revealed that 75.00% of them contained intact prophages that exist as integrated and/or plasmid forms. Comparative genomics revealed the rich diversity of the prophages distributed in 13 main clusters and four distinct singletons. The analysis of the putative prophage proteins indicated the involvement of prophages in encoding the proteins linked to the immunity, bacteriocin production, sporulation, arsenate, and arsenite resistance of the host, enhancing its adaptability to diverse environments. An induction study in 91 B. subtilis collections demonstrated that UV-light treatment was instrumental in producing infective phages in 18.68% of them, showing a wide range of host specificity. The high prevalence and inducibility potential of the prophages observed in this study implies that prophages may play vital roles in the B. subtilis host.

Curtobacterium spp. and Curtobacterium flaccumfaciens: Phylogeny, Genomics-Based Taxonomy, Pathogenicity, and Diagnostics

The genus of Curtobacterium, belonging to the Microbacteriaceae family of the Actinomycetales order, includes economically significant pathogenic bacteria of soybeans and other agricultural crops. Thorough phylogenetic and full-genome analysis using the latest genomic data has demonstrated a complex and contradictory taxonomic picture within the group of organisms classified as the Curtobacterium species. Based on these data, it is possible to delineate about 50 new species and to reclassify a substantial part of the Curtobacterium strains. It is suggested that 53 strains, including most of the Curtobacterium flaccumfaciens pathovars, can compose a monophyletic group classified as C. flaccumfaciens. A genomic analysis using the most recent inventory of bacterial chromosomal and plasmid genomes deposited to GenBank confirmed the possible role of Microbacteriaceae plasmids in pathogenicity and demonstrated the existence of a group of related plasmids carrying virulence factors and possessing a gene distantly related to DNA polymerase found in bacteriophages and archaeal and eukaryotic viruses. A PCR diagnostic assay specific to the genus Curtobacterium was developed and tested. The presented results assist in the understanding of the evolutionary relations within the genus and can lay the foundation for further taxonomic updates.

Selective Isolation of Eggerthella lenta from Human Faeces and Characterisation of the Species Prophage Diversity

Eggerthella lenta is an anaerobic, high GC, Gram-positive bacillus commonly found in the human digestive tract that belongs to the class Coriobacteriia of the phylum Actinobacteria. This species has been of increasing interest as an important player in the metabolism of xenobiotics and dietary compounds. However, little is known regarding its susceptibility to bacteriophage predation and how this may influence its fitness. Here, we report the isolation of seven novel E. lenta strains using cefotaxime and ceftriaxone as selective agents. We conducted comparative and pangenome analyses of these strains and those publicly available to investigate the diversity of prophages associated with this species. Prophage gene products represent a minimum of 5.8% of the E. lenta pangenome, comprising at least ten distantly related prophage clades that display limited homology to currently known bacteriophages. All clades possess genes implicated in virion structure, lysis, lysogeny and, to a limited extent, DNA replication. Some prophages utilise tyrosine recombinases and diversity generating retroelements to generate phase variation among targeted genes. The prophages have differing levels of sensitivity to the CRISPR/cas systems of their hosts, with spacers from 44 E. lenta isolates found to target only five out of the ten identified prophage clades. Furthermore, using a PCR-based approach targeting the prophage attP site, we were able to determine that several of these elements can excise from the host chromosome, thus supporting the notion that these are active prophages. The findings of this study provide further insights into the diversity of prophages infecting species of the phylum Actinobacteria.

Isolation, characterization and biocontrol efficacy of a T4-like phage virulent to multidrug-resistant Enterobacter hormaechei

Enterobacter hormaechei is an important emerging pathogen, often exhibiting resistance to multiple clinically important antibiotics. In this study, E. hormaechei was found, for the first time, to be lethal to fish. Bacteriophages are considered potential treatments for bacterial infections. The lytic phage vB_EhoM-IME523 (abbreviated 'IME523') infecting multidrug-resistant E. hormaechei was isolated from hospital sewage. IME523 exhibits T4-like morphology, including a prolate icosahedral head 110 ± 1.89 nm (mean ± SD) long and 82 ± 0.75 nm wide, and a contractile tail of ca. 110 ± 0.91 nm in length. The complete genome length of phage IME523 is 172763 bp, with a G + C content of 39.97%. The whole genome sequence of IME523 has a 93.10% average nucleotide identity (ANI) and a 53.3% in silico DNA-DNA hybridization (isDDH) value with the closest-related Enterobacter phage vB_EclM_CIP9 ('CIP9'). ANI and isDDH values between IME523 and other phages were lower than 78 and 22%, respectively. IME523 and CIP9 formed a monophyletic branch in a phylogenetic tree based on the terminase large subunit, DNA polymerase protein and whole genome phylogenetic analysis. Results suggest that IME523 is a novel species in the subfamily Tevenvirinae and forms a novel genus together with CIP9. No IME523 open reading frame was found to be associated with virulence factors or antibiotic resistance genes. IME523 showed promising protection to zebrafish and brocade carp against E. hormaechei challenge.

Genomic and Phenotypic Characteristics in Geographically Separated Clinical Campylobacter jejuni ST353CC Isolates

Campylobacter jejuni fecal isolates of eight international travelers, 5 of which had traveled to Ecuador and 3 to Bangladesh, were characterized, and the possible relationship between bacterial traits and clinical symptoms was further analyzed. All eight isolates belonged to the same Multi-Locus Sequence Type clonal complex (ST353CC). The three isolates from Bangladesh were all of the same sequence type (ST-9438), and when compared to isolates of various other sequence types, they had a larger quantity of unique genetic content, higher expression levels of some putative virulence genes involved in adhesion and invasion (flpA, ciaB and iamA), and showed higher adhesion levels to human HT-29 colon cancer cells in an in vitro infection model. However, in contrast to the seemingly higher pathogenic potential of these bacterial isolates, travelers infected with the ST-9438 isolates had no or only very mild symptoms, whereas the other individuals, whose bacterial isolates seemed to have less pathogenic potential, generally reported severe symptoms. When studying the 16S rRNA gene-based fecal microbiota in samples collected prior to travel, there was an individual variation in the relative abundance of the three major bacterial phyla Actinobacteria, Bacteroidetes and Firmicutes, but there were no associations between composition and diversity of microbiota and development of severe symptoms from the infection. It remains to be confirmed by larger studies whether an individual's characteristics such as gut microbiota, might be related to the severity of symptoms in Campylobacter infections.

Extensive Genome Exploration of Clostridium botulinum Group III Field Strains

In animals, botulism is commonly sustained by botulinum neurotoxin C, D or their mosaic variants, which are produced by anaerobic bacteria included in Clostridium botulinum group III. In this study, a WGS has been applied to a large collection of C. botulinum group III field strains in order to expand the knowledge on these BoNT-producing Clostridia and to evaluate the potentiality of this method for epidemiological investigations. Sixty field strains were submitted to WGS, and the results were analyzed with respect to epidemiological information and compared to published sequences. The strains were isolated from biological or environmental samples collected in animal botulism outbreaks which occurred in Italy from 2007 to 2016. The new sequenced strains belonged to subspecific groups, some of which were already defined, while others were newly characterized, peculiar to Italian strains and contained genomic features not yet observed. This included, in particular, two new flicC types (VI and VII) and new plasmids which widen the known plasmidome of the species. The extensive genome exploration shown in this study improves the C. botulinum and related species classification scheme, enriching it with new strains of rare genotypes and permitting the highest grade of discrimination among strains for forensic and epidemiological applications.

Effect of Enterococcus faecium as a Water and/or Feed Additive on the Gut Microbiota, Hematologic and Immunological Parameters, and Resistance Against Francisellosis and Streptococcosis in Nile Tilapia (Oreochromis niloticus)

In the present study, we evaluated the effects of administering Enterococcus faecium in food and/or water on the hematological and immunological parameters, intestinal microbiota, resistance to bacterial diseases (streptococcosis and francisellosis) and growth of Nile tilapia. Before the in vivo experiment, probiotic bacteria isolated from Nile tilapia were selected via inhibition tests. Sequencing, annotation, and assembly of the complete genome of the selected bacteria as well as other tests were performed using bioinformatics tools. Three treatments were implemented: G1 (probiotic feeding), G2 (probiotic in water), and G3 (probiotic in food and water); and a negative control (NC) was also employed. Treatment lasted 38 days, and each group consisted of fish and two repetitions. The fish were divided and infected with Streptococcus agalactiae S13 (serotype Ib) and Francisella orientalis. The G1 group had a higher average final weight gain than the G2, G3, and NC groups. Further, a significant increase in the number of thrombocytes was observed in the groups administered probiotics in the diet (G1 and G3). A statistical difference was observed in the mortality of fish infected with S. agalactiae in the NC compared to the treated groups. Cetobacterium was the 43 most abundant genus in the intestinal microbiota of all groups, including the NC group. E. faecium increased the immunity of fish administered the treatment and decreased the mortality caused by S. agalactiae. As an autochtone probiotic, E. faecium does not interfere with the local ecosystem and thus has a great probiotic potential for Nile tilapia in Brazil.

Exploration of the Diversity of Clustered Regularly Interspaced Short Palindromic Repeats-Cas Systems in Clostridium novyi sensu lato

Classified as the genospecies Clostridium novyi sensu lato and distributed into four lineages (I-IV), Clostridium botulinum (group III), Clostridium novyi, and Clostridium haemolyticum are clostridial pathogens that cause animal diseases. Clostridium novyi sensu lato contains a large mobilome consisting of plasmids and circular bacteriophages. Here, we explored clustered regularly interspaced short palindromic repeats (CRISPR) arrays and their associated proteins (Cas) to shed light on the link between evolution of CRISPR-Cas systems and the plasmid and phage composition in a study of 58 Clostridium novyi sensu lato genomes. In 55 of these genomes, types I-B (complete or partial), I-D, II-C, III-B, III-D, or V-U CRISPR-Cas systems were detected in chromosomes as well as in mobile genetic elements (MGEs). Type I-B predominated (67.2%) and was the only CRISPR type detected in the Ia, III, and IV genomic lineages. Putative type V-U CRISPR Cas14a genes were detected in two different cases: next to partial type-IB CRISPR loci on the phage encoding the botulinum neurotoxin (BoNT) in lineage Ia and in 12 lineage II genomes, as part of a putative integrative element related to a phage-inducible chromosomal island (PICI). In the putative PICI, Cas14a was associated with CRISPR arrays and restriction modification (RM) systems as part of an accessory locus. This is the first time a PICI containing such locus has been detected in C. botulinum. Mobilome composition and dynamics were also investigated based on the contents of the CRISPR arrays and the study of spacers. A large proportion of identified protospacers (20.2%) originated from Clostridium novyi sensu lato (p1_Cst, p4_BKT015925, p6_Cst, CWou-2020a, p1_BKT015925, and p2_BKT015925), confirming active exchanges within this genospecies and the key importance of specific MGEs in Clostridium novyi sensu lato.

Genomic analysis revealed conserved acid tolerance mechanisms from native micro-organisms in fermented feed

AIMS

Fermented feed is an agricultural practice used in many regions of the world to improve the growth performance of farm animals. This study aimed to identify and evaluate the lactic acid bacteria and yeast involved in the production of fermented feed.

METHODS AND RESULTS

We isolated and described two micro-organisms from autochthonous microbiota origin present in a regional feed product, Lactobacillus paracasei IBR07 (Lacticaseibacillus paracasei) and Kazachstania unispora IBR014 (Saccharomyces unisporum). Genome sequence analyses were performed to characterize both micro-organisms. Potential pathways involved in the acid response, tolerance and persistence were predicted in both genomes. Although L. paracasei and K. unispora are considered safe for animal feed, we analysed the presence of virulence factors, antibiotic resistance and pathogenicity islands. Furthermore, the Galleria mellonella model was used to support the safety of both isolates.

CONCLUSIONS

We conclude that IBR07 and IBR014 strains are good candidates to be used as starter cultures for feed fermentation.

SIGNIFICANCE AND IMPACT OF THE STUDY

The data presented here will be helpful to explore other biotechnological aspects and constitute a starting point for further studies to establish the consumption benefit of fermented feed in farm animal production.

Isolation and Characterization of Bacteriophages Infecting Burkholderia glumae, the Major Causal Agent of Bacterial Panicle Blight in Rice

Bacterial panicle blight (BPB), caused by Burkholderia glumae, is one of the most severe seed-borne bacterial diseases of rice in the world, which can decrease rice production by ≤75%. Nevertheless, there are few effective measures to manage this disease. In an attempt to develop an alternative management tool for BPB, we isolated and characterized phages from soil and water that are effective to lyse several strains of B. glumae. After tests of host ranges, the phages NBP1-1, NBP4-7, and NBP4-8 were selected for further comprehensive characterization, all of which could lyse B. glumae BGLa14-8 (phage sensitive) but not B. glumae 336gr-1 (phage insensitive). This result indicates that the phages killing B. glumae cells have specific host ranges at the strain level within the bacterial species. In the greenhouse condition of this study, foliar application of the phage NBP4-7 reduced the severity of BPB caused by B. glumae BGLa14-8 ≤62% but did not cause any significant effect on the infection by B. glumae 336gr-1. Electron microscopy and whole-genome sequencing were also performed to characterize the three selected phages. Transmission electron microscopy revealed that the selected phages belong to the family Myoviridae. Furthermore, whole-genome sequence analysis indicated that the three phages belong to a same species and are closely related to the Burkholderia phage KL3, a member of the Myoviridae family.

Tailed Lytic Bacteriophages of Soft Rot Pectobacteriaceae

The study of the ecological and evolutionary traits of Soft Rot Pectobacteriaceae (SRP) comprising genera Pectobacterium and Dickeya often involves bacterial viruses (bacteriophages). Bacteriophages are considered to be a prospective tool for the ecologically safe and highly specific protection of plants and harvests from bacterial diseases. Information concerning bacteriophages has been growing rapidly in recent years, and this has included new genomics-based principles of taxonomic distribution. In this review, we summarise the data on phages infecting Pectobacterium and Dickeya that are available in publications and genomic databases. The analysis highlights not only major genomic properties that assign phages to taxonomic families and genera, but also the features that make them potentially suitable for phage control applications. Specifically, there is a discussion of the molecular mechanisms of receptor recognition by the phages and problems concerning the evolution of phage-resistant mutants.

Phylogenomic Characterization of a Novel Corynebacterium Species Associated with Fatal Diphtheritic Stomatitis in Endangered Yellow-Eyed Penguins

Yellow-eyed penguins, Megadyptes antipodes, are an endangered species that are endemic to New Zealand. Outbreaks of diphtheritic stomatitis have caused significant mortality for this species, especially among young chicks. In this study, we isolated 16 Corynebacterium sp. isolates from the oral cavities of 2- to 14-day-old chicks at a range of infection stages and sequenced the genomes to understand their virulence mechanisms. Phylogenomic and matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) characterization indicate that these strains belong to a novel Corynebacterium species. A simple multiplex PCR-based diagnostic assay has been developed to identify these strains rapidly and reliably. Similar to other corynebacteria, genomic islands and prophages introduced significant diversity among these strains that has potentially led to minor functional variations between the two lineages. Despite the presence of multiple corynebacterial virulence genes and a spaDEF-type pilus gene cluster among these strains, the survival rate was much higher in Galleria mellonella larvae than in those inoculated with Corynebacterium ulcerans NZRM 818 and Corynebacterium pseudotuberculosis NZRM 3004. Therefore, these strains are opportunistic pathogens causing high mortality among young penguin chicks due to a less-developed immune system.

IMPORTANCE Yellow-eyed penguins, Megadyptes antipodes, are endangered species with a sharp decline in the numbers of breeding pairs over the last 2 decades. Diphtheritic stomatitis, characterized by a thick fibrinopurulent exudate in the oral cavities and symptoms, including inanition and significant weight loss, is responsible for significant mortality among the young chicks. These chicks are treated with antibiotics, amoxicillin-clavulanic acid or enrofloxacin, but do not always recover from the infection. The pathogen causing these infections and the mechanism of pathogenesis are unclear. This study has identified a novel Corynebacterium species to be associated with diphtheritic stomatitis in yellow-eyed penguins with potential virulence genes that are likely involved in pathogenesis. Importantly, a gene encoding an exotoxin, phospholipase D, is present among these strains. The inactivated form of this enzyme could potentially be used as an effective vaccine to protect these penguins from infection.

Closing Clostridium botulinum Group III Genomes Using Long-Read Sequencing

Clostridium botulinum group III is the anaerobic Gram-positive bacterium producing the deadly neurotoxin responsible for animal botulism. Here, we used long-read sequencing to produce four complete genomes from Clostridium botulinum group III neurotoxin types C, D, C/D, and D/C. The protocol for obtaining high-molecular-weight DNA from C. botulinum group III is described.

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.

Genomic Insights Into the Antifungal Activity and Plant Growth-Promoting Ability in Bacillus velezensis CMRP 4490

The main objective of this study was to evaluate Bacillus velezensis strain CMRP 4490 regarding its ability to inhibit soil-borne plant pathogens and to increase plant growth. The study included evaluation of in vitro antifungal control, sequencing the bacterial genome, mining genes responsible for the synthesis of secondary metabolites, root colonization ability, and greenhouse studies for the assessment of plant growth–promoting ability. The strain was obtained from soil samples in the north of Paraná in Brazil and was classified as a B. velezensis, which is considered a promising biological control agent. In vitro assay showed that B. velezensis CMRP 4490 presented antagonistic activity against Sclerotinia sclerotiorum, Macrophomina phaseolina, Botrytis cinerea, and Rhizoctonia solani with a mycelial growth inhibition of approximately 60%, without any significant difference among them. To well understand this strain and to validate its effect on growth-promoting rhizobacteria, it was decided to explore its genetic content through genome sequencing, in vitro, and greenhouse studies. The genome of CMRP 4490 was estimated at 3,996,396 bp with a GC content of 46.4% and presents 4,042 coding DNA sequences. Biosynthetic gene clusters related to the synthesis of molecules with antifungal activity were found in the genome. Genes linked to the regulation/formation of biofilms, motility, and important properties for rhizospheric colonization were also found in the genome. Application of CMRP 4490 as a coating film on soybean increased from 55.5 to 64% on germination rates when compared to the control; no differences were observed among treatments for the maize germination. The results indicated that B. velezensis CMRP 4490 could be a potential biocontrol agent with plant growth–promoting ability.

Antimicrobial activity screening of rhizosphere soil bacteria from tomato and genome-based analysis of their antimicrobial biosynthetic potential

Background

Tomato plant growth is frequently hampered by a high susceptibility to pests and diseases. Traditional chemical control causes a serious impact on both the environment and human health. Therefore, seeking environment-friendly and cost-effective green methods in agricultural production becomes crucial nowadays. Plant Growth Promoting Rhizobacteria (PGPR) can promote plant growth through biological activity. Their use is considered to be a promising sustainable approach for crop growth. Moreover, a vast number of biosynthetic gene clusters (BGCs) for secondary metabolite production are being revealed in PGPR, which helps to find potential anti-microbial activities for tomato disease control.

Results

We isolated 181 Bacillus-like strains from healthy tomato, rhizosphere soil, and tomato tissues. In vitro antagonistic assays revealed that 34 Bacillus strains have antimicrobial activity against Erwinia carotovora, Pseudomonas syringae; Rhizoctonia solani; Botrytis cinerea; Verticillium dahliae and Phytophthora infestans. The genomes of 10 Bacillus and Paenibacillus strains with good antagonistic activity were sequenced. Via genome mining approaches, we identified 120 BGCs encoding NRPs, PKs-NRPs, PKs, terpenes and bacteriocins, including known compounds such as fengycin, surfactin, bacillibactin, subtilin, etc. In addition, several novel BGCs were identified. We discovered that the NRPs and PKs-NRPs BGCs in Bacillus species are encoding highly conserved known compounds as well as various novel variants.

Conclusions

This study highlights the great number of varieties of BGCs in Bacillus strains. These findings pave the road for future usage of Bacillus strains as biocontrol agents for tomato disease control and are a resource arsenal for novel antimicrobial discovery.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07346-8.

Enterococcus phage Nonaheksakonda infecting clinical isolates of Enterococcus faecalis represents a new lineage in the family Siphoviridae

Enterococcus phage Nonaheksakonda was isolated from wastewater, using a vancomycin-resistant strain of the opportunistic pathogen Enterococcus faecalis (VRE) as a host. Nonaheksakonda is a lytic phage infecting E. faecalis V583 and clinical isolates with at least four different multi-locus sequence types (MLSTs). The genome is a 41.9-kb double-stranded DNA molecule (34.6% GC) with 74 coding sequences. Comparative analysis revealed only one close relative, Enterococcus phage heks. All other phages had low protein similarity and shared less than 54% nucleotide sequence identity with phage Nonaheksakonda. The most similar phages were all classified and unclassified efquatroviruses. We propose that the phages Nonaheksakonda and heks represent a novel genus within the family Siphoviridae, order Caudovirales, for which we propose the name "Nonaheksakondavirus".

Comparative genomics with a multidrug-resistant Klebsiella pneumoniae isolate reveals the panorama of unexplored diversity in Northeast Brazil

The emergence of community acquired infections increases the public health concern on K. pneumoniae and closely related bacteria among which antimicrobial resistance spreads. We report a multidrug-resistant K. pneumoniae isolate, B31, of a patient infected in the community and admitted to an intensive care unit in Northeast Brazil. Antimicrobial susceptibility and genome information were thoroughly investigated to characterize B31 in front of 172 sequenced strains of different countries. Assigned to the Sequence Type 15, which is globally spread, B31 presented extended spectrum beta-lactamase, tigecycline and ciprofloxacin resistance. Genome sequencing revealed most resistance genes being carried by plasmids with high dissemination potential. The absence of main virulence factors, like yersiniabactin and colibactin, apparently suggests a mild pathogenic strain which, on the contrary, persisted and caused severe infection in a previously healthy patient. The present study contributes to unveil the unclear genomic scenario of virulent and multidrug-resistant K. pneumoniae in Brazil.

Gut Microbiota and Metabolic Specificity in Ulcerative Colitis and Crohn's Disease

Background: Inflammatory bowel disease (IBD) represents multifactorial chronic inflammatory conditions in the gastrointestinal tract and includes Crohn's disease (CD) and ulcerative colitis (UC). Despite similarities in pathobiology and disease symptoms, UC and CD represent distinct diseases and exhibit diverse therapeutic responses. While studies have now confirmed that IBD is associated with dramatic changes in the gut microbiota, specific changes in the gut microbiome and associated metabolic effects on the host due to CD and UC are less well-understood. Methods: To address this knowledge gap, we performed an extensive unbiased meta-analysis of the gut microbiome data from five different IBD patient cohorts from five different countries using QIIME2, DIAMOND, and STAMP bioinformatics platforms. In-silico profiling of the metabolic pathways and community metabolic modeling were carried out to identify disease-specific association of the metabolic fluxes and signaling pathways. Results: Our results demonstrated a highly conserved gut microbiota community between healthy individuals and IBD patients at higher phylogenetic levels. However, at or below the order level in the taxonomic rank, we found significant disease-specific alterations. Similarly, we identified differential enrichment of the metabolic pathways in CD and UC, which included enriched pathways related to amino acid and glycan biosynthesis and metabolism, in addition to other metabolic pathways. Conclusions: In conclusion, this study highlights the prospects of harnessing the gut microbiota to improve understanding of the etiology of CD and UC and to develop novel prognostic, and therapeutic approaches.

MALDI-TOF MS and genomic analysis can make the difference in the clarification of canine brucellosis outbreaks

Brucellosis is one of the most common bacterial zoonoses worldwide affecting not only livestock and wildlife but also pets. Canine brucellosis is characterized by reproductive failure in dogs. Human Brucella canis infections are rarely reported but probably underestimated due to insufficient diagnostic surveillance. To improve diagnostics, we investigated dogs in a breeding kennel that showed clinical manifestations of brucellosis and revealed positive blood cultures. As an alternative to the time-consuming and hazardous classical identification procedures, a newly developed species-specific intact-cell matrix-assisted laser desorption/ionization–time of flight mass spectrometry analysis was applied, which allowed for rapid identification of B. canis and differentiation from closely related B. suis biovar 1. High-throughput sequencing and comparative genomics using single nucleotide polymorphism analysis clustered our isolates together with canine and human strains from various Central and South American countries in a distinct sub-lineage. Hence, molecular epidemiology clearly defined the outbreak cluster and demonstrated the endemic situation in South America. Our study illustrates that MALDI-TOF MS analysis using a validated in-house reference database facilitates rapid B. canis identification at species level. Additional whole genome sequencing provides more detailed outbreak information and leads to a deeper understanding of the epidemiology of canine brucellosis.

VIRIDIC-A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses

Nucleotide-based intergenomic similarities are useful to understand how viruses are related with each other and to classify them. Here we have developed VIRIDIC, which implements the traditional algorithm used by the International Committee on Taxonomy of Viruses (ICTV), Bacterial and Archaeal Viruses Subcommittee, to calculate virus intergenomic similarities. When compared with other software, VIRIDIC gave the best agreement with the traditional algorithm, which is based on the percent identity between two genomes determined by BLASTN. Furthermore, VIRIDIC proved best at estimating the relatedness between more distantly-related phages, relatedness that other tools can significantly overestimate. In addition to the intergenomic similarities, VIRIDIC also calculates three indicators of the alignment ability to capture the relatedness between viruses: the aligned fractions for each genome in a pair and the length ratio between the two genomes. The main output of VIRIDIC is a heatmap integrating the intergenomic similarity values with information regarding the genome lengths and the aligned genome fraction. Additionally, VIRIDIC can group viruses into clusters, based on user-defined intergenomic similarity thresholds. The sensitivity of VIRIDIC is given by the BLASTN. Thus, it is able to capture relationships between viruses having in common even short genomic regions, with as low as 65% similarity. Below this similarity level, protein-based analyses should be used, as they are the best suited to capture distant relationships. VIRIDIC is available at viridic.icbm.de, both as a web-service and a stand-alone tool. It allows fast analysis of large phage genome datasets, especially in the stand-alone version, which can be run on the user's own servers and can be integrated in bioinformatics pipelines. VIRIDIC was developed having viruses of Bacteria and Archaea in mind; however, it could potentially be used for eukaryotic viruses as well, as long as they are monopartite.