HAG Effector Evolution in Pyricularia Species and Plant Cell Death Suppression by HAG4.

Seventy host-adapted gene (HAG) effector family members from Pyricularia species are found in P. oryzae and three closely related species (isolates LS and 18-2 from an unknown Pyricularia sp., P. grisea, and P. pennisetigena) that share at least eight orthologous HAG family members with P. oryzae. The genome sequence of a more distantly related species, P. penniseti, lacks HAG genes, suggesting a time frame for the origin of the gene family in the genus. In P. oryzae, HAG4 is uniquely found in the genetic lineage that contains populations adapted to Setaria and Oryza hosts. We find a nearly identical HAG4 allele in a P. grisea isolate, suggesting transfer of HAG4 from P. grisea to P. oryzae. HAG4 encodes a suppressor of plant cell death. Yeast two-hybrid screens with several HAG genes independently identify common interacting clones from a rice complementary DNA library, suggesting conservation of protein surface motifs between HAG homologs with as little as 40% protein sequence identity. HAG family orthologs have diverged rapidly and HAG15 orthologs display unusually high rates of sequence divergence compared with adjacent genes suggesting gene-specific accelerated divergence. The sequence diversity of the HAG homologs in Pyricularia species provides a resource for examining mechanisms of gene family evolution and the relationship to structural and functional evolution of HAG effector family activity. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Evolution and Regulation of a Large Effector Family of Pyricularia oryzae.

Plant pathogen effectors play important roles in parasitism, including countering plant immunity. However, investigations of the emergence and diversification of fungal effectors across host-adapted populations has been limited. We previously identified a gene encoding a suppressor of plant cell death in Pyricularia oryzae (syn. Magnaporthe oryzae). Here, we report the gene is one of a 21-member gene family and we characterize sequence diversity in different populations. Within the rice pathogen population, nucleotide diversity is low, however; the majority of gene family members display presence-absence polymorphism or other null alleles. Gene family allelic diversity is greater between host-adapted populations and, thus, we named them host-adapted genes (HAGs). Multiple copies of HAGs were found in some genome assemblies and sequence divergence between the alleles in two cases suggested they were the result of repeat-induced point mutagenesis. Transfer of family members between populations and novel HAG haplotypes resulting from apparent recombination were observed. HAG family transcripts were induced in planta and a subset of HAGs are dependent on a key regulator of pathogenesis, PMK1. We also found differential intron splicing for some HAGs that would prevent ex planta protein expression. For some genes, spliced transcript was expressed in antiphase with an overlapping antisense transcript. Characterization of HAG expression patterns and allelic diversity reveal novel mechanisms for HAG regulation and mechanisms generating sequence diversity and novel allele combinations. This evidence of strong in planta-specific expression and selection operating on the HAG family is suggestive of a role in parasitism.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Emergence of a hybrid PKS-NRPS secondary metabolite cluster in a clonal population of the rice blast fungus Magnaporthe oryzae.

Secondary metabolites (SMs) are crucial for fungi and vary in function from beneficial antibiotics to pathogenicity factors. To generate diversified SMs that enable different functions, SM-coding regions rapidly evolve in fungal genomes. However, the driving force and genetic mechanism of fungal SM diversification in the context of host-pathogen interactions remain largely unknown. Previously, we grouped field populations of the rice blast fungus Magnaporthe oryzae (syn: Pyricularia oryzae) into three major globally distributed clades based on population genomic analyses. Here, we characterize a recent duplication of an avirulent gene-containing SM cluster, ACE1, in a clonal M. oryzae population (Clade 2). We demonstrate that the ACE1 cluster is specifically duplicated in Clade 2, a dominant clade in indica rice-growing areas. With long-read sequencing, we obtained chromosome-level genome sequences of four Clade 2 isolates, which displayed differences in genomic organization of the ACE1 duplication process. Comparative genomic analyses suggested that the original ACE1 cluster experienced frequent rearrangement in Clade 2 isolates and revealed that the new ACE1 cluster is located in a newly formed and transposable element-rich region. Taken together, these results highlight the frequent mutation and expansion of an avirulent gene-containing SM cluster through transposable element-mediated whole-cluster duplication in the context of host-pathogen interactions.

Retromer Is Essential for Autophagy-Dependent Plant Infection by the Rice Blast Fungus.

The retromer mediates protein trafficking through recycling cargo from endosomes to the trans-Golgi network in eukaryotes. However, the role of such trafficking events during pathogen-host interaction remains unclear. Here, we report that the cargo-recognition complex (MoVps35, MoVps26 and MoVps29) of the retromer is essential for appressorium-mediated host penetration by Magnaporthe oryzae, the causal pathogen of the blast disease in rice. Loss of retromer function blocked glycogen distribution and turnover of lipid bodies, delayed nuclear degeneration and reduced turgor during appressorial development. Cytological observation revealed dynamic MoVps35-GFP foci co-localized with autophagy-related protein RFP-MoAtg8 at the periphery of autolysosomes. Furthermore, RFP-MoAtg8 interacted with MoVps35-GFP in vivo, RFP-MoAtg8 was mislocalized to the vacuole and failed to recycle from the autolysosome in the absence of the retromer function, leading to impaired biogenesis of autophagosomes. We therefore conclude that retromer is essential for autophagy-dependent plant infection by the rice blast fungus.

WD40-repeat protein MoCreC is essential for carbon repression and is involved in conidiation, growth and pathogenicity of Magnaporthe oryzae.

Carbon catabolite repression (CCR) is a common regulatory mechanism used by microorganisms to prioritize use of a preferred carbon source (usually glucose). The CreC WD40-repeat protein is a major component of the CCR pathway in Aspergillus nidulans. To clarify the function of the CreC ortholog from Magnaporthe oryzae in regulating gene expression important for pathogenesis, MoCreC was identified and genetically characterized. The vegetative growth rate of the MoCreC deletion mutant on various carbon sources was reduced. The MoCreC mutant produced fewer conidia and with about 60% of conidia having septation defects. Appressorium formation was impaired in the MoCreC mutant. Although some appressoria of the mutant could penetrate the leaf surface successfully, the efficiency of penetration and invasive growth of infection hyphae was reduced, resulting in attenuated virulence toward host plants. The CCR was defective as the mutant was more sensitive to allyl alcohol in the presence of glucose, and 2-deoxyglucose was unable to fully repress utilization of secondary carbon sources. qRT-PCR results indicated that the genes encoding cell wall degradation enzymes, such as ß-glucosidase, feruloyl esterase and exoglucanase, were upregulated in MoCreC mutant. Taken together, we conclude that MoCreC is a major regulator of CCR and plays significant roles in regulating growth, conidiation, and pathogenicity of M. oryzae.

Neurospora crassa ASM-1 complements the conidiation defect in a stuA mutant of Aspergillus nidulans.

Aspergillus nidulans StuA and Neurospora crassa ASM-1 are orthologous APSES (ASM-1, PHD1, SOK2, Efg1, StuA) transcription factors conserved across a diverse group of fungi. StuA and ASM-1 have roles in asexual (conidiation) and sexual (ascospore formation) development in both organisms. To address the hypothesis that the last common ancestor of these diverse fungi regulated conidiation with similar genes, asm-1 was introduced into the stuA1 mutant of A. nidulans. Expression of asm-1 complemented defective conidiophore morphology and restored conidia production to wild type levels in stuA1. Expression of asm-1 in the stuA1 strain did not rescue the defect in sexual development. When the conidiation regulator AbaA was tagged at its C-terminus with GFP in A. nidulans, it localized to nuclei in phialides. When expressed in the stuA1 mutant, AbaA::GFP localized to nuclei in conidiophores but no longer was confined to phialides, suggesting that expression of AbaA in specific cell types of the conidiophore was conditioned by StuA. Our data suggest that the function in conidiation of StuA and ASM-1 is conserved and support the view that, despite the great morphological and ontogenic diversity of their condiphores, the last common ancestor of A. nidulans and N. crassa produced an ortholog of StuA that was involved in conidiophore development.

acon-3, the Neurospora crassa ortholog of the developmental modifier, medA, complements the conidiation defect of the Aspergillus nidulans mutant.

Aspergillus nidulans and Neurospora crassa are ascomycetes that produce asexual spores through morphologically distinct processes. MedA, a protein with unknown function, is required for normal asexual and sexual development in A. nidulans. We determined that the N. crassa ortholog of medA is acon-3, a gene required for early conidiophore development and female fertility. To test hypotheses about the evolutionary origins of asexual development in distinct fungal lineages it is important to understand the degree of conservation of developmental regulators. The amino acid sequences of A. nidulans MedA and N. crassa ACON-3 shared 37% identity and 51% similarity. acon-3 is induced at late time points of conidiation. In contrast, medA is constitutively expressed and MedA protein localizes to nuclei in all tissue types. Nonetheless, expression of acon-3 using its native promoter complemented the conidiation defects of the A. nidulans ΔmedA and medA15 mutants. We conclude that the biochemical activity of the medA orthologs is conserved for conidiation.

The genome sequence of the rice blast fungus Magnaporthe grisea.

Magnaporthe grisea is the most destructive pathogen of rice worldwide and the principal model organism for elucidating the molecular basis of fungal disease of plants. Here, we report the draft sequence of the M. grisea genome. Analysis of the gene set provides an insight into the adaptations required by a fungus to cause disease. The genome encodes a large and diverse set of secreted proteins, including those defined by unusual carbohydrate-binding domains. This fungus also possesses an expanded family of G-protein-coupled receptors, several new virulence-associated genes and large suites of enzymes involved in secondary metabolism. Consistent with a role in fungal pathogenesis, the expression of several of these genes is upregulated during the early stages of infection-related development. The M. grisea genome has been subject to invasion and proliferation of active transposable elements, reflecting the clonal nature of this fungus imposed by widespread rice cultivation.

Structure and dynamics of the microbial communities underlying the carboxylate platform for biofuel production.

The carboxylate platform utilizes a mixed microbial community to convert lignocellulosic biomass into chemicals and fuels. While much of the platform is well understood, little is known about its microbiology. Mesophilic (40 degrees C) and thermophilic (55 degrees C) fermentations employing a sorghum feedstock and marine sediment inoculum were profiled using 16S rRNA tag-pyrosequencing over the course of a 30-day incubation. The contrasting fermentation temperatures converted similar amounts of biomass, but the mesophilic community was significantly more productive, and the two temperatures differed significantly with respect to propionic and butyric acid production. Pyrotag sequencing revealed the presence of dynamic communities that responded rapidly to temperature and changed substantially over time. Both temperatures were dominated by bacteria resembling Clostridia, but they shared few taxa in common. The species-rich mesophilic community harbored a variety of Bacteroidetes, Actinobacteria, and gamma-Proteobacteria, whereas the thermophilic community was composed mainly of Clostridia and Bacilli. Despite differences in composition and productivity, similar patterns of functional class dynamics were observed. Over time, organisms resembling known cellulose degraders decreased in abundance, while organisms resembling known xylose degraders increased. Improved understanding of the carboxylate platform's microbiology will help refine platform performance and contribute to our growing knowledge regarding biomass conversion and biofuel production processes.

Gene Ontology annotation of the rice blast fungus, Magnaporthe oryzae.

BACKGROUND: Magnaporthe oryzae, the causal agent of blast disease of rice, is the most destructive disease of rice worldwide. The genome of this fungal pathogen has been sequenced and an automated annotation has recently been updated to Version 6 http://www.broad.mit.edu/annotation/genome/magnaporthe_grisea/MultiDownloads.html. However, a comprehensive manual curation remains to be performed. Gene Ontology (GO) annotation is a valuable means of assigning functional information using standardized vocabulary. We report an overview of the GO annotation for Version 5 of M. oryzae genome assembly. METHODS: A similarity-based (i.e., computational) GO annotation with manual review was conducted, which was then integrated with a literature-based GO annotation with computational assistance. For similarity-based GO annotation a stringent reciprocal best hits method was used to identify similarity between predicted proteins of M. oryzae and GO proteins from multiple organisms with published associations to GO terms. Significant alignment pairs were manually reviewed. Functional assignments were further cross-validated with manually reviewed data, conserved domains, or data determined by wet lab experiments. Additionally, biological appropriateness of the functional assignments was manually checked. RESULTS: In total, 6,286 proteins received GO term assignment via the homology-based annotation, including 2,870 hypothetical proteins. Literature-based experimental evidence, such as microarray, MPSS, T-DNA insertion mutation, or gene knockout mutation, resulted in 2,810 proteins being annotated with GO terms. Of these, 1,673 proteins were annotated with new terms developed for Plant-Associated Microbe Gene Ontology (PAMGO). In addition, 67 experiment-determined secreted proteins were annotated with PAMGO terms. Integration of the two data sets resulted in 7,412 proteins (57%) being annotated with 1,957 distinct and specific GO terms. Unannotated proteins were assigned to the 3 root terms. The Version 5 GO annotation is publically queryable via the GO site http://amigo.geneontology.org/cgi-bin/amigo/go.cgi. Additionally, the genome of M. oryzae is constantly being refined and updated as new information is incorporated. For the latest GO annotation of Version 6 genome, please visit our website http://scotland.fgl.ncsu.edu/smeng/GoAnnotationMagnaporthegrisea.html. The preliminary GO annotation of Version 6 genome is placed at a local MySql database that is publically queryable via a user-friendly interface Adhoc Query System. CONCLUSION: Our analysis provides comprehensive and robust GO annotations of the M. oryzae genome assemblies that will be solid foundations for further functional interrogation of M. oryzae.

Magnaporthe as a model for understanding host-pathogen interactions.

The rice blast pathosystem has been the subject of intense interest in part because of the importance of the disease to world agriculture, but also because both Magnaporthe oryzae and its host are amenable to advanced experimental approaches. The goal of this review is to provide an overview of the system and to point out recent significant studies that update our understanding of the biology of M. oryzae. The genome sequence of M. oryzae has provided insight into how genome structure and pathogen population genetic variability has been shaped by transposable elements. The sequence allows systematic approaches to long-standing areas of investigation, including pathogen development and the molecular basis of compatible and incompatible interactions with its host. Rice blast provides an integrated system to illustrate most of the important concepts governing fungal/plant interactions and serves as an excellent starting point for gaining a broad perspective of issues in plant pathology.

Magnaporthe oryzae CK2 Accumulates in Nuclei, Nucleoli, at Septal Pores and Forms a Large Ring Structure in Appressoria, and Is Involved in Rice Blast Pathogenesis.

Magnaporthe oryzae (Mo) is a model pathogen causing rice blast resulting in yield and economic losses world-wide. CK2 is a constitutively active, serine/threonine kinase in eukaryotes, having a wide array of known substrates, and involved in many cellular processes. We investigated the localization and role of MoCK2 during growth and infection. BLAST search for MoCK2 components and targeted deletion of subunits was combined with protein-GFP fusions to investigate localization. We found one CKa and two CKb subunits of the CK2 holoenzyme. Deletion of the catalytic subunit CKa was not possible and might indicate that such deletions are lethal. The CKb subunits could be deleted but they were both necessary for normal growth and pathogenicity. Localization studies showed that the CK2 holoenzyme needed to be intact for normal localization at septal pores and at appressorium penetration pores. Nuclear localization of CKa was however not dependent on the intact CK2 holoenzyme. In appressoria, CK2 formed a large ring perpendicular to the penetration pore and the ring formation was dependent on the presence of all CK2 subunits. The effects on growth and pathogenicity of deletion of the b subunits combined with the localization indicate that CK2 can have important regulatory functions not only in the nucleus/nucleolus but also at fungal specific structures such as septa and appressorial pores.

Lessons from the genome sequence of Neurospora crassa: tracing the path from genomic blueprint to multicellular organism.

We present an analysis of over 1,100 of the approximately 10,000 predicted proteins encoded by the genome sequence of the filamentous fungus Neurospora crassa. Seven major areas of Neurospora genomics and biology are covered. First, the basic features of the genome, including the automated assembly, gene calls, and global gene analyses are summarized. The second section covers components of the centromere and kinetochore complexes, chromatin assembly and modification, and transcription and translation initiation factors. The third area discusses genome defense mechanisms, including repeat induced point mutation, quelling and meiotic silencing, and DNA repair and recombination. In the fourth section, topics relevant to metabolism and transport include extracellular digestion; membrane transporters; aspects of carbon, sulfur, nitrogen, and lipid metabolism; the mitochondrion and energy metabolism; the proteasome; and protein glycosylation, secretion, and endocytosis. Environmental sensing is the focus of the fifth section with a treatment of two-component systems; GTP-binding proteins; mitogen-activated protein, p21-activated, and germinal center kinases; calcium signaling; protein phosphatases; photobiology; circadian rhythms; and heat shock and stress responses. The sixth area of analysis is growth and development; it encompasses cell wall synthesis, proteins important for hyphal polarity, cytoskeletal components, the cyclin/cyclin-dependent kinase machinery, macroconidiation, meiosis, and the sexual cycle. The seventh section covers topics relevant to animal and plant pathogenesis and human disease. The results demonstrate that a large proportion of Neurospora genes do not have homologues in the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. The group of unshared genes includes potential new targets for antifungals as well as loci implicated in human and plant physiology and disease.

Recurrent locus-specific mutation resulting from a cryptic ectopic insertion in Neurospora.

New mutations are found among approximately 20% of progeny when one or both parents carry eas allele UCLA191 (eas(UCLA), easily wettable, hydrophobin-deficient, linkage group II). The mutations inactivate the wild-type allele of cya-8 (cytochrome aa3 deficient, linkage group VII), resulting in thin, "transparent" mycelial growth. Other eas alleles fail to produce cya-8 mutant progeny. The recurrent cya-8 mutations are attributed to repeat-induced point mutation (RIP) resulting from a duplicated copy of cya-8+ that was inserted ectopically at eas when the UCLA191 mutation occurred. As expected for RIP, eas(UCLA)-induced cya-8 mutations occur during nuclear proliferation prior to karyogamy. When only one parent is eas(UCLA), the new mutations arise exclusively in eas(UCLA) nuclei. Mutation of cya-8 is suppressed when a long unlinked duplication is present. Stable cya-8 mutations are effectively eliminated in crosses homozygous for rid, a recessive suppressor of RIP. The eas(UCLA) allele is associated with a long paracentric inversion. A discontinuity is present in eas(UCLA) DNA. The eas promoter is methylated in cya-8 progeny of eas(UCLA), presumably by the spreading of methylation beyond the adjoining RIP-inactivated duplication. These findings support a model in which an ectopic insertion that created a mutation at the target site acts as a locus-specific mutator via RIP.

Population genomic analysis of the rice blast fungus reveals specific events associated with expansion of three main clades.

We examined the genomes of 100 isolates of Magnaporthe oryzae (Pyricularia oryzae), the causal agent of rice blast disease. We grouped current field populations of M. oryzae into three major globally distributed groups. A genetically diverse group, clade 1, which may represent a group of closely related lineages, contains isolates of both mating types. Two well-separated clades, clades 2 and 3, appear to have arisen as clonal lineages distinct from the genetically diverse clade. Examination of genes involved in mating pathways identified clade-specific diversification of several genes with orthologs involved in mating behavior in other fungi. All isolates within each clonal lineage are of the same mating type. Clade 2 is distinguished by a unique deletion allele of a gene encoding a small cysteine-rich protein that we determined to be a virulence factor. Clade 3 isolates have a small deletion within the MFA2 pheromone precursor gene, and this allele is shared with an unusual group of isolates we placed within clade 1 that contain AVR1-CO39 alleles. These markers could be used for rapid screening of isolates and suggest specific events in evolution that shaped these populations. Our findings are consistent with the view that M. oryzae populations in Asia generate diversity through recombination and may have served as the source of the clades 2 and 3 isolates that comprise a large fraction of the global population.

A mitogen-activated protein kinase pathway essential for mating and contributing to vegetative growth in Neurospora crassa.

MAP kinases homologous to Saccharomyces cerevisiae Fus3p/Kss1p have been identified in plant pathogenic fungi and are required for pathogenicity and sexual reproduction. To better understand the role of MAP kinase signaling in Neurospora crassa, and to identify downstream target genes of the pathway, we isolated, cloned, and disrupted the FUS3 homolog mak-2. Ste12p is a transcription factor target of Fus3p that activates genes of the mating pathway in yeast, and we also characterized the N. crassa STE12 homolog pp-1. The mak-2 and pp-1 mutants have reduced growth rate, produce short aerial hyphae, and fail to develop protoperithecia. In addition, ascospores carrying null mutations of either gene are inviable. Subtractive cloning was used to isolate genes having reduced expression in the mak-2 mutant. Expression of some of these genes is protoperithecia specific and three of them are part of a gene cluster potentially involved in the production of a polyketide secondary metabolite. Microarray analysis was used to extend the analysis of gene expression in mak-2 and pp-1 mutants. The role of the MAP kinase pathway in both sexual and asexual development as well as secondary metabolism is consistent with the dual regulation of the mating process and pathogencity observed in fungal pathogens.

Directional Selection from Host Plants Is a Major Force Driving Host Specificity in Magnaporthe Species.

One major threat to global food security that requires immediate attention, is the increasing incidence of host shift and host expansion in growing number of pathogenic fungi and emergence of new pathogens. The threat is more alarming because, yield quality and quantity improvement efforts are encouraging the cultivation of uniform plants with low genetic diversity that are increasingly susceptible to emerging pathogens. However, the influence of host genome differentiation on pathogen genome differentiation and its contribution to emergence and adaptability is still obscure. Here, we compared genome sequence of 6 isolates of Magnaporthe species obtained from three different host plants. We demonstrated the evolutionary relationship between Magnaporthe species and the influence of host differentiation on pathogens. Phylogenetic analysis showed that evolution of pathogen directly corresponds with host divergence, suggesting that host-pathogen interaction has led to co-evolution. Furthermore, we identified an asymmetric selection pressure on Magnaporthe species. Oryza sativa-infecting isolates showed higher directional selection from host and subsequently tends to lower the genetic diversity in its genome. We concluded that, frequent gene loss or gain, new transposon acquisition and sequence divergence are host adaptability mechanisms for Magnaporthe species, and this coevolution processes is greatly driven by directional selection from host plants.

The fluffy gene of Neurospora crassa is necessary and sufficient to induce conidiophore development.

The fl (fluffy) gene of Neurospora crassa encodes a binuclear zinc cluster protein that regulates the production of asexual spores called macroconidia. Two other genes, acon-2 and acon-3, play major roles in controlling development. fl is induced specifically in differentiating tissue during conidiation and acon-2 plays a role in this induction. We examined the function of fl by manipulating its level of expression in wild-type and developmental mutant strains. Increasing expression of fl from a heterologous promoter in a wild-type genetic background is sufficient to induce conidiophore development. Elevated expression of fl leads to induction of development of the acon-2 mutant in nitrogen-starved cultures, but does not bypass the conidiation defect of the acon-3 mutant. These findings indicate that fl acts downstream of acon-2 and upstream of acon-3 in regulating gene expression during development. The eas, con-6, and con-10 genes are induced at different times during development. Morphological changes induced by artificially elevated fl expression in the absence of environmental cues were correlated with increased expression of eas, but not con-6 or con-10. Thus, although inappropriate expression of fl in vegetative hyphae is sufficient to induce conidial morphogenesis, complete reconstitution of development leading to the formation of mature conidia may require environmental signals to regulate fl activity and/or appropriate induction of fl expression in the developing conidiophore.

The exocyst complex: delivery hub for morphogenesis and pathogenesis in filamentous fungi.

Regulated by several small GTPases, the octameric exocyst complex directs the docking and tethering of exocytic vesicles to the destined plasma membrane sites, providing the precise spatiotemporal control of exocytosis. Although the exocyst components are well conserved among various fungal species, the mechanisms for the regulation of its assembly and activity are diverse. Exocytosis is crucial for the generation of cell polarity as well as the delivery of effector proteins in filamentous fungi, and thus plays an important role for fungal morphogenesis and pathogenicity on plant hosts. This review focuses on current findings about the roles of the exocyst complex in the morphogenesis and pathogenesis of filamentous fungi.

Gene discovery and gene expression in the rice blast fungus, Magnaporthe grisea: analysis of expressed sequence tags.

Over 28,000 expressed sequence tags (ESTs) were produced from cDNA libraries representing a variety of growth conditions and cell types. Several Magnaporthe grisea strains were used to produce the libraries, including a nonpathogenic strain bearing a mutation in the PMK1 mitogen-activated protein kinase. Approximately 23,000 of the ESTs could be clustered into 3,050 contigs, leaving 5,127 singleton sequences. The estimate of 8,177 unique sequences indicates that over half of the genes of the fungus are represented in the ESTs. Analysis of EST frequency reveals growth and cell type-specific patterns of gene expression. This analysis establishes criteria for identification of fungal genes involved in pathogenesis. A large fraction of the genes represented by ESTs have no known function or described homologs. Manual annotation of the most abundant cDNAs with no known homologs allowed us to identify a family of metallothionein proteins present in M. grisea, Neurospora crassa, and Fusarium graminearum. In addition, multiply represented ESTs permitted the identification of alternatively spliced mRNA species. Alternative splicing was rare, and in most cases, the alternate mRNA forms were unspliced, although alternative 5' splice sites were also observed.