Significantly low level of small RNA accumulation derived from an encapsidated mycovirus with dsRNA genome

Reexamining the Mycovirome of Botrytis spp

Botrytis species cause gray mold disease in more than 200 crops worldwide. To control this disease, chemical fungicides are usually applied. However, more sustainable control alternatives should be explored, such as the use of hypovirulent mycovirus-infected fungal strains. To determine the mycovirome of two Botrytis species, B. cinerea and B. prunorum, we reanalyzed RNA-Seq and small RNA-Seq data using different assembly programs and an updated viral database, aiming to identify new mycoviruses that were previously not described in the same dataset. New mycoviruses were identified, including those previously reported to infect or be associated with B. cinerea and Plasmopara viticola, such as Botrytis cinerea alpha-like virus 1 and Plasmopara viticola lesion-associated ourmia-like virus 80. Additionally, two novel narnaviruses, not previously identified infecting Botrytis species, have been characterized, tentatively named Botrytis cinerea narnavirus 1 and Botrytis narnavirus 1. The analysis of small RNAs suggested that all identified mycoviruses were targeted by the antiviral fungal mechanism, regardless of the viral genome type. In conclusion, the enlarged list of newly found viruses and the application of different bioinformatics approaches have enabled the identification of novel mycoviruses not previously described in Botrytis species, expanding the already extensive list.

Mycologists and Virologists Align: Proposing Botrytis cinerea for Global Mycovirus Studies

Mycoviruses are highly genetically diverse and can significantly change their fungal host's phenotype, yet they are generally under-described in genotypic and biological studies. We propose Botrytis cinerea as a model mycovirus system in which to develop a deeper understanding of mycovirus epidemiology including diversity, impact, and the associated cellular biology of the host and virus interaction. Over 100 mycoviruses have been described in this fungal host. B. cinerea is an ideal model fungus for mycovirology as it has highly tractable characteristics-it is easy to culture, has a worldwide distribution, infects a wide range of host plants, can be transformed and gene-edited, and has an existing depth of biological resources including annotated genomes, transcriptomes, and isolates with gene knockouts. Focusing on a model system for mycoviruses will enable the research community to address deep research questions that cannot be answered in a non-systematic manner. Since B. cinerea is a major plant pathogen, new insights may have immediate utility as well as creating new knowledge that complements and extends the knowledge of mycovirus interactions in other fungi, alone or with their respective plant hosts. In this review, we set out some of the critical steps required to develop B. cinerea as a model mycovirus system and how this may be used in the future.

Target degradation specificity of phytoplasma effector phyllogen is regulated by the recruitment of host proteasome shuttle protein

Phytoplasmas infect a wide variety of plants and can cause distinctive symptoms including the conversion of floral organs into leaf-like organs, known as phyllody. Phyllody is induced by an effector protein family called phyllogens, which interact with floral MADS-box transcription factors (MTFs) responsible for determining the identity of floral organs. The MTF/phyllogen complex then interacts with the proteasomal shuttle protein RADIATION SENSITIVE23 (RAD23), which facilitates delivery of the MTF/phyllogen complex to the host proteasome for MTF degradation. Previous studies have indicated that the MTF degradation specificity of phyllogens is determined by their ability to bind to MTFs. However, in the present study, we discovered a novel mechanism determining the degradation specificity through detailed functional analyses of a phyllogen homologue of rice yellow dwarf phytoplasma (PHYLRYD ). PHYLRYD degraded a narrower range of floral MTFs than other phyllody-inducing phyllogens, resulting in compromised phyllody phenotypes in plants. Interestingly, PHYLRYD was able to bind to some floral MTFs that PHYLRYD was unable to efficiently degrade. However, the complex of PHYLRYD and the non-degradable MTF could not interact with RAD23. These results indicate that the MTF degradation specificity of PHYLRYD is correlated with the ability to form the MTF/PHYLRYD /RAD23 ternary complex, rather than the ability to bind to MTF. This study elucidated that phyllogen target specificity is regulated by both the MTF-binding ability and RAD23 recruitment ability of the MTF/phyllogen complex.

Interaction of EXA1 and eIF4E Family Members Facilitates Potexvirus Infection in Arabidopsis thaliana

Plant viruses depend on a number of host factors for successful infection. Deficiency of critical host factors confers recessively inherited viral resistance in plants. For example, loss of Essential for poteXvirus Accumulation 1 (EXA1) in Arabidopsis thaliana confers resistance to potexviruses. However, the molecular mechanism of how EXA1 assists potexvirus infection remains largely unknown. Previous studies reported that the salicylic acid (SA) pathway is upregulated in exa1 mutants, and EXA1 modulates hypersensitive response-related cell death during EDS1-dependent effector-triggered immunity. Here, we show that exa1-mediated viral resistance is mostly independent of SA and EDS1 pathways. We demonstrate that Arabidopsis EXA1 interacts with three members of the eukaryotic translation initiation factor 4E (eIF4E) family, eIF4E1, eIFiso4E, and novel cap-binding protein (nCBP), through the eIF4E-binding motif (4EBM). Expression of EXA1 in exa1 mutants restored infection by the potexvirus Plantago asiatica mosaic virus (PlAMV), but EXA1 with mutations in 4EBM only partially restored infection. In virus inoculation experiments using Arabidopsis knockout mutants, EXA1 promoted PlAMV infection in concert with nCBP, but the functions of eIFiso4E and nCBP in promoting PlAMV infection were redundant. By contrast, the promotion of PlAMV infection by eIF4E1 was, at least partially, EXA1 independent. Taken together, our results imply that the interaction of EXA1-eIF4E family members is essential for efficient PlAMV multiplication, although specific roles of three eIF4E family members in PlAMV infection differ. IMPORTANCE The genus Potexvirus comprises a group of plant RNA viruses, including viruses that cause serious damage to agricultural crops. We previously showed that loss of Essential for poteXvirus Accumulation 1 (EXA1) in Arabidopsis thaliana confers resistance to potexviruses. EXA1 may thus play a critical role in the success of potexvirus infection; hence, elucidation of its mechanism of action is crucial for understanding the infection process of potexviruses and for effective viral control. Previous studies reported that loss of EXA1 enhances plant immune responses, but our results indicate that this is not the primary mechanism of exa1-mediated viral resistance. Here, we show that Arabidopsis EXA1 assists infection by the potexvirus Plantago asiatica mosaic virus (PlAMV) by interacting with the eukaryotic translation initiation factor 4E family. Our results imply that EXA1 contributes to PlAMV multiplication by regulating translation.

Random mutagenesis-based screening of the interface of phyllogen, a bacterial phyllody-inducing effector, for interaction with plant MADS-box proteins

To understand protein function deeply, it is important to identify how it interacts physically with its target. Phyllogen is a phyllody-inducing effector that interacts with the K domain of plant MADS-box transcription factors (MTFs), which is followed by proteasome-mediated degradation of the MTF. Although several amino acid residues of phyllogen have been identified as being responsible for the interaction, the exact interface of the interaction has not been elucidated. In this study, we comprehensively explored interface residues based on random mutagenesis using error-prone PCR. Two novel residues, at which mutations enhanced the affinity of phyllogen to MTF, were identified. These residues, and all other known interaction-involved residues, are clustered together at the surface of the protein structure of phyllogen, indicating that they constitute the interface of the interaction. Moreover, in silico structural prediction of the protein complex using ColabFold suggested that phyllogen interacts with the K domain of MTF via the putative interface. Our study facilitates an understanding of the interaction mechanisms between phyllogen and MTF.

Similar Characteristics of siRNAs of Plant Viruses Which Replicate in Plant and Fungal Hosts

Simple Summary

RNA silencing in fungi was shown to confer antiviral defense against plant viruses. In this study, using high-throughput sequencing and bioinformatic analyses, we showed that small interfering RNAs (siRNAs) of cucumber mosaic virus and tobacco mosaic virus (TMV) which replicated in phytopathogenic fungi Rhizoctonia solani and Fusarium graminearum had similarities with viral siRNAs produced in plant hosts in regard to the size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini. Additionally, our results also determined that both F. graminearum DCL1 and DCL2 were involved in the production of TMV siRNAs. Thus, the fungal RNA silencing machineries have adaptive capabilities to recognize and process the genome of invading plant viruses.

Abstract

RNA silencing is a host innate antiviral mechanism which acts via the synthesis of viral-derived small interfering RNAs (vsiRNAs). We have previously reported the infection of phytopathogenic fungi by plant viruses such as cucumber mosaic virus (CMV) and tobacco mosaic virus (TMV). Furthermore, fungal RNA silencing was shown to suppress plant virus accumulation, but the characteristics of plant vsiRNAs associated with the antiviral response in this nonconventional host remain unknown. Using high-throughput sequencing, we characterized vsiRNA profiles in two plant RNA virus–fungal host pathosystems: CMV infection in phytopathogenic fungus Rhizoctonia solani and TMV infection in phytopathogenic fungus Fusarium graminearum. The relative abundances of CMV and TMV siRNAs in the respective fungal hosts were much lower than those in the respective experimental plant hosts, Nicotiana benthamiana and Nicotiana tabacum. However, CMV and TMV siRNAs in fungi had similar characteristics to those in plants, particularly in their size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini of vsiRNAs. The abundance of TMV siRNAs largely decreased in F. graminearum mutants with a deletion in either dicer-like 1 (dcl1) or dcl2 genes which encode key proteins for the production of siRNAs and antiviral responses. However, deletion of both dcl1 and dcl2 restored TMV siRNA accumulation in F. graminearum, indicating the production of dcl-independent siRNAs with no antiviral function in the absence of the dcl1 and dcl2 genes. Our results suggest that fungal RNA silencing recognizes and processes the invading plant RNA virus genome in a similar way as in plants.

A phytoplasma effector acts as a ubiquitin-like mediator between floral MADS-box proteins and proteasome shuttle proteins

Plant pathogenic bacteria have developed effectors to manipulate host cell functions to facilitate infection. A certain number of effectors use the conserved ubiquitin-proteasome system in eukaryotic to proteolyze targets. The proteasome utilization mechanism is mainly mediated by ubiquitin interaction with target proteins destined for degradation. Phyllogens are a family of protein effectors produced by pathogenic phytoplasmas that transform flowers into leaves in diverse plants. Here, we present a noncanonical mechanism for phyllogen action that involves the proteasome and is ubiquitin-independent. Phyllogens induce proteasomal degradation of floral MADS-box transcription factors (MTFs) in the presence of RADIATION-SENSITIVE23 (RAD23) shuttle proteins, which recruit ubiquitinated proteins to the proteasome. Intracellular localization analysis revealed that phyllogen induced colocalization of MTF with RAD23. The MTF/phyllogen/RAD23 ternary protein complex was detected not only in planta but also in vitro in the absence of ubiquitin, showing that phyllogen directly mediates interaction between MTF and RAD23. A Lys-less nonubiquitinated phyllogen mutant induced degradation of MTF or a Lys-less mutant of MTF. Furthermore, the method of sequential formation of the MTF/phyllogen/RAD23 protein complex was elucidated, first by MTF/phyllogen interaction and then RAD23 recruitment. Phyllogen recognized both the evolutionarily conserved tetramerization region of MTF and the ubiquitin-associated domain of RAD23. Our findings indicate that phyllogen functionally mimics ubiquitin as a mediator between MTF and RAD23.

Residual Effects Caused by a Past Mycovirus Infection in Fusarium circinatum

Mycoviruses are known to be difficult to cure in fungi but their spontaneous loss occurs commonly. The unexpected disappearance of mycoviruses can be explained by diverse reasons, from methodological procedures to biological events such as posttranscriptional silencing machinery. The long-term effects of a virus infection on the host organism have been well studied in the case of human viruses; however, the possible residual effect on a fungus after the degradation of a mycovirus is unknown. For that, this study analyses a possible residual effect on the transcriptome of the pathogenic fungus Fusarium circinatum after the loss of the mitovirus FcMV1. The mycovirus that previously infected the fungal isolate was not recovered after a 4-year storage period. Only 14 genes were determined as differentially expressed and were related to cell cycle regulation and amino acid metabolism. The results showed a slight acceleration in the metabolism of the host that had lost the mycovirus by the upregulation of the genes involved in essential functions for fungal development. The analysis also revealed a weak expression in the annotated genes of the RNA silencing machinery. To our knowledge, this is the first time that a potential residual effect on the host transcriptome caused by the past infection of a mycovirus is reported.

The RNA Virome and Its Dynamics in an Invasive Fruit Fly, Bactrocera dorsalis, Imply Interactions Between Host and Viruses

The oriental fruit fly, Bactrocera dorsalis (Hendel), is an important invasive agricultural insect pest with a wide host range, and has spread around the world over the last century. This evolutionary trait may have arisen primarily from interactions between B. dorsalis and other invertebrates that share the same ecological niches. The invasive behavior of B. dorsalis also frequently exposes them to diverse species of viruses. Thereby, RNA viromes may be useful microbial markers to understand the ecological evolution of B. dorsalis as well as to investigate virus-host interactions. Here, we reported eight novel RNA viruses in B. dorsalis of a lab colony, including four positive-strand RNA viruses, two negative-strand RNA viruses, and two double-stranded RNA viruses using high-throughput sequencing technology. Analysis of the virus-derived small RNAs suggested that most of these viruses may be active and trigger the host antiviral RNAi responses. The viruses were also detected in various geographical populations of B. dorsalis, implying that there is a strong association between the viromes and host. In addition, these viruses infected specific fly tissues, predominately the central nervous system and gut. Furthermore, we explored the dynamics of the viruses when hosts were exposed to short- or long-term stressors, which showed that titers of some viruses were responsively altered in the stressed B. dorsalis. The discovery of these viruses may enrich our understanding of the species diversity of RNA viruses and also provide information on viruses in association with host adaptation in insects.

Analysis of an Intrinsic Mycovirus Associated With Reduced Virulence of the Human Pathogenic Fungus Aspergillus fumigatus

Aspergillus fumigatus is an airborne fungal pathogen that causes severe infections with invasive growth in immunocompromised patients. Several mycoviruses have recently been isolated from A. fumigatus strains, but there are presently no reports of mycoviral-mediated reduction or elimination of fungal pathogenicity in vertebrate models. Here, we report the biological features of a novel mycovirus, A. fumigatus chrysovirus 41362 (AfuCV41362), isolated from the hypovirulent A. fumigatus strain IFM 41362. The AfuCV41362 genome is comprised of four dsRNAs, each of which contains a single ORF (ORF1-4). ORF1 encodes a protein with sequence similarity to RNA-dependent RNA polymerases of viruses in the family Chrysoviridae, while ORF3 encodes a putative capsid protein. Viral RNAs are expressed primarily during the germination stage, and RNA-seq analysis of virus-infected A. fumigatus at the germination stage suggested that the virus suppressed expression of several pathogenicity-associated host genes, including hypoxia adaptation and nitric oxide detoxification genes. In vitro functional analysis revealed that the virus-infected strain had reduced tolerance to environmental stressors. Virus-infected A. fumigatus strain IFM 41362 had reduced virulence in vivo compared to the virus-free strain in a mouse infection model. Furthermore, introduction of the mycovirus to a natively virus-free KU A. fumigatus strain induced virus-infected phenotypes. To identify mycovirus genes responsible for the reduced virulence of A. fumigatus, each viral ORF was ectopically expressed in the virus-free KU strain. Ectopic expression of the individual ORFs only nominally reduced virulence of the host fungus in a mouse infection model. However, we found that ORF3 and ORF4 reduced tolerance to environmental stresses in in vitro analysis. Based on these results, we suggest that the AfuCV41362 mycovirus ORF3 and ORF4 reduce fungal virulence by suppressing stress tolerance together with other viral genes, rather than alone.

Functional conservation of EXA1 among diverse plant species for the infection by a family of plant viruses

Since the propagation of plant viruses depends on various host susceptibility factors, deficiency in them can prevent viral infection in cultivated and model plants. Recently, we identified the susceptibility factor Essential for poteXvirus Accumulation 1 (EXA1) in Arabidopsis thaliana, and revealed that EXA1-mediated resistance was effective against three potexviruses. Although EXA1 homolog genes are found in tomato and rice, little is known about which viruses depend on EXA1 for their infection capability and whether the function of EXA1 homologs in viral infection is conserved across multiple plant species, including crops. To address these questions, we generated knockdown mutants using virus-induced gene silencing in two Solanaceae species, Nicotiana benthamiana and tomato. In N. benthamiana, silencing of an EXA1 homolog significantly compromised the accumulation of potexviruses and a lolavirus, a close relative of potexviruses, whereas transient expression of EXA1 homologs from tomato and rice complemented viral infection. EXA1 dependency for potexviral infection was also conserved in tomato. These results indicate that EXA1 is necessary for effective accumulation of potexviruses and a lolavirus, and that the function of EXA1 in viral infection is conserved among diverse plant species.

Chrysoviruses in Magnaporthe oryzae

Magnaporthe oryzae, the fungus that causes rice blast, is the most destructive pathogen of rice worldwide. A number of M. oryzae mycoviruses have been identified. These include Magnaporthe oryzae. viruses 1, 2, and 3 (MoV1, MoV2, and MoV3) belonging to the genus, Victorivirus, in the family, Totiviridae; Magnaporthe oryzae. partitivirus 1 (MoPV1) in the family, Partitiviridae; Magnaporthe oryzae. chrysovirus 1 strains A and B (MoCV1-A and MoCV1-B) belonging to cluster II of the family, Chrysoviridae; a mycovirus related to plant viruses of the family, Tombusviridae (Magnaporthe oryzae. virus A); and a (+)ssRNA mycovirus closely related to the ourmia-like viruses (Magnaporthe oryzae. ourmia-like virus 1). Among these, MoCV1-A and MoCV1-B were the first reported mycoviruses that cause hypovirulence traits in their host fungus, such as impaired growth, altered colony morphology, and reduced pigmentation. Recently we reported that, although MoCV1-A infection generally confers hypovirulence to fungi, it is also a driving force behind the development of physiological diversity, including pathogenic races. Another example of modulated pathogenicity caused by mycovirus infection is that of Alternaria alternata chrysovirus 1 (AaCV1), which is closely related to MoCV1-A. AaCV1 exhibits two contrasting effects: Impaired growth of the host fungus while rendering the host hypervirulent to the plant, through increased production of the host-specific AK-toxin. It is inferred that these mycoviruses might be epigenetic factors that cause changes in the pathogenicity of phytopathogenic fungi.

Diverse, Novel Mycoviruses From the Virome of a Hypovirulent Sclerotium rolfsii Strain

Sclerotium rolfsii, which causes southern blight in a wide variety of crops, is a devastating plant pathogen worldwide. Mycoviruses that induce hypovirulence in phytopathogenic fungi are potential biological control resources against fungal plant diseases. However, in S. rolfsii, mycoviruses are rarely reported. In a previous study, we found a hypovirulent strain carrying a diverse pattern of dsRNAs. Here, we utilized the RNA_Seq technique to detect viral sequences. Deep sequencing, RT-PCR and Sanger sequencing validation analyses revealed that this strain harbors various new viral species that show affinity to the distinctly established and proposed families Benyviridae, Endornaviridae, Fusariviridae, Hypoviridae, and Fusagraviridae. Moreover, some viral sequences that could not be assigned to any of the existing families or orders were also identified and showed similarities to the Alphavirus, Ourmiavirus, phlegivirus-like and Curvularia thermal tolerance virus-like groups. In addition, we also conducted deep sequencing analysis of small RNAs in the virus-infecting fugal strain. The results indicated that the Dicer-mediated gene silencing mechanism was present in S. rolfsii. This is the first report of viral diversity in a single S. rolfsii fungal strain, and the results presented herein might provide insight into the taxonomy and evolution of mycoviruses and be useful for the exploration of mycoviruses as biocontrol agents.

Novel Partitivirus Enhances Virulence of and Causes Aberrant Gene Expression in Talaromyces marneffei

Talaromyces marneffei is the most important thermal dimorphic fungus causing systemic mycosis in Southeast Asia. We report the discovery of a novel partitivirus, Talaromyces marneffeipartitivirus-1 (TmPV1). TmPV1 was detected in 7 (12.7%) of 55 clinical T. marneffei isolates. Complete genome sequencing of the seven TmPV1 isolates revealed two double-stranded RNA (dsRNA) segments encoding RNA-dependent RNA polymerase (RdRp) and capsid protein, respectively. Phylogenetic analysis showed that TmPV1 occupied a distinct clade among the members of the genus Gammapartitivirus Transmission electron microscopy confirmed the presence of isometric, nonenveloped viral particles of 30 to 45 nm in diameter, compatible with partitiviruses, in TmPV1-infected T. marneffei Quantitative reverse transcription-PCR (qRT-PCR) demonstrated higher viral load of TmPV1 in the yeast phase than in the mycelial phase of T. marneffei Two virus-free isolates, PM1 and PM41, were successfully infected by purified TmPV1 using protoplast transfection. Mice challenged with TmPV1-infected T. marneffei isolates showed significantly shortened survival time (P < 0.0001) and higher fungal burden in organs than mice challenged with isogenic TmPV1-free isolates. Transcriptomic analysis showed that TmPV1 causes aberrant expression of various genes in T. marneffei, with upregulation of potential virulence factors and suppression of RNA interference (RNAi)-related genes. This is the first report of a mycovirus in a thermally dimorphic fungus. Further studies are required to ascertain the mechanism whereby TmPV1 enhances the virulence of T. marneffei in mice and the potential role of RNAi-related genes in antiviral defense in T. marneffeiIMPORTANCETalaromyces marneffei (formerly Penicillium marneffei) is the most important thermal dimorphic fungus in Southeast Asia, causing highly fatal systemic penicilliosis in HIV-infected and immunocompromised patients. We discovered a novel mycovirus, TmPV1, in seven clinical isolates of T. marneffei TmPV1 belongs to the genus Gammapartitivirus of the family Partitiviridae We showed that TmPV1 enhanced the virulence of T. marneffei in mice, with shortened survival time and higher fungal burden in the organs of mice challenged with TmPV1-infected T. marneffei isolates than in those of mice challenged with virus-free isogenic isolates. Transcriptomics analysis showed that TmPV1 altered the expression of genes involved in various cellular processes in T. marneffei, with upregulation of potential virulence factors and suppression of RNAi machinery which may be involved in antiviral defense. This is the first report of a mycovirus in a thermal dimorphic fungus. The present results offer insights into mycovirus-fungus interactions and pathogenesis of thermal dimorphic fungi.

Differential Contribution of RNA Interference Components in Response to Distinct Fusarium graminearum Virus Infections

The mechanisms of RNA interference (RNAi) as a defense response against viruses remain unclear in many plant-pathogenic fungi. In this study, we used reverse genetics and virus-derived small RNA profiling to investigate the contributions of RNAi components to the antiviral response against Fusarium graminearum viruses 1 to 3 (FgV1, -2, and -3). Real-time reverse transcription-quantitative PCR (qRT-PCR) indicated that infection of Fusarium graminearum by FgV1, -2, or -3 differentially induces the gene expression of RNAi components in F. graminearum Transcripts of the DICER-2 and AGO-1 genes of F. graminearum (FgDICER-2 and FgAGO-1) accumulated at lower levels following FgV1 infection than following FgV2 or FgV3 infection. We constructed gene disruption and overexpression mutants for each of the Argonaute and dicer genes and for two RNA-dependent RNA polymerase (RdRP) genes and generated virus-infected strains of each mutant. Interestingly, mycelial growth was significantly faster for the FgV1-infected FgAGO-1 overexpression mutant than for the FgV1-infected wild type, while neither FgV2 nor FgV3 infection altered the colony morphology of the gene deletion and overexpression mutants. FgV1 RNA accumulation was significantly decreased in the FgAGO-1 overexpression mutant. Furthermore, the levels of induction of FgAGO-1, FgDICER-2, and some of the FgRdRP genes caused by FgV2 and FgV3 infection were similar to those caused by hairpin RNA-induced gene silencing. Using small RNA sequencing analysis, we documented different patterns of virus-derived small interfering RNA (vsiRNA) production in strains infected with FgV1, -2, and -3. Our results suggest that the Argonaute protein encoded by FgAGO-1 is required for RNAi in F. graminearum, that FgAGO-1 induction differs in response to FgV1, -2, and -3, and that FgAGO-1 might contribute to the accumulation of vsiRNAs in FgV1-infected F. graminearumIMPORTANCE To increase our understanding of how RNAi components in Fusarium graminearum react to mycovirus infections, we characterized the role(s) of RNAi components involved in the antiviral defense response against Fusarium graminearum viruses (FgVs). We observed differences in the levels of induction of RNA silencing-related genes, including FgDICER-2 and FgAGO-1, in response to infection by three different FgVs. FgAGO-1 can efficiently induce a robust RNAi response against FgV1 infection, but FgDICER genes might be relatively redundant to FgAGO-1 with respect to antiviral defense. However, the contribution of this gene in the response to the other FgV infections might be small. Compared to previous studies of Cryphonectria parasitica, which showed dicer-like protein 2 and Argonaute-like protein 2 to be important in antiviral RNA silencing, our results showed that F. graminearum developed a more complex and robust RNA silencing system against mycoviruses and that FgDICER-1 and FgDICER-2 and FgAGO-1 and FgAGO-2 had redundant roles in antiviral RNA silencing.

Deep sequencing of mycovirus-derived small RNAs from Botrytis species

RNA silencing is an ancient regulatory mechanism operating in all eukaryotic cells. In fungi, it was first discovered in Neurospora crassa, although its potential as a defence mechanism against mycoviruses was first reported in Cryphonectria parasitica and, later, in several fungal species. There is little evidence of the antiviral potential of RNA silencing in the phytopathogenic species of the fungal genus Botrytis. Moreover, little is known about the RNA silencing components in these fungi, although the analysis of public genome databases identified two Dicer-like genes in B. cinerea, as in most of the ascomycetes sequenced to date. In this work, we used deep sequencing to study the virus-derived small RNA (vsiRNA) populations from different mycoviruses infecting field isolates of Botrytis spp. The mycoviruses under study belong to different genera and species, and have different types of genome [double-stranded RNA (dsRNA), (+)single-stranded RNA (ssRNA) and (-)ssRNA]. In general, vsiRNAs derived from mycoviruses are mostly of 21, 20 and 22 nucleotides in length, possess sense or antisense orientation, either in a similar ratio or with a predominance of sense polarity depending on the virus species, have predominantly U at their 5' end, and are unevenly distributed along the viral genome, showing conspicuous hotspots of vsiRNA accumulation. These characteristics reveal striking similarities with vsiRNAs produced by plant viruses, suggesting similar pathways of viral targeting in plants and fungi. We have shown that the fungal RNA silencing machinery acts against the mycoviruses used in this work in a similar manner independent of their viral or fungal origin.

Profile and functional analysis of small RNAs derived from Aspergillus fumigatus infected with double-stranded RNA mycoviruses

Background

Mycoviruses are viruses that naturally infect and replicate in fungi. Aspergillus fumigatus, an opportunistic pathogen causing fungal lung diseases in humans and animals, was recently shown to harbour several different types of mycoviruses. A well-characterised defence against virus infection is RNA silencing. The A. fumigatus genome encodes essential components of the RNA silencing machinery, including Dicer, Argonaute and RNA-dependent RNA polymerase (RdRP) homologues. Active silencing of double-stranded (ds)RNA and the generation of small RNAs (sRNAs) has been shown for several mycoviruses and it is anticipated that a similar mechanism will be activated in A. fumigatus isolates infected with mycoviruses.

Results

To investigate the existence and nature of A. fumigatus sRNAs, sRNA-seq libraries of virus-free and virus-infected isolates were created using Scriptminer adapters and compared. Three dsRNA viruses were investigated: Aspergillus fumigatus partitivirus-1 (AfuPV-1, PV), Aspergillus fumigatus chrysovirus (AfuCV, CV) and Aspergillus fumigatus tetramycovirus-1 (AfuTmV-1, NK) which were selected because they induce phenotypic changes such as coloration and sectoring. The dsRNAs of all three viruses, which included two conventionally encapsidated ones PV and CV and one unencapsidated example NK, were silenced and yielded characteristic vsiRNAs together with co-incidental silencing of host fungal genes which shared sequence homology with the viral genomes.

Conclusions

Virus-derived sRNAs were detected and characterised in the presence of virus infection. Differentially expressed A. fumigatus microRNA-like (miRNA-like) sRNAs and small interfering RNAs (siRNAs) were detected and validated. Host sRNA loci which were differentially expressed as a result of virus infection were also identified. To our knowledge, this is the first study reporting the sRNA profiles of A. fumigatus isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3773-8) contains supplementary material, which is available to authorized users.

Phytoplasma-conserved phyllogen proteins induce phyllody across the Plantae by degrading floral MADS domain proteins

ABCE-class MADS domain transcription factors (MTFs) are key regulators of floral organ development in angiosperms. Aberrant expression of these genes can result in abnormal floral traits such as phyllody. Phyllogen is a virulence factor conserved in phytoplasmas, plant pathogenic bacteria of the class Mollicutes. It triggers phyllody in Arabidopsis thaliana by inducing degradation of A- and E-class MTFs. However, it is still unknown whether phyllogen can induce phyllody in plants other than A. thaliana, although phytoplasma-associated phyllody symptoms are observed in a broad range of angiosperms. In this study, phyllogen was shown to cause phyllody phenotypes in several eudicot species belonging to three different families. Moreover, phyllogen can interact with MTFs of not only angiosperm species including eudicots and monocots but also gymnosperms and a fern, and induce their degradation. These results suggest that phyllogen induces phyllody in angiosperms and inhibits MTF function in diverse plant species.

Deficiency of the eIF4E isoform nCBP limits the cell-to-cell movement of a plant virus encoding triple-gene-block proteins in Arabidopsis thaliana

One of the important antiviral genetic strategies used in crop breeding is recessive resistance. Two eukaryotic translation initiation factor 4E family genes, eIF4E and eIFiso4E, are the most common recessive resistance genes whose absence inhibits infection by plant viruses in Potyviridae, Carmovirus, and Cucumovirus. Here, we show that another eIF4E family gene, nCBP, acts as a novel recessive resistance gene in Arabidopsis thaliana toward plant viruses in Alpha- and Betaflexiviridae. We found that infection by Plantago asiatica mosaic virus (PlAMV), a potexvirus, was delayed in ncbp mutants of A. thaliana. Virus replication efficiency did not differ between an ncbp mutant and a wild type plant in single cells, but viral cell-to-cell movement was significantly delayed in the ncbp mutant. Furthermore, the accumulation of triple-gene-block protein 2 (TGB2) and TGB3, the movement proteins of potexviruses, decreased in the ncbp mutant. Inoculation experiments with several viruses showed that the accumulation of viruses encoding TGBs in their genomes decreased in the ncbp mutant. These results indicate that nCBP is a novel member of the eIF4E family recessive resistance genes whose loss impairs viral cell-to-cell movement by inhibiting the efficient accumulation of TGB2 and TGB3.

EXA1, a GYF domain protein, is responsible for loss-of-susceptibility to plantago asiatica mosaic virus in Arabidopsis thaliana

One of the plant host resistance machineries to viruses is attributed to recessive alleles of genes encoding critical host factors for virus infection. This type of resistance, also referred to as recessive resistance, is useful for revealing plant-virus interactions and for breeding antivirus resistance in crop plants. Therefore, it is important to identify a novel host factor responsible for robust recessive resistance to plant viruses. Here, we identified a mutant from an ethylmethane sulfonate (EMS)-mutagenized Arabidopsis population which confers resistance to plantago asiatica mosaic virus (PlAMV, genus Potexvirus). Based on map-based cloning and single nucleotide polymorphism analysis, we identified a premature termination codon in a functionally unknown gene containing a GYF domain, which binds to proline-rich sequences in eukaryotes. Complementation analyses and robust resistance to PlAMV in a T-DNA mutant demonstrated that this gene, named Essential for poteXvirus Accumulation 1 (EXA1), is indispensable for PlAMV infection. EXA1 contains a GYF domain and a conserved motif for interaction with eukaryotic translation initiation factor 4E (eIF4E), and is highly conserved among monocot and dicot species. Analysis using qRT-PCR and immunoblotting revealed that EXA1 was expressed in all tissues, and was not transcriptionally responsive to PlAMV infection in Arabidopsis plants. Moreover, accumulation of PlAMV and a PlAMV-derived replicon was drastically diminished in the initially infected cells by the EXA1 deficiency. Accumulation of two other potexviruses also decreased in exa1-1 mutant plants. Our results provided a functional annotation to GYF domain-containing proteins by revealing the function of the highly conserved EXA1 gene in plant-virus interactions.

Small RNA-Based Antiviral Defense in the Phytopathogenic Fungus Colletotrichum higginsianum

Even though the fungal kingdom contains more than 3 million species, little is known about the biological roles of RNA silencing in fungi. The Colletotrichum genus comprises fungal species that are pathogenic for a wide range of crop species worldwide. To investigate the role of RNA silencing in the ascomycete fungus Colletotrichum higginsianum, knock-out mutants affecting genes for three RNA-dependent RNA polymerase (RDR), two Dicer-like (DCL), and two Argonaute (AGO) proteins were generated by targeted gene replacement. No effects were observed on vegetative growth for any mutant strain when grown on complex or minimal media. However, Δdcl1, Δdcl1Δdcl2 double mutant, and Δago1 strains showed severe defects in conidiation and conidia morphology. Total RNA transcripts and small RNA populations were analyzed in parental and mutant strains. The greatest effects on both RNA populations was observed in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains, in which a previously uncharacterized dsRNA mycovirus [termed Colletotrichum higginsianum non-segmented dsRNA virus 1 (ChNRV1)] was derepressed. Phylogenetic analyses clearly showed a close relationship between ChNRV1 and members of the segmented Partitiviridae family, despite the non-segmented nature of the genome. Immunoprecipitation of small RNAs associated with AGO1 showed abundant loading of 5’U-containing viral siRNA. C. higginsianum parental and Δdcl1 mutant strains cured of ChNRV1 revealed that the conidiation and spore morphology defects were primarily caused by ChNRV1. Based on these results, RNA silencing involving ChDCL1 and ChAGO1 in C. higginsianum is proposed to function as an antiviral mechanism.

Colletotrichum sp. comprises a diverse group of fungal pathogens that attack over 3000 plant species worldwide. Understanding the underlying mechanisms that govern fungal development and pathogenicity may enable more effective and sustainable approaches to crop disease management and control. In most organisms, RNA silencing is an important mechanism to control endogenous and exogenous RNA. RNA silencing utilizes small regulatory molecules (small RNAs) produced by proteins called Dicer (DCL), and exercise their function though effector proteins named Argonaute (AGO). Here, we investigated the role of RNA silencing machinery in the fungus Colletotrichum higginsianum, by generating deletions in genes encoding RNA silencing components. Severe defects were observed in both conidiation and conidia morphology in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains. Analysis of transcripts and small RNAs revealed an uncharacterized dsRNA virus persistently infecting C. higginsianum. The virus was shown (1) to be de-repressed in the Δdcl1, Δdcl1Δdcl2 and Δago1 strains, and (2) to cause the conidiation and spore mutant phenotypes. Our results indicate that C. higginsianum employs RNA silencing as an antiviral mechanism to suppress viruses and their debilitating effects.

Differential Inductions of RNA Silencing among Encapsidated Double-Stranded RNA Mycoviruses in the White Root Rot Fungus Rosellinia necatrix

RNA silencing acts as a defense mechanism against virus infection in a wide variety of organisms. Here, we investigated inductions of RNA silencing against encapsidated double-stranded RNA (dsRNA) fungal viruses (mycoviruses), including a partitivirus (RnPV1), a quadrivirus (RnQV1), a victorivirus (RnVV1), a mycoreovirus (RnMyRV3), and a megabirnavirus (RnMBV1) in the phytopathogenic fungus Rosellinia necatrix Expression profiling of RNA silencing-related genes revealed that a dicer-like gene, an Argonaute-like gene, and two RNA-dependent RNA polymerase genes were upregulated by RnMyRV3 or RnMBV1 infection but not by other virus infections or by constitutive expression of dsRNA in R. necatrix Massive analysis of viral small RNAs (vsRNAs) from the five mycoviruses showed that 19- to 22-nucleotide (nt) vsRNAs were predominant; however, their ability to form duplexes with 3' overhangs and the 5' nucleotide preferences of vsRNAs differed among the five mycoviruses. The abundances of 19- to 22-nt vsRNAs from RnPV1, RnQV1, RnVV1, RnMyRV3, and RnMBV1 were 6.8%, 1.2%, 0.3%, 13.0%, and 24.9%, respectively. Importantly, the vsRNA abundances and accumulation levels of viral RNA were not always correlated, and the origins of the vsRNAs were distinguishable among the five mycoviruses. These data corroborated diverse interactions between encapsidated dsRNA mycoviruses and RNA silencing. Moreover, a green fluorescent protein (GFP)-based sensor assay in R. necatrix revealed that RnMBV1 infection induced silencing of the target sensor gene (GFP gene and the partial RnMBV1 sequence), suggesting that vsRNAs from RnMBV1 activated the RNA-induced silencing complex. Overall, this study provides insights into RNA silencing against encapsidated dsRNA mycoviruses.

IMPORTANCE

Encapsidated dsRNA fungal viruses (mycoviruses) are believed to replicate inside their virions; therefore, there is a question of whether they induce RNA silencing. Here, we investigated inductions of RNA silencing against encapsidated dsRNA mycoviruses (a partitivirus, a quadrivirus, a victorivirus, a mycoreovirus, and a megabirnavirus) in Rosellinia necatrix We revealed upregulation of RNA silencing-related genes in R. necatrix infected with a mycoreovirus or a megabirnavirus but not with other viruses, which was consistent with the relatively high abundances of vsRNAs from the two mycoviruses. We also showed common and different molecular features and origins of the vsRNAs from the five mycoviruses. Furthermore, we demonstrated the activation of RNA-induced silencing complex by mycoviruses in R. necatrix Taken together, our data provide insights into an RNA silencing pathway against encapsidated dsRNA mycoviruses which is differentially induced among encapsidated dsRNA mycoviruses; that is, diverse replication strategies exist among encapsidated dsRNA mycoviruses.

Highly activated RNA silencing via strong induction of dicer by one virus can interfere with the replication of an unrelated virus

Viruses often coinfect single host organisms in nature. Depending on the combination of viruses in such coinfections, the interplay between them may be synergistic, apparently neutral with no effect on each other, or antagonistic. RNA silencing is responsible for many cases of interference or cross-protection between viruses, but such antagonistic interactions are usually restricted to closely related strains of the same viral species. In this study, we present an unprecedented example of RNA silencing-mediated one-way interference between unrelated viruses in a filamentous model fungus, Cryphonectria parasitica. The replication of Rosellinia necatrix victorivirus 1 (RnVV1; Totiviridae) was strongly impaired by coinfection with the prototypic member of the genus Mycoreovirus (MyRV1) or a mutant of the prototype hypovirus (Cryphonectria hypovirus 1, CHV1) lacking the RNA silencing suppressor (CHV1-Δp69). This interference was associated with marked transcriptional induction of key genes in antiviral RNA silencing, dicer-like 2 (dcl2) and argonaute-like 2 (agl2), following MyRV1 or CHV1-Δp69 infection. Interestingly, the inhibition of RnVV1 replication was reproduced when the levels of dcl2 and agl2 transcripts were elevated by transgenic expression of a hairpin construct of an endogenous C. parasitica gene. Disruption of dcl2 completely abolished the interference, whereas that of agl2 did not always lead to its abolishment, suggesting more crucial roles of dcl2 in antiviral defense. Taken altogether, these results demonstrated the susceptible nature of RnVV1 to the antiviral silencing in C. parasitica activated by distinct viruses or transgene-derived double-stranded RNAs and provide insight into the potential for broad-spectrum virus control mediated by RNA silencing.

Cell Death Triggered by a Putative Amphipathic Helix of Radish mosaic virus Helicase Protein Is Tightly Correlated With Host Membrane Modification

Systemic necrosis is one of the most severe symptoms caused by plant RNA viruses. Recently, systemic necrosis has been suggested to have similar features to a defense response referred to as the hypersensitive response (HR), a form of programmed cell death. In virus-infected plant cells, host intracellular membrane structures are changed dramatically for more efficient viral replication. However, little is known about whether this replication-associated membrane modification is the cause of the symptoms. In this study, we identified an amino-terminal amphipathic helix of the helicase encoded by Radish mosaic virus (RaMV) (genus Comovirus) as an elicitor of cell death in RaMV-infected plants. Cell death caused by the amphipathic helix had features similar to HR, such as SGT1-dependence. Mutational analyses and inhibitor assays using cerulenin demonstrated that the amphipathic helix-induced cell death was tightly correlated with dramatic alterations in endoplasmic reticulum (ER) membrane structures. Furthermore, the cell death-inducing activity of the amphipathic helix was conserved in Cowpea mosaic virus (genus Comovirus) and Tobacco ringspot virus (genus Nepovirus), both of which are classified in the family Secoviridae. Together, these results indicate that ER membrane modification associated with viral intracellular replication may be recognized to prime defense responses against plant viruses.

Comparative analysis among the small RNA populations of source, sink and conductive tissues in two different plant-virus pathosystems

Background

In plants, RNA silencing plays a fundamental role as defence mechanism against viruses. During last years deep-sequencing technology has allowed to analyze the sRNA profile of a large variety of virus-infected tissues. Nevertheless, the majority of these studies have been restricted to a unique tissue and no comparative analysis between phloem and source/sink tissues has been conducted. In the present work, we compared the sRNA populations of source, sink and conductive (phloem) tissues in two different plant virus pathosystems. We chose two cucurbit species infected with two viruses very different in genome organization and replication strategy; Melon necrotic spot virus (MNSV) and Prunus necrotic ringspot virus (PNRSV).

Results

Our findings showed, in both systems, an increase of the 21-nt total sRNAs together with a decrease of those with a size of 24-nt in all the infected tissues, except for the phloem where the ratio of 21/24-nt sRNA species remained constant. Comparing the vsRNAs, both PNRSV- and MNSV-infected plants share the same vsRNA size distribution in all the analyzed tissues. Similar accumulation levels of sense and antisense vsRNAs were observed in both systems except for roots that showed a prevalence of (+) vsRNAs in both pathosystems. Additionally, the presence of overrepresented discrete sites along the viral genome, hot spots, were identified and validated by stem-loop RT-PCR. Despite that in PNRSV-infected plants the presence of vsRNAs was scarce both viruses modulated the host sRNA profile.

Conclusions

We compare for the first time the sRNA profile of four different tissues, including source, sink and conductive (phloem) tissues, in two plant-virus pathosystems. Our results indicate that antiviral silencing machinery in melon and cucumber acts mainly through DCL4. Upon infection, the total sRNA pattern in phloem remains unchanged in contrast to the rest of the analyzed tissues indicating a certain tissue-tropism to this polulation. Independently of the accumulation level of the vsRNAs both viruses were able to modulate the host sRNA pattern.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1327-5) contains supplementary material, which is available to authorized users.

In Planta Recognition of a Double-Stranded RNA Synthesis Protein Complex by a Potexviral RNA Silencing Suppressor

RNA silencing plays an important antiviral role in plants and invertebrates. To counteract antiviral RNA silencing, most plant viruses have evolved viral suppressors of RNA silencing (VSRs). TRIPLE GENE BLOCK PROTEIN1 (TGBp1) of potexviruses is a well-characterized VSR, but the detailed mechanism by which it suppresses RNA silencing remains unclear. We demonstrate that transgenic expression of TGBp1 of plantago asiatica mosaic virus (PlAMV) induced developmental abnormalities in Arabidopsis thaliana similar to those observed in mutants of SUPPRESSOR OF GENE SILENCING3 (SGS3) and RNA-DEPENDENT RNA POLYMERASE6 (RDR6) required for the trans-acting small interfering RNA synthesis pathway. PlAMV-TGBp1 inhibits SGS3/RDR6-dependent double-stranded RNA synthesis in the trans-acting small interfering RNA pathway. TGBp1 interacts with SGS3 and RDR6 and coaggregates with SGS3/RDR6 bodies, which are normally dispersed in the cytoplasm. In addition, TGBp1 forms homooligomers, whose formation coincides with TGBp1 aggregation with SGS3/RDR6 bodies. These results reveal the detailed molecular function of TGBp1 as a VSR and shed new light on the SGS3/RDR6-dependent double-stranded RNA synthesis pathway as another general target of VSRs.

A mycoreovirus suppresses RNA silencing in the white root rot fungus, Rosellinia necatrix

RNA silencing is a fundamental antiviral response in eukaryotic organisms. We investigated the counterdefense strategy of a fungal virus (mycovirus) against RNA silencing in the white root rot fungus, Rosellinia necatrix. We generated an R. necatrix strain that constitutively induced RNA silencing of the exogenous green fluorescent protein (GFP) gene, and infected it with each of four unrelated mycoviruses, including a partitivirus, a mycoreovirus, a megabirnavirus, and a quadrivirus. Infection with a mycoreovirus (R. necatrix mycoreovirus 3; RnMyRV3) suppressed RNA silencing of GFP, while the other mycoviruses did not. RnMyRV3 reduced accumulation of GFP-small interfering (si) RNAs and increased accumulation of GFP-double-stranded (ds) RNA; suggesting that the virus interferes with the dicing of dsRNA. Moreover, an agroinfiltration assay in planta revealed that the S10 gene of RnMyRV3 has RNA silencing suppressor activity. These data corroborate the counterdefense strategy of RnMyRV3 against host RNA silencing.

Physiological stressors and invasive plant infections alter the small RNA transcriptome of the rice blast fungus, Magnaporthe oryzae

Background

The rice blast fungus, Magnaporthe oryzae is a destructive pathogen of rice and other related crops, causing significant yield losses worldwide. Endogenous small RNAs (sRNAs), including small interfering RNAs (siRNAs) and microRNAs (miRNAs) are critical components of gene regulation in many eukaryotic organisms. Recently several new species of sRNAs have been identified in fungi. This fact along with the availability of genome sequence makes M. oryzae a compelling target for sRNA profiling. We have examined sRNA species and their biosynthetic genes in M. oryzae, and the degree to which these elements regulate fungal stress responses. To this end, we have characterized sRNAs under different physiological stress conditions, which had not yet been examined in this fungus.

Results

The resulting libraries are composed of more than 37 million total genome matched reads mapping to intergenic regions, coding sequences, retrotransposons, inverted, tandem, and other repeated regions of the genome with more than half of the small RNAs arising from intergenic regions. The 24 nucleotide (nt) size class of sRNAs was predominant. A comparison to transcriptional data of M. oryzae undergoing the same physiological stresses indicates that sRNAs play a role in transcriptional regulation for a small subset of genes. Support for this idea comes from generation and characterization of mutants putatively involved in sRNAs biogenesis; our results indicate that the deletion of Dicer-like genes and an RNA-Dependent RNA Polymerase gene increases the transcriptional regulation of this subset of genes, including one involved in virulence.

Conclusions

Various physiological stressors and in planta conditions alter the small RNA profile of the rice blast fungus. Characterization of sRNA biosynthetic mutants helps to clarify the role of sRNAs in transcriptional control.

A novel victorivirus from a phytopathogenic fungus, Rosellinia necatrix, is infectious as particles and targeted by RNA silencing

A novel victorivirus, termed Rosellinia necatrix victorivirus 1 (RnVV1), was isolated from a plant pathogenic ascomycete, white root rot fungus Rosellinia necatrix, coinfected with a partitivirus. The virus was molecularly and biologically characterized using the natural and experimental hosts (chestnut blight fungus, Cryphonectria parasitica). RnVV1 was shown to have typical molecular victorivirus attributes, including a monopartite double-stranded RNA genome with two open reading frames (ORFs) encoding capsid protein (CP) and RNA-dependent RNA polymerase (RdRp), a UAAUG pentamer presumed to facilitate the coupled termination/reinitiation for translation of the two ORFs, a spherical particle structure ~40 nm in diameter, and moderate levels of CP and RdRp sequence identity (34 to 58%) to those of members of the genus Victorivirus within the family Totiviridae. A reproducible transfection system with purified RnVV1 virions was developed for the two distinct fungal hosts. Transfection assay with purified RnVV1 virions combined with virus elimination by hyphal tipping showed that the effects of RnVV1 on the phenotype of the natural host were negligible. Interestingly, comparison of the RNA silencing-competent (standard strain EP155) and -defective (Δdcl-2) strains of C. parasitica infected with RnVV1 showed that RNA silencing acted against the virus to repress its replication, which was restored by coinfection with hypovirus or transgenic expression of an RNA silencing suppressor, hypovirus p29. Phenotypic changes were observed in the Δdcl-2 strain but not in EP155. This is the first reported study on the host range expansion of a Totiviridae member that is targeted by RNA silencing.

Identification of three MAPKKKs forming a linear signaling pathway leading to programmed cell death in Nicotiana benthamiana

BACKGROUND

The mitogen-activated protein kinase (MAPK) cascade is an evolutionarily ancient mechanism of signal transduction found in eukaryotic cells. In plants, MAPK cascades are associated with responses to various abiotic and biotic stresses such as plant pathogens. MAPK cascades function through sequential phosphorylation: MAPK kinase kinases (MAPKKKs) phosphorylate MAPK kinases (MAPKKs), and phosphorylated MAPKKs phosphorylate MAPKs. Of these three types of kinase, the MAPKKKs exhibit the most divergence in the plant genome. Their great diversity is assumed to allow MAPKKKs to regulate many specific signaling pathways in plants despite the relatively limited number of MAPKKs and MAPKs. Although some plant MAPKKKs, including the MAPKKKα of Nicotiana benthamiana (NbMAPKKKα), are known to play crucial roles in plant defense responses, the functional relationship among MAPKKK genes is poorly understood. Here, we performed a comparative functional analysis of MAPKKKs to investigate the signaling pathway leading to the defense response.

RESULTS

We cloned three novel MAPKKK genes from N. benthamiana: NbMAPKKKβ, NbMAPKKKγ, and NbMAPKKKε2. Transient overexpression of full-length NbMAPKKKβ or NbMAPKKKγ or their kinase domains in N. benthamiana leaves induced hypersensitive response (HR)-like cell death associated with hydrogen peroxide production. This activity was dependent on the kinase activity of the overexpressed MAPKKK. In addition, virus-induced silencing of NbMAPKKKβ or NbMAPKKKγ expression significantly suppressed the induction of programmed cell death (PCD) by viral infection. Furthermore, in epistasis analysis of the functional relationships among NbMAPKKKβ, NbMAPKKKγ, and NbMAPKKKα (previously shown to be involved in plant defense responses) conducted by combining transient overexpression analysis and virus-induced gene silencing, silencing of NbMAPKKKα suppressed cell death induced by the overexpression of the NbMAPKKKβ kinase domain or of NbMAPKKKγ, but silencing of NbMAPKKKβ failed to suppress cell death induced by the overexpression of NbMAPKKKα or NbMAPKKKγ. Silencing of NbMAPKKKγ suppressed cell death induced by the NbMAPKKKβ kinase domain but not that induced by NbMAPKKKα.

CONCLUSIONS

These results demonstrate that in addition to NbMAPKKKα, NbMAPKKKβ and NbMAPKKKγ also function as positive regulators of PCD. Furthermore, these three MAPKKKs form a linear signaling pathway leading to PCD; this pathway proceeds from NbMAPKKKβ to NbMAPKKKγ to NbMAPKKKα.