Identification and analysis of in planta expressed genes of Magnaporthe oryzae

Comparative Transcriptomics of Fusarium graminearum and Magnaporthe oryzae Spore Germination Leading up To Infection

For fungal plant pathogens, the germinating spore provides the first interaction with the host. Spore germlings move across the plant surface and use diverse penetration strategies for ingress into plant surfaces. Penetration strategies include pressurized melanized appressoria, which facilitate physically punching through the plant cuticle, and nonmelanized appressoria, which penetrate with the help of enzymes or cuticular damage to breach the plant surface. Two well-studied plant pathogens, Fusarium graminearum and Magnaporthe oryzae, are typical of these two modes of penetration. We applied comparative transcriptomics to Fusarium graminearum and Magnaporthe oryzae to characterize the genetic programming of the early host-pathogen interface. Four sequential stages of development following spore localization on the plant surface, from spore swelling to appressorium formation, were sampled for each species on culture medium and on barley sheaths, and transcriptomic analyses were performed. Gene expression in the prepenetration stages in both species and under both conditions was similar. In contrast, gene expression in the final stage was strongly influenced by the environment. Appressorium formation involved the greatest number of differentially expressed genes. Laser-dissection microscopy was used to perform detailed transcriptomics of initial infection points by F. graminearum. These analyses revealed new and important aspects of early fungal ingress in this species. Expression of the trichothecene genes involved in biosynthesis of deoxynivalenol by F. graminearum implies that toxisomes are not fully functional until after penetration and indicates that deoxynivalenol is not essential for penetration under our conditions. The use of comparative gene expression of divergent fungi promises to advance highly effective targets for antifungal strategies. IMPORTANCE Fusarium graminearum and Magnaporthe oryzae are two of the most important pathogens of cereal grains worldwide. Despite years of research, strong host resistance has not been identified for F. graminearum, so other methods of control are essential. The pathogen takes advantage of multiple entry points to infect the host, including breaches in the florets due to senescence of flower parts and penetration of the weakened trichome bases to breach the epidermis. In contrast, M. oryzae directly punctures leaves that it infects, and resistant cultivars have been characterized. The threat of either pathogen causing a major disease outbreak is ever present. Comparative transcriptomics demonstrated its potential to reveal novel and effective disease prevention strategies that affect the initial stages of disease. Shedding light on the basis of this diversity of infection strategies will result in development of increasingly specific control strategies.

A genome-wide survey of interaction between rice and Magnaporthe oryzae via microarray analysis

The main aim of the work is to study the regulation of gene expression in the interaction between rice and Magnaporthe oryzae by gene chip technology. In this study, we mainly focused on changes of gene expression at 24, 48, and 72 hours post-inoculation (hpi), through which we could conduct a more comprehensive analysis of rice blast-related genes in the process of infection. The results showed that the experimental groups contained 460, 1227, and 3937 significant differentially expressed genes at 24, 48, and 72 hpi, respectively. Furthermore, 115 significantly differentially expressed genes were identified in response to rice blast infection at all three time points. By annotating these 115 genes, they were divided into three categories: metabolic pathways, proteins or enzymes, and organelle components. As expected, many of these genes were known rice blast-related genes; however, we discovered new genes with high fold changes. Most of them encoded conserved hypothetical proteins, and some were hypothetically conserved genes. Our study may contribute to finding new resistance genes and understanding the mechanism of rice blast development.

Comparative genome analyses reveal sequence features reflecting distinct modes of host-adaptation between dicot and monocot powdery mildew

BACKGROUND

Powdery mildew (PM) is one of the most important and widespread plant diseases caused by biotrophic fungi. Notably, while monocot (grass) PM fungi exhibit high-level of host-specialization, many dicot PM fungi display a broad host range. To understand such distinct modes of host-adaptation, we sequenced the genomes of four dicot PM biotypes belonging to Golovinomyces cichoracearum or Oidium neolycopersici.

RESULTS

We compared genomes of the four dicot PM together with those of Blumeria graminis f.sp. hordei (both DH14 and RACE1 isolates), B. graminis f.sp. tritici, and Erysiphe necator infectious on barley, wheat and grapevine, respectively. We found that despite having a similar gene number (6620-6961), the PM genomes vary from 120 to 222 Mb in size. This high-level of genome size variation is indicative of highly differential transposon activities in the PM genomes. While the total number of genes in any given PM genome is only about half of that in the genomes of closely related ascomycete fungi, most (~ 93%) of the ascomycete core genes (ACGs) can be found in the PM genomes. Yet, 186 ACGs were found absent in at least two of the eight PM genomes, of which 35 are missing in some dicot PM biotypes, but present in the three monocot PM genomes, indicating remarkable, independent and perhaps ongoing gene loss in different PM lineages. Consistent with this, we found that only 4192 (3819 singleton) genes are shared by all the eight PM genomes, the remaining genes are lineage- or biotype-specific. Strikingly, whereas the three monocot PM genomes possess up to 661 genes encoding candidate secreted effector proteins (CSEPs) with families containing up to 38 members, all the five dicot PM fungi have only 116-175 genes encoding CSEPs with limited gene amplification.

CONCLUSIONS

Compared to monocot (grass) PM fungi, dicot PM fungi have a much smaller effectorome. This is consistent with their contrasting modes of host-adaption: while the monocot PM fungi show a high-level of host specialization, which may reflect an advanced host-pathogen arms race, the dicot PM fungi tend to practice polyphagy, which might have lessened selective pressure for escalating an with a particular host.

Genome analysis of the ubiquitous boxwood pathogen Pseudonectria foliicola

Boxwood (Buxus spp.) are broad-leaved, evergreen landscape plants valued for their longevity and ornamental qualities. Volutella leaf and stem blight, caused by the ascomycete fungi Pseudonectria foliicola and P. buxi, is one of the major diseases affecting the health and ornamental qualities of boxwood. Although this disease is less severe than boxwood blight caused by Calonectria pseudonaviculata and C. henricotiae, its widespread occurrence and disfiguring symptoms have caused substantial economic losses to the ornamental industry. In this study, we sequenced the genome of P. foliicola isolate ATCC13545 using Illumina technology and compared it to other publicly available fungal pathogen genomes to better understand the biology of this organism. A de novo assembly estimated the genome size of P. foliicola at 28.7 Mb (425 contigs; N50 = 184,987 bp; avg. coverage 188×), with just 9,272 protein-coding genes. To our knowledge, P. foliicola has the smallest known genome within the Nectriaceae. Consistent with the small size of the genome, the secretome, CAzyme and secondary metabolite profiles of this fungus are reduced relative to two other surveyed Nectriaceae fungal genomes: Dactylonectria macrodidyma JAC15-245 and Fusarium graminearum Ph-1. Interestingly, a large cohort of genes associated with reduced virulence and loss of pathogenicity was identified from the P. foliicola dataset. These data are consistent with the latest observations by plant pathologists that P. buxi and most likely P. foliicola, are opportunistic, latent pathogens that prey upon weak and stressed boxwood plants.

Increase of Fungal Pathogenicity and Role of Plant Glutamine in Nitrogen-Induced Susceptibility (NIS) To Rice Blast

Highlight  Modifications in glutamine synthetase OsGS1-2 expression and fungal pathogenicity underlie nitrogen-induced susceptibility to rice blast. Understanding why nitrogen fertilization increase the impact of many plant diseases is of major importance. The interaction between Magnaporthe oryzae and rice was used as a model for analyzing the molecular mechanisms underlying Nitrogen-Induced Susceptibility (NIS). We show that our experimental system in which nitrogen supply strongly affects rice blast susceptibility only slightly affects plant growth. In order to get insights into the mechanisms of NIS, we conducted a dual RNA-seq experiment on rice infected tissues under two nitrogen fertilization regimes. On the one hand, we show that enhanced susceptibility was visible despite an over-induction of defense gene expression by infection under high nitrogen regime. On the other hand, the fungus expressed to high levels effectors and pathogenicity-related genes in plants under high nitrogen regime. We propose that in plants supplied with elevated nitrogen fertilization, the observed enhanced induction of plant defense is over-passed by an increase in the expression of the fungal pathogenicity program, thus leading to enhanced susceptibility. Moreover, some rice genes implicated in nitrogen recycling were highly induced during NIS. We further demonstrate that the OsGS1-2 glutamine synthetase gene enhances plant resistance to M. oryzae and abolishes NIS and pinpoint glutamine as a potential key nutrient during NIS.

Quantitative real-time PCR normalization for gene expression studies in the plant pathogenic fungi Lasiodiplodia theobromae

Lasiodiplodia theobromae is a highly virulent plant pathogen. It has been suggested that heat stress increases its virulence. The aim of this work was to evaluate, compare, and recommend normalization strategies for gene expression analysis of the fungus growing with grapevine wood under heat stress. Using RT-qPCR-derived data, reference gene stability was evaluated through geNorm, NormFinder and Bestkeeper applications. Based on the geometric mean using the ranking position obtained for each independent analysis, genes were ranked from least to most stable as follows: glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin (ACT), β-tubulin (TUB) and elongation factor-1α (EF1α). Using RNAseq-derived data based on the calculated tagwise dispersion these genes were ordered by increasing stability as follows: GAPDH, ACT, TUB, and EF1α. The correlation between RNAseq and RTqPCR results was used as criteria to identify the best RT-qPCR normalization approach. The gene TUB is recommended as the best option for normalization among the commonly used reference genes, but alternative fungal reference genes are also suggested.

Large-scale rewiring of innate immunity circuitry and microRNA regulation during initial rice blast infection

Rice blast is a recurrent fungal disease, and resistance to fungal infection is a complex trait. Therefore, a comprehensive examination of rice transcriptome and its variation during fungal infection is necessary to understand the complex gene regulatory networks. In this study, adopting Next-Generation Sequencing we profiled the transcriptomes and microRNAomes of rice varieties, one susceptible and the other resistant to M. oryzae, at multiple time points during the fungal infection. Our results revealed a substantial variation in the plant transcriptome and microRNAome as well as change to rice innate immunity during fungal infection. A number of putative R gene candidates were identified from a perturbed rice transcriptome analysis. The expression of genes and non-coding RNA molecules changed in both fungal resistant and susceptible plants during M. oryzae invasion discovered distinct pathways triggered in the susceptible and resistant plants. In addition, a number of fungus genes in the susceptible and resistant plants were constantly expressed at different time points, suggesting that they were likely to be the potential AVR genes. Our results revealed large-scale rewiring of innate immunity circuitry and microRNA regulation during initial rice blast infection, which would help to develop more robust blast-resistant rice plants.

Kingdom-Wide Analysis of Fungal Small Secreted Proteins (SSPs) Reveals their Potential Role in Host Association

Fungal secretome consists of various functional groups of proteins, many of which participate in nutrient acquisition, self-protection, or manipulation of the environment and neighboring organisms. The least characterized component of the secretome is small secreted proteins (SSPs). Some SSPs have been reported to function as effectors, but most remain to be characterized. The composition of major secretome components, such as carbohydrate-active enzymes, proteases, lipases, and oxidoreductases, appear to reflect the lifestyle and ecological niche of individual species. We hypothesize that many SSPs participate in manipulating plants as effectors. Obligate biotrophs likely encode more and diverse effector-like SSPs to suppress host defense compared to necrotrophs, which generally use cell wall degrading enzymes and phytotoxins to kill hosts. Because different secretome prediction workflows have been used in different studies, available secretome data are difficult to integrate for comprehensive comparative studies to test this hypothesis. In this study, SSPs encoded by 136 fungal species were identified from data archived in Fungal Secretome Database (FSD) via a refined secretome workflow. Subsequently, compositions of SSPs and other secretome components were compared in light of taxa and lifestyles. Those species that are intimately associated with host cells, such as biotrophs and symbionts, usually have higher proportion of species-specific SSPs (SSSPs) than hemibiotrophs and necrotrophs, but the latter groups displayed higher proportions of secreted enzymes. Results from our study established a foundation for functional studies on SSPs and will also help understand genomic changes potentially underpinning different fungal lifestyles.

Protein kinase C is essential for viability of the rice blast fungus Magnaporthe oryzae

Protein kinase C constitutes a family of serine–threonine kinases found in all eukaryotes and implicated in a wide range of cellular functions, including regulation of cell growth, cellular differentiation and immunity. Here, we present three independent lines of evidence which indicate that protein kinase C is essential for viability of Magnaporthe oryzae. First, all attempts to generate a target deletion of PKC1, the single copy protein kinase C‐encoding gene, proved unsuccessful. Secondly, conditional gene silencing of PKC1 by RNA interference led to severely reduced growth of the fungus, which was reversed by targeted deletion of the Dicer2‐encoding gene, MDL2. Finally, selective kinase inhibition of protein kinase C by targeted allelic replacement with an analogue‐sensitive PKC1^AS allele led to specific loss of fungal viability in the presence of the PP1 inhibitor. Global transcriptional profiling following selective PKC inhibition identified significant changes in gene expression associated with cell wall re‐modelling, autophagy, signal transduction and secondary metabolism. When considered together, these results suggest protein kinase C is essential for growth and development of M. oryzae with extensive downstream targets in addition to the cell integrity pathway. Targeting protein kinase C signalling may therefore prove an effective means of controlling rice blast disease.

Genome plasticity in filamentous plant pathogens contributes to the emergence of novel effectors and their cellular processes in the host

Plant diseases cause extensive yield loss of crops worldwide, and secretory 'warfare' occurs between plants and pathogenic organisms all the time. Filamentous plant pathogens have evolved the ability to manipulate host processes and facilitate colonization through secreting effectors inside plant cells. The stresses from hosts and environment can drive the genome dynamics of plant pathogens. Remarkable advances in plant pathology have been made owing to these adaptable genome regions of several lineages of filamentous phytopathogens. Characterization new effectors and interaction analyses between pathogens and plants have provided molecular insights into the plant pathways perturbed during the infection process. In this mini-review, we highlight promising approaches of identifying novel effectors based on the genome plasticity. We also discuss the interaction mechanisms between plants and their filamentous pathogens and outline the possibilities of effector gene expression under epigenetic control that will be future directions for research.

Global genome and transcriptome analyses of Magnaporthe oryzae epidemic isolate 98-06 uncover novel effectors and pathogenicity-related genes, revealing gene gain and lose dynamics in genome evolution

Genome dynamics of pathogenic organisms are driven by pathogen and host co-evolution, in which pathogen genomes are shaped to overcome stresses imposed by hosts with various genetic backgrounds through generation of a variety of isolates. This same principle applies to the rice blast pathogen Magnaporthe oryzae and the rice host; however, genetic variations among different isolates of M. oryzae remain largely unknown, particularly at genome and transcriptome levels. Here, we applied genomic and transcriptomic analytical tools to investigate M. oryzae isolate 98-06 that is the most aggressive in infection of susceptible rice cultivars. A unique 1.4 Mb of genomic sequences was found in isolate 98-06 in comparison to reference strain 70-15. Genome-wide expression profiling revealed the presence of two critical expression patterns of M. oryzae based on 64 known pathogenicity-related (PaR) genes. In addition, 134 candidate effectors with various segregation patterns were identified. Five tested proteins could suppress BAX-mediated programmed cell death in Nicotiana benthamiana leaves. Characterization of isolate-specific effector candidates Iug6 and Iug9 and PaR candidate Iug18 revealed that they have a role in fungal propagation and pathogenicity. Moreover, Iug6 and Iug9 are located exclusively in the biotrophic interfacial complex (BIC) and their overexpression leads to suppression of defense-related gene expression in rice, suggesting that they might participate in biotrophy by inhibiting the SA and ET pathways within the host. Thus, our studies identify novel effector and PaR proteins involved in pathogenicity of the highly aggressive M. oryzae field isolate 98-06, and reveal molecular and genomic dynamics in the evolution of M. oryzae and rice host interactions.

Genetic variations in pathogens, such as the causal agent of rice blast Magnaporthe oryzae, often lead to circumvention of disease-resistance cultivars. Previous genome-wide analyses of model organisms suggest that pathogen effectors are also rapidly evolving, especially in regions with high genome plasticity. However, genetic variations among different isolates remain largely unknown in M. oryzae, particularly at the genome and transcriptome levels. In this study, we provided a systematic genomic and interaction transcriptome profile for a dominant rice blast field isolate, resulting in identification of 134 candidate effectors. Two effectors, Iug6 and Iug9, and one pathogenicity-related (PaR) gene product, Iug18, were subjected to functional characterization. We found that Iug6 and Iug9 are located in the biotrophic interfacial complex (BIC) and their overexpression leads to suppression of defense-related gene expression in rice, while Iug18 appears to be a novel PaR protein. Our studies support the hypothesis that isolate-unique genes may serve as a source of genetic variability in the M. oryzae population encountering different environments. Our studies also facilitate further understanding of effectors and genomic variations in pathogenicity of M. oryzae.

Systematic characterization of the bZIP transcription factor gene family in the rice blast fungus, Magnaporthe oryzae

Regulatory roles of the basic leucine zipper (bZIP) transcription factors (TFs) in fungi have been identified in diverse cellular processes such as development, nutrient utilization and various stress responses. In this study, the 22 Magnaporthe oryzae genes encoding bZIP TFs were systematically characterized. Phylogenetic analysis of fungal bZIP TFs revealed that seven MobZIPs are Magnaporthe-specific, while others belongs to 15 clades of orthologous Ascomycota genes. Expression patterns of MobZIPs under various conditions showed that they are highly stress responsive. We generated deletion mutants for 13 MobZIPs: nine with orthologues in other fungal species and four Magnaporthe-specific ones. Seven of them exhibited defects in mycelial growth, development and/or pathogenicity. Consistent with the conserved functions of the orthologues, MobZIP22 and MobZIP13 played a role in sulfur metabolism and iron homeostasis respectively. Along with MobZIP22 and MobZIP13, one Magnaporthe-specific gene, MobZIP11 is essential for pathogenicity in a reactive oxygen species-dependent manner. Taken together, our results will contribute to understanding the regulatory mechanisms of the bZIP TF gene family in fungal development, adaptation to environmental stresses and pathogenicity in the rice blast fungus.

Analysis of in planta Expressed Orphan Genes in the Rice Blast Fungus Magnaporthe oryzae

Genomes contain a large number of unique genes which have not been found in other species. Although the origin of such "orphan" genes remains unclear, they are thought to be involved in species-specific adaptive processes. Here, we analyzed seven orphan genes (MoSPC1 to MoSPC7) prioritized based on in planta expressed sequence tag data in the rice blast fungus, Magnaporthe oryzae. Expression analysis using qRT-PCR confirmed the expression of four genes (MoSPC1, MoSPC2, MoSPC3 and MoSPC7) during plant infection. However, individual deletion mutants of these four genes did not differ from the wild-type strain for all phenotypes examined, including pathogenicity. The length, GC contents, codon adaptation index and expression during mycelial growth of the four genes suggest that these genes formed during the evolutionary history of M. oryzae. Synteny analyses using closely related fungal species corroborated the notion that these genes evolved de novo in the M. oryzae genome. In this report, we discuss our inability to detect phenotypic changes in the four deletion mutants. Based on these results, the four orphan genes may be products of de novo gene birth processes, and their adaptive potential is in the course of being tested for retention or extinction through natural selection.

Regulation of cellular diacylglycerol through lipid phosphate phosphatases is required for pathogenesis of the rice blast fungus, Magnaporthe oryzae

Considering implication of diacylglycerol in both metabolism and signaling pathways, maintaining proper levels of diacylglycerol (DAG) is critical to cellular homeostasis and development. Except the PIP2-PLC mediated pathway, metabolic pathways leading to generation of DAG converge on dephosphorylation of phosphatidic acid catalyzed by lipid phosphate phosphatases. Here we report the role of such enzymes in a model plant pathogenic fungus, Magnaporthe oryzae. We identified five genes encoding putative lipid phosphate phosphatases (MoLPP1 to MoLPP5). Targeted disruption of four genes (except MoLPP4) showed that MoLPP3 and MoLPP5 are required for normal progression of infection-specific development and proliferation within host plants, whereas MoLPP1 and MoLPP2 are indispensable for fungal pathogenicity. Reintroduction of MoLPP3 and MoLPP5 into individual deletion mutants restored all the defects. Furthermore, exogenous addition of saturated DAG not only restored defect in appressorium formation but also complemented reduced virulence in both mutants. Taken together, our data indicate differential roles of lipid phosphate phosphatase genes and requirement of proper regulation of cellular DAGs for fungal development and pathogenesis.

Roles of Forkhead-box Transcription Factors in Controlling Development, Pathogenicity, and Stress Response in Magnaporthe oryzae

Although multiple transcription factors (TFs) have been characterized via mutagenesis to understand their roles in controlling pathogenicity and infection-related development in Magnaporthe oryzae, the causal agent of rice blast, if and how forkhead-box (FOX) TFs contribute to these processes remain to be characterized. Four putative FOX TF genes were identified in the genome of M. oryzae, and phylogenetic analysis suggested that two of them (MoFKH1 and MoHCM1) correspond to Ascomycota-specific members of the FOX TF family while the others (MoFOX1 and MoFOX2) are Pezizomycotina-specific members. Deletion of MoFKH1 (ΔMofkh1) resulted in reduced mycelial growth and conidial germination, abnormal septation and stress response, and reduced virulence. Similarly, ΔMohcm1 exhibited reduced mycelial growth and conidial germination. Conidia of ΔMofkh1 and ΔMohcm1 were more sensitive to one or both of the cell cycle inhibitors hydroxyurea and benomyl, suggesting their role in cell cycle control. On the other hand, loss of MoFOX1 (ΔMofox1) did not show any noticeable changes in development, pathogenicity, and stress response. Deletion of MoFOX2 was not successful even after repeated attempts. Taken together, these results suggested that MoFKH1 and Mo-HCM1 are important in fungal development and that MoFKH1 is further implicated in pathogenicity and stress response in M. oryzae.

Evolutionary analysis of hydrophobin gene family in two wood-degrading basidiomycetes, Phlebia brevispora and Heterobasidion annosum s.l

Background

Hydrophobins are small secreted cysteine-rich proteins that play diverse roles during different phases of fungal life cycle. In basidiomycetes, hydrophobin-encoding genes often form large multigene families with up to 40 members. The evolutionary forces driving hydrophobin gene expansion and diversification in basidiomycetes are poorly understood. The functional roles of individual genes within such gene families also remain unclear. The relationship between the hydrophobin gene number, the genome size and the lifestyle of respective fungal species has not yet been thoroughly investigated. Here, we present results of our survey of hydrophobin gene families in two species of wood-degrading basidiomycetes, Phlebia brevispora and Heterobasidion annosum s.l. We have also investigated the regulatory pattern of hydrophobin-encoding genes from H. annosum s.s. during saprotrophic growth on pine wood as well as on culture filtrate from Phlebiopsis gigantea using micro-arrays. These data are supplemented by results of the protein structure modeling for a representative set of hydrophobins.

Results

We have identified hydrophobin genes from the genomes of two wood-degrading species of basidiomycetes, Heterobasidion irregulare, representing one of the microspecies within the aggregate H. annosum s.l., and Phlebia brevispora. Although a high number of hydrophobin-encoding genes were observed in H. irregulare (16 copies), a remarkable expansion of these genes was recorded in P. brevispora (26 copies). A significant expansion of hydrophobin-encoding genes in other analyzed basidiomycetes was also documented (1–40 copies), whereas contraction through gene loss was observed among the analyzed ascomycetes (1–11 copies). Our phylogenetic analysis confirmed the important role of gene duplication events in the evolution of hydrophobins in basidiomycetes. Increased number of hydrophobin-encoding genes appears to have been linked to the species’ ecological strategy, with the non-pathogenic fungi having increased numbers of hydrophobins compared with their pathogenic counterparts. However, there was no significant relationship between the number of hydrophobin-encoding genes and genome size. Furthermore, our results revealed significant differences in the expression levels of the 16 H. annosum s.s. hydrophobin-encoding genes which suggest possible differences in their regulatory patterns.

Conclusions

A considerable expansion of the hydrophobin-encoding genes in basidiomycetes has been observed. The distribution and number of hydrophobin-encoding genes in the analyzed species may be connected to their ecological preferences. Results of our analysis also have shown that H. annosum s.l. hydrophobin-encoding genes may be under positive selection. Our gene expression analysis revealed differential expression of H. annosum s.s. hydrophobin genes under different growth conditions, indicating their possible functional diversification.

Investigating the Ustilago maydis/Zea mays pathosystem: transcriptional responses and novel functional aspects of a fungal calcineurin regulatory B subunit

The sustainable control of basidiomycete biotrophic plant pathogenesis requires an understanding of host responses to infection, as well as the identification and functional analysis of fungal genes involved in disease development. The creation and analysis of a suppressive subtractive hybridization (SSH) cDNA library from Ustilago maydis-infected Zea mays seedlings enabled the identification of fungal and plant genes expressed during disease development, and uncovered new insights into the interactions of this model system. Candidate U. maydis pathogenesis genes were identified by using the current SSH cDNA library analysis, and by knowledge generated from previous cDNA microarray and comparative genomic analyses. These identifications were supported by the independent determination of transcript level changes in different cell-types and during pathogenic development. The basidiomycete specific um01632, the highly in planta expressed um03046 (zig1), and the calcineurin regulatory B subunit (um10226, cnb1), were chosen for deletion experiments. um01632 and zig1 mutants showed no difference in morphology and did not have a statistically significant impact on pathogenesis. cnb1 mutants had a distinct cell division phenotype and reduced virulence in seedling assays. Infections with reciprocal wild-type×Δcnb1 haploid strain crosses revealed that the wild-type allele was unable to fully compensate for the lack of a second cnb1 allele. This haploinsufficiency was undetected in other fungal cnb1 mutational analyses. The reported data improves U. maydis genome annotation and expands on the current understanding of pathogenesis genes in this model basidiomycete.

Comparative genomics and evolutionary analysis of hydrophobins from three species of wood-degrading fungi

Hydrophobins are small, secreted proteins playing important roles at different stages of fungal life cycles. Their characteristic feature is the presence of eight highly conserved cysteine residues. Here we present an inventory and evolutionary analysis of hydrophobin genes from three wood-degrading basidiomycetes, Phlebia brevispora, Ganoderma sp. and Bjerkandera adusta. The genomes of the three analyzed species are characterized by the presence of high copy numbers of hydrophobin genes. Results of the phylogenetic analysis of the identified proteins revealed that many of them share a high degree of sequence similarity and probably originated from a series of duplication events. The presence of several clusters of adjacent copies of the hydrophobin gene in a particular location in the genome further supports the interpretation that gene duplication has played a role in the evolution of hydrophobins in the analyzed species.

Temporal analysis of the magnaporthe oryzae proteome during conidial germination and cyclic AMP (cAMP)-mediated appressorium formation

Rice blast disease caused by Magnaporthe oryzae is one of the most serious threats to global rice production. During the earliest stages of rice infection, M. oryzae conidia germinate on the leaf surface and form a specialized infection structure termed the appressorium. The development of the appressorium represents the first critical stage of infectious development. A total of 3200 unique proteins were identified by nanoLC-MS/MS in a temporal study of conidial germination and cAMP-induced appressorium formation in M. oryzae. Using spectral counting based label free quantification, observed changes in relative protein abundance during the developmental process revealed changes in the cell wall biosynthetic machinery, transport functions, and production of extracellular proteins in developing appressoria. One hundred and sixty-six up-regulated and 208 down-regulated proteins were identified in response to cAMP treatment. Proteomic analysis of a cAMP-dependent protein kinase A mutant that is compromised in the ability to form appressoria identified proteins whose developmental regulation is dependent on cAMP signaling. Selected reaction monitoring was used for absolute quantification of four regulated proteins to validate the global proteomics data and confirmed the germination or appressorium specific regulation of these proteins. Finally, a comparison of the proteome and transcriptome was performed and revealed little correlation between transcript and protein regulation. A subset of regulated proteins were identified whose transcripts show similar regulation patterns and include many of the most strongly regulated proteins indicating a central role in appressorium formation. A temporal quantitative RT-PCR analysis confirmed a strong correlation between transcript and protein abundance for some but not all genes. Collectively, the data presented here provide the first comprehensive view of the M. oryzae proteome during early infection-related development and highlight biological processes important for pathogenicity.

Global expression profiling of transcription factor genes provides new insights into pathogenicity and stress responses in the rice blast fungus

Because most efforts to understand the molecular mechanisms underpinning fungal pathogenicity have focused on studying the function and role of individual genes, relatively little is known about how transcriptional machineries globally regulate and coordinate the expression of a large group of genes involved in pathogenesis. Using quantitative real-time PCR, we analyzed the expression patterns of 206 transcription factor (TF) genes in the rice blast fungus Magnaporthe oryzae under 32 conditions, including multiple infection-related developmental stages and various abiotic stresses. The resulting data, which are publicly available via an online platform, provided new insights into how these TFs are regulated and potentially work together to control cellular responses to a diverse array of stimuli. High degrees of differential TF expression were observed under the conditions tested. More than 50% of the 206 TF genes were up-regulated during conidiation and/or in conidia. Mutations in ten conidiation-specific TF genes caused defects in conidiation. Expression patterns in planta were similar to those under oxidative stress conditions. Mutants of in planta inducible genes not only exhibited sensitive to oxidative stress but also failed to infect rice. These experimental validations clearly demonstrated the value of TF expression patterns in predicting the function of individual TF genes. The regulatory network of TF genes revealed by this study provides a solid foundation for elucidating how M. oryzae regulates its pathogenesis, development, and stress responses.

Identification and characterization of in planta-expressed secreted effector proteins from Magnaporthe oryzae that induce cell death in rice

Interactions between rice and Magnaporthe oryzae involve the recognition of cellular components and the exchange of complex molecular signals from both partners. How these interactions occur in rice cells is still elusive. We employed robust-long serial analysis of gene expression, massively parallel signature sequencing, and sequencing by synthesis to examine transcriptome profiles of infected rice leaves. A total of 6,413 in planta-expressed fungal genes, including 851 genes encoding predicted effector proteins, were identified. We used a protoplast transient expression system to assess 42 of the predicted effector proteins for the ability to induce plant cell death. Ectopic expression assays identified five novel effectors that induced host cell death only when they contained the signal peptide for secretion to the extracellular space. Four of them induced cell death in Nicotiana benthamiana. Although the five effectors are highly diverse in their sequences, the physiological basis of cell death induced by each was similar. This study demonstrates that our integrative genomic approach is effective for the identification of in planta-expressed cell death-inducing effectors from M. oryzae that may play an important role facilitating colonization and fungal growth during infection.

Simultaneous RNA-seq analysis of a mixed transcriptome of rice and blast fungus interaction

A filamentous fungus, Magnaporthe oryzae, is a causal agent of rice blast disease, which is one of the most serious diseases affecting cultivated rice, Oryza sativa. However, the molecular mechanisms underlying both rice defense and fungal attack are not yet fully understood. Extensive past studies have characterized many infection-responsive genes in the pathogen and host plant, separately. To understand the plant-pathogen interaction comprehensively, it is valuable to monitor the gene expression profiles of both interacting organisms simultaneously in the same infected plant tissue. Although the host-pathogen interaction during the initial infection stage is important for the establishment of infection, the detection of fungal gene expression in infected leaves at the stage has been difficult because very few numbers of fungal cells are present. Using the emerging RNA-Seq technique, which has a wide dynamic range for expression analyses, we analyzed the mixed transcriptome of rice and blast fungus in infected leaves at 24 hours post-inoculation, which is the point when the primary infection hyphae penetrate leaf epidermal cells. We demonstrated that our method detected the gene expression of both the host plant and pathogen simultaneously in the same infected leaf blades in natural infection conditions without any artificial treatments. The upregulation of 240 fungal transcripts encoding putative secreted proteins was observed, suggesting that these candidates of fungal effector genes may play important roles in initial infection processes. The upregulation of transcripts encoding glycosyl hydrolases, cutinases and LysM domain-containing proteins were observed in the blast fungus, whereas pathogenesis-related and phytoalexin biosynthetic genes were upregulated in rice. Furthermore, more drastic changes in expression were observed in the incompatible interactions compared with the compatible ones in both rice and blast fungus at this stage. Our mixed transcriptome analysis is useful for the simultaneous elucidation of the tactics of host plant defense and pathogen attack.

Comparative analysis of zinc finger proteins involved in plant disease resistance

A meta-analysis was performed to understand the role of zinc finger domains in proteins of resistance (R) genes cloned from different crops. We analyzed protein sequences of seventy R genes of various crops in which twenty six proteins were found to have zinc finger domains along with nucleotide binding sites - leucine rice repeats (NBS-LRR) domains. We identified thirty four zinc finger domains in the R proteins of nine crops and were grouped into 19 types of zinc fingers. The size of individual zinc finger domain within the R genes varied from 11 to 84 amino acids, whereas the size of proteins containing these domains varied from 263 to 1305 amino acids. The biophysical analysis revealed that molecular weight of Pi54 zinc finger was lowest whereas the highest one was found in rice Pib zinc finger named as Transposes Transcription Factor (TTF). The instability (R(2) =0.95) and the aliphatic (R(2) =0.94) indices profile of zinc finger domains follows the polynomial distribution pattern. The pairwise identity analysis showed that the Lin11, Isl-1 & Mec-3 (LIM) zinc finger domain of rice blast resistance protein pi21 have 12.3% similarity with the nuclear transcription factor, X-box binding-like 1 (NFX) type zinc finger domain of Pi54 protein. For the first time, we reported that Pi54 (Pi-k(h)-Tetep), a rice blast resistance (R) protein have a small zinc finger domain of NFX type located on the C-terminal in between NBS and LRR domains of the R-protein. Compositional analysis depicted by the helical wheel diagram revealed the presence of a hydrophobic region within this domain which might help in exposing the LRR region for a possible R-Avr interaction. This domain is unique among all other cloned plant disease resistance genes and might play an important role in broad-spectrum nature of rice blast resistance gene Pi54.

Secretome analysis of Magnaporthe oryzae using in vitro systems

Magnaporthe oryzae is a devastating blast fungal pathogen of rice (Oryza sativa L.) that causes dramatic decreases in seed yield and quality. During the early stages of infection by this pathogen, the fungal spore senses the rice leaf surface, germinates, and penetrates the cell via an infectious structure known as an appressorium. During this process, M. oryzae secretes several proteins; however, these proteins are largely unknown mainly due to the lack of a suitable method for isolating secreted proteins during germination and appressoria formation. We examined the secretome of M. oryzae by mimicking the early stages of infection in vitro using a glass plate (GP), PVDF membrane, and liquid culture medium (LCM). Microscopic observation of M. oryzae growth revealed appressorium formation on the GP and PVDF membrane resembling natural M. oryzae-rice interactions; however, appresorium formation was not observed in the LCM. Secreted proteins were collected from the GP (3, 8, and 24 h), PVDF membrane (24 h), and LCM (48 h) and identified by two-dimensional gel electrophoresis (2DE) followed by tandem mass spectrometry. The GP, PVDF membrane, and LCM-derived 2D gels showed distinct protein patterns, indicating that they are complementary approaches. Collectively, 53 nonredundant proteins including previously known and novel secreted proteins were identified. Six biological functions were assigned to the proteins, with the predominant functional classes being cell wall modification, reactive oxygen species detoxification, lipid modification, metabolism, and protein modification. The in vitro system using GPs and PVDF membranes applied in this study to survey the M. oryzae secretome, can be used to further our understanding of the early interactions between M. oryzae and rice leaves.

[Transcriptome analysis of early interaction between rice and Magnaporthe oryzae using next-generation sequencing technology]

Magnaporthe oryzae is speculated to express and translocate a set of effector proteins to interfere with rice innate immunity and to promote its invasion. In this study, we carried out a transcriptome analysis of early interaction between rice and M. oryzae using a next generation sequencing technology to identify putative effector protein genes. We totally obtained 338 942 M. oryzae derived sequence tags from around 12.5 M gross mixed tags by referring the genome sequence of reference strain 70-15. These M. oryzae sequences were finally mapped to 779 predicted genes. Of these 799 genes, 108 were predicted to be involved in the interaction between rice and M. oryzae and 42 encoded proteins predicted to be secreted. Expression of these 42 putative secreted protein genes were further assessed using a reverse transcription polymerase chain reaction (RT-PCR) method, which revealed that 12 of them were expressed during the infection. Four of these 12 expressed genes were further confirmed to be expressed specifically during the infection. This study demonstrated an attempt to isolate M. oryzae genes, especially those putative effector protein genes, which are expressed specifically at the early interaction stage using a next-generation sequencing technology.

Plant defense mechanisms are activated during biotrophic and necrotrophic development of Colletotricum graminicola in maize

Hemibiotrophic plant pathogens first establish a biotrophic interaction with the host plant and later switch to a destructive necrotrophic lifestyle. Studies of biotrophic pathogens have shown that they actively suppress plant defenses after an initial microbe-associated molecular pattern-triggered activation. In contrast, studies of the hemibiotrophs suggest that they do not suppress plant defenses during the biotrophic phase, indicating that while there are similarities between the biotrophic phase of hemibiotrophs and biotrophic pathogens, the two lifestyles are not analogous. We performed transcriptomic, histological, and biochemical studies of the early events during the infection of maize (Zea mays) with Colletotrichum graminicola, a model pathosystem for the study of hemibiotrophy. Time-course experiments revealed that mRNAs of several defense-related genes, reactive oxygen species, and antimicrobial compounds all begin to accumulate early in the infection process and continue to accumulate during the biotrophic stage. We also discovered the production of maize-derived vesicular bodies containing hydrogen peroxide targeting the fungal hyphae. We describe the fungal respiratory burst during host infection, paralleled by superoxide ion production in specific fungal cells during the transition from biotrophy to a necrotrophic lifestyle. We also identified several novel putative fungal effectors and studied their expression during anthracnose development in maize. Our results demonstrate a strong induction of defense mechanisms occurring in maize cells during C. graminicola infection, even during the biotrophic development of the pathogen. We hypothesize that the switch to necrotrophic growth enables the fungus to evade the effects of the plant immune system and allows for full fungal pathogenicity.

Eukaryotic extracellular catalase-peroxidase from Magnaporthe grisea - Biophysical/chemical characterization of the first representative from a novel phytopathogenic KatG group

All phytopathogenic fungi have two catalase-peroxidase paralogues located either intracellularly (KatG1) or extracellularly (KatG2). Here, for the first time a secreted bifunctional, homodimeric catalase-peroxidase (KatG2 from the rice blast fungus Magnaporthe grisea) has been produced heterologously with almost 100% heme occupancy and comprehensively investigated by using a broad set of methods including UV-Vis, ECD and resonance Raman spectroscopy (RR), thin-layer spectroelectrochemistry, mass spectrometry, steady-state & presteady-state spectroscopy. RR spectroscopy reveals that MagKatG2 shows a unique mixed-spin state, non-planar heme b, and a proximal histidine with pronounced imidazolate character. At pH 7.0 and 25 °C, the standard reduction potential E°' of the Fe(III)/Fe(II) couple for the high-spin native protein was found to fall in the range typical for the KatG family. Binding of cyanide was relatively slow at pH 7.0 and 25 °C and with a K(d) value significantly higher than for the intracellular counterpart. Demonstrated by mass spectrometry MagKatG2 has the typical Trp118-Tyr251-Met277 adduct that is essential for its predominantly catalase activity at the unique acidic pH optimum. In addition, MagKatG2 acts as a versatile peroxidase using both one- and two-electron donors. Based on these data, structure-function relationships of extracellular eukaryotic KatGs are discussed with respect to intracellular KatGs and possible role(s) in host-pathogen interaction.

Identification of a sugar beet BvM14-MADS box gene through differential gene expression analysis of monosomic addition line M14

Monosomic addition line M14 carrying an additional chromosome 9 from Beta corolliflora Zosimovic ex Buttler was obtained through hybridization between the wild species B. corolliflora and a cultivated species Beta vulgaris L. var Saccharifera Alef. The M14 line showed diplosporic reproduction and stress tolerance. To identify differentially expressed genes in M14, a subtractive cDNA library was prepared by suppression subtractive hybridization (SSH) between M14 (2n=18+1) and B. vulgaris (2n=18). A total of 190 unique sequences were identified in the library and their putative functions were analyzed using Gene Ontology (GO). One of the genes, designated as BvM14-MADS box, encodes a MADS box transcription factor. It was cloned from M14 and over-expressed in transgenic tobacco plants. Interestingly, this gene was located on chromosome 2 of B. vulgaris, not on the additional chromosome 9. Overexpression of BvM14-MADS box led to significant phenotypic changes in tobacco. The differential expression of BvM14-MADS box gene in M14 may be caused by the interaction between the additional chromosome 9 from B. corolliflora and the B. vulgaris chromosomes in M14.

Comparative secretome investigation of Magnaporthe oryzae proteins responsive to nitrogen starvation

Magnaporthe oryzae is a fungal pathogen that causes blast disease in rice. During its early infection process, during which starvation of nutrients, including nitrogen, prevails before establishment of successful infection, the fungally secreted proteins play an important role in the pathogenicity and stress response. In this study, M. oryzae-secreted proteins were investigated in an N-deficient minimal medium using two-dimensional gel electrophoresis (2-DGE) coupled with mass spectrometry analysis (MALDI-TOF-MS and μLC-ESI-MS/MS). The 2-DGE analysis of secreted proteins detected 89 differentially expressed protein spots (14 downregulated and 75 upregulated) responsive to N starvation. Eighty five of the protein spots were identified by mass spectrometry analyses. Identified proteins were mainly cell wall hydrolase enzymes (22.4%), protein and lipid hydrolases (24.7%), reactive oxygen species detoxifying proteins (22.4%), and proteins with unknown function (14.1%), suggesting early production of prerequisite proteins for successful infection of the host. SignalP analysis predicted the presence of signal peptides in 67% of the identified proteins, suggesting that in addition to the classical Golgi/endoplasmic reticulum secretory pathway, M. oryzae might possess other, as yet undefined, secretory pathways. Those nonclassical or leaderless secretion proteins accounted for 25.9% of the total identified proteins by TatP and SecretomeP predictions. Semiquantitative reverse transcriptase polymerase chain reaction of seven randomly selected N-responsive secreted proteins also revealed a good correlation between RNA and protein levels. Taken together, the establishment of the M. oryzae secretome that is responsive to N starvation provides the first evidence of the secretion of 60 unreported and 25 previously known proteins. This developed protein inventory could be exploited to improve our understanding of the secretory mechanisms of M. oryzae and its invasive growth process in rice tissue.

Unifying themes in microbial associations with animal and plant hosts described using the gene ontology

Microbes form intimate relationships with hosts (symbioses) that range from mutualism to parasitism. Common microbial mechanisms involved in a successful host association include adhesion, entry of the microbe or its effector proteins into the host cell, mitigation of host defenses, and nutrient acquisition. Genes associated with these microbial mechanisms are known for a broad range of symbioses, revealing both divergent and convergent strategies. Effective comparisons among these symbioses, however, are hampered by inconsistent descriptive terms in the literature for functionally similar genes. Bioinformatic approaches that use homology-based tools are limited to identifying functionally similar genes based on similarities in their sequences. An effective solution to these limitations is provided by the Gene Ontology (GO), which provides a standardized language to describe gene products from all organisms. The GO comprises three ontologies that enable one to describe the molecular function(s) of gene products, the biological processes to which they contribute, and their cellular locations. Beginning in 2004, the Plant-Associated Microbe Gene Ontology (PAMGO) interest group collaborated with the GO consortium to extend the GO to accommodate terms for describing gene products associated with microbe-host interactions. Currently, over 900 terms that describe biological processes common to diverse plant- and animal-associated microbes are incorporated into the GO database. Here we review some unifying themes common to diverse host-microbe associations and illustrate how the new GO terms facilitate a standardized description of the gene products involved. We also highlight areas where new terms need to be developed, an ongoing process that should involve the whole community.