Gene expression during Ustilago maydis diploid filamentous growth: EST library creation and analyses

Transcriptome analysis of smut fungi reveals widespread intergenic transcription and conserved antisense transcript expression

Background

Biotrophic fungal plant pathogens cause billions of dollars in losses to North American crops annually. The model for functional investigation of these fungi is Ustilago maydis. Its 20.5 Mb annotated genome sequence has been an excellent resource for investigating biotrophic plant pathogenesis. Expressed-sequence tag libraries and microarray hybridizations have provided insight regarding the type of transcripts produced by U. maydis but these analyses were not comprehensive and there were insufficient data for transcriptome comparison to other smut fungi. To improve transcriptome annotation and enable comparative analyses, comprehensive strand-specific RNA-seq was performed on cell-types of three related smut species: U. maydis (common smut of corn), Ustilago hordei (covered smut of barley), and Sporisorium reilianum (head smut of corn).

Results

In total, >1 billion paired-end sequence reads were obtained from haploid cell, dikaryon and teliospore RNA of U. maydis, haploid cell RNA of U. hordei, and haploid and dikaryon cell RNA of S. reilianum. The sequences were assembled into transfrags using Trinity, and updated gene models were created using PASA and categorized with Cufflinks Cuffcompare. Representative genes that were predicted for the first time with these RNA-seq analyses and genes with novel annotation features were independently assessed by reverse transcriptase PCR. The analyses indicate hundreds more predicted proteins, relative to the previous genome annotation, could be produced by U. maydis from altered transcript forms, and that the number of non-coding RNAs produced, including transcribed intergenic sequences and natural antisense transcripts, approximately equals the number of mRNAs. This high representation of non-coding RNAs appears to be a conserved feature of the smut fungi regardless of whether they have RNA interference machinery. Approximately 50% of the identified NATs were conserved among the smut fungi.

Conclusions

Overall, these analyses revealed: 1) smut genomes encode a number of transcriptional units that is twice the number of annotated protein-coding genes, 2) a small number of intergenic transcripts may encode proteins with characteristics of fungal effectors, 3) the vast majority of intergenic and antisense transcripts do not contain ORFs, 4) a large proportion of the identified antisense transcripts were detected at orthologous loci among the smut fungi, and 5) there is an enrichment of functional categories among orthologous loci that suggests antisense RNAs could have a genome-wide, non-RNAi-mediated, influence on gene expression in smut fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3720-8) contains supplementary material, which is available to authorized users.

Investigating the Ustilago maydis/Zea mays pathosystem: transcriptional responses and novel functional aspects of a fungal calcineurin regulatory B subunit

The sustainable control of basidiomycete biotrophic plant pathogenesis requires an understanding of host responses to infection, as well as the identification and functional analysis of fungal genes involved in disease development. The creation and analysis of a suppressive subtractive hybridization (SSH) cDNA library from Ustilago maydis-infected Zea mays seedlings enabled the identification of fungal and plant genes expressed during disease development, and uncovered new insights into the interactions of this model system. Candidate U. maydis pathogenesis genes were identified by using the current SSH cDNA library analysis, and by knowledge generated from previous cDNA microarray and comparative genomic analyses. These identifications were supported by the independent determination of transcript level changes in different cell-types and during pathogenic development. The basidiomycete specific um01632, the highly in planta expressed um03046 (zig1), and the calcineurin regulatory B subunit (um10226, cnb1), were chosen for deletion experiments. um01632 and zig1 mutants showed no difference in morphology and did not have a statistically significant impact on pathogenesis. cnb1 mutants had a distinct cell division phenotype and reduced virulence in seedling assays. Infections with reciprocal wild-type×Δcnb1 haploid strain crosses revealed that the wild-type allele was unable to fully compensate for the lack of a second cnb1 allele. This haploinsufficiency was undetected in other fungal cnb1 mutational analyses. The reported data improves U. maydis genome annotation and expands on the current understanding of pathogenesis genes in this model basidiomycete.

Ustilago maydis natural antisense transcript expression alters mRNA stability and pathogenesis

Ustilago maydis infection of Zea mays leads to the production of thick-walled diploid teliospores that are the dispersal agent for this pathogen. Transcriptome analyses of this model biotrophic basidiomycete fungus identified natural antisense transcripts (NATs) complementary to 247 open reading frames. The U. maydis NAT cDNAs were fully sequenced and annotated. Strand-specific RT-PCR screens confirmed expression and identified NATs preferentially expressed in the teliospore. Targeted screens revealed four U. maydis NATs that are conserved in a related fungus. Expression of NATs in haploid cells, where they are not naturally occurring, resulted in increased steady-state levels of some complementary mRNAs. The expression of one NAT, as-um02151, in haploid cells resulted in a twofold increase in complementary mRNA levels, the formation of sense-antisense double-stranded RNAs, and unchanged Um02151 protein levels. This led to a model for NAT function in the maintenance and expression of stored teliospore mRNAs. In testing this model by deletion of the regulatory region, it was determined that alteration in NAT expression resulted in decreased pathogenesis in both cob and seedling infections. This annotation and functional analysis supports multiple roles for U. maydis NATs in controlling gene expression and influencing pathogenesis.

Ustilago maydis transcript features identified through full-length cDNA analysis

Ustilago maydis is the model for investigating basidiomycete biotrophic plant pathogens. To further the annotation of its genome, 12,943 full-length cDNA sequences were used to construct databases for the promoter and untranslated regions of U. maydis genes. A subset of clones was sequenced to determine full cDNA sequences. These and the original ESTs were assembled into contigs representing 3,058, or 45%, of the predicted U. maydis genes. The new sequencing allowed the confirmation of 2,842 gene models, 690 of which contain an intron. The use of full-length cDNA clone sequences ensured that untranslated regions were physically linked to the open reading frames (ORFs), not merely aligned upstream of the start of transcription. Identified sequence features include: (1) over 500 potential short upstream ORFs, (2) 95 gene models that require further annotation, (3) one new potential ORF, (4) varying GC content in different gene regions, (5) a WebLogo motif for the start of translation, (6) the correlation of UTR length with transcript representation in cDNA libraries and with gene function categories, (7) a relationship between natural antisense transcripts and UTR length that differs from that of Saccharomyces cerevisiae, (8) a potential relationship between DNA replication and the control of transcription, and (9) new insights regarding mechanisms for the control of transcription and mRNA maturation in U. maydis.

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

Background

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.

Results

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.

Conclusions

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.

Detection of antisense RNA transcripts by strand-specific RT-PCR

Comprehensive genome annotation requires extensive cDNA analysis. This analysis has identified natural antisense transcripts (NATs), which are distinct from the microRNAs, siRNAs, and piRNAs, in a number of diverse eukaryotes. This wide conservation supports the possibility of an important role for NATs in regulating cellular processes. Investigating their roles requires the confirmation of expressed sequence tag (EST) data and the detection of antisense transcripts in distinct cellular backgrounds. This chapter describes the use of a reverse transcription polymerase chain reaction (RT-PCR) method for the detection of antisense transcripts. The protocol was designed to reduce the number of first strand synthesis reactions during screening for antisense transcripts through the utilization of antisense directed primers and oligo dT to prime first strand synthesis. These results are further confirmed using sense and antisense directed primers in first strand synthesis. Results indicate that optimization of the screens requires proper controls to confirm removal of gDNA contamination and to rule out self-priming as a source of first strand products.

Comparative gene expression analysis of susceptible and resistant near-isogenic lines in common wheat infected by Puccinia triticina

Gene expression after leaf rust infection was compared in near-isogenic wheat lines differing in the Lr10 leaf rust resistance gene. RNA from susceptible and resistant plants was used for cDNA library construction. In total, 55 008 ESTs were sequenced from the two libraries, then combined and assembled into 14 268 unigenes for further analysis. Of these ESTs, 89% encoded proteins similar to (E value of < or =10(-5)) characterized or annotated proteins from the NCBI non-redundant database representing diverse molecular functions, cellular localization and biological processes based on gene ontology classification. Further, the unigenes were classified into susceptible and resistant classes based on the EST members assembled from the respective libraries. Several genes from the resistant sample (14-3-3 protein, wali5 protein, actin-depolymerization factor and ADP-ribosylation factor) and the susceptible sample (brown plant hopper resistance protein, caffeic acid O-methyltransferase, pathogenesis-related protein and senescence-associated protein) were selected and their differential expression in the resistant and susceptible samples collected at different time points after leaf rust infection was confirmed by RT-PCR analysis. The molecular pathogenicity of leaf rust in wheat was studied and the EST data generated made a foundation for future studies.

Analyses of expressed sequence tags from the maize foliar pathogen Cercospora zeae-maydis identify novel genes expressed during vegetative, infectious, and reproductive growth

Background

The ascomycete fungus Cercospora zeae-maydis is an aggressive foliar pathogen of maize that causes substantial losses annually throughout the Western Hemisphere. Despite its impact on maize production, little is known about the regulation of pathogenesis in C. zeae-maydis at the molecular level. The objectives of this study were to generate a collection of expressed sequence tags (ESTs) from C. zeae-maydis and evaluate their expression during vegetative, infectious, and reproductive growth.

Results

A total of 27,551 ESTs was obtained from five cDNA libraries constructed from vegetative and sporulating cultures of C. zeae-maydis. The ESTs, grouped into 4088 clusters and 531 singlets, represented 4619 putative unique genes. Of these, 36% encoded proteins similar (E value ≤ 10^-05) to characterized or annotated proteins from the NCBI non-redundant database representing diverse molecular functions and biological processes based on Gene Ontology (GO) classification. We identified numerous, previously undescribed genes with potential roles in photoreception, pathogenesis, and the regulation of development as well as Zephyr, a novel, actively transcribed transposable element. Differential expression of selected genes was demonstrated by real-time PCR, supporting their proposed roles in vegetative, infectious, and reproductive growth.

Conclusion

Novel genes that are potentially involved in regulating growth, development, and pathogenesis were identified in C. zeae-maydis, providing specific targets for characterization by molecular genetics and functional genomics. The EST data establish a foundation for future studies in evolutionary and comparative genomics among species of Cercospora and other groups of plant pathogenic fungi.

Gene discovery and transcript analyses in the corn smut pathogen Ustilago maydis: expressed sequence tag and genome sequence comparison

Background

Ustilago maydis is the basidiomycete fungus responsible for common smut of corn and is a model organism for the study of fungal phytopathogenesis. To aid in the annotation of the genome sequence of this organism, several expressed sequence tag (EST) libraries were generated from a variety of U. maydis cell types. In addition to utility in the context of gene identification and structure annotation, the ESTs were analyzed to identify differentially abundant transcripts and to detect evidence of alternative splicing and anti-sense transcription.

Results

Four cDNA libraries were constructed using RNA isolated from U. maydis diploid teliospores (U. maydis strains 518 × 521) and haploid cells of strain 521 grown under nutrient rich, carbon starved, and nitrogen starved conditions. Using the genome sequence as a scaffold, the 15,901 ESTs were assembled into 6,101 contiguous expressed sequences (contigs); among these, 5,482 corresponded to predicted genes in the MUMDB (MIPS Ustilago maydis database), while 619 aligned to regions of the genome not yet designated as genes in MUMDB. A comparison of EST abundance identified numerous genes that may be regulated in a cell type or starvation-specific manner. The transcriptional response to nitrogen starvation was assessed using RT-qPCR. The results of this suggest that there may be cross-talk between the nitrogen and carbon signalling pathways in U. maydis. Bioinformatic analysis identified numerous examples of alternative splicing and anti-sense transcription. While intron retention was the predominant form of alternative splicing in U. maydis, other varieties were also evident (e.g. exon skipping). Selected instances of both alternative splicing and anti-sense transcription were independently confirmed using RT-PCR.

Conclusion

Through this work: 1) substantial sequence information has been provided for U. maydis genome annotation; 2) new genes were identified through the discovery of 619 contigs that had previously escaped annotation; 3) evidence is provided that suggests the regulation of nitrogen metabolism in U. maydis differs from that of other model fungi, and 4) Alternative splicing and anti-sense transcription were identified in U. maydis and, amid similar observations in other basidiomycetes, this suggests these phenomena may be widespread in this group of fungi. These advances emphasize the importance of EST analysis in genome annotation.

Generation of a wheat leaf rust, Puccinia triticina, EST database from stage-specific cDNA libraries

Thirteen cDNA libraries constructed from small amounts of leaf rust mRNA using optimized methods served as the source for the generation of 25 558 high-quality DNA sequence reads. Five life-cycle stages were sampled: resting urediniospores, urediniospores germinated over water or plant extract, compatible, interactive stages during appressorium or haustorium formation just before sporulation, and an incompatible interaction. mRNA populations were subjected to treatments such as full-length cDNA production, subtractive and normalizing hybridizations, and size selection methods combined with PCR amplification. Pathogen and host sequences from interactive libraries were differentiated in silico using cereal and fungal sequences, codon usage analyses, and by means of a partial prototype cDNA microarray hybridized with genomic DNAs. This yielded a non-redundant unigene set of 9760 putative fungal sequences consisting of 6616 singlets and 3144 contigs, representing 4.7 Mbp. At an E-value 10(-5), 3670 unigenes (38%) matched sequences in various databases and collections but only 694 unigenes (7%) were similar to genes with known functions. In total, 296 unigenes were identified as most probably wheat and ten as rRNA sequences. Annotation rates were low for germinated urediniospores (4%) and appressoria (2%). Gene sets obtained from the various life-cycle stages appear to be remarkably different, suggesting drastic reprogramming of the transcriptome during these major differentiation processes. Redundancy within contigs yielded information about possible expression levels of certain genes among stages. Many sequences were similar to genes from other rusts such as Uromyces and Melampsora species; some of these genes have been implicated in pathogenicity and virulence.

Analysis of expressed sequence tags from the fungus Aspergillus oryzae cultured under different conditions

We performed random sequencing of cDNAs from nine biologically or industrially important cultures of the industrially valuable fungus Aspergillus oryzae to obtain expressed sequence tags (ESTs). Consequently, 21 446 raw ESTs were accumulated and subsequently assembled to 7589 non-redundant consensus sequences (contigs). Among all contigs, 5491 (72.4%) were derived from only a particular culture. These included 4735 (62.4%) singletons, i.e. lone ESTs overlapping with no others. These data showed that consideration of culture grown under various conditions as cDNA sources enabled efficient collection of ESTs. BLAST searches against the public databases showed that 2953 (38.9%) of the EST contigs showed significant similarities to deposited sequences with known functions, 793 (10.5%) were similar to hypothetical proteins, and the remaining 3843 (50.6%) showed no significant similarity to sequences in the databases. Culture-specific contigs were extracted on the basis of the EST frequency normalized by the total number for each culture condition. In addition, contig sequences were compared with sequence sets in eukaryotic orthologous groups (KOGs), and classified into the KOG functional categories.

Generation, annotation, and analysis of ESTs from four different Trichoderma strains grown under conditions related to biocontrol

The functional genomics project "TrichoEST" was developed focused on different taxonomic groups of Trichoderma with biocontrol potential. Four cDNA libraries were constructed, using similar growth conditions, from four different Trichoderma strains: Trichoderma longibrachiatum T52, Trichoderma asperellum T53, Trichoderma virens T59, and Trichoderma sp. T78. In this study, we present the analysis of the 8,160 expressed sequence tags (ESTs) generated. Each EST library was independently assembled and 1,000-1,300 unique sequences were identified in each strain. First, we queried our collection of ESTs against the NCBI nonredundant database using the BLASTX algorithm. Moreover, using the Gene Ontology hierarchy, we performed the annotation of 40.9% of the unique sequences. Later, based on the EST abundance, we examined the highly expressed genes in the four strains. A hydrophobin was found as the gene expressed at the highest level in two of the strains, but we also found that other unique sequences similar to the HEX1, QID3, and NMT1 proteins were highly represented in at least two of the Trichoderma strains.

Genetics of morphogenesis and pathogenic development of Ustilago maydis

Ustilago maydis has emerged as an important model system for the study of fungi. Like many fungi, U. maydis undergoes remarkable morphological transitions throughout its life cycle. Fusion of compatible, budding, haploid cells leads to the production of a filamentous dikaryon that penetrates and colonizes the plant, culminating in the production of diploid teliospores within fungal-induced plant galls or tumors. These dramatic morphological transitions are controlled by components of various signaling pathways, including the pheromone-responsive MAP kinase and cAMP/PKA (cyclic AMP/protein kinase A) pathways, which coregulate the dimorphic switch and sexual development of U. maydis. These signaling pathways must somehow cooperate with the regulation of the cytoskeletal and cell cycle machinery. In this chapter, we provide an overview of these processes from pheromone perception and mating to gall production and sporulation in planta. Emphasis is placed on the genetic determinants of morphogenesis and pathogenic development of U. maydis and on the fungus-host interaction. Additionally, we review advances in the development of tools to study U. maydis, including the recently available genome sequence. We conclude with a brief assessment of current challenges and future directions for the genetic study of U. maydis.

EST analysis of cDNA libraries from the entomopathogenic fungus Beauveria (Cordyceps) bassiana. I. Evidence for stage-specific gene expression in aerial conidia, in vitro blastospores and submerged conidia

The entomopathogenic fungus Beauveria (Cordyceps) bassiana holds much promise as a pest biological control agent. B. bassiana produces at least three in vitro single cell infectious propagules, including aerial conidia, vegetative cells termed blastospores and submerged conidia, that display different morphological, biochemical and virulence properties. Populations of aerial conidia, blastospores and submerged conidia were produced on agar plates, rich liquid broth cultures and under conditions of nutrient limitation in submerged cultures, respectively. cDNA libraries were generated from mRNA isolated from each B. bassiana cell type and approximately 2,500 5' end sequences were determined from each library. Sequences derived from aerial conidia clustered into 284 contigs and 963 singlets, with those derived from blastospores and submerged conidia forming 327 contigs with 788 singlets, and 303 contigs and 1,079 contigs, respectively. Almost half (40-45 %) of the sequences in each library displayed either no significant similarity (e value >10(-4)) or similarity to hypothetical proteins found in the NCBI database. The expressed sequence tag dataset also included sequences representing a significant portion of proteins in cellular metabolism, information storage and processing, transport and cell processes, including cell division and posttranslational modifications. Transcripts encoding a diverse array of pathogenicity-related genes, including proteases, lipases, esterases, phosphatases and enzymes producing toxic secondary metabolites, were also identified. Comparative analysis between the libraries identified 2,416 unique sequences, of which 20-30 % were unique to each library, and only approximately 6 % of the sequences were shared between all three libraries. The unique and divergent representation of the B. bassiana transcriptome in the cDNA libraries from each cell type suggests robust differential gene expression profiles in response to environmental conditions.

Analysis of expressed sequence tags from the anamorphic basidiomycetous yeast, Pseudozyma antarctica, which produces glycolipid biosurfactants, mannosylerythritol lipids

Pseudozyma antarctica T-34 secretes a large amount of biosurfactants (BS), mannosylerythritol lipids (MEL), from different carbon sources such as hydrocarbons and vegetable oils. The detailed biosynthetic pathway of MEL remained unknown due to lack of genetic information on the anamorphic basidiomycetous yeasts, including the genus Pseudozyma. Here, in order to obtain genetic information on P. antarctica T-34, we constructed a cDNA library from yeast cells producing MEL from soybean oil and identified the genes expressed through the creation of an expressed sequence tags (EST) library. We generated 398 ESTs, assembled into 146 contiguous sequences. Based upon a BLAST search similarity cut-off of E<or=10(-5), 21.4% of all contigs were orphan, while 78.6% showed similarity to sequences in the protein database; 60.3% of all contiguous sequences shared significant identities to hypothetical protein of Ustilago maydis, which is a smut fungus and BS producer. Based on the gene expression study using real-time reverse transcriptase-PCR, the predicted genes, such as mannosyltranferase and acyltransferase, were demonstrated to be highly involved in MEL biosynthesis in soybean oil-grown cells.

Generation, annotation and analysis of ESTs from Trichoderma harzianum CECT 2413

BACKGROUND

The filamentous fungus Trichoderma harzianum is used as biological control agent of several plant-pathogenic fungi. In order to study the genome of this fungus, a functional genomics project called "TrichoEST" was developed to give insights into genes involved in biological control activities using an approach based on the generation of expressed sequence tags (ESTs).

RESULTS

Eight different cDNA libraries from T. harzianum strain CECT 2413 were constructed. Different growth conditions involving mainly different nutrient conditions and/or stresses were used. We here present the analysis of the 8,710 ESTs generated. A total of 3,478 unique sequences were identified of which 81.4% had sequence similarity with GenBank entries, using the BLASTX algorithm. Using the Gene Ontology hierarchy, we performed the annotation of 51.1% of the unique sequences and compared its distribution among the gene libraries. Additionally, the InterProScan algorithm was used in order to further characterize the sequences. The identification of the putatively secreted proteins was also carried out. Later, based on the EST abundance, we examined the highly expressed genes and a hydrophobin was identified as the gene expressed at the highest level. We compared our collection of ESTs with the previous collections obtained from Trichoderma species and we also compared our sequence set with different complete eukaryotic genomes from several animals, plants and fungi. Accordingly, the presence of similar sequences in different kingdoms was also studied.

CONCLUSION

This EST collection and its annotation provide a significant resource for basic and applied research on T. harzianum, a fungus with a high biotechnological interest.

Genetic linkage map and expression analysis of genes expressed in the lamellae of the edible basidiomycete Pleurotus ostreatus

Pleurotus ostreatus is an industrially cultivated basidiomycete with nutritional and environmental applications. Its genome contains 35 Mbp organized in 11 chromosomes. There is currently available a genetic linkage map based predominantly on anonymous molecular markers complemented with the mapping of QTLs controlling growth rate and industrial productivity. To increase the saturation of the existing linkage maps, we have identified and mapped 82 genes expressed in the lamellae. Their manual annotation revealed that 34.1% of the lamellae-expressed and 71.5% of the lamellae-specific genes correspond to previously unknown sequences or to hypothetical proteins without a clearly established function. Furthermore, the expression pattern of some genes provides an experimental basis for studying gene regulation during the change from vegetative to reproductive growth. Finally, the identification of various differentially regulated genes involved in protein metabolism suggests the relevance of these processes in fruit body formation and maturation.

Serial analysis of gene expression reveals conserved links between protein kinase A, ribosome biogenesis, and phosphate metabolism in Ustilago maydis

The switch from budding to filamentous growth is a key aspect of invasive growth and virulence for the fungal phytopathogen Ustilago maydis. The cyclic AMP (cAMP) signaling pathway regulates dimorphism in U. maydis, as demonstrated by the phenotypes of mutants with defects in protein kinase A (PKA). Specifically, a mutant lacking the regulatory subunit of PKA encoded by the ubc1 gene displays a multiple-budded phenotype and fails to incite disease symptoms, although proliferation does occur in the plant host. A mutant with a defect in a catalytic subunit of PKA, encoded by adr1, has a constitutively filamentous phenotype and is nonpathogenic. We employed serial analysis of gene expression to examine the transcriptomes of a wild-type strain and the ubc1 and adr1 mutants to further define the role of PKA in U. maydis. The mutants displayed changes in the transcript levels for genes encoding ribosomal proteins, genes regulated by the b mating-type proteins, and genes for metabolic functions. Importantly, the ubc1 mutant displayed elevated transcript levels for genes involved in phosphate acquisition and storage, thus revealing a connection between cAMP and phosphate metabolism. Further experimentation indicated a phosphate storage defect and elevated acid phosphatase activity for the ubc1 mutant. Elevated phosphate levels in culture media also enhanced the filamentous growth of wild-type cells in response to lipids, a finding consistent with PKA regulation of morphogenesis in U. maydis. Overall, these findings extend our understanding of cAMP signaling in U. maydis and reveal a link between phosphate metabolism and morphogenesis.

Comparative genomic analysis of phytopathogenic fungi using expressed sequence tag (EST) collections

SUMMARY We describe the analysis of 57 727 unique expressed sequence tags (ESTs) from 15 species of phytopathogenic and three species of saprophytic fungi. This resource is held within the COGEME phytopathogen EST database (http://cogeme.ex.ac.uk/). Comparative analysis was performed to investigate the differences between pathogenic and free-living fungi based on a substantial collection of expressed gene sequences and available, completed fungal genome sequences. We report that the expressed gene inventories of pathogenic fungi were not significantly more similar to each other than to those of free-living filamentous fungi. As expected, however, filamentous fungi as a group share more sequences in common than with the free-living yeast species Saccharomyces cerevisiae. Interestingly, ESTs of the obligate biotrophic fungus Blumeria graminis f. sp. hordei were more dissimilar to those of all other fungal species assessed, having a lower number of sequences in common with filamentous ascomycetes studied to date and also possessing a larger proportion of unisequences of unknown function. Our analysis of ESTs in the COGEME database enabled identification of a set of functional groups of genes that are more highly represented in the genomes of pathogenic fungi than non-pathogenic species.

Expressed sequence tag analysis of the soybean rust pathogen Phakopsora pachyrhizi

Soybean rust is caused by the obligate fungal pathogen Phakopsora pachyrhizi Sydow. A unidirectional cDNA library was constructed using mRNA isolated from germinating P. pachyrhizi urediniospores to identify genes expressed at this physiological stage. Single pass sequence analysis of 908 clones revealed 488 unique expressed sequence tags (ESTs, unigenes) of which 107 appeared as multiple copies. BLASTX analysis identified 189 unigenes with significant similarities (Evalue<10(-5)) to sequences deposited in the NCBI non-redundant protein database. A search against the NCBI dbEST using the BLASTN algorithm revealed 32 ESTs with high or moderate similarities to plant and fungal sequences. Using the Expressed Gene Anatomy Classification, 31.7% of these ESTs were involved in primary metabolism, 14.3% in gene/protein expression, 7.4% in cell structure and growth, 6.9% in cell division, 4.8% in cell signaling/cell communication, and 4.8% in cell/organism defense. Approximately 29.6% of the identities were to hypothetical proteins and proteins with unknown function.

Analysis of expressed sequence tags from the wheat leaf blotch pathogen Mycosphaerella graminicola (anamorph Septoria tritici)

Mycosphaerella graminicola is a major fungal pathogen of wheat as the causal agent of Septoria leaf blotch disease. As a first step toward a greater understanding of the mechanism of host infection we have generated, sequenced, and analyzed three M. graminicola EST libraries from conditions predicted to resemble independent phases of the host infection process, including one library generated from the fungus during interaction with its host. A total of 5180 ESTs were sequenced and clustered into 886 contigs and 2039 singletons to give a set of 2925 unique sequences (unisequences). BLASTX analysis revealed 33% of the unknown M. graminicola unisequences to be orphans. Very limited inter-library overlap of expression was seen with the majority of unisequences (contigs and singletons) being library-specific. Analysis of EST redundancy between libraries demonstrated a significant difference in gene expression in the three conditions. Comparisons made against fully sequenced genomes revealed most M. graminicola sequences to be homologous to genes present in both pathogenic and non-pathogenic Ascomycete filamentous fungi. A range of sequences having significant homology to verified pathogenicity/virulence genes (HvPV-genes) of either plant or mammalian fungal and Oomycete pathogens were also identified (<1e-20). The generation of, and the diversity present within, this EST collection will facilitate future efforts aimed at a more detailed study of the transcriptome of the fungus during host infection.

Differential gene expression during teliospore germination in Ustilago maydis

Ustilago maydis is a model fungal pathogen that induces the formation of tumors in maize. The tumor provides an environment for hyphal differentiation, leading to the formation of thick-walled, diploid teliospores. Such spores serve as a dispersal agent for smut and rust fungi, and their germination leads to new rounds of infection. The morphological changes that occur during teliospore germination in U. maydis have been described in detail. However, the specific molecular events that facilitate this process have not been identified. Through the construction and hybridization of microarrays containing a set of 3918 non-redundant cDNAs, we have identified genes that are differentially regulated during teliospore germination. Teliospores induced to germinate for 4 and 11 h were selected for comparison with dormant teliospores. Genes identified as differentially expressed included many that are presumably involved in as yet undescribed molecular events during teliospore germination, as well as characterized genes previously shown to be required for the process. This study represents the first large-scale investigation of changes in gene expression during teliospore germination.

Molecular cloning of developmentally specific genes by representational difference analysis during the fruiting body formation in the basidiomycete Lentinula edodes

To understand molecular mechanisms of the fruiting body development in basidiomycetes, we attempted to isolate developmentally regulated genes expressed specifically during the fruiting body formation of Lentinula edodes (Shiitake-mushroom). cDNA representational difference analysis (cDNA-RDA) between vegetatively growing mycelium and two developmental substages, primordium and mature fruiting body, resulted in an isolation of 105 individual genes (51 in primordium and 54 in mature fruiting body, respectively). A search of homology with the protein databases and two basidiomycetous genomes in Phanerochaete chrysosporium and Coprinopsis cinerea revealed that the obtained genes encoded various proteins similar to those involved in general metabolism, cell structure, signal transduction, and responses to stress; in addition, there were apparently several metabolic pathways and signal transduction cascades that could be involved in the fruiting body development. The expression products of several genes revealed no significant homologies to those in the databases, implying that those genes are unique in L. edodes and the encoding products may possess possible functions in the course of fruiting body development. RT-PCR analyses revealed that 20 candidates of the obtained genes were specifically or abundantly transcribed in the course of the fruiting body formation, suggesting that the obtained genes in this work play roles in fruiting body development in L. edodes.

Differential gene expression in filamentous cells of Ustilago maydis

When fungi interact with plants as pathogens or as symbionts, there are often changes in fungal cell morphology and nuclear state. This study establishes the use of cDNA microarrays to detect gene expression changes in Ustilago maydis cells that differ in structure and nuclear content. Categorizing differentially expressed genes on the basis of function indicated that U. maydis cell types vary most in the expression of genes related to metabolism. We also observed that more genes are up-regulated in the filamentous dikaryon than in the filamentous diploid, relative to non-pathogenic budding cells. Our comparison of pathogenic development indicated that the dikaryon is more virulent than the diploid. Other identified expression patterns suggest a cell-specific difference in nutrient acquisition, cell metabolism and signal transduction. The relevance of gene expression change to cell type biology is discussed.