Transcriptomic and Metabolomic Analyses Provide Insights into the Pathogenic Mechanism of the Rice False Smut Pathogen Ustilaginoidea virens.

Rice false smut, caused by the fungal pathogen Ustilaginoidea virens, is a worldwide rice fungal disease. However, the molecular mechanism of the pathogenicity of the fungus U. virens remains unclear. To understand the molecular mechanism of pathogenesis of the fungus U. virens, we performed an integrated analysis of the transcriptome and metabolome of strongly (S) and weakly (W) virulent strains both before and after the infection of panicles. A total of 7932 differential expressed genes (DEGs) were identified using transcriptome analysis. Gene ontology (GO) and metabolic pathway enrichment analysis indicated that amino acid metabolism, autophagy-yeast, MAPK signaling pathway-yeast, and starch and sucrose metabolism were closely related to the pathogenicity of U. virens. Genes related to pathogenicity were significantly upregulated in the strongly virulent strain, and were ATG, MAPK, STE, TPS, and NTH genes. However, genes involved in the negative regulation of pathogenesis were significantly downregulated and contained TOR kinase, TORC1, and autophagy-related protein genes. Metabolome analysis identified 698 differentially accumulated metabolites (DAMs), including 13 categories of organic acids and derivatives, lipids and lipid-like molecules, organoheterocyclic compounds. The significantly enriched pathways of DAMs mainly included amino acids and carbohydrates, and they accumulated after infection by the S strain. To understand the relevance of DEGs and DAMs in the pathogenicity of U. virens, transcriptomic and metabolomic data were integrated and analyzed. These results further confirmed that the pathogenesis of U. virens was regulated by DEGs and DAMs related to these four pathways, involving arginine and proline metabolism, lysine biosynthesis, alanine, aspartate and glutamate metabolism, and starch and sugar metabolism. Therefore, we speculate that the pathogenicity of U. virens is closely related to the accumulation of amino acids and carbohydrates, and to the changes in the expression of related genes.

Conjunctive Analysis of BSA-Seq and SSR Markers Unveil the Candidate Genes for Resistance to Rice False Smut.

Rice false smut (RFS) caused by the fungus Ustilaginoidea virens (Cook) leads to serious yield losses in rice. Identification of the gene or quantitative trait loci (QTLs) is crucial to resistance breeding and mitigation of RFS damage. In this study, we crossed a resistant variety, IR77298-14-1-2::IRGC117374-1, with a susceptible indica cultivar, 9311, and evaluated recombinant inbred lines in a greenhouse. The genetic analysis showed that the RFS resistance of IR77298-14-1-2::IRGC117374-1 was controlled by multiple recessive loci. We identified a novel QTL, qRFS12.01, for RFS resistance in IR77298-14-1-2::IRGC117374-1 by combining bulked segregant analysis with whole genome resequencing (BSA-seq) and simple sequence repeat (SSR) marker mapping approaches. The phenotypic effect of qRFS12.01 on RFS resistance reached 28.74%, suggesting that SSR markers linked to qRFS12.01 are valuable for marker-assisted breeding of RFS resistance in rice. The prediction of putative candidate genes within qRFS12.01 revealed five disease resistance proteins containing NB-ARC domains. In conclusion, our findings provide a new rice chromosome region carrying genes/QTLs for resistance to RFS.

Transcription Profiling of Rice Panicle in Response to Crude Toxin Extract of Ustilaginoidea virens.

Ustilaginoidea virens infects rice, causing rice false smut disease and reduced yields. During its growth, U. virens can also produce some toxins but less is known about the response mechanisms of the plant to U. virens toxins. U. virens toxins can inhibit the accumulation of total sugar in rice panicles. We used RNA sequencing to analyze the differential expression profile induced by infiltrating crude toxins into early growth-stage rice panicles. We compared the transcriptomes of the control and crude toxin-treated rice panicles and determined variable transcriptional responses under the action of the crude toxins. A total of 6,127 differentially expressed genes (DEGs) were identified. Among these genes, 3,150 were upregulated and 2,977 were downregulated. Gene Ontology (GO) and metabolic pathway enrichment analyses indicated that U. virens toxins mainly influenced glycometabolism, amino acid metabolism, and secondary metabolism of rice panicles. DEG analysis showed that the gene expression levels of 10 transcription factor families were significantly changed. Genes involved in phenylpropanoid biosynthesis, flavonoid biosynthesis, sugar transporters, and starch synthesis-related were significantly downregulated, including cytochrome P450, beta-glucosidase, CHS1, sucrose transporters, SWEETs, starch-branching enzymes, and UDP-glucose pyrophosphorylase. However, genes involved in programmed cell death (PCD) were significantly upregulated and contained cytochrome c, metacaspase, and protein kinase genes. The results indicate that U. virens toxins may act as the pathogenic factors to reduce stress resistance, disrupt total sugar accumulation and starch formation, and induce PCD.

Diversity Analysis of the Rice False Smut Pathogen Ustilaginoidea virens in Southwest China.

Rice false smut caused by Ustilaginoidea virens is a destructive disease in rice cropping areas of the world. The present study is focused on the morphology, pathogenicity, mating-type loci distribution, and genetic characterization of different isolates of U. virens. A total of 221 strains of U. virens were collected from 13 rice-growing regions in southwest China. The morphological features of these strains exhibited high diversity, and the pathogenicity of the smut fungus showed significant differentiation. There was no correlation between pathogenicity and sporulation. Mating-type locus (MAT) analysis revealed that all 221 isolates comprised heterothallic and homothallic forms, wherein 204 (92.31%) and 17 (7.69%) isolates belonged to heterothallic and homothallic mating types, respectively. Among 204 strains of heterothallic mating types, 62 (28.05%) contained MAT1-1-1 idiomorphs, and 142 isolates (64.25%) had the MAT1-2-1 idiomorph. Interestingly, strains isolated from the same fungus ball had different mating types. The genetic structure of the isolates was analyzed using simple sequence repeats (SSRs) and single-nucleotide polymorphisms (SNPs). All isolates were clustered into five genetic groups. The values of Nei's gene diversity (H) and Shannon's information index (I) indicated that all strains as a group had higher genetic diversity than strains from a single geographical population. The pairwise population fixation index (FST) values also indicated significant genetic differentiation among all compared geographical populations. The analysis of molecular variation (AMOVA) indicated greater genetic variation within individual populations and less genetic variation among populations. The results showed that most of the strains were not clustered according to their geographical origin, showing the rich genetic diversity and the complex and diverse genetic background of U. virens in southwest China. These results should help to better understand the biological and genetic diversity of U. virens in southwest China and provide a theoretical basis for building effective management strategies.

Genome sequence and transcriptome profiles of pathogenic fungus Paecilomyces penicillatus reveal its interactions with edible fungus Morchella importuna.

Paecilomyces penicillatus is one of the pathogens of morels, which greatly affects the yield and quality of Morchella spp.. In the present study, we de novo assembled the genome sequence of the fungus P. penicillatus SAAS_ppe1. We analyzed the transcriptional profile of P. penicillatus SAAS_ppe1 infection of Morchella importuna at different stages (3 days and 6 days after infection) and the response of M. importuna using the transcriptome. The assembled genome sequence of P. penicillatus SAAS_ppe1 was 39.78 Mb in length (11 scaffolds; scaffold N50, 6.50 Mb), in which 99.7% of the expected genes were detected. A total of 7.48% and 19.83% clean transcriptional reads from the infected sites were mapped to the P. penicillatus genome at the early and late stages of infection, respectively. There were 3,943 genes differently expressed in P. penicillatus at different stages of infection, of which 24 genes had increased expression with the infection and infection stage, including diphthamide biosynthesis, aldehyde reductase, and NAD (P)H-hydrate epimerase (P < 0.05). Several genes had variable expression trends at different stages of infection, indicating P. penicillatus had diverse regulation patterns to infect M. importuna. GO function, involving cellular components, and KEGG pathways, involving glycerolipid metabolism, and plant-pathogen interaction were significantly enriched during infection by P. penicillatus. The expression of ten genes in M. importuna increased during the infection and infection stage, and these may regulate the response of M. importuna to P. penicillatus infection. This is the first comprehensive study on P. penicillatus infection mechanism and M. importuna response mechanism, which will lay a foundation for understanding the fungus-fungus interactions, gene functions, and variety breeding of pathogenic and edible fungi.

Comparative mitochondrial genome analysis reveals intron dynamics and gene rearrangements in two Trametes species.

Trametes species are efficient wood decomposers that are widespread throughout the world. Mitogenomes have been widely used to understand the phylogeny and evolution of fungi. Up to now, two mitogenomes from the Trametes genus have been revealed. In the present study, the complete mitogenomes of two novel Trametes species, Trametes versicolor and T. coccinea, were assembled and compared with other Polyporales mitogenomes. Both species contained circular DNA molecules, with sizes of 67,318 bp and 99,976 bp, respectively. Comparative mitogenomic analysis indicated that the gene number, length and base composition varied between the four Trametes mitogenomes we tested. In addition, all of the core protein coding genes in Trametes species were identified and subjected to purifying selection. The mitogenome of T. coccinea contained the largest number of introns among the four Trametes species tested, and introns were considered the main factors contributing to size variations of Polyporales. Several novel introns were detected in the Trametes species we assembled, and introns identified in Polyporales were found to undergo frequent loss/gain events. Large-scale gene rearrangements were detected between closely related Trametes species, including gene inversions, insertions, and migrations. A well-supported phylogenetic tree for 77 Basidiomycetes was obtained based on the combined mitochondrial gene set using 2 phylogenetic inference methods. The results showed that mitochondrial genes are effective molecular markers for understanding the phylogeny of Basidiomycetes. This study is the first to report the mitogenome rearrangement and intron dynamics of Trametes species, which shed light on the evolution of Trametes and other related species.

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny.

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

The complete mitochondrial genome of Cladobotryum mycophilum (Hypocreales: Sordariomycetes).

Cladobotryum mycophilum is the causal agent of cobweb disease in many important mushroom crops. In this study, we report the complete mitochondrial genome of C. mycophilum for the first time. The genome is 78,729 bp long and comprises 52 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNA), and 26 transfer RNA (tRNA) genes. The nucleotide composition of C. mycophilum mitochondrial genome is as follows: A (38.06%), T (34.68%), C (12.19%), and G (15.07%). Phylogenetic analysis revealed that C. mycophilum had a close relationship with Cladobotryum varium from Hypocreaceae. This study provided a basis for studies of the mitochondrial evolution of Hypocreaceae.

Characterization of the complete mitochondrial genome of Corynespora cassiicola (Pleosporales: Dothideomycetes), with its phylogenetic analysis.

Corynespora cassiicola is a well-known plant pathogen with a broad host range and diverse lifestyles. In this study, we presented the complete mitochondrial genome (mitogenome) of C. cassiicola for the first time. It has a total length of 40,752 bp, which encodes 17 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNA), and 27 transfer RNA (tRNA) genes. The nucleotide composition of the mitogenome is: A (36.24%), T (34.62%), G (15.74%), and C (13.39%). Phylogenetic analysis revealed that C. cassiicola has a close relationship with Didymella pinodes from Didymellaceae.

Characterization of the mitochondrial genome of the pathogenic fungus Scytalidium auriculariicola (Leotiomycetes) and insights into its phylogenetics.

Scytalidium auriculariicola is the causative pathogen of slippery scar disease in the cultivated cloud ear fungus, Auricularia polytricha. In the present study, the mitogenome of S. auriculariicola was sequenced and assembled by next-generation sequencing technology. The circular mitogenome is 96,857 bp long and contains 56 protein-coding genes, 2 ribosomal RNA genes, and 30 transfer RNA genes (tRNAs). The high frequency of A and T used in codons contributed to the high AT content (73.70%) of the S. auriculariicola mitogenome. Comparative analysis indicated that the base composition and the number of introns and protein-coding genes in the S. auriculariicola mitogenome varied from that of other Leotiomycetes mitogenomes, including a uniquely positive AT skew. Five distinct groups were found in the gene arrangements of Leotiomycetes. Phylogenetic analyses based on combined gene datasets (15 protein-coding genes) yielded well-supported (BPP = 1) topologies. A single-gene phylogenetic tree indicated that the nad4 gene may be useful as a molecular marker to analyze the phylogenetic relationships of Leotiomycetes species. This study is the first report on the mitochondrial genome of the genus Scytalidium, and it will contribute to our understanding of the population genetics and evolution of S. auriculariicola and related species.

Transcriptional profiling provides new insights into the role of nitric oxide in enhancing Ganoderma oregonense resistance to heat stress.

Ganoderma is well known for its use in traditional Chinese medicine and is widely cultivated in China, Korea, and Japan. Increased temperatures associated with global warming are negatively influencing the growth and development of Ganoderma. Nitric oxide is reported to play an important role in alleviating fungal heat stress (HS). However, the transcriptional profiling of Ganoderma oregonense in response to HS, as well as the transcriptional response regulated by NO to cope with HS has not been reported. We used RNA-Seq technology to generate large-scale transcriptome data from G. oregonense mycelia subjected to HS (32 °C) and exposed to concentrations of exogenous NO. The results showed that heat shock proteins (HSPs), "probable stress-induced proteins", and unigenes involved in "D-amino-acid oxidase activity" and "oxidoreductase activity" were significantly up-regulated in G. oregonense subjected to HS (P < 0.05). The significantly up-regulated HSPs, "monooxygenases", "alcohol dehydrogenase", and "FAD/NAD(P)-binding domain-containing proteins" (P < 0.05) regulated by exogenous NO may play important roles in the enhanced HS tolerance of G. oregonense. These results provide insights into the transcriptional response of G. oregonense to HS and the mechanism by which NO enhances the HS tolerance of fungi at the gene expression level.

The evolution and pathogenic mechanisms of the rice sheath blight pathogen.

Rhizoctonia solani is a major fungal pathogen of rice (Oryza sativa L.) that causes great yield losses in all rice-growing regions of the world. Here we report the draft genome sequence of the rice sheath blight disease pathogen, R. solani AG1 IA, assembled using next-generation Illumina Genome Analyser sequencing technologies. The genome encodes a large and diverse set of secreted proteins, enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, which probably reflect an exclusive necrotrophic lifestyle. We find few repetitive elements, a closer relationship to Agaricomycotina among Basidiomycetes, and expand protein domains and families. Among the 25 candidate pathogen effectors identified according to their functionality and evolution, we validate 3 that trigger crop defence responses; hence we reveal the exclusive expression patterns of the pathogenic determinants during host infection.

Morphological structure of propagules and electrophoretic karyotype analysis of false smut Villosiclava virens in rice.

The target pathogen Villosiclava virens (teleomorph: claviceps oryzae-sativae) was isolated from the infected rice, where it caused false smut. In our study, the forming processes of the chlamydospores, chlamydospore balls, conidiospores, and secondary conidiospores during the asexual reproduction were observed more precisely and in greater detail than previous descriptions. The microstructure of the infected rice kernel showed that the outer dense chlamydospores piled around the false smut balls grown on XBZ medium; moreover the sclerotia consisting of dense mycelium were found. The different morphology was observed across the different growing conditions. In addition, we observed the nuclear numbers of both the conidiospores and hyphae using 4',6-diamidino-2-phenylindole (DAPI) staining. Because the fungus has small chromosomes and the numbers were not previously known, we analyzed the electrophoretic karyotype using a pulsed field gel electrophoresis (PFGE) technique. The results showed that V. virens has at least 10 chromosomes ranging in size from 0.6 kb to 6 Mb. The V. virens genome size is estimated to be 23 Mb. Here, we report the morphological characteristics of the fungus and the process of asexual spores forming asexual propagules, along with the first analyze the molecular karyotype of V. virens. These results supply a foundation for further study of the pathogenicity and biology of this devastating pathogen.