Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life

OrthoPhyl-streamlining large-scale, orthology-based phylogenomic studies of bacteria at broad evolutionary scales

There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis.

Genome-wide identification of MATE and ALMT genes and their expression profiling in mungbean (Vigna radiata L.) under aluminium stress

The Multidrug and toxic compound extrusion (MATE) and aluminium activated malate transporter (ALMT) gene families are involved in response to aluminium (Al) stress. In this study, we identified 48 MATE and 14 ALMT gene families in Vigna radiata genome and classified into 5 (MATE) and 3 (ALMT) clades by phylogenetic analysis. All the VrMATE and VrALMT genes were distributed across mungbean chromosomes. Tandem duplication was the main driving force for evolution and expansion of MATE gene family. Collinearity of mungbean with soybean indicated that MATE gene family is closely linked to Glycine max. Eight MATE transporters in clade 2 were found to be associated with previously characterized Al tolerance related MATEs in various plant species. Citrate exuding motif (CEM) was present in seven VrMATEs of clade 2. Promoter analysis revealed abundant plant hormone and stress responsive cis-elements. Results from quantitative real time-polymerase chain reaction (qRT-PCR) revealed that VrMATE19, VrMATE30 and VrALMT13 genes were markedly up-regulated at different time points under Al stress. Overall, this study offers a new direction for further molecular characterization of the MATE and ALMT genes in mungbean for Al tolerance.

First molecular characterization of Burkholderia mallei strains isolated from horses in Mongolia

Glanders, a highly contagious and often fatal disease affecting equids, is caused by Burkholderia mallei. Although sporadic cases of equine glanders have recently been documented in Mongolia, genome sequencing and molecular studies of the bacteria within this region are lacking. This study provided the first molecular characterization of B. mallei isolated from four native Mongolian horses from two different provinces in 2019 and 2022 by applying whole-genome sequencing with two SNP types (previously developed genotyping with 15 SNP markers that provide global coverage of the B. mallei population and the core genome coding SNP typing developed in this study). The Mongolian isolates were located within the L3B1 cluster, which was previously associated with the V-120 strain from Russia. Within the L3B1 cluster shared by neighboring countries, they were in a unique subbranch. In this study, specific SNP markers unique to the Mongolian strains were identified to track these strains using a high-resolution melting analysis (HRMA). This study revealed the unique phylogenetic background of Mongolian strains isolated from the eastern part of Mongolia. HRMA specific to the Mongolian subbranch may contribute to the molecular epidemiological monitoring of glanders in Mongolia and surrounding countries.

A novel variant in Salmonella genomic island 1 of multidrug-resistant Salmonella enterica serovar Kentucky ST198

Salmonella enterica serovar Kentucky ST198 is a major health threat due to its resistance to ciprofloxacin and several other drugs, including third-generation cephalosporins. Many drug-resistant genes have been identified in the Salmonella genomic island 1 variant K (SGI1-K). In this study, we investigated the antimicrobial resistance (AMR) profile and genotypic relatedness of two isolates of ciprofloxacin-resistant (CIPR) S. Kentucky ST198 from poultry in Northeastern Thailand. We successfully assembled the complete genomes of both isolates, namely SSSE-01 and SSSE-03, using hybrid de novo assembly of both short- and long-read sequence data. The complete genomes revealed their highly similar genomic structures and a novel variant of SGI1-K underlying multidrug-resistant (MDR) patterns, including the presence of blaTEM-1b, which confers resistance to beta-lactams, including cephalosporins and lnu(F) which confers resistance to lincomycin and other lincosamides. In addition, the chromosomal mutations in the quinolone resistance-determining region (QRDR) were found at positions 83 (Ser83Phe) and 87 (Asp87Asn) of GyrA and at positions 57 (Thr57Ser) and 80 (Ser80Ile) of ParC suggesting high resistance to ciprofloxacin. We also compared SSSE-01 and SSSE-03 with publicly available complete genome data and revealed significant variations in SGI1-K genetic structures and variable relationships to antibiotic resistance. In comparison to the other isolates, SGI1-K of SSSE-01 and SSSE-03 had a relatively large deletion in the backbone, spanning from S011 (traG∆) to S027 (resG), and the inversion of the IS26-S044∆-yidY segment. Their MDR region was characterized by the inversion of a large segment, including the mer operon and the relocation of IntI1 and several resistance genes downstream of the IS26-S044∆-yidY segment. These structural changes were likely mediated by the recombination of IS26. The findings broaden our understanding of the possible evolution pathway of SGI1-K in fostering drug resistance, which may provide opportunities to control these MDR strains.IMPORTANCEThe emergence of ciprofloxacin-resistant (CIPR) Salmonella Kentucky ST198 globally has raised significant concerns. This study focuses on two poultry isolates from Thailand, revealing a distinct Salmonella genomic island 1 variant K (SGI1-K) genetic structure. Remarkably, multiple antibiotic resistance genes (ARGs) were identified within the SGI1-K as well as other locations in the chromosome, but not in plasmids. Comparing the SGI1-K genetic structures among global and even within-country isolates unveiled substantial variations. Intriguingly, certain isolates lacked ARGs within the SGI1-K, while others had ARGs relocated outside. The presence of chromosomal extended-spectrum β-lactamase (ESBL) genes and lincosamide resistance, lnu(F), gene, could potentially inform the choices of the treatment of CIPRS. Kentucky ST198 infections in humans. This study highlights the importance of understanding the diverse genetic structures of SGI1-K and emphasizes the role of animals and humans in the emergence of antimicrobial resistance.

Genome sequence of the plant-growth-promoting bacterium Bacillus velezensis EU07

Many Gram-positive spore-forming rhizobacteria of the genus Bacillus show potential as biocontrol biopesticides that promise improved sustainability and ecological safety in agriculture. Here, we present a draft-quality genome sequence for Bacillus velezensis EU07, which shows growth-promotion in tomato plants and biocontrol against Fusarium head blight. We found that the genome of EU07 is almost identical to that of the commercially used strain QST713, but identified 46 single-nucleotide differences that distinguish these strains from each other. The availability of this genome sequence will facilitate future efforts to unravel the genetic and molecular basis for EU07's beneficial properties.

Genomic Characterization of Carbapenemase-Producing Enterobacter hormaechei, Serratia marcescens, Citrobacter freundii, Providencia stuartii, and Morganella morganii Clinical Isolates from Bulgaria

Carbapenemase-producing Enterobacter spp. Serratia marcescens, Citrobacter freundii, Providencia spp., and Morganella morganii (CP-ESCPM) are increasingly identified as causative agents of nosocomial infections but are still not under systematic genomic surveillance. In this study, using a combination of whole-genome sequencing and conjugation experiments, we sought to elucidate the genomic characteristics and transferability of resistance genes in clinical CP-ESCPM isolates from Bulgaria. Among the 36 sequenced isolates, NDM-1 (12/36), VIM-4 (11/36), VIM-86 (8/36), and OXA-48 (7/36) carbapenemases were identified; two isolates carried both NDM-1 and VIM-86. The majority of carbapenemase genes were found on self-conjugative plasmids. IncL plasmids were responsible for the spread of OXA-48 among E. hormaechei, C. freundii, and S. marcescens. IncM2 plasmids were generally associated with the spread of NDM-1 in C. freundii and S. marcescens, and also of VIM-4 in C. freundii. IncC plasmids were involved in the spread of the recently described VIM-86 in P. stuartii isolates. IncC plasmids carrying blaNDM-1 and blaVIM-86 were observed too. blaNDM-1 was also detected on IncX3 in S. marcescens and on IncT plasmid in M. morganii. The significant resistance transfer rates we observed highlight the role of the ESCPM group as a reservoir of resistance determinants and stress the need for strengthening infection control measures.

Isolation of Burkholderia pseudomallei from a goat in New Caledonia: implications for animal and human health monitoring and serological tool comparison

BACKGROUND

Melioidosis is a serious bacterial infection caused by Burkholderia pseudomallei, a gram-negative bacterium commonly found in soil and water. It can affect both humans and animals, and is endemic in regions such as Southeast Asia and Northern Australia. In recent years, there have been reports of an emergence of human melioidosis in other areas, including New Caledonia.

RESULTS

During standard laboratory analysis in New Caledonia in 2021, a strain of B. pseudomallei was isolated from a goat. The strain was characterized using both MLST and WGS techniques and was found to cluster with previously described local human strains from the area. In parallel, several serological tests (CFT, ELISA, Luminex (Hcp1, GroEL, BPSS1840), arrays assay and a latex agglutination test) were performed on animals from the farm where the goat originated, and/or from three other neighboring farms. Using two commercial ELISA kits, seropositive animals were found only on the farm where the infected goat originated and tests based on recombinant proteins confirmed the usefulness of the Hcp1 protein for the diagnosis of melioidosis in animals.

CONCLUSIONS

Despite the regular reports of human cases, this is the first confirmed case of melioidosis in an animal in New Caledonia. These results confirm the presence of the bacterium in the region and highlight the importance of vigilance for both animal and human health. It is critical that all health partners, including breeders, veterinarians, and biologists, work together to monitor and prevent the spread of the disease.

Update on the proposed minimal standards for the use of genome data for the taxonomy of prokaryotes

The field of microbial taxonomy is dynamic, aiming to provide a stable and contemporary classification system for prokaryotes. Traditionally, reliance on phenotypic characteristics limited the comprehensive understanding of microbial diversity and evolution. The introduction of molecular techniques, particularly DNA sequencing and genomics, has transformed our perception of prokaryotic diversity. In the past two decades, advancements in genome sequencing have transitioned from traditional methods to a genome-based taxonomic framework, not only to define species, but also higher taxonomic ranks. As technology and databases rapidly expand, maintaining updated standards is crucial. This work seeks to revise the 2018 guidelines for applying genome sequencing data in microbial taxonomy, adapting minimal standards and recommendations to reflect technological progress during this period.

VBCG: 20 validated bacterial core genes for phylogenomic analysis with high fidelity and resolution

BACKGROUND

Phylogenomic analysis has become an inseparable part of studies of bacterial diversity and evolution, and many different bacterial core genes have been collated and used for phylogenomic tree reconstruction. However, these genes have been selected based on their presence and single-copy ratio in all bacterial genomes, leaving out the gene's 'phylogenetic fidelity' unexamined.

RESULTS

From 30,522 complete genomes covering 11,262 species, we examined 148 bacterial core genes that have been previously used for phylogenomic analysis. In addition to the gene presence and single-copy rations, we evaluated the gene's phylogenetic fidelity by comparing each gene's phylogeny with its corresponding 16S rRNA gene tree. Out of the 148 bacterial genes, 20 validated bacterial core genes (VBCG) were selected as the core gene set with the highest bacterial phylogenetic fidelity. Compared to the larger gene set, the 20-gene core set resulted in more species having all genes present and fewer species with missing data, thereby enhancing the accuracy of phylogenomic analysis. Using Escherichia coli strains as examples of prominent bacterial foodborne pathogens, we demonstrated that the 20 VBCG produced phylogenies with higher fidelity and resolution at species and strain levels while 16S rRNA gene tree alone could not.

CONCLUSION

The 20 validated core gene set improves the fidelity and speed of phylogenomic analysis. Among other uses, this tool improves our ability to explore the evolution, typing and tracking of bacterial strains, such as human pathogens. We have developed a Python pipeline and a desktop graphic app (available on GitHub) for users to perform phylogenomic analysis with high fidelity and resolution. Video Abstract.

Building Phylogenetic Trees From Genome Sequences With kSNP4

Performing phylogenetic analysis with genome sequences maximizes the information used to estimate phylogenies and the resolution of closely related taxa. The use of single-nucleotide polymorphisms (SNPs) permits estimating trees without genome alignments and permits the use of data sets of hundreds of microbial genomes. kSNP4 is a program that identifies SNPs without using a reference genome, estimates parsimony, maximum likelihood, and neighbor-joining trees, and is able to annotate the discovered SNPs. kSNP4 is a command-line program that does not require any additional programs or dependencies to install or use. kSNP4 does not require any programming experience or bioinformatics experience to install and use. It is suitable for use by students through senior investigators. It includes a detailed user guide that explains all of the many features of kSNP4. In this study, we provide a detailed step-by-step protocol for downloading, installing, and using kSNP4 to build phylogenetic trees from genome sequences.

Arbuscular mycorrhizal fungi heterokaryons have two nuclear populations with distinct roles in host-plant interactions

Arbuscular mycorrhizal fungi (AMF) are prominent root symbionts that can carry thousands of nuclei deriving from two parental strains in a large syncytium. These co-existing genomes can also vary in abundance with changing environmental conditions. Here we assemble the nuclear genomes of all four publicly available AMF heterokaryons using PacBio high-fidelity and Hi-C sequencing. We find that the two co-existing genomes of these strains are phylogenetically related but differ in structure, content and epigenetics. We confirm that AMF heterokaryon genomes vary in relative abundance across conditions and show this can lead to nucleus-specific differences in expression during interactions with plants. Population analyses also reveal signatures of genetic exchange indicative of past events of sexual reproduction in these strains. This work uncovers the origin and contribution of two nuclear genomes in AMF heterokaryons and opens avenues for the improvement and environmental application of these strains.

Application of Next-Generation Sequencing to Enterobacter Hormaechei Subspecies Analysis during a Neonatal Intensive Care Unit Outbreak

INTRODUCTION

The Enterobacter cloacae complex (ECC) species are potential neonatal pathogens, and ECC strains are among the most commonly encountered Enterobacter spp. associated with nosocomial bloodstream infections. Outbreaks caused by ECC can lead to significant morbidity and mortality in susceptible neonates. At the molecular level, ECC exhibits genomic heterogeneity, with six closely related species and subspecies. Genetic variability poses a challenge in accurately identifying outbreaks by determining the clonality of ECC isolates. This difficulty is further compounded by the limitations of the commonly used molecular typing methods, such as pulsed field gel electrophoresis, which do not provide reliable accuracy in distinguishing between ECC strains and can lead to incorrect conclusions. Next-generation sequencing (NGS) offers superior resolution in determining strain relatedness. Therefore, we investigated the clinical pertinence of incorporating NGS into existing bundle measures to enhance patient management during an outbreak of ECC in a level-3 neonatal intensive care unit (NICU) in Germany.

METHODS

As the standard of care, all neonates on the NICU received weekly microbiological swabs (nasopharyngeal and rectal) and analysis of endotracheal secretion, where feasible. During the 2.5-month outbreak, colonisation with ECC was detected in n = 10 neonates. The phylogenetic relationship and potential antimicrobial resistance genes as well as mobile genetic elements were identified via bacterial whole-genome sequencing (WGS) using Illumina MiSeq followed by in silico data analysis.

RESULTS

Although all ECC isolates exhibited almost identical antimicrobial susceptibility patterns, the WGS data revealed the involvement of four different ECC clones. The isolates could be characterised as Enterobacter hormaechei subspecies steigerwaltii (n = 6, clonal), subsp. hoffmannii (n = 3, two clones) and subsp. oharae (n = 1). Despite the collection of environmental samples, no source of this diffuse outbreak could be identified. A new standardised operating procedure was implemented to enhance the management of neonates colonised with MRGN. This collaborative approach involved both parents and medical professionals and successfully prevented further transmission of ECC.

CONCLUSIONS

Initially, it was believed that the NICU outbreak was caused by a single ECC clone due to the similarity in antibiotic resistance. However, our findings show that antibiotic susceptibility patterns can be misleading in investigating outbreaks of multi-drug-resistant ECC. In contrast, bacterial WGS accurately identified ECC at the clonal level, which significantly helped to delineate the nature of the observed outbreak.

Spatial and phylogenetic patterns reveal hidden infection sources of Bacillus anthracis in an anthrax outbreak in Son La province, Vietnam

Bacillus anthracis, the bacterial cause of anthrax, is a zoonosis affecting livestock and wildlife often spilling over into humans. In Vietnam, anthrax has been nationally reportable since 2015 with cases occurring annually, mostly in the northern provinces. In April 2022, an outbreak was reported in Son La province following the butchering of a water buffalo, Bubalus bubalis. A total of 137 humans from three villages were likely exposed to contaminated meat from the animal. Early epidemiological investigations suggested a single animal was involved in all exposures. Five B. anthracis isolates were recovered from human clinical cases along with one from the buffalo hide, another from associated maggots, and one from soil at the carcass site. The isolates were whole genome sequenced, allowing global, regional, and local molecular epidemiological analyses of the outbreak strains. All recovered B. anthracis belong to the A.Br.001/002 lineage based on canonical single nucleotide polymorphism analysis (canSNP). Although not previously identified in Vietnam, this lineage has been identified in the nearby countries of China, India, Indonesia, Thailand, as well as Australia. A twenty-five marker multi-locus variable number tandem repeat analysis (MLVA-25) was used to investigate the relationship between human, soil, and buffalo strains. Locally, four MLVA-25 genotypes were identified from the eight isolates. This level of genetic diversity is unusual for the limited geography and timing of cases and differs from past literature using MLVA-25. The coupled spatial and phylogenetic data suggest this outbreak originated from multiple, likely undetected, animal sources. These findings were further supported by local news reports that identified at least two additional buffalo deaths beyond the initial animal sampled in response to the human cases. Future outbreak response should include intensive surveillance for additional animal cases and additional molecular epidemiological traceback to identify pathogen sources.

The core genome evolution of Lactobacillus crispatus as a driving force for niche competition in the human vaginal tract

The lower female reproductive tract is notoriously dominated by Lactobacillus species, among which Lactobacillus crispatus emerges for its protective and health-promoting activities. Although previous comparative genome analyses highlighted genetic and phenotypic diversity within the L. crispatus species, most studies have focused on the presence/absence of accessory genes. Here, we investigated the variation at the single nucleotide level within protein-encoding genes shared across a human-derived L. crispatus strain selection, which includes 200 currently available human-derived L. crispatus genomes as well as 41 chromosome sequences of such taxon that have been decoded in the framework of this study. Such data clearly pointed out the presence of intra-species micro-diversities that could have evolutionary significance contributing to phenotypical diversification by affecting protein domains. Specifically, two single nucleotide variations in the type II pullulanase gene sequence led to specific amino acid substitutions, possibly explaining the substantial differences in the growth performances and competition abilities observed in a multi-strain bioreactor culture simulating the vaginal environment. Accordingly, L. crispatus strains display different growth performances, suggesting that the colonisation and stable persistence in the female reproductive tract between the members of this taxon is highly variable.

Gene flow and species boundaries of the genus Salmonella

The genus Salmonella comprises two species, Salmonella bongori and Salmonella enterica, which are infectious to a wide variety of animal hosts. The diversity within S. enterica has been further partitioned into 6-10 subspecies based on such features as host range, geography, and most recently, genetic relatedness and phylogenetic affiliation. Although Salmonella pathogenicity is attributable to large numbers of acquired virulence factors, the extent of homologous exchange in the species at large is apparently constrained such that the species and subspecies form distinct clusters of strains. To explore the extent of gene flow within and among subspecies, and to ultimately define true biological species, we evaluated patterns of recombination in over 1,000 genomes currently assigned to the genus. Those Salmonella subspecies containing sufficient numbers of sequenced genomes to allow meaningful analysis-i.e., subsp. enterica and diarizonae-were found to be reproductively isolated from one another and from all other subspecies. Based on the configuration of genomic sequence divergence among subspecies, it is expected that each of the other Salmonella subspecies will also represent a biological species. Our findings argue against the application of prescribed nucleotide-identity thresholds to delineate bacterial species and contend that the Biological Species Concept should not be disregarded for bacteria, even those, like Salmonella, that demonstrate complex patterns of species and subspecies divergence. IMPORTANCE The Biological Species Concept (BSC), which defines species boundaries based on the capacity for gene exchange, is widely used to classify sexually reproducing eukaryotes but is generally thought to be inapplicable to bacteria due to their completely asexual mode of reproduction. We show that the genus Salmonella, whose thousands of described serovars were formerly considered to be strictly clonal, undergoes sufficient levels of homologous recombination to be assigned to species according to the BSC. Aside from the two recognized species, Salmonella enterica and Salmonella bongori, several (and likely all) of the subspecies within S. enterica are reproductively isolated from one another and should each be considered a separate biological species. These findings demonstrate that species barriers in bacteria can form despite high levels of nucleotide identity and that commonly applied thresholds of genomic sequence identity are not reliable indicators of bacterial species status.

Genomic and Phylogenetic Analysis of Bacillus cereus Biovar anthracis Isolated from Archival Bone Samples Reveals Earlier Natural History of the Pathogen

(1) Background: Bacillus cereus biovar anthracis (Bcbva) was the causative agent of an anthrax-like fatal disease among wild chimpanzees in 2001 in Côte d'Ivoire. Before this, there had not been any description of an anthrax-like disease caused by typically avirulent Bacillus cereus. Genetic analysis found that B. cereus had acquired two anthrax-like plasmids, one a pXO1-like toxin producing plasmid and the other a pXO2-like plasmid encoding capsule. Bcbva caused animal fatalities in Cameroon, Democratic Republic of Congo, and the Central African Republic between 2004 and 2012. (2) Methods: The pathogen had acquired plasmids in the wild and that was discovered as the cause of widespread animal fatalities in the early 2000s. Primate bones had been shipped out of the endemic zone for anthropological studies prior to the realized danger of contamination with Bcbva. Spores were isolated from the bone fragments and positively identified as Bcbva. Strains were characterized by classical microbiological methods and qPCR. Four new Bcbva isolates were whole-genome sequenced. Chromosomal and plasmid phylogenomic analysis was performed to provide temporal and spatial context to these new strains and previously sequenced Bcbva. Tau and principal component analyses were utilized to identify genetic and spatial case patterns in the Taï National Park anthrax zone. (3) Results: Preliminary studies positively identified Bcbva presence in several archival bone fragments. The animals in question died between 1994 and 2010. Previously, the earliest archival strains of Bcbva were identified in 1996. Though the pathogen has a homogeneous genome, spatial analyses of a subset of mappable isolates from Taï National Park revealed strains found closer together were generally more similar, with strains from chimpanzees and duikers having the widest distribution. Ancestral strains were located mostly in the west of the park and had lower spatial clustering compared to more recent isolates, indicating a local increase in genetic diversity of Bcbva in the park over space and time. Global clustering analysis indicates patterns of genetic diversity and distance are shared between the ancestral and more recently isolated type strains. (4) Conclusions: Our strains have the potential to unveil historical genomic information not available elsewhere. This information sheds light on the evolution and emergence of a dangerous anthrax-causing pathogen.

Genome sequence data of the strongly antagonistic yeast Pichia kluyveri isolate APC 11.10 B as a foundation for analysing biocontrol mechanisms

Pichia kluyveri strain APC 11.10 B was isolated from apple bark in Switzerland and exhibited strong antagonistic activity against plant pathogenic fungi in vitro (e.g., Botrytis, Fusarium or Monilinia isolates). In order to identify the mechanisms underlying this antagonism, we have sequenced the genome of this isolate by long- and short-read sequencing technologies. The sequence data were de novo assembled into nine scaffolds and a fully resolved circularized mitogenome. The total genome size was 10.9 Mbp and 7451 potential open reading frames (ORFs) and 202 tRNA genes were predicted. In comparison to two P. kluyveri genomes deposited at the NCBI (of strains X31-10 and CBA6002), the APC 11.10 B strain seemed to represent a hybrid because backmapping of sequencing reads resulted in a high rate of heterozygous and structural variants in the nuclear genome (this was not observed for the mitochondrial genome). The P. kluyveri (APC 11.10 B) draft genome represents a first step and resource for genome mining, comparative and functional genomics (e.g., identifying the biocontrol mode of action), and evolutionary studies. Since the genus Pichia comprises many biotechnologically relevant yeasts, the genome data may be used in a variety of fields and disciplines.

Draft genome sequences for ten strains of Xanthomonas species that have phylogenomic importance

Here we report draft-quality genome sequences for pathotype strains of eight plant-pathogenic bacterial pathovars: Xanthomonas campestris pv. asclepiadis, X. campestris pv. cannae, X. campestris pv. esculenti, X. campestris pv. nigromaculans, X. campestris pv. parthenii, X. campestris pv. phormiicola, X. campestris pv. zinniae and X. dyei pv. eucalypti (= X. campestris pv. eucalypti). We also sequenced the type strain of species X. melonis and the unclassified Xanthomonas strain NCPPB 1067. These data will be useful for phylogenomic and taxonomic studies, filling some important gaps in sequence coverage of Xanthomonas phylogenetic diversity. We include representatives of previously under-sequenced pathovars and species-level clades. Furthermore, these genome sequences may be useful in elucidating the molecular basis for important phenotypes, such as biosynthesis of coronatine-related toxins and degradation of fungal toxin cercosporin.

Salmonella gallinarum strains from outbreaks of fowl typhoid fever in Southern Africa closely related to SG9R vaccines

Introduction

Salmonella enterica subspecies enterica serovar Gallinarum biovar Gallinarum (SG) is associated with fowl typhoid fever, and the attenuated rough strain SG9R is widely used as a vaccine in many regions. Reversion to virulence of vaccine strains was suspected as the cause during recent fowl typhoid fever outbreaks in poultry in South Africa and Eswatini.

Methods

To compare nine field isolates with global wild-type SG9 strains and the two commercial SG9R vaccines in use, Nobilis^® SG9R and Cevac^®-SG, we used whole-genome comparison with single-nucleotide polymorphism (SNP) detection.

Results

SNP phylogenic analysis showed that all the southern African field isolates were more closely related to the vaccine strains than wild-type SG9 strains. Furthermore, SNPs in the pyruvate dehydrogenase (aceE) and/or lipopolysaccharide 1,2-glucosyltransferase (rfaJ) genes, which are known markers of attenuation, were found in four of the field isolates along with intact spv, SPI-1, and SPI-2 gene clusters, providing conclusive evidence that these four isolates were originally vaccine strains that reverted to virulence. Five other field isolates lacked the SG9R attenuation markers, but variant analysis identified an SNP in the yihX gene, insertions in the ybjX and hydH genes, and deletions in the ftsK and sadA genes that were shared between the field isolates and vaccine strains but absent in wild-type SG9, indicating that these field isolates were also likely revertant vaccines.

Discussion

Overall, this study highlights different mechanisms of reversion of two commercial vaccines, where virulence caused by field isolates closely related to the Nobilis^® SG9R vaccine was associated with the restoration of intact virulence gene clusters, and those derived from the Cevac^®-SG vaccine were characterized by point mutations resulting in restored aceE and rfaJ genes. A possible new marker of attenuation was identified as a point mutation in the yihX gene, as well as four new candidate genes that could potentially be used to distinguish current vaccine strains from wild-type strains using PCR assays.

Phylogenomic analysis supports Mycobacterium tuberculosis transmission between humans and elephants

Introduction

Tuberculosis is an infectious disease caused by a group of acid-fast bacilli known as Mycobacterium tuberculosis complex (MTC), which has a major impact on humans. Transmission of MTC across the human-animal interface has been demonstrated by several studies. However, the reverse zoonotic transmission from humans to animals (zooanthroponosis) has often been neglected.

Methods

In this study, we used Nanopore MinION and Illumina MiSeq approaches to sequence the whole genome of M. tuberculosis strains isolated from two deceased Asian elephants (Elephas maximus) and one human in Chitwan, Nepal. The evolutionary relationships and drug resistance capacity of these strains were assessed using the whole genome data generated by the stand-alone tool Tb-Profiler. Phylogenomic trees were also constructed using a non-synonymous SNP alignment of 2,596 bp, including 94 whole genome sequences representative of the previously described M. tuberculosis lineages from elephants worldwide (lineages 1 and 4) and from humans in Nepal (lineages 1, 2 and 3).

Results and Discussion

The new genomes achieved an average coverage of 99.6%, with an average depth of 55.67x. These M. tuberculosis strains belong to lineage 1 (elephant DG), lineage 2 (elephant PK) and lineage 4 (human), and none of them were found to have drug-resistant variants. The elephant-derived isolates were evolutionarily closely related to human-derived isolates previously described in Nepal, both in lineages 1 and 2, providing additional support for zooanthroponosis or bidirectional transmission between humans and elephants. The human-derived isolate clustered together with other published human isolates from Argentina, Russia and the United Kingdom in the lineage 4 clade. This complex multi-pathogen, multi-host system is challenging and highlights the need for a One Health approach to tuberculosis prevention and control at human-animal interface, particularly in regions where human tuberculosis is highly endemic.

Assembly and comparative genome analysis of a Patagonian Aureobasidium pullulans isolate reveals unexpected intraspecific variation

Aureobasidium pullulans is a yeast-like fungus with remarkable phenotypic plasticity widely studied for its importance for the pharmaceutical and food industries. So far, genomic studies with strains from all over the world suggest they constitute a genetically unstructured population, with no association by habitat. However, the mechanisms by which this genome supports so many phenotypic permutations are still poorly understood. Recent works have shown the importance of sequencing yeast genomes from extreme environments to increase the repertoire of phenotypic diversity of unconventional yeasts. In this study, we present the genomic draft of A. pullulans strain from a Patagonian yeast diversity hotspot, re-evaluate its taxonomic classification based on taxogenomic approaches, and annotate its genome with high-depth transcriptomic data. Our analysis suggests this isolate could be considered a novel variant at an early stage of the speciation process. The discovery of divergent strains in a genomically homogeneous group, such as A. pullulans, can be valuable in understanding the evolution of the species. The identification and characterization of new variants will not only allow finding unique traits of biotechnological importance, but also optimize the choice of strains whose phenotypes will be characterized, providing new elements to explore questions about plasticity and adaptation.

Beyond the spore, the exosporium sugar anthrose impacts vegetative Bacillus anthracis gene regulation in cis and trans

The Bacillus anthracis exosporium nap is the outermost portion of spore that interacts with the environment and host systems. Changes to this layer have the potential to impact wide-ranging physiological and immunological processes. The unique sugar, anthrose, normally coats the exosporium nap at its most distal points. We previously identified additional mechanisms rendering B. anthracis anthrose negative. In this work, several new ant ^- B. anthracis strains are identified and the impact of anthrose negativity on spore physiology is investigated. We demonstrate that live-attenuated Sterne vaccines as well as culture filtrate anthrax vaccines generate antibodies targeting non-protein components of the spore. The role of anthrose as a vegetative B. anthracis Sterne signaling molecule is implicated by luminescent expression strain assays, RNA-seq experiments, and toxin secretion analysis by western blot. Pure anthrose and the sporulation-inducing nucleoside analogue decoyinine had similar effects on toxin expression. Co-culture experiments demonstrated gene expression changes in B. anthracis depend on intracellular anthrose status (cis) in addition to anthrose status of extracellular interactions (trans). These findings provide a mechanism for how a unique spore-specific sugar residue affects physiology, expression and genetics of vegetative B. anthracis with impacts on the ecology, pathogenesis, and vaccinology of anthrax.

Clinical Aureobasidium Isolates Are More Fungicide Sensitive than Many Agricultural Isolates

Fungicide applications in agriculture and medicine can promote the evolution of resistant, pathogenic fungi, which is a growing problem for disease management in both settings. Nonpathogenic mycobiota are also exposed to fungicides, may become tolerant, and could turn into agricultural or medical problems, for example, due to climate change or in immunocompromised individuals. However, quantitative data about fungicide sensitivity of environmental fungi is mostly lacking. Aureobasidium species are widely distributed and frequently isolated yeast-like fungi. One species, A. pullulans, is used as a biocontrol agent, but is also encountered in clinical samples, regularly. Here, we compared 16 clinical and 30 agricultural Aureobasidium isolates based on whole-genome data and by sensitivity testing with the 3 fungicides captan, cyprodinil, and difenoconazole. Our phylogenetic analyses determined that 7 of the 16 clinical isolates did not belong to the species A. pullulans. These isolates clustered with other Aureobasidium species, including A. melanogenum, a recently separated species that expresses virulence traits that are mostly lacking in A. pullulans. Interestingly, the clinical Aureobasidium isolates were significantly more fungicide sensitive than many isolates from agricultural samples, which implies selection for fungicide tolerance of non-target fungi in agricultural ecosystems. IMPORTANCE Environmental microbiota are regularly found in clinical samples and can cause disease, in particular, in immunocompromised individuals. Organisms of the genus Aureobasidium belonging to this group are highly abundant, and some species are even described as pathogens. Many A. pullulans isolates from agricultural samples are tolerant to different fungicides, and it seems inevitable that such strains will eventually appear in the clinics. Selection for fungicide tolerance would be particularly worrisome for species A. melanogenum, which is also found in the environment and exhibits virulence traits. Based on our observation and the strains tested here, clinical Aureobasidium isolates are still fungicide sensitive. We, therefore, suggest monitoring fungicide sensitivity in species, such as A. pullulans and A. melanogenum, and to consider the development of fungicide tolerance in the evaluation process of fungicides.

Phylogenomic Analysis Supports the Transfer of 20 Pathovars from Xanthomonas campestris into Xanthomonas euvesicatoria

The Gram-negative bacterial genus Xanthomonas includes numerous infra-specific taxa known as pathovars, which are defined primarily on host range and disease symptoms. With the advent of molecular sequence data, many pathovars have been transferred from X. campestris into other Xanthomonas species to better harmonise taxonomy and phylogeny. We performed whole-genome shotgun sequencing on pathotype strains of the following X. campestris pathovars: blepharidis, carissae, clerodendri, convolvuli, coriandri, daturae, euphorbiae, fici, heliotropii, ionidii, lawsoniae, mirabilis, obscurae, paulliniae, pennamericanum, spermacoces, uppalii, vernoniae, viegasii and zingibericola. These genomes showed more than 98% average nucleotide identity with the type-strain of X. euvesicatoria and less than 88% with the type-strain of X. campestris. We propose the transfer of these pathovars into X. euvesicatoria and present an emended species description for X. euvesicatoria.

Escherichia Coli: What Is and Which Are?

Escherichia coli have served as important model organisms for over a century-used to elucidate key aspects of genetics, evolution, molecular biology, and pathogenesis. However, defining which strains actually belong to this species is erratic and unstable due to shifts in the characters and criteria used to distinguish bacterial species. Additionally, many isolates designated as E. coli are genetically more closely related to strains of Shigella than to other E. coli, creating a situation in which the entire genus of Shigella and its four species are encompassed within the single species E. coli. We evaluated all complete genomes assigned to E. coli and its closest relatives according to the biological species concept (BSC), using evidence of reproductive isolation and gene flow (i.e., homologous recombination in the case of asexual bacteria) to ascertain species boundaries. The BSC establishes a uniform, consistent, and objective principle that allows species-level classification across all domains of life and does not rely on either phenotypic or genotypic similarity to a defined type-specimen for species membership. Analyzing a total of 1,887 sequenced genomes and comparing our results to other genome-based classification methods, we found few barriers to gene flow among the strains, clades, phylogroups, or species within E. coli and Shigella. Due to the utility in recognizing which strains constitute a true biological species, we designate genomes that form a genetic cohesive group as members of E. coliBIO.

Large-scale comparative genomics to refine the organization of the global Salmonella enterica population structure

The White-Kauffmann-Le Minor (WKL) scheme is the most widely used Salmonella typing scheme for reporting the disease prevalence of the enteric pathogen. With the advent of whole-genome sequencing (WGS), in silico methods have increasingly replaced traditional serotyping due to reproducibility, speed and coverage. However, despite integrating genomic-based typing by in silico serotyping tools such as SISTR, in silico serotyping in certain contexts remains ambiguous and insufficiently informative. Specifically, in silico serotyping does not attempt to resolve polyphyly. Furthermore, in spite of the widespread acknowledgement of polyphyly from genomic studies, the prevalence of polyphyletic serovars is not well characterized. Here, we applied a genomics approach to acquire the necessary resolution to classify genetically discordant serovars and propose an alternative typing scheme that consistently reflect natural Salmonella populations. By accessing the unprecedented volume of bacterial genomic data publicly available in GenomeTrakr and PubMLST databases (>180 000 genomes representing 723 serovars), we characterized the global Salmonella population structure and systematically identified putative non-monophyletic serovars. The proportion of putative non-monophyletic serovars was estimated higher than previous reports, reinforcing the inability of antigenic determinants to depict the complexity of Salmonella evolutionary history. We explored the extent of genetic diversity masked by serotyping labels and found significant intra-serovar molecular differences across many clinically important serovars. To avoid false discovery due to incorrect in silico serotyping calls, we cross-referenced reported serovar labels and concluded a low error rate in in silico serotyping. The combined application of clustering statistics and genome-wide association methods demonstrated effective characterization of stable bacterial populations and explained functional differences. The collective methods adopted in our study have practical values in establishing genomic-based typing nomenclatures for an entire microbial species or closely related subpopulations. Ultimately, we foresee an improved typing scheme to be a hybrid that integrates both genomic and antigenic information such that the resolution from WGS is leveraged to improve the precision of subpopulation classification while preserving the common names defined by the WKL scheme.

Increased zinc levels facilitate phenotypic detection of ceftazidime-avibactam resistance in metallo-β-lactamase-producing Gram-negative bacteria

Ceftazidime-avibactam is one of the last resort antimicrobial agents for the treatment of carbapenem-resistant, Gram-negative bacteria. Metallo-β-lactamase-producing bacteria are considered to be ceftazidime-avibactam resistant. Here, we evaluated a semi-automated antimicrobial susceptibility testing system regarding its capability to detect phenotypic ceftazidime-avibactam resistance in 176 carbapenem-resistant, metallo-β-lactamase-producing Enterobacterales and Pseudomonas aeruginosa isolates. Nine clinical isolates displayed ceftazidime-avibactam susceptibility in the semi-automated system and six of these isolates were susceptible by broth microdilution, too. In all nine isolates, metallo-β-lactamase-mediated hydrolytic activity was demonstrated with the EDTA-modified carbapenemase inactivation method. As zinc is known to be an important co-factor for metallo-β-lactamase activity, test media of the semi-automated antimicrobial susceptibility testing system and broth microdilution were supplemented with zinc. Thereby, the detection of phenotypic resistance was improved in the semi-automated system and in broth microdilution. Currently, ceftazidime-avibactam is not approved as treatment option for infections by metallo-β-lactamase-producing, Gram-negative bacteria. In infections caused by carbapenem-resistant Gram-negatives, we therefore recommend to rule out the presence of metallo-β-lactamases with additional methods before initiating ceftazidime-avibactam treatment.

Metagenome-Assembled Genomes of Komagataeibacter from Kombucha Exposed to Mars-Like Conditions Reveal the Secrets in Tolerating Extraterrestrial Stresses

Kombucha mutualistic community (KMC) is composed by acetic acid bacteria and yeasts, producing fermented tea with health benefits. As part of the BIOlogy and Mars EXperiment (BIOMEX) project, the effect of Mars-like conditions on the KMC was analyzed. Here, we analyzed metagenome-assembled genomes (MAGs) of the Komagataeibacter, which is a predominant genus in KMC, to understand their roles in the KMC after exposure to Mars-like conditions (outside the International Space Station) based on functional genetic elements. We constructed three MAGs: K. hansenii, K. rhaeticus, and K. oboediens. Our results showed that (i) K. oboediens MAG functionally more complex than K. hansenii, (ii) K. hansenii is a keystone in KMCs with specific functional features to tolerate extreme stress, and (iii) genes related to the PPDK, betaine biosynthesis, polyamines biosynthesis, sulfate-sulfur assimilation pathway as well as type II toxin-antitoxin (TA) system, quorum sensing (QS) system, and cellulose production could play important roles in the resilience of KMC after exposure to Mars-like stress. Our findings show the potential mechanisms through which Komagataeibacter tolerates the extraterrestrial stress and will help to understand minimal microbial composition of KMC for space travelers.

Whole-genome single nucleotide polymorphism analysis for typing the pandemic pathogen Fusarium graminearum sensu stricto

Recent improvements in microbiology and molecular epidemiology were largely stimulated by whole- genome sequencing (WGS), which provides an unprecedented resolution in discriminating highly related genetic backgrounds. WGS is becoming the method of choice in epidemiology of fungal diseases, but its application is still in a pioneer stage, mainly due to the limited number of available genomes. Fungal pathogens often belong to complexes composed of numerous cryptic species. Detecting cryptic diversity is fundamental to understand the dynamics and the evolutionary relationships underlying disease outbreaks. In this study, we explore the value of whole-genome SNP analyses in identification of the pandemic pathogen Fusarium graminearum sensu stricto (F.g.). This species is responsible for cereal diseases and negatively impacts grain production worldwide. The fungus belongs to the monophyletic fungal complex referred to as F. graminearum species complex including at least sixteen cryptic species, a few among them may be involved in cereal diseases in certain agricultural areas. We analyzed WGS data from a collection of 99 F.g. strains and 33 strains representing all known cryptic species belonging to the FGSC complex. As a first step, we performed a phylogenomic analysis to reveal species-specific clustering. A RAxML maximum likelihood tree grouped all analyzed strains of F.g. into a single clade, supporting the clustering-based identification approach. Although, phylogenetic reconstructions are essential in detecting cryptic species, a phylogenomic tree does not fulfill the criteria for rapid and cost-effective approach for identification of fungi, due to the time-consuming nature of the analysis. As an alternative, analysis of WGS information by mapping sequence data from individual strains against reference genomes may provide useful markers for the rapid identification of fungi. We provide a robust framework for typing F.g. through the web-based PhaME workflow available at EDGE bioinformatics. The method was validated through multiple comparisons of assembly genomes to F.g. reference strain PH-1. We showed that the difference between intra- and interspecies variability was at least two times higher than intraspecific variation facilitating successful typing of F.g. This is the first study which employs WGS data for typing plant pathogenic fusaria.

The emergence of a multidrug resistant Salmonella Muenchen in Israel is associated with horizontal acquisition of the epidemic pESI plasmid

OBJECTIVES

Horizontal acquisition of mobile genetic elements is a powerful evolutionary driving force that can profoundly affect pathogens epidemiology and their interactions with the environment and host. In the last decade, the role of the epidemic megaplasmid, pESI was demonstrated in the global emergence of multi-drug resistant (MDR) Salmonella enterica serovar Infantis strains, but it was unknown if this was a one-time phenomenon, or that pESI can drive the emergence of other pathogens.

METHODS

Epidemiological, molecular, whole genome sequencing, de-novo assembly, bioinformatics and genetic approaches were used to analyze the emergence of a pESI-positive Salmonella enterica serovar Muenchen strain in Israel.

RESULTS

Since 2018, we report the emergence and high prevalence of S. Muenchen in Israel, which consisted at 2020, 40% (1055/2671) of all clinical Salmonella isolates. We show that the emergence of S. Muenchen is dominated by a clonal MDR strain, report its complete assembled genome sequence, and demonstrate that in contrast to preemergent strains, it harbors the epidemic megaplasmid, pESI, which can be self-mobilized into E. coli and other Salmonella serovars. Additionally, we identified bioinformatically highly similar genomes of clinical isolates that were recently collected in South Africa, UK and USA.

CONCLUSIONS

This is a second documented case of a pathogen emergence associated with pESI acquisition. Considering the genetic cargo of pESI that enhances resistance, stress tolerance and virulence, and its ability to conjugate into prevalent Salmonella serovars, we provide further support that pESI facilities the emergence and spreading of new Salmonella strains.

Genomic signatures of high-altitude adaptation and chromosomal polymorphism in geladas

Primates have adapted to numerous environments and lifestyles, but very few species are native to high elevations. Here, we investigated high-altitude adaptations in the gelada (Theropithecus gelada), a monkey endemic to the Ethiopian Plateau. We examined genome-wide variation in conjunction with measurements of hematological and morphological traits. Our new gelada reference genome is highly intact and assembled at chromosome-length levels. Unexpectedly, we identified a chromosomal polymorphism in geladas that could potentially contribute to reproductive barriers between populations. Compared to baboons at low altitude, we found that high-altitude geladas exhibit significantly expanded chest circumferences, potentially allowing for greater lung surface area for increased oxygen diffusion. We identified gelada-specific amino acid substitutions in the alpha-chain subunit of adult hemoglobin but found that gelada hemoglobin does not exhibit markedly altered oxygenation properties compared to lowland primates. We also found that geladas at high altitude do not exhibit elevated blood hemoglobin concentrations, in contrast to the normal acclimatization response to hypoxia in lowland primates. The absence of altitude-related polycythemia suggests that geladas are able to sustain adequate tissue-oxygen delivery despite environmental hypoxia. Finally, we identified numerous genes and genomic regions exhibiting accelerated rates of evolution, as well as gene families exhibiting expansions in the gelada lineage, potentially reflecting altitude-related selection. Our findings lend insight into putative mechanisms of high-altitude adaptation while suggesting promising avenues for functional hypoxia research.

Integration of Molecular Docking Analysis and Molecular Dynamics Simulations for Studying Food Proteins and Bioactive Peptides

In silico tools, such as molecular docking, are widely applied to study interactions and binding affinity of biological activity of proteins and peptides. However, restricted sampling of both ligand and receptor conformations and use of approximated scoring functions can produce results that do not correlate with actual experimental binding affinities. Molecular dynamics simulations (MDS) can provide valuable information in deciphering functional mechanisms of proteins/peptides and other biomolecules, overcoming the rigid sampling limitations in docking analysis. This review will discuss the information related to the traditional use of in silico models, such as molecular docking, and its application for studying food proteins and bioactive peptides, followed by an in-depth introduction to the theory of MDS and description of why these molecular simulation techniques are important in the theoretical prediction of structural and functional dynamics of food proteins and bioactive peptides. Applications, limitations, and future prospects of MDS will also be discussed.

Development of Single Nucleotide Polymorphism (SNP)-Based Triplex PCR Marker for Serotype-Specific Escherichia coli Detection

Single-nucleotide polymorphisms (SNPs) are one of the most common forms of genetic variation and as such are powerful tools for the identification of bacterial strains, their genetic diversity, phylogenetic analysis, and outbreak surveillance. In this study, we used 15 sets of SNP-containing primers to amplify and sequence the target Escherichia coli. Based on the combination of the 15-sequence primer sets, each SNP site encompassing forward and reverse primer sequences (620–919 bp) were aligned and an SNP-based marker was designed. Each SNP marker exists in at least two SNP sites at the 3′ end of each primer; one natural and the other artificially created by transition or transversion mutation. Thus, 12 sets of SNP primers (225–488 bp) were developed for validation by amplifying the target E. coli. Finally, a temperature gradient triplex PCR kit was designed to detect target E. coli strains. The selected primers were amplified in three genes (ileS, thrB, and polB), with fragment sizes of 401, 337, and 232 bp for E. coli O157:H7, E. coli, and E. coli O145:H28, respectively. This allele-specific SNP-based triplex primer assay provides serotype-specific detection of E. coli strains in one reaction tube. The developed marker would be used to diagnose, investigate, and control food-borne E. coli outbreaks.

Morphological and Molecular Characterization of Paragonimus Species Isolated from Freshwater Crabs in Southern Yunnan, China

Paragonimus species are highly prevalent in various regions of China. The study's objective is to isolate and identify Paragonimus from natural habitats and compare the phylogenetic diversity of Paragonimus in southern Yunnan province, China. Metacercariae of Paragonimus was isolated from crabs, and morphologic identification was performed by microscopy. Metacercariae were injected into experimental Paragonimus free Sprague Dawley rats. After 114 days, adult worms and eggs were isolated from multiple organs. Morphologic identification confirmed the initial identification. DNA was extracted from 5 adult worms, and molecular characterization was performed by amplification and sequencing of CO1 and ITS2 regions, followed by phylogenetic analysis. Out of 447 crabs captured, 186 crabs were found to be infected. A total of 4 species of Paragonimus was observed from naturally infected crabs. Paragonimus microrchis (2), Paragonimus heterotremus (1), Paragonimus proliferus (1), and Paragonimus skrjabini (1) were isolated and identified. A total of 32 sequences were downloaded from the National Center for Biotechnology Information, and 5 sequences generated in the study were used for phylogenetic analysis. In the phylogenetic tree of the CO1 gene, Paragonimus proliferus, Paragonimus heterotremus, and Paragonimus skrjabini were clustered with the same species, and the confidence values of their branches were >95%. A congruent phylogenetic relationship was observed with the ITS2 phylogenetic tree. In the phylogenetic tree constructed with the combined dataset of CO1 and ITS2 datasets, Paragonimus proliferus, Paragonimus heterotremus, and Paragonimus skrjabini clustered with the same species, and their branch confidence values were >94%. Paragonimus microrchis clustered with Paragonimus bangkokensis in both datasets. Phylogenetic analysis revealed robustness of the double loci method as against the single-locus method with either CO1 or ITS2 alone. Paragonimus species isolated from the southern Yunnan province, China, was phylogenetically diverse, and the analysis revealed the clustering of multiple species of Paragonimus isolated from different geographic locations.

Antimicrobial Susceptibility of Lactobacillus delbrueckii subsp. lactis from Milk Products and Other Habitats

As components of many cheese starter cultures, strains of Lactobacillus delbrueckii subsp. lactis (LDL) must be tested for their antimicrobial susceptibility to avoid the potential horizontal transfer of antibiotic resistance (ABR) determinants in the human body or in the environment. To this end, a phenotypic test, as well as a screening for antibiotic resistance genes (ARGs) in genome sequences, is commonly performed. Historically, microbiological cutoffs (MCs), which are used to classify strains as either 'sensitive' or 'resistant' based on the minimal inhibitory concentrations (MICs) of a range of clinically-relevant antibiotics, have been defined for the whole group of the obligate homofermentative lactobacilli, which includes LDL among many other species. This often leads to inaccuracies in the appreciation of the ABR status of tested LDL strains and to false positive results. To define more accurate MCs for LDL, we analyzed the MIC profiles of strains originating from various habitats by using the broth microdilution method. These strains' genomes were sequenced and used to complement our analysis involving a search for ARGs, as well as to assess the phylogenetic proximity between strains. Of LDL strains, 52.1% displayed MICs that were higher than the defined MCs for kanamycin, 9.9% for chloramphenicol, and 5.6% for tetracycline, but no ARG was conclusively detected. On the other hand, all strains displayed MICs below the defined MCs for ampicillin, gentamycin, erythromycin, and clindamycin. Considering our results, we propose the adaptation of the MCs for six of the tested clinically-relevant antibiotics to improve the accuracy of phenotypic antibiotic testing.

Cassava diseases caused by Xanthomonas phaseoli pv. manihotis and Xanthomonas cassavae

Xanthomonas phaseoli pv. manihotis (Xpm) and X. cassavae (Xc) are two bacterial pathogens attacking cassava. Cassava bacterial blight (CBB) is a systemic disease caused by Xpm, which might have dramatic effects on plant growth and crop production. Cassava bacterial necrosis is a nonvascular disease caused by Xc with foliar symptoms similar to CBB, but its impacts on the plant vigour and the crop are limited. In this review, we describe the epidemiology and ecology of the two pathogens, the impacts and management of the diseases, and the main research achievements for each pathosystem. Because Xc data are sparse, our main focus is on Xpm and CBB.

Widespread bacterial diversity within the bacteriome of fungi

Knowledge of associations between fungal hosts and their bacterial associates has steadily grown in recent years as the number and diversity of examinations have increased, but current knowledge is predominantly limited to a small number of fungal taxa and bacterial partners. Here, we screened for potential bacterial associates in over 700 phylogenetically diverse fungal isolates, representing 366 genera, or a tenfold increase compared with previously examined fungal genera, including isolates from several previously unexplored phyla. Both a 16 S rDNA-based exploration of fungal isolates from four distinct culture collections spanning North America, South America and Europe, and a bioinformatic screen for bacterial-specific sequences within fungal genome sequencing projects, revealed that a surprisingly diverse array of bacterial associates are frequently found in otherwise axenic fungal cultures. We demonstrate that bacterial associations with diverse fungal hosts appear to be the rule, rather than the exception, and deserve increased consideration in microbiome studies and in examinations of microbial interactions.

Gulosibacter hominis sp. nov.: a novel human microbiome bacterium that may cause opportunistic infections

We present genomic, phylogenomic, and phenotypic taxonomic data to demonstrate that three human ear isolates represent a novel species within the genus Gulosibacter. These isolates could not be identified reliably using MALDI-TOF mass spectrometry during routine diagnostic work, but partial 16S rRNA gene sequence analysis revealed that they belonged to the genus Gulosibacter. Overall genomic relatedness indices between the draft genome sequences of the three isolates and of the type strains of established Gulosibacter species confirmed that the three isolates represented a single novel Gulosibacter species. A biochemical characterisation yielded differential tests between the novel and established Gulosibacter species, which could also be differentiated using MALDI-TOF mass spectrometry. We propose to formally classify these three isolates into Gulosibacter hominis sp. nov., with 401352-2018 ^T (= LMG 31778 ^T, CCUG 74795 ^T) as the type strain. The whole-genome sequence of strain 401352-2018 ^T has a size of 2,340,181 bp and a G+C content of 62.04 mol%. A Gulosibacter pangenome analysis revealed 467 gene clusters that were exclusively present in G. hominis genomes. While these G. hominis specific gene clusters were enriched in several COG functional categories, this analysis did not reveal functions that suggested a role in the human microbiome, nor did it explain the occurrence of G. hominis in ear infections. The absence of acquired antimicrobial resistance determinants and virulence factors in the G. hominis genomes, and an analysis of publicly available 16S rRNA gene sequences and 16S rRNA amplicon sequencing data sets suggested that G. hominis is a member of the human skin microbiota that may occasionally be involved in opportunistic infections.

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.

Genomic Analysis of SARS-CoV-2 Variants of Concern Circulating in Hawai'i to Facilitate Public-Health Policies

Using genomics, bioinformatics and statistics, herein we demonstrate the effect of statewide and nationwide quarantine on the introduction of SARS-CoV-2 variants of concern (VOC) in Hawai'i. To define the origins of introduced VOC, we analyzed 260 VOC sequences from Hawai'i, and 301,646 VOC sequences worldwide, deposited in the GenBank and global initiative on sharing all influenza data (GISAID), and constructed phylogenetic trees. The trees define the most recent common ancestor as the origin. Further, the multiple sequence alignment used to generate the phylogenetic trees identified the consensus single nucleotide polymorphisms in the VOC genomes. These consensus sequences allow for VOC comparison and identification of mutations of interest in relation to viral immune evasion and host immune activation. Of note is the P71L substitution within the E protein, the protein sensed by TLR2 to produce cytokines, found in the B.1.351 VOC may diminish the efficacy of some vaccines. Based on the phylogenetic trees, the B.1.1.7, B.1.351, B.1.427, and B.1.429 VOC have been introduced in Hawai'i multiple times since December 2020 from several definable geographic regions. From the first worldwide report of VOC in GenBank and GISAID, to the first arrival of VOC in Hawai'i, averages 320 days with quarantine, and 132 days without quarantine. As such, the effect of quarantine is shown to significantly affect the time to arrival of VOC in Hawai'i, both during and following quarantine. Further, the collective 2020 quarantine of 43-states in the United States demonstrates a profound impact in delaying the arrival of VOC in states that did not practice quarantine, such as Utah. Our data demonstrates that at least 76% of all definable SARS-CoV-2 VOC have entered Hawai'i from California, with the B.1.351 variant in Hawai'i originating exclusively from the United Kingdom. These data provide a foundation for policy-makers and public-health officials to apply precision public health genomics to real-world policies such as mandatory screening and quarantine.

Antibiotic Susceptibility Profiles of Pediococcus pentosaceus from Various Origins and Their Implications for the Safety Assessment of Strains with Food-Technology Applications

ABSTRACT

In the fight against the spread of antibiotic resistance, authorities usually require that strains "intentionally added into the food chain" be tested for their antibiotic susceptibility. This applies to strains used in starter or adjunct cultures for the production of fermented foods, such as many strains of Pediococcus pentosaceus. The European Food Safety Authority recommends testing strains for their antibiotic susceptibility based on both genomic and phenotypic approaches. Furthermore, it proposes a set of antibiotics to assess as well as a list of microbiological cutoffs (MCs), allowing classification of lactic acid bacteria as susceptible or resistant. Accurate MCs are essential not only to avoid false-negative strains, which may carry antibiotic resistance genes and remain unnoticed, but also to avoid false-positive strains, which may be discarded while screening potential candidates for food-technology applications. Because of relatively scarce data, MCs have been defined for the whole Pediococcus genus, although differences between species should be expected. In this study, we investigated the antibiotic susceptibility of 35 strains of P. pentosaceus isolated from various matrices in the past 70 yr. MICs were determined using a standard protocol, and MIC distributions were established. Phenotypic analyses were complemented with genome sequencing and by seeking known antibiotic resistance genes. The genomes of all the strains were free of known antibiotic resistance genes, but most displayed MICs above the currently defined MCs for chloramphenicol, and all showed excessive MICs for tetracycline. Based on the distributions, we calculated and proposed new MCs for chloramphenicol (16 instead of 4 mg/L) and tetracycline (256 instead of 8 mg/L).

HIGHLIGHTS

Streptococcus vicugnae sp. nov., isolated from faeces of alpacas (Vicugna pacos) and cattle (Bos taurus), Streptococcus zalophi sp. nov., and Streptococcus pacificus sp. nov., isolated from respiratory tract of California sea lions (Zalophus californianus)

Four novel independent strains of Streptococcus spp. were isolated from faeces of alpaca (SL1232^T), cattle (KCJ4950), and from respiratory tract of wild California sea lions (CSL7508^T, CSL7591^T). The strains were indole-, oxidase- and catalase-negative, non-spore-forming, non-motile Gram-positive cocci in short and long chains, facultative anaerobes. The 16S rRNA gene of SL1232^T and KCJ4950 shared 99.40-99.60% nucleotide similarity to strains of S. equinus, S. lutetiensis, S. infantarius, and the 16S rRNA gene of CSL7508^T and CSL7591^T demonstrated 98.72 and 98.92% similarity, respectively, to S. marimammalium. All other known Streptococcus species had the 16S rRNA gene sequence similarities of ≤95%. The genomes were sequenced for the novel strains. Average nucleotide identity (ANI) analysis for strains SL1232^T and KCJ4950, showed the highest similarity to S. equinus, S. lutetiensis, and S. infantarius with 85.21, 87.17, 88.47, 85.54, 87.47 and 88.89%, respectively, and strains CSL7508^T and CSL7591^T to S. marimammalium with 87.16 and 83.97%, respectively. Results of ANI were confirmed by pairwise digital DNA-DNA hybridization and phylogeny, which also revealed that the strains belong to three novel species of the genus Streptococcus. Phenotypical features of the novel species were in congruence with closely related members of the genus Streptococcus and gave negative reactions with the tested Lancefield serological groups (A-D, F and G). MALDI-TOF mass spectrometry supported identification of the species. Based on these data, we propose three novel species of the genus Streptococcus, for which the name Streptococcus vicugnae sp. nov. is proposed with the type strain SL1232^T (=NCTC 14341^T=DSM 110741^T=CCUG 74371^T), Streptococcus zalophi sp. nov. is proposed with the type strain CSL7508^T (=NCTC 14410^T=DSM 110742^T=CCUG 74374^T) and Streptococcus pacificus sp. nov. is proposed with the type strain CSL7591^T (=NCTC 14455^T=DSM 111148^T=CCUG 74655^T). The genome G+C content is 36.89, 34.85, and 35.34 % and draft genome sizes are 1906993, 1581094 and 1656080 bp for strains SL1232^T, CSL7508^T, and CSL7591^T, respectively.

A public website for the automated assessment and validation of SARS-CoV-2 diagnostic PCR assays

Summary

Polymerase chain reaction-based assays are the current gold standard for detecting and diagnosing SARS-CoV-2. However, as SARS-CoV-2 mutates, we need to constantly assess whether existing PCR-based assays will continue to detect all known viral strains. To enable the continuous monitoring of SARS-CoV-2 assays, we have developed a web-based assay validation algorithm that checks existing PCR-based assays against the ever-expanding genome databases for SARS-CoV-2 using both thermodynamic and edit-distance metrics. The assay-screening results are displayed as a heatmap, showing the number of mismatches between each detection and each SARS-CoV-2 genome sequence. Using a mismatch threshold to define detection failure, assay performance is summarized with the true-positive rate (recall) to simplify assay comparisons.

Availability and implementation

The assay evaluation website and supporting software are Open Source and freely available at https://covid19.edgebioinformatics.org/#/assayValidation, https://github.com/jgans/thermonucleotide BLAST and https://github.com/LANL-Bioinformatics/assay_validation.

Supplementary information

Supplementary data are available at Bioinformatics online.

The ovarian carcinoma risk with the polymorphisms of CYP1B1 come from the positive selection

Ovarian carcinoma is one of the major causes of gynecological cancer. This study aimed to evaluate the association of CYP1 family polymorphism with the risk of ovarian carcinoma and chemotherapy resistance. Positive selection was detected among human CYP1A1, CYP1A2, and CYP1B1, and other species. Several positive sites were detected by site models and brach-site models. Meta-analysis was conducted for the sites rs1056836 (MAF 0.39) and rs1056827 (MAF 0.36) of CYP1B1 to clarify the association between gene polymorphisms and ovarian carcinoma risk. Subgroup analysis showed the association of rs1056836 polymorphism with ovarian cancer risk among Caucasians and Asians, while all the six genetic models showed no association among African-Americans. All the six genetic models showed no association of rs1056827 polymorphism with ovarian cancer risk. The polymorphisms of rs1056836 associated with ovarian cancer risk were detected in chemotherapy-sensitive and drug-resistant ovarian cancer patients. DNA was extracted from 62 chemotherapy resistance Ovarian carcinoma tissue samples and 137 chemotherapy-sensitive ovarian carcinoma tissue samples as controls. Gene polymorphisms were genotyped using the Sequenom MassARRAY SNP approach. There was no significant association between the CYP1B1 rs1056836 polymorphism and chemotherapy resistance of ovarian cancer in all genetic models. The results suggest that rs1056836 polymorphism of gene CYP1B1 under obvious selection pressure had a significantly increased risk for ovarian carcinoma. However, it had no significant correlation with chemotherapy resistance of ovarian cancer.

A hijack mechanism of Indian SARS-CoV-2 isolates for relapsing contemporary antiviral therapeutics

Coronavirus disease (COVID-19) rapidly expands to a global pandemic and its impact on public health varies from country to country. It is caused by a new virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It is imperative for relapsing current antiviral therapeutics owing to randomized genetic drift in global SARS-CoV-2 isolates. A molecular mechanism behind the emerging genomic variants is not yet understood for the prioritization of selective antivirals. The present computational study was aimed to repurpose existing antivirals for Indian SARS-CoV-2 isolates by uncovering a hijack mechanism based on structural and functional characteristics of protein variants. Forty-one protein mutations were identified in 12 Indian SARS-CoV-2 isolates by analysis of genome variations across 460 genome sequences obtained from 30 geographic sites in India. Two unique mutations such as W6152R and N5928H found in exonuclease of Surat (GBRC275b) and Gandhinagar (GBRC239) isolates. We report for the first time the impact of folding rate on stabilizing/retaining a sequence-structure-function-virulence link of emerging protein variants leading to accommodate hijack ability from current antivirals. Binding affinity analysis revealed the effect of point mutations on virus infectivity and the drug-escaping efficiency of Indian isolates. Emodin and artinemol suggested herein as repurposable antivirals for the treatment of COVID-19 patients infected with Indian isolates. Our study concludes that a protein folding rate is a key structural and evolutionary determinant to enhance the receptor-binding specificity and ensure hijack ability from the prevalent antiviral therapeutics.

Telomeres and a repeat-rich chromosome encode effector gene clusters in plant pathogenic Colletotrichum fungi

Members of the Colletotrichum gloeosporioides species complex are causal agents of anthracnose in many commercially important plants. Closely related strains have different levels of pathogenicity on hosts despite their close phylogenetic relationship. To gain insight into the genetics underlying these differences, we generated and annotated whole-genome assemblies of multiple isolates of C. fructicola (Cf) and C. siamense (Cs), as well as three previously unsequenced species, C. aenigma (Ca), C. tropicale and C. viniferum with different pathogenicity on strawberry. Based on comparative genomics, we identified accessory regions with a high degree of conservation in strawberry-pathogenic Cf, Cs and Ca strains. These regions encode homologs of pathogenicity-related genes known as effectors, organized in syntenic gene clusters, with copy number variations in different strains of Cf, Cs and Ca. Analysis of highly contiguous assemblies of Cf, Cs and Ca revealed the association of related accessory effector gene clusters with telomeres and repeat-rich chromosomes and provided evidence of exchange between these two genomic compartments. In addition, expression analysis indicated that orthologues in syntenic gene clusters showed a tendency for correlated gene expression during infection. These data provide insight into mechanisms by which Colletotrichum genomes evolve, acquire and organize effectors.

Pathoadaptation of the passerine-associated Salmonella enterica serovar Typhimurium lineage to the avian host

Salmonella enterica is a diverse bacterial pathogen and a primary cause of human and animal infections. While many S. enterica serovars present a broad host-specificity, several specialized pathotypes have been adapted to colonize and cause disease in one or limited numbers of host species. The underlying mechanisms defining Salmonella host-specificity are far from understood. Here, we present genetic analysis, phenotypic characterization and virulence profiling of a monophasic S. enterica serovar Typhimurium strain that was isolated from several wild sparrows in Israel. Whole genome sequencing and complete assembly of its genome demonstrate a unique genetic signature that includes the integration of the BTP1 prophage, loss of the virulence plasmid, pSLT and pseudogene accumulation in multiple T3SS-2 effectors (sseJ, steC, gogB, sseK2, and sseK3), catalase (katE), tetrathionate respiration (ttrB) and several adhesion/ colonization factors (lpfD, fimH, bigA, ratB, siiC and siiE) encoded genes. Correspondingly, this strain demonstrates impaired biofilm formation, intolerance to oxidative stress and compromised intracellular replication within non-phagocytic host cells. Moreover, while this strain showed attenuated pathogenicity in the mouse, it was highly virulent and caused an inflammatory disease in an avian host. Overall, our findings demonstrate a unique phenotypic profile and genetic makeup of an overlooked S. Typhimurium sparrow-associated lineage and present distinct genetic signatures that are likely to contribute to its pathoadaptation to passerine birds.

During Salmonella enterica evolution, many different ecological niches have been effectively occupied by this highly diverse bacterial pathogen. While many S. enterica serovars successfully maintained their ability to infect and colonize in a wide-array of host species, a few biotypes have evolved to colonize and cause a disease in only one or a small group of hosts. The evolutionary dynamic and the mechanisms shaping the host-specificity of Salmonella adapted strains are important to better understand Salmonella pathogenicity and its ecology, but still not fully understood. Here, we report genetic and phenotypic characterization of a S. Typhimurium strain that was isolated from several wild sparrows in Israel. This strain presented unique phenotypic profile that included impaired biofilm formation, high sensitivity to oxidative stress and reduced intracellular replication in non-phagocytic cells. In addition, while this strain was able to cause high inflammatory disease in an avian host, it was highly attenuated in the mouse model. Genome analysis identified that specific genetic signatures found in the sparrow strain are more frequently associated with poultry isolates than clinical isolates of S. Typhimurium. These genetic features are expected to accumulatively contribute toward the adaptation of this strain to birds.