Analyses of expressed sequence tags from the maize foliar pathogen Cercospora zeae-maydis identify novel genes expressed during vegetative, infectious, and reproductive growth
© Bluhm et al; licensee BioMed Central Ltd. 2008
Received: 17 January 2008
Accepted: 04 November 2008
Published: 04 November 2008
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.
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.
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.
The fungal genus Cercospora represents a large and diverse group of plant pathogens that are distributed worldwide and infect numerous host species. Individual species of Cercospora are usually host specific, but collectively they infect remarkably diverse hosts. More than 3,000 species of Cercospora have been named , and they often are classified according to host association, e.g., C. beticola infects sugar beet (Beta vulgaris), C. oryzae infects rice (Oryza sativa), and C. sorghi infects sorghum (Sorghum bicolor). Most plant-pathogenic species of Cercospora enter host leaves through stomata, a process facilitated in part by the ability of elongating germ tubes to sense nearby stomata and reorient their direction of growth accordingly . Upon reaching stomata, germ tubes differentiate into multilobed infection structures similar to appressoria, from which infectious hyphae penetrate mesophyll tissues. After a period of colonization, the fungus presumably adopts a necrotrophic growth habit, leading to the formation of expanding, necrotic lesions that coalesce in severe outbreaks, leading to a significant reduction in photosynthetic tissue, defoliation, and potentially premature death of the host plant. Reproduction and formation of secondary inocula occur in colonized tissue through the production of asexual spores (conidia) that infect neighboring plants after dispersal by wind and/or rain splash. Many diseases caused by Cercospora species occur periodically throughout the world as epidemics singly or as components of disease complexes [e.g., [3–5]], and for crops such as sugar beet, are major limitations for production . Additionally, the possibility that Cercospora pathogens influence the distribution of plant species in natural ecosystems is a plausible but largely unexplored hypothesis.
Cercospora zeae-maydis is a foliar pathogen causing gray leaf spot of maize. Substantial economic losses from this disease occur annually throughout the Western Hemisphere. First discovered in 1924 in Illinois , C. zeae-maydis did not become an important pathogen of maize until the 1980s; by the mid-1990s, the fungus caused significant losses throughout the corn belt of the U.S. and it is now the most devastating foliar pathogen of maize in much of the world . Colonization of leaves by the fungus causes distinctive rectangular lesions delineated by the major veins. When the incidence of infection is high before grain filling, the impaired photosynthetic capability of diseased leaves results in severe reductions in yield . Management of C. zeae-maydis is especially difficult because commercial hybrids of maize lack effective resistance to gray leaf spot  and the fungus can survive between growing seasons in plant debris . Exactly why C. zeae-maydis has ascended so rapidly as a pathogen of maize during the past two decades is not known, but speculation has linked the phenomenon to global climate change, the emergence of more virulent strains, and the increased practice of conservation tillage in maize production [8, 10, 11].
During pathogenesis, C. zeae-maydis and many other species of Cercospora produce the host non-specific phytotoxin cercosporin, a photosensitizing perylenequinone that causes lipid peroxidation and alters membrane permeability through the action of reactive oxygen species . Cercospora pathogens protect themselves against the toxic effects of cercosporin through the functions of CFP1, which encodes an ABC transporter required for secretion , and PDX1 (formerly SOR1), a gene involved in the biosynthesis of pyridoxine (vitamin B6), which quenches singlet oxygen produced during the interaction of cercosporin with membranes . Consistent with the production of many fungal secondary metabolites, cercosporin biosynthesis was recently demonstrated to result from the expression of genes organized in a cluster . However, a molecular understanding is lacking to explain how Cercospora species integrate diverse environmental inputs to regulate cercosporin biosynthesis and the extent to which regulation is conserved throughout the genus. Typically, cultures of C. zeae-maydis producing asexual spores (conidia) do not produce cercosporin, suggesting that fungal development and secondary metabolism are antagonistic at some level. In culture, cercosporin biosynthesis is repressed by the presence of preferred nitrogen sources ; presumably, this regulation is a component of global changes in gene expression during nitrogen metabolite repression resulting from the actions of a homolog of the nitrogen-responsive transcription factor areA . Additionally, the biosynthesis and activation of cercosporin require light , thus establishing an intriguing link between light and pathogenesis among Cercospora species.
The initial infection of maize leaves by C. zeae-maydis occurs in spring or early summer when propagules of the fungus that survived the winter in plant debris give rise to conidia that are dispersed onto leaves of young plants . During colonization of leaf tissue, the fungus produces stromata that give rise to erumpent conidiophores bearing conidia that serve as secondary inocula . Multiple cycles of secondary infection can occur when environmental conditions are favorable, leading to epidemic levels of infection. Somewhat surprisingly, C. zeae-maydis has not been demonstrated to reproduce sexually in laboratory conditions, and field populations appear to be largely clonal , although recent analyses of the distribution of mating type loci suggest the possibility of cryptic sex .
Despite the impact of C. zeae-maydis and other Cercospora species on global agriculture, very little information is available at the molecular level regarding how members of this genus regulate growth, development, and pathogenesis. The focus of this research was to generate a collection of expressed sequence tags (ESTs) from C. zeae-maydis, and to analyze their expression during defined stages of growth and development. To this end, we generated distinct cDNA libraries from vegetative cultures of C. zeae-maydis (vegetative libraries) as well as cultures producing conidia (sporulation libraries). Among 27,551 ESTs sequenced from both conditions, we identified 4619 unique sequences representing a broad range of molecular functions and biological processes. Of 4088 clusters containing two or more ESTs, 1436 were comprised of ESTs found exclusively in the sporulation libraries, whereas 1744 were unique to vegetative libraries. At least eight clusters encode putative photoreceptors and light-responsive genes, six are similar to genes regulating morphogenesis in other fungi, and 20 are implicated in host/pathogen interactions. The expression profiles of 15 clusters were characterized by real-time quantitative PCR, which largely confirmed their proposed roles in photoreception, conidiation, and pathogenesis. Furthermore, we identified Zephyr, a novel, highly transcribed member of the Ty3/Gypsy family of transposable elements. This research represents the first comprehensive EST sequencing project for C. zeae-maydis, and provides specific targets for subsequent studies in molecular genetics as well as a framework for future investigations into the evolution of pathogenesis among species of Cercospora and closely related genera.
Fungal strain and culture conditions
Wild-type C. zeae-maydis strain SCOH1-5, isolated from infected maize plants near South Charleston, Ohio in 1999, was used in all experiments. Cultures were maintained on V8 agar in constant darkness to provide conidia for inoculations. For library construction, the fungus was grown at 24°C on V8 agar, 0.2× potato dextrose agar (PDA; BD Biosciences, Sparks, MD), or 0.2× PDA supplemented with 10 mM ammonium phosphate. Cultures grown in constant light received 8–10 μE m-2s-1 of illumination. To facilitate collection of fungal tissue from agar plates, conidial suspensions were inoculated onto cellophane membranes placed on the surface of the medium. Maize inbred line B73, which is highly susceptible to infection by C. zeae-maydis, was grown in a greenhouse and inoculated with conidia (105/ml) with a fine-mist atomizer.
RNA isolation, cDNA library construction and sequencing
EST libraries constructed for this study
Insert size selection
Total unique ESTs
Passing to cluster
Small insert (100–700 bp)
Total number with BLAST hits2
Average diversity (passing clones)
Library quality was assessed first by randomly selecting 24 clones and amplifying the cDNA inserts by PCR with the primers M13-F (5'-GTAAAACGACGGCCAGT) andM13-R (5'-AGGAAACAGCTATGACCAT). The number of clones without inserts was determined and 384 clones for each library were picked, inoculated into 384-well plates (Nunc; Nalge Nunc International, Rochester, NY) and grown for 18 hr at 37°C. After amplification by rolling-circle amplification (RCA), the 5' and 3' ends of each insert were sequenced using vector-specific primers (FW: 5'- ATTTAGGTGACACTA TAGAA and RV 5' – TAATACGACTCACTATAGGG) and Big Dye chemistry (Applied Biosystems; Foster City, CA). For each insert, the clone identification information was retained for the 3' and 5' sequence reads. An additional sporulation library was generated from the cultures producing conidia (library CCAW) although the inserts were not size selected, and they were directionally ligated into the Sfi IA/B sites of the vector pDNR-Lib (BD Biosciences). In total, bidirectional sequencing of each library generated 9888 ESTs from CBYB, 3072 from CBYC, 3072 from CBYF, 2304 from CBYG and 9216 from CCAW. All sequences were deposited into the GenBank dbEST database; accession numbers are provided for each EST in the Additional Materials.
EST analyses and clustering
To trim vector sequences, common sequence patterns at the ends of ESTs were identified and removed. Clones were determined to lack inserts if ≥ 200 bases from the 5' end of the EST were identified as vector or if the insert was comprised of fewer than 100 bases of non-vector masked sequence. ESTs were then trimmed for quality with a sliding window trimmer (window = 11 bases). Once the average quality score in the window was below the quality threshold (Q15), the EST was split and the longest remaining sequence segment was retained as the trimmed EST. EST sequences with fewer than 100 bases of high-quality sequence were removed. ESTs were screened for the presence of polyA- or polyT-tails (which, if present, were deleted) and re-evaluated for length; ESTs with fewer than 100 bases were removed. ESTs consisting of more than 50% low-complexity sequence were removed from the final set of usable ESTs. If an EST required re-sequencing, the longest high-quality EST was retained. Sister ESTs or end-pair reads were categorized as follows: if one EST was insertless or a contaminant, then, by default, the sister EST was categorized as the same. However, each sister EST was treated separately for complexity and quality scores. Finally, EST sequences were compared against the GenBank nucleotide database by BLAST  to identify contaminants; undesirable ESTs such as those matching non-cellular sequences were removed.
For clustering, ESTs were evaluated with MALIGN , a kmer-based alignment tool that clusters ESTs based on sequence overlap (kmer = 16, seed length requirement = 32, alignment ID >= 98%). Clusters of ESTs were further merged based on sister reads using double linkage, which requires that two or more matching sister ESTs are in each cluster to be merged. EST clusters were then assembled using CAP3 to form consensus sequences. Clusters may contain more than one consensus sequence for various reasons (e.g., clone has long insert, clones are splice variants, consensus sequences are erroneously assembled). Cluster singlets are clusters of one EST, whereas CAP3 singlets are single ESTs that had joined a cluster but during cluster assembly were isolated into a separate consensus sequence. ESTs from each separate cDNA library were clustered and assembled separately, and subsequently the entire set of ESTs from all five cDNA libraries was clustered and assembled together with an external cDNA library (designated EXTA) obtained in an earlier study . A file containing all clusters, cluster singlets, and CAP3 singlets is available in the Additional Materials.
Annotation of ESTs with GO terms was done with Blast2Go . First, sequences were evaluated with BLASTx against the NCBI nr (non-redundant) database with an E-value threshold of 10-5. From a total of 7120 clusters, cluster singlets, and CAP3 singlets, 2526 sequences had no blast hits. Out of the remaining 4594 sequences, 2208 (48.1%) were categorized into different gene ontology (GO) classes at level three organization. For most clusters containing multiple consensus sequences or cluster singlets, a single Blast hit was selected for annotation. For a few clusters, some consensus sequences and/or cluster singlets corresponded to distinctly different genes, possibly due to overclustering, and thus were included in the final analysis. A file containing the annotation data for each cluster is provided in the Additional Materials.
Real-time quantitative PCR (qPCR)
Sequences examined by qPCR
Putative function or identity
Cluster or singlet ID
Primers for qPCR (5' to 3')
Glycine-rich cell wall protein
Conserved hypothetical protein
Blue light-induced gene 3 (bli-3)
Polyketide synthase (CTB1)
Regulator of xylanase activity (xlnR)
Regulator of cutinase activity (cf1α)
Transposable element (Zephyr)
Results and discussion
Construction and sequencing of cDNA libraries
During the interaction between C. zeae-maydis and maize, two key aspects of the disease cycle are colonization of host tissue and the production of conidia for secondary inocula. Although we are especially interested in identifying genes underlying host/pathogen interactions, a major drawback of constructing cDNA libraries from inoculated leaves is that a high percentage of ESTs are likely to correspond to plant rather than fungal genes. To circumvent this problem, we created cDNA libraries from sporulating cultures in early and late stages of conidiation (sporulation libraries) as well as vegetative cultures grown under a variety of conditions that support or repress cercosporin biosynthesis (vegetative libraries).
ESTs from the five cDNA libraries described above were combined for clustering analysis. Among the 4088 clusters, we identified 7120 consensus identification sequences that primarily reflected non-overlapping sequencing reads due to a high percentage of large inserts in the CBYB library (data not shown). Fifteen clusters (0.4%) contained three consensus identification sequences, and six (0.1%) contained four or more, reflecting a combination of alternative transcript splicing as well as erroneous grouping of sequences (overclustering). To determine the distribution of cluster sequences between the two sets of conditions, we performed a cluster overlap analysis with the 4088 clusters containing two or more ESTs. A total of 1744 clusters were comprised of ESTs found exclusively in the vegetative libraries, whereas 1436 clusters were comprised of ESTs found exclusively in the sporulation libraries. Only 908 (18%) of the clusters were comprised of ESTs found in both the vegetative and sporulation libraries, thus indicating that the conditions selected for library construction have a substantial impact on the transcriptome of C. zeae-maydis.
Sequence annotation and analysis
Fungal tissue from which cDNA libraries were constructed was obtained from cultures grown under a variety of conditions representing multiple stages of fungal development with the goal of obtaining a diverse collection of ESTs representing a range of molecular functions. ESTs were annotated according to Gene Ontology (GO)  guidelines with Blast2Go, a universal, web-based annotation application . To ensure the highest recovery of GO terms, we submitted all 7120 consensus identification sequences derived from the 4088 clusters and 531 singletons for Blast2Go analysis. In total, 2526 sequences had no blast hits with an E value ≤ 10-05. Out of these, 2515 sequences contained one or more predicted open reading frames of at least 100 amino acids. The sequences that have coding potential but do not share significant homology to deposited sequences could represent conserved genes that are not yet described in other fungi or genes that are unique to C. zeae-maydis. Of the 4594 sequences with BLAST hits, 2208 sequences (48.1%) were assigned GO terms. To eliminate over-representation of GO terms, a single BLAST hit was included in the final analysis for each cluster unless multiple consensus sequences for a given cluster corresponded to remarkably different proteins. In total, 1471 clusters were assigned GO terms.
Hydrolases and oxidoreductases comprised over 45% of the total number of molecular functions identified by GO analysis (Fig 3B). Hydrolases, which utilize water molecules to break chemical bonds, perform a broad range of functions in fungi, including the extracellular digestion of complex carbon sources such as cellulose and other components of plant cell walls. Oxidoreductases catalyze the transfer of electrons between molecules and in fungi are involved in primary and secondary metabolism (including cercosporin biosynthesis) as well as the detoxification of compounds such as reactive oxygen species, superoxide and hydrogen peroxide. Intriguingly, these same compounds are frequently associated with the oxidative burst component of plant defense . Although it is reasonable to propose that oxidoreductases of fungal foliar pathogens could be involved in detoxification of reactive oxygen species during pathogenesis, such a relationship has not been demonstrated. Sequences involved in signal transduction comprised 2% of the molecular functions identified (Fig 3B). Most of these genes were predicted to encode protein kinases, including 13 genes predicted to encode histidine kinases analogous to the two-component sensor histidine kinase family and three genes predicted to encode mitogen-activated protein kinases (MAPKs). The role of MAPKs in regulating morphology and virulence is well established in many fungi, including C. zeae-maydis [27, 28]. In filamentous fungi and yeasts, histidine kinases trigger phosphorelay signaling mechanisms that interact with various MAPKs to regulate growth, differentiation, and virulence [29, 30].
The vast majority of annotated sequences are predicted to encode intracellular proteins (Fig 3C). Considering the pathogenic lifestyle of C. zeae-maydis, we anticipated identifying a substantial number of secreted proteins, but found that limitations inherent to EST sequencing projects (e.g., 3' bias of sequence data, clones not corresponding to full-length transcripts) made predictions regarding secretion unreliable. However, nearly 1% of sequences were categorized by Blast2Go analysis as comprising external encapsulating structures, defined as any constituent of a structure that lies outside the plasma membrane and surrounds the entire cell . Unlike bacteria, filamentous fungi generally produce highly hydrophobic proteins (collectively referred to as hydrophobins) rather than polysaccharide capsules as a protective barrier against the environment. The extent to which hydrophobins are involved in pathogenesis among filamentous fungi is not clear, but in Magnaporthe grisea, a hydrophobin encoded by MPG1 is required for the efficient induction of appressoria, possibly by mediating aspects of surface recognition .
Consistent with many fungal EST projects, a substantial number of sequences could not be annotated due to either a lack of BLAST hits or hits to uncharacterized fungal sequences [e.g., [32, 33]]. Of the sequences with no BLAST hits, some fraction could be unique to C. zeae-maydis, whereas a significant percentage is likely to be too short to yield BLAST hits or correspond to untranslated regions of the mRNA (such as the 5' or 3' UTR). Of the sequences with BLAST hits, well over half could not be annotated due to a general lack of knowledge regarding the specific molecular functions of many fungal genes. For example, the genome of the closely related fungus Mycosphaerella graminicola is predicted to contain 11–12,000 genes, but to date, only ~30% have been annotated as to biological process, ~15% by cellular component, and ~40% by molecular function http://genome.jgi-psf.org/Mycgr1/Mycgr1.home.html.
Highly differentially expressed sequences
Consensus sequences highly enriched during either sporulation or vegetative growth
Cluster ID_consensus sequence
BLAST hit (gi #)1
No. of ESTs in libraries
Enriched during sporulation
No hits found
No hits found
No hits found
No hits found
Conserved hypothetical protein
No hits found
Conserved hypothetical protein
Conserved hypothetical protein
Conserved hypothetical protein
No hits found
ATP synthase protein 9
Ribosomal protein L30
ATP synthase subunit Atp18
60S ribosomal protein L36
60S ribosomal protein L29
Blue light induced gene-3
Cytochrome c oxidase
Conserved hypothetical protein
No hits found
Conserved hypothetical protein
Enriched during vegetative growth
No hits found
No hits found
No hits found
Photoreceptors and light-responsive genes
Putative photoreceptors and light-regulated genes found in EST libraries of Cercospora zeae-maydis
Cluster or singleton #
GenBank gi #1
BLAST e value
No. of ESTs in libraries
Clock-controlled gene 6
Glyceraldehyde 3-phosphate dehydrogenase (clock-controlled gene 7)
Clock-controlled gene 8
To investigate the transcriptional regulation of putative photoreceptors and light-regulated genes, we analyzed the expression of selected genes in response to light and growth medium. The sequence similar to members of the cryptochrome/photolyase family was more highly expressed in light than dark (Fig 4B) and, consistent with the distribution of corresponding ESTs (Table 4), the expression of the phytochrome-like gene did not appear to be affected by light (Fig 4B). Additionally, expression of the sequence similar to bli-3 appeared to be regulated primarily by growth medium (Fig 4B), thus explaining its enrichment in the sporulation rather than vegetative libraries.
Sequences implicated in the regulation of development
Sequences from EST libraries of Cercospora zeae-maydis implicated in fungal development
GenBank gi #1
Putative function (fungal homolog)2
No. of ESTs in libraries
Conidiophore development protein (hymA)
Regulator of conidiation (flbA)
Regulator of fruiting body formation (nosA)
Regulator of conidiation (nrc-2)
Regulator of development in response to nutrient availability (rcd1)
Pleiotropic developmental regulator (fst12)
Because conidia play a key role in the propagation of diseases caused by Cercospora species, we are particularly interested in identifying genes involved in the regulation of asexual development. In C. zeae-maydis and many other filamentous fungi, conidia are borne on specialized structures termed conidiophores (Fig. 1A, C, D). However, the morphological characteristics of conidia and conidiophores vary widely among fungi, often to the extent that the size and shape of conidia and/or conidiophores form a basis for taxonomic identification of genera or species. Given the structural complexity of conidiophores and conidia as well as the extent to which conidiation is regulated by environmental cues, asexual development presumably requires the coordinated expression of many genes. However, relatively little is known at the molecular level regarding how fungi regulate conidiation. Much of the existing knowledge is derived from model fungi, such as Aspergillus nidulans and Neurospora crassa, which are only distantly related to C. zeae-maydis.
We identified several consensus sequences from the EST libraries corresponding to genes known to regulate conidiation in filamentous fungi (Table 5). In A. nidulans, the regulator of G-protein signaling flbA is required for asexual sporulation , and the Zn(II)2Cys6 transcription factor encoded by nosA that is required for the induction of sexual development is also transcriptionally upregulated during asexual development . The exact molecular function of the protein encoded by hymA, also required for conidiophore formation in A. nidulans, is unknown . In Neurospora crassa, an insertional mutant that constitutively initiated, but failed to complete, conidial development arose from disruption of nrc-2, a gene encoding a serine-threonine protein kinase . Also, the putative green-light photoreceptor encoded by nop-1 regulates conidiation-specific gene expression in N. crassa, thus implicating the gene in fungal development . Because of the complexity of conidiophore and conidial development and their relatively poor evolutionary conservation among taxonomic classes of fungi, further characterization of candidate genes involved in asexual development in C. zeae-maydis will require functional characterization such as targeted disruption.
Pathogenesis-related sequences identified in EST libraries of Cercospora zeae-maydis
GenBank gi #1
No. of ESTs in libraries
Polyketide synthase involved in Cercosporin biosynthesis (CTB1)
Hybrid non-ribosomal peptide synthetase/polyketide synthase
Non-ribosomal peptide synthetase
Regulator of xylanase activity (xlnR)
Regulator of cutinase activity (cf1α)
Nitrogen response factor (areA)
Many fungal secondary metabolites, including cercosporin, are polyketide compounds formed by the head-to-tail condensation of acetate molecules as catalyzed by polyketide synthases. Recently, a gene cluster encoding a group of biosynthetic genes required for cercosporin biosynthesis was identified in Cercospora nicotianae, a foliar pathogen of tobacco . The cluster contains a polyketide synthase (CTB1), disruption of which abolishes cercosporin biosynthesis , as well as other coordinately regulated genes such as oxidoreductases hypothesized to catalyze specific steps in the biosynthesis of cercosporin . Consistent with established patterns of cercosporin biosynthesis in culture, the C. zeae-maydis homolog of CTB1 was induced by light on 0.2× PDA and was repressed by V8 agar irrespective of exposure to light (Fig 4B).
Although cercosporin is known to function as a virulence/pathogenicity factor in many Cercospora species, the dynamics of cercosporin biosynthesis during pathogenesis are largely unknown. To explore this question, we monitored CTB1 expression during colonization of leaf tissue. Although expression of CTB1 increased two-fold by 14 days after inoculation, it was somewhat surprising that expression of CTB1 changed little from 3–10 days after inoculation (Fig 5B). During this time, the fungus makes its initial penetration of mesophyll tissue and, as reflected by the visible development of lesions, commences necrotrophic growth. The absence of CTB1 induction during these stages of pathogenesis suggests that other virulence/pathogenicity factors may play a greater role in the initial colonization of leaf tissue.
Somewhat surprisingly, no other sequence similar to genes in the cercosporin biosynthesis (CTB) cluster was found among the ESTs obtained in this study. Among fungi that produce a given secondary metabolite, the underlying gene clusters are generally highly conserved, making it unlikely that C. zeae-maydis possesses a fundamentally different mechanism responsible for cercosporin biosynthesis. Rather, the most likely explanation for the absence of other CTB homologs from the EST dataset is that the cultures from which the vegetative libraries were produced represented a variety of growth conditions, not all of which supported cercosporin biosynthesis; therefore, the relative concentration of mRNAs corresponding to CTB genes was diluted. We hypothesize that more extensive sequencing of the vegetative libraries would lead to the identification of homologs of CTB genes such as those identified in C. nicotianae.
Identification and characterization of Zephyr, a novel transposable element
Among the ESTs highly represented in the vegetative library compared to the sporulation library, we identified a sequence highly similar to members of the Ty3/Gypsy family of long terminal repeat (LTR) transposons, including Grasshopper from Magnaporthe grisea , REAL from Alternaria alternata , and Skippy from Fusarium oxysporum . Members of the family typically contain two long, partially overlapping open reading frames encoding a protein similar to retroviral structural proteins and a poly protein containing protease, reverse transcriptase, RNaseH, and integrase domains . The retroelement identified in this study, designated Zephyr, is comprised of four clusters of 1749, 3664, 1028, and 1707 bp as well as four cluster singlets consisting of 110, 766, 745, and 252 bp. A conceptual translation of the 3664-bp cluster results in a protein of 1221 amino acids that corresponds to the poly protein of the element and is highly similar to Ty3/gypsy elements found in Magnaporthe grisea and other filamentous fungi, including the closely related fungus Mycosphaerella graminicola (data not shown). To date, only one other retroelement has been identified in C. zeae-maydis: Malazy, a degenerate, presumably non-coding member of the gypsy family  that shares substantial identity with Zephyr at the nucleotide level. Because the numerous premature stop codons found in Malazy are absent from Zephyr and EST evidence indicates Zephyr is an active element, we hypothesize that Malazy represents a defective/inactivated descendent of Zephyr.
Transposition of retroelements in fungi can be induced by a variety of biotic and abiotic stresses as well as morphological changes such as sexual reproduction. However, activation of transposable elements is relatively rare during the normal growth and development of most organisms, including filamentous fungi . Therefore, it is somewhat surprising that ESTs representing the polyprotein-encoding region of Zephyr are highly enriched in the vegetative libraries (34 ESTs in vegetative libraries compared to two in sporulation libraries). To further characterize the regulation of Zephyr, we profiled its expression during growth of C. zeae-maydis under a variety of environmental conditions. When evaluated by qPCR, cycle threshold values were low (< 20), thus indicating high levels of expression. Consistent with the distribution of ESTs in vegetative and sporulation libraries, the pol polyprotein-encoding region of Zephyr was expressed nearly 4-fold greater during growth on 0.2× PDA in light than on V8 agar in darkness. These results suggest that Zephyr is an actively transcribed element that is regulated by growth medium and possibly by light. Although EST and qPCR data indicate that Zephyr is highly expressed, further studies will be required to verify the transposition of the element.
Currently, little is known regarding the molecular mechanisms controlling the activity of transposable elements in fungi. Transposable elements have been implicated as a driving force behind genetic diversity; their activation in response to environmental stress is hypothesized to be a mechanism of adaptation, and consequently, genomic evolution . Because Cercospora is believed to be a largely asexual genus, transposition of elements such as Zephyr could be a driving force behind the remarkably high level of host-specific speciation that has evolved among Cercospora species.
By generating ESTs from vegetative and sporulating cultures of C. zeae-maydis, we identified novel genes involved in a wide range of biological processes. Functional annotation and expression profiling implicated subsets of genes in pathogenesis and conidiation. Consistent with the crucial role light plays in host-pathogen interactions between C. zeae-maydis and maize, we identified a large number of photoreceptors and light-regulated genes, plus Zephyr, a novel, highly expressed transposable element. We conclude that light plays a key role in the dichotomy between vegetative and reproductive growth in C. zeae-maydis and that future characterization of the underlying molecular mechanisms will contribute significantly to the fundamental understanding of how fungi respond to light.
We thank Corie Shaner and Kaila Zink for technical assistance and Charles Crane for assistance with bioinformatic analyses. This report constitutes ARP 2008-18287 of the Purdue University Agriculture Experiment Station. This research was supported by USDA CRIS project 3602-22000-013-00D and the Community Sequencing Program of the Joint Genome Institute. This work was performed under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program and by the University of California, Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48, Lawrence Berkeley National Laboratory under contract No. DE-AC03-76SF00098 and Los Alamos National Laboratory under contract No. W-7405-ENG-36.
- Pollack FG: An annotated compilation of Cercospora names. Mycolog Mem. 1987, 12: 1-212.Google Scholar
- Beckman PM, Payne GA: External growth, penetration, and development of Cercospora zeae-maydis and lesion development in corn. Phytopathology. 1982, 73: 286-289. 10.1094/Phyto-73-286.View ArticleGoogle Scholar
- Berger RD: Early blight of celery: analysis of disease spread in Florida. Phytopathology. 1973, 63: 1161-65.View ArticleGoogle Scholar
- de Nazareno NRX, Madden LV, Lipps PE: Characterization of gray leaf spot epidemics of maize. J Plant Dis Prot. 1993, 100: 410-425.Google Scholar
- Nelson SC, Campbell CL: Disease progress, defoliation, and spatial pattern in a multiple-pathogen disease complex on white clover. Phytopathology. 1993, 83: 419-429. 10.1094/Phyto-83-419.View ArticleGoogle Scholar
- Cooke DA, Scott RK: The Sugar Beet Crop. 1993, London: Chapman and HallView ArticleGoogle Scholar
- Tehon LR, Daniels E: Notes on the parasitic fungi of Illinois: II. Mycologia. 1925, 17: 240-249. 10.2307/3753890.View ArticleGoogle Scholar
- Ward JMJ, Stromberg EL, Nowell DC, Nutter FW: Gray leaf spot: a disease of global importance in maize production. Plant Dis. 1999, 83: 884-895. 10.1094/PDIS.19126.96.36.1994.View ArticleGoogle Scholar
- de Nazareno NRX, Lipps PE, Madden LV: Survival of Cercospora zeae-maydis in corn residues in Ohio. Plant Dis. 1992, 76: 560-564.View ArticleGoogle Scholar
- Latterell FM, Rossi AE: Gray leaf spot of maize: A disease on the move. Plant Dis. 1983, 67: 842-847. 10.1094/PD-67-842.View ArticleGoogle Scholar
- Rosenzweig C, Iglesias A, Yang XB, Epstein PR, Chivan E: Climate change and U.S. agriculture: The impacts of warming and extreme weather events on productivity, plant diseases, and pests. 2000, Boston: Center for Health and the Global Environment, Harvard Medical SchoolGoogle Scholar
- Daub ME, Ehrenshaft M: The photoactivated Cercospora toxin cercosporin: Contributions to plant disease and fundamental biology. Annu Rev Phytopathol. 2000, 38: 461-490. 10.1146/annurev.phyto.38.1.461.PubMedView ArticleGoogle Scholar
- Callahan TM, Rose MS, Meade MJ, Ehrenshaft M, Upchurch RG: CFP, the putative cercosporin transporter of Cercospora kikuchii, is required for wild type cercosporin production, resistance, and virulence on soybean. Mol Plant-Microbe Interact. 1999, 12: 901-910. 10.1094/MPMI.19188.8.131.521.PubMedView ArticleGoogle Scholar
- Ehrenshaft M, Bilski P, Li M, Chignell CF, Daub ME: A highly conserved sequence is a novel gene involved in de novo vitamin B6 biosynthesis. Proc Natl Acad Sci USA. 1999, 96: 9374-9378. 10.1073/pnas.96.16.9374.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen H, Lee MH, Daub ME, Chung KR: Molecular analysis of the cercosporin biosynthetic gene cluster in Cercospora nicotianae. Mol Microbiology. 2007, 64: 755-770. 10.1111/j.1365-2958.2007.05689.x.View ArticleGoogle Scholar
- Shim WB, Dunkle LD: Identification of genes expressed during cercosporin biosynthesis in Cercospora zeae-maydis. Physiol Mol Plant Pathol. 2002, 61: 237-248. 10.1006/pmpp.2002.0437.View ArticleGoogle Scholar
- Kudla B, Caddick MX, Langdon T, Martinez-Rossi NM, Bennett CF, Sibley S, Davies RW, Arst HN: The regulatory gene areA mediating nitrogen metabolite repression in Aspergillus nidulans. Mutations affecting specificity of gene activation alter a loop residue of a putative zinc finger. EMBO J. 1990, 9: 1355-1364.PubMedPubMed CentralGoogle Scholar
- Ehrenshaft M, Upchurch RG: Isolation of light-enhanced complementary DNAs of Cercospora kikuchii. Appl Environ Microbiol. 1991, 57: 2671-2676.PubMedPubMed CentralGoogle Scholar
- Wang J, Levy M, Dunkle LD: Sibling species of Cercospora associated with gray leaf spot of maize. Phytopathology. 1998, 88: 1269-1275. 10.1094/PHYTO.19184.108.40.2069.PubMedView ArticleGoogle Scholar
- Groenewald M, Groenewald JZ, Harrington TC, Abeln ECA, Crous PW: Mating type gene analysis in apparently asexual Cercospora species is suggestive of cryptic sex. Fungal Genet Biol. 2006, 43: 813-825. 10.1016/j.fgb.2006.05.008.PubMedView ArticleGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor, New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Wheeler WC, Gladstein DS: MALIGN: a multiple sequence alignment program. J Hered. 1994, 85: 417-418.Google Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMedPubMed CentralView ArticleGoogle Scholar
- Low PS, Merida JR: The oxidative burst in plant defense: function and signal transduction. Physiol Plant. 1996, 96: 533-542. 10.1111/j.1399-3054.1996.tb00469.x.View ArticleGoogle Scholar
- Roman E, Arana DM, Nombela C, Alonso-Monge R, Pla J: MAP kinase pathways as regulators of fungal virulence. Trends Microbiol. 2007, 15: 181-190. 10.1016/j.tim.2007.02.001.PubMedView ArticleGoogle Scholar
- Shim WB, Dunkle LD: CZK3, a MAP kinase kinase kinase homolog in Cercospora zeae-maydis, regulates cercosporin biosynthesis, fungal development, and pathogenesis. Mol Plant Microbe Int. 2003, 16: 760-768. 10.1094/MPMI.2003.16.9.760.View ArticleGoogle Scholar
- Santos JL, Shiozaki K: Fungal histidine kinases. Science's STKE. 2001, 98: RE1-Google Scholar
- Vargas-Perez I, Sanchez O, Kawasaki L, Georgellis D, Aguirre J: Response regulators SrrA and SskA are central components of a phosphorelay system involved in stress signal transduction and asexual sporulation in Aspergillus nidulans. Eukary Cell. 2007, 6: 1570-1583. 10.1128/EC.00085-07.View ArticleGoogle Scholar
- Beckerman JL, Ebbole DJ: MPG1, a gene encoding a fungal hydrophobin of Magnaporthe grisea, is involved in surface recognition. Mol Plant-Microbe Interact. 1996, 9: 450-456.PubMedView ArticleGoogle Scholar
- Nugent KG, Choffe K, Saville BJ: Gene expression during Ustilago maydis diploid filamentous growth: EST library creation and analyses. Fungal Genet Biol. 2004, 41: 349-360. 10.1016/j.fgb.2003.11.006.PubMedView ArticleGoogle Scholar
- Yockteng R, Marthey S, Chiapello H, Gendrault A, Hood ME, Rodolphe F, Devier B, Wincker P, Dossat C, Giraud T: Expressed sequences tags of the anther smut fungus, Microbotryum violaceum, identify mating and pathogenicity genes. BMC Genomics. 2007, 8: 272-10.1186/1471-2164-8-272.PubMedPubMed CentralView ArticleGoogle Scholar
- McClung JK, Jupe ER, Liu X-T, Dell'Orco RT: Prohibitin: Potential role in senescence, development, and tumor suppression. Exper Gerontol. 1995, 30: 99-124. 10.1016/0531-5565(94)00069-7.View ArticleGoogle Scholar
- Shaffer PM, Arst HN, Estberg L, Fernando L, Ly T, Sitter M: An asparaginase of Aspergillus nidulans is subject to oxygen repression in addition to nitrogen metabolite repression. Mol Gen Genet. 1988, 212: 337-341. 10.1007/BF00334704.PubMedView ArticleGoogle Scholar
- Herrera-Estrella A, Horwitz BA: Looking through the eyes of fungi: molecular genetics of photoreception. Mol Microbiol. 2007, 64: 5-15. 10.1111/j.1365-2958.2007.05632.x.PubMedView ArticleGoogle Scholar
- Blumenstein A, Vienken K, Tasler R, Purschwitz J, Veith D, Frankenberg-Dinkel N, Fischer R: The Aspergillus nidulans phytochrome FphA represses sexual development in red light. Curr Biol. 2005, 15: 1833-1838. 10.1016/j.cub.2005.08.061.PubMedView ArticleGoogle Scholar
- Yakura M, Ozoe F, Ishida H, Nakagawa T, Tanaka K, Matsuda H, Kawamukai M: Zds1, a novel gene encoding an ortholog of Zds1 and Zds2, controls sexual differentiation, cell wall integrity and cell morphology in fission yeast. Genetics. 2006, 172: 811-825. 10.1534/genetics.105.050906.PubMedPubMed CentralView ArticleGoogle Scholar
- Okazaki N, Okazaki K, Watanabe Y, Kato-Hayashi M, Yamamoto M, Okayama H: Novel factor highly conserved among eukaryotes controls sexual development in fission yeast. Mol Cell Biol. 1998, 18: 887-895.PubMedPubMed CentralView ArticleGoogle Scholar
- Goldar MM, Jeong HT, Tanaka K, Matsuda H, Kawamukai M: Moc3, a novel Zn finger type protein involved in sexual development, ascus formation, and stress response of Schizosaccharomyces pombe. Curr Genet. 2005, 48: 345-355. 10.1007/s00294-005-0028-z.PubMedView ArticleGoogle Scholar
- Herskowitz I: MAP kinase pathways in yeast: for mating and more. Cell. 1995, 80: 187-197. 10.1016/0092-8674(95)90402-6.PubMedView ArticleGoogle Scholar
- Deng FY, Allen TD, Hillman BI, Nuss DL: Comparative analysis of alterations in host phenotype and transcript accumulation following hypovirus and mycoreovirus infections of the chestnut blight fungus Cryphonectria parasitica. Eukary Cell. 2007, 6: 1286-1298. 10.1128/EC.00166-07.View ArticleGoogle Scholar
- Hoi JWS, Herbert C, Bacha N, O'Connell R, Lafitte C, Borderies G, Rossignol M, Rouge P, Dumas B: Regulation and role of a STE12-like transcription factor from the plant pathogen Colletotrichum lindemuthianum. Mol Microbiol. 2007, 64: 68-82. 10.1111/j.1365-2958.2007.05639.x.View ArticleGoogle Scholar
- Morita H, Hatamoto O, Masuda T, Sato T, Takeuchi M: Function analysis of a steA homolog in Aspergillus oryzae. Fungal Genet Biol. 2007, 44: 330-338. 10.1016/j.fgb.2006.10.009.PubMedView ArticleGoogle Scholar
- Yu JH, Wieser J, Adams TH: The Aspergillus FlbA RGS domain protein antagonizes G protein signaling to block proliferation and allow development. EMBO J. 1996, 15: 5184-5190.PubMedPubMed CentralGoogle Scholar
- Vienken K, Rishcher R: The Zn(II)2Cys6 putative transcription factor NosA controls fruiting body formation in Aspergillus nidulans. Mol Microbiol. 2006, 61: 544-554. 10.1111/j.1365-2958.2006.05257.x.PubMedView ArticleGoogle Scholar
- Karos M, Fischer R: Molecular characterization of HymA, an evolutionarily conserved and highly expressed protein of Aspergillus nidulans. Mol Gen Genet. 1999, 260: 510-521. 10.1007/s004380050924.PubMedView ArticleGoogle Scholar
- Kothe GO, Free SJ: The isolation and characterization of nrc-1 and nrc-2, two genes encoding protein kinases that control growth and development in Neurospora crassa. Genetics. 1998, 149: 117-130.PubMedPubMed CentralGoogle Scholar
- Bieszke JA, Li L, Borkovich KA: The fungal opsin gene nop-1 is negatively regulated by a component of the blue light sensing pathway and influences conidiation-specific gene expression in Neurospora crassa. Curr Genet. 2007, 52: 149-157. 10.1007/s00294-007-0148-8.PubMedView ArticleGoogle Scholar
- Scazzocchio C: The fungal GATA factors. Curr Opin Microbiol. 2000, 3: 126-131. 10.1016/S1369-5274(00)00063-1.PubMedView ArticleGoogle Scholar
- Li D, Kolattukudy PE: Cloning of cutinase transcription factor 1, a transactivating protein containing Cys6Zn2 binuclear cluster DNA-binding motif. J Biol Chem. 1997, 272: 12462-12467. 10.1074/jbc.272.19.12462.PubMedView ArticleGoogle Scholar
- van Peij NNME, Visser J, de Graaff LH: Isolation and analysis of xlnR, encoding a transcriptional activator co-ordinating xylanolytic expression in Aspergillus niger. Mol Microbiol. 1998, 27: 131-142. 10.1046/j.1365-2958.1998.00666.x.PubMedView ArticleGoogle Scholar
- Choquer M, Lahey KA, Chen HL, Cao L, Ueng PP, Daub ME, Chung KR: The CTB1 gene encoding a fungal polyketide synthase is required for cercosporin toxin biosynthesis and fungal virulence in Cercospora nicotianae. Mol Plant Microbe Interact. 2005, 18: 468-476. 10.1094/MPMI-18-0468.PubMedView ArticleGoogle Scholar
- Chen HQ, Lee MH, Chung KR: Functional characterization of three genes encoding putative oxidoreductases required for cercosporin toxin biosynthesis in the fungus Cercospora nicotianae. Microbiol. 2007, 153: 2781-2790. 10.1099/mic.0.2007/007294-0.View ArticleGoogle Scholar
- Dobinson KF, Harris RE, Hamer JE: Grasshopper, a long terminal repeat (LTR) retroelement in the phytopathogenic fungus Magnaporthe grisea. Mol Plant Microbe Interact. 1993, 6: 114-126.PubMedView ArticleGoogle Scholar
- Kaneko I, Tanaka A, Tsuge T: REAL, an LTR retrotransposon from the plant pathogenic fungus Alternaria alternata. Mol Gen Genet. 2000, 263: 625-634. 10.1007/s004380051210.PubMedView ArticleGoogle Scholar
- Anaya N, Roncero MI: Skippy, a retrotransposon from the fungal plant pathogen Fusarium oxysporum. Mol Gen Genet. 1995, 249: 637-647. 10.1007/BF00418033.PubMedView ArticleGoogle Scholar
- Lesage P, Todeshini AL: Happy together: the life and times of Ty retrotransposons and their hosts. Cytogenet Genome Res. 2005, 110: 70-90. 10.1159/000084940.PubMedView ArticleGoogle Scholar
- Shim WB, Dunkle LD: Malazy, a degenerate, species-specific transposable element in Cercospora zeae-maydis. Mycologia. 2005, 97: 349-355. 10.3852/mycologia.97.2.349.PubMedView ArticleGoogle Scholar
- Daboussi MJ, Capy P: Transposable elements in filamentous fungi. Annu Rev Microbiol. 2003, 57: 275-299. 10.1146/annurev.micro.57.030502.091029.PubMedView ArticleGoogle Scholar
- McClintock B: The significance of responses of the genome to challenge. Science. 1984, 226: 792-801. 10.1126/science.15739260.PubMedView ArticleGoogle Scholar
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.