Gene discovery in EST sequences from the wheat leaf rust fungus Puccinia triticina sexual spores, asexual spores and haustoria, compared to other rust and corn smut fungi
- Junhuan Xu†1, 10,
- Rob Linning†1,
- John Fellers2,
- Matthew Dickinson3,
- Wenhan Zhu4,
- Ivan Antonov4,
- David L Joly1,
- Michael E Donaldson5,
- Tamar Eilam6,
- Yehoshua Anikster6,
- Travis Banks7,
- Sarah Munro8,
- Michael Mayo8,
- Brian Wynhoven8,
- Johar Ali8,
- Richard Moore8,
- Brent McCallum7,
- Mark Borodovsky4,
- Barry Saville9 and
- Guus Bakkeren1Email author
© Xu et al; licensee BioMed Central Ltd. 2011
Received: 6 October 2010
Accepted: 24 March 2011
Published: 24 March 2011
Rust fungi are biotrophic basidiomycete plant pathogens that cause major diseases on plants and trees world-wide, affecting agriculture and forestry. Their biotrophic nature precludes many established molecular genetic manipulations and lines of research. The generation of genomic resources for these microbes is leading to novel insights into biology such as interactions with the hosts and guiding directions for breakthrough research in plant pathology.
To support gene discovery and gene model verification in the genome of the wheat leaf rust fungus, Puccinia triticina (Pt), we have generated Expressed Sequence Tags (ESTs) by sampling several life cycle stages. We focused on several spore stages and isolated haustorial structures from infected wheat, generating 17,684 ESTs. We produced sequences from both the sexual (pycniospores, aeciospores and teliospores) and asexual (germinated urediniospores) stages of the life cycle. From pycniospores and aeciospores, produced by infecting the alternate host, meadow rue (Thalictrum speciosissimum), 4,869 and 1,292 reads were generated, respectively. We generated 3,703 ESTs from teliospores produced on the senescent primary wheat host. Finally, we generated 6,817 reads from haustoria isolated from infected wheat as well as 1,003 sequences from germinated urediniospores. Along with 25,558 previously generated ESTs, we compiled a database of 13,328 non-redundant sequences (4,506 singlets and 8,822 contigs). Fungal genes were predicted using the EST version of the self-training GeneMarkS algorithm. To refine the EST database, we compared EST sequences by BLASTN to a set of 454 pyrosequencing-generated contigs and Sanger BAC-end sequences derived both from the Pt genome, and to ESTs and genome reads from wheat. A collection of 6,308 fungal genes was identified and compared to sequences of the cereal rusts, Puccinia graminis f. sp. tritici (Pgt) and stripe rust, P. striiformis f. sp. tritici (Pst), and poplar leaf rust Melampsora species, and the corn smut fungus, Ustilago maydis (Um). While extensive homologies were found, many genes appeared novel and species-specific; over 40% of genes did not match any known sequence in existing databases. Focusing on spore stages, direct comparison to Um identified potential functional homologs, possibly allowing heterologous functional analysis in that model fungus. Many potentially secreted protein genes were identified by similarity searches against genes and proteins of Pgt and Melampsora spp., revealing apparent orthologs.
The current set of Pt unigenes contributes to gene discovery in this major cereal pathogen and will be invaluable for gene model verification in the genome sequence.
Puccinia triticina (Pt) has a complex life cycle which includes five different spore types and two hosts: wheat (Triticum aestivum L.) and meadow rue (Thalictrum speciosissimum L.). The latter plant is the so-called alternate host on which the fungus completes its sexual stage [1–3]. Sex is not essential and infection of and spread on wheat through re-infection constitutes the asexual cycle. The brown-coloured urediniospores from which the rust got its name, are the asexual infectious propagules. They are easily carried long distances by prevailing winds and can lead to epidemics. Early processes in infection include urediniospore attachment, germination and the formation of a germtube. The tips of germtubes differentiate into appressoria which develop over stomatal lips and entry into the substomatal cavity is gained forcibly by turgor pressure. Within 24 hrs after spore germination, a haustorial mother cell is formed adjacent to a plant cell within the cavity and cell wall penetration takes place. Subsequent invagination of the host plasmalemma results in the first intimate contact. Thereafter, a microscopically visible haustorial interface surrounding the mature feeding structure, is produced, likely made up of both fungal and host material . This interface is critical in governing protein and metabolite traffic [5, 6] since haustoria are thought to secrete a suite of proteins, some of which are aimed at suppressing host defence responses that may be triggered by the fungus when it penetrates the plant cell wall or at establishing the feeding interaction. In compatible interactions, the fungus colonizes the plant and within 7 days can produce uredinia (pustules) containing asexual urediniospores which are released and give rise to new rounds of infection.
On senescing wheat plants, uredinia respond to cues and switch to producing black teliospores. These are survival propagules with a complex, multi-layered wall and no vacuoles. They contain lipid droplets and glycogen-like material [7, 8]. Teliospores are primarily 2-celled with each cell containing two haploid nuclei that have paired, if not fused to form the diploid state . They often appear in low numbers on the lower parts of the plant, including the stem. Under the proper conditions, the teliospore germinates and a metabasidium forms, generally from both cells, in which meiosis occurs and on which four haploid basidiospores develop. A third mitotic division occurs resulting in basidiospores having two nuclei though being monokaryons, containing two nuclei of the same type . These basidiospores are ephemeral and can be dispersed to infect the alternate host. This infection occurs via direct penetration of the plant epidermis and the subsequent production of monokaryotic (M-) haustoria. These structures are extensions of the intercellular hyphae which penetrate the plant cells in a manner reminiscent of that described for the corn smut fungus, Ustilago maydis. Neither infection shows the morphological specialization seen for the P. triticina dikaryotic (D-) haustorium . Upon establishing a feeding relationship and presumably suppressing alternate host defence responses, the monokaryotic hyphae yield specialized pycnia which generate pycniospores embedded in nectar. Because they originated from haploid meiotic products, the pycniospores represent different mating types and can cross-fertilize, often through the action of insects attracted by the nectar. After fertilization, that is, the fusion of one pycniospore to a receptive hypha in the pycnium of a different mating type followed by nuclear transfer, the newly formed dikaryon undergoes developmental reprogramming. The resultant mycelium traverses the leaf and forms aecia on the underside in which dikaryotic aeciospores develop . Aeciospores are dispersal propagules which will infect the primary wheat host.
The generation of genomic resources for the cereal rust fungi is gearing up. A draft genome sequence for the related stem rust fungus, P. graminis f. sp. tritici (Pgt) was released in 2007, a rough draft for Pt was released in 2009, and genome sequences are being generated for the related wheat stripe rust fungus, P. striiformis f. sp. tritici (Pst; http://www.broadinstitute.org/annotation/genome/puccinia_group/MultiHome.html). The generation of Expressed Sequence Tags (ESTs) is essential for proper gene prediction in genomes. Previous Pt EST libraries contributed to gene discovery and stage-specific expression analyses [12–14]. Similar EST collections were generated from other Puccinia species [15–21]. However, all these studies were focused on urediniospores and the wheat infection cycle. There is no molecular data on genes involved in the sexual stages in the rusts. Unique to this study, we surveyed three other spore types representing the sexual stage: teliospores, pycniospores and aeciospores. We found that each spore type yielded rather specific EST sequences and that pycniospores and teliospores in particular seemed to express a unique set, when compared to all other sampled stages.
In addition, we generated and analysed ESTs from two sets of isolated haustoria and from the germinated urediniospore stage. Haustoria ESTs included a large set of unique sequences of unknown function. In total, over 17,000 new ESTs, of which many appeared spore-specific, were added to our existing collection of 25,558 sequences and compiled into a single unigene set. Subsequently, fungal genes were predicted using bioinformatic algorithms trained on Pgt sequences, and were then compared to Pgt and model basidiomycete plant pathogen, U. maydis genomic resources. Notably, comparison to U. maydis revealed interesting leads for research into spore biology.
Results and Discussion
Life cycle stage-specific cDNA libraries
Three different sexual stage spore types - Minute amounts of nectar, containing pycniospores of different mating types, as well as dikaryotic aeciospores, were collected from Thalictrum pycnia and aecia residing on the upper and abaxial side of leaves, respectively. Both life cycle stages were from the same isolate, but different from race 1. Race 1 teliospores were induced on senescing wheat and small amounts of dormant spores cleanly dissected from telial pustules. From these three stages, only nanogram-amounts of total RNA could be obtained. From the aeciospore material, cDNA library PT029 was constructed and the pycniospore RNA resulted in library PT030. The RNA from the teliospores produced cDNA library PT031.
Isolated haustoria -The haustorium is an important structure involved in bidirectional traffic of proteins and uptake of nutrients and ESTs might reveal expression of genes involved in host defence suppression and fungal feeding process. Moreover, haustoria are expected to produce effectors, some of which might have avirulence functions as has been shown to occur in other systems [22–24]. We therefore isolated haustoria from heavily-infected wheat leaves and produced two cDNA libraries from two different isolates: library designation PT033 (isolate PRTUS3) and PTDH (isolate WBRS-97-3; see Methods).
Germtube - Rust isolate WBRS-97-3 was also used to produce a cDNA library from urediniospores germinated over water, library designation PTDG, from which bi-directional reads were generated.
Compilation of EST sequences and Pt gene finding
cDNA libraries produced and number of generated quality EST reads
predicted fungal c
Several of the cDNA libraries were generated from RNA isolated from infected plants and would contain host cDNA; other (bacterial) contaminants are also fairly abundant when such approaches are used. Therefore, we predicted the putative fungal genes using two approaches, ab initio and comparative genomics; we compared the results and now present a consolidated prediction set of 6,308 genes.
Ab initio gene prediction approach using the EST version of GeneMarkS
Second, to further improve estimates of the HMM parameters, we proceeded with validation of the automatically derived training set. We used BLASTP to search for homologs of predicted Pt proteins among 20,567 proteins encoded in the genome of Pgt. At a cut-off e-value of 1e-10, this procedure identified a set of 2,093 Pt genes with Pgt homologs. This reduced set was used as a refined Pt training set to derive parameters for the 4th order Markov model of a protein-coding region, a critically important submodel of the gHMM.
The 4th order model was more specific; it omitted 3,579 genes identified by the 2nd order model. A comparative G + C content analysis of the genes predicted by the 2nd and 4th order models showed that the genes predicted by the 4th order model had a slightly higher G + C content than those predicted by the 2nd order model (51% versus 48%). Notably, a number of genes now absent from the predictions could be attributed to contaminant sequences, possibly from the wheat host (see Contaminants section below).
Also, we observed that the use of the 4th order model led to a significant increase, by 1,503, in the number of fragments with no genes predicted. A total of 5,350 unigenes out of 10,576 analyzed unigenes longer than 300 nt, were identified as containing no fungal genes. Still, these 5,350 unigenes did contain protein-coding sequences. When we switched gene prediction model in GeneMark.hmm to a generic 2nd order heuristic model specific for genomic DNA with a given G + C content , we identified 5,156 protein coding regions (Additional file 1). However, at e < 10-5, BLASTP searches revealed similarity to Pgt proteins for only 14 out of these 5,156 predicted proteins and 19 had similarity to proteins in the NCBI nr database. Of the 19 "non-Pgt" proteins, at least 9 had similarity to proteins in plants, human and E. coli thus indicating possible contaminants (Additional file 2). The majority of the 5,156 predicted protein coding regions (3,501 or 68%) had lengths of less than 150 nt. It is likely that these EST sequences represent 5', or 3'- UTRs and predicted short coding regions are mere artifacts or represent very small proteins.
Comparative genomics approach using BLASTN
As part of a Pt genome sequencing project (Cuomo, Fellers, Szabo and Bakkeren), an initial set of short genome sequences generated through '454 pyrosequencing' (see Methods) was available. The 13,328 EST unigene sequences were compared by BLASTN to this genome collection and 4,749 yielded significant matches at e < 10-5. This number differed from the Pt unigene set of 5,007 that we arrived at upon using GeneMarkS; a BLASTN comparison showed the two sets had 3,448 unigenes in common. At this early stage of Pt genome analysis it is uncertain how much coverage has been obtained with the 454-genomic reads. In addition, 13 out of the 14 unigene sequences not called by GeneMark.hmm but matching Pgt proteins by BLAST (see above), had BLASTN hits to the Pt genomic reads. We therefore assembled a final set of 6,308 putative Pt unigenes composed of the union of the two predicted sets. Of these, 3,383 (54%) were contributed through this new study and constitute over 85% of the total 3,990 newly-added unigene sequences. The percentage of called fungal sequences was even higher for individual spore stages (89 - 95%; Table 1) indicating the expected purity of collected spore samples.
Using the BLASTN algorithm, out of the total set of 13,328 EST unigenes, 6,095 did not match any sequences in the surveyed databases although 1,452 still matched the preliminary Pt genome reads in a BLASTN search and 1,656 the more-inclusive predicted fungal unigene set. Whether these sequences represent Pt-specific gene sequences or transcribed genome regions where no genes have (yet) been predicted, remains to be investigated. Although such sequences do not currently contribute to gene discovery, they will be valuable for the verification of gene models or transcribed regions in the completed genome sequence.
Although one would expect to see just one continuous protein-coding region per unigene, in the output of the EST version of GeneMarkS with the Pt-specific 4th order model, we could identify up to nine predicted genes per unigene and sometimes no genes (Figure 2). Multiple predictions could be explained by the presence of sequencing errors. To identify possible frameshifts caused by sequencing errors or chimeric cDNA inserts, a new frameshift prediction program GeneTack  was applied to 3,851 ESTs with several genes predicted in the same strand. Separate genes with no detected frameshifts were identified in 2,645 ESTs (69%) whereas a single frameshift was predicted in 905 ESTs (24%). Protein-coding regions with several frameshifts were identified in the remaining 301 ESTs (Additional file 3). Coding sequences with a single predicted frameshift were conceptually translated. This approach produced 805 sequences used to validate the GeneTack frameshift predictions by a BLASTP search against the Pgt proteome (20,567 sequences) and the NCBI non-redundant database (8,924,078 sequences). At e < 10-5, database hits overlapping the sites of predicted frameshifts confirmed 325 frameshift predictions; in most cases the frameshift was caused by misinterpretation of single-pass reads. When inspection of sequencing electropherograms allowed, the corresponding EST sequences were corrected.
Several cDNA libraries were constructed from host infection stages. In order to not discard any potential fungal sequences from our unigene set, we employed a rather conservative strategy which most likely left in the database a number of "contaminating" wheat host sequences. To further understand the discrepancy between the total set of 13,328 EST unigenes and the 6,308 predicted fungal unigenes, all unigenes were also compared by BLASTN to a large collection of available wheat EST sequences to identify possible host contaminants. Some 3,080 unigenes matched wheat ESTs at e < 10-5 and just over 3,300 had similarity to a 1x wheat genome coverage of '454-generated' sequences (http://www.cerealsdb.uk.net/index.htm). Still, 42% of unigenes matching wheat sequences (1,392 unigenes) also matched the preliminary Pt genomic reads. Since e-values are related to the size of the database, care has to be taken when directly comparing values obtained from databases of significantly different sizes. Nevertheless, in a smaller EST unigene set , we had also observed unigenes that matched fungal and wheat sequences with e-values in a similar range. Whether possibly horizontal gene transfer during co-evolution is involved might be solved once the genomes of both organisms are available.
Some 7,020 sequences were not included in the 'putative Pt fungal' set; they did not match the preliminary, still partial Pt genome reads. Of these 63% (4,439 unigenes) did not match any sequences in the databases surveyed. However, 26 sequences matched Pgt proteins at e < 10-5, 2% matched sequences in the nr database, 11% matched wheat ESTs, 22% matched fungal and oomycete ESTs in COGEME , and 30% matched sequences in EST-other. Those sequences without database matches (4,439 unigenes or 63%) were mainly derived from one appressorial library ('pta': 1,727 unigenes), from one specific germinated urediniospore stage library PT002 ('ptg': 1,598 unigenes), and from a wheat-infection stage library library PT009 ('pth': 520 unigenes; ). What these three cDNA libraries had in common was that cDNA inserts were generated by PCR which can reduce tag size. When the size distribution of these 4,439 unigene sequences was evaluated, almost 45% had lengths of less than 500 bp (compared to less than 29% among all 13,328 unigenes). When a comprehensive Pt genome is available, these sequences will be re-evaluated.
Annotation and classification of ESTs
The 13,328 Pt EST unigene sequences and the predicted fungal subset were compared by TBLASTX or BLASTX to various databases. Using a cut-off value of e < 10-5, approximately 40% of the putative fungal subset of 6,308 unigenes matched NCBI's nr, with a similar number matching the UniProt database. The corresponding percentages were 27% for the fungal and oomycete-specific EST database COGEME and almost 60% for dbEST (Additional file 4). When the total 13,328-member unigene set was compared to dbEST, roughly a third more (2,147 unigenes for a total of 5,854 or 44%) matched, but most of these (1,628 unigenes) scored with low significance at e < 10-10.
Preliminary annotation was achieved for 2,556 unigenes by choosing the most informative header arising from the various database BLAST searches (at e < 10-5) which was not necessarily the most statistically significant return; 2,286 (90%) were judged to be of fungal origin and 48 had similarity to ribosomal sequences of various rusts (Additional file 5). It is notoriously difficult to assign genes to categories since the proteins they encode can have multiple functions and/or be associated with various aspects of their biology. However, to achieve a better understanding of their role(s), we analysed all unigenes with the program 'annot8r'  which assigned various GO (Gene Ontology), EC (Enzyme Commission number) and/or KEGG (Kyoto Encyclopedia of Genes and Genomes) annotations to almost 1,500 unigenes.
Similarity searches to related fungi
Similarity searches of Pt unigene sequences to related fungal databases
Pgt genome a
Pgt proteins b
Pt genome c
Pst ESTs d
Mlp proteins e
U. maydis ESTs f
U. maydis genome g
< = -100
-100 < × ≤ -50
-50 < × ≤ -20
-20 < × ≤ -5
The Pt unigene set was also compared to a set of 3,297 published, putative wheat stripe rust (Pst) ESTs from comparable life cycle stages (urediniospores and haustoria, Table 2). At e < 10-5, roughly 20% of the Pt fungal unigenes matched this Pst collection, whereas almost half matched the Pgt resources (see above) which themselves have an estimated genome coverage of approximately 92%. The current number of gene calls for Pgt is 20,567, 18,241 when disregarding many transposable element- (TE) related calls, and will likely be similar for Pt (and Pst). We have only compared a partial Pt gene set, based on the EST unigenes, but, assuming we can extrapolate the numbers, these surprisingly low percentages suggest substantial species-specific gene complements and could indicate considerable genome evolution. Indeed, even though these rust fungi are considered to belong to the same genus, they fall into well-separated phylogenetic sub-clades . Gene divergence could occur due to strong selection of these biotrophs. By comparison, two related biotrophic smut fungi from different genera, U. maydis and Sporisorium reilianum, still share an average 74% amino acid identity among predicted proteins . A closer look at the set of unigenes that matched Pt genome sequences but not Pgt resources (2,157 unigenes or 34%), revealed that 725 unigenes had representation in the 'in planta infection' and 'isolated haustoria' cDNA libraries (out of a pool of 1,344 generated unigenes matching Pt genomic sequences), 719 in urediniospores (out of 1,627), 374 in pycniospores (out of 1,101), 247 in teliospores (out of 790) and less than 100 in other stages. As a percentage of their respective pools (54%, 44%, 34% and 31%), this suggests that genes expressed in these four life cycle stages might be under different selection pressures. Using this rough estimate, the highest discrepancy between Pt and Pgt is seen among genes expressed during plant infection and it could be argued that many of those represent virulence factors likely under selection pressure as has been seen among poplar rust, Melampsora species .
A particular focus of this study was to compare the Pt unigenes to U. maydis genomic resources since this fungus is the best-studied, closest-related model representing cereal-biotrophic basidiomycete interactions. We have previously illustrated the possibility of using the U. maydis model system for performing functional analyses of Pt genes . A BLASTX search against the 6,846 predicted U. maydis proteins  yielded 1,535 matches at e < 10-5, 907 (13% of the total Um complement) with significant similarity at e < 10-20 (Additional file 6). Most matching Pt unigenes (1,506) were predicted 'fungal' as might be expected from genes having orthologs in basidiomycetes belonging to different genera (Table 2). Indeed, a large majority of 1,416 Pt-Um "pairs" (92%) also matched Pgt proteins suggesting the existence of a common set of proteins/functions among these basidiomycetes. However, 119 Pt-Um "pairs" (8%) did not match Pgt proteins. In a complementary approach, the complete set of 20,567 predicted Pgt proteins was compared to the same set of Um proteins, yielding 5,624 matches of which 4,012 (59% of the total Um complement and 20% of the Pgt complement) revealed significant similarities at e < 10-20. Of these Pgt-Um "pairs", 2,604 (46%) also matched Pt unigenes; this smaller percentage likely reflects the large comprehensive Pgt genome-predicted protein set compared to the Pt EST subset. The identification of many homologs and the existence of many Um gene deletion mutants or the ease with which Um gene deletions can be obtained , could allow for functional Pt gene analyses through complementation in U. maydis.
Several unigenes revealed interesting putative functions, in particular among the ESTs from the pycniospore library, a stage not yet covered in the literature. PtContig7448 matched a pheromone transporter from Cryptococcus neoformans, a human opportunistic basidiomycete pathogen, and was covered by sequences from 4 different cDNA clones from only that cell type. Contig7830 (one cDNA clone) matched pheromone receptor sequences in many fungi (e.g., the A2 pheromone receptor from Microbotryum violaceum, the anther smut fungus) and was similar to a 357-amino acid protein in Pgt. Pt and related rusts are thought to possess bipolar mating systems (+ or -)  although some might possess multiple-allelomorphic tetrapolar systems similar to the higher mushrooms . Because pycniospores are haploid cells in which mating type has segregated and are primed for mating to form the dikaryotic aeciospores, they are expected to express mating-type specific genes. In rusts, no information on mating-type genes has been published. It would be tremendously valuable to uncover the molecular basis of the mating system and reveal the structure of MAT loci in cereal rusts. The revealed Pt pheromone receptor gene, pra, could represent a component of the mating-type locus. Contig7901 (one cDNA clone) has similarity to a protein belonging to the oligopeptide OPT transporter family found in many fungi and annotated as 'sexual differentiation process protein isp4'. The pycniospore stage also revealed several genes involved in sterol biosynthesis with possible roles in mating (Contigs 7279, 7412, 7754 and 7777). In S. cerevisiae, (ergo-)sterols have been shown to be involved in cell fusion during mating, for example, sphingolipid and ergosterol biosynthetic mutants fail to polarize proteins to the tip of the deforming yeast cell ("shmoo") and are therefore deficient in mating [39, 40]. However, another ergosterol biosynthesis gene, Erg28 (Contig7716) and an oxysterol-binding protein (PT0318.P01.C21.ptt), were found in the teliospores library, whereas a sterol O-acyltransferase (Contig107) had one representative EST in urediniospores and one in germinated urediniospores. Among the whole unigene set, a few Contigs, 7490 (2 cDNA clones), 7853 (9 cDNA clones) and 7905 (2 cDNA clones), were obtained only from the pycniospore library and were most similar to Pgt predicted protein PGTG_01034.2, a carbon catabolite-derepressing protein kinase. A comparative multiple sequence alignment of DNA sequences in AlignX confirmed all three Pt unigenes were closely related but had four approximately 80 bp-stretches of DNA that were shared in different combinations. Whether these represent unique genes or are different splice products will be solved once the Pt genome is available; however, the unique stretches did not match predicted introns in various Pgt homologs of which there were many (at least 3 at e = 0). We are interested in developmental-specific kinases and intriguingly, these Pt contigs also bore great resemblance to the catalytic subunit of cAMP-dependent protein kinase, known to be involved in responsiveness to mating pheromone in several fungi [41, 42].
Teliospores are survival propagules and are melanized for protection. Contig7966 (one teliospores-specific cDNA clone), matched a diphenol/urishiol (multi-copper) oxidase or laccase (EC:22.214.171.124), an enzyme found in fungi and plants which oxidizes different phenols and diamines and is implicated in the production of melanin pigment.
Several unigenes had EST representation mainly in all spore stages: Contig5605 (TAR1, regulation of mitochondrial gene expression), Contig8158 (a P450 monooxygenase) and Contigs 279, 5253 and 7316 (no hits). It is possible that certain genes are involved in general spore-forming processes despite their very different morphology and physiology.
Correlating elevated EST levels of spore-specific libraries between P. triticina and U. maydis
Spore production is critical to the dispersal of rusts and smuts; therefore, knowledge of the development and biology of spores is essential to understanding disease spread. Since spore development in these fungi requires growth in the host plant, it is reasonable to assume that their development is triggered in response to a signal(s) received from the host. In this context, identifying genes with elevated transcript levels in spores of both fungi may lead to the discovery of common responses to plant signals. Of the four Pt rust spore stages we generated ESTs from in this study, only teliospores have a biological equivalent in the life cycle of Um (represented by ESTs from the 'dormant teliospore' cDNA library TDO ). Comparing expression among libraries of these two species is complicated by variation in library preparation methods, notably the requirement to include an amplification step to clone cDNAs from the low levels of Pt teliospores, and the need to normalize libraries being compared to correct for EST pool size. With these limitations, we expected to find few commonly expressed genes. Additional file 6 shows that Pt contig5874, a putative mitochondrial inner membrane protein involved in protein import, has elevated numbers of ESTs in both teliospore derived libraries. Elevated numbers of ESTs in Pt contig7329, a NADH-ubiquinone oxidoreductase, originated from Pt pycniospore and teliospore libraries. The similar um10989 (e = 10-34) has elevated EST numbers in Um teliospore and haploid cell libraries (Additional file 6). Pt pycniospores and Um haploid cells are both descendents of basidiospores, the direct products of meiosis. The discovery of these genes is interesting and they deserve further investigation. However, the fact that we found any commonly expressed genes given the limitations in the analysis is perhaps more important and supports a deeper analysis of the teliospore transcriptome in these fungi.
Comparing libraries from other spore (cell) types that are functionally equivalent between Pt and Um revealed further similarities in elevated EST levels. Four Pt contigs with elevated numbers in the pycniospore had homologs with elevated numbers in Um haploid cells or germinating teliospores. Two represent conserved hypothetical proteins and the others represent cofilin, an actin binding and severing protein, and a 26S proteasome non-ATPase regulatory subunit (Additional file 8). Six Pt contigs having elevated EST counts in the germinated urediniospore stage had Um homologs with higher number of ESTs, relative to the other stages, in the functionally equivalent Um dikaryon (Additional file 8). It is notable that four of these were conserved hypothetical proteins. These may represent functionally conserved genes. The function of these genes can be investigated in the tractable Um. Further expression analysis is required to confirm these results, but finding similarly high EST counts between similar cell types supports the benefit of continuing these comparative approaches in more depth using RNA sequencing techniques.
Stage-specific Pt unigene counts and representation in Pgt
# unique unigenes a
predicted fungal b
# matching Pgt homologs d
The combined 13,328 unigene set was queried using InterProScan to identify protein sequences predicted to be secreted. To complement this, 11,638 proteins predicted from the partial Pt genome were queried using the method by Joly et al.  to yield 758 potentially secreted Pt proteins. Reciprocal BLAST searches were performed to generate a list of Pt-specific, potentially secreted proteins. This list was complemented by searches against a subset of 1,699 potentially secreted Pgt proteins derived similarly from the 20,567 genome-predicted Pgt proteins, and additionally against collected sets of various predicted "secretome" proteins: 689 assembled EST sequences from various poplar leaf rust fungi, Melampsora spp. (S + category in ), 28 identified flax rust, Melampsora lini haustoria-specific secreted proteins , 100 bean rust, Uromyces fabae, secretome proteins , and 386 described secreted proteins in U. maydis. Additional file 9 shows which Pt unigenes are predicted to be secreted or are similar by BLAST searches to predicted secreted proteins in other (related) rust or basidiomycete plant biotrophic fungi. The representation of the ESTs among the various cDNA libraries is indicated, representing expression in the respective life cycle stages. Striking were the large groupings revealed; a given Pt unigene would match various different sequences, up to 27 Pgt proteins in the case of Contig75, with various e-values. These are likely family members of related secreted proteins, some of which seem to belong to clusters (such as PGTG_8431 to PGTG_8436, and PGTG_04883 to PGTG04887). Although many such groups were found matching Pgt, large poplar rust groups were also evident, and often a Pt unigene would match similar large groups in both organisms; sometimes they would match uniquely in poplar rust. A TBLASTX search of 583 Pt unigenes predicted to be secreted (based on the searches mentioned above) against the complete Pt unigene set ('self-BLAST') yielded related Pt groups which in general correlated well with the Pgt and poplar rust groups (Additional file 9, 'Pt secr prots self-BLAST' tab).
Overall, 2,421 predicted fungal unigenes, of which 1,344 matched Pt genomic 454 reads, had at least one EST generated from plant-infected material, represented by the stages pth, pti or ptih. Of the 1,402 unigenes representing isolated haustoria, 337 did not match any sequence in various databases including Pgt proteins, but of these, still 196 were predicted to be fungal and 157 matched Pt genomic reads. These could represent species-specific genes that are potentially preferentially expressed in haustoria, among which could be effectors. Transcript profiling or gene-specific quantitative PCR will be needed to verify whether these genes are expressed uniquely or at elevated levels during plant infection and hence might play important roles.
The Pt EST sequences generated, collected and analysed in the current study, enriched by sampling spore stages not covered in other studies, will provide essential information for proper gene calling and annotation of the Pt genome currently being assembled. For the first time, all spore stages from a macrocyclic cereal rust fungus, except for basidiospores, have been investigated for expressed genes through this EST approach. Many obtained unigene sequences from the various spore stages seemed specific suggesting radically different developmental gene sets are involved. Unfortunately, many could not be annotated due to lack of similar sequences in the various databases. Comparing potential spore-specific gene sequences between Pt and the model basidiomycete plant biotroph U. maydis, seems feasible and has yielded some interesting candidates for follow-up studies.
Although substantial homology exists between the Pt unigene sequences and near-comprehensive genomic resources from the related stem rust fungus, P. graminis f. sp. tritici, many genes seem to have diverged substantially from those from the putative common ancestor and many others appear to be species specific. A preliminary comparative sequence analysis of several BAC clones harbouring large Pt genomic inserts with the syntenic regions in Pgt, seems to corroborate this (JF and GB, unpublished). Future comprehensive comparisons once the genomes of several cereal rust species are available will shed light on the molecular evolution of their genes and genomes.
The generated ESTs have also been invaluable in generating molecular markers. The current gene sequences have been used to reveal SSRs (Simple Sequence Repeats or microsatellites) which were used in genetic diversity studies , and the variation among compiled contigs, potentially being derived from two parental nuclei in case of the dikaryotic material or the 4 different isolates used, revealed SNPs (Single Nucleotide Polymorphisms) which are being used to saturate a Pt genetic map (unpublished).
The computational prediction of a comprehensive secretome for Pt is one of the goals once the genome sequence is completed. Many secreted proteins will function in general fungal development such as cell wall biogenesis, but a significant subset will be crucial to support the biotrophic life style and hence will be considered virulence factors. Some will function in generating and maintaining the matrix surrounding haustoria and possibly hyphae roaming the plant tissue and contribute to feeding, but others, effectors among them, will be needed to subdue the host defence machinery. It is anticipated that many of the effectors will represent the many avirulence gene products that trigger resistance in the numerous wheat cultivars. Unravelling this complement is crucial in revealing old and new interactions which will aid in breeding durable resistance in our cereal crops.
Fungal strains and plant materials
To produce teliospores, Pt race1 (BBBD; ) was inoculated on susceptible wheat (Triticum aestivum L.) cultivar Thatcher and the infection was allowed to proceed for months until senescence of the plants was induced by reducing watering and temperature. Dormant teliospores were then dissected from senescent wheat stalks with a fine scalpel; this material was used to generate cDNA library PT031. Aeciospores and pycniospores were collected from the alternate host, meadow rue or Thalictrum speciosissimum, infected with germinating and basidiospores-producing telial material (isolate 99193) in the greenhouse (Y. Anikster and T. Eilam). Pycnial nectar with spores was collected and dried in a lyophilizer for 24 h; this material was used to construct cDNA library PT030. Aeciospores were dissected from infected T. speciosissimum leaves resulting in cDNA library PT029. Haustoria were isolated from susceptible wheat, cv. Wichita, heavily-infected with Pt urediniospores from isolate PRTUS6 (race PBJL; ), using conA-loaded columns essentially as described ; this material generated cDNA library PT033 (J. Fellers). The latter method was also employed to isolate haustoria from wheat cv. Vuka heavily infected with isolate WBRS-97-3 (material for cDNA library PTDH, M. Dickinson). The same isolate was also used to generate urediniospores which were then germinated over water for 8 hours and then collected as described  yielding material for cDNA library PTDG (M. Dickinson).
Isolation of total RNA and cDNA library construction
Dissected teliospores, haustorial preparations and germinated urediniospores were ground to a fine powder in liquid nitrogen using a mortar and pestle. The minute amounts of aeciospore and pycniospore material were ground in Eppendorf tubes with sterile sand. Total RNA was extracted in a solution of phenol-guanidine isothiocyanate (TRI Reagent, Molecular Research Center, Inc., Cincinnati, OH) according to instructions provided by the manufacturer. Due to the very small sample sizes, only 0.05 μg of total RNA could be obtained from the pycnio-, and aeciospores, and 1 μg from the teliospores. From these materials, cDNA libraries were constructed from total RNA using the SMART cDNA library kit (Clontech) according to the manufacturer's instructions with minor changes. The cycle number of LD PCR for cDNA synthesis was adapted; 26 cycles for pycnio- and teliospores, and 2 times 26 cycles (using an aliquot of the first amplification for the next 26 cycles) for aeciospores. The cDNA PCR products were digested with SfiI, then with proteinase K, and subsequently size-fractionated on 1% agarose gels in 0.5xTBE (45 mM Tris, pH 7.0, 45 mM boric acid, 1 mM EDTA). Fragments larger than 600 bp were excised and purified with the QiaQuick Gel Extraction Kit (Qiagen, Mississauga, ON). DNA fractions were ligated directionally into vector pDNR-LIB and transformed into E. coli DH10B/r. Total RNA from haustorial and germinated urediniospore material was converted to cDNA using oligo-dT primers containing a Xho1 site and Reverse Transcriptase (ZAP cDNA synthesis kit from Strategene which contains StrataScript™ Reverse Transcriptase). After adding EcoRI linkers to blunted 5'-ends, cDNA was cloned into vector lambda-zap; conversion of the library resulted in cDNA clones in pBluescript SK- in E. coli DH10B/r. RNA isolation and construction of cDNA library PT033 (isolated haustoria) was according to Li et al. .
DNA sequencing of the cDNA clones from the teliospore -, aeciospore -, pycniospore - and isolated haustoria (PT033) stages was performed at the Michael Smith Genome Sciences Centre (Vancouver, BC, Canada). Clones were picked on a Genetix Q-Pix into 384-well plates containing 80 μL 2× YT medium with 7.5% glycerol and plasmids were isolated from 96-well plate cultures using a modified alkaline lysis method. Inserts representing the 3'-end of the mRNA were sequenced using the standard M13 forward primer (F: 5'-GTAAAACGACGGCCAG, or B21/C21: 5'-TGTAAAACGACGGCCAGT), the T7 promoter primer (B7/C7: 5'-AATACGACTCACTATAG) or primer TB (5'-(T)23 A/C/G) used to overcome stuttering and slippage during sequencing due to long polyT stretches. Inserts representing the 5'-end of the mRNA were sequenced using M13 reverse primers (BR/CR: 5'-CAGGAAACAGCTATGAC, or R/CPTR: 5'-AACAGCTATGACCATG). The BigDye Terminator Cycle Sequence chemistry was used for all reactions (Applied Biosystems, Foster City, CA). Sequences were run on the ABI3700 platform (Applied Biosystems). The isolated haustoria (PTDH) and germinated urediniospores (PTDG) cDNA sequences contributed by M. Dickinson were end-sequenced by Syngenta UK and run through the Stackpack program, which aligns identical sequences to give a non-redundant output. The sequences were also "trimmed" to remove any contaminating vector sequences that may interfere with database searching.
ESTs database analysis
Generated ESTs were trimmed via the Consed suite of programs , including the 'Cross-Match' algorithm  using the UniVec database (http://www.ncbi.nlm.nih.gov/VecScreen/UniVec.html) to remove vector sequences and low quality reads. Sequences shorter than 50 bp were discarded and all were visually inspected. The resulting quality sequences, including the previously generated Pt ESTs  and a small set of unigene sequences from infected wheat material predicted to be fungal , were assembled into a non-redundant unigene set via the CAP3 assembly program . Databases for the various comparisons were comprised of publicly available and assembled sequences from wheat where care was taken to include only sequences derived from "clean", uninfected material, and from various fungi and oomycetes such as the COGEME db, http://cogeme.ex.ac.uk/. More specific collections were from the closely related Pgt genomic resources, such as the 20,567 predicted protein sequences or the 18,241 when disregarding many TE-related protein calls (release Jan 2010) and over 80,000 ESTs ; http://www.broadinstitute.org/annotation/genome/puccinia_group/MultiHome.html), as well as 30,000 Pt genomic BAC clone end-reads and initial pyrosequencing (454)-generated Pt genome fragments (http://www.ebi.ac.uk/ena/data/view/Taxon:208348). A collection of 3,297 putative fungal Pst EST sequences was collected from NCBI dbEST [17, 20, 21]. Mlp genomic resources were collected from JGI (http://genome.jgi-psf.org/Mellp1/Mellp1.home.html). U. maydis genomic resources were retrieved from the Munich Information Center for Protein Sequences (MIPS), MUMDB (http://mips.helmholtz-muenchen.de/genre/proj/ustilago).
The 17,684 new Pt EST sequences presented in this article can be retrieved from dbEST at NCBI (http://www.ncbi.nlm.nih.gov/sites/entrez?db=nucest) under accession numbers GR487994 to GR505442 and GR911120 to GR911355. Contig (assembled) unigene sequences are available from the Transcriptome Shotgun Assembly Sequence Database, TSA, at NCBI (http://www.ncbi.nlm.nih.gov/Genbank/TSA.html). Accession numbers are also given in Additional file 5.
This work was supported by a Strategic Grant from the National Science and Engineering Research Council of Canada to BS and GB (Grant No. STPGP 270337-2003) and by funding from the Agriculture & Agri-Food Canada 'Canadian Crop Genomics Initiative', the USDA-ARS-CRIS (#5430-21000-006-00D to JF). The research of MD was supported by a BBSRC-CASE studentship with Syngenta UK, which funded Helen Meakin, whose contribution is gratefully acknowledged. The authors would like to thank Dr. Li Huang for her assistance in the construction of a haustorial library used in this work. The authors also greatly acknowledge the funding from NSF/USDA grant 2008-35600-04693 to support the work done at the Broad Institute, Cambridge, MA, in generating the initial P. triticina genome reads (C. Cuomo, L. Szabo, J. Fellers and G. Bakkeren).
- Samborski DJ: Wheat leaf rust. The Cereal Rusts Vol II: Diseases, distribution, epidemiology and control. Edited by: Roelfs AP, Bushnell WR. 1985, Orlando, Fl.: Academic Press, 39-59.View Article
- Horton JS, Bakkeren G, Klosterman SJ, Garcia-Pedrajas M, Gold SE: Genetics of morphogenesis in Basidiomycetes. Applied Mycology and Biotechnology Genes and Genomics. Edited by: Arora DK, Berka R. 2005, Dordrecht: Elsevier, 5: 353-422. full_text.
- Bolton MD, Kolmer JA, Garvin DF: Wheat leaf rust caused by Puccinia triticina. Mol Plant Pathol. 2008, 9 (5): 563-575. 10.1111/j.1364-3703.2008.00487.x.PubMedView Article
- Chong J, Harder DE, Rohringer R: Cytochemical studies on Puccinia graminis f.sp. tritici in a compatible wheat host. II. Haustorium mother cell walls at the host cell penetration site, haustorial walls, and the extrahaustorial matrix. Can J Bot. 1986, 64 (11): 2561-2575. 10.1139/b86-339.View Article
- Samborski DJ, Dyck PL: Inheritance of virulence in wheat leaf rust on the standard differential wheat varieties. Can J Genet Cytol. 1968, 10: 24-32.View Article
- Mendgen K, Struck C, Voegele RT, Hahn M: Biotrophy and rust haustoria. Physiol Mol Plant Pathol. 2000, 56 (4): 141-145. 10.1006/pmpp.2000.0264.View Article
- Harder DE: Developmental ultrastructure of hyphae and spores. The cerreal rusts Vol I Origins, specificity, structure and physiology. Edited by: Bushnell WR, Roelfs AP. 1984, Orlando, FL: Academic Press, Inc, I: 333-373.
- Mendgen K: Development and physiology of teliospores. The cereal rusts. Edited by: Bushnell WR, Roelfs AP. 1984, New York: Academic Press, Inc, 1: 375-398.
- Anikster Y: Binucleate basidiospores; a general rule in rust fungi. Trans Br Mycol Soc. 1983, 81 (3): 624-626. 10.1016/S0007-1536(83)80136-3.View Article
- Snetselaar KM, Mims CW: Light and electron microscopy of Ustilago maydis hyphae in maize. Mycol Res. 1994, 98 (3): 347-355. 10.1016/S0953-7562(09)80463-2.View Article
- Gold RE, Mendgen K: Rust basidiospore germlings and disease initiation. The Fungal spore and disease initiation in plants and animals. Edited by: Garry T Cole, Harvey C. 1991, Hoch New York : Plenum Press, 67-99. c1991 pView Article
- Zhang L, Meakin H, Dickinson M: Isolation of genes expressed during compatible interactions between leaf rust (Puccinia triticina) and wheat using cDNA-AFLP. Mol Plant Pathol. 2003, 4 (6): 469-477. 10.1046/j.1364-3703.2003.00192.x.PubMedView Article
- Thara VK, Fellers JP, Zhou JM: In planta induced genes of Puccinia triticina. Mol Plant Pathol. 2003, 4 (1): 51-56. 10.1046/j.1364-3703.2003.00142.x.PubMedView Article
- Hu G, Linning ROB, McCallum B, Banks T, Cloutier S, Butterfield Y, Liu J, Kirkpatrick R, Stott J, Yang G, Smailus D, Jones S, Marra M, Schein J, Bakkeren G: Generation of a wheat leaf rust, Puccinia triticina, EST database from stage-specific cDNA libraries. Mol Plant Pathol. 2007, 8 (4): 451-467. 10.1111/j.1364-3703.2007.00406.x.PubMedView Article
- Broeker K, Bernard F, Moerschbacher BM: An EST library from Puccinia graminis f. sp. tritici reveals genes potentially involved in fungal differentiation. FEMS Microbiol Lett. 2006, 256 (2): 273-281. 10.1111/j.1574-6968.2006.00127.x.PubMedView Article
- Dracatos PM, Dumsday JL, Olle RS, Cogan NOI, Dobrowolski MP, Fujimori M, Roderick H, Stewart AV, Smith KF, Forster JW: Development and characterization of EST-SSR markers for the crown rust pathogen of ryegrass (Puccinia coronata f.sp. lolii). Genome. 2006, 49 (6): 572-583. 10.1139/G06-006.PubMedView Article
- Zhang Y, Qu Z, Zheng W, Liu B, Wang X, Xue X, Xu L, Huang L, Han Q, Zhao J, Kang Z: Stage-specific gene expression during urediniospore germination in Puccinia striiformis f. sp tritici. BMC Genomics. 2008, 9: 203-10.1186/1471-2164-9-203.PubMed CentralPubMedView Article
- Ling P, Wang M, Chen X, Campbell KG: Construction and characterization of a full-length cDNA library for the wheat stripe rust pathogen (Puccinia striiformis f. sp. tritici). BMC Genomics. 2007, 8 (1): 145-10.1186/1471-2164-8-145.PubMed CentralPubMedView Article
- Zhong S, Leng Y, Friesen TL, Faris JD, Szabo LJ: Development and Characterization of Expressed Sequence Tag-Derived Microsatellite Markers for the Wheat Stem Rust Fungus Puccinia graminis f. sp. tritici. Phytopathology. 2009, 99 (3): 282-289. 10.1094/PHYTO-99-3-0282.PubMedView Article
- Ma J, Huang X, Wang X, Chen X, Qu Z, Huang L, Kang Z: Identification of expressed genes during compatible interaction between stripe rust (Puccinia striiformis) and wheat using a cDNA library. BMC Genomics. 2009, 10 (1): 586-10.1186/1471-2164-10-586.PubMed CentralPubMedView Article
- Yin C, Chen X, Wang X, Han Q, Kang Z, Hulbert S: Generation and analysis of expression sequence tags from haustoria of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici. BMC Genomics. 2009, 10 (1): 626-10.1186/1471-2164-10-626.PubMed CentralPubMedView Article
- Jakupovic M, Heintz M, Reichmann P, Mendgen K, Hahn M: Microarray analysis of expressed sequence tags from haustoria of the rust fungus Uromyces fabae. Fungal Genet Biol. 2006, 43 (1): 8-19. 10.1016/j.fgb.2005.09.001.PubMedView Article
- 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.PubMed CentralPubMedView Article
- Puthoff DP, Neelam A, Ehrenfried ML, Scheffler BE, Ballard L, Song Q, Campbell KB, Cooper B, Tucker ML: Analysis of Expressed Sequence Tags from Uromyces appendiculatus Hyphae and Haustoria and Their Comparison to Sequences from Other Rust Fungi. Phytopathology. 2008, 98 (10): 1126-1135. 10.1094/PHYTO-98-10-1126.PubMedView Article
- Besemer J, Lomsadze A, Borodovsky M: GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001, 29 (12): 2607-2618. 10.1093/nar/29.12.2607.PubMed CentralPubMedView Article
- Lukashin AV, Borodovsky M: GeneMark.hmm: new solutions for gene finding. Nucleic Acids Res. 1998, 26 (4): 1107-1115. 10.1093/nar/26.4.1107.PubMed CentralPubMedView Article
- Besemer J, Borodovsky M: Heuristic approach to deriving models for gene finding. Nucleic Acids Res. 1999, 27 (19): 3911-3920. 10.1093/nar/27.19.3911.PubMed CentralPubMedView Article
- Antonov I, Borodovsky M: GeneTack: Frameshift identification in protein coding sequences by the Viterbi algorithm. J Bioinformatics Comp Biol. 2010, 8 (3): 1-17.
- Soanes DM, Skinner W, Keon J, Hargreaves J, Talbot NJ: Genomics of phytopathogenic fungi and the development of bioinformatic resources. Mol Plant-Microb Interact. 2002, 15 (5): 421-427. 10.1094/MPMI.2002.15.5.421.View Article
- Schmid R, Blaxter ML: annot8r: GO, EC and KEGG annotation of EST datasets. BMC Bioinformatics. 2008, 9: 180-10.1186/1471-2105-9-180.PubMed CentralPubMedView Article
- van der Merwe MM, Walker J, Ericson L, Burdon JJ: Coevolution with higher taxonomic host groups within the Puccinia/Uromyces rust lineage obscured by host jumps. Mycol Res. 2008, 112 (12): 1387-1408. 10.1016/j.mycres.2008.06.027.PubMedView Article
- Schirawski J, Mannhaupt G, Münch K, Brefort T, Schipper K, Doehlemann G, Di Stasio M, Rössel N, Mendoza-Mendoza A, Pester D, Müller O, Winterberg B, Meyer E, Ghareeb H, Wollenberg T, Münsterkötter M, Wong P, Walter M, Stukenbrock E, Güldener U, Kahmann R: Pathogenicity Determinants in Smut Fungi Revealed by Genome Comparison. Science. 2010, 330 (6010): 1546-1548. 10.1126/science.1195330.PubMedView Article
- Joly D, Feau N, Tanguay P, Hamelin R: Comparative analysis of secreted protein evolution using expressed sequence tags from four poplar leaf rusts (Melampsora spp.). BMC Genomics. 2010, 11 (1): 422-10.1186/1471-2164-11-422.PubMed CentralPubMedView Article
- Hu G, Kamp A, Linning R, Naik S, Bakkeren G: Complementation of Ustilago maydis MAPK mutants by a wheat leaf rust, Puccinia triticina homolog: potential for functional analyses of rust genes. Mol Plant-Microbe Interact. 2007, 20 (6): 637-647. 10.1094/MPMI-20-6-0637.PubMedView Article
- Mewes HW, Dietmann S, Frishman D, Gregory R, Mannhaupt G, Mayer KF, Münsterkötter M, Ruepp A, Spannagl M, Stumpflen V, Rattei T: MIPS: analysis and annotation of genome information in 2007. Nucleic Acids Res. 2008, 36: D196-201. 10.1093/nar/gkm980.PubMed CentralPubMedView Article
- García-Pedrajas MD, Nadal M, Denny T, Baeza-Montañez L, Paz Z, Gold SE: DelsGate: A Robust and Rapid Method for Gene Deletion. Methods Mol Biol. 2010, 638: 55-76.PubMedView Article
- Anikster Y, Eilam T, Mittelman L, Szabo LJ, Bushnell WR: Pycnial nectar of rust fungi induces cap formation on pycniospores of opposite mating type. Mycologia. 1999, 91 (5): 858-870. 10.2307/3761539.View Article
- Narisawa K, Yamaoka Y, Katsuya K: Mating type of isolates derived from the spermogonial state of Puccinia coronata var. coronata. Mycoscience. 1994, 35 (2): 131-135. 10.1007/BF02318489.View Article
- Bagnat M, Simons K: Cell surface polarization during yeast mating. Proc Natl Acad Sci USA. 2002, 99 (22): 14183-14188. 10.1073/pnas.172517799.PubMed CentralPubMedView Article
- Aguilar PS, Heiman MG, Walther TC, Engel A, Schwudke D, Gushwa N, Kurzchalia T, Walter P: Structure of sterol aliphatic chains affects yeast cell shape and cell fusion during mating. Proc Natl Acad Sci USA. 2010, 107 (9): 4170-4175. 10.1073/pnas.0914094107.PubMed CentralPubMedView Article
- Lee N, D'Souza CA, Kronstad JW: Of Smuts, Blasts, Mildews, and Blights:cAMP Signaling in Phytopathogenic Fungi. Annu Rev Phytopathol. 2003, 41: 399-427. 10.1146/annurev.phyto.41.052002.095728.PubMedView Article
- Cervantes-Chavez JA, Kronberg F, Passeron S, Ruiz-Herrera J: Regulatory role of the PKA pathway in dimorphism and mating in Yarrowia lipolytica. Fungal Genet Biol. 2009, 46 (5): 390-399. 10.1016/j.fgb.2009.02.005.PubMedView Article
- Ho EC, Cahill MJ, Saville BJ: Gene discovery and transcript analyses in the corn smut pathogen Ustilago maydis: expressed sequence tag and genome sequence comparison. BMC Genomics. 2007, 8 (1): 334-10.1186/1471-2164-8-334.PubMed CentralPubMedView Article
- Link TI, Voegele RT: Secreted proteins of Uromyces fabae: similarities and stage specificity. Mol Plant Pathol. 2008, 9 (1): 59-66.PubMed
- Müller 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: S63-S70.View Article
- Wang X, Bakkeren G, McCallum B: Virulence and molecular polymorphisms of the wheat leaf rust fungus Puccinia triticina in Canada from 1997 to 2007. Botany. 2010, 88 (6): 575-589. 10.1139/B10-034.View Article
- Long DL, Kolmer JA: A North American system of nomenclature for Puccinia recondita f. sp. tritici. Phytopathology. 1989, 79 (5): 525-529. 10.1094/Phyto-79-525.View Article
- Webb CA, Szabo LJ, Bakkeren G, Garry C, Staples RC, Eversmeyer M, Fellers JP: Transient expression and insertional mutagenesis of Puccinia triticina using biolistics. Funct Integr Genomics. 2006, 6 (3): 250-260. 10.1007/s10142-005-0009-9.PubMedView Article
- Hahn M, Mendgen K: Isolation of ConA binding of haustoria from different rust fungi and comparison of their surface qualities. Protoplasma. 1992, 170 (3-4): 95-103. 10.1007/BF01378785.View Article
- Li WL, Faris JD, Muthukrishnan S, Liu DJ, Chen PD, Gill BS: Isolation and characterization of novel cDNA clones of acidic chitinases and β-1,3-glucanases from wheat spikes infected by Fusarium graminearum. Theoretical Appl Genet. 2001, 102 (2): 353-362. 10.1007/s001220051653.View Article
- Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8 (3): 195-202.PubMedView Article
- Gordon D, Desmarais C, Green P: Automated finishing with autofinish. Genome Res. 2001, 11 (4): 614-625. 10.1101/gr.171401.PubMed CentralPubMedView Article
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView Article
- Soanes DM, Talbot NJ: Comparative genomic analysis of phytopathogenic fungi using expressed sequence tag (EST) collections. Mol Plant Pathol. 2006, 7 (1): 61-70. 10.1111/j.1364-3703.2005.00317.x.PubMedView Article
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.PubMedView Article
- Sacadura NT, Saville BJ: Gene expression and EST analyses of Ustilago maydis germinating teliospores. Fungal Genet Biol. 2003, 40 (1): 47-64. 10.1016/S1087-1845(03)00078-1.PubMedView Article
- Nugent KG, Choffe K, Saville BJ: Gene expression during Ustilago maydis diploid filamentous growth: EST library creation and analyses. Fungal Genet Biol. 2004, 41 (3): 349-360. 10.1016/j.fgb.2003.11.006.PubMedView Article
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.