Distinct dynamics of the nucleolus in response to nutrient availability and during development in the rice blast fungus

IMPORTANCE

The nucleolus is a dynamic subnuclear structure that is involved in many fundamental processes of the nucleus. In higher eukaryotic cells, the size and shape of nucleoli correlate with nucleolar activities. For fungi, knowledge of the nucleolus and its functions is primarily gleaned from budding yeast. Whether such correlation is conserved and how nucleolar functions are regulated in filamentous fungi including important human and crop pathogens are largely unknown. Our observations reveal that the dynamics of nucleolus in a model plant pathogenic fungus, Magnaporthe oryzae, is distinct from those of animal and yeast nucleoli under low nutrient availability and during pathogenic development. Our data not only provide new insight into the nucleoli in filamentous fungi but also highlight the need for investigating how nucleolar dynamics is regulated in comparison to other eukaryotes.

Tissues and mechanisms associated with Verticillium wilt resistance in tomato using bi-grafted near-isogenic lines

Host resistance is the primary means to control Verticillium dahliae, a soil-borne pathogen causing major losses on a broad range of plants, including tomato. The tissues and mechanisms responsible for resistance remain obscure. In the field, resistant tomato used as rootstocks does not confer resistance. Here, we created bi-grafted plants with near-isogenic lines (NILs) exhibiting (Ve1) or lacking (ve1) resistance to V. dahliae race 1. Ten days after inoculation, scion and rootstock tissues were subjected to differential gene expression and co-expression network analyses. Symptoms only developed in susceptible scions regardless of the rootstock. Infection caused more dramatic alteration of tomato gene expression in susceptible compared with resistant tissues, including pathogen receptor, signaling pathway, pathogenesis-related protein, and cell wall modification genes. Differences were observed between scions and rootstocks, primarily related to physiological processes in these tissues. Gene expression in scions was influenced by the rootstock genotype. A few genes were associated with the Ve1 genotype, which was independent of infection or tissue type. Several were physically clustered, some near the Ve1 locus on chromosome 9. Transcripts mapped to V. dahliae were dominated by secreted candidate effector proteins. These findings advance knowledge of molecular mechanisms underlying the tomato-V. dahliae interaction.

Host induced gene silencing of Magnaporthe oryzae by targeting pathogenicity and development genes to control rice blast disease

Rice blast disease caused by the hemi-biotrophic fungus Magnaporthe oryzae is the most destructive disease of rice world-wide. Traditional disease resistance strategies for the control of rice blast disease have not proved durable. HIGS (host induced gene silencing) is being developed as an alternative strategy. Six genes (CRZ1, PMC1, MAGB, LHS1, CYP51A, CYP51B) that play important roles in pathogenicity and development of M. oryzae were chosen for HIGS. HIGS vectors were transformed into rice calli through Agrobacterium-mediated transformation and T0, T1 and T2 generations of transgenic rice plants were generated. Except for PMC1 and LHS1, HIGS transgenic rice plants challenged with M. oryzae showed significantly reduced disease compared with non-silenced control plants. Following infection with M. oryzae of HIGS transgenic plants, expression levels of target genes were reduced as demonstrated by Quantitative RT-PCR. In addition, treating M. oryzae with small RNA derived from the target genes inhibited fungal growth. These findings suggest RNA silencing signals can be transferred from host to an invasive fungus and that HIGS has potential to generate resistant rice against M. oryzae.

Soybean Cyst Nematodes Influence Aboveground Plant Volatile Signals Prior to Symptom Development

Soybean cyst nematode (SCN), Heterodera glycines, is one of the most destructive soybean pests worldwide. Unlike many diseases, SCN doesn't show above ground evidence of disease until several weeks after infestation. Knowledge of Volatile Organic Compounds (VOCs) related to pests and pathogens of foliar tissue is extensive, however, information related to above ground VOCs in response to root damage is lacking. In temporal studies, gas chromatography-mass spectrometry analysis of VOCs from the foliar tissues of SCN infested plants yielded 107 VOCs, referred to as Common Plant Volatiles (CPVs), 33 with confirmed identities. Plants showed no significant stunting until 10 days after infestation. Total CPVs increased over time and were significantly higher from SCN infested plants compared to mock infested plants post 7 days after infestation (DAI). Hierarchical clustering analysis of expression ratios (SCN: Mock) across all time points revealed 5 groups, with the largest group containing VOCs elevated in response to SCN infestation. Linear projection of Principal Component Analysis clearly separated SCN infested from mock infested plants at time points 5, 7, 10 and 14 DAI. Elevated Styrene (CPV11), D-Limonene (CPV32), Tetradecane (CPV65), 2,6-Di-T-butyl-4-methylene-2,5-cyclohexadiene-1-one (CPV74), Butylated Hydroxytoluene (CPV76) and suppressed Ethylhexyl benzoate (CPV87) levels, were associated with SCN infestation prior to stunting. Our findings demonstrate that SCN infestation elevates the release of certain VOCs from foliage and that some are evident prior to symptom development. VOCs associated with SCN infestations prior to symptom development may be valuable for innovative diagnostic approaches.

Comparative Genome Analyses of 18 Verticillium dahliae Tomato Isolates Reveals Phylogenetic and Race Specific Signatures

Host resistance is one of the few strategies available to combat the soil borne pathogenic fungus Verticillium dahliae. Understanding pathogen diversity in populations is key to successfully deploying host resistance. In this study the genomes of 18 V. dahliae isolates of races 1 (n = 2), 2 (n = 4), and 3 (n = 12) from Japan, California, and North Carolina were sequenced and mapped to the reference genome of JR2 (from tomato). The genomes were analyzed for phylogenetic and pathogen specific signatures to classify specific strains or genes for future research. Four highly clonal lineages/groups were discovered, including a lineage unique to North Carolina isolates, which had the rare MAT1-1 mating type. No evidence for recombination between isolates of different mating types was observed, even in isolates of different mating types discovered in the same field. By mapping these 18 isolates genomes to the JR2 reference genome, 193 unique candidate effectors were found using SignalP and EffectorP. Within these effectors, 144 highly conserved effectors, 42 mutable effectors (truncated or present in some isolates but absent in others), and 7 effectors present in highly variable regions of the chromosomes were discovered. Of the 144 core effectors, 21 were highly conserved in V. alfalfae and V. longisporum, 7 of which have no known function. Within the non-core effectors 30 contained large numbers of non-synonymous mutations, while 15 of them contained indels, frameshift mutations, or were present on highly variable regions of the chromosome. Two of these highly variable region effectors (HVREs) were only present in race 2 isolates, but not in race 3 isolates. The race 1 effector Ave1 was also present in a highly variable region. These data may suggest that these highly variable regions are enriched in race determinant genes, consistent with the two-speed genome hypothesis.

Opportunities and Challenges in Studies of Host-Pathogen Interactions and Management of Verticillium dahliae in Tomatoes

Tomatoes (Solanum lycopersicum L.) are a valuable horticultural crop that are grown and consumed worldwide. Optimal production is hindered by several factors, among which Verticillium dahliae, the cause of Verticillium wilt, is considered a major biological constraint in temperate production regions. V. dahliae is difficult to mitigate because it is a vascular pathogen, has a broad host range and worldwide distribution, and can persist in soil for years. Understanding pathogen virulence and genetic diversity, host resistance, and plant-pathogen interactions could ultimately inform the development of integrated strategies to manage the disease. In recent years, considerable research has focused on providing new insights into these processes, as well as the development and integration of environment-friendly management approaches. Here, we discuss the current knowledge on the race and population structure of V. dahliae, including pathogenicity factors, host genes, proteins, enzymes involved in defense, and the emergent management strategies and future research directions for managing Verticillium wilt in tomatoes.

Dynamic Changes in the Microbiome of Rice During Shoot and Root Growth Derived From Seeds

Microbes form close associations with host plants including rice as both surface (epiphytes) and internal (endophytes) inhabitants. Yet despite rice being one of the most important cereal crops agriculturally and economically, knowledge of its microbiome, particularly core inhabitants and any functional properties bestowed is limited. In this study, the microbiome in rice seedlings derived directly from seeds was identified, characterized and compared to the microbiome of the seed. Rice seeds were sourced from two different locations in Arkansas, USA of two different rice genotypes (Katy, M202) from two different harvest years (2013, 2014). Seeds were planted in sterile media and bacterial as well as fungal communities were identified through 16S and ITS sequencing, respectively, for four seedling compartments (root surface, root endosphere, shoot surface, shoot endosphere). Overall, 966 bacterial and 280 fungal ASVs were found in seedlings. Greater abundance and diversity were detected for the microbiome associated with roots compared to shoots and with more epiphytes than endophytes. The seedling compartments were the driving factor for microbial community composition rather than other factors such as rice genotype, location and harvest year. Comparison with datasets from seeds revealed that 91 (out of 296) bacterial and 11 (out of 341) fungal ASVs were shared with seedlings with the majority being retained within root tissues. Core bacterial and fungal microbiome shared across seedling samples were identified. Core bacteria genera identified in this study such as Rhizobium, Pantoea, Sphingomonas, and Paenibacillus have been reported as plant growth promoting bacteria while core fungi such as Pleosporales, Alternaria and Occultifur have potential as biocontrol agents.

Comparative Proteomic Analysis of Wild Type and Mutant Lacking an SCF E3 Ligase F-Box Protein in Magnaporthe oryzae

Magnaporthe oryzae (M. oryzae) is a pathogenic, filamentous fungus that is a primary cause of rice blast disease. The M. oryzae protein MGG_13065, SCF E3 ubiquitin ligase complex F-box protein, has been identified as playing a crucial role in the infection process, specifically, as part of the ubiquitin mediated proteolysis pathway. Proteins targeted by MGG_13065 E3 ligase are first phosphorylated and then ubiquitinated by E3 ligase. In this study, we used a label-free quantitative global proteomics technique to probe the role of ubiquitination and phosphorylation in the mechanism of how E3 ligase regulates change in virulence of M. oryzae. To do this, we compared the WT M. oryzae 70-15 strain with a gene knock out (E3 ligase KO) strain. After applying a ≥ 5 normalized spectral count cutoff, a total of 4432 unique proteins were identified comprised of 4360 and 4372 in the WT and E3 ligase KO samples, respectively. Eighty proteins drastically increased in abundance, while 65 proteins decreased in abundance in the E3 ligase KO strain. Proteins (59) were identified only in the WT strain; 13 of these proteins had both phosphorylation and ubiquitination post-translational modifications. Proteins (71) were revealed to be only in the E3 ligase KO strain; 23 of the proteins have both phosphorylation and ubiquitination post-translational modifications. Several of these proteins were associated with key biological processes. These data greatly assist in the selection of future genes for functional studies and enable mechanistic insight related to virulence.

Movement of small RNAs in and between plants and fungi

RNA interference is a biological process whereby small RNAs inhibit gene expression through neutralizing targeted mRNA molecules. This process is conserved in eukaryotes. Here, recent work regarding the mechanisms of how small RNAs move within and between organisms is examined. Small RNAs can move locally and systemically in plants through plasmodesmata and phloem, respectively. In fungi, transportation of small RNAs may also be achieved by septal pores and vesicles. Recent evidence also supports bidirectional cross-kingdom communication of small RNAs between host plants and adapted fungal pathogens to affect the outcome of infection. We discuss several mechanisms for small RNA trafficking and describe evidence for transport through naked form, combined with RNA-binding proteins or enclosed by vesicles.

Comparative proteomic analysis between nitrogen supplemented and starved conditions in Magnaporthe oryzae

BACKGROUND

Fungi are constantly exposed to nitrogen limiting environments, and thus the efficient regulation of nitrogen metabolism is essential for their survival, growth, development and pathogenicity. To understand how the rice blast pathogen Magnaporthe oryzae copes with limited nitrogen availability, a global proteome analysis under nitrogen supplemented and nitrogen starved conditions was completed.

METHODS

M. oryzae strain 70-15 was cultivated in liquid minimal media and transferred to media with nitrate or without a nitrogen source. Proteins were isolated and subjected to unfractionated gel-free based liquid chromatography-tandem mass spectrometry (LC-MS/MS). The subcellular localization and function of the identified proteins were predicted using bioinformatics tools.

RESULTS

A total of 5498 M. oryzae proteins were identified. Comparative analysis of protein expression showed 363 proteins and 266 proteins significantly induced or uniquely expressed under nitrogen starved or nitrogen supplemented conditions, respectively. A functional analysis of differentially expressed proteins revealed that during nitrogen starvation nitrogen catabolite repression, melanin biosynthesis, protein degradation and protein translation pathways underwent extensive alterations. In addition, nitrogen starvation induced accumulation of various extracellular proteins including small extracellular proteins consistent with observations of a link between nitrogen starvation and the development of pathogenicity in M. oryzae.

CONCLUSION

The results from this study provide a comprehensive understanding of fungal responses to nitrogen availability.

Identification and characterization of suppressors of plant cell death (SPD) effectors from Magnaporthe oryzae

Phytopathogenic microorganisms, including the fungal pathogen Magnaporthe oryzae, secrete a myriad of effector proteins to facilitate infection. Utilizing the transient expression of candidate effectors in the leaves of the model plant Nicotiana benthamiana, we identified 11 suppressors of plant cell death (SPD) effectors from M. oryzae that were able to block the host cell death reaction induced by Nep1. Ten of these 11 were also able to suppress BAX-mediated plant cell death. Five of the 11 SPD genes have been identified previously as either essential for the pathogenicity of M. oryzae, secreted into the plant during disease development, or as suppressors or homologues of other characterized suppressors. In addition, of the remaining six, we showed that SPD8 (previously identified as BAS162) was localized to the rice cytoplasm in invaded and surrounding uninvaded cells during biotrophic invasion. Sequence analysis of the 11 SPD genes across 43 re-sequenced M. oryzae genomes revealed that SPD2, SPD4 and SPD7 have nucleotide polymorphisms amongst the isolates. SPD4 exhibited the highest level of nucleotide diversity of any currently known effector from M. oryzae in addition to the presence/absence polymorphisms, suggesting that this gene is potentially undergoing selection to avoid recognition by the host. Taken together, we have identified a series of effectors, some of which were previously unknown or whose function was unknown, that probably act at different stages of the infection process and contribute to the virulence of M. oryzae.

Examining ubiquitinated peptide enrichment efficiency through an epitope labeled protein

Ubiquitination is a dynamic process that is responsible for regulation of cellular responses to stimuli in a number of biological systems. Previous efforts to study this post-translational modification have focused on protein enrichment; however, recent research utilizes the presence of the di-glycine (Gly-Gly) remnants following trypsin digestion to immuno-enrich ubiquitinated peptides. Monoclonal antibodies developed to the cleaved ubiquitin modification epitope, (tert-butoxycarbonyl) glycylglycine (Boc-Gly-Gly-NHS)(1), are used to identify the Gly-Gly signature. Here, we have successfully generated the Boc-Gly-Gly-NHS modification and showed that when conjugated to a lysine containing protein, such as lysozyme, it can be applied as a standard protein to examine ubiquitinated peptide enrichment within a complex background.

Comparative genome analysis and genome evolution of members of the magnaporthaceae family of fungi

BACKGROUND

Magnaporthaceae, a family of ascomycetes, includes three fungi of great economic importance that cause disease in cereal and turf grasses: Magnaporthe oryzae (rice blast), Gaeumannomyces graminis var. tritici (take-all disease), and Magnaporthe poae (summer patch disease). Recently, the sequenced and assembled genomes for these three fungi were reported. Here, the genomes were compared for orthologous genes in order to identified genes that are unique to the Magnaporthaceae family of fungi. In addition, ortholog clustering was used to identify a core proteome for the Magnaporthaceae, which was examined for diversifying and purifying selection and evidence of two-speed genome evolution.

RESULTS

A genome-scale comparative study was conducted across 74 fungal genomes to identify clusters of orthologous genes unique to the three Magnaporthaceae species as well as species specific genes. We found 1149 clusters that were unique to the Magnaporthaceae family of fungi with 295 of those containing genes from all three species. Gene clusters involved in metabolic and enzymatic activities were highly represented in the Magnaporthaceae specific clusters. Also highly represented in the Magnaporthaceae specific clusters as well as in the species specific genes were transcriptional regulators. In addition, we examined the relationship between gene evolution and distance to repetitive elements found in the genome. No correlations between diversifying or purifying selection and distance to repetitive elements or an increased rate of evolution in secreted and small secreted proteins were observed.

CONCLUSIONS

Taken together, these data show that at the genome level, there is no evidence to suggest multi-speed genome evolution or that proximity to repetitive elements play a role in diversification of genes.

Phosphoproteome Analysis Links Protein Phosphorylation to Cellular Remodeling and Metabolic Adaptation during Magnaporthe oryzae Appressorium Development

The rice pathogen, Magnaporthe oryzae, undergoes a complex developmental process leading to formation of an appressorium prior to plant infection. In an effort to better understand phosphoregulation during appressorium development, a mass spectrometry based phosphoproteomics study was undertaken. A total of 2924 class I phosphosites were identified from 1514 phosphoproteins from mycelia, conidia, germlings, and appressoria of the wild type and a protein kinase A (PKA) mutant. Phosphoregulation during appressorium development was observed for 448 phosphosites on 320 phosphoproteins. In addition, a set of candidate PKA targets was identified encompassing 253 phosphosites on 227 phosphoproteins. Network analysis incorporating regulation from transcriptomic, proteomic, and phosphoproteomic data revealed new insights into the regulation of the metabolism of conidial storage reserves and phospholipids, autophagy, actin dynamics, and cell wall metabolism during appressorium formation. In particular, protein phosphorylation appears to play a central role in the regulation of autophagic recycling and actin dynamics during appressorium formation. Changes in phosphorylation were observed in multiple components of the cell wall integrity pathway providing evidence that this pathway is highly active during appressorium development. Several transcription factors were phosphoregulated during appressorium formation including the bHLH domain transcription factor MGG_05709. Functional analysis of MGG_05709 provided further evidence for the role of protein phosphorylation in regulation of glycerol metabolism and the metabolic reprogramming characteristic of appressorium formation. The data presented here represent a comprehensive investigation of the M. oryzae phosphoproteome and provide key insights on the role of protein phosphorylation during infection-related development.

Genome-wide profiling of DNA methylation provides insights into epigenetic regulation of fungal development in a plant pathogenic fungus, Magnaporthe oryzae

DNA methylation is an important epigenetic modification that regulates development of plants and mammals. To investigate the roles of DNA methylation in fungal development, we profiled genome-wide methylation patterns at single-nucleotide resolution during vegetative growth, asexual reproduction, and infection-related morphogenesis in a model plant pathogenic fungus, Magnaporthe oryzae. We found that DNA methylation occurs in and around genes as well as transposable elements and undergoes global reprogramming during fungal development. Such reprogramming of DNA methylation suggests that it may have acquired new roles other than controlling the proliferation of TEs. Genetic analysis of DNA methyltransferase deletion mutants also indicated that proper reprogramming in methylomes is required for asexual reproduction in the fungus. Furthermore, RNA-seq analysis showed that DNA methylation is associated with transcriptional silencing of transposable elements and transcript abundance of genes in context-dependent manner, reinforcing the role of DNA methylation as a genome defense mechanism. This comprehensive approach suggests that DNA methylation in fungi can be a dynamic epigenetic entity contributing to fungal development and genome defense. Furthermore, our DNA methylomes provide a foundation for future studies exploring this key epigenetic modification in fungal development and pathogenesis.

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.

Experimental evolution reveals genome-wide spectrum and dynamics of mutations in the rice blast fungus, Magnaporthe oryzae

Knowledge on mutation processes is central to interpreting genetic analysis data as well as understanding the underlying nature of almost all evolutionary phenomena. However, studies on genome-wide mutational spectrum and dynamics in fungal pathogens are scarce, hindering our understanding of their evolution and biology. Here, we explored changes in the phenotypes and genome sequences of the rice blast fungus Magnaporthe oryzae during the forced in vitro evolution by weekly transfer of cultures on artificial media. Through combination of experimental evolution with high throughput sequencing technology, we found that mutations accumulate rapidly prior to visible phenotypic changes and that both genetic drift and selection seem to contribute to shaping mutational landscape, suggesting the buffering capacity of fungal genome against mutations. Inference of mutational effects on phenotypes through the use of T-DNA insertion mutants suggested that at least some of the DNA sequence mutations are likely associated with the observed phenotypic changes. Furthermore, our data suggest oxidative damages and UV as major sources of mutation during subcultures. Taken together, our work revealed important properties of original source of variation in the genome of the rice blast fungus. We believe that these results provide not only insights into stability of pathogenicity and genome evolution in plant pathogenic fungi but also a model in which evolution of fungal pathogens in natura can be comparatively investigated.

In-depth analysis of the Magnaporthe oryzae conidial proteome

The filamentous fungus Magnaporthe oryzae (M. oryzae) is the causative agent of rice blast disease and presents a significant threat to worldwide rice production. To establish the groundwork for future research on the pathogenic development of M. oryzae, a global proteomic study of conidia was performed. The filter aided sample preparation method (FASP) and anion StageTip fractionation combined with long, optimized shallow 210 min nanoLC gradients prior to mass spectrometry analysis on an Orbitrap XL was applied, which resulted in a doubling of protein identifications in comparison to our previous GeLC analysis. Herein, we report the identification of 2912 conidial proteins at a 1% protein false discovery rate (FDR) and we present the most extensive study performed on M. oryzae conidia to date. A similar distribution between identified proteins and the predicted proteome was observed when subcellular localization analysis was performed, suggesting the detected proteins build a representative portion of the predicted proteome. A higher percentage of cytoplasmic proteins (associated with translation, energy, and metabolism) were observed in the conidial proteome relative to the whole predicted proteome. Conversely, nuclear and extracellular proteins were less well represented in the conidial proteome. Further analysis by gene ontology revealed biological insights into identified proteins important for central metabolic processes and the physiology of conidia.

Accurate discrimination of bHLH domains in plants, animals, and fungi using biologically meaningful sites

Background

The highly conserved bHLH (basic Helix-Loop-Helix) domain, found in many transcription factors, has been well characterized separately in Plants, Animals, and Fungi. While conserved, even functionally constrained sites have varied since the Eukarya split. Our research identifies those slightly variable sites that were highly characteristic of Plants, Animals, or Fungi.

Results

Through discriminant analysis, we identified five highly discerning DNA-binding amino acid sites. Additionally, by incorporating Kingdom specific HMMs, we were able to construct a tool to quickly and accurately identify and classify bHLH sequences using these sites.

Conclusions

We conclude that highly discerning sites identified through our analysis were likely under functional constraints specific to each Kingdom. We also demonstrated the utility of our tool by identifying and classifying previously unknown bHLH domains in both characterized genomes and from sequences in a large environmental sample.

Polyubiquitin is required for growth, development and pathogenicity in the rice blast fungus Magnaporthe oryzae

Protein ubiquitination, which is highly selective, regulates many important biological processes including cellular differentiation and pathogenesis in eukaryotic cells. Here, we integrated pharmacological, molecular and proteomic approaches to explore the role of ubiquitination in Magnaporthe oryzae, the leading fungal disease of rice world-wide. Inhibition of ubiquitin-mediated proteolysis using the 26S proteasome inhibitor, Bortezomib, significantly attenuated conidia germination, appressorium formation and pathogenicity in M. oryzae. Gene expression analysis revealed that many genes associated with protein ubiquitination were developmentally regulated during conidia germination. Only a few, including a polyubiquitin encoding gene, MGG_01282, were more abundantly expressed during appressorium formation and under nitrogen starvation. Targeted gene deletion of MGG_01282, in addition to a significant reduction in protein ubiquitination as determined by immuno blot assays, resulted in pleiotropic effects on M. oryzae including reduced growth and sporulation, abnormal conidia morphology, reduced germination and appressorium formation, and the inability to cause disease. Mutants were also defective in sexual development and were female sterile. Using mass spectrometry, we identified 63 candidate polyubiquitinated proteins under nitrogen starvation, which included overrepresentation of proteins involved in translation, transport and protein modification. Our study suggests that ubiquitination of target proteins plays an important role in nutrient assimilation, development and pathogenicity of M. oryzae.

Comparative analysis of the genomes of two field isolates of the rice blast fungus Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae is one of the most destructive diseases of rice worldwide. The fungal pathogen is notorious for its ability to overcome host resistance. To better understand its genetic variation in nature, we sequenced the genomes of two field isolates, Y34 and P131. In comparison with the previously sequenced laboratory strain 70-15, both field isolates had a similar genome size but slightly more genes. Sequences from the field isolates were used to improve genome assembly and gene prediction of 70-15. Although the overall genome structure is similar, a number of gene families that are likely involved in plant-fungal interactions are expanded in the field isolates. Genome-wide analysis on asynonymous to synonymous nucleotide substitution rates revealed that many infection-related genes underwent diversifying selection. The field isolates also have hundreds of isolate-specific genes and a number of isolate-specific gene duplication events. Functional characterization of randomly selected isolate-specific genes revealed that they play diverse roles, some of which affect virulence. Furthermore, each genome contains thousands of loci of transposon-like elements, but less than 30% of them are conserved among different isolates, suggesting active transposition events in M. oryzae. A total of approximately 200 genes were disrupted in these three strains by transposable elements. Interestingly, transposon-like elements tend to be associated with isolate-specific or duplicated sequences. Overall, our results indicate that gain or loss of unique genes, DNA duplication, gene family expansion, and frequent translocation of transposon-like elements are important factors in genome variation of the rice blast fungus.

Comparative analysis of the genomes of two field isolates of the rice blast fungus Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae is one of the most destructive diseases of rice worldwide. The fungal pathogen is notorious for its ability to overcome host resistance. To better understand its genetic variation in nature, we sequenced the genomes of two field isolates, Y34 and P131. In comparison with the previously sequenced laboratory strain 70-15, both field isolates had a similar genome size but slightly more genes. Sequences from the field isolates were used to improve genome assembly and gene prediction of 70-15. Although the overall genome structure is similar, a number of gene families that are likely involved in plant-fungal interactions are expanded in the field isolates. Genome-wide analysis on asynonymous to synonymous nucleotide substitution rates revealed that many infection-related genes underwent diversifying selection. The field isolates also have hundreds of isolate-specific genes and a number of isolate-specific gene duplication events. Functional characterization of randomly selected isolate-specific genes revealed that they play diverse roles, some of which affect virulence. Furthermore, each genome contains thousands of loci of transposon-like elements, but less than 30% of them are conserved among different isolates, suggesting active transposition events in M. oryzae. A total of approximately 200 genes were disrupted in these three strains by transposable elements. Interestingly, transposon-like elements tend to be associated with isolate-specific or duplicated sequences. Overall, our results indicate that gain or loss of unique genes, DNA duplication, gene family expansion, and frequent translocation of transposon-like elements are important factors in genome variation of the rice blast fungus.

Magnaporthe oryzae is the causal agent of rice blast that is mainly controlled with resistance cultivars. However, genetic variations in the pathogen often lead to overcoming R gene-mediated resistance in rice cultivars. In this study we sequenced two field isolates from China and Japan. In comparison with the laboratory strain that was previously sequenced, the field isolates have a similar genome size and overall genome structure. However, they have slightly more genes and contain a number of expanded gene families that are likely involved in plant-fungal interactions. Each of the isolates has specific genes, some of which affect virulence and some others are important for asexual development. The three strains differ noticeably in the distribution of transposon-like elements. Many of the transposable elements tend to be associated with isolate-specific or duplicated sequences. This study revealed genetic factors involved in genome variation of the rice blast fungus.

Host-induced gene silencing: a tool for understanding fungal host interaction and for developing novel disease control strategies

Recent discoveries regarding small RNAs and the mechanisms of gene silencing are providing new opportunities to explore fungal pathogen-host interactions and potential strategies for novel disease control. Plant pathogenic fungi are a constant and major threat to global food security; they represent the largest group of disease-causing agents on crop plants on the planet. An initial understanding of RNA silencing mechanisms and small RNAs was derived from model fungi. Now, new knowledge with practical implications for RNA silencing is beginning to emerge from the study of plant-fungus interactions. Recent studies have shown that the expression of silencing constructs in plants designed on fungal genes can specifically silence their targets in invading pathogenic fungi, such as Fusarium verticillioides, Blumeria graminis and Puccinia striiformis f.sp. tritici. Here, we highlight the important general aspects of RNA silencing mechanisms and emphasize recent findings from plant pathogenic fungi. Strategies to employ RNA silencing to investigate the basis of fungal pathogenesis are discussed. Finally, we address important aspects for the development of fungal-derived resistance through the expression of silencing constructs in host plants as a powerful strategy to control fungal disease.

Evaluation of normalization methods on GeLC-MS/MS label-free spectral counting data to correct for variation during proteomic workflows

Normalization of spectral counts (SpCs) in label-free shotgun proteomic approaches is important to achieve reliable relative quantification. Three different SpC normalization methods, total spectral count (TSpC) normalization, normalized spectral abundance factor (NSAF) normalization, and normalization to selected proteins (NSP) were evaluated based on their ability to correct for day-to-day variation between gel-based sample preparation and chromatographic performance. Three spectral counting data sets obtained from the same biological conidia sample of the rice blast fungus Magnaporthe oryzae were analyzed by 1D gel and liquid chromatography-tandem mass spectrometry (GeLC-MS/MS). Equine myoglobin and chicken ovalbumin were spiked into the protein extracts prior to 1D-SDS- PAGE as internal protein standards for NSP. The correlation between SpCs of the same proteins across the different data sets was investigated. We report that TSpC normalization and NSAF normalization yielded almost ideal slopes of unity for normalized SpC versus average normalized SpC plots, while NSP did not afford effective corrections of the unnormalized data. Furthermore, when utilizing TSpC normalization prior to relative protein quantification, t-testing and fold-change revealed the cutoff limits for determining real biological change to be a function of the absolute number of SpCs. For instance, we observed the variance decreased as the number of SpCs increased, which resulted in a higher propensity for detecting statistically significant, yet artificial, change for highly abundant proteins. Thus, we suggest applying higher confidence level and lower fold-change cutoffs for proteins with higher SpCs, rather than using a single criterion for the entire data set. By choosing appropriate cutoff values to maintain a constant false positive rate across different protein levels (i.e., SpC levels), it is expected this will reduce the overall false negative rate, particularly for proteins with higher SpCs.

Phylogenetic analysis and classification of the fungal bHLH domain

The basic Helix-Loop-Helix (bHLH) domain is an essential highly conserved DNA-binding domain found in many transcription factors in all eukaryotic organisms. The bHLH domain has been well studied in the Animal and Plant Kingdoms but has yet to be characterized within Fungi. Herein, we obtained and evaluated the phylogenetic relationship of 490 fungal-specific bHLH containing proteins from 55 whole genome projects composed of 49 Ascomycota and 6 Basidiomycota organisms. We identified 12 major groupings within Fungi (F1-F12); identifying conserved motifs and functions specific to each group. Several classification models were built to distinguish the 12 groups and elucidate the most discerning sites in the domain. Performance testing on these models, for correct group classification, resulted in a maximum sensitivity and specificity of 98.5% and 99.8%, respectively. We identified 12 highly discerning sites and incorporated those into a set of rules (simplified model) to classify sequences into the correct group. Conservation of amino acid sites and phylogenetic analyses established that like plant bHLH proteins, fungal bHLH-containing proteins are most closely related to animal Group B. The models used in these analyses were incorporated into a software package, the source code for which is available at www.fungalgenomics.ncsu.edu.

Comparison of stable-isotope labeling with amino acids in cell culture and spectral counting for relative quantification of protein expression

Protein quantification is one of the principal goals of mass spectrometry (MS)-based proteomics, and many strategies exist to achieve it. Several approaches involve the incorporation of a stable-isotope label using either chemical derivatization, enzymatically catalyzed incorporation of (18)O, or metabolic labeling in a cell or tissue culture. These techniques can be cost or time prohibitive or not amenable to the biological system of interest. Label-free techniques including those utilizing integrated ion abundance and spectral counting offer an alternative to stable-isotope-based methodologies. Herein, we present the comparison of stable-isotope labeling of amino acids in cell culture (SILAC) with spectral counting for the quantification of human embryonic stem cells as they differentiate toward the trophectoderm at three time points. Our spectral counting experimental strategy resulted in the identification of 2641 protein groups across three time points with an average sequence coverage of 30.3%, of which 1837 could be quantified with more than five spectral counts. SILAC quantification was able to identify 1369 protein groups with an average coverage of 24.7%, of which 1027 could be quantified across all time points. Within this context we further explore the capacity of each strategy for proteome coverage, variation in quantification, and the relative sensitivity of each technique to the detection of change in relative protein expression.

Diverse and tissue-enriched small RNAs in the plant pathogenic fungus, Magnaporthe oryzae

Background

Emerging knowledge of the impact of small RNAs as important cellular regulators has prompted an explosion of small transcriptome sequencing projects. Although significant progress has been made towards small RNA discovery and biogenesis in higher eukaryotes and other model organisms, knowledge in simple eukaryotes such as filamentous fungi remains limited.

Results

Here, we used 454 pyrosequencing to present a detailed analysis of the small RNA transcriptome (~ 15 - 40 nucleotides in length) from mycelia and appressoria tissues of the rice blast fungal pathogen, Magnaporthe oryzae. Small RNAs mapped to numerous nuclear and mitochondrial genomic features including repetitive elements, tRNA loci, rRNAs, protein coding genes, snRNAs and intergenic regions. For most elements, small RNAs mapped primarily to the sense strand with the exception of repetitive elements to which small RNAs mapped in the sense and antisense orientation in near equal proportions. Inspection of the small RNAs revealed a preference for U and suppression of C at position 1, particularly for antisense mapping small RNAs. In the mycelia library, small RNAs of the size 18 - 23 nt were enriched for intergenic regions and repetitive elements. Small RNAs mapping to LTR retrotransposons were classified as LTR retrotransposon-siRNAs (LTR-siRNAs). Conversely, the appressoria library had a greater proportion of 28 - 35 nt small RNAs mapping to tRNA loci, and were classified as tRNA-derived RNA fragments (tRFs). LTR-siRNAs and tRFs were independently validated by 3' RACE PCR and northern blots, respectively.

Conclusions

Our findings suggest M. oryzae small RNAs differentially accumulate in vegetative and specialized-infection tissues and may play an active role in genome integrity and regulating growth and development.

MoSfl1 is important for virulence and heat tolerance in Magnaporthe oryzae

The formation of appressoria, specialized plant penetration structures of Magnaporthe oryzae, is regulated by the MST11-MST7-PMK1 MAP kinase cascade. One of its downstream transcription factor, MST12, is important for penetration and invasive growth but dispensable for appressorium formation. To identify additional downstream targets that are regulated by Pmk1, in this study we performed phosphorylation assays with a protein microarray composed of 573 M. oryzae transcription factor (TF) genes. Three of the TF genes phosphorylated by Pmk1 in vitro were further analyzed by coimmunoprecipitation assays. One of them, MoSFL1, was found to interact with Pmk1 in vivo. Like other Sfl1 orthologs, the MoSfl1 protein has the HSF-like domain. When expressed in yeast, MoSFL1 functionally complemented the flocculation defects of the sfl1 mutant. In M. oryzae, deletion of MoSFl1 resulted in a significant reduction in virulence on rice and barley seedlings. Consistent with this observation, the Mosfl1 mutant was defective in invasive growth in penetration assays with rice leaf sheaths. In comparison with that of vegetative hyphae, the expression level of MoSFL1 was increased in appressoria and infected rice leaves. The Mosfl1 mutant also had increased sensitivity to elevated temperatures. In CM cultures of the Mosfl1 and pmk1 mutants grown at 30°C, the production of aerial hyphae and melanization were reduced but their growth rate was not altered. When assayed by qRT-PCR, the transcription levels of the MoHSP30 and MoHSP98 genes were reduced 10- and 3-fold, respectively, in the Mosfl1 mutant. SFL1 orthologs are conserved in filamentous ascomycetes but none of them have been functionally characterized in non-Saccharomycetales fungi. MoSfl1 has one putative MAPK docking site and three putative MAPK phosphorylation sites. Therefore, it may be functionally related to Pmk1 in the regulation of invasive growth and stress responses in M. oryzae.

Direct comparison of stable isotope labeling by amino acids in cell culture and spectral counting for quantitative proteomics

Numerous experimental strategies exist for relative protein quantification, one of the primary objectives of mass spectrometry based proteomics analysis. These strategies mostly involve the incorporation of a stable isotope label via either metabolic incorporation in cell or tissue culture (¹⁵N/¹⁴N metabolic labeling, stable isotope labeling by amino acids in cell culture (SILAC)), chemical derivatization (ICAT, iTRAQ, TMT), or enzymatically catalyzed incorporation (¹⁸O labeling). Also, these techniques can be cost or time prohibitive or not amenable to the biological system of interest (i.e., metabolic labeling of clinical samples, most animals, or fungi). This is the case with the quantification of fungal proteomes, which often require auxotroph mutants to fully metabolically label. Alternatively, label-free strategies for protein quantification such as using integrated ion abundance and spectral counting have been demonstrated for quantification affording over 2 orders of magnitude of dynamic range which is comparable to metabolic labeling strategies. Direct comparisons of these quantitative techniques are largely lacking in the literature but are highly warranted in order to evaluate the capabilities, limitations, and analytical variability of available quantitative strategies. Here, we present the direct comparison of SILAC to label-free quantification by spectral counting of an identical set of data from the bottom-up proteomic analysis of human embryonic stem cells, which are readily able to be quantified using both strategies, finding that both strategies result in a similar number of protein identifications. We also discuss necessary constraints for accurate quantification using spectral counting and assess the potential of this label-free strategy as a viable alternative for quantitative proteomics.

Improving proteome coverage on a LTQ-Orbitrap using design of experiments

Design of experiments (DOE) was used to determine improved settings for a LTQ-Orbitrap XL to maximize proteome coverage of Saccharomyces cerevisiae. A total of nine instrument parameters were evaluated with the best values affording an increase of approximately 60% in proteome coverage. Utilizing JMP software, 2 DOE screening design tables were generated and used to specify parameter values for instrument methods. DOE 1, a fractional factorial design, required 32 methods fully resolving the investigation of six instrument parameters involving only half the time necessary for a full factorial design of the same resolution. It was advantageous to complete a full factorial design for the analysis of three additional instrument parameters. Measured with a maximum of 1% false discovery rate, protein groups, unique peptides, and spectral counts gauged instrument performance. Randomized triplicate nanoLC-LTQ-Orbitrap XL MS/MS analysis of the S. cerevisiae digest demonstrated that the following five parameters significantly influenced proteome coverage of the sample: (1) maximum ion trap ionization time; (2) monoisotopic precursor selection; (3) number of MS/MS events; (4) capillary temperature; and (5) tube lens voltage. Minimal influence on the proteome coverage was observed for the remaining four parameters (dynamic exclusion duration, resolving power, minimum count threshold to trigger a MS/MS event, and normalized collision energy). The DOE approach represents a time- and cost-effective method for empirically optimizing MS-based proteomics workflows including sample preparation, LC conditions, and multiple instrument platforms.

Increasing proteome coverage with offline RP HPLC coupled to online RP nanoLC-MS

Fractionation prior to mass spectrometry is an indispensable step in proteomics. In this paper we report the success of performing offline reversed phase high pressure liquid chromatography (HPLC) fractionation on a C18 2.0 mm×150 mm column at the peptide level with microliter per minute flow rates prior to online nano-flow reversed phase liquid chromatography mass spectrometry (nanoLC-MS) using the well-studied fungus Saccharomyces cerevisiae. A C18 75 μm×150 mm column was used online and the online elution gradients for each fraction were adjusted in order to obtain well resolved separation. Comparing this method directly to only performing nanoLC-MS we observed a 61.6% increase in the number of identified proteins. At a 1% false discovery rate 1028 proteins were identified using two dimensions of RPLC versus 636 proteins identified in a single nano-flow separation. The majority of proteins identified by one dimension of nano-LC were present in the proteins identified in our two dimensional strategy. Although increasing analysis time, this non-orthogonal and facile pre-fractionation method affords a more comprehensive examination of the proteome.

Genome-wide characterization of methylguanosine-capped and polyadenylated small RNAs in the rice blast fungus Magnaporthe oryzae

Small RNAs are well described in higher eukaryotes such as mammals and plants; however, knowledge in simple eukaryotes such as filamentous fungi is limited. In this study, we discovered and characterized methylguanosine-capped and polyadenylated small RNAs (CPA-sRNAs) by using differential RNA selection, full-length cDNA cloning and 454 transcriptome sequencing of the rice blast fungus Magnaporthe oryzae. This fungus causes blast, a devastating disease on rice, the principle food staple for over half the world’s population. CPA-sRNAs mapped primarily to the transcription initiation and termination sites of protein-coding genes and were positively correlated with gene expression, particularly for highly expressed genes including those encoding ribosomal proteins. Numerous CPA-sRNAs also mapped to rRNAs, tRNAs, snRNAs, transposable elements and intergenic regions. Many other 454 sequence reads could not be mapped to the genome; however, inspection revealed evidence for non-template additions and chimeric sequences. CPA-sRNAs were independently confirmed using a high affinity variant of eIF-4E to capture 5′-methylguanosine-capped RNA followed by 3′-RACE sequencing. These results expand the repertoire of small RNAs in filamentous fungi.

Combining ChIP-chip and expression profiling to model the MoCRZ1 mediated circuit for Ca/calcineurin signaling in the rice blast fungus

Significant progress has been made in defining the central signaling networks in many organisms, but collectively we know little about the downstream targets of these networks and the genes they regulate. To reconstruct the regulatory circuit of calcineurin signal transduction via MoCRZ1, a Magnaporthe oryzae C2H2 transcription factor activated by calcineurin dephosphorylation, we used a combined approach of chromatin immunoprecipitation - chip (ChIP-chip), coupled with microarray expression studies. One hundred forty genes were identified as being both a direct target of MoCRZ1 and having expression concurrently differentially regulated in a calcium/calcineurin/MoCRZ1 dependent manner. Highly represented were genes involved in calcium signaling, small molecule transport, ion homeostasis, cell wall synthesis/maintenance, and fungal virulence. Of particular note, genes involved in vesicle mediated secretion necessary for establishing host associations, were also found. MoCRZ1 itself was a target, suggesting a previously unreported autoregulation control point. The data also implicated a previously unreported feedback regulation mechanism of calcineurin activity. We propose that calcium/calcineurin regulated signal transduction circuits controlling development and pathogenicity manifest through multiple layers of regulation. We present results from the ChIP-chip and expression analysis along with a refined model of calcium/calcineurin signaling in this important plant pathogen.

All organisms have the innate ability to perceive their environment and respond to it, largely through controlling gene expression. Tailored specificity of a response is primarily achieved through signal cascades involving unique receptors, downstream transcription factors (proteins that bind to DNA to regulate gene expression), and the genes these transcription factors regulate. For fungal plant pathogens, signal transduction cascades are involved in perception of hosts, transgression of physical barriers, suppression or elicitation of host defenses, in vivo nutrient acquisition, and completion of their life cycle. We know that the Ca²⁺/calcineurin signaling pathway is a central conduit regulating these aspects of the life cycle for fungal pathogens of plants and animals. In this study, we used advanced ChIP-chip and microarray gene expression technologies to identify the genes that the Ca²⁺/calcineurin responsive transcription factor MoCRZ1 directly binds to and regulates the expression of. Our findings show conservations and divergence in this pathway within the fungal kingdom. It also identifies points of control in the pathway that were previously unidentified. Most importantly, this study implicates this pathway in the establishment of host associations and virulence for the causal agent of rice blast disease, Magnaporthe oryzae, the most important disease of rice worldwide.

Aspergillus flavus genomics as a tool for studying the mechanism of aflatoxin formation

Aspergillus flavus is a weak pathogen that infects plants, animals and humans. When it infects agricultural crops, however, it produces one of the most potent carcinogens known (aflatoxins). To devise strategies to control aflatoxin contamination of pre-harvest agricultural crops and post-harvest grains during storage, we launched the A. flavus genomics program. The major objective of this program is the identification of genes involved in aflatoxin biosynthesis and regulation, as well as in pathogenicity, to gain a better understanding of the mechanism of aflatoxin formation. The sequencing of A. flavus whole genome has been completed. Initial annotation of the sequence revealed that there are about 13,071 genes in the A. flavus genome. Genes which potentially encode for enzymes involved in secondary metabolite production in the A. flavus genome have been identified. Preliminary comparative genome analysis of A. flavus with A. oryzae is summarized here.

Infection strategies of filamentous microbes described with the Gene Ontology

Filamentous microbes that form highly developed symbiotic associations (ranging from pathogenesis to mutualism) with their hosts include fungi, oomycetes and actinomycete bacteria. These organisms share many common features in growth, development and infection and have evolved similar strategies for neutralizing host defense responses to establish symbioses. Recent advances in sequencing technologies have led to a remarkable increase in the number of sequenced genomes of filamentous organisms. Analysis of the available genomes has provided useful information about genes that might be important for host infection and colonization. However, because many functional similarities among these organisms have arisen by convergent evolution, sequence-based genomic comparisons will miss many genes that are functionally analogous. In the absence of sequence similarity, annotating genes with standardized terms from the Gene Ontology (GO) can facilitate functional comparisons. Here, we review common strategies employed by filamentous organisms during colonization of their hosts, with reference to GO terms that best describe the processes involved.

The influence of matrix type, diurnal rhythm and sample collection and processing on the measurement of plasma beta-amyloid isoforms using the INNO-BIA plasma Abeta forms multiplex assay

OBJECTIVE

The aim of the study was to determine the extent to which plasma matrix types, diurnal rhythm and sample collection and processing procedures contribute to overall variability of measurements with the INNO-BIA plasma Abeta forms assay.

METHODS

Plasma samples from healthy volunteers were collected at BARC-CRI. Analyte concentrations from various plasma matrix types (EDTA, heparin, fluoride) were compared to serum after collection of blood in commercial plastic and glass tubes. Sample processing variables including time and temperature before and after centrifugation, centrifugal force and plasma dilution factor were also investigated. Diurnal variability in plasma Abeta isoforms was determined in 29 healthy volunteers by analysis of EDTA plasma specimens serially collected over 24 hours and stored frozen following oral administration of a placebo treatment. All plasma samples from a given individual and experiment were analyzed in a single analytical run.

RESULTS

Highest Abeta levels were obtained using EDTA-plasma samples (in contrast to serum, heparin, citrate, or fluoride). Addition of aprotinin to EDTA plasma had no effect on Abeta peptide recovery. The elapsed time and temperature exposure, before and after sample processing affects the recovery of Abeta isoforms. Analyte recovery was not significantly affected by the presence of platelets in plasma samples. At the subject level, analysis of serially collected EDTA plasma specimens from healthy volunteers revealed no evidence of diurnal variation in any of the Abeta isoforms investigated and results from samples collected on a monthly basis showed only very limited intra-individual variation.

CONCLUSIONS

Optimal recovery of Abeta peptides was obtained from blood drawn into EDTA tubes and processed within 4 h. Plasma that was refrigerated after separation and analysed within 4 h gave comparable results to samples immediately processed and frozen at -70 degrees C.

Common processes in pathogenesis by fungal and oomycete plant pathogens, described with Gene Ontology terms

Plant diseases caused by fungi and oomycetes result in significant economic losses every year. Although phylogenetically distant, the infection processes by these organisms share many common features. These include dispersal of an infectious particle, host adhesion, recognition, penetration, invasive growth, and lesion development. Previously, many of these common processes did not have corresponding Gene Ontology (GO) terms. For example, no GO terms existed to describe processes related to the appressorium, an important structure for infection by many fungi and oomycetes. In this mini-review, we identify common features of the pathogenic processes of fungi and oomycetes and create a pathogenesis model using 256 newly developed and 38 extant GO terms, with an emphasis on the appressorium and signal transduction. This set of standardized GO terms provides a solid base to further compare and contrast the molecular underpinnings of fungal and oomycete pathogenesis.

Gene Ontology annotation of the rice blast fungus, Magnaporthe oryzae

Background

Magnaporthe oryzae, the causal agent of blast disease of rice, is the most destructive disease of rice worldwide. The genome of this fungal pathogen has been sequenced and an automated annotation has recently been updated to Version 6 . However, a comprehensive manual curation remains to be performed. Gene Ontology (GO) annotation is a valuable means of assigning functional information using standardized vocabulary. We report an overview of the GO annotation for Version 5 of M. oryzae genome assembly.

Methods

A similarity-based (i.e., computational) GO annotation with manual review was conducted, which was then integrated with a literature-based GO annotation with computational assistance. For similarity-based GO annotation a stringent reciprocal best hits method was used to identify similarity between predicted proteins of M. oryzae and GO proteins from multiple organisms with published associations to GO terms. Significant alignment pairs were manually reviewed. Functional assignments were further cross-validated with manually reviewed data, conserved domains, or data determined by wet lab experiments. Additionally, biological appropriateness of the functional assignments was manually checked.

Results

In total, 6,286 proteins received GO term assignment via the homology-based annotation, including 2,870 hypothetical proteins. Literature-based experimental evidence, such as microarray, MPSS, T-DNA insertion mutation, or gene knockout mutation, resulted in 2,810 proteins being annotated with GO terms. Of these, 1,673 proteins were annotated with new terms developed for Plant-Associated Microbe Gene Ontology (PAMGO). In addition, 67 experiment-determined secreted proteins were annotated with PAMGO terms. Integration of the two data sets resulted in 7,412 proteins (57%) being annotated with 1,957 distinct and specific GO terms. Unannotated proteins were assigned to the 3 root terms. The Version 5 GO annotation is publically queryable via the GO site . Additionally, the genome of M. oryzae is constantly being refined and updated as new information is incorporated. For the latest GO annotation of Version 6 genome, please visit our website . The preliminary GO annotation of Version 6 genome is placed at a local MySql database that is publically queryable via a user-friendly interface Adhoc Query System.

Conclusion

Our analysis provides comprehensive and robust GO annotations of the M. oryzae genome assemblies that will be solid foundations for further functional interrogation of M. oryzae.

GT-Miner: a graph-theoretic data miner, viewer, and model processor

Inexpensive computational power combined with high-throughput experimental platforms has created a wealth of biological information requiring analytical tools and techniques for interpretation. Graph-theoretic concepts and tools have provided an important foundation for information visualization, integration, and analysis of datasets, but they have often been relegated to background analysis tasks. GT-Miner is designed for visual data analysis and mining operations, interacts with other software, including databases, and works with diverse data types. It facilitates a discovery-oriented approach to data mining wherein exploration of alterations of the data and variations of the visualization is encouraged. The user is presented with a basic iterative process, consisting of loading, visualizing, transforming, and then storing the resultant information. Complex analyses are built-up through repeated iterations and user interactions. The iterative process is optimized by automatic layout following transformations and by maintaining a current selection set of interest for elements modified by the transformations. Multiple visualizations are supported including hierarchical, spring, and force-directed self-organizing layouts. Graphs can be transformed with an extensible set of algorithms or manually with an integral visual editor. GT-Miner is intended to allow easier access to visual data mining for the non-expert.

AVAILABILITY

The GT-Miner program and supplemental materials, including example uses and a user guide, are freely available from http://www.cifr.ncsu.edu/bioinformatics/downloads/

Transcriptome analysis reveals new insight into appressorium formation and function in the rice blast fungus Magnaporthe oryzae

BACKGROUND

Rice blast disease is caused by the filamentous Ascomycetous fungus Magnaporthe oryzae and results in significant annual rice yield losses worldwide. Infection by this and many other fungal plant pathogens requires the development of a specialized infection cell called an appressorium. The molecular processes regulating appressorium formation are incompletely understood.

RESULTS

We analyzed genome-wide gene expression changes during spore germination and appressorium formation on a hydrophobic surface compared to induction by cAMP. During spore germination, 2,154 (approximately 21%) genes showed differential expression, with the majority being up-regulated. During appressorium formation, 357 genes were differentially expressed in response to both stimuli. These genes, which we refer to as appressorium consensus genes, were functionally grouped into Gene Ontology categories. Overall, we found a significant decrease in expression of genes involved in protein synthesis. Conversely, expression of genes associated with protein and amino acid degradation, lipid metabolism, secondary metabolism and cellular transportation exhibited a dramatic increase. We functionally characterized several differentially regulated genes, including a subtilisin protease (SPM1) and a NAD specific glutamate dehydrogenase (Mgd1), by targeted gene disruption. These studies revealed hitherto unknown findings that protein degradation and amino acid metabolism are essential for appressorium formation and subsequent infection.

CONCLUSION

We present the first comprehensive genome-wide transcript profile study and functional analysis of infection structure formation by a fungal plant pathogen. Our data provide novel insight into the underlying molecular mechanisms that will directly benefit efforts to identify fungal pathogenicity factors and aid the development of new disease management strategies.

Expression of a harpin-encoding gene in rice confers durable nonspecific resistance to Magnaporthe grisea

Engineering durable nonspecific resistance to phytopathogens is one of the ultimate goals of plant breeding. However, most attempts to reach this goal fail as a result of rapid changes in pathogen populations and the sheer diversity of pathogen infection mechanisms. In this study, we show that the expression of a harpin-encoding gene (hrf1), derived from Xanthomonas oryzae pv. oryzae, confers nonspecific resistance in rice to the blast fungus Magnaporthe grisea. Transgenic plants and their T1-T7 progenies were highly resistant to all major M. grisea races in rice-growing areas along the Yangtze River, China. The expression of defence-related genes was activated in resistant transgenic plants, and the formation of melanized appressoria, which is essential for foliar infection, was inhibited on plant leaves. These results suggest that harpins may offer new opportunities for generating broad-spectrum disease resistance in other crops.

A systematic analysis of T-DNA insertion events in Magnaporthe oryzae

We describe here the analysis of random T-DNA insertions that were generated as part of a large-scale insertional mutagenesis project for Magnaporthe oryzae. Chromosomal regions flanking T-DNA insertions were rescued by inverse PCR, sequenced and used to search the M. oryzae genome assembly. Among the 175 insertions for which at least one flank was rescued, 137 had integrated in single-copy regions of the genome, 17 were in repeated sequences, one had no match to the genome, and the remainder were unassigned due to illegitimate T-DNA integration events. These included in order of abundance: head-to-tail tandem insertions, right border excision failures, left border excision failures and insertion of one T-DNA into another. The left borders of the T-DNA were frequently truncated and inserted in sequences with micro-homology to the left terminus. By contrast the right borders were less prone to degradation and appeared to have been integrated in a homology-independent manner. Gross genome rearrangements rarely occurred when the T-DNAs integrated in single-copy regions, although most insertions did cause small deletions at the target site. Significant insertion bias was detected, with promoters receiving two times more T-DNA hits than expected, and open reading frames receiving three times fewer. In addition, we found that the distribution of T-DNA inserts among the M. oryzae chromosomes was not random. The implications of these findings with regard to saturation mutagenesis of the M. oryzae genome are discussed.

Development of a high throughput transformation system for insertional mutagenesis in Magnaporthe oryzae

Towards the goal of disrupting all genes in the genome of Magnaporthe oryzae and identifying their function, a collection of >55,000 random insertion lines of M. oryzae strain 70-15 were generated. All strains were screened to identify genes involved in growth rate, conidiation, pigmentation, auxotrophy, and pathogenicity. Here, we provide a description of the high throughput transformation and analysis pipeline used to create our library. Transformed lines were generated either by CaCl(2)/PEG treatment of protoplasts with DNA or by Agrobacterium tumefaciens-mediated transformation (ATMT). We describe the optimization of both approaches and compare their efficiency. ATMT was found to be a more reproducible method, resulting in predominantly single copy insertions, and its efficiency was high with up to 0.3% of conidia being transformed. The phenotypic data is accessible via a public database called MGOS and all strains are publicly available. This represents the most comprehensive insertional mutagenesis analysis of a fungal pathogen.

Transcript profiling of a conifer pathosystem: response of Pinus sylvestris root tissues to pathogen (Heterobasidion annosum) invasion

The mechanisms underlying defence reactions to a pathogen attack, though well studied in crop plants, are poorly understood in conifers. To analyze changes in gene transcript abundance in Pinus sylvestris L. root tissues infected by Heterobasidion annosum (Fr.) Bref. s.l., a cDNA microarray containing 2109 ESTs from P. taeda L. was used. Mixed model statistical analysis identified 179 expressed sequence tags differentially expressed at 1, 5 or 15 days post inoculation. In general, the total number of genes differentially expressed during the infection increased over time. The most abundant group of genes up-regulated upon infection coded for enzymes involved in metabolism (phenylpropanoid pathway) and defence-related proteins with antimicrobial properties. A class III peroxidase responsible for lignin biosynthesis and cell wall thickening had increased transcript abundance at all measurement times. Real-time RT-PCR verified the microarray results with high reproducibility. The similarity of the expression profiling pattern observed in this pathosystem to those documented in crop pathology suggests that angiosperms and gymnosperms use similar genetic programs in responding to invasive growth by microbial pathogens.