Search

The seed microbiota is an important component given by nature to plants, protecting seeds from damage by other organisms and abiotic stress. However, little is known about the dynamic changes and potential functions of the seed microbiota during seed development. In this study, we investigated the composition and potential functions of the seed microbiota of rapeseed (Brassica napus). A total of 2496 amplicon sequence variants (ASVs) belonging to 504 genera in 25 phyla were identified, and the seed microbiota of all sampling stages were divided into three groups. The microbiota of flower buds, young pods, and seeds at 20 days after flowering (daf) formed the first group; that of seeds at 30 daf, 40 daf and 50 daf formed the second group; that of mature seeds and parental seeds were clustered into the third group. The functions of seed microbiota were identified by using PICRUSt2, and it was found that the substance metabolism of seed microbiota was correlated with those of the seeds. Finally, sixty-one core ASVs, including several potential human pathogens, were identified, and a member of the seed core microbiota, Sphingomonas endophytica, was isolated from seeds and found to promote seedling growth and enhance resistance against Sclerotinia sclerotiorum, a major pathogen in rapeseed. Our findings provide a novel perspective for understanding the composition and functions of microbiota during seed development and may enhance the efficiency of mining beneficial seed microbes.

Phytopathogenic fungi normally secrete large amounts of CWDEs to enhance infection of plants. In this study, we identified and characterized a secreted glycosyl hydrolase 5 family member in Sclerotinia sclerotiorum (SsGH5, Sclerotinia sclerotiorum Glycosyl Hydrolase 5). SsGH5 was significantly upregulated during the early stages of infection. Knocking out SsGH5 did not affect the growth and acid production of S. sclerotiorum but resulted in decreased glucan utilization and significantly reduced virulence. In addition, Arabidopsis thaliana expressing SsGH5 became more susceptible to necrotrophic pathogens and basal immune responses were inhibited in these plants. Remarkably, the lost virulence of the ΔSsGH5 mutants was restored after inoculating onto SsGH5 transgenic Arabidopsis. In summary, these results highlight that S. sclerotiorum suppresses the immune responses of Arabidopsis through secreting SsGH5, and thus exerts full virulence for successful infection.

Sclerotinia sclerotiorum is a fungal pathogen with a broad range of hosts, which can cause diseases and pose a great threat to many crops. Fungal-specific Zn2Cys6 transcription factors (TFs) constitute a large family prevalent among plant pathogens. However, the function of Zn2Cys6 TFs remains largely unknown. In this study, we identified and characterized SsZNC1, a Zn2Cys6 TF in S. sclerotiorum, which is involved in virulence, sclerotial development, and osmotic stress response. The expression of SsZNC1 was significantly up-regulated in the early stages of S. sclerotiorum infection on Arabidopsis leaves. The target deletion of SsZNC1 resulted in reduced virulence on Arabidopsis and oilseed rape. In addition, sclerotial development ability and growth ability under hyperosmotic conditions of SsZNC1 knockout transformants were reduced. A transcriptomic analysis unveiled its regulatory role in key cellular functions, including cellulose catabolic process, methyltransferase activity, and virulence, etc. Together, our results indicated that SsZNC1, a core regulatory gene involved in virulence, sclerotial development and stress response, provides new insight into the transcription regulation and pathogenesis of S. sclerotiorum.

The concentration of selenium (Se) in agricultural products primarily depends on the concentration of Se in soil and the ability of plants to accumulate Se. Selenium deficiency not only leads to decreased body resistance, but also increases the risk of cancer. The form and concentration of bioavailable Se is important for diet. The present study was carried out via field experiment with wheat and broad beans in soil of different Se concentrations (0, 1.12, and 11.2 kg·ha−1), which was determined based on the national standard and the team’s previous experience. Results indicated that the concentration of Se in the edible organs of wheat was higher than in broad bean, while the enriched Se concentration in the root of broad bean was more than twice and three times higher than that of wheat at medium and high levels of Se, respectively. Selenomethionine, which accounted for over half of the total Se speciations, was the dominant species in the edible parts of the two crops, followed by Selenocystine and methylselenocysteine. Through the analysis of the rhizosphere soil, it was found that Fe-Mn oxide-bound Se exceeded 80% of the total Se. Additionally, there was a significant linear correlation between the Se content in the edible parts of the two crops and the Se content in the soil. Findings suggested that wheat was more favorable than broad beans as Se supplement crops in a Se-supplied field.

Rhizosphere microbiota is important for plant growth and health. Domestication is a process to select suitable plants to satisfy the needs of humans, which may have great impacts on the interaction between the host and its rhizosphere microbiota. Rapeseed (Brassica napus) is an important oilseed crop derived from the hybridization between Brassica rapa and Brassica oleracea ~7500 years ago. However, variations in rhizosphere microbiota along with rapeseed domestication remain poorly understood. Here, we characterized the composition and structure of the rhizosphere microbiota among diverse rapeseed accessions, including ten B. napus, two B. rapa, and three B. oleracea accessions through bacterial 16S rRNA gene sequencing. B. napus exhibited a higher Shannon index and different bacterial relative abundance compared with its wild relatives in rhizosphere microbiota. Moreover, artificial synthetic B. napus lines G3D001 and No.2127 showed significantly different rhizosphere microbiota diversity and composition from other B. napus accessions and their ancestors. The core rhizosphere microbiota of B. napus and its wild relatives was also described. FAPROTAX annotation predicted that the synthetic B. napus lines had more abundant pathways related to nitrogen metabolism, and the co-occurrence network results demonstrated that Rhodoplanes acted as hub nodes to promote nitrogen metabolism in the synthetic B. napus lines. This study provides new insights into the impacts of rapeseed domestication on the diversity and community structure of rhizosphere microbiota, which may highlight the contribution of rhizosphere microbiota to plant health.

Sclerotinia sclerotiorum is an ascomycetous fungus and hosts various mycoviruses. In this study, a novel fungal alphaflexivirus with a special genomic structure, named Sclerotinia sclerotiorum alphaflexivirus 1 (SsAFV1), was cloned from a hypovirulent strain, AHS31. Strain AHS31 was also co-infected with two botourmiaviruses and two mitoviruses. The complete genome of SsAFV1 comprised 6939 bases with four open reading frames (ORFs), a conserved 5′-untranslated region (UTR), and a poly(A) tail in the 3′ terminal; the ORF1 and ORF3 encoded a replicase and a coat protein (CP), respectively, while the function of the proteins encoded by ORF2 and ORF4 was unknown. The virion of SsAFV1 was flexuous filamentous 480–510 nm in length and 9–10 nm in diameter. The results of the alignment and the phylogenetic analysis showed that SsAFV1 is related to allexivirus and botrexvirus, such as Garlic virus X of the genus Allexivirus and Botrytis virus X of the genus Botrevirus, both with 44% amino-acid (aa) identity of replicase. Thus, SsAFV1 is a novel virus and a new genus, Sclerotexvirus, is proposed to accommodate this novel alphaflexivirus.

Verticillium wilt is a kind of plant vascular disease caused by the soilborne fungus Verticillium dahliae, which severely limits cotton production. Our previous studies showed that the endophytic fungus Gibellulopsis nigrescens CEF08111 can effectively control Verticillium wilt and induce a defense response in cotton plants. However, the comprehensive molecular mechanism governing this response is not yet clear. To study the signaling mechanism induced by strain CEF08111, the transcriptome of cotton seedlings pretreated with CEF08111 was sequenced. The results revealed 249, 3559 and 33 differentially expressed genes (DEGs) at 3, 12 and 48 h post inoculation with CEF08111, respectively. At 12 h post inoculation with CEF08111, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that the DEGs were enriched mainly in the plant–pathogen interaction, mitogen-activated protein kinase (MAPK) signaling pathway-plant, and plant hormone signal transduction pathways. Gene ontology (GO) analysis revealed that these DEGs were enriched mainly in the following terms: response to external stimulus, systemic acquired resistance, kinase activity, phosphotransferase activity, xyloglucan: xyloglucosyl transferase activity, xyloglucan metabolic process, cell wall polysaccharide metabolic process and hemicellulose metabolic process. Moreover, many genes, such as calcium-dependent protein kinase (CDPK), flagellin-sensing 2 (FLS2), resistance to Pseudomonas syringae pv. maculicola 1(RPM1) and myelocytomatosis protein 2 (MYC2), that regulate crucial points in defense-related pathways were identified and may contribute to V. dahliae resistance in cotton. Seven DEGs of the pathway phenylpropanoid biosynthesis were identified by weighted gene co-expression network analysis (WGCNA), and these genes are related to lignin synthesis. The above genes were compared and analyzed, a total of 710 candidate genes that may be related to the resistance of cotton to Verticillium wilt were identified. These results provide a basis for understanding the molecular mechanism by which the biocontrol fungus CEF08111 increases the resistance of cotton to Verticillium wilt.

Sclerotinia sclerotiorum is a broad-spectrum necrotrophic phytopathogen that can infect many plant species worldwide. The application of fungicides is a common measure for controlling Sclerotinia sclerotiorum. Due to the risk of developing resistance to fungicides, it is imperative to find ways to be environmentally friendly and even effective. Using bioactive compounds in plants to reduce the amounts of fungicides has become a clean and sustainable strategy of controlling Sclerotinia sclerotiorum. Our study found that selenium in soil mediated the phenylacetic acid-related metabolic pathway in oilseed rape and reduced the incidence rate of Sclerotinia sclerotiorum. The growth-inhibition rates of Sclerotinia sclerotiorum were observed at 25.82%, 19.67%, and 52.61% for treatments of 0.8 mg·L−1 dimethachlon, 0.1 mg·mL−1 phenylacetic acid, and dimethachlon (0.8 mg·L−1) + phenylacetic acid (0.1 mg·mL−1), respectively. Phenylacetic acid reduced the application amount of dimethachlon and enhanced the inhibition effect for Sclerotinia sclerotiorum. Results also suggested that phenylacetic acid severely damaged the morphological structure, changed the electrical conductivity, and reduced the capacity of acid production and oxalic acid secretion of Sclerotinia sclerotiorum mycelium. Further studies revealed that phenylacetic acid increased the gene-expression level of Ssodc1, Ssodc2, CWDE2 and CWDE10 in mycelium while decreasing the expression level of SsGgt1, and phenylacetic acid + dimethachlon reduced the relative expression level of SsBil. These findings verified that phenylacetic acid could partially replace the amount of dimethachlon, as well as enhance the prevention of Sclerotinia sclerotiorum by dimethachlon, which provides evidence for developing an environment-friendly method for Sclerotinia sclerotiorum control.

Heredity and variation are inherent characteristics of species and are mainly reflected in the stability and variation of the genome; the former is relative, while the latter is continuous. However, whether life has both stable genomes and extremely diverse genomes at the same time is unknown. In this study, we isolated Sclerotinia sclerotiorum strains from sclerotium samples in Quincy, Washington State, USA, and found that four single-sclerotium-isolation strains (PB4, PB273, PB615, and PB623) had almost identical genomes to the reference strain 1980 isolated in the west of Nebraska 40 years ago. The genome of strain PB4 sequenced by the next-generation sequencing (NGS) and Pacific Biosciences (PacBio) sequencing carried only 135 single nucleotide polymorphisms (SNPs) and 18 structural variations (SVs) compared with the genome of strain 1980 and 48 SNPs were distributed on Contig_20. Based on data generated by NGS, three other strains, PB273, PB615, and PB623, had 256, 275, and 262 SNPs, respectively, against strain 1980, which were much less than in strain PB4 (532 SNPs) and none of them occurred on Contig_20, suggesting much closer genomes to strain 1980 than to strain PB4. All other strains from America and China are rich in SNPs with a range of 34,391–77,618 when compared with strain 1980. We also found that there were 39–79 SNPs between strain PB4 and its sexual offspring, 53.1% of which also occurred on Contig_20. Our discoveries show that there are two types of genomes in S. sclerotiorum, one is very stable and the other tends to change constantly. Investigating the mechanism of such genome stability will enhance our understanding of heredity and variation.

Tomato, which is regarded as an important worldwide crop, is susceptible to gray mold caused by Botrytis cinerea. Selenium and methyl jasmonate can act as antifungal agents against pathogenic infections. To clarify the effect of selenium and methyl jasmonate on the fungal pathogen, the spore germination and mycelial growth of B. cinerea were investigated in vitro using the growth rate method. Additionally, the electrical conductivity, soluble protein content, malondialdehyde content and oxalic acid secretion of B. cinerea mycelium were also determined to further explore the antifungal mechanism of selenium and methyl jasmonate. The results showed that selenium application significantly increased cell membrane permeability and malondialdehyde content, and methyl jasmonate treatment decreased the soluble protein content in mycelium of B. cinerea. Furthermore, supplementation of the medium with both selenium and methyl jasmonate effectively inhibited spore germination and colony growth of B. cinerea by compromising membrane integrity, and significantly reduced soluble protein content and the oxalic acid secretion of hypha. The resulting incidence of postharvest tomato gray mold with the combination of selenium and methyl jasmonate was 34.7%, which was approximately half of that of the control. To sum up, the combined use of selenium and methyl jasmonate inhibited the normal physiological activity and pathogenicity of B. cinerea, which suggests that selenium and methyl jasmonate have the potential for controlling gray mold disease caused by B. cinerea in postharvest fruits and vegetables. These findings may offer a promising and eco-friendly strategy to control gray mold disease in postharvest fruits and vegetables.

Transcriptional activity is mediated by chromatin remodeling, which in turn is affected by post-translational modifications, including histone acetylation. Histone acetyltransferases (HATs) are capable of promoting euchromatin formation and then activating gene transcription. Here, we characterize the Elp3 GNAT family HAT, which is also a subunit of Elongator complex, in the environmentally and economically important fungal insect pathogen, Beauveria bassiana. BbElp3 showed high localization levels to mitochondria, with some nuclear and cytoplasmic localization also apparent. Targeted gene knockout of BbElp3 resulted in impaired asexual development and morphogenesis, reduced tolerances to multiple stress conditions, reduced the ability of the fungus to utilize various carbon/nitrogen sources, increased susceptibility to rapamycin, and attenuated virulence in bioassays using the greater wax moth, Galleria mellonella. The ΔBbElp3 mutant also showed disrupted cell cycle, abnormal hyphal septation patterns, and enlarged autophagosomes in vegetative hyphae. Transcriptome analyses revealed differential expression of 775 genes (DEGs), including 336 downregulated and 438 upregulated genes in the ΔBbElp3 strain as compared to the wild type. Downregulated genes were mainly enriched in pathways involved in DNA processing and transcription, cell cycle control, cellular transportation, cell defense, and virulence, including hydrophobins, cellular transporters (ABC and MFS multidrug transporters), and insect cuticular degrading enzymes, while upregulated genes were mainly enriched in carbohydrate metabolism and amino acid metabolism. These data indicate pleiotropic effects of BbElp3 in impacting specific cellular processes related to asexual development, cell cycle, autophagy, and virulence.

Tomato cultivation is seriously affected by infection from Botrytis cinerea. The safe and effective control of tomato gray mold remains elusive. Plant-related microbial communities regulate not only plant metabolism but also plant immune systems. In this study, we observed that Selenium application in soil combined with foliar spraying of methyl jasmonate could reduce Botrytis cinerea infection in tomato fruits and leaves and improve tomato fruit quality. The infection rate of leaves decreased from 42.19% to 25.00%, and the vitamin C content increased by 22.14%. The bacterial community structure of the tomato was studied by using amplicon sequencing technology. The leaf bacterial alpha diversity of tomatoes treated with Se plus methyl jasmonate was significantly higher than that of the control. Then we isolated five strains antagonistic to Botrytis cinerea in vitro from tomato leaves in the treatment of Se plus methyl jasmonate. The antagonistic strains were identified as Bacillus subtilis and Bacillus velezensis. Spraying mixed antagonistic strain suspension significantly inhibited the diameter of Botrytis cinerea with an inhibition rate of 40.99%. This study revealed the key role of plant-beneficial bacteria recruited by Se combined with methyl jasmonate in improving tomato plant disease resistance. These findings may benefit our understanding of the new regulation of microorganisms on Botrytis cinerea.

Clubroot, a soil-infective disease caused by Plasmodiophora brassicae, is a serious disease affecting cruciferous plants around the world. There is no effective control measure to completely remove this pathogen from fields after infection. Here, we screened and identified two strains (Hz36, Trichoderma guizhouense; Hk37, Trichoderma koningiopsis) of Trichoderma from the gall of clubroot in rapeseed fields with biocontrol potential for clubroot. The fermentation broth of Hz36 could significantly inhibit the germination of resting spores of P. brassicae, and promote the seed germination and root growth of rapeseed. The biocontrol efficiency of Hz36 strain on clubroot for rapeseed and Arabidopsis thaliana was 44.29% and 52.18%, respectively. The qPCR results revealed that strain Hz36 treatment could significantly reduce the content of P. brassicae in root cells, and paraffin section analysis revealed that it could delay the development of P. brassicae. Strain Hk37 showed similar effects to strain Hz36, whose biocontrol efficiency of clubroot could reach 57.30% in rapeseed and 68.01% in A. thaliana. These results indicate that strains Hz36 and Hk37 have the potential for the biocontrol of clubroot.

Codon usage bias (CUB) could reflect co-evolutionary changes between viruses and hosts in contrast to plant and animal viruses, and the systematic analysis of codon usage among the mycoviruses that infect plant pathogenic fungi is limited. We performed an extensive analysis of codon usage patterns among 98 characterized RNA mycoviruses from eight phytopathogenic fungi. The GC and GC3s contents of mycoviruses have a wide variation from 29.35% to 64.62% and 24.32% to 97.13%, respectively. Mycoviral CUB is weak, and natural selection plays a major role in the formation of mycoviral codon usage pattern. In this study, we demonstrated that the codon usage of mycoviruses is similar to that of some host genes, especially those involved in RNA biosynthetic process and transcription, suggesting that CUB is a potential evolutionary mechanism that mycoviruses adapt to in their hosts.

Mycoviruses are ubiquitous in pathogenic fungi including Sclerotinia sclerotiorum. Using RNA sequencing, more mycoviruses have been identified in individual strains, which were previously reported to be infected by a single mycovirus. A hypovirulent strain of S. sclerotiorum, HC025, was previously thought to harbor a single mitovirus, Sclerotinia sclerotiorum mitovirus 1 (SsMV1), based on the analysis of the conventional dsRNA extraction method. We found HC025 to be co-infected by five mycoviruses. In addition to SsMV1, four mycoviruses were identified: Sclerotinia sclerotiorum narnavirus 4 (SsNV4), Sclerotinia sclerotiorum negative-stranded RNA virus 1 (SsNSRV1), Sclerotinia sclerotiorum ourmia-like virus 14 (SsOLV14), and SsOLV22. Three mycoviruses including SsNV4, SsNSRV1, and SsOLV14 share high replicase identities (more than 95%) with the previously reported corresponding mycoviruses, and SsOLV22 shows lower identity to the known viruses. The complete genome of SsOLV22 is 3987 nt long and contains a single ORF-encoded RdRp, which shares 24.84% identity with the RNA-dependent RNA polymerase (RdRp) of Hubei narna-like virus 10 (query coverage: 26%; e-value: 8 × 10−19). The phylogenetic tree of RdRp suggests that SsOLV22 is a new member within the family Botourmiaviridae. All of the mycoviruses except for SsNSRV1 could horizontally co-transfer from HC025 to the virulent strain Ep-1PNA367 with hypovirulent phenotypes, and converted a later strain into a hypovirulent strain. In summary, we molecularly characterized the hypovirulent strain HC025 and identified five RNA mycoviruses including a new member within Botourmiaviridae.

The plant cuticle is one of the protective layers of the external surface of plant tissues. Plants use the cuticle layer to reduce water loss and resist pathogen infection. Fungi release cell wall-degrading enzymes to destroy the epidermis of plants to achieve the purpose of infection. Sclerotinia sclerotiorum secretes a large amount of cutinase to disrupt the cuticle layer of plants during the infection process. In order to further understand the role of cutinase in the pathogenic process of S. sclerotiorum, the S. sclerotiorum cutinsae 1 (SsCut1) gene was cloned and analyzed. The protein SsCut1 contains the conserved cutinase domain and a fungal cellulose-binding domain. RT-qPCR results showed that the expression of SsCut1 was significantly upregulated during infection. Split-Marker recombination was utilized for the deletion of the SsCut1 gene, ΔSsCut1 mutants showed reduced cutinase activity and virulence, but the deletion of the SsCut1 gene had no effect on the growth rate, colony morphology, oxalic acid production, infection cushion formation and sclerotial development. Complementation with the wild-type SsCut1 allele restored the cutinase activity and virulence to the wild-type level. Interestingly, expression of SsCut1 in plants can trigger defense responses, but it also enhanced plant susceptibility to SsCut1 gene knock-out mutants. Taken together, our finding demonstrated that the SsCut1 gene promotes the virulence of S. sclerotiorum by enhancing its cutinase activity.

In pursuing higher rice production, we have often jeopardized soil at an alarming rate. It is hypothesized that intensive farming practices degrade soil health and increase the abundance of rice diseases while fallowing increases the abundance of biocontrol agents. In this study, the bacterial community was monitored in the paddy and fallow soil. Proteobacteria, Bacteroidetes, and Actinobacteria were abundant in the fallow soil, whereas Acidobacteria, Chloroflexi, and Gemmatimonadetes were more abundant in the paddy soil. Among the potential biocontrol agents, Bacillus, Thiobacillus, Rhizobium, Massilia, Rhizobacter, Streptomyces, Micromonospora, and Pseudonocardia were more abundant in the fallow soil, while Pseudomonas and Burkholderia were more abundant in the paddy soil. The possible rice pathogens, i.e., Xanthomonas and Erwinia, were more abundant in the paddy soil. The alpha diversity was higher in paddy soil than in fallow soil. Additionally, the principal coordinate analysis based on UniFrac distances revealed distinct clusters of the soils. Moreover, the functional analysis suggested that the fallow soil was abundant in genes associated with the biosynthesis of siderophores and secondary metabolites. In contrast, the paddy soil was abundant in genes related to plant-pathogen interactions. In conclusion, these results highlight the significance of fallowing to improve soil health.

Endornaviruses are capsidless linear (+) ssRNA viruses in the family Endornaviridae. In this study, Scelrotinia sclerotiorum endornavirus 11 (SsEV11), a novel endornavirus infecting hypovirulent Sclerotinia sclerotiorum strain XY79, was identified and cloned using virome sequencing analysis and rapid amplification of cDNA ends (RACE) techniques. The full-length genome of SsEV11 is 11906 nt in length with a large ORF, which encodes a large polyprotein of 3928 amino acid residues, containing a viral methyltransferase domain, a cysteine-rich region, a putative DEADc, a viral helicase domain, and an RNA-dependent RNA polymerase (RdRp) 2 domain. The 5’ and 3’ untranslated regions (UTR) are 31 nt and 90 nt, respectively. According to the BLAST result of the nucleotide sequence, SsEV11 shows the highest identity (45%) with Sclerotinia minor endornavirus 1 (SmEV1). Phylogenetic analysis based on amino acid sequence of RdRp demonstrated that SsEV11 clusters to endornavirus and has a close relationship with Betaendornavirus. Phylogenetic analysis based on the sequence of endornaviral RdRp domain indicated that there were three large clusters in the phylogenetic tree. Combining the results of alignment analysis, Cluster I at least has five subclusters including typical members of Alphaendornavirus and many unclassified endornaviruses that isolated from fungi, oomycetes, algae, and insects; Cluster II also has five subclusters including typical members of Betaendornavirus, SsEV11, and other unclassified viruses that infected fungi; Cluster III includes many endorna-like viruses that infect nematodes, mites, and insects. Viruses in Cluster I and Cluster II are close to each other and relatively distant to those in Cluster III. Our study characterized a novel betaendornavirus, SsEV11, infected fungal pathogen S. sclerotiorum, and suggested that notable phylogenetic diverse exists in endornaviruses. In addition, at least, one novel genus, Gammaendornavirus, should be established to accommodate those endorna-like viruses in Cluster III.

A novel mycovirus named Fusarium oxysporum alternavirus 1(FoAV1) was identified as infecting Fusarium oxysporum strain BH19, which was isolated from a fusarium wilt diseased stem of Lilium brownii. The genome of FoAV1 contains four double-stranded RNA (dsRNA) segments (dsRNA1, dsRNA 2, dsRNA 3 and dsRNA 4, with lengths of 3.3, 2.6, 2.3 and 1.8 kbp, respectively). Additionally, dsRNA1 encodes RNA-dependent RNA polymerase (RdRp), and dsRNA2- dsRNA3- and dsRNA4-encoded hypothetical proteins (ORF2, ORF3 and ORF4), respectively. A homology BLAST search, along with multiple alignments based on RdRp, ORF2 and ORF3 sequences, identified FoAV1 as a novel member of the proposed family “Alternaviridae”. Evolutionary relation analyses indicated that FoAV1 may be related to alternaviruses, thus dividing the family “Alternaviridae” members into four clades. In addition, we determined that dsRNA4 was dispensable for replication and may be a satellite-like RNA of FoAV1—and could perhaps play a role in the evolution of alternaviruses. Our results provided evidence for potential genera establishment within the proposed family “Alternaviridae”. Additionally, FoAV1 exhibited biological control of Fusarium wilt. Our results also laid the foundations for the further study of mycoviruses within the family “Alternaviridae”, and provide a potential agent for the biocontrol of diseases caused by F. oxysporum.

The infection by a single-stranded DNA virus, Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1 (SsHADV-1), causes hypovirulence, a reduced growth rate, and other colony morphological changes in its host Sclerotinia sclerotiorum strain DT-8. However, the mechanisms of the decline are still unclear. Using digital RNA sequencing, a transcriptome analysis was conducted to elucidate the phenotype-related genes with expression changes in response to SsHADV-1 infection. A total of 3110 S. sclerotiorum differentially expressed genes (DEGs) were detected during SsHADV-1 infection, 1741 of which were up-regulated, and 1369 were down-regulated. The identified DEGs were involved in several important pathways. DNA replication, DNA damage response, carbohydrate and lipid metabolism, ribosomal assembly, and translation were the affected categories in S. sclerotiorum upon SsHADV-1 infection. Moreover, the infection of SsHADV-1 also suppressed the expression of antiviral RNA silencing and virulence factor genes. These results provide further detailed insights into the effects of SsHADV-1 infection on the whole genome transcription in S. sclerotiorum.

Considering the huge economic loss caused by postharvest diseases, the identification and prevention of citrus postharvest diseases is vital to the citrus industry. In 2018, 16 decayed citrus fruit from four citrus varieties—Satsuma mandarin (Citrus unshiu), Ponkan (Citrus reticulata Blanco cv. Ponkan), Nanfeng mandarin (Citrus reticulata cv. nanfengmiju), and Sugar orange (Citrus reticulata Blanco)—showing soft rot and sogginess on their surfaces and covered with white mycelia were collected from storage rooms in seven provinces. The pathogens were isolated and the pathogenicity of the isolates was tested. The fungal strains were identified as Lasiodiplodia pseudotheobromae based on their morphological characteristics and phylogenetic analyses using the internal transcribed spacer regions (ITS), translation elongation factor 1-α gene (TEF), and beta-tubulin (TUB) gene sequences. The strains could infect wounded citrus fruit and cause decay within two days post inoculation, but could not infect unwounded fruit. To our knowledge, this is the first report of citrus fruit decay caused by L. pseudotheobromae in China.

Via virome sequencing, six viruses were detected from Magnaporthe oryzae strains YC81-2, including one virus in the family Tombusviridae, one virus in the family Narnaviridae and four viruses in the family Botourmiaviridae. Since the RNA-dependent RNA polymerase (RdRp) of one botourmiavirus show the highest identity (79%) with Magnaporthe oryzae ourmia-like virus 1 (MOLV1), the virus that was grouped into the genus Magoulivirus was designated as Magnaporthe oryzae botourmiavirus 2 (MOBV2). The three other novel botourmiaviruses were selected for further study. The complete nucleotide sequences of the three botourmiaviruses were determined. Sequence analysis showed that virus 1, virus 2, and virus 3 were 2598, 2385, and 2326 nts in length, respectively. The variable 3′ untranslated region (3′-UTR) and 5′-UTR of each virus could be folded into a stable stem-loop secondary structure. Each virus consisted of a unique ORF encoding a putative RdRp. The putative proteins with a conserved GDD motif of RdRp showed the highest sequence similarity to RdRps of viruses in the family Botourmiaviridae. Phylogenetic analysis demonstrated that these viruses were three distinct novel botourmiaviruses, clustered into the Botourmiaviridae family but not belonging to any known genera of this family. Thus, virus 1, virus 2, and virus 3 were designated as Magnaporthe oryzae botourmiavirus 5, 6, and 7 (MOBV5, MOBV6, and MOBV7), respectively. Our results suggest that four distinct botourmiaviruses, MOBV2, MOBV5, MOBV6, and MOBV7, co-infect a single strain of Magnaporthe oryzae, and MOBV5, MOBV6, and MOBV7 are members of three unclassified genera in the family Botourmiaviridae.

We previously identified Sclerotinia sclerotiorum negative-stranded virus 1 (SsNSRV-1), the first (−) ssRNA mycovirus, associated with hypovirulence of its fungal host Sclerotinia sclerotiorum. In this study, functional analysis of Open Reading Frame Ι (ORF Ι) of SsNSRV-1 was performed. The integration and expression of ORF Ι led to defects in hyphal tips, vegetative growth, and virulence of the mutant strains of S. sclerotiorum. Further, differentially expressed genes (DEGs) responding to the expression of ORF Ι were identified by transcriptome analysis. In all, 686 DEGs consisted of 267 up-regulated genes and 419 down-regulated genes. DEGs reprogramed by ORF Ι were relevant to secretory proteins, pathogenicity, transcription, transmembrane transport, protein biosynthesis, modification, and metabolism. Alternative splicing was also detected in all mutant strains, but not in hypovirulent strain AH98, which was co-infected by SsNSRV-1 and Sclerotinia sclerotiorum hypovirus 1 (SsHV-1). Thus, the integrity of SsNSRV-1 genome may be necessary to protect viral mRNA from splicing and inactivation by the host. Taken together, the results suggested that protein ORF Ι could regulate the transcription, translation, and modification of host genes in order to facilitate viral proliferation and reduce the virulence of the host. Therefore, ORF Ι may be a potential gene used for the prevention of S. sclerotiorum.

Botrytis cinerea, a ubiquitous necrotrophic plant-pathogenic fungus, is responsible for grey mold and rot disease in a very wide range of plant species. Subtilisin-like proteases (or subtilases) are a very diverse family of serine proteases present in many organisms and are reported to have a broad spectrum of biological functions. Here, we identified two genes encoding subtilisin-like proteases (Bcser1 and Bcser2) in the genome of B. cinerea, both of which contain an inhibitor I9 domain and a peptidase S8 domain. The expression levels of Bcser1 and Bcser2 increased during the sclerotial forming stage, as well as during a later stage of hyphal infection on Arabidopsis thaliana leaves, but the up-regulation of Bcser1 was significantly higher than that of Bcser2. Interestingly, deletion of Bcser1 had no effect on the fungal development or virulence of B. cinerea. However, deletion of Bcser2 or double deletion of Bcser1 and Bcser2 severely impaired the hyphal growth, sclerotial formation and conidiation of B. cinerea. We also found that ∆Bcser2 and ∆Bcser1/2 could not form complete infection cushions and then lost the ability to infect intact plant leaves of Arabidopsis and tomato but could infect wounded plant tissues. Taken together, our results indicate that the subtilisin-like protease Bcser2 is crucial for the sclerotial formation, conidiation, and virulence of B. cinerea.

Recent studies have demonstrated that a diverse array of mycoviruses infect the plant pathogenic fungus Sclerotinia sclerotiorum. Here, we report the molecular characterization of a newly identified mycovirus, Sclerotinia sclerotiorum fusarivirus 1 (SsFV1), which was isolated from a sclerotia-defective strain JMTJ14 of S. sclerotiorum. Excluding a poly (A) tail, the genome of SsFV1 comprises 7754 nucleotides (nts) in length with 83 and 418 nts for 5'- and 3'-untranslated regions, respectively. SsFV1 has four non-overlapping open reading frames (ORFs): ORF1 encodes a 191 kDa polyprotein (1664 amino acid residues in length) containing conserved RNA-dependent RNA polymerase (RdRp) and helicase domains; the other three ORFs encode three putative hypothetical proteins of unknown function. Phylogenetic analysis, based on RdRp and Helicase domains, indicated that SsFV1 is phylogenetically related to Rosellinia necatrix fusarivirus 1 (RnFV1), Fusarium graminearum virus-DK21 (FgV1), and Penicillium roqueforti RNA mycovirus 1 (PrRV1), a cluster of an independent group belonging to a newly proposed family Fusarividae. However, SsFV1 is markedly different from FgV1 and RnFV1 in genome organization and nucleotide sequence. SsFV1 was transmitted successfully to two vegetatively incompatible virus-free strains. SsFV1 is not responsible for the abnormal phenotype of strain JMTJ14.

Infection by diverse mycoviruses is a common phenomenon in Sclerotinia sclerotiorum. In this study, the full genome of a single-stranded RNA mycovirus, tentatively named Hubei sclerotinia RNA virus 1 (HuSRV1), was determined in the hypovirulent strain 277 of S. sclerotiorum. The HuSRV1 genome is 4492 nucleotides (nt) long and lacks a poly (A) tail at the 3ˊ- terminus. Sequence analyses showed that the HuSRV1 genome contains four putative open reading frames (ORFs). ORF1a was presumed to encode a protein with a conserved protease domain and a transmembrane domain. This protein is 27% identical to the P2a protein encoded by the subterranean clover mottle virus. ORF1b encodes a protein containing a conserved RNA-dependent RNA polymerase (RdRp) domain, which may be translated into a fusion protein by a -1 ribosome frameshift. This protein is 45.9% identical to P2b encoded by the sowbane mosaic virus. ORF2 was found to encode a putative coat protein, which shares 23% identical to the coat protein encoded by the olive mild mosaic virus. ORF3 was presumed to encode a putative protein with an unknown function. Evolutionary relation analyses indicated that HuSRV1 is related to members within Sobemovirus, but forms a unique phylogenetic branch, suggesting that HuSRV1 represents a new member within Solemoviridae. HuSRV1 virions, approximately 30 nm in diameter, were purified from strain 277. The purified virions were successfully introduced into virulent strain Ep-1PNA367, resulting in a new hypovirulent strain, which confirmed that HuSRV1 confers hypovirulence on S. sclerotiorum.

We previously determined that virions of Sclerotinia sclerotiorum hypovirulence associated DNA virus 1 (SsHADV-1) could directly infect hyphae of Sclerotinia sclerotiorum, resulting in hypovirulence of the fungal host. However, the molecular mechanisms of SsHADV-1 virions disruption of the fungal cell wall barrier and entrance into the host cell are still unclear. To investigate the early response of S. sclerotiorum to SsHADV-1 infection, S. sclerotiorum hyphae were inoculated with purified SsHADV-1 virions. The pre- and post-infection hyphae were collected at one–three hours post-inoculation for transcriptome analysis. Further, bioinformatic analysis showed that differentially expressed genes (DEGs) regulated by SsHADV-1 infection were identified in S. sclerotiorum. In total, 187 genes were differentially expressed, consisting of more up-regulated (114) than down-regulated (73) genes. The identified DEGs were involved in several important pathways. Metabolic processes, biosynthesis of antibiotics, and secondary metabolites were the most affected categories in S. sclerotiorum upon SsHADV-1 infection. Cell structure analysis suggested that 26% of the total DEGs were related to membrane tissues. Furthermore, 10 and 27 DEGs were predicted to be located in the cell membrane and mitochondria, respectively. Gene ontology enrichment analyses of the DEGs were performed, followed by functional annotation of the genes. Interestingly, one third of the annotated functional DEGs could be involved in the Ras-small G protein signal transduction pathway. These results revealed that SsHADV-1 virions may be able to bind host membrane proteins and influence signal transduction through Ras-small G protein-coupled receptors during early infection, providing new insight towards the molecular mechanisms of virions infection in S. sclerotiorum.

Various mycoviruses have been isolated from Sclerotinia sclerotiorum. Here, we identified a viral RNA sequence contig, representing a novel virus, Sclerotinia sclerotiorum deltaflexivirus 2 (SsDFV2), from an RNA_Seq database. We found that SsDFV2 was harbored in the hypovirulent strain, 228, which grew slowly on potato dextrose agar, produced a few sclerotia, and could not induce typical lesions on detached rapeseed (Brassica napus) leaves. Strain 228 was also infected by Botrytis porri RNA Virus 1 (BpRV1), a virus originally isolated from Botrytis porri. The genome of SsDFV2 comprised 6711 nucleotides, excluding the poly (A) tail, and contained a single large predicted open reading frame encoding a putative viral RNA replicase. Phylogenetic analysis demonstrated that SsDFV2 is closely related to viruses in the family Deltaflexiviridae; however, it also differs significantly from members of this family, suggesting that it may represent a new species. Further we determined that SsDFV2 could be efficiently transmitted to host vegetative incompatible individuals by dual culture. To our best knowledge, this is the first report that a (+) ssRNA mycovirus can overcome the transmission limitations of the vegetative incompatibility system, a phenomenon that may facilitate the potential use of mycoviruses for the control of crop fungal diseases.