PHI-base: a new database for pathogen host interactions

Comprehensive genome analysis of two Cytospora (Cytosporaceae, Diaporthales) species associated with canker disease of spruce: C.piceae and C.piceicola sp. nov

Cytospora canker (CC) is among the most important diseases in conifer trees (Picea spp., mainly). This disease poses a significant risk factor for forest health, potentially leading to economic losses for wood producers. To provide a genomic basis of the CC pathogenesis, the genomes of two Cytospora species associated with the disease were sequenced and further analyzed within a set of Diaporthales species. The first species was identified as C.piceae. The second was described as C.piceicola sp. nov. based on morphological characteristics and multi-gene phylogenetic analysis. The novel species is sister to other Cytospora species isolated from conifers. Here, we report 39.7 and 43.8 Mb highly contiguous genome assemblies of C.piceae EI-19(A) and C.piceicola EI-20, respectively, obtained using Illumina sequencing technology. Despite notably different genome sizes, these species share the main genome characteristics, such as predicted gene number (10,862 and 10,742) and assembly completeness (97.6% and 98.1%). A wide range of genes encoding carbohydrate-active enzymes, secondary metabolite biosynthesis clusters, and secreted effectors were found. Multiple experimentally validated virulence genes were also identified in the studied species. The defined arsenals of enzymes and effectors generally relate to the hemibiotrophic lifestyle with a capability to switch to biotrophy. The obtained evidence also supports that C.piceae EI-19(A) and C.piceicola EI-20 can cause severe canker disease symptoms in Picea spp. specifically. It was additionally observed that the strains of C.piceae may have different pathogenicity and virulence characteristics based on the analyses of predicted secondary metabolite complements, effectomes, and virulence-related genes. Phylogenomic analysis and timetree estimations indicated that divergence of the studied species may have occurred relatively late, 11-10 million years ago. Compared to other members of Diaporthales, C.piceae EI-19(A) and C.piceicola EI-20 implied a moderate rate of gene contraction, but the latter experienced significant gene loss that can additionally support host specificity attributed to these species. But uncovered gene contraction events may point out potential lifestyle differentiation and host shift of the studied species. It was revealed that EI-19(A) and C.piceicola EI-20 carry distinct secretomes and effectomes among Diaporthales species. This feature can indicate a species lifestyle and pathogenicity potential. These findings highlight potential targets for identification and/or detection of pathogenic Cytospora in conifers. The introduced draft genome sequences of C.piceae and C.piceicola can be employed as tools to understand basic genetics and pathogenicity mechanisms of fungal species causing canker disease in woody plants. The identified pathogenicity and virulence-related genes would serve as potential candidates for host-induced gene silencing aimed at making plant hosts more resistant to pathogenic species. Furthermore, the comparative genomics component of the study will facilitate the functional analysis of the genes of unknown function in all fungal pathogens.

Genome sequencing of a novel Verticillium dahliae strain (huangweibingjun)

Verticillium dahliae is a soilborne pathogenic fungus that causes vascular discoloration and wilting in a broad spectrum of plant hosts, affecting about 400 species, such as cotton, potatoes, watermelon, cucumber, spinach, etc. In 2021, V. dahliae was estimated to cause about 15-20% reduction in cotton in China. Here, we report the genome sequencing of a novel strain named huangweibingjun, isolated from diseased cotton roots in the Henan province of China. The huangweibingjun genome consists of a total size of 35.84 Mb, GC content of 59.835%, and harbors six chromosomes (scaffold7561, scaffold7329, scaffold7795, scaffold5491, scaffold5473, and scaffold4511). The genome architecture showed a high diversity of cell wall-degrading secretory proteins that might influence the pathogenicity of the fungal strain. Moreover, preliminary metabolic pathway prediction showed that this novel strain synthesizes polyketide, terpenoids, shikimic acid-derived compounds and could also be aflatoxigenic. Consistent with other pathogenic microbes, the huangweibingjun genome comprises several virulent-associated genes. This genome assembly lays the foundation for further investigation of the pathogenicity of huangweibingjun.

Unveiling the complete genome sequence of Paenibacillus taichungensis: genomic features and biocontrol potential

OBJECTIVES

The genus Paenibacillus encompasses a diverse group of Gram-positive bacteria with significant biotechnological potential. However, the research data and application cases of Paenibacillus taichungensis were still poorly understood. In this study, we isolated a P. taichungensis strain BB507, which demonstrated antibacterial effect on Ralstonia solanacearum species complex, and provided data and analysis of its complete genome.

DATA DESCRIPTION

Strain BB507 was isolated from a pine rhizosphere in Meizhou city, Guangdong province of China, and showed antibacterial activity against Ralstonia solanacearum species complex. Complete genome was sequenced using Illumina NovaSeq (second-generation) and Oxford Nanopore (third-generation) platforms. The genome of BB507 comprised of a 6,974,531 bp circular chromosome and a 352,197 bp circular plasmid, encoding a total of 6,649 gene with an average gene length of 950 bp, 103 transfer RNAs, 2 sRNAs, and 36 rRNAs. Three candidate CRISPRs, 6 genomic islands and 14 prophages were predicted. The bacterial orthologous average nucleotide identity (OAT) and the type genome server (TYGS) analysis highlighted the strain BB507 was clustered into a subgroup with P. taichungensis. antiSMASH v7.0 predicted the presence of 10 secondary metabolite gene clusters in the genome. These findings will serve as a useful resource for further research in industrial and agricultural biotechnology.

Characterization and Genome Analysis of Fusarium oxysporum Provides Insights into the Pathogenic Mechanisms of the Pokkah Boeng Disease in China

Pokkah Boeng Disease (PBD) is a severe and devastating disease that causes significant damage and yield losses in China. The pathogenic fungus Fusarium oxysporum is responsible for the rapid onset of top rot symptoms in sugarcane. In this study, we selected a representative strain, BS2-6, to perform morphological observations of colonies and determine pathogenicity. We examined the effects of BS2-6 infestation on the ultrastructure of sugarcane leaves. Moreover, we sequenced the whole genome of BS2-6 and examined the effects of various nitrogen sources and chemical reagents on its growth and pathogenicity. Our results indicate that sugarcane leaves inoculated with BS2-6 quickly succumb to heart leaf and growing rot. Ultrastructural analysis revealed that the surface tissues of the diseased leaves were destroyed with mycelium, and conidia blocked leaf stomata, which ultimately led to the degradation of leaf tissues. Ammoniacal nitrogen significantly promoted mycelial growth, pigment secretion, and the expression of genes related to secondary metabolite synthesis, thereby accelerating the development of PBD. In addition, we found that carbendazim effectively inhibited the growth of BS2-6 at various concentrations. These findings provide important insights for the effective prevention and control of PBD during sugarcane production.

Genomic Analysis of Penicillium griseofulvum CF3 Reveals Potential for Plant Growth Promotion and Disease Resistance

Penicillium griseofulvum CF3 is a fungus isolated from healthy strawberry soil, with the potential to promote the growth of plants and enhance their resistance to diseases. However, the genome sequence of P. griseofulvum CF3 remains unclear. Therefore, we performed the whole-genome CCS sequencing of P. griseofulvum CF3 using the PacBio Sequel II platform. The assembled genome comprised 104 contigs, with a total length of 37,564,657 bp, encoding 13,252 protein-coding genes. Comprehensive functional annotation was performed using various BLAST databases, including the non-redundant (Nr) protein sequence database, Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), EuKaryotic Orthologous Groups (KOG), and the Carbohydrate-Active enZymes (CAZy) database, to identify and predict protein-coding genes, tRNAs, and rRNAs. The Antibiotics and Secondary Metabolites Analysis Shell (Antismash) analysis identified 50 biosynthetic gene clusters involved in secondary metabolite production within the P. griseofulvum CF3 genome. The whole-genome sequencing of P. griseofulvum CF3 helps us to understand its potential mechanisms in promoting plant growth and enhancing disease resistance, paving the way for the application of the CF3 strain in sustainable crop production.

Design and regulation of engineered bacteria for environmental release

Emerging products of biotechnology involve the release of living genetically modified microbes (GMMs) into the environment. However, regulatory challenges limit their use. So far, GMMs have mainly been tested in agriculture and environmental cleanup, with few approved for commercial purposes. Current government regulations do not sufficiently address modern genetic engineering and limit the potential of new applications, including living therapeutics, engineered living materials, self-healing infrastructure, anticorrosion coatings and consumer products. Here, based on 47 global studies on soil-released GMMs and laboratory microcosm experiments, we discuss the environmental behaviour of released bacteria and offer engineering strategies to help improve performance, control persistence and reduce risk. Furthermore, advanced technologies that improve GMM function and control, but lead to increases in regulatory scrutiny, are reviewed. Finally, we propose a new regulatory framework informed by recent data to maximize the benefits of GMMs and address risks.

How do various strategies for returning residues change microbiota modulation: potential implications for soil health

Introduction

Residue incorporation is a crucial aspect of anthropogenic land management practices in agricultural fields. However, the effects of various returning strategies on the soil microbiota, which play an essential vital role in maintaining soil health, remains largely unexplored.

Methods

In a study conducted, different residue management strategies were implemented, involving the application of chemical fertilizers and residues that had undergone chopping (SD), composting (SC), and pyrolysis (BC) processes, with conventional fertilization serving as the control (CK).

Results and discussion

Using metagenomic sequencing, the analysis revealed that while all residue returning strategies had minimal effects on the diversity (both α and β) of microbiota, they did significantly alter microbial functional genes related to carbon (C), nitrogen (N), phosphorus (P), and sulfur (S) cycling, as well as the presence of antibiotic resistance genes (ARGs) and pathogens. Specifically, chopped residues were found to enhance microbial genes associated with C, N, P, and S cycling, while composted residues primarily stimulated C and S cycling. Furthermore, all residue treatments resulted in a disruption of relationships among nutrient cycles, with varying degrees of impact observed across the different management strategies, with the sequence of impact being SD < SC < BC. Moreover, the residue additions resulted in the accumulation of ARGs, while only SC caused an increase in certain pathogens. Finally, through analyzing the correlation network among indices that exhibited active responses to residue additions, potential indicators for functional changes in response to residue additions were identified. This study further offered recommendations for future cropland management practices aimed at enhancing soil health through microbiomes.

Host-Pathogen Interaction Databases: Tools for Rapid Understanding of Microbial Pathogenesis

Understanding of microbial pathogenesis has greatly revolutionized after conventional culture-based techniques are replaced by molecular methods. This technological shift is generating huge host-pathogen interactions (HPIs) data. Moreover, computational predictions of biological interactions are also adding to HPI understanding. Recently, several dedicated databases are developed for exclusively cataloging HPIs. Present article covers about some available HPI databases, types, and evolution of this area, along with recent trends in the application of these databases for biological research. As per the recent understanding in microbial pathogenesis, HPIs are considered highly dynamic in nature with multiple outcomes, which goes beyond simple microbes-disease association. Therefore, careful cataloging of complete information about HPIs can open several avenues to understand microbial pathogenesis considering their multifaceted effects on host system. HPI databases are indispensable tools for understanding microbial pathogenesis, and this article provides comprehensive information about their uses in the field of microbial pathogenesis research.

Comparative secretome analysis unveils species-specific virulence factors in Elsinoe perseae, the causative agent of the scab disease of avocado (Persea americana)

The scab disease, caused by Elsinoe perseae, poses a significant risk to avocado (Persea americana) production in countries with warm and humid climates. Although the genome has been published, the precise virulence factors accountable for the pathogenicity of E. perseae have not yet been determined. The current study employed an in silico approach to identify and functionally characterize the secretory proteins of E. perseae. A total of 654 potential secretory proteins were identified, of which 190 were classified as carbohydrate-active enzymes (CAZymes), 49 as proteases, and 155 as potential effectors. A comparison to six other closely related species identified 40 species-specific putative effectors in E. perseae, indicating their specific involvement in the pathogenicity of E. perseae on avocado. The data presented in this study might be valuable for further research focused on understanding the molecular mechanisms that contribute to the pathogenicity of E. perseae on avocado.

Whole-Genome Sequencing of the Entomopathogenic Fungus Fusarium solani KMZW-1 and Its Efficacy Against Bactrocera dorsalis

Fusarium solani KMZW-1 is recognized for its potential as a biocontrol agent against agricultural and forestry pests, particularly due to its compatibility with integrated pest management (IPM) strategies. This study aimed to investigate the complete genome of F. solani KMZW-1 and assess its pathogenicity against Bactrocera dorsalis. Whole-genome sequencing revealed a genome size of 47,239,278 bp, comprising 27 contigs, with a GC content of 51.16% and fungus identified as F. solani KMZW-1. The genome completeness was assessed as 97.93% using BUSCO analysis, the DFVF sequence identifier was Fusarium 0G092560.1, and AntiSMASH analysis identified 35 gene clusters associated with secondary metabolite biosynthesis, providing insights into the genetic basis of its pathogenic mechanisms and biocontrol potential. Comparative genomic analysis found 269 unique genes for F. solani KMZW-1, and collinearity analysis exhibited a high degree of synteny with Fusarium solani-melongenae. The pathogenicity of F. solani KMZW-1 was assessed using concentrations ranging from 1 × 104 to 1 × 1011 conidia/mL. Higher concentrations (1 × 1010 to 1 × 1011 conidia/mL) resulted in significantly increased cumulative mortality rates of B. dorsalis adults compared to the control group. Notably, the pathogenicity was higher in male adults than in females. Probit analysis yielded LC50 (50% lethal concentration) values of 5.662 for female and 4.486 for male B. dorsalis adults. In summary, F. solani, KMZW-1 exhibits strong insecticidal activity against B. dorsalis and shows potential as a biocontrol agent with IPM strategies. These findings provide robust genomic evidence supporting the use of F. solani KMZW-1 in managing against B. dorsalis populations.

Complete genome sequences of two Pantoea stewartii strains ATCC 8199 from maize and PSCN1 from sugarcane

OBJECTIVES

The pathogen of Pantoea stewartii (Ps) is the causal agent of bacterial disease in corn and various graminaceous plants. Ps has two subspecies, Pantoea stewartii subsp. stewartia (Pss) and Pantoea stewartii subsp. indologenes (Psi). This study presents two complete genomes of Ps strains including ATCC 8199 isolated from maize and PSCN1 causing bacterial wilt in sugarcane. The two bacterial genomes information will be helpful for taxonomy analysis in this genus Pantoea at whole-genome levels and accurately discriminated the two subspecies of Pss and Psi.

DATA DESCRIPTION

The reference strain ATCC 8199 isolated from maize was purchased from Beijing Biobw Biotechnology Co., Ltd. (China) and the strain of PSCN1 was isolated from sugarcane cultivar YZ08-1095 in Zhanjiang, Guangdong province of China. Two complete genomes were sequenced using Illumina Hiseq (second-generation) and Oxford Nanopore (third-generation) platforms. The genome of the strain ATCC 8199 comprised of 4.78 Mb with an average GC content of 54.03%, along with five plasmids, encoding a total of 4,846 gene with an average gene length of 827 bp. The genome of PSCN1 comprised of 5.03 Mb with an average GC content of 53.78%, along with two plasmids, encoding a total of 4,725 gene with an average gene length of 913 bp. The bacterial pan-genome analysis highlighted the strain ATCC 8199 was clustered into a subgroup with a Pss strain CCUG 26,359 from USA, while the strain PSCN1 was clustered into another subgroup with a Ps strain NRRLB-133 from USA. These findings will serve as a useful resource for further analyses of the evolution of Ps strains and corresponding disease epidemiology worldwide.

Comparative genomics of Metarhizium brunneum strains V275 and ARSEF 4556: unraveling intraspecies diversity

Entomopathogenic fungi belonging to the Order Hypocreales are renowned for their ability to infect and kill insect hosts, while their endophytic mode of life and the beneficial rhizosphere effects on plant hosts have only been recently recognized. Understanding the molecular mechanisms underlying their different lifestyles could optimize their potential as both biocontrol and biofertilizer agents, as well as the wider appreciation of niche plasticity in fungal ecology. This study describes the comprehensive whole genome sequencing and analysis of one of the most effective entomopathogenic and endophytic EPF strains, Metarhizium brunneum V275 (commercially known as Lalguard Met52), achieved through Nanopore and Illumina reads. Comparative genomics for exploring intraspecies variability and analyses of key gene sets were conducted with a second effective EPF strain, M. brunneum ARSEF 4556. The search for strain- or species-specific genes was extended to M. brunneum strain ARSEF 3297 and other species of genus Metarhizium, to identify molecular mechanisms and putative key genome adaptations associated with mode of life differences. Genome size differed significantly, with M. brunneum V275 having the largest genome amongst M. brunneum strains sequenced to date. Genome analyses revealed an abundance of plant-degrading enzymes, plant colonization-associated genes, and intriguing intraspecies variations regarding their predicted secondary metabolic compounds and the number and localization of Transposable Elements. The potential significance of the differences found between closely related endophytic and entomopathogenic fungi, regarding plant growth-promoting and entomopathogenic abilities, are discussed, enhancing our understanding of their diverse functionalities and putative applications in agriculture and ecology.

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.

A Genomic Sequence Resource of Diaporthe mahothocarpus GZU-Y2 Causing Leaf Spot Blight in Camellia oleifera

Diaporthe mahothocarpus GZU-Y2, a new pathogen responsible for leaf spot blight disease, leads to significant damage and economic losses in some Camellia oleifera plantations. The current study annotated the genome of the D. mahothocarpus isolate GZU-Y2 to advance our knowledge of the pathogen and facilitate improved disease management of leaf spot blight. The initial Pacbio-Illumina hybrid draft genome for GZU-Y2 resulted in a high-quality assembly with 62 contigs, characterized by an N50 length of 7.07 Mb. The complete genome of isolate GZU-Y2 was 58.97 Mbp, with a GC content of 50.65%. Importantly, the assembly exhibits remarkable integrity, with 97.93% of complete BUSCO validating genome completeness. The prediction results showed that a total of 15,918 protein-coding genes were annotated using multiple bioinformatics databases. The genome assembly and annotation resource reported here will be useful for the further study of fungal infection mechanisms and pathogen-host interaction.

Potential Probiotic Properties and Complete Genome Analysis of Limosilactobacillus reuteri LRA7 from Dogs

This study aimed to isolate and screen canine-derived probiotics with excellent probiotic properties. Strain characterization was conducted using a combination of in vitro and in vivo probiotic characterization and safety assessments, as well as complete genome analysis. The results showed that Limosilactobacillus reuteri LRA7 exhibited excellent bacteriostatic and antioxidant activities. The survival rate at pH 2.5 was 79.98%, and the viable counts after exposure to gastrointestinal fluid and 0.5% bile salts were 7.77 log CFU/mL and 5.29 log CFU/mL, respectively. The bacterium also exhibited high hydrophobicity, self-coagulation, and high temperature tolerance, was negative for hemolysis, and was sensitive to clindamycin. In vivo studies in mice showed that the serum superoxide dismutase activity level was 53.69 U/mL higher in the MR group of mice compared to that of the control group, the malondialdehyde content was 0.53 nmol/mL lower in the HR group, and the highest jejunal V/C value was 4.11 ± 1.05 in the HR group (p < 0.05). The L. reuteri LRA7 gene is 2.021 megabases in size, contains one chromosome and one plasmid, and is annotated with 1978 functional genes. In conclusion, L. reuteri LRA7 has good probiotic potential and is safe. It can be used as an ideal probiotic candidate strain of canine origin.

Whole-genome sequence of Sclerotium delphinii, a pathogenic fungus of Dendrobium officinale southern blight

Dendrobium officinale is a rare and precious medicinal plant. Southern blight is a destructive disease in the artificial cultivation of D. officinale, and one of its pathogens is Sclerotium delphinii. S. delphinii is a phytopathogenic fungus with a wide host range with extremely strong pathogenicity. In this study, S. delphinii was isolated from D. officinale with southern blight. Subsequently, this specific strain underwent thorough whole-genome sequencing using the PacBio Sequel II platform, which employed single-molecule real-time (SMRT) technology. Comprehensive annotations were obtained through functional annotation of protein sequences using various publicly available databases. The genome of S. delphinii measures 73.66 Mb, with an N90 contig size of 2,707,110 bp, and it contains 18,506 putative predictive genes. This study represents the first report on the genome size assembly and annotation of S. delphinii, making it the initial species within the Sclerotium genus to undergo whole-genome sequencing, which can provide solid data and a theoretical basis for further research on the pathogenesis, omics of S. delphinii.

Association of pathogenic determinants of Fusobacterium necrophorum with bacteremia, and Lemierre's syndrome

Fusobacterium necrophorum is a Gram-negative anaerobic bacterium responsible for localized infections of the oropharynx that can evolve into bacteremia and/or septic thrombophlebitis of the jugular vein or peritonsillar vein, called Lemierre's syndrome. To identify microbial genetic determinants associated with the severity of this life-threatening disease, 70 F. necrophorum strains were collected and grouped into two categories according to the clinical presentation: (i) localized infection, (ii) bacteremia with/without Lemierre's syndrome. Comparative genomic analyses revealed two clades with distinct genetic content, one clade being significantly enriched with isolates from subjects with bacteremia. To identify genetic determinants contributing to F. necrophorum pathogenicity, genomic islands and virulence factor orthogroups (OVFs) were predicted. The presence/absence profiles of OVFs did not group isolates according to their clinical category, but rather according to their phylogeny. However, a variant of lktA, a key virulence factor, with a frameshift deletion that results in two open reading frames, was associated with bacteremia. Moreover, a genome-wide association study identified three orthogroups associated with bacteremic strains: (i) cas8a1, (ii) a sodium/solute symporter, and (iii) a POP1 domain-containing protein. Further studies must be performed to assess the functional impact of lktA mutation and of these orthogroups on the physiopathological mechanisms of F. necrophorum infections.

Dual-RNA-sequencing to elucidate the interactions between sorghum and Colletotrichum sublineola

In warm and humid regions, the productivity of sorghum is significantly limited by the fungal hemibiotrophic pathogen Colletotrichum sublineola, the causal agent of anthracnose, a problematic disease of sorghum (Sorghum bicolor (L.) Moench) that can result in grain and biomass yield losses of up to 50%. Despite available genomic resources of both the host and fungal pathogen, the molecular basis of sorghum-C. sublineola interactions are poorly understood. By employing a dual-RNA sequencing approach, the molecular crosstalk between sorghum and C. sublineola can be elucidated. In this study, we examined the transcriptomes of four resistant sorghum accessions from the sorghum association panel (SAP) at varying time points post-infection with C. sublineola. Approximately 0.3% and 93% of the reads mapped to the genomes of C. sublineola and Sorghum bicolor, respectively. Expression profiling of in vitro versus in planta C. sublineola at 1-, 3-, and 5-days post-infection (dpi) indicated that genes encoding secreted candidate effectors, carbohydrate-active enzymes (CAZymes), and membrane transporters increased in expression during the transition from the biotrophic to the necrotrophic phase (3 dpi). The hallmark of the pathogen-associated molecular pattern (PAMP)-triggered immunity in sorghum includes the production of reactive oxygen species (ROS) and phytoalexins. The majority of effector candidates secreted by C. sublineola were predicted to be localized in the host apoplast, where they could interfere with the PAMP-triggered immunity response, specifically in the host ROS signaling pathway. The genes encoding critical molecular factors influencing pathogenicity identified in this study are a useful resource for subsequent genetic experiments aimed at validating their contributions to pathogen virulence. This comprehensive study not only provides a better understanding of the biology of C. sublineola but also supports the long-term goal of developing resistant sorghum cultivars.

Characteristics of Corynespora cassiicola, the causal agent of tobacco Corynespora leaf spot, revealed by genomic and metabolic phenomic analysis

Corynespora cassiicola is a highly diverse fungal pathogen that commonly occurs in tropical, subtropical, and greenhouse environments worldwide. In this study, the isolates were identified as C. cassiicola, and the optimum growth and sporulation were studied. The phenotypic characteristics of C. cassiicola, concerning 950 different growth conditions, were tested using Biolog PM plates 1-10. In addition, the strain of C. cassiicola DWZ from tobacco hosts was sequenced for the using Illumina PE150 and Pacbio technologies. The host resistance of tobacco Yunyan 87 with different maturity levels was investigated. In addition, the resistance evaluation of 10 common tobacco varieties was investigated. The results showed that C. cassiicola metabolized 89.47% of the tested carbon source, 100% of the nitrogen source, 100% of the phosphorus source, and 97.14% of the sulfur source. It can adapt to a variety of different osmotic pressure and pH environments, and has good decarboxylase and deaminase activities. The optimum conditions for pathogen growth and sporulation were 25-30 °C, and the growth was better on AEA and OA medium. The total length of the genome was 45.9 Mbp, the GC content was 51.23%, and a total of 13,061 protein-coding genes, 202 non-coding RNAs and 2801 and repeat sequences were predicted. Mature leaves were more susceptible than proper mature and immature leaves, and the average diameter of diseased spots reached 17.74 mm at 12 days. None of the tested ten cultivars exhibited obvious resistance to Corynespora leaf spot of tobacco, whereby all disease spot diameters reached > 10 mm and > 30 mm when at 5 and 10 days after inoculation, respectively. The phenotypic characteristics, genomic analysis of C. cassiicola and the cultivar resistance assessment of this pathogen have increased our understanding of Corynespora leaf spot of tobacco.

Biodegradation of atrazine and nicosulfuron by Streptomyces nigra LM01: Performance, degradative pathway, and possible genes involved

Microbial herbicide degradation is an efficient bioremediation method. In this study, a strain of Streptomyces nigra, LM01, which efficiently degrades atrazine and nicosulfuron, was isolated from a corn field using a direct isolation method. The degradation effects of the identified strain on two herbicides were investigated and optimized using an artificial neural network. The maximum degradation rates of S. nigra LM01 were 58.09 % and 42.97 % for atrazine and nicosulfuron, respectively. The degradation rate of atrazine in the soil reached 67.94 % when the concentration was 108 CFU/g after 5 d and was less effective than that of nicosulfuron. Whole genome sequencing of strain LM01 helped elucidate the possible degradation pathways of atrazine and nicosulfuron. The protein sequences of strain LM01 were aligned with the sequences of the degraded proteins of the two herbicides by using the National Center for Biotechnology Information platform. The sequence (GE005358, GE001556, GE004212, GE005218, GE004846, GE002487) with the highest query cover was retained and docked with the small-molecule ligands of the herbicides. The results revealed a binding energy of - 6.23 kcal/mol between GE005358 and the atrazine ligand and - 6.66 kcal/mol between GE002487 and the nicosulfuron ligand.

Genome-Wide Identification and Characterization of Homeobox Transcription Factors in Phoma sorghina var. saccharum Causing Sugarcane Twisted Leaf Disease

A homeobox transcription factor is a conserved transcription factor, ubiquitous in eukaryotes, that regulates the tissue formation of structure, cell differentiation, proliferation, and cancer. This study identified the homeobox transcription factor family and its distribution in Phoma sorghina var. saccharum at the whole genome level. It elucidated the gene structures and evolutionary characteristics of this family. Additionally, knockout experiments were carried out and the preliminary function of these transcription factors was studied. Through bioinformatics approaches, nine homeobox transcription factors (PsHOX1-PsHOX9) were identified in P. sorghina var. saccharum, and these contained HOX-conserved domains and helix-turn-helix secondary structures. Nine homeobox gene deletion mutants were obtained using the homologous recombinant gene knockout technique. Protoplast transformation was mediated by polyethylene glycol (PEG) and the transformants were identified using PCR. The knockouts of PsHOX1, PsHOX2, PsHOX3, PsHOX4, PsHOX6, PsHOX8, and PsHOX9 genes resulted in a smaller growth diameter in P. sorghina var. saccharum. In contrast, the knockouts of the PsHOX3, PsHOX6, and PsHOX9 genes inhibited the formation of conidia and led to a significant decrease in the pathogenicity. This study's results will provide insights for understanding the growth and development of P. sorghina var. saccharum. The pathogenic mechanism of the affected sugarcane will provide an essential theoretical basis for preventing and controlling sugarcane twisted leaf disease.

A polylactic acid degrading lipase from Bacillus safensis: Characterization and structural analysis

A polylactic acid degrading triacylglycerol lipase (TGL) was identified from Bacillus safensis based on genome annotation and validated by real-time quantitative PCR. TGL displayed optimal activity at pH 9.0 and 55 °C. It maintained stability at pH 9.0 and temperatures 45 °C. The activity of TGL was found to benefit from the presence of potassium sodium ions, and low concentrations of Triton X-100. The TGL could erode the surface of polylactic acid films and increase its hydrophilicity. The hydrolysis products of polylactic acid by TGL were lactate monomer and dimer. TGL contains a classical catalytic triad structure of lipase (Ser77, Asp133, and His156) and an Ala-X-Ser-X-Gly sequence. Compared with some lipases produced by the same genus Bacillus, TGL is highly conserved in its amino acid sequence, mainly reflected in the amino acid residues that exercise the enzyme activity, including the catalytic activity center and the substrate binding sites.

Genome sequencing and comparative genome analysis of Rhizoctonia solani AG-3

Rhizoctonia solani AG-3 is a plant pathogenic fungus that belongs to the group of multinucleate Rhizoctonia. According to its internal transcribed spacer (ITS) cluster analysis and host range, it is divided into TB, PT, and TM subgroups. AG-3 TB mainly causes tobacco target spots, AG-3 PT mainly causes potato black scurf, and AG-3 TM mainly causes tomato leaf blight. In our previous study, we found that all 36 tobacco target spot strains isolated from Yunnan (Southwest China) were classified into AG-3 TB subgroup, while only two of the six tobacco target spot strains isolated from Liaoning (Northeast China) were classified into AG-3 TB subgroup, and the remaining four strains were classified into AG-3 TM subgroup, which had a unique taxonomic status, and there was no previous report on the whole genome information of AG-3 TM subgroup. In this study, the whole genomes of R. solani AG-3 strains 3T-1 (AG-3 TM isolated from Liaoning) and MJ-102 (AG-3 TB isolated from Yunnan) isolated from tobacco target spot in Liaoning and Yunnan were sequenced by IIumina and PacBio sequencing platforms. Comparative genomic analysis was performed with the previously reported AG-3 PT strain Rhs1AP, revealing their differences in genomes and virulence factors. The results indicated that the genome size of 3T-1 was 42,103,597 bp with 11,290 coding genes and 49.74% GC content, and the genome size of MJ-102 was 41,908,281 bp with 10,592 coding genes and 48.91% GC content. Through comparative genomic analysis with the previously reported strain Rhs1AP (AG-3 PT), it was found that the GC content between the genomes was similar, but the strains 3T-1 and MJ-102 contained more repetitive sequences. Similarly, there are similarities between their virulence factors, but there are also some differences. In addition, the results of collinearity analysis showed that 3T-1 and MJ-102 had lower similarity and longer evolutionary distance with Rhs1AP, but the genetic relationship between 3T-1 and MJ-102 was closer. This study can lay a foundation for studying the molecular pathogenesis and virulence factors of R. solani AG-3, and revealing its genomic composition will also help to develop more effective disease control strategies.

Interaction of tetracycline and copper co-intake in inducing antibiotic resistance genes and potential pathogens in mouse gut

The widespread use of copper and tetracycline as growth promoters in the breeding industry poses a potential threat to environmental health. Nevertheless, to the best of our knowledge, the potential adverse effects of copper and tetracycline on the gut microbiota remain unknown. Herein, mice were fed different concentrations of copper and/or tetracycline for 6 weeks to simulate real life-like exposure in the breeding industry. Following the exposure, antibiotic resistance genes (ARGs), potential pathogens, and other pathogenic factors were analyzed in mouse feces. The co-exposure of copper with tetracycline significantly increased the abundance of ARGs and enriched more potential pathogens in the gut of the co-treated mice. Copper and/or tetracycline exposure increased the abundance of bacteria carrying either ARGs, metal resistance genes, or virulence factors, contributing to the widespread dissemination of potentially harmful genes posing a severe risk to public health. Our study provides insights into the effects of copper and tetracycline exposure on the gut resistome and potential pathogens, and our findings can help reduce the risks associated with antibiotic resistance under the One Health framework.

Whole-genome sequencing and analysis of Chryseobacterium arthrosphaerae from Rana nigromaculata

Chryseobacterium arthrosphaerae strain FS91703 was isolated from Rana nigromaculata in our previous study. To investigate the genomic characteristics, pathogenicity-related genes, antimicrobial resistance, and phylogenetic relationship of this strain, PacBio RS II and Illumina HiSeq 2000 platforms were used for the whole genome sequencing. The genome size of strain FS91703 was 5,435,691 bp and GC content was 37.78%. A total of 4,951 coding genes were predicted; 99 potential virulence factors homologs were identified. Analysis of antibiotic resistance genes revealed that strain FS91703 harbored 10 antibiotic resistance genes in 6 categories and 2 multidrug-resistant efflux pump genes, including adeG and farA. Strain FS91703 was sensitive to β-lactam combination drugs, cephem, monobactam and carbapenems, intermediately resistant to phenicol, and resistant to penicillin, aminoglycosides, tetracycline, fluoroquinolones, and folate pathway inhibitors. Phylogenetic analysis revealed that strain FS91703 and C. arthrosphaerae CC-VM-7T were on the same branch of the phylogenetic tree based on 16 S rRNA; the ANI value between them was 96.99%; and the DDH values were 80.2, 72.2 and 81.6% by three default calculation formulae. These results suggested that strain FS91703 was a species of C. arthrosphaerae. Pan-genome analysis showed FS91703 had 566 unique genes compared with 13 other C. arthrosphaerae strains, and had a distant phylogenetic relationship with the other C. arthrosphaerae strains of the same branch in phylogenetic tree based on orthologous genes. The results of this study suggest that strain FS91703 is a multidrug-resistant and highly virulent bacterium, that differs from other C. arthrosphaerae strains at the genomic level. The knowledge about the genomic characteristics and antimicrobial resistance of strain FS91703 provides valuable insights into this rare species, as well as guidance for the treatment of the disease caused by FS91703 in Rana nigromaculata.

Integrated Genome Sequencing and Transcriptome Analysis Identifies Candidate Pathogenicity Genes from Ustilago crameri

Ustilago crameri is a pathogenic basidiomycete fungus that causes foxtail millet kernel smut (FMKS), a devastating grain disease in most foxtail-millet-growing regions of the world. Here, we report an assembled high-quality genome sequence of U. crameri strain SCZ-6 isolated from the diseased grains of foxtail millet in Changzhi, Shanxi Province, China. The genome size is 19.55 Mb, consisting of 73 contigs (N50 = 840,209 bp) with a G + C content of 54.09%, and encoding 6576 predicted genes and 6486 genes supported by RNA-seq. Evolutionarily, U. crameri lies close to the barley smut U. hordei, and an obvious co-linearity was observed between these two smut fungi. We annotated the genome of U. crameri strain SCZ-6 using databases, identifying 1827 pathogen-host interaction (PHI)-associated genes, 1324 genes encoding fungal virulence factors, 259 CAZy-related genes, 80 genes encoding transporters, and 206 putative cytochrome P450 genes; their expression profiles at different inoculation time points were also detected. Additionally, 70 candidate pathogen effectors were identified according to their expression patterns and predicted functions. In summary, our results provide important insights into the pathogenic mechanisms of the pathogenesis-related genes of U. crameri and a robust foundation for further investigation.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

Understanding the mechanism of pathogenicity through interactome studies between Arachis hypogaea L. and Aspergillus flavus

Aspergillus flavus (A. flavus) infects the peanut seeds during pre-and post-harvest stages, causing seed quality destruction for humans and livestock consumption. Even though many resistant varieties were developed, the molecular mechanism of defense interactions of peanut against A. flavus still needs further investigation. Hence, an interologous host-pathogen protein interaction (HPPI) network was constructed to understand the subcellular level interaction mechanism between peanut and A. flavus. Out of the top 10 hub proteins of both organisms, protein phosphatase 2C and cyclic nucleotide-binding/kinase domain-containing protein and different ribosomal proteins were identified as candidate proteins involved in defense. Functional annotation and subcellular localization based characterization of HPPI identified protein SGT1 homolog, calmodulin and Rac-like GTP-binding proteins to be involved in defense response against fungus. The relevance of HPPI in infectious conditions was assessed using two transcriptome data which identified the interplay of host kinase class R proteins, bHLH TFs and cell wall related proteins to impart resistance against pathogen infection. Further, the pathogenicity analysis identified glycogen phosphorylase and molecular chaperone and allergen Mod-E/Hsp90/Hsp1 as potential pathogen targets to enhance the host defense mechanism. Hence, the computationally predicted host-pathogen PPI network could provide valuable support for molecular biology experiments to understand the host-pathogen interaction. SIGNIFICANCE: Protein-protein interactions execute significant cellular interactions in an organism and are influenced majorly by stress conditions. Here we reported the host-pathogen protein-protein interaction between peanut and A. flavus, and a detailed network analysis based on function, subcellular localization, gene co-expression, and pathogenicity was performed. The network analysis identified key proteins such as host kinase class R proteins, calmodulin, SGT1 homolog, Rac-like GTP-binding proteins bHLH TFs and cell wall related to impart resistance against pathogen infection. We observed the interplay of defense related proteins and cell wall related proteins predominantly, which could be subjected to further studies. The network analysis described in this study could be applied to understand other host-pathogen systems generally.

Whole Genome Sequence of Lactiplantibacillus plantarum HOM3204 and Its Antioxidant Effect on D-Galactose-Induced Aging in Mice

Lactiplantibacillus plantarum, previously named Lactobacillus plantarum, is a facultative, homofermentative lactic acid bacterium widely distributed in nature. Several Lpb. plantarum strains have been demonstrated to possess good probiotic properties, and Lpb. plantarum HOM3204 is a potential probiotic strain isolated from homemade pickled cabbage plants. In this study, whole-genome sequencing was performed to acquire genetic information and predict the function of HOM3204, which has a circular chromosome of 3,232,697 bp and two plasmids of 48,573 and 17,060 bp, respectively. Moreover, various oxidative stress-related genes were identified in the strain, and its antioxidant activity was evaluated in vitro and in vivo. Compared to reference strains, the intracellular cell-free extracts of Lpb. plantarum HOM3204 at a dose of 1010 colony-forming units (CFU)/ml in vitro exhibited stronger antioxidant properties, such as total antioxidant activity, 2,2-diphenyl-1-picrylhydrazyl radical scavenging rate, superoxide dismutase activity, and glutathione (GSH) content. Daily administration of 109 CFU Lpb. plantarum HOM3204 for 45 days significantly improved the antioxidant function by increasing the glutathione peroxidase activity in the whole blood and GSH concentration in the livers of D-galactose-induced aging mice. These results suggest that Lpb. plantarum HOM3204 can potentially be used as a food ingredient with good antioxidant properties.

Multiomics analysis reveals the molecular mechanisms underlying virulence in Rhizoctonia and jasmonic acid-mediated resistance in Tartary buckwheat (Fagopyrum tataricum)

Rhizoctonia solani is a devastating soil-borne pathogen that seriously threatens the cultivation of economically important crops. Multiple strains with a very broad host range have been identified, but only 1 (AG1-IA, which causes rice sheath blight disease) has been examined in detail. Here, we analyzed AG4-HGI 3 originally isolated from Tartary buckwheat (Fagopyrum tataricum), but with a host range comparable to AG1-IA. Genome comparison reveals abundant pathogenicity genes in this strain. We used multiomic approaches to improve the efficiency of screening for disease resistance genes. Transcriptomes of the plant-fungi interaction identified differentially expressed genes associated with virulence in Rhizoctonia and resistance in Tartary buckwheat. Integration with jasmonate-mediated transcriptome and metabolome changes revealed a negative regulator of jasmonate signaling, cytochrome P450 (FtCYP94C1), as increasing disease resistance probably via accumulation of resistance-related flavonoids. The integration of resistance data for 320 Tartary buckwheat accessions identified a gene homolog to aspartic proteinase (FtASP), with peak expression following R. solani inoculation. FtASP exhibits no proteinase activity but functions as an antibacterial peptide that slows fungal growth. This work reveals a potential mechanism behind pathogen virulence and host resistance, which should accelerate the molecular breeding of resistant varieties in economically essential crops.

A framework for community curation of interspecies interactions literature

The quantity and complexity of data being generated and published in biology has increased substantially, but few methods exist for capturing knowledge about phenotypes derived from molecular interactions between diverse groups of species, in such a way that is amenable to data-driven biology and research. To improve access to this knowledge, we have constructed a framework for the curation of the scientific literature studying interspecies interactions, using data curated for the Pathogen-Host Interactions database (PHI-base) as a case study. The framework provides a curation tool, phenotype ontology, and controlled vocabularies to curate pathogen-host interaction data, at the level of the host, pathogen, strain, gene, and genotype. The concept of a multispecies genotype, the 'metagenotype,' is introduced to facilitate capturing changes in the disease-causing abilities of pathogens, and host resistance or susceptibility, observed by gene alterations. We report on this framework and describe PHI-Canto, a community curation tool for use by publication authors.

VmPacC-mediated pH regulation of Valsa mali confers to host acidification identified by comparative proteomics analysis

Apple valsa canker caused by the Ascomycete fungus Valsa mali is one of the most serious diseases of apple, resulting in huge economic losses in the apple-growing area of China. Previous study found that the pathogen could acidify the infected tissues to make lower ambient pH (from 6.0 to 3.5) for their successfully colonization. The pH signaling transcription factor VmPacC is required for acidification of its environment and for full virulence in V. mali. It is known that the functional cooperation of proteins secreted by V. mali plays pivotal role in its successful colonization of host plants. In this study, we used tandem mass tag (TMT) labeling coupled with LC-MS/MS-based quantitative proteomics to analyze the VmPacC-mediated pH regulation in V. mali, focusing on differentially expressed proteins (DEPs). We identified 222 DEPs specific to VmPacC deletion, and 921 DEPs specific to different pH conditions (pH 6.0 and 3.4). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that these DEPs were mainly involved in pathways associated with carbon metabolism, biosynthesis of antibiotics, citrate cycle (TCA cycle), glycolysis/gluconeogenesis, glutathione metabolism, ribosomes, and pentose phosphate pathways. Additionally, we identified 119 DEPs that were shared among the VmPacC deletion mutant and different pH conditions, which were mainly related to energy metabolism pathways, providing the energy required for the hyphal growth and responses to environmental stresses. A protein-protein interaction (PPI) network analysis indicated that most of the shared proteins were mapped to an interaction network with a medium confidence score of 0.4. Notably, one uncharacterized protein (KUI69106.1), and two known proteins (heat shock protein 60 (KUI73579.1), aspartate aminotransferase (KUI73864.1)) located in the core of the network were highly connected (with ≥ 38 directed edges) with the other shared DEPs. Our results suggest that VmPacC participates in the pathogen's regulation to ambient pH through the regulation of energy metabolism pathways such as the glycolysis/gluconeogenesis pathway and TCA cycle. Finally, we proposed a sophisticated molecular regulatory network to explain pH decrease in V. mali. Our study, by providing insights into V. mali regulating pH, helps to elucidate the mechanisms of host acidification during pathogen infection.

High-quality genome assembly of Verticillium dahliae VD991 allows for screening and validation of pathogenic genes

Verticillium dahliae (V. dahliae) is a notorious soil-borne pathogen causing Verticillium wilt in more than 400 dicotyledonous plants, including a wide range of economically important crops, such as cotton, tomato, lettuce, potato, and romaine lettuce, which can result in extensive economic losses. In the last decade, several studies have been conducted on the physiological and molecular mechanisms of plant resistance to V. dahliae. However, the lack of a complete genome sequence with a high-quality assembly and complete genomic annotations for V. dahliae has limited these studies. In this study, we produced a full genomic assembly for V. dahliae VD991 using Nanopore sequencing technology, consisting of 35.77 Mb across eight pseudochromosomes and with a GC content of 53.41%. Analysis of the genome completeness assessment (BUSCO alignment: 98.62%; Illumina reads alignment: 99.17%) indicated that our efforts resulted in a nearly complete and high-quality genomic assembly. We selected 25 species closely related to V. dahliae for evolutionary analysis, confirming the evolutionary relationship between V. dahliae and related species, and the identification of a possible whole genome duplication event in V. dahliae. The interaction between cotton and V. dahliae was investigated by transcriptome sequencing resulting in the identification of many genes and pathways associated with cotton disease resistance and V. dahliae pathogenesis. These results will provide new insights into the pathogenic mechanisms of V. dahliae and contribute to the cultivation of cotton varieties resistant to Verticillium wilt.

A small cysteine-rich fungal effector, BsCE66 is essential for the virulence of Bipolaris sorokiniana on wheat plants

The Spot Blotch (SB) caused by hemibiotrophic fungal pathogen Bipolaris sorokiniana is one of the most devastating wheat diseases leading to 15-100% crop loss. However, the biology of Triticum-Bipolaris interactions and host immunity modulation by secreted effector proteins remain underexplored. Here, we identified a total of 692 secretory proteins including 186 predicted effectors encoded by B. sorokiniana genome. Gene Ontology categorization showed that these proteins belong to cellular, metabolic and signaling processes, and exhibit catalytic and binding activities. Further, we functionally characterized a cysteine-rich, B. sorokiniana Candidate Effector 66 (BsCE66) that was induced at 24-96 hpi during host colonization. The Δbsce66 mutant did not show vegetative growth defects or stress sensitivity compared to wild-type, but developed drastically reduced necrotic lesions upon infection in wheat plants. The loss-of-virulence phenotype was rescued upon complementing the Δbsce66 mutant with BsCE66 gene. Moreover, BsCE66 does not form homodimer and conserved cysteine residues form intra-molecular disulphide bonds. BsCE66 localizes to the host nucleus and cytosol, and triggers a strong oxidative burst and cell death in Nicotiana benthamiana. Overall, our findings demonstrate that BsCE66 is a key virulence factor that is necessary for host immunity modulation and SB disease progression. These findings would significantly improve our understanding of Triticum-Bipolaris interactions and assist in the development of SB resistant wheat varieties.

Intraspecific Comparative Analysis Reveals Genomic Variation of Didymella arachidicola and Pathogenicity Factors Potentially Related to Lesion Phenotype

Didymella arachidicola is one of the most important fungal pathogens, causing foliar disease and leading to severe yield losses of peanuts (Arachis hypogaea L.) in China. Two main lesion phenotypes of peanut web blotch have been identified as reticulation type (R type) and blotch type (B type). As no satisfactory reference genome is available, the genomic variations and pathogenicity factors of D. arachidicola remain to be revealed. In the present study, we collected 41 D. arachidicola isolates from 26 geographic locations across China (33 for R type and 8 for B type). The chromosome-scale genome of the most virulent isolate (YY187) was assembled as a reference using PacBio and Hi-C technologies. In addition, we re-sequenced 40 isolates from different sampling sites. Genome-wide alignments showed high similarity among the genomic sequences from the 40 isolates, with an average mapping rate of 97.38%. An average of 3242 SNPs and 315 InDels were identified in the genomic variation analysis, which revealed an intraspecific polymorphism in D. arachidicola. The comparative analysis of the most and least virulent isolates generated an integrated gene set containing 512 differential genes. Moreover, 225 genes individually or simultaneously harbored hits in CAZy-base, PHI-base, DFVF, etc. Compared with the R type reference, the differential gene sets from all B type isolates identified 13 shared genes potentially related to lesion phenotype. Our results reveal the intraspecific genomic variation of D. arachidicola isolates and pathogenicity factors potentially related to different lesion phenotypes. This work sets a genomic foundation for understanding the mechanisms behind genomic diversity driving different pathogenic phenotypes of D. arachidicola.

Interactomics in plant defence: progress and opportunities

Interactomics is a branch of systems biology that deals with the study of protein-protein interactions and how these interactions influence phenotypes. Identifying the interactomes involved during host-pathogen interaction events may bring us a step closer to deciphering the molecular mechanisms underlying plant defence. Here, we conducted a systematic review of plant interactomics studies over the last two decades and found that while a substantial progress has been made in the field, plant-pathogen interactomics remains a less-travelled route. As an effort to facilitate the progress in this field, we provide here a comprehensive research pipeline for an in planta plant-pathogen interactomics study that encompasses the in silico prediction step to the validation step, unconfined to model plants. We also highlight four challenges in plant-pathogen interactomics with plausible solution(s) for each.

Genome Resource of Rhizoctonia solani Anastomosis Group 4 Strain AG4-JY, a Pathomycete of Sheath Blight of Foxtail Millet

The basidiomycetous fungus Rhizoctonia solani Kühn (teleomorph Thanatephorus cucumeris [Frank] Donk) is a fungal pathogen that causes various diseases on economically important crops, such as foxtail millet, maize, and rice. Using the PacBio Sequel platform, we assembled a draft genome of an R. solani strain AG4-JY that was isolated from foxtail millet with sheath blight at the stem. The genome was approximately 43.43 Mb on 53 scaffolds, with a scaffold N50 length of 2.10 Mb. In all, 10,545 genes and 179 noncoding RNAs were predicted, and 10,488 genes had at least one database annotation. In addition, the proteins encoded by 709 genes were predicted as secretory proteins. The AG4-JY genome sequence provides a valuable resource for understanding the interactions between R. solani and foxtail millet and controls sheath blight in the world.

Metagenomics revealing molecular profiles of microbial community structure and metabolic capacity in Bamucuo lake, Tibet

Microorganisms play critical ecological roles in the global biogeochemical cycles. However, extensive information on the microbial communities in Qinghai-Tibet Plateau (QTP), which is the highest plateau in the world, is still lacking, particularly in high elevation locations above 4500 m. Here, we performed a survey of th e soil and water microbial communities in Bamucuo Lake, Tibet, by using shotgun metagenomic methods. In the soil and water samples, we reconstructed 75 almost complete metagenomic assembly genomes, and 74 of the metagenomic assembly genomes from the water sample represented novel species. Proteobacteria and Actinobacteria were found to be the dominant bacterial phyla, while Euryarchaeota was the dominant archaeal phylum. The largest virus, Pandoravirus salinus, was found in the soil microbial community. We concluded that the microorganisms in Bamucuo Lake are most likely to fix carbon mainly through the 3-hydroxypropionic bi-cycle pathway. This study, for the first time, characterized the microbial community composition and metabolic capacity in QTP high-elevation locations with 4555 m, confirming that QTP is a vast and valuable resource pool, in which many microorganisms can be used to develop new bioactive substances and new antibiotics to which pathogenic microorganisms have not yet developed resistance.

Virulence and pathogenicity determinants in whole genome sequence of Fusarium udum causing wilt of pigeon pea

The present study deals with whole genome analysis of Fusarium udum, a wilt causing pathogen of pigeon pea. The de novo assembly identified a total of 16,179 protein-coding genes, of which 11,892 genes (73.50%) were annotated using BlastP and 8,928 genes (55.18%) from KOG annotation. In addition, 5,134 unique InterPro domains were detected in the annotated genes. Apart from this, we also analyzed genome sequence for key pathogenic genes involved in virulence, and identified 1,060 genes (6.55%) as virulence genes as per the PHI-BASE database. The secretome profiling of these virulence genes indicated the presence of 1,439 secretory proteins. Of those, an annotation of 506 predicted secretory proteins through CAZyme database indicated maximum abundance of Glycosyl hydrolase (GH, 45%) family proteins followed by auxiliary activity (AA) family proteins. Interestingly, the presence of effectors for cell wall degradation, pectin degradation, and host cell death was found. The genome comprised approximately 895,132 bp of repetitive elements, which includes 128 long terminal repeats (LTRs), and 4,921 simple sequence repeats (SSRs) of 80,875 bp length. The comparative mining of effector genes among different Fusarium species revealed five common and two specific effectors in F. udum that are related to host cell death. Furthermore, wet lab experiment validated the presence of effector genes like SIX (for Secreted in Xylem). We conclude that deciphering the whole genome of F. udum would be instrumental in understanding evolution, virulence determinants, host-pathogen interaction, possible control strategies, ecological behavior, and many other complexities of the pathogen.

Whole Genome Sequencing and Comparative Genomics of Indian Isolates of Wheat Spot Blotch Pathogen Bipolaris sorokiniana Reveals Expansion of Pathogenicity Gene Clusters

Spot blotch is a highly destructive disease in wheat caused by the fungal pathogen Bipolaris sorokiniana (teleomorph, Cochliobolus sativus). It is prevalent in warm and humid areas, including Africa, Asia, Latin America, and the USA. In the present study, twelve isolates of B. sorokiniana were collected from wheat fields in three different geographical locations in India. The pathogenicity of seven sporulating isolates was assessed on 'DDK 1025', a spot blotch-susceptible wheat variety under greenhouse conditions. The isolate 'D2' illustrated the highest virulence, followed by 'SI' and 'BS52'. These three isolates were sequenced using the Illumina HiSeq1000 platform. The estimated genome sizes of the isolates BS52, D2, and SI were 35.19 MB, 39.32 MB, and 32.76 MB, with GC contents of 48.48%, 50.43%, and 49.42%, respectively. The numbers of pathogenicity genes identified in BS52, D2, and SI isolates were 2015, 2476, and 2018, respectively. Notably, the isolate D2 exhibited a relatively larger genome with expanded arsenals of Biosynthetic Gene Clusters (BGCs), CAZymes, secretome, and pathogenicity genes, which could have contributed to its higher virulence among the tested isolates. This study provides the first comparative genome analysis of the Indian isolates of B. sorokiniana using whole genome sequencing.

Comparative genomic analyses reveal genetic characteristics and pathogenic factors of Bacillus pumilus HM-7

Bacillus pumilus plays an important role in industrial application and biocontrol activities, as well as causing humans and plants disease, leading to economic losses and biosafety concerns. However, until now, the pathogenesis and underlying mechanisms of B. pumilus strains remain unclear. In our previous study, one representative isolate of B. pumilus named HM-7 has been recovered and proved to be the causal agent of fruit rot on muskmelon (Cucumis melo). Herein, we present a complete and annotated genome sequence of HM-7 that contains 4,111 coding genes in a single 3,951,520 bp chromosome with 41.04% GC content. A total of 3,481 genes were functionally annotated with the GO, COG, and KEGG databases. Pan-core genome analysis of HM-7 and 20 representative B. pumilus strains, as well as six closely related Bacillus species, discovered 740 core genes and 15,205 genes in the pan-genome of 21 B. pumilus strains, in which 485 specific-genes were identified in HM-7 genome. The average nucleotide identity (ANI), and whole-genome-based phylogenetic analysis revealed that HM-7 was most closely related to the C4, GR8, MTCC-B6033, TUAT1 and SH-B11 strains, but evolutionarily distinct from other strains in B. pumilus. Collinearity analysis of the six similar B. pumilus strains showed high levels of synteny but also several divergent regions for each strains. In the HM-7 genome, we identified 484 genes in the carbohydrate-active enzymes (CAZyme) class, 650 genes encoding virulence factors, and 1,115 genes associated with pathogen-host interactions. Moreover, three HM-7-specific regions were determined, which contained 424 protein-coding genes. Further investigation of these genes showed that 19 pathogenesis-related genes were mainly associated with flagella formation and secretion of toxic products, which might be involved in the virulence of strain HM-7. Our results provided detailed genomic and taxonomic information for the HM-7 strain, and discovered its potential pathogenic mechanism, which lay a foundation for developing effective prevention and control strategies against this pathogen in the future.

Streptomyces: Derived Active Extract Inhibits Candida albicans Biofilm Formation

Candida albicans is an opportunistic pathogen that causes biofilm-associated infections. C. albicans biofilms are known to display reduced susceptibility to antimicrobials and high rates of acquired antibiotic resistance, and biofilm forming in C. albicans further hampers treatment options and highlights the need for new antibiofilm strategies. Identifying active components from desert actinomycetes strains to inhibit the formation of C. albicans biofilms represents an effective treatment strategy. In this study, actinomycetes that can inhibit C. albicans biofilm formation were isolated from the Taklimakan Desert, and the underlying mechanisms were explored. After screening the anti-C.albicans biofilm activities of culture supernatants from 170 Actinomycete strains, six strains showed significant inhibition of C. albicans biofilm formation. Microscopic examination showed a reduction in biofilm formation of C. albicans treated with supernatants from actinomycetes. Scanning electron microscopy showed that the morphological changes in biofilm cells were caused by cell membrane rupture and cell material leakage. Then, C.albicans biofilms were destroyed by changing the content of extracellular polysaccharides or degrading extracellular DNA. Finally, a preliminary study on active substances extracted from a new species (TRM43335) showed that the substances that inhibited the formation of biofilms might be peptides. This study provides preliminary evidence that desert actinomyces strains have inhibitory effects on the biofilm development of C. albicans.

Virulence-Associated Genes of Calonectria ilicola, Responsible for Cylindrocladium Black Rot

The Cylindrocladium black rot caused by Calonectria ilicicola is a destructive disease affecting a broad range of crops. Herein, we study virulence-associated genes of C. ilicicolaCi14017 isolated from diseased peanut roots (Arachis hypogaea L.). Ci14017 was identified via phylogenetic analysis of the internal transcribed spacer region and standard Koch’s postulate testing. Virulence-associated genes were based on genome analyses and comparative analysis of transcriptome and proteome profiles of sensitive and resistant peanut cultivars. Ci14017 identified as C. ilicicola has a 66 Mb chromosome with 18,366 predicted protein-coding genes. Overall, 46 virulence-associated genes with enhanced expression levels in the sensitive cultivars were identified. Sequence analysis indicated that the 46 gene products included two merops proteins, eight carbohydrate-active enzymes, seven cytochrome P450 enzymes, eight lipases, and 20 proteins with multi-conserved enzyme domains. The results indicate a complex infection mechanism employed by Ci14017 for causing Cylindrocladium black rot in peanuts.

Comparative genomic analysis reveals cellulase plays an important role in the pathogenicity of Setosphaeria turcica f. sp. zeae

Setosphaeria turcica f. sp. zeae and S. turcica f. sp. sorghi, the two formae speciales of S. turcica, cause northern leaf blight disease of corn and sorghum, respectively, and often cause serious economic losses. They have obvious physiological differentiation and show complete host specificity. Host specificity is often closely related to pathogen virulence factors, including secreted protein effectors and secondary metabolites. Genomic sequencing can provide more information for understanding the virulence mechanisms of pathogens. However, the complete genomic sequence of S. turcica f. sp. sorghi has not yet been reported, and no comparative genomic information is available for the two formae speciales. In this study, S. turcica f. sp. zeae was predicted to have fewer secreted proteins, pathogen-host interaction (PHI) genes and carbohydrate-active enzymes (CAZys) than S. turcica f. sp. sorghi. Fifteen and 20 polyketide synthase (PKS) genes were identified in S. turcica f. sp. zeae and S. turcica f. sp. sorghi, respectively, which maintained high homology. There were eight functionally annotated effector protein-encoding genes specifically in S. turcica f. sp. zeae, among which the encoding gene StCEL2 of endo-1, 4-β-D-glucanase, an important component of cellulase, was significantly up-regulated during the interaction process. Finally, gluconolactone inhibited cellulase activity and decreased infection rate and pathogenicity, which indicates that cellulase is essential for maintaining virulence. These findings demonstrate that cellulase plays an important role in the pathogenicity of S. turcica f. sp. zeae. Our results also provide a theoretical basis for future research on the molecular mechanisms underlying the pathogenicity of the two formae speciales and for identifying any associated genes.

Comparative Analysis of Transcriptomes of Ophiostoma novo-ulmi ssp. americana Colonizing Resistant or Sensitive Genotypes of American Elm

The Ascomycete Ophiostoma novo-ulmi threatens elm populations worldwide. The molecular mechanisms underlying its pathogenicity and virulence are still largely uncharacterized. As part of a collaborative study of the O. novo-ulmi-elm interactome, we analyzed the O. novo-ulmi ssp. americana transcriptomes obtained by deep sequencing of messenger RNAs recovered from Ulmus americana saplings from one resistant (Valley Forge, VF) and one susceptible (S) elm genotypes at 0 and 96 h post-inoculation (hpi). Transcripts were identified for 6424 of the 8640 protein-coding genes annotated in the O. novo-ulmi nuclear genome. A total of 1439 genes expressed in planta had orthologs in the PHI-base curated database of genes involved in host-pathogen interactions, whereas 472 genes were considered differentially expressed (DEG) in S elms (370 genes) and VF elms (102 genes) at 96 hpi. Gene ontology (GO) terms for processes and activities associated with transport and transmembrane transport accounted for half (27/55) of GO terms that were significantly enriched in fungal genes upregulated in S elms, whereas the 22 GO terms enriched in genes overexpressed in VF elms included nine GO terms associated with metabolism, catabolism and transport of carbohydrates. Weighted gene co-expression network analysis identified three modules that were significantly associated with higher gene expression in S elms. The three modules accounted for 727 genes expressed in planta and included 103 DEGs upregulated in S elms. Knockdown- and knockout mutants were obtained for eight O. novo-ulmi genes. Although mutants remained virulent towards U. americana saplings, we identified a large repertoire of additional candidate O. novo-ulmi pathogenicity genes for functional validation by loss-of-function approaches.

Transcriptome profiles of genes related to growth and virulence potential in Vibrio alginolyticus treated with modified clay

Vibrio alginolyticus is a globally distributed opportunistic pathogen that causes different degrees of disease in various marine organisms, such as fish, shrimp and shellfish. At present, vibriosis caused by V. alginolyticus has a wide epidemic range and causes frequent outbreaks, resulting in substantial losses in aquaculture. According to previous studies, modified clay (MC) could effectively flocculate and reduce the density of Vibrio in water, but it is still unknown whether MC inhibits growth and how it affects virulence in bottom flocs. Here, we studied the response mechanism of V. alginolyticus in flocs treated with MC at the transcriptome level and verified the transcriptomic data combined with relevant physiological experiments and reverse transcription quantitative real-time PCR (RT-qPCR) for the first time. It was found that the morphology of Vibrio in the MC flocs changed, the membrane function was damaged, the antioxidant system was activated, and the material and energy metabolism also changed. In addition, MC could inhibit the expression of virulence factors of V. alginolyticus; for example, flagella, pilus, siderophores, quorum sensing, and other related genes were significantly downregulated. In general, MC effectively inhibited the growth of Vibrio and reduced its virulence potential in flocs, which could provide theoretical support for a new model of healthy aquaculture.

Comparative Genomic Analyses Provide Insight Into the Pathogenicity of Metschnikowia bicuspidata LNES0119

Metschnikowia bicuspidata is a globally distributed pathogenic yeast with a wide range of aquatic hosts. A new strain, M. bicuspidata LNES0119, isolated from the Chinese mitten crab Eriocheir sinensis, has caused a serious reduction in production and marked economic loss for the aquaculture industry in China. Therefore, the whole-genome sequence of M. bicuspidata LNES0119 was sequenced using Illumina and Oxford Nanopore technology; whole-genome annotation and comparative genomic analyses of this pathogen were performed as well. A high-quality genome of M. bicuspidata LNES0119 was 16.13 Mb in size, with six scaffolds and six contigs, and encoded 5,567 putative predicted genes. Of these, 1,467 genes shared substantial homology with genes in the pathogen–host interactions database. Comparative genomic analyses of three M. bicuspidata strains and one non-pathogenic yeast, M. aff. pulcherrima, showed 331 unique genes in M. bicuspidata LNES0119, 30 of which were putatively related to pathogenicity. Overall, we identified several meaningful characteristics related to pathogenicity and virulence that may play essential roles in the infection and pathogenicity of M. bicuspidata LNES0119. Our study will aid in identifying potential targets for further exploration of the molecular basis of the pathogenicity of M. bicuspidata as well as the therapeutic intervention of M. bicuspidata infection.

Seed-Borne Erwinia persicina Affects the Growth and Physiology of Alfalfa (Medicago sativa L.)

Seed-borne Erwinia persicina can be transmitted globally via alfalfa (Medicago sativa L.) seed trade, but there is limited information about the impact of this plant-pathogenic bacterium on alfalfa plants. In this study, strain Cp2, isolated from alfalfa seeds, was confirmed by whole-genome sequencing to belong to E. persicina. Subsequently, the effects of Cp2 on alfalfa growth and physiology were evaluated by constructing a rhizosphere infection model. Strain Cp2 had a strong inhibitory effect on the elongation and growth of alfalfa roots, which was very unfavorable to these perennial plants. Furthermore, an increased number of leaf spots and yellowing symptoms were observed in plants of the Cp2 group from day 10 to day 21 and the strain Cp2 was re-isolated from these leaves. Correlation between growth and photosynthetic parameters was analyzed and the significant decreases in fresh weight and root and plant lengths in the Cp2 group were related to the marked reduction of chlorophyll b, carotenoid, transpiration rate, and stomatal conductance of leaves (r &gt; 0.75). In addition, nine physiological indicators of root, stem, and leaf were measured in the plants 21 days after treatment with Cp2. The physiological response of root and leaf to Cp2 treatment was stronger than that of stem. The physiological indicators with the greatest response to Cp2 infection were further explored through principal component analysis, and superoxide dismutase, peroxidase, ascorbate peroxidase, and soluble protein showed the greatest changes in roots, stems, and leaves (P &lt; 0.001). Among tissues, the commonality was the change of soluble protein. Therefore, soluble protein is speculated to be a physiological marker during alfalfa–E. persicina interactions. These findings indicate that once E. persicina spreads from alfalfa seeds to the rhizosphere, it can invade alfalfa roots and cause disease. This study demonstrates that this plant pathogenic bacterium may be a potential threat to new environment when they spread via seed trade and these “dangerous hitchhikers” warrant further attention, especially in the study of bacterial diseases in pasture-based production systems.