Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Stage-specific gene expression during urediniospore germination in Puccinia striiformis f. sp tritici

Yonghong Zhang1, Zhipeng Qu1, Wenming Zheng2, Bo Liu1, Xiaojie Wang1, Xiaodan Xue1, Liangsheng Xu3, Lili Huang1, Qingmei Han1, Jie Zhao1 and Zhensheng Kang1*

Author Affiliations

1 College of Plant Protection and Shaanxi Key Laboratory of Molecular Biology for Agriculture, Northwest A&F University, Yangling, Shaanxi, 712100, P.R. China

2 College of Life Sciences, Henan Agriculture University, Zhengzhou, Henan, 450002, P.R. China

3 College of Life Sciences, Northwest A&F University, Yangling, Shaanxi, 712100, P.R. China

For all author emails, please log on.

BMC Genomics 2008, 9:203  doi:10.1186/1471-2164-9-203


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/9/203


Received:26 September 2007
Accepted:1 May 2008
Published:1 May 2008

© 2008 Zhang et al; licensee 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.

Abstract

Background

Puccinia striiformis f. sp. tritici is an obligate biotrophic pathogen that causes leaf stripe rust on wheat. Although it is critical to understand molecular mechanisms of pathogenesis in the wheat stripe rust fungus for developing novel disease management strategies, little is known about its genome and gene functions due to difficulties in molecular studies with this important pathogen. To identify genes expressed during early infection stages, in this study we constructed a cDNA library with RNA isolated from urediniospores of P. striiformis f. sp. tritici germinated for 10 h.

Results

A total of 4798 ESTs were sequenced from the germinated urediniospore library and assembled into 315 contigs and 803 singletons. About 23.9% and 13.3% of the resulting 1118 unisequences were homologous to functionally characterized proteins and hypothetical proteins, respectively. The rest 62.8% unisequences had no significant homologs in GenBank. Several of these ESTs shared significant homology with known fungal pathogenicity or virulence factors, such as HESP767 of the flax rust and PMK1, GAS1, and GAS2 of the rice blast fungus. We selected six ESTs (Ps28, Ps85, Ps87, Ps259, Ps261, and Ps159) for assaying their expression patterns during urediniospore germination and wheat infection by quantitative real-time PCR. All of them had the highest transcript level in germinated urediniospores and a much less transcript level in un-germinated urediniospores and infected wheat tissues (1–7 dpi). The transcript level of Ps159 increased at later infection stages (6–7 dpi). Our data indicated that these genes were highly expressed in germinated urediniospores and may play important roles in fungal-plant interactions during early infection stages in the wheat stripe rust fungus.

Conclusion

Genes expressed in germinated urediniospores of P. striiformis f. sp. tritici were identified by EST analysis. Six of them were confirmed by quantitative real-time PCR assays to be highly expressed in germinated urediniospores.

Background

Wheat stripe rust is one of the most important diseases of wheat throughout the world. It is a major constraint in wheat production and is a serious threat to food security worldwide [1,2]. Puccinia striiformis Westend f. sp. tritici Eriks is the causal agent of wheat stripe rust. It is an obligate biotrophic basidiomycete with an incomplete life cycle. Because of the high variability in the pathogen population, race-specific resistance in many newly developed cultivars often fails within a few years of cultivation and results in severe yield losses. In the past century, most of the P. striiformis f. sp. tritici studies had been focused on the identification of physiological races, virulence variation, and ultrastructural and histological examinations [1,3,4].

Unlike the wheat stem rust fungus, the sexual stage of P. striiformis f. sp. tritici has not been identified. Urediniospore is the most common spore form that has been observed in the wheat stripe rust fungus, which is strictly dependent on living host cells for growth and reproduction. To date, no stable transformation system has been established for P. striiformis f. sp. tritici and other Puccinia species. These biological characteristics make molecular and genetic studies of this important fungus relatively more difficult. Although much progresses have been achieved in researches on its genetic diversity, population structure, and evolution [5-8], there is only very limited knowledge about genes involved in the initial infection and biotrophic growth stages of the wheat stripe rust fungus. Such knowledge is critical for understanding infection mechanisms of this important pathogen and developing better disease management strategies.

Various genomic approaches, such as expressed sequence tag (EST) [9], serial analysis of gene expression (SAGE) [10], and massive parallel signature sequencing (MPSS) [11], have been widely used in genome-wide gene expression studies in various organisms. EST analysis was the first method used for rapid identification of expressed genes [9]. It has been employed to identify genes that are expressed in various tissues, cell types, or developmental stages in different organisms [12-14]. The availability of EST sequences has accelerated molecular characterization of genes of interest and provided sequence information for microarray design.

EST analyses have been conducted in a few rust fungi. To examine gene expression during infection of broad bean by Uromyces fabae, Hahn and his colleagues sequenced ESTs from purified haustoria. A major shift in gene expression was observed between urediniospore germination and the biotrophic growth stage [15,16]. A total of 25,558 ESTs have been generated from 13 cDNA libraries representing various stages of Puccinia triticina, including resting and germinated urediniospores, appressorium and haustorium formation stages during infection of a susceptible wheat cultivar, and infected leaves of a resistant wheat cultivar. While 38% unigenes matched sequences in various databases and collections, the annotation rates were low for ESTs from germinated urediniospores (4%) and appressoria (2%). Gene sets obtained from these different libraries appeared to be remarkably different, suggesting drastic reprogramming of the transcriptome during these major differentiation processes [17]. In this study, we generated a cDNA library from germinated urediniospores of P. striiformis f. sp. tritici. A total of 4798 ESTs were sequenced to generate 1118 unisequences (uniseqs). The majority of these ESTs (over 60%) had no significant homologs in GenBank, indicating that many of them may represent genes unique to P. striiformis f. sp. tritici. Several of these ESTs share significant homology with known fungal pathogenicity or virulence factors, such as HESP767 of the flax rust and PMK1 of the rice blast fungus. The high transcript level of six selected genes in germinated urediniospores was confirmed by quantitative real-time PCR (qRT-PCR) assays. Some of these genes highly expressed in germinated urediniospores may be important for early infection processes in P. striiformis f. sp. tritici.

Results

Generating ESTs of germinated urediniospores

To identify genes related to early events of urediniospore germination and differentiation, we incubated freshly harvested urediniospores of P. striiformis f. sp tritici (race CY32) in sterile distilled water in plastic plates. After incubating for 10 h at 9°C, about 60% of urediniospores attached to the plastic surface and produced long, un-branched germ tubes (Fig. 1). Some of these germ tubes displayed various degree of swelling at the tip (Fig. 1). Germinated and un-germinated urediniospores were collected at 10 h and used for RNA isolation. A directional cDNA library consisting of 6.05 × 105 primary clones was constructed with the λTripEx2 vector (Clontech). On average, about 95% of the cDNA clones had inserts longer than 0.2 kb. The insert size varied from 200–3500 bp, with an average of 750 bp. A total of 5500 random cDNA clones were sequenced from the 5'-end to obtain 4810 quality reads or ESTs.

thumbnailFigure 1. Urediniospores and germ tubes of P. striiformis f. sp tritici (race CY32) after incubating for 10 h on plastic surface. Solid arrowheads marked germinated and un-germinated urediniospores. Dashed arrowhead indicated the swelling of germ tubes at the tip. Bar = 100 μm.

Data analysis and functional classification

Before clustering analysis, ESTs containing no insert sequence or insert shorter than 100 bp were removed. The remaining 4798 ESTs were aligned and assembled into 315 contigs (containing two or more ESTs) and 803 singletons. A total of 1118 unisequences (uniseqs) were submitted to dbEST at NCBI [18] under GenBank accession numbers ES321929ES323046. Most of the uniseqs had the insert size of 200–900 bp. Only 36 uniseqs were longer than 1000 bp. The G+C content of these ESTs ranged from 40.74% to 53.73%, with an average of 45.08%, which was similar to the G+C content of ESTs in P. graminis [12].

The number of individual ESTs belonging to each contig ranged from 2 to 608 [see Additional file 1]. Nine most abundant contigs contained more than 100 ESTs each, suggesting a high transcript level of the corresponding genes. Contig Ps303 was represented by 608 clones, more than any other contigs [see Additional file 2]. It is homologous to a Saccharomyces cerevisiae gene encoding a hypothetical protein. Five of these most abundant contigs [see Additional file 2] had no homologs in EST or protein database bases. Three other abundant contigs displayed limited homology with entries in GenBank. Contig Ps253 had weak homology with a putative secreted protein in Ixodes scapularis. Contig Ps28 consisting of 128 ESTs was homologous to a differentiation-related protein Inf24 from Uromyces appendiculatus. Inf24 was identified as infection structure protein that was highly expressed during urediniospore germination [19]. For contig Ps314, it consisted of 154 ESTs and was homologous to a predicted protein of Kluyveromyces lactis.

Additional file 1. The redundancy of P. striiformis f. sp tritici ESTs derived from the germinated urediniospore cDNA library. The number of contigs consisting of 1, 2, 3–5, 6–10, 11–35, and more than 36 ESTs was presented by the columns.

Format: PDF Size: 56KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 2. The most abundant ESTs identified in the P. striiformis f. sp. tritici germinated urediniospore library. These data provided represent the nine most abundant expressed ESTs.

Format: DOC Size: 42KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

All uniseqs were subjected to similarity searches against sequences in the non-redundant protein (nr) database at GenBank using the BLASTX algorithm. InterProScan was used to analyze uniseqs with no BLASTX matches for known protein domains. The vast majority (703) of the uniseqs had no significant homolog in GenBank by BLASTX search. Among 415 uniseqs displayed similarity (E-value < 10e-5) to entries in the nr database, about 56% of them were similar to genes coding for proteins with unknown functions (not classified in Gene Ontology). The frequency of orphan sequences in P. striiformis f. sp. tritici ESTs is similar to that has been observed in other fungal EST projects [20].

Among the 415 uniseqs that exhibited similarity to entries in the non-redundant protein database at GenBank, about 40.72% shared homology to proteins from filamentous fungi while the rest were homologous to proteins from a wide variety of organisms, including yeast, bacteria, nematodes, plants, insects, and animals [see Additional file 3]. The uniseqs with significant similarity to hypothetical proteins were placed in the unclassified protein category. Based on the results from BlastX and InterProScan searches, ESTs with significant matches were categorized according to their putative functions (Table 1, for detail [see Additional file 4]).

Additional file 3. Uniseqs with significant homology to genes from different organisms. These data provided display the number of uniseqs have homology to genes from diversity organisms.

Format: DOC Size: 38KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 4. Uniseqs with similarity (BLASTX, E-value < 10-5 or InterProScan, E-value < 10-5) to proteins in public databases were grouped into functional categories according to Gene Ontology. These data provided represent the original EST number and best hit.

Format: DOC Size: 2.2MB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Table 1. Major functional categories of P. striiformis ESTs from the germinated urediniospore library.

A total of 267 uniseqs were identified to be homologous to proteins with known cellular functions. Figure 2 listed the putative functions of these uniseqs and their occurrence. Over 70% of them were functionally related to primary metabolism (42.59%) and protein (21.58%) or (11.11%) RNA synthesis. About 10% were involved in cellular signaling and defense responses. The ESTs with no hits in GenBank were searched against the genome sequence of the wheat stem rust P. graminis [21]. Many of them had homologous sequences in the genome of the stem rust fungus (data not shown). A total of 195 uniseqs had no homologous sequences in the stem rust fungus genome. These ESTs may represent genes unique to P. striiformis f. sp. tritici.

thumbnailFigure 2. Classification of the P. striiformis f. sp. tritici uniseqs homologous to proteins of known functions. A total of 267 uniseqs (23.9% of 1118 uniseqs) with significant homologs in GenBank were categorized according to the functions of their homologs. Over 60% of them were involved in primary metabolism (42.59%) and protein synthesis (21.58%).

Putative pathogenicity or virulence factors identified in the germinated urediniospore ESTs

To date, many genes important for fungal pathogenesis have been identified in plant and human pathogens [22-24]. Homology searches revealed that several ESTs were homologous to known fungal pathogenicity or virulence factors (E-value < 1e-5), including key components of signal transduction pathways, transporters, and genes involved in infection-related morphogenesis (Table 2). For three genes related to signal transduction, Ps1728 encodes a putative adenylate cyclase that has been shown to be essential for pathogenicity in Magnaporthe grisea, Ustilago maydis, and several other fungi [25-28]. Ps259 encodes a putative Gα subunit that is similar to MagB of M. grisea and Gpa1 of Cryptococcus neoformans [29,30]. The MAP kinase encoded by Ps261 is homologous to Pmk1 of M. grisea, which is important for regulating plant infection processes in a number of phytopathogenic fungi [31-33].

Ps8 and Ps28 may be related to the infection structure differentiation (Table 2). Both of them were homologous to Inf24 of U. appendiculatus and appeared to be members of a multigene family. Inf24 is highly expressed during urediniospore germination. Microinjection of antisense oligonucleotides of Inf24 inhibits appressorium development in U. appendiculatus [19]. Both Ps28 and Ps228 were contigs represented by over 100 ESTs, indicating that these two genes are highly expressed in germinated urediniospores, and may have similar functions during early stages of plant infection in P. striiformis f. sp. tritici.

Contigs Ps6 and Ps159 encoded proteins homologous to gEgh16 of Blumeria graminis [22], which is highly expressed at early infection stages (16 h). In M. grisea, its homologs, GAS1 and GAS2, are important virulence factors that are involved in appressorial penetration [34]. Both GAS1 and GAS2 are specifically expressed during appressorium formation. We also identified two contigs, Ps300 and Ps5709, that were putative ATP-binding cassette (ABC) transporter genes. In M. grisea and other plant pathogens, several ABC transporter genes have been implicated in fungal-plant interactions, possibly by controlling the efflux of phytotoxic fungal metabolites or plant defensive compounds [35].

Table 2. Uniseqs with significant homology to known fungal pathogenicity factors

We also identified an EST, Ps238 (Table 2), that encodes a putative copper- and zinc-superoxide dismutase. In some fungal pathogens, such as Candida albicans, superoxide dismutase is required for tolerance to oxidative stress and full virulence [36]. Contig Ps85 encoded a protein consisting of 169 amino acid residues. It is homologous a putative cell surface antigen gene in C. albicans [37] and contains a CFEM (

    C
ommon in several
    F
ungal
    E
xtracellular
    M
embrane proteins) domain that is unique to filamentous ascomycetes [38]. In M. grisea, the CFEM domain-containing protein Pth11 is known to play an important role in pathogenesis as a receptor [39]. Pth11 is not required for appressorium morphogenesis in vitro but is involved in host surface recognition. The protein encoded by contig Ps85 may play a similar role in P. striiformis f. sp. tritici, and is involved in appressorium differentiation.

ESTs homologous to known rust effectors or avirulence genes

The first Avr protein identified in rust fungi was AvrL567 of the flax rust fungus Melampsora lini [40]. Recently, a number of additional avirulence elicitors have been identified in the flax rust fungus [41]. In the ESTs of the wheat stripe rust fungus, we identified a candidate avirulence elicitor encoded by contig Ps87, which had significant homology with HESP767. In the flax rust fungus, HESP767 is expressed in haustoria and may function as an effector in interactions with the host plant [41]. Ps87 contained the full length open reading frame (ORF) that was predicted to encode a protein of 85 amino acid residues. It had a typical signal peptide at the N-terminus (data not shown) as predicted by SignalP 3.0 [42]. Similar to HESP767, Ps87 may function as an effector or avirulence factor during wheat infection in P. striiformis f. sp. tritici.

Assaying the transcript level of six selected ESTs

Six uniseqs, Ps28 (a differentiation-related protein), Ps85 (a CFEM domain-containing protein), Ps159 (gEgh16 homolog), Ps259 (G-alpha subunit), Ps261 (Pmk1 homolog), and Ps87 (HESP767 homolog), were selected for assaying their transcript levels during different developmental and infection stages by quantitative real-time PCR (qRT-PCR). All these six genes were expressed in urediniospores, germinated urediniospores, and infected wheat tissues harvested at 1 through 7 days post inoculation (dpi). However, their transcription levels were much higher in germinated urediniospores than in un-germinated urediniospores or infected wheat tissues (Fig. 3; Fig. 4).

thumbnailFigure 3. Assays for the transcript levels of ESTs Ps28, Ps85, Ps159, Ps259, and Ps261 in P. striiformis f. sp. tritici. RNA samples were isolated from urediniospores, germinated urediniospores, or leaves of wheat cultivar Huixianhong inoculated with CY32 and collected at the indicated time points. The expression level of these genes was calculated by the comparative Ct method with the actin gene of P. striiformis f. sp. tritici as the endogenous reference for normalization. Relative quantification was computed with their expression levels in different stages in comparison to that in urediniospores. The average and standard error were calculated from three biological replicates. S, urediniospores; GS, germinated urediniospores; dpi, days post inoculation.

thumbnailFigure 4. Assays for the transcript level of Ps87 during incompatible interactions. RNA samples were isolated from urediniospores, germinated urediniospores, or leaves of wheat cultivar Shuiyuan 11 inoculated with CY23 and collected at the indicated time points. The expression level of Ps87 was calculated by the comparative Ct method with the actin gene of P. striiformis f. sp. tritici as the endogenous reference for normalization. Relative quantification was computed with its expression levels in different stages in comparison to that in urediniospores. The average and standard error were calculated from three biological replicates. S, urediniospores; GS, germinated urediniospores; dpi, days post inoculation.

Although the transcript level of these genes in infected wheat tissues was relatively low, they had different expression patterns during plant infection. In comparison with the others, Ps159 had a relatively high transcript level in infected wheat plants. Its expression was up-regulated from 2 to 6 dpi and then slightly decreased at 7 dpi. The transcript level of Ps85 and Ps28 increased at early infection stages but decreased after 5 dpi. The transcript level of Ps259 and Ps261 in wheat tissues was lower than that of the other ESTs examined. Those two genes may be specifically and highly expressed during urediniospore germination. These data suggest that a considerable reprogramming of gene expression occurs during urediniospore germination, early stages of plant infection, and the biotrophic growth in P. striiformis f. sp. tritici.

During incompatible interactions, Ps87 also displayed an expression pattern that decreased at early stages but slightly increased after 3 dpi (Fig. 4). The high transcript level of Ps87 in germinated urediniospores suggests that it may be involved in germ tube growth and appressorium differentiation. It is possible that Ps87 also plays a role in late stages of plant infection, probably for host-pathogen recognitions.

Discussion

For plant pathogens with little or no history of genetic research, single-pass sequencing of random cDNA clones as in EST projects represents a relatively inexpensive and rapid procedure for finding novel genes and information about their expression. To date, extensive EST databases have been established for various plant pathogenic fungi, including M. grisea [43], Phakopsora pachyrhizi [52], and F. oxysporum [44]. For wheat rust fungi, only small EST libraries have been described for P. triticina [45] and P. graminis [12]. Recently, a full-length cDNA library has been constructed with RNA isolated from un-germinated urediniospores of a race PST-78 isolate of P. striiformis f. sp. tritici [46]. Among 196 random cDNA clones sequenced from this library, only 73 of them (37.2%) have homologs of known functions. Most of them are involved in various housekeeping functions. A few ESTs encoding hypothetical proteins have homologs in other plant pathogenic fungi [46]. Different from the cDNA library constructed in this study, the library constructed by Ling and colleagues may have high percentage of storage transcripts in un-germinated urediniospores.

Urediniospore germination represents an early stage of the interaction of P. striiformis f. sp. tritici with host plants. As an obligate biotroph, germ tubes of the wheat stripe rust fungus fail to progress beyond this stage of development in vitro [47,48]. In this study, we constructed and sequenced a cDNA library with RNA isolated from germinated urediniospores. Different from other filamentous fungi, a relatively high level of redundancy was found in this P. striiformis f. sp. tritici library. Among the 4798 ESTs generated in this study, 3995 of them could be assembled into contigs. A few contigs consisting of over 100 EST clones may represent genes that are highly expressed during urediniospore germination. When the uniseqs were queried against the nr protein database, about 76.2% of them had no significant homology with proteins of known functions, which may be related to the difficulty of functional characterization of genes in rust fungi. Among these contigs encoding proteins of known functions, 18 of them are ribosomal proteins. ESTs corresponding to genes encoding ribosomal proteins also are abundant in other fungal ESTs, such as those of B. graminis [49] and M. grisea [50].

The 1118 uniseqs identified in this study contained a wide range of genes involved in different cellular functions. The most abundant genes during urediniospore germination were those involved in metabolic activities as well as those responsible for protein biosynthesis, which accounted for 10.29% and 5.28% of the uniseqs, respectively. A few genes are involved in RNA synthesis, cell signal and communication, cell structure and growth, and cell/organism defense, indicating that active metabolism and protein synthesis is important for urediniospore germination and germ tube growth. Several ESTs, including Ps8, Ps28, and Ps228, shared similarity to differentiation-related protein Inf24 from U. appendiculatus, which is specifically expressed during urediniospore germination [51]. Ps88 has similarity to a chitin synthase from P. graminis f. sp. tritici. Ps55 and Ps1010 were homologous to a deacetylase and a cell wall organization and biogenesis-related protein from C. neoformans, respectively [52]. These three genes may be involved in the modification of cell wall during urediniospore germination. The contig Ps5712 encodes a putative calcium/calmodulin-dependent protein kinase. The calcium-signaling pathway has been implicated in regulating appressorium formation in M. grisea and other plant pathogens [53,54]. In the wheat stripe rust fungus, it may regulate germ tube emergence and infection structure differentiation.

Expression patterns of six selected uniseqs, Ps28, Ps85, Ps87, Ps159, Ps261, and Ps295, were examined by qRT-PCR analysis. The transcript level of these genes in germinated urediniospores was several-folds higher than in un-germinated urediniospores or infected wheat tissues, indicating that gene expression in the wheat stripe rust fungus changes dramatically during urediniospore germination and parasitic growth in wheat plants (Fig. 3; Fig. 4). These observations are consistent with what have been reported in U. fabae [16].

A number of genes that are important for plant infection have been identified in various phytopathogenic fungi [55]. Several of them are key components of conserved signaling pathways. In the ESTs generated in this study, we identified several genes involved in signal transduction, such as Ps259 and Ps261 that encode a G-alpha subunit and a Pmk1 homolog, respectively. Similarly, several signaling components are identified in P. triticina ESTs, including the Pmk1 homolog PtMAPK1 [17]. The PtMAPK1 gene has increased transcript levels during urediniospore germination and plant infection. When expressed in U. maydis, it can complement the defects of the kpp2 mutant in mating and plant infection [56].

Other putative fungal virulence or pathogenicity factors identified in this EST analysis included contigs Ps8, Ps28, Ps6, Ps159, and Ps87. Contigs Ps8 and Ps28 shared similarity to differentiation-related protein Inf24 from U. appendiculatus. In U. appendiculatus, microinjection of Inf24 antisense fragment into germinated urediniospores blocked its transcription and appressorium formation [19]. Injection with sense fragments has no effect on responses of germ tubes to the topographical stimuli and development of appressoria. Injection of antisense fragments into mature appressoria has no inhibitory effects on the development of subsequent infection structures. The Inf24 protein appears to play a critical role in the germ tube before the formation of appressoria. The presence of Inf24 homologs in this EST library of germinated urediniospores suggested that they may have similar functions in the wheat stripe rust fungus.

Ps6 and Ps159 have similarity to GAS1 and GAS2 of M. grisea [34], respectively. GAS1 and GAS2 are specifically expressed during appressorium formation and important for appressorial penetration in M. grisea. They are homologous to gEgh16 of B. graminis and members of a small protein family unique to filamentous fungi. In several fungal pathogens, gEgh16 homologs are expressed at early plant infection stages [23,34]. In the wheat stripe rust and leaf rust fungi [17], gEgh16 homologs may be involved in early stages of appressorium development. Although Ps159 was transcribed at a level higher than those of other ESTs in infected plant tissues, its transcript was much more abundant in germinated urediniospores (Fig. 3). After germinating for 10 h, the germ tube tips tend to swell in P. striiformis f. sp. tritici (Fig. 1). Genes involved in appressorium formation may be highly expressed at this stage.

When reliable transformation systems become available for P. striiformis f. sp. tritici in the future, it will be important to determine the functions of these putative pathogenicity factors identified in this EST library by generating gene knockout or silencing mutants. Some of these genes highly expressed in germinated urediniospores may play important roles in fungal-plant interactions during early infection stages in the wheat stripe rust fungus.

Conclusion

A cDNA library was constructed from germinated urediniospores of P. striiformis f. sp. tritici. A total of 4798 ESTs were sequenced and assembled into 315 contigs and 803 singletons. About 62.8% of the resulting 1118 uniseqs had no significant homologs in GenBank. Among the uniseqs with assigned functions, over 70% of them were functionally related to primary metabolism and protein or RNA synthesis. The rest were associated with various cellular functions. Several of them were homologous to known fungal pathogenicity factors or effector proteins. The high transcript level of six selected ESTs in germinated urediniospores was confirmed by qRT-PCR. Genes identified in this study to be highly expressed in germinated urediniospores may be important for early infection processes in P. striiformis f. sp. tritici.

Methods

Strains and culture conditions

P. striiformis f. sp. tritici strain CY32 was inoculated and propagated on wheat cultivar Huixianhong as described previously [57]. Fresh urediniospores were harvested from infected wheat plants and resuspended in sterile distilled water (6 mg/200 ml). After incubating in a 50 × 125 cm dish at 9°C for 10 h, germinated urediniospores were collected with a spatula, frozen in liquid nitrogen, and stored at -80°C. For isolating RNA from infected plants, wheat leaves of susceptible cultivar Huixianhong and resistant cultivar Shuiyuan 11 were inoculated with CY32 urediniospores and harvested at 1, 2, 3, 4, 5, 6, and 7 days post inoculation (dpi).

cDNA synthesis, library construction, and DNA sequencing

Total RNA was isolated from 100 mg of germinated urediniospores with the RNeasy Plant Mini RNA purification kit (QIAGEN, Germany) following the instruction provided by the manufacturer. The SMARTTM cDNA library construction kit (Clontech, USA) was used for cDNA synthesis and library construction. The MaxPlaxTM Lambda Packaging Extract (Epicentre, USA) was used for in vitro packaging and transfection of Escherichia coli strain XL1-Blue. The resulting directional library consisting of 6.05 × 105 primary clones was amplified and stored in 15% glycerol at -80°C. Randomly selected cDNA clones were sequenced with primer seq1 (5' CGACTCTAGACTCGAGCAAG 3') from the 5'-end with an ABI 3130-XL DNA sequencer.

Sequence analysis and bioinformatics

Sequence reads longer than 100 base pairs were processed by cross_match [58] that allows the removal of contaminated sequences. Repeats and low complexity sequences were masked using RepeatMasker [59]. The resulting quality trimmed sequences were extracted and assembled with the CAP3 assembler, and viewed with Consed [60]. Statistics of the assemblies were generated by perl scripts using BioPerl modules.

For functional classification, the resulting unisequences were searched against the NCBI non-redundant protein database using the BLASTX program [61,62]. The unisequences with significant BLASTX matches were classified according to their likely cellular functions following the general categories outlined by the Gene Ontology Consortium [63]. InterProScan was used to search for protein domains in ESTs with no significant homologs in BLASTX searches.

Isolating RNA from urediniospores, germinated urediniospores, and infected wheat leaves

Standard protocols [64] were used to extract total RNA from urediniospores, germinated urediniospores, and wheat leaves inoculated with the stripe rust fungus. Two micrograms of total RNA each was used for cDNA synthesis with the SuperScript First-strand Synthesis System (Invitrogen) following the instruction provided by the manufacture. The resulting 1st-strand cDNA products synthesized with RNA from urediniospores, germinated urediniospores, and infected wheat leaves were used as the templates for qRT-PCR assays.

Quantitative real-time PCR analysis

Four sets of primers and fluorescently labeled TaqMan probes [see Additional file 5] were used to assay the transcript levels of ESTs Ps28, Ps85, Ps159, Ps259, Ps261, and Ps87. The expected sizes of resulting PCR products were 107, 65, 81, 73, 103, and 179 bp for Ps28, Ps85, Ps87, Ps159, Ps259, and Ps261, respectively. The fluorogenic probes were labeled at the 5'end with the fluorescent reporter dye FAM (6-carboxy-fluorescein) and modified at the 3'end with the quencher dye TAMRA (6-carboxy-tetramethylrhodamine). Transcript abundance was assessed with three independent biological replicates. Each real-time PCR mixture (25 μl) contained 12.5 μl premix (Premix Ex TaqTM, Takara), 7.5 μl distilled H2O, 1 μl 10 μM forward primer, 1 μl 10 μM reverse primer, 1 μl TaqMan probe, and 2 μl the cDNA template. The Smart CyclerSystem (Cepheid, USA) was used for PCR reactions that consisted of 10 s at 95°C (1 cycle) and 40 cycles of 95°C for 5 s and 60°C for 20 s. The transcript level of selected ESTs was calculated by the 2-ΔΔCT method [65] with the actin gene of P. striiformis f. sp. tritici as the endogenous reference for normalization. Relative quantification of these ESTs was computed with their transcript levels in different stages in comparison to that in urediniospores.

Additional file 5. Sequences of primers and probes used for qRT-PCR. The primer pairs, probes and reference sequences used for qRT-PCR were listed.

Format: DOC Size: 51KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Authors' contributions

YZ constructed the cDNA library, participated in EST sequencing and sequence analysis, conducted qRT-PCR assays, and prepared the manuscript; ZQ, LX and XW contributed to EST sequencing and Blast searches. WZ participated in various experiments; BL and XX assisted in qRT-PCR assays; LH, JZ and QH coordinated the project; ZK conceived of this study, contributed to experimental designs, provided supervisions, and revised the manuscript; All authors read and approved the final manuscript.

Acknowledgements

We thank Drs. Larry Dunkle and Jin-Rong Xu at Purdue University for critical reading of this manuscript. This study was supported by grants from the National Basic Research Program of China (No.2006CB101901); the Cultivation Fund of the Key Scientific and Technical Innovation Project, Ministry of Education of China (No. 2004-295); Nature Science Foundation of China (No.30671350); the Program for Changjiang Scholars and Innovative Research Team in University, Ministry of Education of China (No.200558); and the 111 Project from Ministry of Education of China (B07049).

References

  1. Li ZQ, Zeng SM: Wheat stripe rust in China. Beijing, Agriculture Press; 2002. OpenURL

  2. Line RF: Stripe rust of wheat and barley in North America: A retrospective historical review.

    Annu Rev Phytopathol 2002, 40:75-118. PubMed Abstract | Publisher Full Text OpenURL

  3. Kang ZS: Ultrastructure of Plant Pathogenic Fungi. Beijing, China Science and Technology Press; 1995. OpenURL

  4. Staples RC: Research on the rust fungi during the twentieth century.

    Annu Rev Phytopathol 2000, 38:49-69. PubMed Abstract | Publisher Full Text OpenURL

  5. Chen X, Line RF, Leung H: Relationship between virulence variation and DNA polymorphism in Puccinia striiformis.

    Phytopathology 1993, 83:1489-1496. Publisher Full Text OpenURL

  6. Steele KA, Humphreys E, Wellings CR, Dickinson MJ: Support for a stepwise mutation model for pathogen evolution in Australasian Puccinia striiformis f. sp. tritici by use of molecular markers.

    Plant Pathol 2001, 50:174-180. Publisher Full Text OpenURL

  7. Zheng WM, Liu F, Kang Z, Chen SY, Li ZQ, Wu L: AFLP fingerprinting of Chinese epidemic strains of Puccinia striiformis f. sp. Tritici.

    Prog Nat Sci 2001, 11:587-593. OpenURL

  8. Lorys MMA, Villaréal CL, Claude de Vallavieille-Pope CN: Genetic Variability in Puccinia striiformis f. sp. tritici Populations Sampled on a Local Scale during Natural Epidemics.

    Appl Environ Microbiol 2002, 68:6138-6145. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Adams M, Kerlavage AR, Fleischmann RD, Fuldner RA, Bult CJ, Lee NH, Kirkness EF, Weinstock KG, Gocayne JD, White O: Initial assessment of human gene diversity and expression patterns based upon 83 million nucleotides of cDNA sequence.

    Nature 1995, 377:173-174. PubMed Abstract | Publisher Full Text OpenURL

  10. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: Serial analysis of gene expression.

    Science 1995, 270:484-487. PubMed Abstract | Publisher Full Text OpenURL

  11. Brenner S, Johnson M, Bridgham J, Golda G, Lloyd DH, Johnson D, Luo S, McCurdy S, Foy M, Ewan M: Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays.

    Nat Biotechnol 2000, 18:630-634. PubMed Abstract | Publisher Full Text OpenURL

  12. 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:273-281. PubMed Abstract | Publisher Full Text OpenURL

  13. Ronning CM, Stegalkina SS, Ascenzi RA, Bougri O, Hart AL, Utterbach TR, Vanaken SE, Riedmuller SB, White JA, Cho J: Comparative analyses of potato expressed sequence tag libraries.

    Plant Physiol 2003, 131:419-429. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Ogihara Y, Mochida K, Nemoto Y, Murai K, Yamazaki Y, Shin-IT KY: Correlated clustering and virtual display of gene expression patterns in the wheat life cycle by large-scale statistical analyses of expressed sequence tags.

    Plant J 2003, 33:1001-1011. PubMed Abstract | Publisher Full Text OpenURL

  15. Hahn M, Mendgen K: Characterization of in planta-induced rust genes isolated from a haustorium-specific cDNA library.

    Mol Plant Microbe Interact 1997, 10:427-437. PubMed Abstract | Publisher Full Text OpenURL

  16. Jakupoviæ 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:8-19. PubMed Abstract | Publisher Full Text OpenURL

  17. Hu GG, Linning R, 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:451-467. Publisher Full Text OpenURL

  18. National Center for Biotechnology Information [http://www.ncbi.nlm.nih.gov] webcite

  19. Baria F, Correa AJ, Staples RC, Hoch HC: Microinjected antisense Inf24 oligo-nucleotides inhibit appressorium development in Uromyces.

    Mycol Res 1998, 102:1513-1518. Publisher Full Text OpenURL

  20. Yoder OC, Turgeon BG: Fungal genomics and pathogenicity.

    Curr Opin Plant Biol 2001, 4:315-321. PubMed Abstract | Publisher Full Text OpenURL

  21. Broad Institute of MIT and Harvard [http://www.broad.mit.edu/] webcite

  22. Ahn N, Kin S, Choi W, Im KH, Lee YH: Extracellular matrix protein gene, EMP1, is required for appressorium formation and pathogenicity of the rice blast fungus, Magnaporthe grisea.

    Mol Cells 2003, 17:166-173. OpenURL

  23. Justesen A, Somerville S, Christiansen S, Giese H: Isolation and characterization of two novel genes expressed in germinating conidia of the obligate biotroph Erysiphe graminis.

    Gene 1996, 170:131-135. PubMed Abstract | Publisher Full Text OpenURL

  24. Xue CY, Bahn YS, Cox GM, Heitman J: G Protein-coupled Receptor Gpr4 Senses Amino Acids and Activates the cAMP-PKA Pathway in Cryptococcus neoformans.

    Mol Biol Cell 2006, 17:667-679. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Adachi K, Hamer JE: Divergent cAMP signaling pathways regulate growth and pathogenesis in the rice blast fungus Magnaporthe grisea.

    Plant Cell 1998, 10:1361-1374. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Barhoom S, Sharon A: cAMP regulation of "pathogenic" and "saprophytic" fungal spore germination.

    Fungal Genet Biol 2004, 41:317-326. PubMed Abstract | Publisher Full Text OpenURL

  27. Gold SE, Brogdon SM, Mayorga ME, Kronstad JW: The Ustilago maydis regulatory subunit of a cAMP-dependent protein kinase is required for gall formation in maize.

    Plant Cell 1997, 9:1585-1594. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Liebmann B, Müller M, Braun A, Brakhage AA: The Cyclic AMP-Dependent Protein Kinase A Network Regulates Development and Virulence in Aspergillus fumigatus.

    Infect Immun 2004, 72:5193-5203. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Bahn YS, Hicks JK, Giles SS, Cox GM, Heitman J: Adenylyl cyclase-associated protein Aca1 regulates virulence and differentiation of Cryptococcus neoformans via the cyclic AMP-protein kinase A cascade.

    Eukaryot Cell 2004, 3:1476-1491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Liu S, Dean RA: G protein α subunit genes control growth, development, and pathogenicity of Magnaporthe grisea.

    Mol Plant Microbe In 1997, 10:1075-1086. Publisher Full Text OpenURL

  31. Mey G, Oeser B, Lebrun MH, Tudzynski P: The biotrophic, non-appressoria forming grass pathogen Claviceps purpurea needs a Fus3/Pmk1 homologous MAP kinase for colonization of rye ovarian tissue.

    Mol Plant-Microbe Interact 2002, 15:303-312. PubMed Abstract | Publisher Full Text OpenURL

  32. Xu JR: MAP kinases in fungal pathogens.

    Fungal Genet Biol 2000, 31:137-152. PubMed Abstract | Publisher Full Text OpenURL

  33. Xu JR, Hamer JE: MAP kinase and cAMP signaling regulate infection structure formation and pathogenic growth in the rice blast fungus Magnaporthe grisea.

    Genes Dev 1996, 10:2696-2706. PubMed Abstract | Publisher Full Text OpenURL

  34. Xue C, Park G, Choi W, Zheng L, Dean RA, Xu JR: Two novel fungal virulence genes specifically expressed in appressoria of the rice blast fungus.

    Plant Cell 2002, 14:2107-2119. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Urban M, Bhargava T, Hamer JE: An ATP-driven efflux pump is a novel pathogenicity factor in rice blast disease.

    EMBO J 1999, 18:512-521. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Hwang CS, Rhie GE, Oh JH, Huh WK, Yim HS, Kang SO: Copper- and zinc-containing superoxide dismutase (Cu/ZnSOD) is required for the protection of Candida albicans against oxidative stresses and the expression of its full virulence.

    Microbiology 2002, 148:3705-3713. PubMed Abstract | Publisher Full Text OpenURL

  37. Jones T, Federspiel NA, Chibana H, Dungan J, Kalman S, Magee BB, Newport G, Thorstenson YR, Agabian N, Magee PT, Davis RW, Scherer S: The diploid genome sequence of Candida albicans.

    Proc Natl Acad Sci USA 2004, 101:7329-7334. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. 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:118-121. PubMed Abstract | Publisher Full Text OpenURL

  39. DeZwaan TM, Carroll AM, Valent B, Sweigard JA: Magnaporthe grisea pth11p is a novel plasma membrane protein that mediates appressorium differentiation in response to inductive substrate cues.

    Plant Cell 1999, 11:2013-2030. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Dodds PN, Lawrence GJ, Catanzariti AM, Ayliffe MA, Ellis JG: The Melampsora lini AvrL567 avirulence genes are expressed in haustoria and their products are recognized inside plant cells.

    Plant Cell 2004, 16:755-768. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. 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:243-256. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. SignalP 3.0 Server [http://www.cbs.dtu.dk/services/SignalP/] webcite

  43. Ebbole DJ, Jin Y, Thon M, Pan H, Bhattarai E, Thomas T, Dean R: Gene discovery and gene expression in the rice blast fungus, Magnaporthe grisea: analysis of expressed sequence tags.

    Mol Plant-Microbe Interact 2004, 17:1337-1347. PubMed Abstract | Publisher Full Text OpenURL

  44. Iida Y, Ohara T, Tsuge T: Identification of genes up-regulated during conidiation of Fusarium oxysporum through expressed sequence tag analysis.

    Fungal Genet Biol 2006, 43:179-189. PubMed Abstract | Publisher Full Text OpenURL

  45. Thara VK, Fellers JP, Zhou JM: In Planta induced genes of Puccinia triticina.

    Mol Plant Pathol 2003, 4:51-56. Publisher Full Text OpenURL

  46. Ling P, Wang MN, Chen XM, 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:145. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  47. Mendgen K, Struck C, Voegele RT, Hahn M: Biotrophy and rust haustoria.

    Physiol Mol Plant P 2000, 56:141-145. Publisher Full Text OpenURL

  48. Mendgen K, Hahn M: Plant infection and the establishment of fungal biotrophy.

    Trends Plant Sci 2002, 7:352-356. PubMed Abstract | Publisher Full Text OpenURL

  49. Thomas SW, Rasmussen SW, Glaring MA, Rouster JA, Christiansen SK, Oliver RP: Gene identification in the obligate fungal pathogen Blumeria graminis by expressed sequence tag analysis.

    Fungal Genet Biol 2001, 33:195-211. PubMed Abstract | Publisher Full Text OpenURL

  50. Rauyaree P, Choi W, Fang E, Blackmon B, Dean R: Genes expressed during early stages of rice infection with the rice blast fungus Magnaporthe grisea.

    Mol Plant Pathol 2001, 2:347-354. Publisher Full Text OpenURL

  51. Xuei XL, bhairi S, Staples RC, Yoder OC: INF56 represents a family of differentiation-specific genes from Uromyces appendiculatus.

    Curr Genet 1993, 24:84-88. PubMed Abstract | Publisher Full Text OpenURL

  52. Posada-buitrago ML, Frederick RD: Expressed sequence tag analysis of the soybean rust pathogen Phakopsora pachyrhizi.

    Fungal Genet Biol 2005, 42:949-962. PubMed Abstract | Publisher Full Text OpenURL

  53. Ahn IP, Suh SC: Calcium/calmodulin-dependent signaling for prepenetration development in Cochliobolus miyabeanus infecting rice.

    J Gen Plant Pathol 2007, 73:113-120. Publisher Full Text OpenURL

  54. Karunakaran M, Nair NV, Rho HS, Lee YH: Effects of ca2+/calmodulin mediated signaling on spore germination and appressorium development in Colletotrichum falcatum went.

    Acta Phytopathol Entomol Hungarica 2006, 41:249-260. Publisher Full Text OpenURL

  55. Idnurm A, Howlett BJ: Pathogenicity genes of phytopathogenic fungi.

    Mol Plant Pathol 2001, 2:241-255. Publisher Full Text OpenURL

  56. 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:637-647. PubMed Abstract | Publisher Full Text OpenURL

  57. Kang Z, Li Z, Shang H, CHONG J, Rohringer R: Ultrastructure and cytochemistry of intercellular hyphae of wheat stripe rust.

    Acta Mycol Sinica 1993, 12:208-213. OpenURL

  58. Green Group [http://www.phrap.org/] webcite

  59. Repeat Masker Server Home Page [http://repeatmasker.genome.washington.edu/] webcite

  60. Huang X, Madan A: CAP 3: A DNA sequence assembly program.

    Genome Res 1999, 9:868-877. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic Local Alignment Search Tool.

    J Mol Biol 1990, 215:403-410. PubMed Abstract | Publisher Full Text OpenURL

  62. Altschul SF, Madden TL, Schaffer AA, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

    Nucl Acid Res 1997, 25:3389-3402. Publisher Full Text OpenURL

  63. The Gene Ontology [http://www.geneontology.org/] webcite

  64. Sambrook J, Fitsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, Cold Spring Harbor Press; 1989. OpenURL

  65. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method.

    Methods 2001, 25:402-408. PubMed Abstract | Publisher Full Text OpenURL