Skip to main content
  • Research article
  • Open access
  • Published:

A comparative hidden Markov model analysis pipeline identifies proteins characteristic of cereal-infecting fungi

Abstract

Background

Fungal pathogens cause devastating losses in economically important cereal crops by utilising pathogen proteins to infect host plants. Secreted pathogen proteins are referred to as effectors and have thus far been identified by selecting small, cysteine-rich peptides from the secretome despite increasing evidence that not all effectors share these attributes.

Results

We take advantage of the availability of sequenced fungal genomes and present an unbiased method for finding putative pathogen proteins and secreted effectors in a query genome via comparative hidden Markov model analyses followed by unsupervised protein clustering. Our method returns experimentally validated fungal effectors in Stagonospora nodorum and Fusarium oxysporum as well as the N-terminal Y/F/WxC-motif from the barley powdery mildew pathogen. Application to the cereal pathogen Fusarium graminearum reveals a secreted phosphorylcholine phosphatase that is characteristic of hemibiotrophic and necrotrophic cereal pathogens and shares an ancient selection process with bacterial plant pathogens. Three F. graminearum protein clusters are found with an enriched secretion signal. One of these putative effector clusters contains proteins that share a [SG]-P-C-[KR]-P sequence motif in the N-terminal and show features not commonly associated with fungal effectors. This motif is conserved in secreted pathogenic Fusarium proteins and a prime candidate for functional testing.

Conclusions

Our pipeline has successfully uncovered conservation patterns, putative effectors and motifs of fungal pathogens that would have been overlooked by existing approaches that identify effectors as small, secreted, cysteine-rich peptides. It can be applied to any pathogenic proteome data, such as microbial pathogen data of plants and other organisms.

Background

Fungal pathogens represent one of the most important restraints of food production worldwide and can also threaten natural plant populations [1, 2]. Any one plant species can come under attack from a diverse range of fungal pathogens that may not be related. These different pathogens may also have different modes of infection: biotrophs which require living tissue to obtain nutrients, necrotrophs which kill the plant cells to obtain nutrition and hemibiotrophs which show phasic interactions with plant hosts that are a combination of both biotrophy and necrotrophy. The outcome of the interaction between any of these infection modes will depend on both the ability of the plant to defend itself and the ability of the pathogen to cause disease. Understanding how diverse pathogens cause disease on specific hosts will provide information for the development of new crop protection strategies.

Pathogens utilise a wide range of virulence strategies to bring about disease. These include the production of small secreted proteins (or effectors) that interfere with host functions, the production of small secondary metabolites that may have effector-like function or be toxic to host cells and the production of reactive oxygen species. Pathogens also possess mechanisms to defend themselves against the host response to infection such as proteins to mask other molecules that may otherwise be recognised by the host, enzymes to metabolise toxic plant derived compounds and enzymes to detoxify reactive oxygen species. However, our understanding of the virulence arsenal of most pathogens is fragmentary. This is in part because most virulence strategies possessed by a single pathogen contribute in a quantitative manner to the outcome of an interaction, making their detection difficult. Proteinaceous effectors that act as avirulence determinants in specific host genotypes are an exception to a quantitative contribution to an interaction, but the virulence function of these in compatible hosts are also poorly defined.

De novo prediction of pathogen virulence mechanisms in fungi based on common sequence features has to date been difficult. In contrast, proteinaceous effectors of bacteria and oomycete pathogens can be predicted from conserved targeting signals in the N-terminal regions with comparative ease. Using whole-genome sequence and deduced protein-coding gene sequences of these organisms, proteins with certain conserved amino-acid motifs in their N-terminal region have been observed to be strongly associated with predicted secreted proteins. In bacteria, this is the signal that targets proteins to the type III secretion system which directly delivers proteins to the host cytoplasm [35]. Effectors in oomycete secretomes have conserved RXLR-dEER [6], Crinkler [7], CHXC [8] and YXSL[RK] motifs [9] allowing the identification of new Oomycete effectors. The RXLR-dEER motif has been demonstrated to facilitate effector uptake by the host cell [1012]. oomycete genomes encode hundreds of such potential effectors [13], a growing number of which have been functionally validated.

In the fungi, no functional uptake signal has been described yet and thus robust effector prediction based on signature sequence motifs are not feasible. Effector detection in fungal pathogens has been based largely on host specific responses mediated by host recognition of the effectors [14]. Bioinformatics approaches utilize known effector features for analyzing the secretome of fungal pathogens to find proteins with possible roles in virulence. The secretome is typically predicted using a combination of signal peptide prediction methods [15, 16] and accuracy relies heavily on the reliability of the secretion prediction methods used. Effector prediction in the literature commonly involves selecting small cysteine-rich proteins from the predicted secretome as these have been associated with properties of known effectors [15, 17, 18]. Small, cysteine-rich proteins are typically defined as being shorter than 150 amino acids [19] and having more than four cysteine residues [20], however there is no consensus in the literature about how these thresholds are set. Therefore, this type of selective effector prediction approach is dependent on the validity of the thresholds used and also makes a priori assumptions about effector protein properties.

These a priori assumptions about effectors are beginning to be challenged with the discovery of effector proteins that are larger in size, do not contain a high number of cysteines and are not predicted to be secreted based on a signal peptide sequence. For example, a 290 amino acid secreted chorismate mutase enzyme of the Basidiomycete fungus Ustilago maydis is taken up by maize cells where it interferes with the production of the defence-stimulating plant hormone, salicylic acid [21]. In the bacterium Xanthomonas campestris a 536 amino acid uridine monophosphate transferase acts as a virulence effector inibiting host kinases [22]. The SIX3 gene in Fusarium oxysporum encodes an effector protein with only two cysteine residues [23] and the AvrM gene encodes an effector without cysteines in the haustoria-forming pathogen Melampsora lini[24]. The barley powdery mildew fungus genes AVR-a10 and AVRk1 encode effector proteins without a secretion signal peptide [25]. Furthermore, in a proteomic study of the interaction between F. graminearum and wheat, a large proportion of proteins of fungal origin isolated from the apoplastic space between host plant cells were also not predicted to be secreted [26]. The reverse approach of functionally analysing large sets of secreted proteins for roles in virulence has to date also had a relatively low success rates. For example, in an analysis of 78 secreted proteins from Magnaporthe oryzae only one was shown to have a role in virulence [27]. In summary, we propose that fungal pathogens encode many more effectors than those currently known and with additional or alternate molecular properties than those that are predicted by current approaches.

An unbiased and powerful way to predict the virulence arsenal of fungal pathogens is to compare and contrast fungal pathogens via comparative genomic analyses. Fungal pathogens of the same or similar host plant can be expected to have evolved, to some degree, common virulence mechanisms that target the same host components to bring about disease. The relatively small size of fungal genomes (compared to many plant hosts) makes them amenable to whole genome sequencing and the number of fungal genome sequences available provides an unprecedented opportunity to use comparative analyses to associate genes with virulence functions. Cross species comparative analysis has been successfully applied to the identification of completely novel virulence mechanisms in wilt pathogens of dicot hosts [28] and in the cereal-infecting fungus F. pseudograminearum[29]. In both of these examples, proteins uniquely present in organisms sharing a common biological feature (i.e wilt-causing fungi or cereal-infecting fungi, respectively) were shown to be involved in virulence. Both of these studies relied on the sequence similarity search tool BLAST [30]. BLAST is a position-dependent local alignment tool that works well for regions of high sequence similarity. However, more sensitive methods are likely to be required in fast-evolving genome regions where mutations might lead to little sequence identity amongst a protein family, yet the three-dimensional structure, i.e. the functionality, will be preserved. For example, RXLR effector proteins might share as little as 15% sequence identity yet they can have almost identical three-dimensional structures [31]. Sequence similarity search tools such as BLAST will fail to detect those important functional relationships.

More sensitive interrogation techniques for screening for putative virulence and effector molecules in fungal pathogens are hidden Markov models (HMMs). Profile HMMs are an alternative to BLAST and are statistical models of multiple sequence alignments with the ability to capture position-specific information. Certain positions can have specific scores for amino acids as well as specific gap penalties for insertions or deletions. This is a clear mathematical and biological advantage over position-independent local alignment tools such as BLAST. Profile HMMs can group sequences into evolutionarily related families and find remote homologues. The power of profile HMMs becomes obvious when a family of evolutionarily related sequences is known and a carefully curated multiple sequence alignment is available. Instead of scanning a large number of weak BLAST hits by eye for remote homologues, a profile HMM can be calculated from the alignment and can be used to search a database for remote family members. Profile HMM searches are now essentially as fast as BLAST and available in the HMMER3 software package [32]. If no initial multiple sequence alignment is available, a profile HMM can be built from a single sequence and can be used in a database search using phmmer. Despite these advantages, to our knowledge, single sequence profile HMM methods such as phmmer have not been applied to the de novo identification of putative virulence-related proteins and effectors in pathogens.

In this work, we present a sensitive and accurate HMM-based pipeline for predicting proteins with shared functions, in this particular study with putative virulence or effector functions, based on a comparative analysis and a subsequent protein clustering step (Figure 1). We use 174 fungal genomes from the Joint Genome Institute (JGI) MycoCosm [33], which covers diverse pathogenic and non-pathogenic fungi and is a well-established resource in fungal genomics [34]. Instead of restricting ourselves to a small number of genomes from a particular lineage (and thus creating taxonomic bias), we intentionally used 174 publicly available fungal genomes which cover the fungal kingdom in a phylogenetically and phenotypically diverse manner. The initial step of the pipeline is to find proteins in a query fungal pathogen genome which have HMM sequence similarity hits specific to fungal pathogens and at the same time predominantly absent from non-pathogenic fungi. These pathogen-associated proteins are ranked according to their degree of conservation across other fungal pathogens, with and without a cereal host, to find infection mechanisms common to pathogenic fungi. Unsupervised clustering of the pathogen-associated protein set according to 35 sequence-derived features (e.g. signal peptide prediction score, molecular weight, amino acid composition) leads to protein clusters which share similar characteristics. From this we identify clusters that have an enriched secretion signal or conserved potential signalling domain as putative effector protein groups. The major advantage of the unsupervised clustering of proteins for finding effector genes is that no heuristic cut-offs such as protein length or the number of cysteines is used. However, if these are features common to a group of effectors, they will show up in the clusters naturally. Furthermore, the clustering of pathogen-associated proteins using sequence-derived features is a first step towards finding structural homologies in the twilight zone of sequence similarity (20–30% pairwise sequence identity). Our method has revealed promising novel candidate virulence factors as well as putative effector molecules, some of which appear to have unusual evolutionary origins suggesting they have been subjected to specific selection and have important functions in pathogenesis. Furthermore, we have identified protein groups with related sequence-derived features, which may indicate common functionality. Our method can have wide application in effector discovery in microbial pathogens of plants and other organisms.

Figure 1
figure 1

The comparative pipeline for predicting candidate virulence genes and effectors. Our comparative analysis pipeline predicts proteins in a query pathogen genome which have HMM sequence similarity hits predominantly in fungal pathogens. For a protein to be associated with pathogenicity, the list of phmmer hits must include at least 80% hits to proteins from pathogen species (including the query genome itself). These pathogen-associated proteins are then analyzed in terms of two criteria: (1) their degree of conservation across other fungal pathogens with and without a cereal host and (2) their potential to act as effectors outside the fungal cell. Proteins which are highly conserved across a diverse range of other pathogens as identified by our F-measure ranking are prime candidates for virulence-related proteins. To identify putative effector candidates in an unbiased way, an unsupervised clustering technique based on 35 sequence-derived protein features is used to look for protein clusters with an enriched secretion signal. Further investigation of the clusters with regards to predicted sequence motifs, novel secretion signals, amino acid composition and functional annotation is conducted to find and characterise novel putative effector families.

Results and discussion

Development of a candidate virulence gene selection pipeline

To identify genes with potential roles in virulence on plant hosts a prediction pipeline was developed that utilises distribution across the fungal kingdom as a means of identifying candidates. The hypothesis underlying this pipeline is that genes involved in virulence processes are more likely to be retained or uniquely present in organisms that share a biological trait such as pathogenicity towards plants. The comparative analysis pipeline queries a genome of interest by a systematic and sensitive HMM comparison to a large number of publicly available fungal genomes. Here, 174 genome sequences were used, but as new genomes are published they can be added to the pipeline. Our approach intentionally avoids taxonomic groupings of the 174 fungal genomes as pathogens can be phylogenetically closely related yet at the same time be phenotypically diverse, e.g. by infecting different hosts or tissues. For example, the group species Fusarium oxysporum contains plant and opportunistic animal pathogens, host-specific pathogens, non-pathogens and biocontrol agents.

The initial step identifies proteins predicted in the query genome that are enriched in pathogen genomes and are predominantly absent in the genomes of non-pathogens (Figure 1). Thus, this step identifies putative virulence- and pathogenicity-associated proteins. The next step in the pipeline is to rank these pathogen-associated proteins according to their respective frequency of association across pathogens, and this was achieved using the F- measure, which is the harmonic mean of precision and recall. A perfect F-measure of 1 indicates presence uniquely in fungal pathogens and across all fungal pathogens (Figure 2). To investigate host-specific patterns, we report the F-measure with respect to two backgrounds, the frequency of association across pathogens in general (pathogen F-measure) and the frequency of association across pathogens with a cereal host (cereal F-measure). To dissect the nature of the pathogen-associated proteins we next used an unsupervised clustering method based on 35 features (e.g. signal peptide prediction score, molecular weight, amino acid composition), deduced from the sequences. This leads to a suite of clusters of related pathogen-associated proteins. One of our aims was to identify potential effector proteins, which we define as proteins active outside the fungal cell that interact with host defences as well as pathogen recognition and signalling by the host, and these would be expected to represent a subset of the pathogen-associated protein clusters. Traditionally, effector candidates have been identified as small secreted peptides, whereas our model makes no assumption on size, and also includes searches for novel potential secretion signatures, as well as conventional secretion signals. If a family of effectors has sequence-derived features in common, such as the amino acid composition or the predicted signal peptide score, and if these features are distinct from the remaining proteins, it can be expected to form a natural group or cluster. This cluster will then have an enrichment or depletion in its characteristic features from background, which is the set of all pathogen-associated proteins. At this stage, proteins for specific functional and evolutionary analysis can be identified to further elucidate roles in pathogenesis.

Figure 2
figure 2

A visualization of the F -measure for finding highly conserved virulence genes in pathogens, with and without a cereal host. The set of 174 fungal proteomes from the JGI MycoCosm tree [33] is shown with the subsets of pathogenic fungi and cereal-infecting fungi. The genera of cereal-infecting fungi are colored according to fungal lifestyle (biotrophs, hemibiotrophs, necrotrophs). The F-measure is the harmonic mean of precision and recall and ranges from 0 to 1. For a query pathogen protein of interest, a F-measure of 0 means that the protein is only found in the query pathogen genome. A perfect pathogen F-measure of 1 means that this protein is found in all pathogen genomes and at the same time only in pathogenic fungi, thus indicating a protein highly relevant to pathogenicity. In the middle of the scale sit pathogen proteins which can be found mostly in pathogens and at the same time in some non-pathogens. We also introduce the more specific cereal F-measure to identify proteins which are highly conserved and exclusive to cereal-infecting fungi to identify host-specific infection patterns.

Application to S. nodorum and F. oxysporum returns known effectors

F. graminearum is one of the most important pathogens of cereals globally and was the focus for our study. However, no effectors have been functionally validated so far in the literature for this pathogen. Therefore, we initially tested our pipeline on two fungal genomes which have a number of characterized effector proteins, i.e. Stagonospora nodorum SN15 [35] and Fusarium oxysporum f. sp. lycopersici strain 4287 [36]. Note that unlike many bacterial pathogens, no individual fungal pathogen has sufficient numbers of known virulence and non-virulence proteins to allow benchmarking by identifying true positives and false positive rates, which prohibits a large-scale validation of the pipeline.

S. nodorum SN15 is a necrotrophic pathogen of wheat with a number of characterized effectors or host-selective toxins, called SnToxA, SnTox1 and SnTox3. We applied our gene-selection pipeline to the predicted protein set of the S. nodorum genome and found 3528 proteins (28% of gene models) which are enriched in pathogen genomes and are predominantly absent in the genomes of non-pathogens (Additional file 1). The proteins with the highest combined degree of conservation and exclusivity across fungal pathogens according to the pathogen F-measure are all predicted to be secreted (SNOG_12346, SNOG_01958, SNOG_09549, SNOG_13500, SNOG_05633). Two of the proteins have a functional annotation, one as a putative secreted phosphorylcholine phosphatase (SNOG_09549) and another one as a predicted cystatin-like fold (SNOG_13500). The S. nodorum proteins with the highest cereal F-measure are a zinc finger domain containing protein (SNOG_09637), two small secreted cysteine-rich proteins with no functional annotation (SNOG_20011, SNOG_05982) and a putative secreted phosphorylcholine phosphatase (SNOG_09549).

Our predicted set of pathogen-associated S. nodorum proteins includes the known effectors SnToxA (SNOG_16571), SnTox1 (SNOG_20078) and SnTox3 (SNOG_08981). Overall, SnToxA was ranked 65th in the cereal F-measure list. Due to some proteins sharing the same F-measure, the top 65 rank positions account for 636 of the 3528 proteins enriched in pathogen genomes. For the three known effectors SnToxA, SnTox1 and SnTox3, SnToxA has the highest cereal F-measure of 0.35 due to HMM sequence similarity hits to Pyrenophora tritici-repentis and Cochliobolus heterostrophus isolates C4 and C5 [37]. For both C. heterostrophus isolates C4 and C5 proteins a ToxA domain is predicted using the 3D structure prediction software Phyre2 (100% confidence, 44% sequence identity). Interestingly, despite the high level of sequence similarity, the RGD-containing loop [38] is not conserved in C. heterostrophus isolates C4 and C5 (SGN-containing loop, Figure 3). The other two host-selective toxins SnTox1 and SnTox3 are only found in S.nodorum and thus have the lowest possible F- measure. In general, the unsupervised clustering of pathogen-associated proteins returns three clusters with a significantly higher secretion signal than the background distribution for all clusters (Mann–Whitney U test, p-value < 2.2e-16). Interestingly, the three effectors SnToxA, SnTox1 and SnTox3 all belong to the same cluster, which has 191 members in total and is enriched in small and non-polar amino acids (data not shown). The 191 proteins in this cluster share sequence-derived similarites and based on the membership of known effectors there is good potential that this cluster includes proteins that are worth investigating in terms of their effector potential.

Figure 3
figure 3

HMM search returns ToxA-like proteins in Cochliobolus heterostrophus isolates C4 and C5. A phmmer search with SnToxA (SNOG_16571) returns sequence similarity hits to Pyrenophora tritici-repentis (PTRG_04889) and hits to potential ToxA-like domains in two proteins in Cochliobolus heterostrophus isolates C4 and C5. The corresponding multiple alignment of the phmmer domain hits is visualized by Jalview [39].

F. oxysporum f. sp. lycopersici strain 4287 is a hemibiotroph pathogen that causes vascular wilt in tomato plants. Our pipeline returns 4613 proteins (26% of gene models) which are enriched in pathogen genomes and are predominantly absent in the genomes of non-pathogens (Additional file 2). As F. oxysporum is not a cereal-infecting pathogen, we investigated the proteins with highest pathogen F-measure. These include a putative secreted phosphorylcholine phosphatase (FOXG_09820), a putative secreted peptidase C69 domain protein (FOXG_05778) and a putative synaptotagmin-1 protein fold (FOXG_13700). F. oxysporum is known to secrete several proteins during infection of the host and these are referred to as secreted-in-xylem (SIX) proteins, which include several functionally proven effectors [23, 4042]. Our set of predicted pathogen-associated proteins contains all the SIX proteins for this strain, i.e. SIX1 (FOXG_16418), SIX2 (FOXG_16416), SIX3 (FOXG_16398), SIX6 (FOXG_14246), SIX8 (FOXG_16464, FOXG_17445), SIX9 (FOXG_14223) and SIX10 (FOXG_17457). These effectors have lowest possible F- measure due to the lack of hits to any other fungal pathogen. The unsupervised clustering returns three clusters with a significantly higher secretion signal than the background distribution for all clusters (data not shown). For the SIX2 effector, the software SIGNALP 4.0 does not predict a signal peptide, yet in our analysis it is still placed in the same cluster as the secreted SIX1, SIX6, SIX9 and SIX10 effectors based on its sequence-derived features. This might indicate that our clustering approach successfully recognizes structural similarities amongst the pathogen-associated proteins. This particular cluster has 221 members and is enriched in tiny, small and non-polar amino acids and in particular in the amino acids cysteine, glycine and asparagine. Due to the occurence of known effectors in this cluster and the traditional association of some of the cluster characteristics (small, secreted, cysteine-rich) with effectors, this cluster might contain other potential effector candidates. The SIX3 and SIX8 proteins are also not predicted to be secreted by SIGNALP 4.0. SIX3 sits in a cluster with high average positive protein charge and a significantly higher percentage of prolines, serines and threonines than the other clusters. SIX8 is part of a cluster with a strong secretion signal and a significantly higher percentage of aliphatic amino acids than the background distribution for all clusters.

Taken together, the blind study of S. nodorum and F. oxysporum using comparative analysis to the larger fungal genome collection was able to highlight pathogen-associated proteins highly conserved across fungal pathogens using the F-measure, as well as to recover all the proteins that are known effectors and to predict novel effector candidates due to cluster membership. Pathogen-associated proteins with high F-measure, i.e. proteins that score highly for pathogen or cereal-infecting specificity and conservation, are worthy candidates for putative virulence or pathogenicity factors common to a large set of pathogens. On the other hand, the SIX effectors and SnTox1 and SnTox3 have the lowest possible F-measure due to being exclusive to F. oxysporum and S. nodorum, respectively. However, for both genomes we find that these known effectors sit in clusters of proteins with strong signals for being small and secreted. This observation fits with known signatures for effectors and other proteins clustered with these known effectors based on similar sequence features may represent candidates worthy of further investigation. Similarly, other clusters of proteins that contain known effectors but are defined by different sets of sequence-derived features potentially represent uncharacterised effector properties. In summary, the F-measure and clustering methods are complementary approaches to analyze the space of pathogen-associated proteins. Fungal pathogens acting on the same host plant can be expected to utilize a common set of invasion strategies reflected in proteins highly conserved across fungal pathogens. At the same time, there will be effector proteins which are exclusively found in a particular fungal pathogen and as shown in this section, these can be found in an unbiased way using a clustering approach and subsequent protein selection step.

Application to F. graminearum: pathogen-associated proteins highly conserved across fungal pathogens

After validation of the analysis pipeline on two relatively well characterized fungal pathogens (S. nodorum and F. oxysporum) with known effectors, the pipeline was applied to investigate the set of pathogen-associated proteins in one of the most important pathogens of cereals globally, the hemibiotrophic ascomycete F. graminearum[43]. As set in our pipeline, pathogen-associated proteins are returned initially and these are defined by having hits with more than 80% pathogen precision. Therefore, pathogen-associated proteins are proteins which are predominantly absent from the genomes of non-pathogens.

For F. graminearum proteins we find 2830 pathogen-associated proteins which represents 21% of all F. graminearum gene models (Additional file 3). We find that the vast majority (99%) of the 2830 F. graminearum proteins belong to the unclassified category when using MIPS Functional Catalogue Database (FunCatDB) for searching the functional distribution [44]. Only 241 of the 2830 proteins have a protein domain hit when sequence-based tools such as Pfam search are used. 3D structure similarities for 487 proteins of the 2830 proteins are reported by Phyre2 with confidence > 95% [45] and interestingly, 298 of these predicted protein structures have no Pfam domain hit. Thus, predicted 3D structure similarities can in many cases complement sequence-similarity based functional annotations.

The majority of the 2830 pathogen-associated proteins from F. graminearum (90%) have no significant sequence similarity hits to a fungal pathogen outside the Fusarium genus. The remaining 10% of the 2830 pathogen-associated proteins have sequence similarity hits to a diverse range of fungal pathogens. The 2830 pathogen-associated proteins from F. graminearum were ranked according to their degree of conservation across other fungal pathogens, with and without a cereal host, to reveal proteins highly relevant to pathogenicity and virulence. Table 1 shows the top 10 proteins in terms of cereal F-measure and pathogen F-measure, respectively, with their corresponding Pfam domain hits and predicted 3D structure similarities. We also report the taxonomical distribution of sequence similarity hits outside the MycoCosm tree by using an analysis of the phmmer web search with the non-redundant protein databank (NR) [32].

Table 1 The top 10 pathogen-associated F. graminearum proteins in terms of cereal and pathogen F-measure

A protein with a haloacid dehalogenase (HAD)-like hydrolase domain hit (FGSG_03333) has the highest cereal F-measure in our analysis. Using the Phyre2 structure prediction program FGSG_03333 is predicted to encode a phosphorylcholine phosphatase which is a subfamily of the HAD superfamily. Phosphorylcholine phosphatases have been associated with breakdown of the host cell membrane in bacteria [46]. FGSG_03333 is also predicted to be secreted and is orthologous to SNOG_009549 and FOXG_09820 identified in the analysis of the S. nodorum and F. oxysporum genomes (see above). A NCBI BLAST search and phylogenetic tree analysis reveals strong sequence-based matches to bacteria many of which are pathogenic on plants, e.g. Pseudomonas syringae (Figure 4). The I-TASSER server [47] also returns a confident 3D structure prediction (TM-score 0.61 ± 0.14) with the top structural analogue again a phosphorylcholine phosphatase from P. aeruginosa[46]. This P. aeruginosa protein has not previously been reported to be associated with cereal pathogenicity, however the phylogenetic analysis for FGSG_03333 suggests that its presence is characteristic of hemibiotrophic and necrotrophic cereal pathogens and that it may have undergone an ancient selection process with bacterial plant pathogens. Interestingly, we also find two proteins in close genomic proximity (FGSG_03338 and FGSG_03339) in the top 10 cereal F-measure table, which both have hits restricted to the Ascomycota (Table 1).

Figure 4
figure 4

FGSG_03333 encodes a putative phosphorylcholine phosphatase with a shared selection process with bacterial plant pathogens. For protein FGSG_03333, a BLASTp search returns 231 hits with E-value < 1.0e-05 and these were merged with our results from the MycoCosm tree. A MUSCLE alignment [48] and maximum likelihood phylogenetic tree estimation by PhyML [49] returns this phylogenetic tree with branch support values, visualized by FigTree [50]. The tree indicates that this protein has been selectively retained in hemibiotrophic and necrotrophic fungal pathogens of cereals, plants and insects. Furthermore, this protein is closely related to a highly conserved protein in a large number of Pseudomonads, indicating a shared selection process and suggesting a common role in pathogenicity.

The highest pathogen F-measure of 0.65 is shared by two proteins and both are predicted to be secreted (FGSG_04060 and FGSG_09841). Protein FGSG_04060 is confidently predicted to be a beta-expansin 1a-like protein by Phyre2 as well as by the I-TASSER server (TM-score 0.73 ± 0.11), with the top structural analogue being a beta-expansin and group-1 pollen allergen from maize [51]. A public database search with FGSG_04060 returns a large number of significant hits to a wide range of pathogenic fungi including ascomycetes and basidiomycetes as well as to oomycetes and proteobacteria pathogenic on plants (Table 1). Bacterial and fungal expansin-like proteins have been associated with plant pathogenesis or plant cell wall degradation [52, 53] and thus this protein might constitute a virulence factor that is worth functionally characterizing. Protein FGSG_09841 has an equally high pathogen F- measure with hits mainly to ascomycetes, however it has no significant Pfam domain hits or confident 3D structure predictions.

In summary, the ranking of pathogen-associated proteins in terms of their F- measure is a technique that is able to identify proteins in the query genome which are highly conserved across other fungal pathogens and at the same time predominantly absent from non-pathogenic fungi. Proteins with these properties are prime candidates for key factors in pathogenicity or virulence. The use of two different types of rankings, namely the cereal and pathogen F- measure aims to identify proteins that might be specific to pathogens acting on cereal hosts or on plants in general. All of the pathogen-associated proteins reported by us are predominantly absent from saprophyte genomes, however they vary in terms of how their sequence similarity hits are distributed in other kingdoms. For example, we find a large number of hits for FGSG_03333 to the bacterial kingdom or more specifically to bacterial plant pathogens, supporting a role in pathogenicity or virulence. On the other hand, for several of the proteins in Table 1 the taxonomical distribution outside the MycoCosm tree might not give a clear indication as to whether they have a functional role in virulence. For example, some proteins have sequence similarity hits restricted to pathogenic ascomycetes (FGSG_03338, FGSG_03339, FGSG_04015, FGSG_03549, FGSG_11152) whereas others contain domains that are widespread across the bacterial and eukaryotic kingdoms (FGSG_04507, FGSG_07909, FGSG_07846, FGSG_04060, FGSG_11496). Proteins are generally composed of several domains and the combination of the domains determines the function of the protein. Therefore, the presence of one domain that is widespread across pathogenic and non-pathogenic organisms, e.g. the F-box domain, and thus results in a large number of hits to diverse species, is not enough evidence to dismiss a protein as not having a role in pathogenicity or virulence.

Clustering of the pathogen-associated proteins in F. graminearum reveals groups of putative effectors

To find putative F. graminearum effectors in the set of pathogen-associated proteins, we applied unsupervised clustering based on 35 sequence-derived features. As is required by the k-means approach, the number of clusters has to be predefined and it was set at 12, based on estimates using the elbow plot method (see Methods). The feature attributes of each cluster were compared to the full background distribution to identify those features that statistically discriminated that cluster from background. Clusters were examined further based on likely or established models of effector characteristics, for example a strong predicted secretion signal, potential novel signalling domain or a striking amino acid composition.

Table 2 shows the 12 clusters (numbered from C1 to C12) and their associated properties that were found using our k- means clustering approach (Additional file 4). We find three clusters of proteins (C2, C11, C12) with an enriched positive protein charge that are distinct from each other in terms of their amino acid composition. Cluster C12 contains proteins with a high percentage of aromatic amino acids, whereas cluster C2 features tiny amino acids and cluster C11 has a large proportion of polar amino acids. On the other hand, five clusters show depletion in positive protein charge and thus contain proteins with low charge (C3, C5, C8, C9, C10). Two of these clusters (C3, C9) contain a significantly higher number of aliphatic amino acids that are very hydrophobic and tend to stay in the interior of a protein structure. Cluster C5 shows depletion in cysteines and cluster C10 is enriched in aromatic amino acids. Another cluster (C6) contains 166 proteins with significantly higher molecular weight than the background distribution for all clusters, however no other discriminative features are evident. Cluster C4 has 210 proteins with enrichment of the amino acid proline. Proline has unusual biochemical properties and has been associated with a variety of functions, such as membrane proteins, cell adhesion or cereal storage proteins. Three clusters (C1, C7, C8) show an enriched secretion signal predicted by SignalP and WoLF PSORT, and are thus likely candidates for proteins which are secreted from the pathogen to the host apoplast.

Table 2 Selected properties for the clusters of pathogen-associated proteins

In general, we find that our k-means clustering approach divides the space of pathogen-associated proteins into groups of proteins with distinct sequence-derived features. As expected, some clusters have more discriminative features than others. We anticipate the set of pathogen-associated proteins to contain protein groups which share a number of features that are clearly distinct from background, as for example seen for small, secreted, cysteine-rich proteins in cluster C7. On the other hand, there will be protein groups which have been assigned to a cluster based on a less discriminative feature set, such as cluster C6 which contains proteins with high molecular weight and no other characteristic feature. Taken together, k-means clustering is an unsupervised method which will find natural groups of proteins with discriminative features in the set of pathogen-associated proteins. Unlike methods for effector prediction which use a number of thresholds and assumptions for selecting putative effectors, k-means clustering will return putative effector clusters in an unbiased way, if discriminative characteristics for effector proteins exist.

Examination of putative effectors reveals extracellular domains and conserved N-terminal sequence motifs

From the clustering of pathogen-associated proteins and the extraction of significant cluster features, we identified three clusters (C1, C7, C8) of proteins with an enriched predicted secretion signal, thus prime candidates for containing putative effector proteins (Table 2). These three clusters differ significantly in terms of their protein characteristics and amino acid composition. Cluster C1 contains 268 proteins with enrichment in the secretion signal as well as in the non-polar, aromatic and aliphatic amino acid content. Aliphatic amino acids have been associated with transmembrane regions and indeed the TMHMM program [54] predicts transmembrane regions for 98 of the 268 proteins. Only 9 of the 268 proteins have a significant Pfam domain hit, e.g. a copper amine oxidase (FGSG_01761), a carbon-nitrogen hydrolase (FGSG_07204), an ion-channel (FGSG_11239) and a ferroportin 1 (FGSG_13446). The two other protein clusters (C7, C8) are enriched in tiny and small amino acids and are at the same time predicted to be secreted, which are features traditionally associated with effector proteins.

Cluster C7 contains 150 proteins and is enriched in amino acids cysteine and glycine. Cysteines are commonly associated with extracellular proteins due to their ability to form stable disulfide bonds. Glycine is the smallest amino acid and has no side chain. Thus glycine often plays a structural role in sterically restrictive turn regions. Significant Pfam domain hits are found for only 14 of the 150 proteins and some have associations with pathogenicity in the literature. FGSG_03573 has a predicted CFEM domain which has been found in extracellular fungal membrane proteins and has been associated with fungal pathogenesis, for example in Magnaporthe grisea[55]. Another protein FGSG_03109 has a cysteine-rich secretory Pfam domain hit and has predicted structural similarity to the plant PR-1 class of pathogen-related proteins. FGSG_00029 has a LysM domain which has been associated with carbohydrate and chitin binding and has been shown to be part of effector proteins in Mycosphaerella graminicola[56] and Cladosporium fulvum[57]. Phyre2 returns confident structure predictions for 31 of the 150 proteins, e.g. for putative beta expansions 1a (FGSG_04060, FGSG_11496) or a putative cellulose 1,4-beta-cellobiosidase fold (FGSG_07728). MEME does not find a motif common to more than five proteins from this set, except for the signal peptide. For example, FGSG_02841, FGSG_07766 and FGSG_08175 share common motifs as well as the three proteins FGSG_00111, FGSG_00112 and FGSG_13443. In summary, predicted extracellular protein domains and the enrichment in cysteines, small amino acids and the secretion signal form strong evidence for the presence of potential effector proteins in this cluster.

Cluster C8 contains 110 proteins and has sequence-derived features that are distinct from Cluster C7. A striking difference to the other clusters is the enriched threonine content and molecular weight as well as the depleted leucine content and protein charge (Table 1). Interestingly, enrichment in serines, threonines and depletion in leucine has been reported for type III secreted effectors of animal and plant pathogens [3]. Significant Pfam domain hits are found for only 6 of the 110 proteins, i. e. for a WSC domain involved in carbohydrate binding (FGSG_09755, FGSG_10335, FGSG_11313, FGSG_11507), a pectate lyase domain (FGSG_13834) and a CFEM domain (FGSG_02109). Phyre2 returns confident structure predictions for 22 of the 110 proteins, e.g. a putative cellulose 1,4-beta-cellobiosidase fold (FGSG_02999) or a putative pectin lyase-like fold (FGSG_13834). For 37 of the proteins in this cluster, MEME returns a [SG]-P-C-[KR]-P motif (e-value 3.4e-091) adjacent to the predicted signal peptide cleavage site followed by a clear stretch of serines/threonines (Figure 5). We ran an iterative HMM search (jackhmmer) on the NR database with a MUSCLE alignment of the first 50 amino acids of a core set of well-aligned 31 sequences to look for the prevalence of this motif. After five iterations, we find 88 significant hits of the motif to proteins in only F. graminearum, F. pseudograminearum, F. oxysporum and F. solani. Subsets of these proteins share a common domain structure, however with no functional annotation for the majority of hits. A “G-P-C-R-P” motif was initially reported for F. graminearum in the Supplementary Material of Ma et al.[36], however it has not been functionally characterized. Zhang et al.[58] refer to this particular F. graminearum protein group as “CPGRP-anchored” (sic) and report transient induction at 16 h after inoculation during wheat coleoptile infection and at 2 to 8 h during in vitro germination. We suggest that this intriguing Fusarium-specific motif should be referred to as the [SG]-P-C-[KR]-P motif and that functional testing is necessary to reveal its role in the host-pathogen interaction.

Figure 5
figure 5

A highly conserved N-terminal [SG]-P-C-[KR]-P motif in a subset of F. graminearum proteins. MEME returns a [SG]-P-C-[KR]-P motif (e-value 3.4e-091) adjacent to the predicted signal peptide cleavage site for 37 of the 110 proteins in cluster C8. A MUSCLE alignment of the N-terminal region of a core set of well-aligned 31 of those proteins is visualized by Jalview. The MEME logos are shown with the [SG]-P-C-[KR]-P motif followed by a serine-/threonine-rich stretch. An iterative public database search of these aligned sequences (jackhmmer) returned additional members of this putative protein family in F. graminearum, F. pseudograminearum, F. oxysporum and F. solani.

For seven of the proteins lacking the [SG]-P-C-[KR]-P motif in Cluster C8, MEME predicts an alternative motif [WYF]-C-x-T-Y-x-S-T-Y-L about ten amino acids after the predicted N-terminal signal peptide (Figure 6). These 7 proteins are around 600 aas in length, have low average charge, and lack a functional annotation and confident 3D structure prediction. An iterative HMM search (jackhmmer) on the NR database with a MUSCLE alignment of the first 60 amino acids converged after three iterations. We collected 20 hits which are restricted to F. graminearum, F. pseudograminearum, F. oxysporum, F. solani and the cereal pathogen Colletotrichum graminicola. MEME predicts a highly conserved domain structure in the N-terminal region for these proteins (Additional file 5).

Figure 6
figure 6

A protein family with a conserved N-terminal [WYF]-C-x-T-Y-x-S-T-Y-L motif in F. graminearum . MEME returns a [WYF]-C-x-T-Y-x-S-T-Y-L motif (e-value 1.2e-028) following the predicted signal peptide cleavage site for seven of the proteins in cluster C8. A MUSCLE alignment of the truncated 7 sequences in cluster C8 is visualized by Jalview.

In summary, we have identified three protein clusters with an enrichment in the signal peptide prediction score for F. graminearum. These three clusters have distinct molecular features. Cluster C1 contains a high number of aliphatic amino acids and transmembrane proteins. The presence of extracellular domains and enriched cysteine content in cluster C7 indicates the occurrence of potential effectors in the set of 150 proteins. Cluster C8 has enriched threonine content, depleted leucine content and features proteins with high molecular weight and low charge. For proteins in this cluster, we predict two distinct conserved N-terminal motifs ([SG]-P-C-[KR]-P and [WYF]-C-x-T-Y-x-S-T-Y-L) which will be our prime candidates for functional characterisation in future work.

F. graminearum putative effector clusters sit in regions of genome innovation

Subtelomeric and ancient sub-telomeric regions of the F. graminearum genome have been previously described as regions of genome innovation and are characterised by higher levels of genetic diversity between strains and other Fusarium compared to the rest of the genome, higher density of secreted and in planta expressed proteins and higher genetic recombination rates per unit length of DNA [43, 59]. Pathogenicity genes from the Pathogen-Host Interactions database (PHI-base [60]) have been negatively associated with these regions [61]. However, these pathogenicity genes are dominated by highly conserved components of fungal biology such as signalling cascades which when mutated render the fungus with many pleiotropic phenotypes. Genes associated with virulence may have different distributions across the genome compared to pathogenicity genes. Analysis of the clusters of genes identified in this work reveals distinct distributions across the F. graminearum chromosomes (Figure 7). Members of cluster C1, which is enriched in secretion signals and features a number of proteins with transmembrane regions, are found throughout the genome (data not shown). In contrast, members of putative effector clusters C7 and C8 and especially proteins containing the [SG]-P-C-[KR]-P motif are found in regions of genome innovation as depicted in Figure 7.

Figure 7
figure 7

Distribution of predicted protein clusters across the F. graminearum genome. A Circos plot [62] is shown which visualizes the four F. graminearum chromosomes with ticks in kilobase (Kbp) units. The following bands for F. graminearum are visualized: (I) Recombination frequency (blue bars) and SNP density (line), (II) gene density, (III) GC content, (IV) set of 2830 pathogen-associated proteins, (V) protein cluster C7, (VI) protein cluster C8, (VII) proteins with the [SG]-P-C-[KR]-P motif. For each gene set, the location of the genes on the chromosomes and a heatmap in red shading for the gene count in 500 Kbp bins are shown. It can be seen that proteins from putative effector clusters C7 and C8 occur predominantly in regions of genome innovation.

De novo motif search reveals potential signalling, targeting or uptake motifs in F. graminearum

The aims of our study were to find the set of pathogen-associated proteins in F. graminearum and to dissect them in terms of their pathogenicity and virulence potential as well as their effector potential in an unbiased way. In this section, we applied de novo motif search to look for novel effector sequence signatures in the set of pathogen-associated proteins. Effector proteins are active outside the fungal cell and can be expected to carry a secretion signal. Bioinformatics signal peptide prediction programs such as SignalP [63] use neural networks and HMMs for predictions and thus rely on the quality of the training data. Especially for fungal pathogens, the data on known secretion signals or host-translocation motifs is scarce and thus, we applied de novo motif search to look for novel sequence signatures. We expect sequence signatures comprising potential secretion signals and host-translocation motifs to occur within the 150 N-terminal residues of a protein. We apply the motif search tool MEME to each predicted cluster and report motifs which occur on average within the first 150 residues of the proteins and are present in more than five proteins.

To validate our approach for finding conserved N-terminal sequence signatures with potential function in signalling, targeting or uptake, we applied our pipeline to the powdery mildew fungus Blumeria graminis f. sp. hordei which infects barley [64]. An N-terminal Y/F/WxC-motif has been described in powdery mildew fungal effector candidates [65]. We applied MEME to the pathogen-associated protein clusters predicted by our pipeline in B. graminis f. sp. hordei (Additional file 6). We find that only one of the clusters has a significantly higher secretion signal than the other clusters and for 98 of the proteins in this cluster, MEME predicts an N-terminal Y/F/WxC-motif (Additional files 7 and 8). Our set of 98 proteins predicted from the B. graminis f. sp. hordei is of similar size to the set of 107 proteins identified by Godfrey et al.[65] from a set of around 3000 genes from sequenced ESTs. After positive validation of our motif search approach in the space of pathogen-associated protein clusters in B. graminis f. sp. hordei, we applied it to F. graminearum.

For F. graminearum, the de novo sequence motif prediction results for the protein N-terminal regions within the 12 clusters are shown in Table 3. For each of the motifs found in the clusters, we built a HMM from the multiple alignment computed by MUSCLE and used an iterative, sensitive jackhmmer search against the non-redundant protein databank (NR), with a maximum of five rounds, to record the prevalence of the motifs. As reported in the previous section, we find two motifs in cluster 8 which are largely restricted to pathogenic Fusarium and occur, on average, within the first 50 residues of the proteins and after a predicted signal peptide. Another motif sits in cluster C12 and shares resemblance to the L-P-x-E motif which occurs in F-box like domains. An iterative jackhammer search finds that this motif is exclusive to Ascomycota with hits predominantly to fungal pathogens, for example most prominently to S. nodorum or Pyrenophora tritici-repentis. In future work we will experimentally test these motifs for their functional role in the host-pathogen interaction.

Table 3 De novo N-terminal motif search results for the pathogen-associated protein clusters in F. graminearum

Conclusions

The bioinformatics search for fungal pathogen proteins involved in attack on the host plant has thus far focussed on finding small and cysteine-rich peptides in the predicted fungal secretome [20]. Therefore, such methods have been limited in their ability to return putative effectors with those previously recognized features. However, there is growing evidence in the literature that effectors can have unconventional characteristics, such as no predicted signal peptide, a low number of cysteine residues or a large size [2123, 25]. Here, we take advantage of the availability of sequenced fungal genomes and present a novel and unbiased method for isolating genes strongly associated with pathogenic species using comparative genomics analysis followed by an unsupervised clustering step using multiple protein properties. This type of approach allows us to predict putative pathogen-associated proteins or effectors with previously unrecognized characteristics. The main contributions of our method are as follows.

Beyond BLAST: protein comparisons in fast-evolving genome regions

For finding proteins in a query pathogen genome which are conserved across other pathogens and predominantly absent from non-pathogens, we use a HMM methodology [32] as a sensitive method for finding pathogen-associated proteins. Sensitive HMM methods are particularly well-suited for comparative analysis in fast-evolving genome regions, such as pathogenicity islands or dispensable chromosomes. Unlike previous studies which made use of the stringent reciprocal best BLAST hit method [28, 29], our pipeline is able to detect remote sequence similarity relationships across a large number of genomes and protein clusters sharing domains and common structural features derived from their peptide sequences. The reciprocal best BLAST hit method and the triangular extension of clusters of orthologous groups (COGs) are efficient tools for deriving orthology across genomic sequences, however they are inadequate for large evolutionary distances as found in fast-evolving genome regions [66]. Orthology is a powerful concept, however ideally phylogenetic methods should be used in fast-evolving pathogen genomes instead of reciprocal best BLAST hits and COGs and for these, a resolved fungal tree of life would be a prerequisite.

F-measure ranking for finding virulence mechanisms common and exclusive to pathogens

We introduce a novel way for ranking pathogen-associated proteins in a query genome according to their degree of conservation across other fungal pathogens, with and without a cereal host, using the F-measure. This allows for the identification of fungal pathogen proteins that are widely distributed across other pathogenic fungi and at the same time predominantly absent from non-pathogenic fungi. Using this methodology, we were able to identify a highly ranked putative secreted phosphorylcholine phosphatase in F. graminearum for which a phylogenetic analysis suggests that it is characteristic of necrotrophic and hemibiotrophic pathogens and has undergone a shared selection process with bacterial plant pathogens. Proteins which are highly ranked according to our F-measure method are prime candidates for proteins essential for a pathogenic lifestyle and virulence.

Clustering approach returns known effectors and novel effector candidates

We use a k-means unsupervised clustering approach to divide the space of pathogen-associated proteins into natural groups with shared sequence-derived features. This is a step towards finding structural similarities within the pathogen-associated proteins in the twilight zone of sequence similarity. By using an unsupervised clustering approach on the set of pathogenic proteins without the use of heuristic thresholds, we can predict putative effectors which would have gone undetected with pipeline approaches that rely on heuristic assumptions. We demonstrate that our pipeline returns known effectors in S. nodorum and F. oxysporum. For F. graminearum, we predict three clusters with a high secretion signal that contain prime candidates for effector proteins. We found that the proteins of two of these putative effector clusters are found predominantly in regions of genome innovation. One of the clusters is enriched in threonines and molecular weight as well as depleted in leucines and protein charge. Enrichment of threonines and depletion of leucines is a characteristic of type III secreted proteins in animals and plants [3]. Within this cluster, we identified two putative protein families that each have a distinct conserved sequence motif adjacent to the predicted signal peptide cleavage site. One of them is a highly conserved N-terminal [SG]-P-C-[KR]-P motif exclusive to Fusarium with unknown function. We found several other N-terminal motifs conserved in proteins lacking a predicted secretion signal, e.g. a L-P-x-E-R motif with remote homology to F-box domains that is highly conserved in pathogenic Ascomycota such as S. nodorum. In future work we will experimentally test these motifs for their functional role, i.e. whether they might be signalling, subcellular targetting or uptake motifs in the host-pathogen interaction.

Quality of genomic data and application to other genomic data sets

Our bioinformatics pipeline can be applied to any pathogenic genome for which a contrasting set of non-pathogenic relatives is available, such as publicly available fungal or bacterial genomes. It should be emphasized that accuracy and sensitivity of any comparative genomics analysis relies heavily on the quality of genome assembly and subsequent gene prediction and annotation. Databases housing a large number of genome projects such as the JGI MycoCosm will inevitably feature entries of diverse annotation quality. Gene boundary prediction errors are of particular concern for effector candidates which are expected to harbour N-terminal motifs such as signal peptides. However, unlike fungal effector prediction pipelines that strictly select effector candidates based on a predicted signal peptide, our clustering approach can to some extent compensate for incorrect gene boundaries. For example, the F. oxysporum SIX2 effector for which the software SIGNALP 4.0 does not predict a signal peptide is still placed in the same cluster as the secreted SIX1, SIX6, SIX9 and SIX10 effectors based on its sequence-derived features. This might indicate that our clustering approach successfully recognizes structural similarities amongst the pathogen-associated proteins and can offset sporadic gene boundary prediction errors. However, improvements in gene assembly tools and automated gene prediction pipelines are of utmost importance and our method will naturally be able to take advantage of the increasing number of fungal genomes available and of improved genome annotation quality.

Before taking promising effector candidates to experimental validation, one can incorporate expression data of in planta infection such as microarray or RNAseq data. The pipeline scripts are available from the authors on request. With regard to the cereal pathogen genome set, already pre-processed datasets are available for further analysis.

Methods

Fungal genomes and proteomes

Predicted protein sets from 172 fungal genomes were downloaded from the DOE Joint Genome Institute MycoCosm web site [33] and supplemented with the B. graminis f. sp. hordei genome [64] and our in-house genome assembly for F. pseudograminearum CS3096 (GenBank: AFNW00000000.1, [29]). The majority of MycoCosm genomes are of non-pathogenic fungi. More specifically, 56 of the 174 fungal genomes belong to pathogenic fungi and 19 are fungal pathogens of a cereal host (Table 4).

Table 4 The set of 19 fungal pathogen genomes used in our comparative study which have a cereal host

One-against-all proteome HMM comparisons

Sensitive proteome comparisons are used to identify proteins in the query genome of interest which are conserved across other fungal genomes and share a similar protein domain structure. One of the pathogenic fungal genomes is set as the query genome of interest. For the protein set of the query genome, the set of significant protein hits across the 174 fungal genomes is calculated using the profile HMM search tool phmmer [32]. This returns a list of hits with full-sequence E-values for each protein from the query genome. Proteins are composed of one or more domains. There can be multiple matches across a query, i.e. protein domain hits, and these are returned by phmmer as a table of per-domain hits. A query protein is said to have a significant match in another fungal genome if there is a full-sequence hit with E-value < 10-5 and if the combined domain hits cover more than 60% of the query and target sequences, respectively. For each query protein, its significant fungal genome matches are recorded.

Comparative analysis: pathogen-associated proteins

As a result of the phmmer searches, each protein in the query genome has a number of protein hits across the list of 174 fungal genomes, including to itself. To decide if a query protein x i is specific to fungal pathogens and highly conserved across fungal pathogens, four measurements of the relevance of the phylogenetic distribution of the phmmer hits with respect to pathogenicity are calculated: cereal-recall Rc (i), pathogen-recall Rp (i), cereal-precision Pc (i) and pathogen-precision Pp (i). These values are defined as follows:

R c i = 100 × hits to cereal pathogens cereal pathogens ,
(1)
R p i = 100 × hits to fungal pathogens fungal pathogens ,
(2)
P c i = 100 × hits to cereal pathogens hits ,
(3)
P p i = 100 × hits to fungal pathogens hits .
(4)

A high recall means that a method returns the majority of relevant results. High precision means that a method returns more relevant results than irrelevant results. If a protein has high cereal- or pathogen-recall it means that its phmmer hits cover the majority of fungal pathogens with or without a cereal host, respectively. If a protein has high cereal- or pathogen-precision it means that its phmmer hits mostly contain fungal pathogens with or without a cereal host, respectively. Precision and recall can be combined into one measurement, the F-measure:

F = 2 × precision × recall precision + recall .
(5)

For each query protein x i , recall and precision for its hits across the 174 fungal genomes are recorded. Proteins x i are kept with pathogen-precision of Pp (i) > 80%, i.e. at least 80% of its phmmer hits are to other pathogenic fungi, and we call these pathogen-associated proteins. These pathogen-associated proteins are ranked according to their F-measure in terms of cereal and fungal pathogen hits.

Clustering of pathogen-associated proteins

Centroid-based clustering with the k-means method was used to find clusters of similar query proteins in the set of pathogen-associated proteins S. Each protein x i S has a 35-dimensional feature vector based on its sequence information. Let the feature vector v i for protein be given as

v i = v i , s i , e i , m i , h i , l i , p i 1 , , p i 9 , a i 1 , , a i 20 ,

where |v i | is the protein length. Let s i be the D-score returned by SIGNALP 4.0, which corresponds to the prediction whether the sequence contains a signal peptide or not [63]. Let e i be the score for extracellular localization site prediction calculated by WoLF PSORT [72]. The molecular weight of a protein is given by m i , the charge of a protein is given by h i and the isoelectric point is given by l i . Furthermore, the chemical properties of the amino acids in the protein are used as a first step to detect structural homologues with little sequence identity [73]. Let p i 1 , , p i 9 be the percentage of amino acids in a protein which are classified as tiny, small, aliphatic, aromatic, non-polar, polar, charged, basic and acidic, respectively. These features are calculated using pepstats from the EMBOSS package [74]. Let a i 1 , , a i 20 be the percentage of amino acid residues in a protein. The entries in the 35-dimensional feature vectors are numerical and were normalized to avoid a bias towards a feature with higher range. Centroid-based clustering with the k-means method requires the number of clusters to be specified in advance and we use the heuristic elbow plot method for this purpose. For each feature in the 35-dimensional feature vector, Mann–Whitney U tests were used to test whether the distribution within a cluster is identical to the full background distribution, i.e. all clusters, and highly significant p-values for both directions (lesser and greater, p-value < 2.2e-16) are recorded.

Functional annotation and motif search

3D protein structure predictions were obtained using the remote homology recognition technique Phyre2 [45]. Phyre2 is able to predict 3D structure similarities in cases when sequence similarity search tools fail. For Phyre2, a confidence threshold of 95% was used. For selected proteins, the computationally more demanding I-TASSER server was used [47]. The TM-score returned by I-TASSER measures the structural similarity between two structures. A TM-score < 0.17 means that the result is of random similarity. A TM-score > 0.5 indicates that the model is of correct topology. The C-score returned by I-TASSER is a confidence score for estimating the quality of predicted models by and is typically in the range of from −5 to 2. A high C-score relates to a model with a high confidence. Protein domain searches were performed using a local Pfam search for each query protein against Pfam’s library of profile HMMs [32]. Hits with E-value < 10-5 are considered as significant hits. The tool MEME was used for sequence motif discovery [75]. MEME was run with a minimum and maximum motif width of 4 and 50, respectively, and set to return the best 30 motifs.

Implementation of the comparative HMM analysis pipeline

All analyses described above were conducted with custom Python and R scripts. The k-means clustering was performed using SciPy and the Mann–Whitney U tests were carried out using R.

Availability of supporting data

The following additional data are available with the online version of this paper. Additional file 1 is a table listing the properties of the predicted pathogen-associated proteins in S. nodorum. Additional file 2 is a table listing the properties of the predicted pathogen-associated proteins in F. oxysporum. Additional file 3 is a table listing the properties of the predicted pathogen-associated proteins in F. graminearum as well as the individual clusters and predicted motifs. Additional file 4 contains box plots for the 12 clusters in F. graminearum and the 35 protein features used for the clustering. Additional file 5 contains the full sequence alignment for the proteins containing the motif [WYF]-C-x-T-Y-x-S-T-Y-L. Additional file 6 is a table listing the properties of the predicted pathogen-associated proteins in B. graminis. Additional files 7 and 8 contain information about the predicted Y/F/WxC-motif proteins B. graminis.

References

  1. Strange RN, Scott PR: Plant disease: a threat to global food security. Annu Rev Phytopathol. 2005, 43: 83-116. 10.1146/annurev.phyto.43.113004.133839.

    Article  CAS  PubMed  Google Scholar 

  2. Dean R, Van Kan JA, Pretorius ZA, Hammond-Kosack KE, Di Pietro A, Spanu PD, Rudd JJ, Dickman M, Kahmann R, Ellis J: The Top 10 fungal pathogens in molecular plant pathology. Mol Plant Pathol. 2012, 13 (4): 414-430. 10.1111/j.1364-3703.2011.00783.x.

    Article  PubMed  Google Scholar 

  3. Arnold R, Brandmaier S, Kleine F, Tischler P, Heinz E, Behrens S, Niinikoski A, Mewes HW, Horn M, Rattei T: Sequence-based prediction of type III secreted proteins. PLoS Pathog. 2009, 5 (4): e1000376-10.1371/journal.ppat.1000376.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Galan JE, Wolf-Watz H: Protein delivery into eukaryotic cells by type III secretion machines. Nature. 2006, 444 (7119): 567-573. 10.1038/nature05272.

    Article  CAS  PubMed  Google Scholar 

  5. Samudrala R, Heffron F, McDermott JE: Accurate prediction of secreted substrates and identification of a conserved putative secretion signal for type III secretion systems. PLoS Pathog. 2009, 5 (4): e1000375-10.1371/journal.ppat.1000375.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Bhattacharjee S, Hiller NL, Liolios K, Win J, Kanneganti TD, Young C, Kamoun S, Haldar K: The malarial host-targeting signal is conserved in the Irish potato famine pathogen. PLoS Pathog. 2006, 2 (5): e50-10.1371/journal.ppat.0020050.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Torto TA, Li S, Styer A, Huitema E, Testa A, Gow NA, van West P, Kamoun S: EST mining and functional expression assays identify extracellular effector proteins from the plant pathogen Phytophthora. Genome Res. 2003, 13 (7): 1675-1685. 10.1101/gr.910003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Kemen E, Gardiner A, Schultz-Larsen T, Kemen AC, Balmuth AL, Robert-Seilaniantz A, Bailey K, Holub E, Studholme DJ, Maclean D: Gene gain and loss during evolution of obligate parasitism in the white rust pathogen of Arabidopsis thaliana. PLoS Biol. 2011, 9 (7): e1001094-10.1371/journal.pbio.1001094.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Levesque CA, Brouwer H, Cano L, Hamilton JP, Holt C, Huitema E, Raffaele S, Robideau GP, Thines M, Win J: Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire. Genome Biol. 2010, 11 (7): R73-10.1186/gb-2010-11-7-r73.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Dou D, Kale SD, Wang X, Jiang RH, Bruce NA, Arredondo FD, Zhang X, Tyler BM: RXLR-mediated entry of Phytophthora sojae effector Avr1b into soybean cells does not require pathogen-encoded machinery. Plant Cell. 2008, 20 (7): 1930-1947. 10.1105/tpc.107.056093.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Whisson SC, Boevink PC, Moleleki L, Avrova AO, Morales JG, Gilroy EM, Armstrong MR, Grouffaud S, van West P, Chapman S: A translocation signal for delivery of oomycete effector proteins into host plant cells. Nature. 2007, 450 (7166): 115-118. 10.1038/nature06203.

    Article  CAS  PubMed  Google Scholar 

  12. Tyler BM, Kale SD, Wang Q, Tao K, Clark HR, Drews K, Antignani V, Rumore A, Hayes T, Plett JM: Microbe-independent entry of oomycete RxLR effectors and fungal RxLR-like effectors into plant and animal cells is specific and reproducible. Mol Plant Microbe Interact. 2013, 26 (6): 611-616. 10.1094/MPMI-02-13-0051-IA.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Haas BJ, Kamoun S, Zody MC, Jiang RH, Handsaker RE, Cano LM, Grabherr M, Kodira CD, Raffaele S, Torto-Alalibo T: Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature. 2009, 461 (7262): 393-398. 10.1038/nature08358.

    Article  CAS  PubMed  Google Scholar 

  14. Stergiopoulos I, de Wit PJ: Fungal effector proteins. Annu Rev Phytopathol. 2009, 47: 233-263. 10.1146/annurev.phyto.112408.132637.

    Article  CAS  PubMed  Google Scholar 

  15. Brown NA, Antoniw J, Hammond-Kosack KE: The predicted secretome of the plant pathogenic fungus Fusarium graminearum: a refined comparative analysis. PLoS One. 2012, 7 (4): e33731-10.1371/journal.pone.0033731.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Mueller O, Kahmann R, Aguilar G, Trejo-Aguilar B, Wu A, de Vries RP: The secretome of the maize pathogen Ustilago maydis. Fungal Genet Biol. 2008, 45 (Suppl 1): S63-S70.

    Article  CAS  PubMed  Google Scholar 

  17. Duplessis S, Cuomo CA, Lin YC, Aerts A, Tisserant E, Veneault-Fourrey C, Joly DL, Hacquard S, Amselem J, Cantarel BL: Obligate biotrophy features unraveled by the genomic analysis of rust fungi. Proc Natl Acad Sci USA. 2011, 108 (22): 9166-9171. 10.1073/pnas.1019315108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Saunders DG, Win J, Cano LM, Szabo LJ, Kamoun S, Raffaele S: Using hierarchical clustering of secreted protein families to classify and rank candidate effectors of rust fungi. PLoS One. 2012, 7 (1): e29847-10.1371/journal.pone.0029847.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Kamoun S: A catalogue of the effector secretome of plant pathogenic oomycetes. Annu Rev Phytopathol. 2006, 44: 41-60. 10.1146/annurev.phyto.44.070505.143436.

    Article  CAS  PubMed  Google Scholar 

  20. Ellis JG, Rafiqi M, Gan P, Chakrabarti A, Dodds PN: Recent progress in discovery and functional analysis of effector proteins of fungal and oomycete plant pathogens. Curr Opin Plant Biol. 2009, 12 (4): 399-405. 10.1016/j.pbi.2009.05.004.

    Article  CAS  PubMed  Google Scholar 

  21. Djamei A, Kahmann R: Ustilago maydis: dissecting the molecular interface between pathogen and plant. PLoS Pathog. 2012, 8 (11): e1002955-10.1371/journal.ppat.1002955.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Feng F, Yang F, Rong W, Wu X, Zhang J, Chen S, He C, Zhou J-M: A Xanthomonas uridine 5’-monophosphate transferase inhibits plant immune kinases. Nature. 2012, 485 (7396): 114-118. 10.1038/nature10962.

    Article  CAS  PubMed  Google Scholar 

  23. Houterman PM, Speijer D, Dekker HL, DEK CG, Cornelissen BJ, Rep M: The mixed xylem sap proteome of Fusarium oxysporum-infected tomato plants. Mol Plant Pathol. 2007, 8 (2): 215-221. 10.1111/j.1364-3703.2007.00384.x.

    Article  CAS  PubMed  Google Scholar 

  24. Catanzariti AM, Dodds PN, Lawrence GJ, Ayliffe MA, Ellis JG: Haustorially expressed secreted proteins from flax rust are highly enriched for avirulence elicitors. Plant Cell. 2006, 18 (1): 243-256. 10.1105/tpc.105.035980.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Ridout CJ, Skamnioti P, Porritt O, Sacristan S, Jones JD, Brown JK: Multiple avirulence paralogues in cereal powdery mildew fungi may contribute to parasite fitness and defeat of plant resistance. Plant Cell. 2006, 18 (9): 2402-2414. 10.1105/tpc.106.043307.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Paper JM, Scott-Craig JS, Adhikari ND, Cuomo CA, Walton JD: Comparative proteomics of extracellular proteins in vitro and in planta from the pathogenic fungus Fusarium graminearum. Proteomics. 2007, 7 (17): 3171-3183. 10.1002/pmic.200700184.

    Article  CAS  PubMed  Google Scholar 

  27. Saitoh H, Fujisawa S, Mitsuoka C, Ito A, Hirabuchi A, Ikeda K, Irieda H, Yoshino K, Yoshida K, Matsumura H: Large-scale gene disruption in Magnaporthe oryzae identifies MC69, a secreted protein required for infection by monocot and dicot fungal pathogens. PLoS Pathog. 2012, 8 (5): e1002711-10.1371/journal.ppat.1002711.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Klosterman SJ, Subbarao KV, Kang S, Veronese P, Gold SE, Thomma BPHJ, Chen Z, Henrissat B, Lee Y-H, Park J: Comparative genomics yields insights into niche adaptation of plant vascular wilt pathogens. PLoS Pathog. 2011, 7 (7): e1002137-10.1371/journal.ppat.1002137.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Gardiner DM, McDonald MC, Covarelli L, Solomon PS, Rusu AG, Marshall M, Kazan K, Chakraborty S, McDonald BA, Manners JM: Comparative pathogenomics reveals horizontally acquired novel virulence genes in fungi infecting cereal hosts. PLoS Pathog. 2012, 8 (9): e1002952-10.1371/journal.ppat.1002952.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.

    Article  CAS  PubMed  Google Scholar 

  31. Boutemy LS, King SR, Win J, Hughes RK, Clarke TA, Blumenschein TM, Kamoun S, Banfield MJ: Structures of Phytophthora RXLR effector proteins: a conserved but adaptable fold underpins functional diversity. J Biol Chem. 2011, 286 (41): 35834-35842. 10.1074/jbc.M111.262303.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Finn RD, Clements J, Eddy SR, HMMER web server: interactive sequence similarity searching: Nucleic Acids Res. 2011, 39 (Web Server issue): W29-W37.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Grigoriev IV, Nordberg H, Shabalov I, Aerts A, Cantor M, Goodstein D, Kuo A, Minovitsky S, Nikitin R, Ohm RA: The genome portal of the Department of Energy Joint Genome Institute. Nucleic Acids Res. 2012, 40 (Database issue): D26-D32.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Triplet T, Butler G: A review of genomic data warehousing systems. Brief Bioinform. 2013, doi: 10.1093/bib/bbt031

    Google Scholar 

  35. Hane JK, Lowe RG, Solomon PS, Tan KC, Schoch CL, Spatafora JW, Crous PW, Kodira C, Birren BW, Galagan JE: Dothideomycete plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum. Plant Cell. 2007, 19 (11): 3347-3368. 10.1105/tpc.107.052829.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Ma LJ, van der Does HC, Borkovich KA, Coleman JJ, Daboussi MJ, Di Pietro A, Dufresne M, Freitag M, Grabherr M, Henrissat B: Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature. 2010, 464 (7287): 367-373. 10.1038/nature08850.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Manning VA, Pandelova I, Dhillon B, Wilhelm LJ, Goodwin SB, Berlin AM, Figueroa M, Freitag M, Hane JK, Henrissat B: Comparative genomics of a plant-pathogenic fungus, Pyrenophora tritici-repentis, reveals transduplication and the impact of repeat elements on pathogenicity and population divergence. G3 (Bethesda). 2013, 3 (1): 41-63. 2013.

    Article  CAS  Google Scholar 

  38. Sarma GN, Manning VA, Ciuffetti LM, Karplus PA: Structure of Ptr ToxA: an RGD-containing host-selective toxin from Pyrenophora tritici-repentis. Plant Cell. 2005, 17 (11): 3190-3202. 10.1105/tpc.105.034918.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ: Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25 (9): 1189-1191. 10.1093/bioinformatics/btp033.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Houterman PM, Cornelissen BJ, Rep M: Suppression of plant resistance gene-based immunity by a fungal effector. PLoS Pathog. 2008, 4 (5): e1000061-10.1371/journal.ppat.1000061.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Rep M, van der Does HC, Meijer M, van Wijk R, Houterman PM, Dekker HL, de Koster CG, Cornelissen BJ: A small, cysteine-rich protein secreted by Fusarium oxysporum during colonization of xylem vessels is required for I-3-mediated resistance in tomato. Mol Microbiol. 2004, 53 (5): 1373-1383. 10.1111/j.1365-2958.2004.04177.x.

    Article  CAS  PubMed  Google Scholar 

  42. Houterman PM, Ma L, van Ooijen G, de Vroomen MJ, Cornelissen BJ, Takken FL, Rep M: The effector protein Avr2 of the xylem-colonizing fungus Fusarium oxysporum activates the tomato resistance protein I-2 intracellularly. Plant J. 2009, 58 (6): 970-978. 10.1111/j.1365-313X.2009.03838.x.

    Article  CAS  PubMed  Google Scholar 

  43. Cuomo CA, Guldener U, Xu JR, Trail F, Turgeon BG, Di Pietro A, Walton JD, Ma LJ, Baker SE, Rep M: The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science. 2007, 317 (5843): 1400-1402. 10.1126/science.1143708.

    Article  CAS  PubMed  Google Scholar 

  44. Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, Tetko I, Guldener U, Mannhaupt G, Munsterkotter M: The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004, 32 (18): 5539-5545. 10.1093/nar/gkh894.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Kelley LA, Sternberg MJ: Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc. 2009, 4 (3): 363-371. 10.1038/nprot.2009.2.

    Article  CAS  PubMed  Google Scholar 

  46. Infantes L, Otero LH, Beassoni PR, Boetsch C, Lisa AT, Domenech CE, Albert A: The structural domains of Pseudomonas aeruginosa phosphorylcholine phosphatase cooperate in substrate hydrolysis: 3D structure and enzymatic mechanism. J Mol Biol. 2012, 423 (4): 503-514. 10.1016/j.jmb.2012.07.024.

    Article  CAS  PubMed  Google Scholar 

  47. Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinforma. 2008, 9: 40-10.1186/1471-2105-9-40.

    Article  Google Scholar 

  48. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M: Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36 (Web Server issue): W465-W469.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. FigTree.http://tree.bio.ed.ac.uk/software/figtree/,

  51. Yennawar NH, Li LC, Dudzinski DM, Tabuchi A, Cosgrove DJ: Crystal structure and activities of EXPB1 (Zea m 1), a beta-expansin and group-1 pollen allergen from maize. Proc Natl Acad Sci USA. 2006, 103 (40): 14664-14671. 10.1073/pnas.0605979103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Li Y, Darley CP, Ongaro V, Fleming A, Schipper O, Baldauf SL, McQueen-Mason SJ: Plant expansins are a complex multigene family with an ancient evolutionary origin. Plant Physiol. 2002, 128 (3): 854-864. 10.1104/pp.010658.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Cosgrove DJ, Li LC, Cho HT, Hoffmann-Benning S, Moore RC, Blecker D: The growing world of expansins. Plant Cell Physiol. 2002, 43 (12): 1436-1444. 10.1093/pcp/pcf180.

    Article  CAS  PubMed  Google Scholar 

  54. Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 305 (3): 567-580. 10.1006/jmbi.2000.4315.

    Article  CAS  PubMed  Google Scholar 

  55. Kulkarni RD, Kelkar HS, Dean RA: An eight-cysteine-containing CFEM domain unique to a group of fungal membrane proteins. Trends Biochem Sci. 2003, 28 (3): 118-121. 10.1016/S0968-0004(03)00025-2.

    Article  CAS  PubMed  Google Scholar 

  56. Marshall R, Kombrink A, Motteram J, Loza-Reyes E, Lucas J, Hammond-Kosack KE, Thomma BP, Rudd JJ: Analysis of two in planta expressed LysM effector homologs from the fungus Mycosphaerella graminicola reveals novel functional properties and varying contributions to virulence on wheat. Plant Physiol. 2011, 156 (2): 756-769. 10.1104/pp.111.176347.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. de Jonge R, van Esse HP, Kombrink A, Shinya T, Desaki Y, Bours R, van der Krol S, Shibuya N, Joosten MH, Thomma BP: Conserved fungal LysM effector Ecp6 prevents chitin-triggered immunity in plants. Science. 2010, 329 (5994): 953-955. 10.1126/science.1190859.

    Article  CAS  PubMed  Google Scholar 

  58. Zhang XW, Jia LJ, Zhang Y, Jiang G, Li X, Zhang D, Tang WH: In planta stage-specific fungal gene profiling elucidates the molecular strategies of Fusarium graminearum growing inside wheat coleoptiles. Plant Cell. 2012, 24 (12): 5159-5176. 10.1105/tpc.112.105957.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Gale LR, Bryant JD, Calvo S, Giese H, Katan T, O’Donnell K, Suga H, Taga M, Usgaard TR, Ward TJ: Chromosome complement of the fungal plant pathogen Fusarium graminearum based on genetic and physical mapping and cytological observations. Genetics. 2005, 171 (3): 985-1001. 10.1534/genetics.105.044842.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Winnenburg R, Baldwin TK, Urban M, Rawlings C, Kohler J, Hammond-Kosack KE: PHI-base: a new database for pathogen host interactions. Nucleic Acids Res. 2006, 34 (Database issue): D459-D464.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Andrew M, Beacham J: A genomic fungal foray. Biologist. 2009, 56 (2): 98-

    Google Scholar 

  62. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19 (9): 1639-1645. 10.1101/gr.092759.109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8 (10): 785-786. 10.1038/nmeth.1701.

    Article  CAS  PubMed  Google Scholar 

  64. Spanu PD, Abbott JC, Amselem J, Burgis TA, Soanes DM, Stuber K, Ver Loren Van Themaat E, Brown JK, Butcher SA, Gurr SJ: Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism. Science. 2010, 330 (6010): 1543-1546. 10.1126/science.1194573.

    Article  CAS  PubMed  Google Scholar 

  65. Godfrey D, Bohlenius H, Pedersen C, Zhang Z, Emmersen J, Thordal-Christensen H: Powdery mildew fungal effector candidates share N-terminal Y/F/WxC-motif. BMC Genomics. 2010, 11: 317-10.1186/1471-2164-11-317.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Gabaldon T: Large-scale assignment of orthology: back to phylogenetics?. Genome Biol. 2008, 9 (10): 235-10.1186/gb-2008-9-10-235.

    Article  PubMed Central  PubMed  Google Scholar 

  67. Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434 (7036): 980-986. 10.1038/nature03449.

    Article  CAS  PubMed  Google Scholar 

  68. Kamper J, Kahmann R, Bolker M, Ma LJ, Brefort T, Saville BJ, Banuett F, Kronstad JW, Gold SE, Muller O: Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis. Nature. 2006, 444 (7115): 97-101. 10.1038/nature05248.

    Article  PubMed  Google Scholar 

  69. Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, Barry KW, Condon BJ, Copeland AC, Dhillon B, Glaser F: Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 2012, 8 (12): e1003037-10.1371/journal.ppat.1003037.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  70. Goodwin SB, M’Barek SB, Dhillon B, Wittenberg AH, Crane CF, Hane JK, Foster AJ, Van der Lee TA, Grimwood J, Aerts A: Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet. 2011, 7 (6): e1002070-10.1371/journal.pgen.1002070.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  71. Ellwood SR, Liu Z, Syme RA, Lai Z, Hane JK, Keiper F, Moffat CS, Oliver RP, Friesen TL: A first genome assembly of the barley fungal pathogen Pyrenophora teres f. teres. Genome Biol. 2010, 11 (11): R109-10.1186/gb-2010-11-11-r109.

    Article  PubMed Central  PubMed  Google Scholar 

  72. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K: WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007, 35 (Web Server issue): W585-W587.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Hobohm U, Sander C: A sequence property approach to searching protein databases. J Mol Biol. 1995, 251 (3): 390-399. 10.1006/jmbi.1995.0442.

    Article  CAS  PubMed  Google Scholar 

  74. Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.

    Article  CAS  PubMed  Google Scholar 

  75. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, 37 (Web Server issue): W202-W208.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

JS was supported by the CSIRO Transformational Biology Capability platform. DMG was partially supported by the Australian Grains Research and Development Corporation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jana Sperschneider.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JS developed and implemented the bioinformatics pipeline, performed data analysis and drafted the manuscript. All authors contributed to the design and development of the bioinformatics pipeline. DMG and JMM analysed the biological data and predictions. JMT performed statistical data analysis. JKH assisted in the development of the bioinformatics pipeline and the genomic analysis. JS, DMG, JMT, KBS and JMM conceived the study. All authors contributed to the draft and read and approved the final version.

Electronic supplementary material

12864_2013_7163_MOESM1_ESM.xlsx

Additional file 1:Table listing the properties of the predicted pathogen-associated proteins in S. nodorum .(XLSX 405 KB)

12864_2013_7163_MOESM2_ESM.xlsx

Additional file 2:Table listing the properties of the predicted pathogen-associated proteins in F. oxysporum .(XLSX 509 KB)

12864_2013_7163_MOESM3_ESM.xlsx

Additional file 3:Table listing the properties of the predicted pathogen-associated proteins in F. graminearum as well as the individual clusters and predicted motifs.(XLSX 822 KB)

12864_2013_7163_MOESM4_ESM.pdf

Additional file 4:Box plots for the 12 clusters in F. graminearum and the 35 protein features used for the clustering.(PDF 484 KB)

12864_2013_7163_MOESM5_ESM.pdf

Additional file 5:Contains the full sequence alignment for the proteins containing the motif [WYF]-C-x-T-Y-x-S-T-Y-L.(PDF 64 KB)

12864_2013_7163_MOESM6_ESM.xlsx

Additional file 6:Table listing the properties of the predicted pathogen-associated proteins in B. graminis .(XLSX 151 KB)

Additional file 7:Contain information about the predicted Y/F/WxC-motif proteins B. graminis .(ZIP 2 KB)

Additional file 8:Contain information about the predicted Y/F/WxC-motif proteins B. graminis .(PDF 11 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sperschneider, J., Gardiner, D.M., Taylor, J.M. et al. A comparative hidden Markov model analysis pipeline identifies proteins characteristic of cereal-infecting fungi. BMC Genomics 14, 807 (2013). https://doi.org/10.1186/1471-2164-14-807

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-14-807

Keywords