Draft whole genome sequence of groundnut stem rot fungus Athelia rolfsii revealing genetic architect of its pathogenicity and virulence

Ca2+ affects the hyphal differentiation to sclerotia formation of Athelia rolfsii

RNA-Sequencing (RNA-Seq) and transcriptomic analyses have become powerful tools to study the developmental stages of fungal structures scuh as sclerotia. While RNA-Seq experiments have been set up for many important sclerotia- and microsclerotia-forming fungi, it has not been implemented to study Athelia rolfsii, which is one of the earliest fungi used in literature to uncover the roles of reactive oxygen species (ROS) in stimulating sclerotia formation. This study applied RNA-Seq to profile gene expression in four developmental stages of A. rolfsii sclerotia. Surprisingly, gene ontology and expression patterns suggested that most ROS-scavenging genes were not up-regulated in the stages from hyphal differentiation to the initial sclerotia stage. Using antioxidant and oxidant-amended culture assay, the results suggested none of the ascorbic acid, dithiothreitol (DTT), H2O2, or superoxide dismutase inhibitors [diethyldithiocarbamate (DETC), NaN3, and sodium dodecyl sulfate] affected the sclerotia number. Instead, only glutathione reduced the sclerotia number. Because glutathione has also been suggested to facilitate Ca2+ influx, therefore, glutathione culture assays with the combination of CaCl2, Ca2+-chelator egtazic acid, DETC, and H2O2 were tested on A. rolfsii, as well as two other fungi (Sclerotinia sclerotiorum and Macrophomina phaseolina) for comparison. Although the addition of CaCl2 caused sclerotia or microsclerotia reduction for all three fungi, the CaCl2-ROS interaction was only observed for S. sclerotiorum and M. phaseolina, but not A. rolfsi. Collectively, this study not only pointed out a conserved function of Ca2+ in suppressing fungal sclerotia and microsclerotia formation but also highlighted sclerotia formation of A. rolfsii being only sensitive to Ca2+ and independent of ROS stimuli.IMPORTANCEManagement for plant diseases caused by soil-borne fungal pathogens is challenging because many soil-borne fungal pathogens form sclerotia for long-term survival. Advanced understanding of the molecular and cellular mechanisms of sclerotia formation may provide novel insights to prevent these fungal residues in fields. This study discovered that Ca2+ acts as a negative signal cue to suppress sclerotia and microsclerotia formation in three economically important fungal pathogens. Moreover, the southern blight fungus Athelia rolfsii appears to be only regulated by Ca2+ but not reactive oxygen species. Accordingly, A. rolfsii can be a useful system for studying the detailed mechanism of Ca2+, and the applicability of Ca2+ in reducing sclerotia could be further assessed for disease management.

Bacillus lipopeptide-mediated biocontrol of peanut stem rot caused by Athelia rolfsii

INTRODUCTION

Peanut (Arachis hypogaea L.) is a widespread oilseed crop of high agricultural importance in tropical and subtropical areas. It plays a major role in the food supply in the Democratic Republic of Congo (DRC). However, one major constraint in the production of this plant is the stem rot (white mold or southern blight) disease caused by Athelia rolfsii which is so far controlled mainly using chemicals. Considering the harmful effect of chemical pesticides, the implementation of eco-friendly alternatives such as biological control is required for disease management in a more sustainable agriculture in the DRC as in the other developing countries concerned. Bacillus velezensis is among the rhizobacteria best described for its plant protective effect notably due to the production of a wide range of bioactive secondary metabolites. In this work, we wanted to evaluate the potential of B. velezensis strain GA1 at reducing A. rolfsii infection and to unravel the molecular basis of the protective effect.

RESULTS AND DISCUSSION

Upon growth under the nutritional conditions dictated by peanut root exudation, the bacterium efficiently produces the three types of lipopeptides surfactin, iturin and fengycin known for their antagonistic activities against a wide range of fungal phytopathogens. By testing a range of GA1 mutants specifically repressed in the production of those metabolites, we point out an important role for iturin and another unidentified compound in the antagonistic activity against the pathogen. Biocontrol experiments performed in greenhouse further revealed the efficacy of B. velezensis to reduce peanut disease caused by A. rolfsii both via direct antagonism against the fungus and by stimulating systemic resistance in the host plant. As treatment with pure surfactin yielded a similar level of protection, we postulate that this lipopeptide acts as main elicitor of peanut resistance against A. rolfsii infection.

Phenotypic analysis and genome sequence of Rhizopus oryzae strain Y5, the causal agent of tobacco pole rot

Rhizopus oryzae is a destructive pathogen that frequently causes tobacco pole rot in curing chambers. Phenotypic characterization of the pathogen was conducted to provide basic biological and pathological information using Biolog Phenotype MicroArray (PM). In addition, the Y5 strain of R. oryzae was sequenced using Illumina HiSeq and Pacific Biosciences (PacBio) technologies. Using PM plates 1-8, 758 growth conditions were tested. Results indicated that R. oryzae could metabolize 54.21% of tested carbon sources, 86.84% of nitrogen sources, 100% of sulfur sources, and 98.31% of phosphorus sources. About 37 carbon compounds, including D-xylose, N-acetyl-D-glucosamine, D-sorbitol, β-methyl-D-glucoside, D-galactose, L-arabinose, and D-cellobiose, significantly supported the growth of the pathogen. PM 3 indicated the active nitrogen sources, including Gly-Asn, Ala-Asp., Ala-Gln, and uric acid. PM 6-8 showed 285 different nitrogen pathways, indicating that different combinations of different amino acids support the growth of the pathogen. Genome sequencing results showed that the R. oryzae Y5 strain had raw data assembled into 2,271 Mbp with an N50 value of 10,563 bp. A genome sequence of 50.3 Mb was polished and assembled into 53 contigs with an N50 length of 1,785,794 bp, maximum contig length of 3,223,184 bp, and a sum of contig lengths of 51,182,778 bp. A total of 12,680 protein-coding genes were predicted using the Nonredundant, Gene Ontology, Clusters of Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes, and SWISS-PROT databases. The genome sequence and annotation resources of R. oryzae provided a reference for studying its biological characteristics, trait-specific genes, pathogen-host interaction, pathogen evolution, and population genetic diversity. The phenomics and genome of R. oryzae will provide insights into microfungal biology, pathogen evolution, and the genetic diversity of epidemics.

Genome-wide analysis of basic helix-loop-helix genes in Dendrobium catenatum and functional characterization of DcMYC2 in jasmonate-mediated immunity to Sclerotium delphinii

Dendrobium catenatum, belonging to the Orchidaceae, is a precious Chinese herbal medicine. Sclerotium delphinii (P1) is a broad-spectrum fungal disease, which causes widespread loss in the near-wild cultivation of D. catenatum. Thus, resistance breeding of D. catenatum has become the key to solve this problem. The basic helix-loop-helix (bHLH) gene family is closely related to plant resistance to external stresses, but the related research in D. catenatum is not deep enough yet. Phylogenetic analysis showed that 108 DcbHLH genes could be divided into 23 subgroups. Promoter cis-acting elements revealed that DcbHLHs contain a large number of stress-related cis-acting elements. Transcriptome analysis of MeJA and P1 treatment manifested that exogenous MeJA can change the expression pattern of most bHLH genes, especially the IIIe subgroup, including inhibiting the expression of DcbHLH026 (MYC2a) and promoting the expression of DcbHLH027 (MYC2b). Subcellular localization indicated that they were located in the nucleus. Furthermore, exogenous MeJA treatment significantly delayed disease time and reduced lesion size after infection with P1. DcMYC2b-overexpression Arabidopsis lines showed significantly smaller lesions after being infected with P1 than the wild type, indicating that DcMYC2b functions as an important positive regulator in D. catenatum defense against P1. Our findings shed more insights into the critical role of the DcbHLH family in plants and the resistance breeding of D. catenatum.

Structured Framework and Genome Analysis of Magnaporthe grisea Inciting Pearl Millet Blast Disease Reveals Versatile Metabolic Pathways, Protein Families, and Virulence Factors

Magnaporthe grisea (T.T. Herbert) M.E. Barr is a major fungal phytopathogen that causes blast disease in cereals, resulting in economic losses worldwide. An in-depth understanding of the basis of virulence and ecological adaptation of M. grisea is vital for devising effective disease management strategies. Here, we aimed to determine the genomic basis of the pathogenicity and underlying biochemical pathways in Magnaporthe using the genome sequence of a pearl millet-infecting M. grisea PMg_Dl generated by dual NGS techniques, Illumina NextSeq 500 and PacBio RS II. The short and long nucleotide reads could be draft assembled in 341 contigs and showed a genome size of 47.89 Mb with the N50 value of 765.4 Kb. Magnaporthe grisea PMg_Dl showed an average nucleotide identity (ANI) of 86% and 98% with M. oryzae and Pyricularia pennisetigena, respectively. The gene-calling method revealed a total of 10,218 genes and 10,184 protein-coding sequences in the genome of PMg_Dl. InterProScan of predicted protein showed a distinct 3637 protein families and 695 superfamilies in the PMg_Dl genome. In silico virulence analysis revealed the presence of 51VFs and 539 CAZymes in the genome. The genomic regions for the biosynthesis of cellulolytic endo-glucanase and beta-glucosidase, as well as pectinolytic endo-polygalacturonase, pectin-esterase, and pectate-lyases (pectinolytic) were detected. Signaling pathways modulated by MAPK, PI3K-Akt, AMPK, and mTOR were also deciphered. Multicopy sequences suggestive of transposable elements such as Type LTR, LTR/Copia, LTR/Gypsy, DNA/TcMar-Fot1, and Type LINE were recorded. The genomic resource presented here will be of use in the development of molecular marker and diagnosis, population genetics, disease management, and molecular taxonomy, and also provide a genomic reference for ascomycetous genome investigations in the future.

Structure and Function of Rhizosphere Soil and Root Endophytic Microbial Communities Associated With Root Rot of Panax notoginseng

Panax notoginseng (Burk.) F. H. Chen is a Chinese medicinal plant of the Araliaceae family used for the treatment of cardiovascular and cerebrovascular diseases in Asia. P. notoginseng is vulnerable to root rot disease, which reduces the yield of P. notoginseng. In this study, we analyzed the rhizosphere soil and root endophyte microbial communities of P. notoginseng from different geographical locations using high-throughput sequencing. Our results revealed that the P. notoginseng rhizosphere soil microbial community was more diverse than the root endophyte community. Rhodopseudomonas, Actinoplanes, Burkholderia, and Variovorax paradoxus can help P. notoginseng resist the invasion of root rot disease. Ilyonectria mors-panacis, Pseudomonas fluorescens, and Pseudopyrenochaeta lycopersici are pathogenic bacteria of P. notoginseng. The upregulation of amino acid transport and metabolism in the soil would help to resist pathogens and improve the resistance of P. notoginseng. The ABC transporter and gene modulating resistance genes can improve the disease resistance of P. notoginseng, and the increase in the number of GTs (glycosyltransferases) and GHs (glycoside hydrolases) families may be a molecular manifestation of P. notoginseng root rot. In addition, the complete genomes of two Flavobacteriaceae species and one Bacteroides species were obtained. This study demonstrated the microbial and functional diversity in the rhizosphere and root microbial community of P. notoginseng and provided useful information for a better understanding of the microbial community in P. notoginseng root rot. Our results provide insights into the molecular mechanism underlying P. notoginseng root rot and other plant rhizosphere microbial communities.

Hybrid de novo genome-reassembly reveals new insights on pathways and pathogenicity determinants in rice blast pathogen Magnaporthe oryzae RMg_Dl

Blast disease incited by Magnaporthe oryzae is a major threat to sustain rice production in all rice growing nations. The pathogen is widely distributed in all rice paddies and displays rapid aerial transmissions, and seed-borne latent infection. In order to understand the genetic variability, host specificity, and molecular basis of the pathogenicity-associated traits, the whole genome of rice infecting Magnaporthe oryzae (Strain RMg_Dl) was sequenced using the Illumina and PacBio (RSII compatible) platforms. The high-throughput hybrid assembly of short and long reads resulted in a total of 375 scaffolds with a genome size of 42.43 Mb. Furthermore, comparative genome analysis revealed 99% average nucleotide identity (ANI) with other oryzae genomes and 83% against M. grisea, and 73% against M. poe genomes. The gene calling identified 10,553 genes with 10,539 protein-coding sequences. Among the detected transposable elements, the LTR/Gypsy and Type LINE showed high occurrence. The InterProScan of predicted protein sequences revealed that 97% protein family (PFAM), 98% superfamily, and 95% CDD were shared among RMg_Dl and reference 70-15 genome, respectively. Additionally, 550 CAZymes with high GH family content/distribution and cell wall degrading enzymes (CWDE) such endoglucanase, beta-glucosidase, and pectate lyase were also deciphered in RMg_Dl. The prevalence of virulence factors determination revealed that 51 different VFs were found in the genome. The biochemical pathway such as starch and sucrose metabolism, mTOR signaling, cAMP signaling, MAPK signaling pathways related genes were identified in the genome. The 49,065 SNPs, 3267 insertions and 3611 deletions were detected, and majority of these varinats were located on downstream and upstream region. Taken together, the generated information will be useful to develop a specific marker for diagnosis, pathogen surveillance and tracking, molecular taxonomy, and species delineation which ultimately leads to device improved management strategies for blast disease.

JA signal-mediated immunity of Dendrobium catenatum to necrotrophic Southern Blight pathogen

BACKGROUND

Dendrobium catenatum belongs to the Orchidaceae, and is a precious Chinese herbal medicine. In the past 20 years, D. catenatum industry has developed from an endangered medicinal plant to multi-billion dollar grade industry. The necrotrophic pathogen Sclerotium delphinii has a devastating effection on over 500 plant species, especially resulting in widespread infection and severe yield loss in the process of large-scale cultivation of D. catenatum. It has been widely reported that Jasmonate (JA) is involved in plant immunity to pathogens, but the mechanisms of JA-induced plant resistance to S. delphinii are unclear.

RESULTS

In the present study, the role of JA in enhancing D. catenatum resistance to S. delphinii was investigated. We identified 2 COI1, 13 JAZ, and 12 MYC proteins in D. catenatum genome. Subsequently, systematic analyses containing phylogenetic relationship, gene structure, protein domain, and motif architecture of core JA pathway proteins were conducted in D. catenatum and the newly characterized homologs from its closely related orchid species Phalaenopsis equestris and Apostasia shenzhenica, along with the well-investigated homologs from Arabidopsis thaliana and Oryza sativa. Public RNA-seq data were investigated to analyze the expression patterns of D. catenatum core JA pathway genes in various tissues and organs. Transcriptome analysis of MeJA and S. delphinii treatment showed exogenous MeJA changed most of the expression of the above genes, and several key members, including DcJAZ1/2/5 and DcMYC2b, are involved in enhancing defense ability to S. delphinii in D. catenatum.

CONCLUSIONS

The findings indicate exogenous MeJA treatment affects the expression level of DcJAZ1/2/5 and DcMYC2b, thereby enhancing D. catenatum resistance to S. delphinii. This research would be helpful for future functional identification of core JA pathway genes involved in breeding for disease resistance in D. catenatum.

Genome sequencing and comparative genomic analysis of highly and weakly aggressive strains of Sclerotium rolfsii, the causal agent of peanut stem rot

BACKGROUND

Stem rot caused by Sclerotium rolfsii is a very important soil-borne disease of peanut. S. rolfsii is a necrotrophic plant pathogenic fungus with an extensive host range and worldwide distribution. It can infect peanut stems, roots, pegs and pods, leading to varied yield losses. S. rolfsii strains GP3 and ZY collected from peanut in different provinces of China exhibited a significant difference in aggressiveness on peanut plants by artificial inoculation test. In this study, de-novo genome sequencing of these two distinct strains was performed aiming to reveal the genomic basis of difference in aggressiveness.

RESULTS

Scleotium rolfsii strains GP3 and ZY, with weak and high aggressiveness on peanut plants, exhibited similar growth rate and oxalic acid production in laboratory. The genomes of S. rolfsii strains GP3 and ZY were sequenced by Pacbio long read technology and exhibited 70.51 Mb and 70.61 Mb, with contigs of 27 and 23, and encoded 17,097 and 16,743 gene models, respectively. Comparative genomic analysis revealed that the pathogenicity-related gene repertoires, which might be associated with aggressiveness, differed between GP3 and ZY. There were 58 and 45 unique pathogen-host interaction (PHI) genes in GP3 and ZY, respectively. The ZY strain had more carbohydrate-active enzymes (CAZymes) in its secretome than GP3, especially in the glycoside hydrolase family (GH), the carbohydrate esterase family (CBM), and the polysaccharide lyase family (PL). GP3 and ZY also had different effector candidates and putative secondary metabolite synthetic gene clusters. These results indicated that differences in PHI, secreted CAZymes, effectors and secondary metabolites may play important roles in aggressive difference between these two strains.

CONCLUSIONS

The data provided a further understanding of the S. rolfsii genome. Genomic comparison provided clues to the difference in aggressiveness of S. rolfsii strains.

Optimization of scleroglucan production by Sclerotium rolfsii by lowering pH during fermentation via oxalate metabolic pathway manipulation using CRISPR/Cas9

Background

Sclerotium rolfsii is a potent producer of many secondary metabolites, one of which like scleroglucan is an exopolysaccharide (EPS) appreciated as a multipurpose compound applicable in many industrial fields.

Results

Aspartate transaminase (AAT1) catalyzes the interconversion of aspartate and α-ketoglutarate to glutamate and oxaloacetate. We selected AAT1 in the oxalate metabolic pathway as a target of CRISPR/Cas9. Disruption of AAT1 leads to the accumulation of oxalate, rather than its conversion to α-ketoglutarate (AKG). Therefore, AAT1-mutant serves to lower the pH (pH 3–4) so as to increase the production of the pH-sensitive metabolite scleroglucan to 21.03 g L⁻¹ with a productivity of up to 0.25 g L⁻¹·h⁻¹.

Conclusions

We established a platform for gene editing that could rapidly generate and select mutants to provide a new beneficial strain of S. rolfsii as a scleroglucan hyper-producer, which is expected to reduce the cost of controlling the optimum pH condition in the fermentation industry.

Biocontrol effects of Penicillium griseofulvum against monkshood (Aconitum carmichaelii Debx.) root diseases caused by Sclerotium rolfsiii and Fusarium spp

AIMS

The aims of this study were to investigate the biocontrol effects of Penicillium griseofulvum strain CF3 and its mechanisms against soil-borne root pathogens (Fusarium oxysporum and Sclerotium rolfsii) of the medical plant Aconitum carmichaelii Debx.

METHODS AND RESULTS

The effects of P. griseofulvum strain CF3 were evaluated with regard to the hyphal growth of S. rolfsii and F. oxysporum, the sclerotial formation and germination of S. rolfsii and its expression of sclerotia-formation-related genes. A field experiment was conducted to explore how strain CF3 controls the severity of soil-borne diseases, promotes the growth of A. carmichaelii plants and mediates shifts in the culturable rhizosphere microbial populations. The results showed that treatment with a cell-free culture filtrate of strain CF3 considerably inhibited the hyphal growth of both S. rolfsii and F. oxysporum, in addition to limiting the sclerotial formation and germination of S. rolfsii. Three genes related to sclerotial formation (ArsclR, ArnsdD1 and ArnsdD2) were predicted in S. rolfsii and their expression was found suppressed by the CF3 treatment. Field application of the CF3 biocontrol agent in a powder form (1·9 × 10¹⁰ conidia per gram of substrate) reduced soil-borne disease severity by 15·0%. The shoot and root growth of A. carmichaelii plants was promoted by 61·6 and 83·1% respectively, as the biocontrol strain massively colonized the rhizosphere soil. The CF3 treatment also markedly reduced the density of some known species harmful to plants while increasing the density of some beneficial species in the rhizosphere soil.

SIGNIFICANCE AND IMPACT OF THE STUDY

Genes related to sclerotia formation of S. rolfsii are predicted for the first time and their expression patterns in the presence of P. griseofulvum strain CF3 are evaluated. This comprehensive study provides a candidate fungal biocontrol strain and reveals its potential mechanisms against S. rolfsii and F. oxysporum in A. carmichaelii plants.