Comparative genome analyses of four rice-infecting Rhizoctonia solani isolates reveal extensive enrichment of homogalacturonan modification genes

The genome sequencing and comparative genomics analysis of Rhizoctonia solani reveals a novel effector family owning a uinque domain in Basidiomycetes

Rhizoctonia solani is a soil-borne pathogen with 14 anastomosis groups (AGs), and different subgroups are genetically diverse. However, the genetic factors contributing to the pathogenicity of the fungus have not been well characterized. In this study, the genome of R. solani AG1-ZJ was sequenced. As the result, a 41.57 Mb draft genome containing 12,197 putative coding genes was obtained. Comparative genomic analysis of 11 different AGs revealed conservation and unique characteristics between the AGs. Furthermore, a novel effector family containing a 68 amino acid conserved domain unique in basidiomycetous fungi was characterized. Two effectors containing the conserved domain in AG4-JY were identified, and named as RsUEB1 and RsUEB2. Furthermore, the spray-induced gene silencing strategy was used to generate a dsRNA capable of silencing the conserved domain sequence of RsUEB1 and RsUEB2. This dsRNA can significantly reduce the expression of RsUEB1 and RsUEB2 and the pathogenicity of AG4-JY on foxtail millet, maize, rice and wheat. In conclusion, this study provides significant insights into the pathogenicity mechanisms of R. solani. The identification of the conserved domain and the successful use of dsRNA silencing of the gene containing the conserved domain will offer a new strategy for controlling sheath blight in cereal crops.

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.

Heritability and gene functions associated with sclerotia formation of Rhizoctonia solani AG-7 using whole genome sequencing and genome-wide association study

Sclerotia are specialized fungal structures formed by pigmented and aggregated hyphae, which can survive under unfavourable environmental conditions and serve as the primary inocula for several phytopathogenic fungi including Rhizoctonia solani. Among 154 R. solani anastomosis group 7 (AG-7) isolates collected in fields, the sclerotia-forming capability regarding sclerotia number and sclerotia size varied in the fungal population, but the genetic makeup of these phenotypes remained unclear. As limited studies have focused on the genomics of R. solani AG-7 and the population genetics of sclerotia formation, this study completed the whole genome sequencing and gene prediction of R. solani AG-7 using the Oxford NanoPore and Illumina RNA sequencing. Meanwhile, a high-throughput image-based method was established to quantify the sclerotia-forming capability, and the phenotypic correlation between sclerotia number and sclerotia size was low. A genome-wide association study identified three and five significant SNPs associated with sclerotia number and size in distinct genomic regions, respectively. Of these significant SNPs, two and four showed significant differences in the phenotypic mean separation for sclerotia number and sclerotia size, respectively. Gene ontology enrichment analysis focusing on the linkage disequilibrium blocks of significant SNPs identified more categories related to oxidative stress for sclerotia number, and more categories related to cell development, signalling and metabolism for sclerotia size. These results indicated that different genetic mechanisms may underlie these two phenotypes. Moreover, the heritability of sclerotia number and sclerotia size were estimated for the first time to be 0.92 and 0.31, respectively. This study provides new insights into the heritability and gene functions related to the development of sclerotia number and sclerotia size, which could provide additional knowledge to reduce fungal residues in fields and achieve sustainable disease management.

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.

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.

Evolution of pathogenicity-associated genes in Rhizoctonia solani AG1-IA by genome duplication and transposon-mediated gene function alterations

BACKGROUND

Rhizoctonia solani is a polyphagous fungal pathogen that causes diseases in crops. The fungal strains are classified into anastomosis groups (AGs); however, genomic complexity, diversification into the AGs and the evolution of pathogenicity-associated genes remain poorly understood.

RESULTS

We report a recent whole-genome duplication and sequential segmental duplications in AG1-IA strains of R. solani. Transposable element (TE) clusters have caused loss of synteny in the duplicated blocks and introduced differential structural alterations in the functional domains of several pathogenicity-associated paralogous gene pairs. We demonstrate that the TE-mediated structural variations in a glycosyl hydrolase domain and a GMC oxidoreductase domain in two paralogous pairs affect the pathogenicity of R. solani. Furthermore, to investigate the association of TEs with the natural selection and evolution of pathogenicity, we sequenced the genomes of forty-two rice field isolates of R. solani AG1-IA. The genomic regions with high population mutation rates and with the lowest nucleotide diversity are enriched with TEs. Genetic diversity analysis predicted the genes that are most likely under diversifying and purifying selections. We present evidence that a smaller variant of a glucosamine phosphate N-acetyltransferase (GNAT) protein, predicted to be under purifying selection, and an LPMP_AA9 domain-containing protein, predicted to be under diversifying selection, are important for the successful pathogenesis of R. solani in rice as well as tomato.

CONCLUSIONS

Our study has unravelled whole-genome duplication, TE-mediated neofunctionalization of genes and evolution of pathogenicity traits in R. solani AG1-IA. The pathogenicity-associated genes identified during the study can serve as novel targets for disease control.

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.

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.