Comparative genomics of plant pathogenic Diaporthe species and transcriptomics of Diaporthe caulivora during host infection reveal insights into pathogenic strategies of the genus

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.

Elucidation of mechanisms underlying active oxygen burst in Citrus sinensis after Diaporthe citri infection using transcriptome analysis

Introduction

Reactive oxygen species (ROS) generation is a common disease defense mechanism in plants. However, it is unclear whether Citrus host activates defense response against Diaporthe citri causing citrus melanose disease by producing ROS, and the underlying molecular mechanisms are unknown.

Methods

DAB staining and RNA-Seq technology were used to compare the active oxygen burst and differential gene expression, respectively, in uninfected and infected Citrus sinensis leaves at different time points during D. citri infection in vivo. The functions of CsRBOH (a significant DEG) were confirmed in N. benthamiana through the Agrobacterium-mediated transient expression system.

Results

DAB staining indicated that C. sinensis initiated defense against D. citri infection within 24 h by generating ROS. Illumina sequencing revealed 25,557 expressed genes of C. sinensis. The most upregulated DEGs (n = 1,570) were identified 72 h after fungal inoculation (sample denoted as CD72). In the CD72 vs. Cs (samples at 0 h after fungal inoculation) comparison, the KEGG pathway category with the highest number of genes (n = 62) and most significant enrichment was Protein processing in endoplasmic reticulum, followed by Glutathione metabolism and MAPK signaling pathway-plant. GO analysis revealed that the DEGs of CD72 vs. Cs related to active oxygen burst and chitin recognition were significantly grouped into the regulation of biological processes and molecular functions, with GO terms including response to ROS, response to fungus, and oxidoreductase activity. Remarkably, CsRBOH was significantly enriched in the GO and KEGG analyses, and its expression pattern in qRT-PCR and DAB staining results were consistent. Among the 63 ROS-related DEGs, HSP genes and genes associated with the peroxidase family were highly significant as revealed by protein-protein interaction networks. Furthermore, ROS accumulation, cell death, and upregulation of defense-related genes were observed in N. benthamiana leaves with CsRBOH expressed through the Agrobacterium-mediated transient expression system.

Conclusion

Our findings suggested that C. sinensis activates CsRBOH and ROS-related genes, leading to ROS accumulation to resist the invasion by D. citri. This study laid the foundation for future research on molecular mechanisms and breeding of C. sinensis cultivars resistant to citrus melanose.

A repertoire of candidate effector proteins of the fungus Ceratocystis cacaofunesta

The genus Ceratocystis includes many phytopathogenic fungi that affect different plant species. One of these is Ceratocystis cacaofunesta, which is pathogenic to the cocoa tree and causes Ceratocystis wilt, a lethal disease for the crop. However, little is known about how this pathogen interacts with its host. The knowledge and identification of possible genes encoding effector proteins are essential to understanding this pathosystem. The present work aimed to predict genes that code effector proteins of C. cacaofunesta from a comparative analysis of the genomes of five Ceratocystis species available in databases. We performed a new genome annotation through an in-silico analysis. We analyzed the secretome and effectorome of C. cacaofunesta using the characteristics of the peptides, such as the presence of signal peptide for secretion, absence of transmembrane domain, and richness of cysteine residues. We identified 160 candidate effector proteins in the C. cacaofunesta proteome that could be classified as cytoplasmic (102) or apoplastic (58). Of the total number of candidate effector proteins, 146 were expressed, presenting an average of 206.56 transcripts per million. Our database was created using a robust bioinformatics strategy, followed by manual curation, generating information on pathogenicity-related genes involved in plant interactions, including CAZymes, hydrolases, lyases, and oxidoreductases. Comparing proteins already characterized as effectors in Sordariomycetes species revealed five groups of protein sequences homologous to C. cacaofunesta. These data provide a valuable resource for studying the infection mechanisms of these pathogens in their hosts.

Comparative analysis of soybean transcriptional profiles reveals defense mechanisms involved in resistance against Diaporthe caulivora

Soybean stem canker (SSC) caused by the fungal pathogen Diaporthe caulivora is an important disease affecting soybean production worldwide. However, limited information related to the molecular mechanisms underlying soybean resistance to Diaporthe species is available. In the present work, we analyzed the defense responses to D. caulivora in the soybean genotypes Williams and Génesis 5601. The results showed that compared to Williams, Génesis 5601 is more resistant to fungal infection evidenced by significantly smaller lesion length, reduced disease severity and pathogen biomass. Transcriptional profiling was performed in untreated plants and in D. caulivora-inoculated and control-treated tissues at 8 and 48 h post inoculation (hpi). In total, 2.322 and 1.855 genes were differentially expressed in Génesis 5601 and Williams, respectively. Interestingly, Génesis 5601 exhibited a significantly higher number of upregulated genes compared to Williams at 8 hpi, 1.028 versus 434 genes. Resistance to D. caulivora was associated with defense activation through transcriptional reprogramming mediating perception of the pathogen by receptors, biosynthesis of phenylpropanoids, hormone signaling, small heat shock proteins and pathogenesis related (PR) genes. These findings provide novel insights into soybean defense mechanisms leading to host resistance against D. caulivora, and generate a foundation for the development of resistant SSC varieties within soybean breeding programs.

The Biochemical Response of Soybean Cultivars Infected by Diaporthe Species Complex

Oxidative stress in soybean plants infected with Diaporthe isolates was evaluated in order to select (1) the least aggressive inoculation method, (2) to determine the most aggressive Diaporthe isolate, and (3) to determine the most tolerant soybean cultivar to this isolate. Based on the present malondialdehyde (MDA) content, the main end product of the lipid peroxidation process, and the biomarker for oxidative stress, the mycelium contact method was chosen as the least aggressive inoculation method, compared to the toothpick method and plug method. The activity of the antioxidant enzymes (superoxide-dismutase (SOD), catalase (CAT), and peroxidase (PX)), the reduced glutathione (GSH) content, and the level of lipid peroxidation (LP) were measured in soybean cv. Sava infected by five different Diaporthe species (DPM1F-D. aspalathi, DPC/KR19-D. caulivora, DPC004NY15-D. eres, 18-DIA-SOY-14-D. gulyae, and PL157A-D. longicolla). The most pathogenic Diaporthe species to cv. Sava was D. eres. The screening of the antioxidant enzymes activity in the leaves of 12 different soybean cultivars (Altona, Atlas, Capital, Chico, CX134, Favorit, Lakota, McCall, Morsoy, Strain, Rubin, and Victoria) infected with D. eres by the mycelium contact inoculation method showed that Capital, McCall, and Morsoy were the cultivars with the highest tolerance to D. eres, followed by Chico, Favorit, Lakota, and Rubin. The most sensitive cultivars were Atlas, CX134, Victoria, and Strain.

The genome of a Far Eastern isolate of Diaporthe caulivora, a soybean fungal pathogen

Diaporthe caulivora is an economically important fungal pathogen and a causal agent of soybean stem canker and seed decay. Here, the genome of a Russian Far Eastern isolate of D. caulivora was sequenced, assembled, and announced. Assembly quality was enough for advanced annotation, including prediction of potential disease-related genes encoding virulence factors and molecular determinants contributing to pathogen-host selection, interactions, and adaptation. Comparative analysis of 15 Diaporthe species was conducted regarding general genome properties, collinearity, and proteomes, and included detailed investigation of interspersed repeats. A notable feature of this analysis is a high recombinant variability of Diaporthe genomes, determined by the number and distribution of interspersed repeats, which also proved to be responsible for the diversity of GC content and genome size. This variability is assumed the main determinant of the divergence of Diaporthe genomes. A Bayesian multi-gene phylogeny was inferred for the 15 Diaporthe species on the basis of twenty thousand polymorphic sites of > 100 orthologous genes using independently adjusted evolutionary models. This allowed for the most accurate determination of evolutionary relationships and species boundaries for effective reporting about these plant pathogens. The evidence, obtained by different genome analysis techniques, implies the host-independent evolution of Diaporthe species. KEY POINTS: • The genome of a Far Eastern isolate of D. caulivora was announced. • A high degree of recombinant variability determines genomic divergence in Diaporthe genus. • The multi-gene phylogeny implies host-independent evolution of Diaporthe species.

Pathogenicity Variation in Two Genomes of Cercospora Species Causing Gray Leaf Spot in Maize

The gray leaf spots caused by Cercospora spp. severely affect the yield and quality of maize. However, the evolutionary relation and pathogenicity variation between species of the Cercospora genus is largely unknown. In this study, we constructed high-quality reference genomes by nanopore sequencing two Cercospora species, namely, C. zeae-maydis and C. zeina, with differing pathogenicity, collected from northeast (Liaoning [LN]) and southeast (Yunnan [YN]) China, respectively. The genome size of C. zeae-maydis-LN is 45.08 Mb, containing 10,839 annotated genes, whereas that of Cercospora zeina-YN is 42.18 Mb, containing 10,867 annotated genes, of which approximately 86.58% are common in the two species. The difference in their genome size is largely attributed to increased long terminal repeat retrotransposons of 3.8 Mb in total length in C. zeae-maydis-LN. There are 41 and 30 carbohydrate-binding gene subfamilies identified in C. zeae-maydis-LN and C. zeina-YN, respectively. A higher number of carbohydrate-binding families found in C. zeae-maydis-LN, and its unique CBM4, CBM37, and CBM66, in particular, may contribute to variation in pathogenicity between the two species, as the carbohydrate-binding genes are known to encode cell wall-degrading enzymes. Moreover, there are 114 and 107 effectors predicted, with 47 and 46 having unique potential pathogenicity in C. zeae-maydis-LN and C. zeina-YN, respectively. Of eight effectors randomly selected for pathogenic testing, five were found to inhibit cell apoptosis induced by Bcl-2-associated X. Taken together, our results provide genomic insights into variation in pathogenicity between C. zeae-maydis and C. zeina. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Investigations of specialised metabolites of endophyte Diaporthe destruens hosted in Illigera orbiculata C. Y. Wu

A chemical investigation of the endophytic fungus Diaporthe destruens from the Hernandiaceae plant Illigera orbiculata C. Y. Wu collected from southern Yunnan Province, China, led to the isolation of six undescribed compounds, including two azaphilone analogs, which are a pair of epimers (13R-hydroxy-chermesinone A and 13S-hydroxy-chermesinone A); a pyrrole derivative (1-(4-(methoxymethyl)-1H-pyrrol-3-yl)ethan-1-one); an isoindolone derivative (4-hydroxy-6-methoxyisoindolin-1-one); a benzylbenzene derivative (destruensine A) and a conjectural fragment of polyketide ((2R,4R)-2-(methoxymethyl)pentane-1,4-diol) along with nine known compounds. Their structures were elucidated by spectroscopic methods and HRESIMS, and the absolute configurations were further confirmed by electronic circular dichroism (ECD) and chemical derivatization. The antimicrobial activities, anti-acetylcholinesterase activities, antiproliferation, and NO production inhibitory effects of compounds 1-15 were evaluated.

Genome Analyses of Two Blueberry Pathogens: Diaportheamygdali CAA958 and Diaporthe eres CBS 160.32

The genus Diaporthe includes pathogenic species distributed worldwide and affecting a wide variety of hosts. Diaporthe amygdali and Diaporthe eres have been found to cause cankers, dieback, or twig blights on economically important crops such as soybean, almond, grapevine, and blueberry. Despite their importance as plant pathogens, the strategies of species of Diaporthe to infect host plants are poorly explored. To provide a genomic basis of pathogenicity, the genomes of D. amygdali CAA958 and D. eres CBS 160.32 were sequenced and analyzed. Cellular transporters involved in the transport of toxins, ions, sugars, effectors, and genes implicated in pathogenicity were detected in both genomes. Hydrolases and oxidoreductases were the most prevalent carbohydrate-active enzymes (CAZymes). However, analyses of the secreted proteins revealed that the secretome of D. eres CBS 160.32 is represented by 5.4% of CAZymes, whereas the secreted CAZymes repertoire of D. amygdali CAA958 represents 29.1% of all secretomes. Biosynthetic gene clusters (BGCs) encoding compounds related to phytotoxins and mycotoxins were detected in D. eres and D. amygdali genomes. The core gene clusters of the phytotoxin Fusicoccin A in D. amygdali are reported here through a genome-scale assembly. Comparative analyses of the genomes from 11 Diaporthe species revealed an average of 874 CAZymes, 101 secondary metabolite BGCs, 1640 secreted proteins per species, and genome sizes ranging from 51.5 to 63.6 Mbp. This study offers insights into the overall features and characteristics of Diaporthe genomes. Our findings enrich the knowledge about D. eres and D. amygdali, which will facilitate further research into the pathogenicity mechanisms of these species.