PHI-base in 2022: a multi-species phenotype database for Pathogen-Host Interactions

Spatiotemporal differences induced changes in the structure and function of the gut microbiota in an endangered ungulate

The composition and function of animal gut microbiota are shaped by various factors, among which diet is one of the major factors. Diet is affected by seasonal shifts and geographical differences, which in turn impact the host's nutritional levels. To adapt to these environmental changes, the gut microbiome often produces matching responses. Understanding the relationships among the environment, diet, host and the gut microbiome is helpful for exploring the environmental adaptation of wildlife. Here, we chose wild sika deer (Cervus nippon), which is composed natural allopatric populations, to explore how the environment shapes the gut microbiome and affects the relationship between microbiota composition and function and the mutual adaptation of the seasonal living environment to seasonal dietary changes. To this purpose we used DNA metabarcoding, 16S RNA gene amplification sequencing, metagenomic shotgun sequencing and nutritional analyses to comprehensively examine the relationships among the forage plant, nutrient status and host gut microbiome. Our analyses showed spatiotemporal differences in diet between the Tiebu and Hunchun regions, which ultimately led to varying intakes of protein, cellulose, and soluble sugar. The microbiome composition and function showed unique characteristics in each group, and significant differences were detected at the gene level for the protein absorption and metabolism pathway, the carbohydrate metabolic absorption pathway, and cellulase enzyme function, which are related to nutrition. We also found differences in the pathogenic bacteria and resistance mechanisms genes of the gut microbiota in different groups. Our results showed that the gut microbiome of allopatric populations adapts to changes in food composition and nutrition in different seasons and areas to help the host cope with spatiotemporal changes in the living environment. At the same time, varying levels of human activity can have potential health impacts on wild animals.

Whole-genome sequencing and genomic analysis of four Akkermansia strains newly isolated from human feces

Background

Numerous studies have demonstrated that Akkermansia is closely associated with human health. These bacteria colonize the mucus layer of the gastrointestinal tract and utilize mucin as their sole source of carbon and nitrogen. Akkermansia spp. exhibit potential as probiotics under specific conditions. However, the gene accumulation curve derived from pan-genome analysis suggests that the genome of Akkermansia strains remains open. Consequently, current genome mining efforts are insufficient to fully capture the intraspecific and interspecific characteristics of Akkermansia, necessitating continuous exploration of the genomic and phenotypic diversity of new isolates.

Methods

Based on this finding, we sequenced, assembled, and functionally annotated the whole genomes of four new human isolates from our laboratory: AKK-HX001, AKK-HX002, AKK-HX003, and AKK-HX004.

Results

Phylogenetic analysis revealed that all four isolates belonged to the AmII phylogroup, whereas the type strain DSM 22959 is classified within the AmI phylogroup. Moreover, 2,184 shared homologous genes were identified among the four isolates. Functional annotation using the COG, KEGG, and CAZy databases indicated that the functional genes of the four isolates were primarily associated with metabolism. Two antibiotic resistance genes were identified in AKK-HX001 and AKK-HX002, while three resistance genes were detected in AKK-HX003 and AKK-HX004. Additionally, each of the four isolates possessed two virulence genes and three pathogenicity genes, none of which were associated with pathogenicity. The prediction of mobile genetic elements indicated unequal distributions of GIs among the isolates, and a complete CRISPR system was identified in all isolates except AKK-HX003. Two annotated regions of secondary metabolite biosynthesis genes, both belonging to Terpene, were detected using the antiSMASH online tool.

Conclusion

These findings indicate that the four Akkermansia isolates, which belong to a phylogroup distinct from the model strain DSM 22959, exhibit lower genetic risk and may serve as potential probiotic resources for future research.

Evolution of the Black solider fly larvae gut antibiotic resistome during kitchen waste disposal

Kitchen waste (KW) is an important reservoir of antibiotic resistance genes (ARGs). Black solider fly larvae (BSFL) are extensively employed in KW disposal, closely linking to their robust gut microbes. However, antibiotic resistome in BSFL gut during the KW disposal processes and the mechanism remain unclear. In the present study, the antibiotic resistome in BSFL gut within the 12 days KW disposal processes were investigated. Results showed that, ARGs abundance initially increased and subsequently decreased, the five most prevalent core ARG classes were tetracycline, aminoglycoside, cephalosporin, lincosamide and multidrug. A total of 7 MGE types were observed and the horizontal gene transfer (HGT) of ARGs was predominantly mediated by plasmids. Host microbes were mainly categorized into Proteobacteria (98.12 %) and their assemblies were mainly classified into the deterministic processes. To elucidate the driving mechanisms, the mantel test and the structural equation model (SEM) were developed. Results indicated that microbial functions (0.912, p < 0.0001) and microbial community (1.014, p = 0.036), consistently showed very significant relationships with the patterns of ARGs, which presented higher direct effects than indirect effects. Overall, this study makes an initial contribution to a more deepgoing comprehension of the gut antibiotic resistome of BSFL during KW disposal.

Genome Sequencing and Comparative Genomic Analysis of Attenuated Strain Gibellulopsis nigrescens GnVn.1 Causing Mild Wilt in Sunflower

Gibellulopsis nigrescens, previously classified in the Verticillium genus until 2007, is an attenuated pathogen known to provide cross-protection against Verticillium wilt in various crops. To investigate the potential mechanisms underlying its reduced virulence, we conducted genome sequencing, annotation, and a comparative genome analysis of G. nigrescens GnVn.1 (GnVn.1), an attenuated strain isolated from sunflower. The genome sequencing and annotation results revealed that the GnVn.1 genome consists of 22 contigs, with a total size of 31.79 Mb. We predicted 10,876 genes, resulting in a gene density of 342 genes per Mb. The pathogenicity gene prediction results indicated 1733 high-confidence pathogenicity factors (HCPFs), 895 carbohydrate-active enzymes (CAZys), and 359 effectors. Moreover, we predicted 40 secondary metabolite clusters (SMCs). The comparative genome analysis indicated that GnVn.1 contains more CAZys, SMCs, predicted effectors, and HCPF genes than Verticillium dahliae (VdLs.17) and Verticillium alfalfae (VaMas.102). The core-pan analysis results showed that GnVn.1 had more specific HCPFs, effectors, CAZys, and secreted protein (SP) genes, and lost many critical pathogenic genes compared to VdLs.17 and VaMs.102. Our results indicate that the GnVn.1 genome harbors more pathogenicity-related genes than the VdLs.17 and VaMs.102 genomes. These abundant genes may play critical roles in regulating virulence. The loss of critical pathogenic genes causes weak virulence and confers biocontrol strategies to GnVn.1.

Ensembl 2025

Ensembl (www.ensembl.org) is an open platform integrating publicly available genomics data across the tree of life with a focus on eukaryotic species related to human health, agriculture and biodiversity. This year has seen a continued expansion in the number of species represented, with >4800 eukaryotic and >31 300 prokaryotic genomes available. The new Ensembl site, currently in beta, has continued to develop, currently holding >2700 eukaryotic genome assemblies. The new site provides genome, gene, transcript, homology and variation views, and will replace the current Rapid Release site; this represents a key step towards provision of a single integrated Ensembl site. Additional activities have included developing improved regulatory annotation for human, mouse and agricultural species, and expanding the Ensembl Variant Effect Predictor tool. To learn more about Ensembl, help and documentation are available along with an extensive training program that can be accessed via our training pages.

Comparative Genomic Analysis of an Apiotrichum cacaoliposimilis Strain Isolated from a Patient with Urinary Tract Infection

Opportunistic infections caused by fungi, particularly those occurring in immunocompromised patients, are considered challenging worldwide. Therefore, a comprehensive understanding of pathogenic fungi is necessary. The present study reports the isolation of a strain of Apiotrichum cacaoliposimilis, which is difficult to detect using conventional clinical assays, from the sterile urine samples of a patient with a urinary tract infection. Sanger sequencing of the internal transcribed spacer regions confirmed the genus of the microbe, while whole-genome sequencing yielded the initial genome assembly of A. cacaoliposimilis. A total of 7,161 predicted proteincoding genes were mapped using multiple databases, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, non-redundant protein database, Pathogen-Host Interactions Database, and Comprehensive Antibiotic Resistance Database. The phenotypic data, biochemical reactions, and antimicrobial susceptibility analyses were conducted to reveal the metabolic properties, virulence, and drug resistance profile of the isolated A. cacaoliposimilis. The rank-sum test revealed the differences in the intergeneric distribution of the highly virulent genes UgeB and Pem2. In addition, other genes exhibited significant overlap in terms of virulence factors with the clinical isolate Apiotrichum mycotoxinivorans GMU1709. Fortunately, similar to most fungi belonging to the Apiotrichum genus, the isolate investigated in the present study was also sensitive to the drug voriconazole (MIC = 0.06 μg/ml). In summary, the phylogenetic placement, potential pathogenic genes, drug sensitivity patterns, and morphological characteristics of the isolated A. cacaoliposimilis were determined precisely in the present study.

Mapping an avirulence gene in the sunflower parasitic weed Orobanche cumana and characterization of host selection based on virulence alleles

BACKGROUND

Sunflower broomrape (Orobanche cumana Wallr.) is a holoparasitic plant that jeopardizes sunflower production in most areas of Europe and Asia. Recently, populations with increased virulence, classified as race GGV, have been identified in Southern Spain's Guadalquivir Valley gene pool. These populations overcome resistance genes in hybrids resistant to the predominant race FGV. This study aimed to (i) determine the inheritance and map the avirulence trait segregating in a cross between O. cumana individuals from populations EK23 (FGV) and IN201 (GGV), and (ii) characterize the host effect on the IN201 parental population allelic diversity.

RESULTS

A segregating population consisting of 144 F2:3 families was evaluated for virulence using a differential sunflower genotype (Hybrid 1, resistant to race FGV and susceptible to race GGV) and genotyped with SNP markers. The ratio of avirulent to virulent F2:3 families was not significantly different to 1:3 (χ2 = 0.93; P = 0.34), indicating monogenic control of the avirulence/virulence trait. The AvrG-GV locus was mapped on the upper end of O. cumana chromosome 2, 9.2 cM distal from the SNP markers OS04791 and OS02805. Secretome analysis in the AvrG-GV region revealed a cysteine-rich CAP superfamily- and a glucan 1,3-beta-glucosidase family 3-encoding genes as possible candidates for AvrG-GV. SNP allelic analysis on the IN201 population parasitizing a highly susceptible genotype or the differential genotype Hybrid 1 showed that (i) IN201 structure was shaped towards virulent alleles at SNP loci linked to AvrG-GV (ii) there were significant allelic frequency differences associated with the host genotype at AvrG-GV-linked loci.

CONCLUSIONS

This study mapped for the first time an avirulence gene in parasitic plants using a classical genetic approach, confirmed a gene-for-gene model in the O.cumana -sunflower system, and showed the implication of this single avirulence gene in determining the structure of broomrape populations subjected to selection pressure posed by a resistant genotype. The results will contribute to a better understanding of the interaction between crops and weedy parasitic plants, and to effectively manage evolution of virulence by sustainable control strategies based on host genetic resistance.

Genomic Features of Taiwanofungus gaoligongensis and the Transcriptional Regulation of Secondary Metabolite Biosynthesis

Fungal secondary metabolites (SMs) have broad applications in biomedicine, biocontrol, and the food industry. In this study, whole-genome sequencing and annotation of Taiwanofungus gaoligongensis were conducted, followed by comparative genomic analysis with 11 other species of Polyporales to examine genomic variations and secondary metabolite biosynthesis pathways. Additionally, transcriptome data were used to analyze the differential expression of polyketide synthase (PKS), terpene synthase (TPS) genes, and transcription factors (TFs) under different culture conditions. The results show that T. gaoligongensis differs from other fungal species in genome size (34.58 Mb) and GC content (50.72%). The antibiotics and Secondary Metabolites Analysis Shell (AntiSMASH) analysis reveals significant variation in the number of SM biosynthetic gene clusters (SMBGCs) across the 12 species (12-29), with T. gaoligongensis containing 25 SMBGCs: 4 PKS, 6 non-ribosomal peptide synthetase (NRPS), and 15 TPS clusters. The TgPKS1 gene is hypothesized to be involved in the biosynthesis of orsellinic acid or its derivatives, while TgPKS2 might catalyze the synthesis of 6-methylsalicylic acid (6MSA) and its derivatives. The TgTRI5 genes are suggested to synthesize tetracyclic sesquiterpene type B trichothecene compounds, while TgPentS may be involved in the synthesis of δ-cadinol, β-copaene, and α-murolene analogs or derivatives. Comparative genomic analysis shows that the genome size of T. gaoligongensis is similar to that of T. camphoratus, with comparable SMs. Both species share four types of PKS domains and five distinct types of TPS. Additionally, T. gaoligongensis exhibits a high degree of similarity to Laetiporus sulphureus, despite belonging to a different genus within the same family. Transcriptome analysis reveals significant variation in the expression levels of PKS and TPS genes across different cultivation conditions. The TgPKS1 and TgPKS4 genes, along with nine TgTFs, are significantly upregulated under three solid culture conditions. In contrast, under three different liquid culture conditions, the TgPKS3, TgTRI5-1, and TgTRI5-2 genes, along with twelve TgTFs, exhibit higher activity. Co-expression network analysis and TgTFs binding site prediction in the promoter regions of TgPKS and TgTPS genes suggest that TgMYB9 and TgFTD4 regulate TgPKS4 expression. TgHOX1, TgHSF2, TgHSF3, and TgZnF4 likely modulate TgPKS3 transcriptional activity. TgTRI5-1 and TgTRI5-5 expression is likely regulated by TgbZIP2 and TgZnF15, respectively. This study provides new insights into the regulatory mechanisms of SMs in T. gaoligongensis and offers potential strategies for enhancing the biosynthesis of target compounds through artificial intervention.

PHI-base - the multi-species pathogen-host interaction database in 2025

The Pathogen-Host Interactions Database (PHI-base) has, since 2005, provided manually curated genes from fungal, bacterial and protist pathogens that have been experimentally verified to have important pathogenicity, virulence and/or effector functions during different types of interactions involving human, animal, plant, invertebrate and fungal hosts. PHI-base provides phenotypic annotation and genotypic information for both native and model host interactions, including gene alterations that do not alter the phenotype of the interaction. In this article, we describe major updates to PHI-base. The latest version of PHI-base, 4.17, contains a 19% increase in genes and a 23% increase in interactions relative to version 4.12 (released September 2022). We also describe the unification of data in PHI-base 4 with the data curated from a new curation workflow (PHI-Canto), which forms the first complete release of PHI-base version 5.0. Additionally, we describe adding support for the Frictionless Data framework to PHI-base 4 datasets, new ways of sharing interaction data with the Ensembl database, an analysis of the conserved orthologous genes in PHI-base, and the increasing variety of research studies that make use of PHI-base. PHI-base version 4.17 is freely available at www.phi-base.org and PHI-base version 5.0 is freely available at phi5.phi-base.org.

Natal factors influencing newborn's oral microbiome diversity

The early microbiota of neonates is crucial for developing the postnatal immune system and establishing normal physiological, metabolic, and neurological functions. This study aimed to investigate the factors influencing the diversity of the neonatal oral microbiome, including mother-to-newborn microbial transmission. The study includes a prospective cohort comprising 73 mothers and 87 neonates and a retrospective cohort comprising 991 mothers and 1,121 neonates. Samples from the maternal cervix and neonatal gastric, bronchial, and oral cavities were analyzed using culture-based methods. Neonatal oral swab samples were also analyzed using 16S rRNA gene sequencing to characterize microbial diversity and composition. Similar genera were detected in the neonatal gastric, bronchial, and oral samples, and the neonatal gastric culture was the most similar to the maternal cervical culture. In addition to mother-to-newborn microbial transmission, various natal factors including birth type, birth weight, delivery mode, maternal chorioamnionitis, maternal diabetes and the presence of microbes in other sites influenced neonatal oral microbiome diversity. Among these factors, the birth type was the most significant, and preterm neonates exhibited decreased oral microbiome diversity, with fewer beneficial bacteria and more pathogens. These findings could serve as a baseline for research on the establishment of the oral microbiota in preterm neonates and its health implications.

Discovery of Alanomyces manoharacharyi: A Novel Fungus Identified Using Genome Sequencing and Metabolomic Analysis

In this study, a new species of Alanomyces was isolated as an endophyte from the bark of Azadirachta indica from Mulshi, Maharashtra. The identity of this isolate was confirmed based on the asexual morphological characteristics as well as multi-gene phylogeny based on the internal transcribed spacer (ITS) and large subunit (LSU) nuclear ribosomal RNA (rRNA) regions. As this was the second species to be reported in this genus, we sequenced the genome of this species to increase our knowledge about the possible applicability of this genus to various industries. Its genome length was found to be 35.01 Mb, harboring 7870 protein-coding genes as per Augustus and 8101 genes using GeMoMa. Many genes were annotated using the Clusters of Orthologous Groups (COGs) database, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), Swiss-Prot, NCBI non-redundant nucleotide sequences (NTs), and NCBI non-redundant protein sequences (NRs). The number of repeating sequences was predicted using Proteinmask and RepeatMasker; tRNA were detected using tRNAscan and snRNA were predicted using rfam_scan. The genome was also annotated using the Pathogen-Host Interactions Database (PHI-base) and AntiSMASH. To confirm the evolutionary history, average nucleotide identity (ANIb), phylogeny based on orthologous proteins, and single nucleotide polymorphisms (SNPs) were carried out. Metabolic profiling of the methanolic extract of dried biomass and ethyl acetate extract of the filtrate revealed a variety of compounds of great importance in the pharmaceutical and cosmetic industry. The characterization and genomic analysis of the newly discovered species Alanomyces manoharacharyi highlights its potential applicability across multiple industries, particularly in pharmaceuticals and cosmetics due to its diverse secondary metabolites and unique genetic features it possesses.

Genomic Characterization and Comparative Analysis of Streptococcus zhangguiae sp. nov. Isolated from the Respiratory Tract of Marmota Himalayana

Two Gram-stain-positive, oxidase-negative, non-motile, facultative anaerobic, α-hemolytic, coccus-shaped bacteria (zg-86T and zg-70) were isolated from the respiratory tracts of marmots (Marmota Himalayana) on the Qinghai-Tibet Plateau of China. Phylogenetic analysis of the 16S rRNA gene and 545 core genes revealed that these two strains belong to the Streptococcus genus. These strains were most closely related to Streptococcus respiraculi HTS25T, Streptococcus cuniculi CCUG 65085T, and Streptococcus marmotae HTS5T. The average nucleotide identity (ANI) and digital DNA‒DNA hybridization (dDDH) were below the threshold for species delineation. The predominant cellular fatty acids (CFAs) in this novel species were C16:0, C18:0, and C18:1ω9c, whereas the primary polar lipids were phosphatidylglycerol (PG), phosphatidylethanolamine (PE) and an unknown phosphoglycolipid (PGL). The optimal growth conditions for the strains were 37 °C, pH 7.0, and 0.5% (w/v) NaCl on brain-heart infusion (BHI) agar supplemented with 5% defibrinated sheep blood. Comparative genomics analyses revealed the potential pathogenicity of strain zg-86T through comparisons with suis subclade strains in terms of virulence factors, pathogen-host interactions (PHIs) and mobile genetic factors (MGEs). Based on the phenotypic characteristics and phylogenetic analyses, we propose that these two isolates represent novel species in the genus Streptococcus, for which the names Streptococcus zhangguiae sp. nov. (the type strain zg-86T=GDMCC 1.1758T=JCM 34273T) is proposed.

Genomic insights into fish pathogenic bacteria: A systems biology perspective for sustainable aquaculture

Fish diseases significantly challenge global aquaculture, causing substantial financial losses and impacting sustainability, trade, and socioeconomic conditions. Understanding microbial pathogenesis and virulence at the molecular level is crucial for disease prevention in commercial fish. This review provides genomic insights into fish pathogenic bacteria from a systems biology perspective, aiming to promote sustainable aquaculture. It covers the genomic characteristics of various fish pathogens and their industry impact. The review also explores the systems biology of zebrafish, fish bacterial pathogens, and probiotic bacteria, offering insights into fish production, potential vaccines, and therapeutic drugs. Genome-scale metabolic models aid in studying pathogenic bacteria, contributing to disease management and antimicrobial development. Researchers have also investigated probiotic strains to improve aquaculture health. Additionally, the review highlights bioinformatics resources for fish and fish pathogens, which are essential for researchers. Systems biology approaches enhance understanding of bacterial fish pathogens by revealing virulence factors and host interactions. Despite challenges from the adaptability and pathogenicity of bacterial infections, sustainable alternatives are necessary to meet seafood demand. This review underscores the potential of systems biology in understanding fish pathogen biology, improving production, and promoting sustainable aquaculture.

Distinct species turnover patterns shaped the richness of antibiotic resistance genes on eight different microplastic polymers

Elucidating the formation mechanism of plastisphere antibiotic resistance genes (ARGs) on different polymers is necessary to understand the ecological risks of plastisphere ARGs. Here, we explored the turnover and assembly mechanism of plastisphere ARGs on 8 different microplastic polymers (4 biodegradable (bMPs) and 4 non-biodegradable microplastics (nMPs)) by metagenomic sequencing. Our study revealed the presence of 479 ARGs with abundance ranging from 41.37 to 58.17 copies/16S rRNA gene in all plastispheres. These ARGs were predominantly multidrug resistance genes. The richness of plastisphere ARGs on different polymers had a significant correlation with the contribution of species turnover to plastisphere ARGs β diversity. Furthermore, polymer type was the most critical factor affecting the composition of plastisphere ARGs. More opportunistic pathogens carrying diverse ARGs on BMPs (PBAT, PBS, and PHA) with higher horizontal gene transfer potential may further magnify the ecological risks and human health threats. For example, the opportunistic pathogens Riemerella anatipestifer, Vibrio campbellii, and Vibrio cholerae are closely related to human production and life, which were the important potential hosts of many plastisphere ARGs and mobile genetic elements on BMPs. Thus, we emphasize the urgency of developing the formation mechanism of plastisphere ARGs and the necessity of controlling BMPs and ARG pollution, especially BMPs, with ever-increasing usage in daily life.

First draft genome sequence of a Pectobacterium polaris strain isolated in South Africa from potato tuber affected by soft rot

A phytopathogenic strain of Pectobacterium polaris (designated SRB2) was isolated for the first time in South Africa from a potato tuber affected by soft rot. The draft genome of strain SRB2 encodes various plant cell wall-degrading enzymes and genes associated with biofilm formation and virulence. Antibiotic resistance genes were not detected.

Complete genome identified of clinical isolate Prototheca

Introduction. Prototheca is an opportunistic pathogen that can infect both humans and animals, of which Prototheca wickerhamii (P. wickerhamii) being the most significant pathogenic green algae.Gap statement. The incidence of human diseases caused by Prototheca has been on the rise, yet there is a significant gap in genetic research pertaining to the pathophysiological aspects of these infections.Aim. The aim of this study is to present the whole genome data from the clinical isolate InPu-22_FZ strain and to understand its genomic characteristics through comparative genomic analysis and phylogenetic tree analysis. Functional annotation of protein-coding genes and analysis of their pathogenicity are also conducted.Methodology. We described the high-quality de novo genome assembly of the clinical isolate InPu-22_FZ strain, achieved by combining Nanopore ONT and Illumina NovaSeq sequencing technologies. Phylogenetic tree was constructed to study the evolutionary relationship between the InPu-22_FZ strain and other species. The average nucleotide identity (ANI) analysis was used to assess the similarity between different species. Additionally, the size, distribution and composition of synteny blocks were also analysed to infer the evolutionary relationships of the genomes.Results. The size of the assembled nuclear genome was 18.47 Mb with 48 contigs. Key features of the genome include high overall GC content (63.31%), high number (5478) and proportion (62.24%) of protein-coding genes and more than 96.71% of genes annotated for gene function. Phylogenetic analyses showed that the InPu-22_FZ strain and other P. wickerhamii clustered into one clade with a bootstrap value of 99% and collinearity analysis revealed high levels of collinearity with ATCC 16529. The ANI analysis revealed only a relatively high similarity (89-93%) to available P. wickerhamii genomes, suggesting the overall genomic novelty of InPu-22_FZ strain. Interestingly, the analysis of the pathogen-host interaction database unveiled and demonstrated reduced virulence of this strain, albeit it was isolated from a chronic upper-limb cutaneous infection.Conclusion. The study provides an in-depth insight into the genomic structure and biological function of the InPu-22_FZ strain, revealing the genetic basis of its pathogenicity and virulence.

Comparative Genomic Analyses of Colletotrichum lindemuthianum Pathotypes with Different Virulence Levels and Lifestyles

Colletotrichum lindemuthianum is the most frequent pathogenic fungus of the common bean Phaseolus vulgaris. This filamentous fungus employs a hemibiotrophic nutrition/infection strategy, which is characteristic of many Colletotrichum species. Due to host-pathogen coevolution, C. lindemuthianum includes pathotypes with a diversity of virulence against differential common bean varieties. In this study, we performed comparative genomic analyses on three pathotypes with different virulence levels and a non-pathogenic pathotype, isolated from different geographical areas in Mexico. Our results revealed large genomes with high transposable element contents that have undergone expansions, generating intraspecific diversity. All the pathotypes exhibited a similar number of clusters of orthologous genes (COGs) and Gene Ontology (GO) terms. TFomes contain families that are typical in fungal genomes; however, they show different contents between pathotypes, mainly in transcription factors with the fungal-specific TF and Zn2Cys6 domains. Peptidase families mainly contain abundant serine peptidases, metallopeptidases, and cysteine peptidases. In the secretomes, the number of genes differed between the pathotypes, with a high percentage of candidate effectors. Both the virulence gene and CAZyme gene content for each pathotype was abundant and diverse, and the latter was enriched in hemicellulolytic enzymes. We provide new insights into the nature of intraspecific diversity among C. lindemuthianum pathotypes and the origin of their ability to rapidly adapt to genetic changes in its host and environmental conditions.

Pannonibacter tanglangensis sp. nov., a New Species Isolated from Pond Sediment

Two bacterial strains (XCT-34T and XCT-53) isolated from sediment samples of an artificial freshwater reservoir were analyzed using a polyphasic approach. The two isolates are aerobic, Gram-stain-negative, oxidase-negative, catalase-positive, motile with polar flagella, rod-shaped, and approximately 1.4-3.4 × 0.4-0.9 μm in size. Phylogenetic analyses based on 16S rRNA gene and whole-genome sequences showed that the two strains formed a distinct branch within the evolutionary radiation of the genus Pannonibacter, closest to Pannonibacter carbonis Q4.6T (KCTC 52466). Furthermore, lower than threshold average nucleotide identity values (ANI, 85.7-86.4%) and digital DNA-DNA hybridization values (dDDH, 22.3-30.5%) of the two strains compared to the nearest type strains also confirmed that they represented a novel species. Genomic analyses, including annotation of the KEGG pathways, prediction of the secondary metabolism biosynthetic gene clusters and PHI phenotypes, supported functional inference and differentiation of the strains from the closely related taxa. Results of chemotaxonomic and physiological studies revealed that their distinct phenotypic characteristics distinguished them from existing Pannonibacter species. Thus, the two strains are considered to represent a novel species of Pannonibacter, for which the name of Pannonibacter tanglangensis sp. nov. is proposed, with XCT-34T (= KCTC 82332T = GDMCC 1.1947T) as the respective type strain.

The study of iron- and copper-binding proteome of Fusarium oxysporum and its effector candidates

Fusarium oxysporum f.sp. lycopersici is a phytopathogen which causes vascular wilt disease in tomato plants. The survival tactics of both pathogens and hosts depend on intricate interactions between host plants and pathogenic microbes. Iron-binding proteins (IBPs) and copper-binding proteins (CBPs) play a crucial role in these interactions by participating in enzyme reactions, virulence, metabolism, and transport processes. We employed high-throughput computational tools at the sequence and structural levels to investigate the IBPs and CBPs of F. oxysporum. A total of 124 IBPs and 37 CBPs were identified in the proteome of Fusarium. The ranking of amino acids based on their affinity for binding with iron is Glu > His> Asp > Asn > Cys, and for copper is His > Asp > Cys respectively. The functional annotation, determination of subcellular localization, and Gene Ontology analysis of these putative IBPs and CBPs have unveiled their potential involvement in a diverse array of cellular and biological processes. Three iron-binding glycosyl hydrolase family proteins, along with four CBPs with carbohydrate-binding domains, have been identified as potential effector candidates. These proteins are distinct from the host Solanum lycopersicum proteome. Moreover, they are known to be located extracellularly and function as enzymes that degrade the host cell wall during pathogen-host interactions. The insights gained from this report on the role of metal ions in plant-pathogen interactions can help develop a better understanding of their fundamental biology and control vascular wilt disease in tomato plants.

Annotation of Functions of Sequences of Concern and Its Relevance to the New Biosecurity Regulatory Framework in the United States

Introduction

Recent regulations from United States Government agencies reshape the screening of synthetic nucleic acids. These take a step away from categorizing hazard on the basis of "bad" taxa and invoke the function of the sequence in pathogenesis or intoxication. Ascertaining functions related to pathogenesis and distinguishing these from other molecular abilities that are unproblematic is not simple. Some have suggested that this information can be readily obtained from existing databases of pathogens.

Objectives

We evaluate how virulence factors are described in current databases of pathogens and their adequacy for biothreat data science. We discuss limitations of how virulence factors have been conceived and propose using the sequence of concern (SoC) term to distinguish sequences with biothreat from those without. We discuss ways in which databases of SoCs might be implemented for research and regulatory purposes. We describe ongoing work improving functional descriptions of SoCs.

Methods

We assess the adequacy of descriptions of virulence factors in pathogen databases following extensive engagement with the literature in microbial pathogenesis.

Results/Conclusions

Descriptions of virulence factors in pathogen databases are inadequate for understanding biothreats. Many are not biothreats and would not be concerning if transferred to another pathogen. New gene ontology terms have been authored, and those specific to pathogenic viral processes are being generalized to make them relevant to other pathogenic taxa. This allows better understanding by humans and better recognition by machines. A database of annotated functions of SoCs could benefit the evolving biosecurity regulatory framework in the United States.

Definition of the effector landscape across 13 phytoplasma proteomes with LEAPH and EffectorComb

'Candidatus Phytoplasma' genus, a group of fastidious phloem-restricted bacteria, can infect a wide variety of both ornamental and agro-economically important plants. Phytoplasmas secrete effector proteins responsible for the symptoms associated with the disease. Identifying and characterizing these proteins is of prime importance for expanding our knowledge of the molecular bases of the disease. We faced the challenge of identifying phytoplasma's effectors by developing LEAPH, a machine learning ensemble predictor composed of four models. LEAPH was trained on 479 proteins from 53 phytoplasma species, described by 30 features. LEAPH achieved 97.49% accuracy, 95.26% precision and 98.37% recall, ensuring a low false-positive rate and outperforming available state-of-the-art methods. The application of LEAPH to 13 phytoplasma proteomes yields a comprehensive landscape of 2089 putative pathogenicity proteins. We identified three classes according to different secretion models: 'classical', 'classical-like' and 'non-classical'. Importantly, LEAPH identified 15 out of 17 known experimentally validated effectors belonging to the three classes. Furthermore, to help the selection of novel candidates for biological validation, we applied the Self-Organizing Maps algorithm and developed a Shiny app called EffectorComb. LEAPH and the EffectorComb app can be used to boost the characterization of putative effectors at both computational and experimental levels, and can be employed in other phytopathological models.

Genome sequencing of Elaeocarpus spp. stem blight pathogen Pseudocryphonectria elaeocarpicola reveals potential adaptations to colonize woody bark

BACKGROUND

Elaeocarpus spp. stem blight, caused by Pseudocryphonectria elaeocarpicola, is a destructive disease, which will significantly reduce the productivity and longevity of Elaeocarpus spp. plants, especially in the Guangdong Province of China. However, few information is available for P. elaeocarpicola. To unravel the potential adaptation mechanism of stem adaptation, the whole genome of P. elaeocarpicola was sequenced by using the DNBSEQ and PacBio platforms.

RESULTS

P. elaeocarpicola harbors 44.49 Mb genome with 10,894 predicted coding genes. Genome analysis revealed that the P. elaeocarpicola genome encodes a plethora of pathogenicity-related genes. Analysis of carbohydrate-active enzymes (CAZymes) revealed a rich variety of enzymes participated in plant cell wall degradation, which could effectively degrade cellulose, hemicellulose and xyloglucans in the plant cell wall and promote the invasion of the host plant. There are 213 CAZyme families found in P. elaeocarpicola, among which glycoside hydrolase (GH) family has the largest number, far exceeding other tested fungi by 53%. Besides, P. elaeocarpicola has twice as many genes encoding chitin and cellulose degradation as Cryphonectria parasitica, which belong to the same family. The predicted typical secreted proteins of P. elaeocarpicola are numerous and functional, including many known virulence effector factors, indicating that P. elaeocarpicola has great potential to secrete virulence effectors to promote pathogenicity on host plants. AntiSMASH revealed that the genome encoded 61 secondary metabolic gene clusters including 86 secondary metabolic core genes which was much higher than C. parasitica (49). Among them, two gene cluster of P. elaeocarpicola, cluster12 and cluster52 showed 100% similarity with the mycotoxins synthesis clusters from Aspergillus steynii and Alternaria alternata, respectively. In addition, we annotated cytochrome P450 related enzymes, transporters, and transcription factors in P. elaeocarpicola, which are important virulence determinants of pathogenic fungi.

CONCLUSIONS

Taken together, our study represents the first genome assembly for P. elaeocarpicola and reveals the key virulence factors in the pathogenic process of P. elaeocarpicola, which will promote our understanding of its pathogenic mechanism. The acquired knowledge lays a foundation for further exploration of molecular interactions with the host and provide target for management strategies in future research.

Bile acid-gut microbiota imbalance in cholestasis and its long-term effect in mice

Cholestasis is a common morbid state that may occur in different phases; however, a comprehensive evaluation of the long-term effect post-recovery is still lacking. In the hepatic cholestasis mouse model, which was induced by a temporary complete blockage of the bile duct, the stasis of bile acids and liver damage typically recovered within a short period. However, we found that the temporary hepatic cholestasis had a long-term effect on gut microbiota dysbiosis, including overgrowth of small intestinal bacteria, decreased diversity of the gut microbiota, and an overall imbalance in its composition accompanied by an elevated inflammation level. Additionally, we observed an increase in Escherichia-Shigella (represented by ASV136078), rich in virulence factors, in both small and large intestines following cholestasis. To confirm the causal role of dysregulated gut microbiota in promoting hepatic inflammation and injury, we conducted gut microbiota transplantation into germ-free mice. We found that recipient mice transplanted with feces from cholestasis mice exhibited liver inflammation, damage, and accumulation of hepatic bile acids. In conclusion, our study demonstrates that cholestasis disrupts the overall load and structural composition of the gut microbiota in mice, and these adverse effects persist after recovery from cholestatic liver injury. This finding suggests the importance of monitoring the structural composition of the gut microbiota in patients with cholestasis and during their recovery.

IMPORTANCE

Our pre-clinical study using a mouse model of cholestasis underscores that cholestasis not only disrupts the equilibrium and structural configuration of the gut microbiota but also emphasizes the persistence of these adverse effects even after bile stasis restoration. This suggests the need of monitoring and initiating interventions for gut microbiota structural restoration in patients with cholestasis during and after recovery. We believe that our study contributes to novel and better understanding of the intricate interplay among bile acid homeostasis, gut microbiota, and cholestasis-associated complications. Our pre-clinical findings may provide implications for the clinical management of patients with cholestasis.

Decoding the Arsenal: Protist Effectors and Their Impact on Photosynthetic Hosts

Interactions between various microbial pathogens including viruses, bacteria, fungi, oomycetes, and their plant hosts have traditionally been the focus of phytopathology. In recent years, a significant and growing interest in the study of eukaryotic microorganisms not classified among fungi or oomycetes has emerged. Many of these protists establish complex interactions with photosynthetic hosts, and understanding these interactions is crucial in understanding the dynamics of these parasites within traditional and emerging types of farming, including marine aquaculture. Many phytopathogenic protists are biotrophs with complex polyphasic life cycles, which makes them difficult or impossible to culture, a fact reflected in a wide gap in the availability of comprehensive genomic data when compared to fungal and oomycete plant pathogens. Furthermore, our ability to use available genomic resources for these protists is limited by the broad taxonomic distance that these organisms span, which makes comparisons with other genomic datasets difficult. The current rapid progress in genomics and computational tools for the prediction of protein functions and interactions is revolutionizing the landscape in plant pathology. This is also opening novel possibilities, specifically for a deeper understanding of protist effectors. Tools like AlphaFold2 enable structure-based function prediction of effector candidates with divergent protein sequences. In turn, this allows us to ask better biological questions and, coupled with innovative experimental strategies, will lead into a new era of effector research, especially for protists, to expand our knowledge on these elusive pathogens and their interactions with photosynthetic hosts. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Transcriptome analysis of Phytophthora cactorum infecting strawberry identified RXLR effectors that induce cell death when transiently expressed in Nicotiana benthamiana

Phytophthora cactorum is a plant pathogenic oomycete that causes crown rot in strawberry leading to significant economic losses every year. To invade the host, P. cactorum secretes an arsenal of effectors that can manipulate host physiology and impair its defense system promoting infection. A transcriptome analysis was conducted on a susceptible wild strawberry genotype (Fragaria vesca) 48 hours post inoculation with P. cactorum to identify effectors expressed during the early infection stage. The analysis revealed 4,668 P. cactorum genes expressed during infection of F. vesca. A total of 539 secreted proteins encoded by transcripts were identified, including 120 carbohydrate-active enzymes, 40 RXLRs, 23 proteolytic enzymes, nine elicitins, seven cysteine rich proteins, seven necrosis inducing proteins and 216 hypothetical proteins with unknown function. Twenty of the 40 RXLR effector candidates were transiently expressed in Nicotiana benthamiana using agroinfiltration and five previously unreported RXLR effector genes (Pc741, Pc8318, Pc10890, Pc20813, and Pc22290) triggered cell death when transiently expressed. The identified cell death inducing RXLR effectors showed 31-66% identity to known RXLR effectors in different Phytophthora species having roles in pathogenicity including both activation and suppression of defense response in the host. Furthermore, homology analysis revealed that these cell death inducing RXLR effectors were highly conserved (82 - 100% identity) across 23 different strains of P. cactorum originating from apple or strawberry.

Mobile genetic elements mediate the cross-media transmission of antibiotic resistance genes from pig farms and their risks

Manure composting in traditional small-scale pig farms leads to the migration and diffusion of antibiotics and antibiotics resistance genes (ARGs) along the chain of transmission to the surrounding environment, increasing the risk of environmental resistance. Understanding the transmission patterns, driving factors, and health risks of ARGs on small-scale pig farms is important for effective control of ARGs transmission. This study was conducted on a small pig farm and its surrounding environment. The cross-media transmission of ARGs and their risks in the farming habitat were investigated using Metagenomic annotation and qPCR quantitative detection. The results indicate that ARGs in farms spread with manure pile-soil-channel sediment-mudflat sediment. Pig farm manure contributed 22.49 % of the mudflat sediment ARGs. Mobile genetic elements mediate the spread of ARGs across different media. Among them, tnpA and IS26 have the highest degree. Transmission of high-risk ARGs sul1 and tetM resulted in a 50 % and 116 % increase in host risk for sediment, respectively. This study provides a basis for farm manure management and control of the ARGs spread.

What is new in FungiDB: a web-based bioinformatics platform for omics-scale data analysis for fungal and oomycete species

FungiDB (https://fungidb.org) serves as a valuable online resource that seamlessly integrates genomic and related large-scale data for a wide range of fungal and oomycete species. As an integral part of the VEuPathDB Bioinformatics Resource Center (https://veupathdb.org), FungiDB continually integrates both published and unpublished data addressing various aspects of fungal biology. Established in early 2011, the database has evolved to support 674 datasets. The datasets include over 300 genomes spanning various taxa (e.g. Ascomycota, Basidiomycota, Blastocladiomycota, Chytridiomycota, Mucoromycota, as well as Albuginales, Peronosporales, Pythiales, and Saprolegniales). In addition to genomic assemblies and annotation, over 300 extra datasets encompassing diverse information, such as expression and variation data, are also available. The resource also provides an intuitive web-based interface, facilitating comprehensive approaches to data mining and visualization. Users can test their hypotheses and navigate through omics-scale datasets using a built-in search strategy system. Moreover, FungiDB offers capabilities for private data analysis via the integrated VEuPathDB Galaxy platform. FungiDB also permits genome improvements by capturing expert knowledge through the User Comments system and the Apollo genome annotation editor for structural and functional gene curation. FungiDB facilitates data exploration and analysis and contributes to advancing research efforts by capturing expert knowledge for fungal and oomycete species.

Genomic insights into the cellular specialization of predation in raptorial protists

BACKGROUND

Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited.

RESULTS

In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria.

CONCLUSIONS

Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

Analyzing the safety of the parasiticide fungus Mucor circinelloides: first insights on its virulence profile and interactions with the avian gut microbial community

Parasiticide fungi are considered an accurate, sustainable, and safe solution for the biocontrol of animal gastrointestinal (GI) parasites. This research provides an initial characterization of the virulence of the native parasiticide fungus Mucor circinelloides (FMV-FR1) and an assessment of its impact on birds' gut microbes. The genome of this fungus was sequenced to identify the genes coding for virulence factors. Also, this fungus was checked for the phenotypic expression of proteinase, lecithinase, DNase, gelatinase, hemolysin, and biofilm production. Finally, an in vivo trial was developed based on feeding M. circinelloides spores to laying hens and peacocks three times a week. Bird feces were collected for 3 months, with total genomic DNA being extracted and subjected to long-read 16S and 25S-28S sequencing. Genes coding for an iron permease (FTR1), iron receptors (FOB1 and FOB2), ADP-ribosylation factors (ARFs) (ARF2 and ARF6), and a GTPase (CDC42) were identified in this M. circinelloides genome. Also, this fungus was positive only for lecithinase activity. The field trial revealed a fecal microbiome dominated by Firmicutes and Proteobacteria in laying hens, and Firmicutes and Bacteroidetes in peacocks, whereas the fecal mycobiome of both bird species was mainly composed of Ascomycetes and Basidiomycetes fungi. Bacterial and fungal alpha-diversities did not differ between sampling time points after M. circinelloides administrations (P = 0.62 and P = 0.15, respectively). Although findings from this research suggest the lack of virulence of this M. circinelloides parasiticide isolate, more complementary in vitro and in vivo research is needed to conclude about the safety of its administration to birds, aiming at controlling their GI parasites.IMPORTANCEA previous study revealed that the native Mucor circinelloides isolate (FMV-FR1) can develop parasiticide activity toward coccidia oocysts, one of the most pathogenic GI parasites in birds. However, ensuring its safety for birds is of utmost importance, namely by studying its virulence profile and potential effect on commensal gut microbes. This initial study revealed that although this M. circinelloides isolate had genes coding for four types of virulence factors-iron permease, iron receptors, ADP-ribosylation factors, and GTPase-and only expressed phenotypically the enzyme lecithinase, the administration of its spores to laying hens and peacocks did not interfere with the abundances and diversities of their gut commensal bacteria and fungi. Although overall results suggest the lack of virulence of this M. circinelloides isolate, more complementary research is needed to conclude about the safety of its administration to birds in the scope of parasite biocontrol programs.

Molecular Phylogenetic and Comparative Genomic Analysis of Pleurocordyceps fusiformispora sp. nov. and Perennicordyceps elaphomyceticola in the Family Polycephalomycetaceae

Several Pleurocordyceps species have been reported as hyperparasitic fungi. A new species, Pleurocordyceps fusiformispora, and a known species, Perennicordyceps elaphomyceticola, are described here based on morphology and phylogenetic evidence from six genes (ITS, SSU, LSU, TET1-α, RPB1, and RPB2). Pl. fusiformispora differed from the other Pleurocordyceps species by producing flaky colonies, ovoid or elliptic α-conidia, and fusiform or long fusiform β-conidia. Both full genomes of Pe. elaphomyceticola and Pl. fusiformispora were sequenced, annotated, and compared. The antiSMASH and local BLAST analyses revealed significant differences in the number and types of putative secondary metabolite biosynthetic gene clusters, i.e., NPPS, PKS, and hybrid PKS-NRPS domains, between the two species. In addition, the putative BGCs of six compounds, namely ε-poly lysine, 4-epi-15-epi-brefeldin A, Monorden D/monocillin IV/monocillin VII/pochonin M/monocillin V/monocillin II, Tolypyridone, Piperazine, and Triticone DABFC, were excavated in the present study. This study motivates the use of heterologous expression and gene knockout methods to discover novel biologically active SMs from Polycephalomycetaceae.

Deciphering the Genomic Landscape and Virulence Mechanisms of the Wheat Powdery Mildew Pathogen Blumeria graminis f. sp. tritici Wtn1: Insights from Integrated Genome Assembly and Conidial Transcriptomics

A high-quality genome sequence from an Indian isolate of Blumeria graminis f. sp. tritici Wtn1, a persistent threat in wheat farming, was obtained using a hybrid method. The assembly of over 9.24 million DNA-sequence reads resulted in 93 contigs, totaling a 140.61 Mb genome size, potentially encoding 8480 genes. Notably, more than 73.80% of the genome, spanning approximately 102.14 Mb, comprises retro-elements, LTR elements, and P elements, influencing evolution and adaptation significantly. The phylogenomic analysis placed B. graminis f. sp. tritici Wtn1 in a distinct monocot-infecting clade. A total of 583 tRNA anticodon sequences were identified from the whole genome of the native virulent strain B. graminis f. sp. tritici, which comprises distinct genome features with high counts of tRNA anticodons for leucine (70), cysteine (61), alanine (58), and arginine (45), with only two stop codons (Opal and Ochre) present and the absence of the Amber stop codon. Comparative InterProScan analysis unveiled "shared and unique" proteins in B. graminis f. sp. tritici Wtn1. Identified were 7707 protein-encoding genes, annotated to different categories such as 805 effectors, 156 CAZymes, 6102 orthologous proteins, and 3180 distinct protein families (PFAMs). Among the effectors, genes like Avra10, Avrk1, Bcg-7, BEC1005, CSEP0105, CSEP0162, BEC1016, BEC1040, and HopI1 closely linked to pathogenesis and virulence were recognized. Transcriptome analysis highlighted abundant proteins associated with RNA processing and modification, post-translational modification, protein turnover, chaperones, and signal transduction. Examining the Environmental Information Processing Pathways in B. graminis f. sp. tritici Wtn1 revealed 393 genes across 33 signal transduction pathways. The key pathways included yeast MAPK signaling (53 genes), mTOR signaling (38 genes), PI3K-Akt signaling (23 genes), and AMPK signaling (21 genes). Additionally, pathways like FoxO, Phosphatidylinositol, the two-component system, and Ras signaling showed significant gene representation, each with 15-16 genes, key SNPs, and Indels in specific chromosomes highlighting their relevance to environmental responses and pathotype evolution. The SNP and InDel analysis resulted in about 3.56 million variants, including 3.45 million SNPs, 5050 insertions, and 5651 deletions within the whole genome of B. graminis f. sp. tritici Wtn1. These comprehensive genome and transcriptome datasets serve as crucial resources for understanding the pathogenicity, virulence effectors, retro-elements, and evolutionary origins of B. graminis f. sp. tritici Wtn1, aiding in developing robust strategies for the effective management of wheat powdery mildew.

Genomic insights into the evolution and adaptation of secondary metabolite gene clusters in fungicolous species Cladobotryum mycophilum ATHUM6906

Mycophilic or fungicolous fungi can be found wherever fungi exist since they are able to colonize other fungi, which occupy a diverse range of habitats. Some fungicolous species cause important diseases on Basidiomycetes, and thus, they are the main reason for the destruction of mushroom cultivations. Nonetheless, despite their ecological significance, their genomic data remain limited. Cladobotryum mycophilum is one of the most aggressive species of the genus, destroying the economically important Agaricus bisporus cultivations. The 40.7 Mb whole genome of the Greek isolate ATHUM6906 is assembled in 16 fragments, including the mitochondrial genome and 2 small circular mitochondrial plasmids, in this study. This genome includes a comprehensive set of 12,282 protein coding, 56 rRNA, and 273 tRNA genes. Transposable elements, CAZymes, and pathogenicity related genes were also examined. The genome of C. mycophilum contained a diverse arsenal of genes involved in secondary metabolism, forming 106 biosynthetic gene clusters, which renders this genome as one of the most BGC abundant among fungicolous species. Comparative analyses were performed for genomes of species of the family Hypocreaceae. Some BGCs identified in C. mycophilum genome exhibited similarities to clusters found in the family Hypocreaceae, suggesting vertical heritage. In contrast, certain BGCs showed a scattered distribution among Hypocreaceae species or were solely found in Cladobotryum genomes. This work provides evidence of extensive BGC losses, horizontal gene transfer events, and formation of novel BGCs during evolution, potentially driven by neutral or even positive selection pressures. These events may increase Cladobotryum fitness under various environmental conditions and potentially during host-fungus interaction.

Metabolites of pathogenic microorganisms database (MPMdb) and its seed metabolite applications

Seed metabolites are the combination of essential compounds required by an organism across various potential environmental conditions. The seed metabolites screening framework based on the network topology approach can capture important biological information of species. This study aims to identify comprehensively the relationship between seed metabolites and pathogenic bacteria. A large-scale data set was compiled, describing the seed metabolite sets and metabolite sets of 124,192 pathogenic strains from 34 genera, by constructing genome-scale metabolic models. The enrichment analysis method was used to screen the specific seed metabolites of each species/genus of pathogenic bacteria. The metabolites of pathogenic microorganisms database (MPMdb) (http://qyzhanglab.hzau.edu.cn/MPMdb/) was established for browsing, searching, predicting, or downloading metabolites and seed metabolites of pathogenic microorganisms. Based on the MPMdb, taxonomic and phylogenetic analyses of pathogenic bacteria were performed according to the function of seed metabolites and metabolites. The results showed that the seed metabolites could be used as a feature for microorganism chemotaxonomy, and they could mirror the phylogeny of pathogenic bacteria. In addition, our screened specific seed metabolites of pathogenic bacteria can be used not only for further tapping the nutritional resources and identifying auxotrophies of pathogenic bacteria but also for designing targeted bactericidal compounds by combining with existing antimicrobial agents.IMPORTANCEMetabolites serve as key communication links between pathogenic microorganisms and hosts, with seed metabolites being crucial for microbial growth, reproduction, external communication, and host infection. However, the large-scale screening of metabolites and the identification of seed metabolites have always been the main technical bottleneck due to the low throughput and costly analysis. Genome-scale metabolic models have become a recognized research paradigm to investigate the metabolic characteristics of species. The developed metabolites of pathogenic microorganisms database in this study is committed to systematically predicting and identifying the metabolites and seed metabolites of pathogenic microorganisms, which could provide a powerful resource platform for pathogenic bacteria research.

Draft genome sequence of 'Candidatus Phytoplasma asteris,' strain SW86 associated with sandal spike disease (SSD)

The sandal spike disease (SSD), related to 'Ca. Phytoplasma asteris' (Aster Yellows group), poses a significant threat to Indian sandalwood (Santalum album L.), making it the second most expensive wood globally due to declining population density. The epidemiology of SSD and the nature of the pathogen remain poorly understood. The SW86 isolate, collected from the Marayoor Sandalwood Reserve, was chosen for genome sequencing subsequent to confirming its titer and enriching phytoplasma DNA. Genome sequencing, utilizing Illumina and Oxford Nanopore Technology platforms, enabled a targeted hybrid metagenomic assembly resulting in 20 scaffolds totaling 554,025 bp, housing 436 protein-coding genes, 27 tRNA, and 1 rRNA operon. The genome analysis highlighted specific gene distributions, emphasizing translation, ribosomal structure, and biogenesis, with 352 genes assigned to 18 functional categories. Additionally, 322 proteins received functional assignments in the KEGG database, emphasizing 'Genetic Information Processing' and 'Environmental Information Processing'. Key potential pathogenicity factors, including signal peptide proteins and virulence proteins, were identified. Noteworthy findings include homologs of effectors genes like SAP11 and SAP05 and pathogenesis-related proteins, such as hemolysin III and SodA genes, in the SW86 genome. The duplicated cation-transporting P-type ATPase in the SW86 genome suggests a role in enhancing adaptability and contributing to the severity of SSD symptoms. This genome analysis provides crucial insights into the genomic features and potential virulence factors of 'Ca. Phytoplasma asteris' strain SW86, advancing our understanding of pathogenicity mechanisms and offering avenues for future disease management strategies in Indian sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03952-5.

Genomic and Proteomic Analyses of Extracellular Products Reveal Major Virulence Factors Likely Accounting for Differences in Pathogenicity to Bivalves between Vibrio mediterranei Strains

Vibrio mediterranei, a bacterial pathogen of bivalves, has exhibited strain-dependent virulence. The mechanisms behind the variations in bivalve pathogenicity between V. mediterranei strains have remained unclear. However, a preliminary analysis of the extracellular product (ECP) proteomes has revealed differences in protein compositions between low- and high-virulence strains; in addition to 1265 shared proteins, 127 proteins have been identified to be specific to one low-virulence strain and 95 proteins to be specific to two high-virulence strains. We further studied the ECP proteins of the three V. mediterranei strains from functional perspectives using integrated genomics and proteomics approaches. The results showed that lipid metabolism, transporter activity and membrane transporter pathways were more enriched in the ECPs of the two high-virulence strains than in those of the low-virulence strain. Additionally, 73 of the 95 high-virulence strain-specific proteins were found to have coding genes in the genome but were not expressed in the low-virulence strain. Moreover, comparisons with known virulence factors in the Virulence Factor Database (VFDB) and the Pathogen-Host Interactions Database (PHI-base) allowed us to predict more than 10 virulence factors in the categories of antimicrobial activity/competitive advantage, the effector delivery system and immune modulation, and the high-virulence strain-specific ECP proteins consisted of a greater percentage of known virulence factors than the low-virulence strain. Particularly, two virulence factors, MtrC and KatG, were identified in the ECPs of the two high-virulence strains but not in those of the low-virulence strain. Most coding genes of the ECP proteins including known virulence factors were identified on chromosome 1 of V. mediterranei. Our findings indicate that variations in virulence factor composition in the bacterial ECPs may partially account for the differences in the bivalve pathogenicity between V. mediterranei strains.

Genomic and Comparative Transcriptomic Analyses Reveal Key Genes Associated with the Biosynthesis Regulation of Okaramine B in Penicillium daleae NBP-49626

Restricted production of fungal secondary metabolites hinders the ability to conduct comprehensive research and development of novel biopesticides. Okaramine B from Penicillium demonstrates remarkable insecticidal efficacy; however, its biosynthetic yield is low, and its regulatory mechanism remains unknown. The present study found that the yield difference was influenced by fermentation modes in okaramine-producing strains and performed genomic and comparative transcriptome analysis of P. daleae strain NBP-49626, which exhibits significant features. The NBP-49626 genome is 37.4 Mb, and it encodes 10,131 protein-encoding genes. Up to 5097 differentially expressed genes (DEGs) were identified during the submerged and semi-solid fermentation processes. The oka gene cluster, lacking regulatory and transport genes, displayed distinct transcriptional patterns in response to the fermentation modes and yield of Okaramine B. Although transcription trends of most known global regulatory genes are inconsistent with those of oka, this study identified five potential regulatory genes, including two novel Zn(II)2Cys6 transcription factors, Reg2 and Reg19. A significant correlation was also observed between tryptophan metabolism and Okaramine B yields. In addition, several transporter genes were identified as DEGs. These results were confirmed using real-time quantitative PCR. This study provides comprehensive information regarding the regulatory mechanism of Okaramine B biosynthesis in Penicillium and is critical to the further yield improvement for the development of insecticides.

PBITV3 : A robust and comprehensive tool for screening pathogenic proteomes for drug targets and prioritizing vaccine candidates

Rise of life-threatening superbugs, pandemics and epidemics warrants the need for cost-effective and novel pharmacological interventions. Availability of publicly available proteomes of pathogens supports development of high-throughput discovery platforms to prioritize potential drug-targets and develop testable hypothesis for pharmacological screening. The pipeline builder for identification of target (PBIT) was developed in 2016 and updated in 2021, with the purpose of accelerating the search for drug-targets by integration of methods like comparative and subtractive genomics, essentiality/virulence and druggability analysis. Since then, it has been used for identification of drugs and vaccine targets, safety profiling of multiepitope vaccines and mRNA vaccine construction against a broad-spectrum of pathogens. This tool has now been updated with functionalities related to systems biology and immuno-informatics and validated by analyzing 48 putative antigens of Mycobacterium tuberculosis documented in literature. PBITv3 available as both online and offline tools will enhance drug discovery against emerging drug-resistant infectious agents. PBITv3 can be freely accessed at http://pbit.bicnirrh.res.in/.

Strategies for tailoring functional microbial synthetic communities

Natural ecosystems harbor a huge reservoir of taxonomically diverse microbes that are important for plant growth and health. The vast diversity of soil microorganisms and their complex interactions make it challenging to pinpoint the main players important for the life support functions microbes can provide to plants, including enhanced tolerance to (a)biotic stress factors. Designing simplified microbial synthetic communities (SynComs) helps reduce this complexity to unravel the molecular and chemical basis and interplay of specific microbiome functions. While SynComs have been successfully employed to dissect microbial interactions or reproduce microbiome-associated phenotypes, the assembly and reconstitution of these communities have often been based on generic abundance patterns or taxonomic identities and co-occurrences but have only rarely been informed by functional traits. Here, we review recent studies on designing functional SynComs to reveal common principles and discuss multidimensional approaches for community design. We propose a strategy for tailoring the design of functional SynComs based on integration of high-throughput experimental assays with microbial strains and computational genomic analyses of their functional capabilities.

Whole Genome Sequencing and Comparative Analysis of the First Ehrlichia canis Isolate in China

Ehrlichia canis, a prominent tick-borne pathogen causing canine monocytic ehrlichiosis (CME), is one of the six recognized Ehrlichia species worldwide. Despite its widespread presence in ticks and host dogs in China, comprehensive genomic information about this pathogen remains limited. This study focuses on an in-depth analysis of E. canis YZ-1, isolated and cultured from an infected dog in China. The complete genome of E. canis YZ-1 was sequenced (1,314,789 bp, 1022 genes, 29% GC content, and 73% coding bases), systematically characterizing its genomic elements and functions. Comparative analysis with representative genomes of Ehrlichia species, including E. canis strain Jake, E. chaffeensis, Ehrlichia spp., E. muris, E. ruminantium, and E. minasensis, revealed conserved genes, indicating potential evolutionary connections with E. ruminantium. The observed reduction in virulence-associated genes, coupled with a type IV secretion system (T4SS), suggests an intricate balance between pathogenicity and host adaptation. The close relationship with E. canis Jake and E. chaffeensis, alongside nuanced genomic variations with E. ruminantium and E. mineirensis, underscores the need to explore emerging strains and advancements in sequencing technologies continuously. This genetic insight opens avenues for innovative medications, studies on probiotic resistance, development of new detection markers, and progress in vaccine development for ehrlichiosis. Further investigations into the functional significance of identified genes and their role in host-pathogen interactions will contribute to a more holistic comprehension of Ehrlichia's biology and its implications for pathogenicity and transmission.

Ensembl 2024

Ensembl (https://www.ensembl.org) is a freely available genomic resource that has produced high-quality annotations, tools, and services for vertebrates and model organisms for more than two decades. In recent years, there has been a dramatic shift in the genomic landscape, with a large increase in the number and phylogenetic breadth of high-quality reference genomes, alongside major advances in the pan-genome representations of higher species. In order to support these efforts and accelerate downstream research, Ensembl continues to focus on scaling for the rapid annotation of new genome assemblies, developing new methods for comparative analysis, and expanding the depth and quality of our genome annotations. This year we have continued our expansion to support global biodiversity research, doubling the number of annotated genomes we support on our Rapid Release site to over 1700, driven by our close collaboration with biodiversity projects such as Darwin Tree of Life. We have also strengthened support for key agricultural species, including the first regulatory builds for farmed animals, and have updated key tools and resources that support the global scientific community, notably the Ensembl Variant Effect Predictor. Ensembl data, software, and tools are freely available.

Reference Genome Resource for the Citrus Pathogen Phytophthora citrophthora

Phytophthora citrophthora is an oomycete pathogen that infects citrus. Its occurrence in citrus-growing regions worldwide is considered a major contributor to crop losses. This study presents a high-quality genome resource for P. citrophthora, which was generated using PacBio HiFi long-read high-throughput sequencing technology. We successfully assembled a 48.5 Mb genome containing 16,409 protein-coding genes from high-quality reads. This marks the first complete genome assembly of P. citrophthora, providing a valuable resource to enhance the understanding of pathogenic behaviour and fungicide sensitivity of this destructive citrus pathogen.

Genome-resolved metatranscriptomics reveals conserved root colonization determinants in a synthetic microbiota

The identification of processes activated by specific microbes during microbiota colonization of plant roots has been hampered by technical constraints in metatranscriptomics. These include lack of reference genomes, high representation of host or microbial rRNA sequences in datasets, or difficulty to experimentally validate gene functions. Here, we recolonized germ-free Arabidopsis thaliana with a synthetic, yet representative root microbiota comprising 106 genome-sequenced bacterial and fungal isolates. We used multi-kingdom rRNA depletion, deep RNA-sequencing and read mapping against reference microbial genomes to analyse the in planta metatranscriptome of abundant colonizers. We identified over 3,000 microbial genes that were differentially regulated at the soil-root interface. Translation and energy production processes were consistently activated in planta, and their induction correlated with bacterial strains' abundance in roots. Finally, we used targeted mutagenesis to show that several genes consistently induced by multiple bacteria are required for root colonization in one of the abundant bacterial strains (a genetically tractable Rhodanobacter). Our results indicate that microbiota members activate strain-specific processes but also common gene sets to colonize plant roots.

Virulence factors of the genus Fusarium with targets in plants

Fusarium spp. comprise various species of filamentous fungi that cause severe diseases in plant crops of both agricultural and forestry interest. These plant pathogens produce a wide range of molecules with diverse chemical structures and biological activities. Genetic functional analyses of some of these compounds have shown their role as virulence factors (VF). However, their mode of action and contributions to the infection process for many of these molecules are still unknown. This review aims to analyze the state of the art in Fusarium VF, emphasizing their biological targets on the plant hosts. It also addresses the current experimental approaches to improve our understanding of their role in virulence and suggests relevant research questions that remain to be answered with a greater focus on species of agroeconomic importance. In this review, a total of 37 confirmed VF are described, including 22 proteinaceous and 15 non-proteinaceous molecules, mainly from Fusarium oxysporum and Fusarium graminearum and, to a lesser extent, in Fusarium verticillioides and Fusarium solani.

Insights into the genomic architecture of a newly discovered endophytic Fusarium species belonging to the Fusarium concolor complex from India

In this study, a new species Fusarium indicum belonging to the Fusarium concolor species complex is established to accommodate an endophytic fungus isolated from Bambusa sp. and collected from Himachal Pradesh. The identity of this isolate was confirmed based on the asexual morphs, its cultural characteristics, and phylogenetic analyses. This isolate revealed out to be distinct by showing less similarity with described species in the genus Fusarium based on molecular sequence data, approximately 93.9% similarity based on translation elongation factor 1-alpha, and 94.2% similarity based on RNA polymerase II subunit. Furthermore, to increase knowledge about this novel species, whole-genome sequencing was carried out. The results displayed that Fusarium indicum NFCCI 5145 possesses a 40.2 Mb genome and 48.39% of GC content. Approximately 12,963 functional protein-coding genes were carefully predicted and annotated using different BLAST databases, such as Uniprot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), Pathogen Host Interactions (PHI), Clusters of Orthologous Groups (COG), and Carbohydrate-Active enzymes (CAZy). The orthologous proteins were identified using OrthoFinder and used for the phylogenetic analysis. ANIb confirmed that the isolate is closely related to the F. concolor species complex. It is known that Fusarium strains can produce a wide range of bioactive secondary metabolites. Therefore, in-depth mining for biosynthetic gene clusters for secondary metabolite biosynthesis of Fusarium indicum NFCCI 5145 was investigated using Antibiotics and Secondary Metabolites Analysis Shell (AntiSMASH) annotation. AntiSMASH results displayed that this isolate possesses 45 secondary metabolites of biosynthetic gene clusters (BGCs). These findings significantly improved our understanding of the strain Fusarium indicum NFCCI 5145 and its possible applications in different sectors including industry for the secondary metabolites and enzymes it can produce.

Genomic Comparison of Two Species of Samsoniella with Other Genera in the Family Cordycipitaceae

Whole genomes of Samsoniella hepiali ICMM 82-2 and S. yunnanensis YFCC 1527 were sequenced and annotated, as well as compared with whole genome sequences of other species in the family Cordycipitaceae. S. hepiali ICMM 82-2, S. hepiali FENG and S. yunnanensis YFCC 1527 had 54, 57 and 58 putative secondary metabolite biosynthetic gene clusters, respectively. S. hepiali had one unique domain and S. yunnanensis YFCC 1527 six. Both S. hepiali and S. yunnanensis YFCC 1527 had curvupallide-B, fumosorinone and fujikurin putative biosynthetic gene clusters. C. javanica had biosynthetic gene clusters for fumonisin. The 14 genomes had common domains, namely A-P-C-P-C and KS-AT-DH-ER-KR-ACP. The A-P-C-P-C domain may be involved in the biosynthesis of dimethylcoprogen. The maximum likelihood and the Bayesian inference trees of KS-AT-DH-ER-KR-ACP were highly consistent with the multigene phylogenetic tree for the 13 species of Cordycipitaceae. This study facilitates the discovery of novel biologically active SMs from Cordycipitaceae using heterologous expression and gene knockdown methods.

Application of high-throughput sequencing technologies and analytical tools for pathogen detection in urban water systems: Progress and future perspectives

The ubiquitous presence of pathogenic microorganisms, such as viruses, bacteria, fungi, and protozoa, in urban water systems poses a significant risk to public health. The emergence of infectious waterborne diseases mediated by urban water systems has become one of the leading global causes of mortality. However, the detection and monitoring of these pathogenic microorganisms have been limited by the complexity and diversity in the environmental samples. Conventional methods were restricted by long assay time, high benchmarks of identification, and narrow application sceneries. Novel technologies, such as high-throughput sequencing technologies, enable potentially full-spectrum detection of trace pathogenic microorganisms in complex environmental matrices. This review discusses the current state of high-throughput sequencing technologies for identifying pathogenic microorganisms in urban water systems with a concise summary. Furthermore, future perspectives in pathogen research emphasize the need for detection methods with high accuracy and sensitivity, the establishment of precise detection standards and procedures, and the significance of bioinformatics software and platforms. We have compiled a list of pathogens analysis software/platforms/databases that boast robust engines and high accuracy for preference. We highlight the significance of analyses by combining targeted and non-targeted sequencing technologies, short and long reads technologies, sequencing technologies, and bioinformatic tools in pursuing upgraded biosafety in urban water systems.

LeishMANIAdb: a comparative resource for Leishmania proteins

Leishmaniasis is a detrimental disease causing serious changes in quality of life and some forms can lead to death. The disease is spread by the parasite Leishmania transmitted by sandfly vectors and their primary hosts are vertebrates including humans. The pathogen penetrates host cells and secretes proteins (the secretome) to repurpose cells for pathogen growth and to alter cell signaling via host-pathogen protein-protein interactions). Here, we present LeishMANIAdb, a database specifically designed to investigate how Leishmania virulence factors may interfere with host proteins. Since the secretomes of different Leishmania species are only partially characterized, we collated various experimental evidence and used computational predictions to identify Leishmania secreted proteins to generate a user-friendly unified web resource allowing users to access all information available on experimental and predicted secretomes. In addition, we manually annotated host-pathogen interactions of 211 proteins and the localization/function of 3764 transmembrane (TM) proteins of different Leishmania species. We also enriched all proteins with automatic structural and functional predictions that can provide new insights in the molecular mechanisms of infection. Our database may provide novel insights into Leishmania host-pathogen interactions and help to identify new therapeutic targets for this neglected disease. Database URL  https://leishmaniadb.ttk.hu/.

Host Recognition and Specific Infection of Endomelanconiopsis endophytica during Early Infection

Coevolution between the pathogen and host plant drives pathogenic effector diversity. However, the molecular mechanism behind host-specific pathogenesis remains to be explored. Here, we present a 43 Mb whole-genome sequence of Endomelanconiopsis endophytica strain LS29, a host-specific pathogen of the common subtropical tree Castanopsis fissa. We described its genome annotations and identified its effector candidates. By performing temporal transcriptome sequencing of E. endophytica on C. fissa during early infection, we found that E. endophytica repressed other microbes in order to attack the tissue of the host by producing antibiotics earlier than 24 h post-inoculation (hpi). Simultaneously, a variety of effectors were secreted to recognize the host plant, but most of them showed a significantly opposing expression regulation trend after 24 hpi, indicating that 24 hpi represents a key time point between host recognition and specific infection. Furthermore, a comparison of isoenzymes showed that only a few effectors were identified as specific effectors, which were involved in hydrolyzing the compounds of the plant cell wall and releasing fatty acids during the early infection of C. fissa. Our results determined host recognition timing and identified a specific catalog of effectors, which are crucial for revealing the molecular mechanism of host-specific pathogenesis.

Telluribacter roseus sp. nov., Isolated from the Kumtag Desert Soil

A pink-pigmented bacterium, designated as strain SYSU D00476T, was isolated from sandy soil collected from the Kumtag Desert in China. Colonies were opaque, smooth and of a slight convexity with a clearly defined border. Cells were rod-shaped, Gram-stain-negative, catalase- and oxidase-positive. Growth occurred at 4-45 ℃ (optimum at 28-30 ℃), pH 6.0-8.0 (optimum at 7.0), and with 0-3.0% NaCl (w/v, optimum at 0-2.0%). Major fatty acids (> 10%) were C16:0, summed feature 3 (C16:1 ω7c and/or C16:1 ω6c), iso-C17:0 3-OH and iso-C15:0. Polar lipids comprised of three unidentified polar aminolipids (ALs), two unidentified aminophosphoglycolipids (APLs), one unidentified glycolipid (GL) and three unidentified phospholipids (PLs). The predominant respiratory quinone was MK-7. The genomic DNA G + C content was 50.5%. The low digital DNA-DNA hybridization (dDDH, 27.4%) and average nucleotide identity (ANI, 85%) values between strain SYSU D00476T and Telluribacter humicola KCTC 42819T indicated that SYSU D00476T represent a distinct species. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain SYSU D00476T belonged to the genus Telluribacter, showing 97.5% similarity with T. humicola KCTC 42819T. All these data support that strain SYSU D00476T represent a novel species of the genus Telluribacter within the family Spirosomataceae, named as Telluribacter roseus sp. nov. The type strain is SYSU D00476T (= KCTC 82285T = CGMCC 1.18647T = MCCC 1K04983T).

Comprehensive genomic analysis of Bacillus paralicheniformis strain BP9, pan-genomic and genetic basis of biocontrol mechanism

Many Bacillus species are essential antibacterial agents, but their antibiosis potential still needs to be elucidated to its full extent. Here, we isolated a soil bacterium, BP9, which has significant antibiosis activity against fungal and bacterial pathogens. BP9 improved the growth of wheat seedlings via active colonization and demonstrated effective biofilm and swarming activity. BP9 sequenced genome contains 4282 genes with a mean G-C content of 45.94% of the whole genome. A single copy concatenated 802 core genes of 28 genomes, and their calculated average nucleotide identity (ANI) discriminated the strain BP9 from Bacillus licheniformis and classified it as Bacillus paralicheniformis. Furthermore, a comparative pan-genome analysis of 40 B. paralicheniformis strains suggested that the genetic repertoire of BP9 belongs to open-type genome species. A comparative analysis of a pan-genome dataset using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Cluster of Orthologous Gene groups (COG) revealed the diversity of secondary metabolic pathways, where BP9 distinguishes itself by exhibiting a greater prevalence of loci associated with the metabolism and transportation of organic and inorganic substances, carbohydrate and amino acid for effective inhabitation in diverse environments. The primary secondary metabolites and their genes involved in synthesizing bacillibactin, fencing, bacitracin, and lantibiotics were identified as acquired through a recent Horizontal gene transfer (HGT) event, which contributes to a significant part of the strain`s antimicrobial potential. Finally, we report some genes essential for plant-host interaction identified in BP9, which reduce spore germination and virulence of multiple fungal and bacterial species. The effective colonization, diverse predicted metabolic pathways and secondary metabolites (antibiotics) suggest testing the suitability of strain BP9 as a potential bio-preparation in agricultural fields.