Evolutionary and genomic comparisons of hybrid uninucleate and nonhybrid Rhizoctonia fungi

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.

Intraspecific diploidization of a halophyte root fungus drives heterosis

How organisms respond to environmental stress is a key topic in evolutionary biology. This study focused on the genomic evolution of Laburnicola rhizohalophila, a dark-septate endophytic fungus from roots of a halophyte. Chromosome-level assemblies were generated from five representative isolates from structured subpopulations. The data revealed significant genomic plasticity resulting from chromosomal polymorphisms created by fusion and fission events, known as dysploidy. Analyses of genomic features, phylogenomics, and macrosynteny have provided clear evidence for the origin of intraspecific diploid-like hybrids. Notably, one diploid phenotype stood out as an outlier and exhibited a conditional fitness advantage when exposed to a range of abiotic stresses compared with its parents. By comparing the gene expression patterns in each hybrid parent triad under the four growth conditions, the mechanisms underlying growth vigor were corroborated through an analysis of transgressively upregulated genes enriched in membrane glycerolipid biosynthesis and transmembrane transporter activity. In vitro assays suggested increased membrane integrity and lipid accumulation, as well as decreased malondialdehyde production under optimal salt conditions (0.3 M NaCl) in the hybrid. These attributes have been implicated in salinity tolerance. This study supports the notion that hybridization-induced genome doubling leads to the emergence of phenotypic innovations in an extremophilic endophyte.

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.

Variety-driven rhizosphere microbiome bestows differential salt tolerance to alfalfa for coping with salinity stress

Soil salinization is a global environmental issue and a significant abiotic stress that threatens crop production. Root-associated rhizosphere microbiota play a pivotal role in enhancing plant tolerance to abiotic stresses. However, limited information is available concerning the specific variations in rhizosphere microbiota driven by different plant genotypes (varieties) in response to varying levels of salinity stress. In this study, we compared the growth performance of three alfalfa varieties with varying salt tolerance levels in soils with different degrees of salinization. High-throughput 16S rRNA and ITS sequencing were employed to analyze the rhizosphere microbial communities. Undoubtedly, the increasing salinity significantly inhibited alfalfa growth and reduced rhizosphere microbial diversity. However, intriguingly, salt-tolerant varieties exhibited relatively lower susceptibility to salinity, maintaining more stable rhizosphere bacterial community structure, whereas the reverse was observed for salt-sensitive varieties. Bacillus emerged as the dominant species in alfalfa's adaptation to salinity stress, constituting 21.20% of the shared bacterial genera among the three varieties. The higher abundance of Bacillus, Ensifer, and Pseudomonas in the rhizosphere of salt-tolerant alfalfa varieties is crucial in determining their elevated salt tolerance. As salinity levels increased, salt-sensitive varieties gradually accumulated a substantial population of pathogenic fungi, such as Fusarium and Rhizoctonia. Furthermore, rhizosphere bacteria of salt-tolerant varieties exhibited increased activity in various metabolic pathways, including biosynthesis of secondary metabolites, carbon metabolism, and biosynthesis of amino acids. It is suggested that salt-tolerant alfalfa varieties can provide more carbon sources to the rhizosphere, enriching more effective plant growth-promoting bacteria (PGPB) such as Pseudomonas to mitigate salinity stress. In conclusion, our results highlight the variety-mediated enrichment of rhizosphere microbiota in response to salinity stress, confirming that the high-abundance enrichment of specific dominant rhizosphere microbes and their vital roles play a significant role in conferring high salt adaptability to these varieties.

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.

Assessment of Genetic Variability and Evolutionary Relationships of Rhizoctonia solani Inherent in Legume Crops

Rhizoctonia solani is one of the most common soil-borne fungal pathogens of legume crops worldwide. We collected rDNA-ITS sequences from NCBI GenBank, and the aim of this study was to examine the genetic diversity and phylogenetic relationships of various R. solani anastomosis groups (AGs) that are commonly associated with grain legumes (such as soybean, common bean, pea, peanut, cowpea, and chickpea) and forage legumes (including alfalfa and clover). Soybean is recognized as a host for multiple AGs, with AG-1 and AG-2 being extensively investigated. This is evidenced by the higher representation of sequences associated with these AGs in the NCBI GenBank. Other AGs documented in soybean include AG-4, AG-7, AG-11, AG-5, AG-6, and AG-9. Moreover, AG-4 has been extensively studied concerning its occurrence in chickpea, pea, peanut, and alfalfa. Research on the common bean has been primarily focused on AG-2, AG-4, and AG-1. Similarly, AG-1 has been the subject of extensive investigation in clover and cowpea. Collectively, AG-1, AG-2, and AG-4 have consistently been identified and studied across these diverse legume crops. The phylogenetic analysis of R. solani isolates across different legumes indicates that the distinct clades or subclades formed by the isolates correspond to their specific anastomosis groups (AGs) and subgroups, rather than being determined by their host legume crop. Additionally, there is a high degree of sequence similarity among isolates within the same clade or subclade. Principal coordinate analysis (PCoA) further supports this finding, as isolates belonging to the same AGs and/or subgroups cluster together, irrespective of their host legume. Therefore, the observed clustering of R. solani AGs and subgroups without a direct association with the host legume crop provides additional support for the concept of AGs in understanding the genetic relationships and evolution of R. solani.

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.

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.

Phylogenomic analysis of the Neocallimastigomycota: proposal of Caecomycetaceae fam. nov., Piromycetaceae fam. nov., and emended description of the families Neocallimastigaceae and Anaeromycetaceae

The anaerobic gut fungi (AGF) represent a coherent phylogenetic clade within the Mycota. Twenty genera have been described so far. Currently, the phylogenetic and evolutionary relationships between AGF genera remain poorly understood. Here, we utilized 52 transcriptomic datasets from 14 genera to resolve AGF inter-genus relationships using phylogenomics, and to provide a quantitative estimate (amino acid identity, AAI) for intermediate rank assignments. We identify four distinct supra-genus clades, encompassing all genera producing polyflagellated zoospores, bulbous rhizoids, the broadly circumscribed genus Piromyces, and the Anaeromyces and affiliated genera. We also identify the genus Khoyollomyces as the earliest evolving AGF genus. Concordance between phylogenomic outputs and RPB1 and D1/D2 LSU, but not RPB2, MCM7, EF1α or ITS1, phylogenies was observed. We combine phylogenomic analysis and AAI outputs with informative phenotypic traits to propose accommodating 14/20 AGF genera into four families: Caecomycetaceae fam. nov. (encompassing the genera Caecomyces and Cyllamyces), Piromycetaceae fam. nov. (encompassing the genus Piromyces), emend the description of the family Neocallimastigaceae to encompass the genera Neocallimastix, Orpinomyces, Pecoramyces, Feramyces, Ghazallomyces, Aestipascuomyces and Paucimyces, as well as the family Anaeromycetaceae to include the genera Oontomyces, Liebetanzomyces and Capellomyces in addition to Anaeromyces. We refrain from proposing families for the deeply branching genus Khoyollomyces and for genera with uncertain position (Buwchfawromyces, Joblinomyces, Tahromyces, Agriosomyces and Aklioshbomyces) pending availability of additional isolates and sequence data; and these genera are designated as 'genera incertae sedis' in the order Neocallimastigales. Our results establish an evolutionary-grounded Linnaean taxonomic framework for the AGF, provide quantitative estimates for rank assignments, and demonstrate the utility of RPB1 as an additional informative marker in Neocallimastigomycota taxonomy.

Reduction of Rhizoctonia cerealis Infection on Wheat Through Host- and Spray-Induced Gene Silencing of an Orphan Secreted Gene

Rhizoctonia cerealis is a soilborne fungus that can cause sharp eyespot in wheat, resulting in massive yield losses found in many countries. Due to the lack of resistant cultivars, fungicides have been widely used to control this pathogen. However, chemical control is not environmentally friendly and is costly. Meanwhile, the lack of genetic transformation tools has hindered the functional characterization of virulence genes. In this study, we attempted to characterize the function of virulence genes by two transient methods, host-induced gene silencing (HIGS) and spray-induced gene silencing (SIGS), which use RNA interference to suppress the pathogenic development. We identified ten secretory orphan genes from the genome. After silencing these ten genes, only the RcOSP1 knocked-down plant significantly inhibited the growth of R. cerealis. We then described RcOSP1 as an effector that could impair wheat biological processes and suppress pathogen-associated molecular pattern-triggered immunity in the infection process. These findings confirm that HIGS and SIGS can be practical tools for researching R. cerealis virulence genes. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Morphopathological and Molecular Morphometric Characterization of Waitea circinata var. prodigus Causing a Novel Sheath Spot Disease of Maize in India

Maize brown sheath spot (MBSS), a new disease of maize, was discovered while surveying for maize leaf and sheath blight diseases in the Indian states of Assam, Jharkhand, Meghalaya, Manipur, and Odisha. Maize is the third most important cereal after rice and wheat in India. Unlike banded leaf and sheath blight disease caused by Rhizoctonia solani, MBSS symptoms on maize were discrete and limited to sheaths only. Symptoms of MBSS in the field were initially water-soaked necrotic lesions of 1 to 2 cm in diameter on the lowermost leaf sheaths, which then progressed to the upper sheaths. Lesions coalesced and covered approximately 2 to 5% of the sheath area. Infected dried lower leaves were shed, whereas infected upper leaves remained on the stem. The pathogen was isolated, characterized morphologically, pathologically, and molecularly, and identified as Waitea circinata var. prodigus, a basidiomycete known to cause basal leaf blight of seashore paspalum. The internal transcribed spacer (ITS) sequence 2 (ITS2) of rDNA from MBSS isolates formed a well supported clade with known W. circinata var. prodigus isolates. Molecular morphometric analysis of the ITS2 regions of the five known varieties of W. circinata detected distinguishing variations in GC content, compensatory base changes (CBCs), hemi- CBCs, indels, and altered base-pairing of helices. Variation in these characteristics may indicate that varieties are distinct biological species within W. circinata sensu lato. The geographical distribution and potential impacts of MBSS on the maize crop in India necessitate further investigations of pathogen identification and disease management.

Salicylic Acid Is Required for Broad-Spectrum Disease Resistance in Rice

Rice plants contain high basal levels of salicylic acid (SA), but some of their functions remain elusive. To elucidate the importance of SA homeostasis in rice immunity, we characterized four rice SA hydroxylase genes (OsSAHs) and verified their roles in SA metabolism and disease resistance. Recombinant OsSAH proteins catalyzed SA in vitro, while OsSAH3 protein showed only SA 5-hydroxylase (SA5H) activity, which was remarkably higher than that of other OsSAHs that presented both SA3H and SA5H activities. Amino acid substitutions revealed that three amino acids in the binding pocket affected SAH enzyme activity and/or specificity. Knockout OsSAH2 and OsSAH3 (sahKO) genes conferred enhanced resistance to both hemibiotrophic and necrotrophic pathogens, whereas overexpression of each OsSAH gene increased susceptibility to the pathogens. sahKO mutants showed increased SA and jasmonate levels compared to those of the wild type and OsSAH-overexpressing plants. Analysis of the OsSAH3 promoter indicated that its induction was mainly restricted around Magnaporthe oryzae infection sites. Taken together, our findings indicate that SA plays a vital role in immune signaling. Moreover, fine-tuning SA homeostasis through suppression of SA metabolism is an effective approach in studying broad-spectrum disease resistance in rice.

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.

Genomic Analysis of Sarcomyxa edulis Reveals the Basis of Its Medicinal Properties and Evolutionary Relationships

Yuanmo [Sarcomyxa edulis (Y.C. Dai, Niemelä & G.F. Qin) T. Saito, Tonouchi & T. Harada] is an important edible and medicinal mushroom endemic to Northeastern China. Here we report the de novo sequencing and assembly of the S. edulis genome using single-molecule real-time sequencing technology. The whole genome was approximately 35.65 Mb, with a G + C content of 48.31%. Genome assembly generated 41 contigs with an N50 length of 1,772,559 bp. The genome comprised 9,364 annotated protein-coding genes, many of which encoded enzymes involved in the modification, biosynthesis, and degradation of glycoconjugates and carbohydrates or enzymes predicted to be involved in the biosynthesis of secondary metabolites such as terpene, type I polyketide, siderophore, and fatty acids, which are responsible for the pharmacodynamic activities of S. edulis. We also identified genes encoding 1,3-β-glucan synthase and endo-1,3(4)-β-glucanase, which are involved in polysaccharide and uridine diphosphate glucose biosynthesis. Phylogenetic and comparative analyses of Basidiomycota fungi based on a single-copy orthologous protein indicated that the Sarcomyxa genus is an independent group that evolved from the Pleurotaceae family. The annotated whole-genome sequence of S. edulis can serve as a reference for investigations of bioactive compounds with medicinal value and the development and commercial production of superior S. edulis varieties.

Genome structure and diversity of novel endornaviruses from wheat sharp eyespot pathogen Rhizoctonia cerealis

Rhizoctonia cerealis (teleomorph Ceratobasidium cereale) is a soil-borne plant pathogenic fungus that can cause sharp eyespot in wheat or yellow patch in grasses. In this study, 21 new endornavirus genomes were obtained from five R. cerealis strains through the high-throughput sequencing of viral double-stranded RNA. Eighteen viruses were identified as Alphaendornavirus, and three viruses were identified as new species of Betaendornavirus on the basis of the phylogenetic analysis of the deduced amino acid sequences of RNA-dependent RNA polymerase. Notably, 12 of the new alphaendornaviruses could encode two open reading frames (ORFs), which were a rare feature of Endornaviridae. The amino acid sequences encoded by ORF2 from different endornaviruses had very low identity, and their functions and evolution origins remained unclear. Different endornavirus species with remarkably different genome structures could be found in the same R. cerealis strain. This study indicated that endornaviruses are common in R. cerealis and display wide diversity. Betaendornaviruses were found in R. cerealis, and a new species was proposed. This study is the first to report that the endornaviruses from R. cerealis can encode two ORFs and enhances our understanding of the viruses in the Endornaviridae family.