Proteomic identification of extracellular proteins regulated by the Gna1 Gα subunit in Stagonospora nodorum

Comparative sub-cellular proteome analyses reveals metabolic differentiation and production of effector-like molecules in the dieback phytopathogen Phytophthora cinnamomi

Phytopathogenic oomycetes pose a significant threat to global biodiversity and food security. The proteomes of these oomycetes likely contain important factors that contribute to their pathogenic success, making their discovery crucial for elucidating pathogenicity. Phytophthora cinnamomi is a root pathogen that causes dieback in a wide variety of crops and native vegetation world-wide. Virulence proteins produced by P. cinnamomi are not well defined and a large-scale approach to understand the biochemistry of this pathogen has not been documented. Soluble mycelial, zoospore and secreted proteomes were obtained and label-free quantitative proteomics was used to compare the composition of the three sub-proteomes. A total of 4635 proteins were identified, validating 17.7% of the predicted gene set. The mycelia were abundant in transporters for nutrient acquisition, metabolism and cellular proliferation. The zoospores had less metabolic related ontologies but were abundant in energy generating, motility and signalling associated proteins. Virulence-associated proteins were identified in the secretome such as candidate effector and effector-like proteins, which interfere with the host immune system. These include hydrolases, cell wall degrading enzymes, putative necrosis-inducing proteins and elicitins. The secretome elicited a hypersensitive response on the roots of a model host and thus suggests evidence of effector activity. SIGNIFICANCE: Phytophthora cinnamomi is a phytopathogenic oomycete that causes dieback disease in native vegetation and several horticultural crops such as avocado, pineapple and macadamia. Whilst this pathogen has significance world-wide, its pathogenicity and virulence have not been described in depth. We carried out comparative label-free proteomics of the mycelia, zoospores and secretome of P. cinnamomi. This study highlights the differential metabolism and cellular processes between the sub-proteomes. Proteins associated with metabolism, nutrient transport and cellular proliferation were over represented in the mycelia. The zoospores have a specialised proteome showing increased energy generation geared towards motility. Candidate effectors and effector-like secreted proteins were also identified, which can be exploited for genetic resistance. This demonstrates a better understanding of the biology and pathogenicity of P. cinnamomi infection that can subsequently be used to develop effective methods of disease management.

The G-protein alpha subunit CgGa1 mediates growth, sporulation, penetration and pathogenicity in Colletotrichum gloeosporioides

Colletotrichum gloeosporioides is the main pathogen causing rubber anthracnose, which brings huge economic loss to the natural rubber industry. Heterotrimeric G proteins play a vital role in signal transduction in filamentous fungi, and G alpha subunits are the major component of G proteins. In this study, we characterize a group I Gα subunit CgGa1 in C. gloeosporioides as a homolog of MagB in Pyricularia oryzae. CgGa1 encodes a 353-amino acid protein and has a G_alpha domain. Deletion of CgGa1 results in reduced vegetative growth and conidia yield, and the mutant cannot produce a fruiting body. The CgGa1 deletion mutant also exhibits decreased conidial germination and appressorium formation significantly. Moreover, the mutant has an obvious deficiency in penetration and loses its virulence completely. Transcriptome analysis showed that CgGa1 could affect the expression of many genes related to carbohydrate metabolism, amino acid metabolism and signal transduction, etc. In conclusion, CgGa1 regulates growth, asexual and sexual sporulation, appressorium formation, penetration and pathogenicity of C. gloeosporioides.

Gene Validation and Remodelling Using Proteogenomics of Phytophthora cinnamomi, the Causal Agent of Dieback

Phytophthora cinnamomi is a pathogenic oomycete that causes plant dieback disease across a range of natural ecosystems and in many agriculturally important crops on a global scale. An annotated draft genome sequence is publicly available (JGI Mycocosm) and suggests 26,131 gene models. In this study, soluble mycelial, extracellular (secretome), and zoospore proteins of P. cinnamomi were exploited to refine the genome by correcting gene annotations and discovering novel genes. By implementing the diverse set of sub-proteomes into a generated proteogenomics pipeline, we were able to improve the P. cinnamomi genome annotation. Liquid chromatography mass spectrometry was used to obtain high confidence peptides with spectral matching to both the annotated genome and a generated 6-frame translation. Two thousand seven hundred sixty-four annotations from the draft genome were confirmed by spectral matching. Using a proteogenomic pipeline, mass spectra were used to edit the P. cinnamomi genome and allowed identification of 23 new gene models and 60 edited gene features using high confidence peptides obtained by mass spectrometry, suggesting a rate of incorrect annotations of 3% of the detectable proteome. The novel features were further validated by total peptide support, alongside functional analysis including the use of Gene Ontology and functional domain identification. We demonstrated the use of spectral data in combination with our proteogenomics pipeline can be used to improve the genome annotation of important plant diseases and identify missed genes. This study presents the first use of spectral data to edit and manually annotate an oomycete pathogen.

Comparative quantitative proteomics of osmotic signal transduction mutants in Botrytis cinerea explain mutant phenotypes and highlight interaction with cAMP and Ca2+ signalling pathways

Signal transduction (ST) is essential for rapid adaptive responses to changing environmental conditions. It acts through rapid post-translational modifications of signalling proteins and downstream effectors that regulate the activity and/or subcellular localisation of target proteins, or the expression of downstream genes. We have performed a quantitative, comparative proteomics study of ST mutants in the phytopathogenic fungus Botrytis cinerea during axenic growth under non-stressed conditions to decipher the roles of two kinases of the hyper-osmolarity pathway in B. cinerea physiology. We studied the mutants of the sensor histidine kinase Bos1 and of the MAP kinase Sak1. Label-free shotgun proteomics detected 2425 proteins, 628 differentially abundant between mutants and wild-type, 270 common to both mutants, indicating independent and shared regulatory functions for both kinases. Gene ontology analysis showed significant changes in functional categories that may explain in vitro growth and virulence defects of both mutants (secondary metabolism enzymes, lytic enzymes, proteins linked to osmotic, oxidative and cell wall stress). The proteome data also highlight a new link between Sak1 MAPK, cAMP and Ca²⁺ signalling. This study reveals the potential of proteomic analyses of signal transduction mutants to decipher their biological functions. TEXT-VULGARISATION: The fungus Botrytis cinerea is responsible for grey mold disease of hundreds of plant species. During infection, the fungus has to face important changes of its environment. Adaptation to these changing environmental conditions involves proteins of such called signal transduction pathways that regulate the production, activity or localisation of cellular components, mainly proteins. While the components of such signal transduction pathways are well known, their role globally understood, the precise impact on protein production remains unknown. In this study we have analysed and compared the global protein content of two Botrytis cinerea signal transduction mutants - both avirulent - to the pathogenic parental strain. The data of 628 differential proteins between mutants and wild-type, showed significant changes in proteins related to plant infection (secondary metabolism enzymes, lytic enzymes, proteins linked to osmotic, oxidative and cell wall stress) that may explain the virulence defects of both mutants. Moreover, we observed intracellular accumulation of secreted proteins in one of the mutants suggesting a potential secretion defect.

Analysis of Reproducibility of Proteome Coverage and Quantitation Using Isobaric Mass Tags (iTRAQ and TMT)

This study aimed to compare the depth and reproducibility of total proteome and differentially expressed protein coverage in technical duplicates and triplicates using iTRAQ 4-plex, iTRAQ 8-plex, and TMT 6-plex reagents. The analysis was undertaken because comprehensive comparisons of isobaric mass tag reproducibility have not been widely reported in the literature. The highest number of proteins was identified with 4-plex, followed by 8-plex and then 6-plex reagents. Quantitative analyses revealed that more differentially expressed proteins were identified with 4-plex reagents than 8-plex reagents and 6-plex reagents. Replicate reproducibility was determined to be ≥69% for technical duplicates and ≥57% for technical triplicates. The results indicate that running an 8-plex or 6-plex experiment instead of a 4-plex experiment resulted in 26 or 39% fewer protein identifications, respectively. When 4-plex spectra were searched with three software tools-ProteinPilot, Mascot, and Proteome Discoverer-the highest number of protein identifications were obtained with Mascot. The analysis of negative controls demonstrated the importance of running experiments as replicates. Overall, this study demonstrates the advantages of using iTRAQ 4-plex reagents over iTRAQ 8-plex and TMT 6-plex reagents, provides estimates of technical duplicate and triplicate reproducibility, and emphasizes the value of running replicate samples.

Genome and Transcriptome Sequences Reveal the Specific Parasitism of the Nematophagous Purpureocillium lilacinum 36-1

Purpureocillium lilacinum is a promising nematophagous ascomycete able to adapt diverse environments and it is also an opportunistic fungus that infects humans. A microbial inoculant of P. lilacinum has been registered to control plant parasitic nematodes. However, the molecular mechanism of the toxicological processes is still unclear because of the relatively few reports on the subject. In this study, using Illumina paired-end sequencing, the draft genome sequence and the transcriptome of P. lilacinum strain 36-1 infecting nematode-eggs were determined. Whole genome alignment indicated that P. lilacinum 36-1 possessed a more dynamic genome in comparison with P. lilacinum India strain. Moreover, a phylogenetic analysis showed that the P. lilacinum 36-1 had a closer relation to entomophagous fungi. The protein-coding genes in P. lilacinum 36-1 occurred much more frequently than they did in other fungi, which was a result of the depletion of repeat-induced point mutations (RIP). Comparative genome and transcriptome analyses revealed the genes that were involved in pathogenicity, particularly in the recognition, adhesion of nematode-eggs, downstream signal transduction pathways and hydrolase genes. By contrast, certain numbers of cellulose and xylan degradation genes and a lack of polysaccharide lyase genes showed the potential of P. lilacinum 36-1 as an endophyte. Notably, the expression of appressorium-formation and antioxidants-related genes exhibited similar infection patterns in P. lilacinum strain 36-1 to those of the model entomophagous fungi Metarhizium spp. These results uncovered the specific parasitism of P. lilacinum and presented the genes responsible for the infection of nematode-eggs.

Dissecting the role of histidine kinase and HOG1 mitogen-activated protein kinase signalling in stress tolerance and pathogenicity of Parastagonospora nodorum on wheat

The HOG1 mitogen-activated protein kinase (MAPK) pathway is activated through two-component histidine kinase (HK) signalling. This pathway was first characterized in the budding yeast Saccharomyces cerevisiae as a regulator of osmotolerance. The fungus Parastagonospora nodorum is the causal agent of septoria nodorum blotch of wheat. This pathogen uses host-specific effectors in tandem with general pathogenicity mechanisms to carry out its infection process. Genes showing strong sequence homology to S. cerevisiae HOG1 signalling pathway genes have been identified in the genome of P. nodorum. In this study, we examined the role of the pathway in the virulence of P. nodorum on wheat by disrupting putative pathway component genes: HOG1 (SNOG_13296) MAPK and NIK1 (SNOG_11631) hybrid HK. Mutants deleted in NIK1 and HOG1 were insensitive to dicarboximide and phenylpyrrole fungicides, but not a fungicide that targets ergosterol biosynthesis. Furthermore, both Δnik1 and Δhog1 mutants showed increased sensitivity to hyperosmotic stress. However, HOG1, but not NIK1, is required for tolerance to elevated temperatures. HOG1 deletion conferred increased tolerance to 6-methoxy-2-benzoxazolinone, a cereal phytoalexin. This suggests that the HOG1 signalling pathway is not exclusively associated with NIK1. Both Δnik1 and Δhog1 mutants retained the ability to infect and cause necrotic lesions on wheat. However, we observed that the Δhog1 mutation resulted in reduced production of pycnidia, asexual fruiting bodies that facilitate spore dispersal during late infection. Our study demonstrated the overlapping and distinct roles of a HOG1 MAPK and two-component HK signalling in P. nodorum growth and pathogenicity.

Small GTPase Rac1 and its interaction partner Cla4 regulate polarized growth and pathogenicity in Verticillium dahliae

Rac1 is a small GTPase coordinating diverse cellular functions such as cell polarity, vesicular trafficking, the cell cycle and transcriptional dynamics in many organisms. In this study, we investigate the biological functions of VdRac1, a Rac1 homolog in the soil-borne, wilt-causing fungus Verticillium dahliae. The VdRac1 gene was deleted in a V. dahliae virulence strain Vd8 isolated from a local cotton cultivar. ΔVdrac1 mutants display drastic reduction in colony expansion and form compact, convoluted colonies, show hyper-branching, loss of polarity and ability to penetrate, leading to severely reduced virulence. The p21-activated kinase Cla4 (named as VdCla4 in V. dahliae) null mutants ΔVdcla4 share identical phenotypes with ΔVdrac1. Yeast two-hybrid studies prove that VdCla4 is an effector of VdRac1. Localizations of actin and reactive oxygen species (ROS) in ΔVdrac1 and ΔVdcla4 compared with the corresponding wild-type strain reveal that VdRac1 and VdCla4 play a primary role in polarized hyphal growth via organization of ROS and play only a minor role in the organization of actin. The Vdrac1 and Vdcla4 null mutants are defective in conidiation and trace elements can partially compensate for the defect. Our data demonstrate that VdRac1 regulates polarized growth and pathogenicity by interacting with its effector VdCla4 in V. dahliae.

A comparative analysis of the heterotrimeric G-protein Gα, Gβ and Gγ subunits in the wheat pathogen Stagonospora nodorum

Background

It has been well established that the Gα subunit of the heterotrimeric G-protein in the wheat pathogen Stagonospora nodorum is required for a variety of phenotypes including pathogenicity, melanisation and asexual differentiation. The roles though of the Gγ and Gβ subunits though were unclear. The objective of this study was to identify and understand the role of these subunits and assess their requirement for pathogenicity and development.

Results

G-protein Gγ and Gβ subunits, named Gga1 and Gba1 respectively, were identified in the Stagonospora nodorum genome by comparative analysis with known fungal orthologues. A reverse genetics technique was used to study the role of these and revealed that the mutant strains displayed altered in vitro growth including a differential response to a variety of exogenous carbon sources. Pathogenicity assays showed that Stagonospora nodorum strains lacking Gba1 were essentially non-pathogenic whilst Gga1-impaired strains displayed significantly slower growth in planta. Subsequent sporulation assays showed that like the previously described Gα subunit mutants, both Gba1 and Gga1 were required for asexual sporulation with neither mutant strain being able to differentiate either pycnidia nor pycnidiospores under normal growth conditions. Continued incubation at 4°C was found to complement the mutation in each of the G-protein subunits with nearly wild-type levels of pycnidia recovered.

Conclusion

This study provides further evidence on the significance of cAMP-dependent signal transduction for many aspects of fungal development and pathogenicity. The observation that cold temperatures can complement the G-protein sporulation defect now provides an ideal tool by which asexual differentiation can now be dissected.

Proteomic techniques for plant-fungal interactions

Proteomics is a key technique that is helping elucidate many complex biological processes. The analysis of plant-pathogen interactions using proteomics is complicated by the presence of the proteomes of two species, but is benefiting from the developing maturity and power of these techniques. More and more pathogen genomes are being sequenced, so fungal proteomics is reaching its full potential and remains the chosen technology to unravel the molecular pathways of pathogenicity and resistance. In this chapter, we suggest proteomic strategies that have proved successful on various plant-interacting fungal species. Several protein extraction methods are described. For adequate quantitative analyses of protein abundances, we recommend either separation using two-dimensional gel electrophoresis or labelling with isobaric tags followed by two-dimensional HPLC separation. Proteins of interest are then identified using mass spectrometry. Identified proteins can assist in refining genome annotations, otherwise known as proteogenomics.

The G protein β subunit controls virulence and multiple growth- and development-related traits in Verticillium dahliae

To gain insight into the role of G protein-mediated signaling in virulence and development of the soilborne, wilt causing fungus Verticillium dahliae, the G protein β subunit gene (named as VGB) was disrupted in tomato race 1 strain of V. dahliae. A resulting mutant strain, 70ΔGb15, displayed drastic reduction in virulence, increased microsclerotia formation and conidiation, and decreased ethylene production compared to the corresponding wild type (wt) strain 70wt-r1. Moreover, 70ΔGb15 exhibited an elongated rather than radial growth pattern on agar media. A transformant of 70ΔGb15 (named as 70ΔGbPKAC1) that carries an extra copy of VdPKAC1, a V. dahliae gene encoding the catalytic subunit of the cAMP-dependent protein kinase A, exhibited wt growth pattern and conidiation, was unable to form microsclerotia, produced high amounts of ethylene, and exhibited virulence between that of 70ΔGb15 and 70wt-r1 on tomato plants. Phenotypical changes observed in 70ΔGb15 and 70ΔGbPKAC1 correlated with transcriptional changes in several genes involved in signaling (MAP kinase VMK1) and development (hydrophobin VDH1 and ACC synthase ACS1) of V. dahliae. Results from the present work suggest a linkage between VGB and VdPKAC1 signaling pathways in regulating virulence, hormone production and development in V. dahliae.

Stagonospora nodorum: from pathology to genomics and host resistance

Stagonospora nodorum is a major necrotrophic pathogen of wheat that causes the diseases S. nodorum leaf and glume blotch. A series of tools and resources, including functional genomics, a genome sequence, proteomics and metabolomics, host-mapping populations, and a worldwide collection of isolates, have enabled the dissection of pathogenicity mechanisms. Metabolic and signaling genes required for pathogenicity have been defined. Interaction with the host is dominated by interplay of fungal effectors that induce necrosis on wheat lines carrying specific sensitivity loci. As such, the pathogen has emerged as a model for the Pleosporales group of pathogens.

A novel mode of chromosomal evolution peculiar to filamentous Ascomycete fungi

Background

Gene loss, inversions, translocations, and other chromosomal rearrangements vary among species, resulting in different rates of structural genome evolution. Major chromosomal rearrangements are rare in most eukaryotes, giving large regions with the same genes in the same order and orientation across species. These regions of macrosynteny have been very useful for locating homologous genes in different species and to guide the assembly of genome sequences. Previous analyses in the fungi have indicated that macrosynteny is rare; instead, comparisons across species show no synteny or only microsyntenic regions encompassing usually five or fewer genes. To test the hypothesis that chromosomal evolution is different in the fungi compared to other eukaryotes, synteny was compared between species of the major fungal taxa.

Results

These analyses identified a novel form of evolution in which genes are conserved within homologous chromosomes, but with randomized orders and orientations. This mode of evolution is designated mesosynteny, to differentiate it from micro- and macrosynteny seen in other organisms. Mesosynteny is an alternative evolutionary pathway very different from macrosyntenic conservation. Surprisingly, mesosynteny was not found in all fungal groups. Instead, mesosynteny appears to be restricted to filamentous Ascomycetes and was most striking between species in the Dothideomycetes.

Conclusions

The existence of mesosynteny between relatively distantly related Ascomycetes could be explained by a high frequency of chromosomal inversions, but translocations must be extremely rare. The mechanism for this phenomenon is not known, but presumably involves generation of frequent inversions during meiosis.

Hunting down fungal secretomes using liquid-phase IEF prior to high resolution 2-DE

The secreted proteins (secretome) of fungi play a key role in interactions of pathogenic and symbiotic fungi with plants. Using the plant pathogenic fungus Leptosphaeria maculans and symbiont Laccaria bicolor grown in culture, we have established a proteomic protocol for extraction, concentration and resolution of the fungal secretome. As no proteomic data were available on mycelium tissues from both L. maculans and L. bicolor, mycelial proteins were studied; they also helped verifying the purity of secretome samples. The quality of protein extracts was initially assessed by both 1-DE and 2-DE using first a broad pH range for IEF, and then narrower acidic and basic pH ranges, prior to 2-DE. Compared with the previously published protocols for which only dozens of 2-D spots were recovered from fungal secretome samples, up to approximately 2000 2-D spots were resolved by our method. MS identification of proteins along several pH gradients confirmed this high resolution, as well as the presence of major secretome markers such as endopolygalacturonases, beta-glucanosyltransferases, pectate lyases and endoglucanases. Shotgun proteomic experiments evidenced the enrichment of secreted protein within the liquid medium. This is the first description of the proteome of L. maculans and L. bicolor, and the first application of liquid-phase IEF to any fungal extracts.

Deep proteogenomics; high throughput gene validation by multidimensional liquid chromatography and mass spectrometry of proteins from the fungal wheat pathogen Stagonospora nodorum

Background

Stagonospora nodorum, a fungal ascomycete in the class dothideomycetes, is a damaging pathogen of wheat. It is a model for necrotrophic fungi that cause necrotic symptoms via the interaction of multiple effector proteins with cultivar-specific receptors. A draft genome sequence and annotation was published in 2007. A second-pass gene prediction using a training set of 795 fully EST-supported genes predicted a total of 10762 version 2 nuclear-encoded genes, with an additional 5354 less reliable version 1 genes also retained.

Results

In this study, we subjected soluble mycelial proteins to proteolysis followed by 2D LC MALDI-MS/MS. Comparison of the detected peptides with the gene models validated 2134 genes. 62% of these genes (1324) were not supported by prior EST evidence. Of the 2134 validated genes, all but 188 were version 2 annotations. Statistical analysis of the validated gene models revealed a preponderance of cytoplasmic and nuclear localised proteins, and proteins with intracellular-associated GO terms. These statistical associations are consistent with the source of the peptides used in the study. Comparison with a 6-frame translation of the S. nodorum genome assembly confirmed 905 existing gene annotations (including 119 not previously confirmed) and provided evidence supporting 144 genes with coding exon frameshift modifications, 604 genes with extensions of coding exons into annotated introns or untranslated regions (UTRs), 3 new gene annotations which were supported by tblastn to NR, and 44 potential new genes residing within un-assembled regions of the genome.

Conclusion

We conclude that 2D LC MALDI-MS/MS is a powerful, rapid and economical tool to aid in the annotation of fungal genomic assemblies.

Assessing the impact of transcriptomics, proteomics and metabolomics on fungal phytopathology

SUMMARY Peer-reviewed literature is today littered with exciting new tools and techniques that are being used in all areas of biology and medicine. Transcriptomics, proteomics and, more recently, metabolomics are three of these techniques that have impacted on fungal plant pathology. Used individually, each of these techniques can generate a plethora of data that could occupy a laboratory for years. When used in combination, they have the potential to comprehensively dissect a system at the transcriptional and translational level. Transcriptomics, or quantitative gene expression profiling, is arguably the most familiar to researchers in the field of fungal plant pathology. Microarrays have been the primary technique for the last decade, but others are now emerging. Proteomics has also been exploited by the fungal phytopathogen community, but perhaps not to its potential. A lack of genome sequence information has frustrated proteomics researchers and has largely contributed to this technique not fulfilling its potential. The coming of the genome sequencing era has partially alleviated this problem. Metabolomics is the most recent of these techniques to emerge and is concerned with the non-targeted profiling of all metabolites in a given system. Metabolomics studies on fungal plant pathogens are only just beginning to appear, although its potential to dissect many facets of the pathogen and disease will see its popularity increase quickly. This review assesses the impact of transcriptomics, proteomics and metabolomics on fungal plant pathology over the last decade and discusses their futures. Each of the techniques is described briefly with further reading recommended. Key examples highlighting the application of these technologies to fungal plant pathogens are also reviewed.