Integrative transcriptome analysis discloses the molecular basis of a heterogeneous fungal phytopathogen complex, Rhizoctonia solani AG-1 subgroups

Molecular characterization revealed the role of thaumatin-like proteins of Rhizoctonia solani AG4-JY in inducing maize disease resistance

The gene family of thaumatin-like proteins (TLPs) plays a crucial role in the adaptation of organisms to environmental stresses. In recent years, fungal secreted proteins (SP) with inducing disease resistance activity in plants have emerged as important elicitors in the control of fungal diseases. Identifying SPs with inducing disease resistance activity and studying their mechanisms are crucial for controlling sheath blight. In the present study, 10 proteins containing the thaumatin-like domain were identified in strain AG4-JY of Rhizoctonia solani and eight of the 10 proteins had signal peptides. Analysis of the TLP genes of the 10 different anastomosis groups (AGs) showed that the evolutionary relationship of the TLP gene was consistent with that between different AGs of R. solani. Furthermore, it was found that RsTLP3, RsTLP9 and RsTLP10 were regarded as secreted proteins for their signaling peptides exhibited secretory activity. Prokaryotic expression and enzyme activity analysis revealed that the three secreted proteins possess glycoside hydrolase activity, suggesting they belong to the TLP family. Additionally, spraying the crude enzyme solution of the three TLP proteins could enhance maize resistance to sheath blight. Further analysis showed that genes associated with the salicylic acid and ethylene pathways were up-regulated following RsTLP3 application. The results indicated that RsTLP3 had a good application prospect in biological control.

Mitochondrial-targeting effector RsIA_CtaG/Cox11 in Rhizoctonia solani AG-1 IA has two functions: plant immunity suppression and cell death induction mediated by a rice cytochrome c oxidase subunit

Rhizoctonia solani AG-1 IA causes a necrotrophic rice disease and is a serious threat to rice production. To date, only a few effectors have been characterized in AG-1 IA. We previously identified RsIA_CtaG/Cox11 and showed that infiltration of the recombinant protein into rice leaves caused disease-like symptoms. In the present study, we further characterized the functionality of RsIA_CtaG/Cox11. RsIA_CtaG/Cox11 is an alternative transcript of cytochrome c oxidase copper chaperone Cox11 that starts from the second AUG codon, but contains a functional secretion signal peptide. RNA interference with RsIA_CtaG/Cox11 reduced the pathogenicity of AG-1 IA towards rice and Nicotiana benthamiana without affecting its fitness or mycelial morphology. Transient expression of the RsIA_CtaG/Cox11-GFP fusion protein demonstrated the localization of RsIA_CtaG/Cox11 to mitochondria. Agro-infiltration of RsIA_CtaG/Cox11 into N. benthamiana leaves inhibited cell death by BAX and INF1. In contrast to rice, agro-infiltration of RsIA_CtaG/Cox11 did not induce cell death in N. benthamiana. However, cell death was observed when it was coinfiltrated with Os_CoxVIIa, which encodes a subunit of cytochrome c oxidase. Os_CoxVIIa appeared to interact with RsIA_CtaG/Cox11. The cell death triggered by coexpression of RsIA_CtaG/Cox11 and Os_CoxVIIa is independent of the leucine-rich repeat receptor kinases BAK1/SOBIR1 and enhanced the susceptibility of N. benthamiana to AG-1 IA. Two of the three evolutionarily conserved cysteine residues at positions 25 and 126 of RsIA_CtaG/Cox11 were essential for its immunosuppressive activity, but not for cell death induction. This report suggests that RsIA_CtaG/Cox11 appears to have a dual role in immunosuppression and cell death induction during pathogenesis.

Global distribution, traditional and modern detection, diagnostic, and management approaches of Rhizoctonia solani associated with legume crops

Sustainable development relies heavily on a food system that is both safe and secure. Several approaches may lead to sustainability and food safety. An increase in the cultivation of legume crops is one of the approaches for enhancing agricultural viability and ensuring adequate food supply. Legumes may increase daily intake of fiber, folate, and protein as substitutes for meat and dairy. They are also crucial in various intercropping systems worldwide. However, legume production has been hampered by Rhizoctonia solani due to its destructive lifestyle. R. solani causes blights, damping off, and rotting diseases in legume crops. Our knowledge of the global distribution of R. solani associated with legume crops (alfalfa, soybean, chickpea, pea, lentil, common bean, and peanut), detection, diagnosis, and management of legume crops diseases caused by R. solani is limited. Traditional approaches rely on the incubation of R. solani, visual examination of symptoms on host legume crops, and microscopy identification. However, these approaches are time-consuming, require technical expertise, fail to detect a minimal amount of inoculum, and are unreliable. Biochemical and molecular-based approaches have been used with great success recently because of their excellent sensitivity and specificity. Along with conventional PCR, nested PCR, multiplex PCR, real-time PCR, magnetic-capture hybridization PCR, and loop-mediated isothermal amplification have been widely used to detect and diagnose R. solani. In the future, Next-generation sequencing will likely be used to a greater extent to detect R. solani. This review outlines global distribution, survival, infection and disease cycle, traditional, biochemical, molecular, and next-generation sequencing detection and diagnostic approaches, and an overview of the resistant resources and other management strategies to cope with R. solani.

Identification and characterization of pathogenicity-related genes of Rhizoctonia solani AG3 during tobacco infection

Tobacco target spot disease is caused by a ubiquitous soil-borne phytopathogen Rhizoctonia solani; the pathogenic mechanisms underlying the effects of R. solani remain unclear. Deeper understanding of the functional responses to R. solani during host plant infection would help identify the molecular mechanisms essential for successful host invasion. In this study, we performed global transcriptional analysis of R. solani during various stages (12, 24, 48, 72, 96, and 120 h) of tobacco infection via an RNA sequencing method, while utilizing the pathosystem model R. solani AG3-tobacco (Nicotiana tabacum L.). After R. solani inoculation, the number of differentially expressed genes of R. solani differed at the various time points. Moreover, several gene ontology and Kyoto encyclopedia of genes and genomes pathways were unique in different infection stages, especially with respect to the genes involved in plant cell wall degradation and catalysis of biotransformation reactions, such as the pectin metabolic process and pectin catabolic process. The overexpressing-PD8 N. benthamiana plants enhanced the susceptibility to R. solani. In addition, we found that large amounts of reactive oxygen species (ROS) were generated in tobacco after infected by R. solani. R. solani encoding FAD/NAD binding oxidoreductase and peroxidase gene family to eliminating ROS and counteract oxidative stress. Moreover, Perox3 was validated that can enhance the ability of scavenging ROS by co-injecting. Overall, our findings show that pectin-degrading enzymes and cytochrome P450 genes are critical for plant infection. These results provide comprehensive insights into R. solani AG3 transcriptome responses during tobacco invasion.

Identification of 3-Methoxyphenylacetic Acid as a Phytotoxin, Produced by Rhizoctonia solani AG-3 TB

Tobacco target spot disease is caused by Rhizoctonia solani AG-3 TB, which causes serious harm to the quality and yield of tobacco. In this study, thin layer chromatography (TLC), high performance liquid chromatography (HPLC), infrared absorption spectroscopy (IR), and nuclear magnetic resonance spectroscopy (NMR) were used to purify and identify the potential phytotoxin produced by R. solani AG-3 TB. The result indicated that the purified toxin compound was 3-methoxyphenylacetic acid (3-MOPAA) (molecular formula: C₉H₁₀O₃). The exogenous purified compound 3-MOPAA was tested, and the results revealed that 3-MOPAA can cause necrosis in tobacco leaves. 3-MOPAA is a derivative of phenylacetic acid (PAA), which should be produced by specific enzymes, such as hydroxylase or methylase, in the presence of PAA. These results enrich the research on the pathogenic phytotoxins of R. solani and provide valuable insights into the pathogenic mechanism of AG-3 TB.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

Integrative transcriptome analysis revealed the pathogenic molecular basis of Rhizoctonia solani AG-3 TB at three progressive stages of infection

Rhizoctonia solani has a broad host range and results in significant losses in agricultural production. Here, an integrated transcriptomic analysis was performed to reveal the critical genes responsible for the pathogenesis of R. solani AG-3 TB on Nicotiana tabacum at different infection stages. The results showed that various differential expressed genes (DEGs) were enriched in fatty acid metabolism, amino sugar, carbon metabolism, and cellular carbohydrate biosynthetic process at the early (6–12 hpi), middle (24–36 hpi), and late stage (48–72 hpi) of infection. Specifically, several critical genes such as shikimate kinase that were involved in the biosynthesis of an important fungal toxin, phenylacetic acid (PAA) showed markedly increase at 24 hpi. Additionally, the genes expression levels of carbohydrate-active enzymes (CAZymes) and cell wall degrading enzymes (CWDEs) were significantly increased at the late infection stage. Furthermore, we identified 807 potential secreted proteins and 78 small cysteine-rich proteins, which may function as fungal effectors and involved in the pathogenicity. These results provide valuable insights into critical and potential genes as well as the pathways involved in the pathogenesis of R. solani AG-3 TB.

Secreted Glycosyltransferase RsIA_GT of Rhizoctonia solani AG-1 IA Inhibits Defense Responses in Nicotiana benthamiana

Anastomosis group AG-1 IA of Rhizoctonia solani Khün has a wide host range and threatens crop production. Various glycosyltransferases secreted by phytopathogenic fungi play an essential role in pathogenicity. Previously, we identified a glycosyltransferase RsIA_GT (AG11A_09161) as a secreted protein-encoding gene of R. solani AG-1 IA, whose expression levels increased during infection in rice. In this study, we further characterized the virulence function of RsIA_GT. It is conserved not only in Basidiomycota, including multiple anastomosis groups of R. solani, but also in other primary fungal taxonomic categories. RsIA_GT possesses a signal peptide (SP) for protein secretion, and its functionality was proven using yeast and Nicotiana benthamiana. The SP-truncated form of RsIA_GT (RsIA_GT(ΔS)) expressed in Escherichia coli-induced lesion-like phenotype in rice leaves when applied to punched leaves. However, Agrobacterium-mediated transient expressions of both the full-length RsIA_GT and RsIA_GT(ΔS) did not induce cell death in N. benthamiana leaves. Instead, only RsIA_GT(ΔS) suppressed the cell death induced by two reference cell death factors BAX and INF1 in N.benthamiana. RsIA_GT(ΔS)^{R154A D168A D170A}, a mutant RsIA_GT(ΔS) for the glycosyltransferase catalytic domain, still suppressed the BAX- or INF1-induced cell death, suggesting that the cell death suppression activity of RsIA_GT(ΔS) would be independent from its enzymatic activity. RsIA_GT(ΔS) also suppressed the H₂O₂ production and callose deposition and showed an effect on the induction of defense genes associated with the expression of BAX and INF1. The transient expression of RsIA_GT(ΔS) in N. benthamiana enhanced the lesion area caused by R. solani AG-1 IA. The secreted glycosyltransferase, RsIA_GT, of R. solani AG-1 IA is likely to have a dual role in virulence inside and outside of host cells.

Pangenome Analysis of the Soilborne Fungal Phytopathogen Rhizoctonia solani and Development of a Comprehensive Web Resource: RsolaniDB

Rhizoctonia solani is a collective group of genetically and pathologically diverse basidiomycetous fungi that damage economically important crops. Its isolates are classified into 13 Anastomosis Groups (AGs) and subgroups having distinctive morphology and host ranges. The genetic factors driving the unique features of R. solani pathology are not well characterized due to the limited availability of its annotated genomes. Therefore, we performed genome sequencing, assembly, annotation and functional analysis of 12 R. solani isolates covering 7 AGs and select subgroups (AG1-IA; AG1-IB; AG1-IC; AG2-2IIIB; AG3-PT, isolates Rhs 1AP and the hypovirulent Rhs1A1; AG3-TB; AG4-HG-I, isolates Rs23 and R118-11; AG5; AG6; and AG8), in which six genomes are reported for the first time. Using a pangenome comparative analysis of 12 R. solani isolates and 15 other Basidiomycetes, we defined the unique and shared secretomes, CAZymes, and effectors across the AGs. We have also elucidated the R. solani-derived factors potentially involved in determining AG-specific host preference, and the attributes distinguishing them from other Basidiomycetes. Finally, we present the largest repertoire of R. solani genomes and their annotated components as a comprehensive database, viz. RsolaniDB, with tools for large-scale data mining, functional enrichment and sequence analysis not available with other state-of-the-art platforms.

Contemporary proteomic research on lignocellulosic enzymes and enzymolysis: A review

This review overviewed the current researches on the isolation of novel strains, the development of novel identification protocols, the key enzymes and their synergistic interactions with other functional enzyme systems, and the strategies for enhancing enzymolysis efficiencies. The main obstacle for realizing biorefinery of lignocellulosic biomass to biofuels or biochemicals is the high cost of enzymolysis stage. Therefore, research prospects to reduce the costs for lignocellulose hydrolysis were outlined.

Functional validation of pathogenicity genes in rice sheath blight pathogen Rhizoctonia solani by a novel host-induced gene silencing system

Rice sheath blight, caused by the soilborne fungus Rhizoctonia solani, causes severe yield losses worldwide. Elucidation of the pathogenic mechanism of R. solani is highly desired. However, the lack of a stable genetic transformation system has made it challenging to examine genes' functions in this fungus. Here, we present functional validation of pathogenicity genes in the rice sheath blight pathogen R. solani by a newly established tobacco rattle virus (TRV)-host-induced gene silencing (HIGS) system using the virulent R. solani AG-1 IA strain GD-118. RNA interference constructs of 33 candidate pathogenicity genes were infiltrated into Nicotiana benthamiana leaves with the TRV-HIGS system. Of these constructs, 29 resulted in a significant reduction in necrosis caused by GD-118 infection. For further validation of one of the positive genes, trehalose-6-phosphate phosphatase (Rstps2), stable rice transformants harbouring the double-stranded RNA (dsRNA) construct for Rstps2 were created. The transformants exhibited reduced gene expression of Rstps2, virulence, and trehalose accumulation in GD-118. We showed that the dsRNA for Rstps2 was taken up by GD-118 mycelia and sclerotial differentiation of GD-118 was inhibited. These findings offer gene identification opportunities for the rice sheath blight pathogen and a theoretical basis for controlling this disease by spray-induced gene silencing.

Comparative Mitogenomic Analysis and the Evolution of Rhizoctonia solani Anastomosis Groups

Mitochondria are the major energy source for cell functions. However, for the plant fungal pathogens, mitogenome variations and their roles during the host infection processes remain largely unknown. Rhizoctonia solani, an important soil-borne pathogen, forms different anastomosis groups (AGs) and adapts to a broad range of hosts in nature. Here, we reported three complete mitogenomes of AG1-IA RSIA1, AG1-IB RSIB1, and AG1-IC, and performed a comparative analysis with nine published Rhizoctonia mitogenomes (AG1-IA XN, AG1-IB 7/3/14, AG3, AG4, and five Rhizoctonia sp. mitogenomes). These mitogenomes encoded 15 typical proteins (cox1-3, cob, atp6, atp8-9, nad1-6, nad4L, and rps3) and several LAGLIDADG/GIY-YIG endonucleases with sizes ranging from 109,017 bp (Rhizoctonia sp. SM) to 235,849 bp (AG3). We found that their large sizes were mainly contributed by repeat sequences and genes encoding endonucleases. We identified the complete sequence of the rps3 gene in 10 Rhizoctonia mitogenomes, which contained 14 positively selected sites. Moreover, we inferred a robust maximum-likelihood phylogeny of 32 Basidiomycota mitogenomes, representing that seven R. solani and other five Rhizoctonia sp. lineages formed two parallel branches in Agaricomycotina. The comparative analysis showed that mitogenomes of Basidiomycota pathogens had high GC content and mitogenomes of R. solani had high repeat content. Compared to other strains, the AG1-IC strain had low substitution rates, which may affect its mitochondrial phylogenetic placement in the R. solani clade. Additionally, with the published RNA-seq data, we investigated gene expression patterns from different AGs during host infection stages. The expressed genes from AG1-IA (host: rice) and AG3 (host: potato) mainly formed four groups by k-mean partitioning analysis. However, conserved genes represented varied expression patterns, and only the patterns of rps3-nad2 and nad1-m3g18/mag28 (an LAGLIDADG endonuclease) were conserved in AG1-IA and AG3 as shown by the correlation coefficient analysis, suggesting regulation of gene repertoires adapting to infect varied hosts. The results of variations in mitogenome characteristics and the gene substitution rates and expression patterns may provide insights into the evolution of R. solani mitogenomes.

Genome Analyses of the Less Aggressive Rhizoctonia solani AG1-IB Isolates 1/2/21 and O8/2 Compared to the Reference AG1-IB Isolate 7/3/14

Rhizoctonia solani AG1-IB of the phylum Basidiomycota is known as phytopathogenic fungus affecting various economically important crops, such as bean, rice, soybean, figs, cabbage and lettuce. The isolates 1/2/21 and O8/2 of the anastomosis group AG1-IB originating from lettuce plants with bottom rot symptoms represent two less aggressive R. solani isolates, as confirmed in a pathogenicity test on lettuce. They were deeply sequenced on the Illumina MiSeq system applying the mate-pair and paired-end mode to establish their genome sequences. Assemblies of obtained sequences resulted in 2092 and 1492 scaffolds, respectively, for isolate 1/2/21 and O8/2, amounting to a size of approximately 43 Mb for each isolate. Gene prediction by applying AUGUSTUS (v. 3.2.1.) yielded 12,827 and 12,973 identified genes, respectively. Based on automatic functional annotation, genes potentially encoding cellulases and enzymes involved in secondary metabolite synthesis were identified in the AG1-IB genomes. The annotated genome sequences of the less aggressive AG1-IB isolates were compared with the isolate 7/3/14, which is highly aggressive on lettuce and other vegetable crops such as bean, cabbage and carrot. This analysis revealed the first insights into core genes of AG1-IB isolates and unique determinants of each genome that may explain the different aggressiveness levels of the strains.

Comparative Genomics: Insights on the Pathogenicity and Lifestyle of Rhizoctonia solani

Proper management of agricultural disease is important to ensure sustainable food security. Staple food crops like rice, wheat, cereals, and other cash crops hold great export value for countries. Ensuring proper supply is critical; hence any biotic or abiotic factors contributing to the shortfall in yield of these crops should be alleviated. Rhizoctonia solani is a major biotic factor that results in yield losses in many agriculturally important crops. This paper focuses on genome informatics of our Malaysian Draft R. solani AG1-IA, and the comparative genomics (inter- and intra- AG) with four AGs including China AG1-IA (AG1-IA_KB317705.1), AG1-IB, AG3, and AG8. The genomic content of repeat elements, transposable elements (TEs), syntenic genomic blocks, functions of protein-coding genes as well as core orthologous genic information that underlies R. solani’s pathogenicity strategy were investigated. Our analyses show that all studied AGs have low content and varying profiles of TEs. All AGs were dominant for Class I TE, much like other basidiomycete pathogens. All AGs demonstrate dominance in Glycoside Hydrolase protein-coding gene assignments suggesting its importance in infiltration and infection of host. Our profiling also provides a basis for further investigation on lack of correlation observed between number of pathogenicity and enzyme-related genes with host range. Despite being grouped within the same AG with China AG1-IA, our Draft AG1-IA exhibits differences in terms of protein-coding gene proportions and classifications. This implies that strains from similar AG do not necessarily have to retain similar proportions and classification of TE but must have the necessary arsenal to enable successful infiltration and colonization of host. In a larger perspective, all the studied AGs essentially share core genes that are generally involved in adhesion, penetration, and host colonization. However, the different infiltration strategies will depend on the level of host resilience where this is clearly exhibited by the gene sets encoded for the process of infiltration, infection, and protection from host.

Do Lytic Polysaccharide Monooxygenases Aid in Plant Pathogenesis and Herbivory?

Lytic polysaccharide monooxygenases (LPMOs), copper-dependent enzymes mainly found in fungi, bacteria, and viruses, are responsible for enabling plant infection and degradation processes. Since their discovery 10 years ago, significant progress has been made in understanding the major role these enzymes play in biomass conversion. The recent discovery of additional LPMO families in fungi and oomycetes (AA16) as well as insects (AA15) strongly suggests that LPMOs might also be involved in biological processes such as overcoming plant defenses. In this review, we aim to give a comprehensive overview of the potential role of different LPMO families from the perspective of plant defense and their multiple implications in devising new strategies for achieving crop protection from plant pathogens and insect pests.

Identification of effector candidate genes of Rhizoctonia solani AG-1 IA expressed during infection in Brachypodium distachyon

Rhizoctonia solani is a necrotrophic phytopathogen belonging to basidiomycetes. It causes rice sheath blight which inflicts serious damage in rice production. The infection strategy of this pathogen remains unclear. We previously demonstrated that salicylic acid-induced immunity could block R. solani AG-1 IA infection in both rice and Brachypodium distachyon. R. solani may undergo biotrophic process using effector proteins to suppress host immunity before necrotrophic stage. To identify pathogen genes expressed at the early infection process, here we developed an inoculation method using B. distachyon which enables to sample an increased amount of semi-synchronous infection hyphae. Sixty-one R. solani secretory effector-like protein genes (RsSEPGs) were identified using in silico approach with the publicly available gene annotation of R. solani AG-1 IA genome and our RNA-sequencing results obtained from hyphae grown on agar medium. Expression of RsSEPGs was analyzed at 6, 10, 16, 24, and 32 h after inoculation by a quantitative reverse transcription-polymerase chain reaction and 52 genes could be detected at least on a single time point tested. Their expressions showed phase-specific patterns which were classified into 6 clusters. The 23 RsSEPGs in the cluster 1–3 and 29 RsSEPGs in the cluster 4–6 are expected to be involved in biotrophic and necrotrophic interactions, respectively.