Aflatoxins are toxic and carcinogenic polyketides produced by several Aspergillus species that are known to contaminate agricultural commodities, posing a serious threat to animal and human health. Aflatoxin (AF) biosynthesis is almost fully characterized and involves the coordinated expression of approximately 25 genes clustered in a 70-kb DNA region. Aspergillus parasiticus is an economically important and common agent of AF contamination. Naturally occurring nonaflatoxigenic strains of A. parasiticus are rarely found and generally produce O-methylsterigmatocystin (OMST), the immediate precursor of AF. To elucidate the evolutionary forces acting to retain AF and OMST pathway extrolites (chemotypes), we sequenced 21 intergenic regions spanning the entire cluster in 24 A. parasiticus isolates chosen to represent the genetic diversity within a single Georgia field population. Linkage disequilibrium analyses revealed five distinct recombination blocks in the A. parasiticus cluster. Phylogenetic network analyses showed a history of recombination between chemotype-specific haplotypes, as well as evidence of contemporary recombination. We performed coalescent simulations of variation in recombination blocks and found an approximately twofold deeper coalescence for cluster genealogies compared to noncluster genealogies, our internal standard of neutral evolution. Significantly deeper cluster genealogies are indicative of balancing selection in the AF cluster of A. parasiticus and are further corroborated by the existence of trans-species polymorphisms and common haplotypes in the cluster for several closely related species. Estimates of Ka/Ks for representative cluster genes provide evidence of selection for OMST and AF chemotypes, and indicate a possible role of chemotypes in ecological adaptation and speciation.

BackgroundThe biosynthesis of aflatoxin (AF) involves over 20 enzymatic reactions in a complex polyketide pathway that converts acetate and malonate to the intermediates sterigmatocystin (ST) and O-methylsterigmatocystin (OMST), the respective penultimate and ultimate precursors of AF. Although these precursors are chemically and structurally very similar, their accumulation differs at the species level for Aspergilli. Notable examples are A. nidulans that synthesizes only ST, A. flavus that makes predominantly AF, and A. parasiticus that generally produces either AF or OMST. Whether these differences are important in the evolutionary/ecological processes of species adaptation and diversification is unknown. Equally unknown are the specific genomic mechanisms responsible for ordering and clustering of genes in the AF pathway of Aspergillus.ResultsTo elucidate the mechanisms that have driven formation of these clusters, we performed systematic searches of aflatoxin cluster homologs across five Aspergillus genomes. We found a high level of gene duplication and identified seven modules consisting of highly correlated gene pairs (aflA/aflB, aflR/aflS, aflX/aflY, aflF/aflE, aflT/aflQ, aflC/aflW, and aflG/aflL). With the exception of A. nomius, contrasts of mean Ka/Ks values across all cluster genes showed significant differences in selective pressure between section Flavi and non-section Flavi species. A. nomius mean Ka/Ks values were more similar to partial clusters in A. fumigatus and A. terreus. Overall, mean Ka/Ks values were significantly higher for section Flavi than for non-section Flavi species.ConclusionOur results implicate several genomic mechanisms in the evolution of ST, OMST and AF cluster genes. Gene modules may arise from duplications of a single gene, whereby the function of the pre-duplication gene is retained in the copy (aflF/aflE) or the copies may partition the ancestral function (aflA/aflB). In some gene modules, the duplicated copy may simply augment/supplement a specific pathway function (aflR/aflS and aflX/aflY) or the duplicated copy may evolve a completely new function (aflT/aflQ and aflC/aflW). Gene modules that are contiguous in one species and noncontiguous in others point to possible rearrangements of cluster genes in the evolution of these species. Significantly higher mean Ka/Ks values in section Flavi compared to non-section Flavi species indicate increased positive selection acting in the evolution of genes in OMST and AF gene clusters.

Transcripts for phenylalanine ammonia-lyase, chalcone synthase, chalcone isomerase, and chitinase accumulated in common bean after infiltration with the Pseudomonas syringae pv tabaci Hrp- mutant Pt11528::Hrp1, even though a hypersensitive reaction did not occur. The temporal pattern of this transcript accumulation was similar to that seen after infiltration with wild-type P. s. tabaci Pt11528, which resulted in a hypersensitive reaction. Escherichia coli DH5[alpha], P. fluorescens Pf101, heat-killed Pt11528 cells, and Pt11528 cells treated with protein synthesis inhibitors also induced accumulation of defense transcripts but not a hypersensitive reaction. In contrast, these transcripts were not detected in plants infiltrated with water or P.s. pv phaseolicola NPS3121, a compatible pathogen that causes halo blight. Phytoalexins were produced in bean after infiltration with Pt11528, Pt11528::Hrp1, Pt11528 cells treated with neomycin, or Pf101, but not in plants infiltrated with NPS3121 or water. These results suggest that there are unique biochemical events associated with the expression of a hypersensitive reaction which are distinct from other plant defense responses such as phytoalexin biosynthesis. In addition, our results support the hypothesis that there is a general, nonspecific mechanism for the induction of defense transcripts and phytoalexins by pathogenic and saprophytic bacteria that is distinct from the more specific mechanism associated with the induction of the hypersensitive reaction.

We have developed a model system to examine suppression of defense responses in bean by the compatible bacterium Pseudomonas syringae pv phaseolicola. Previously, we have shown that there is a general mechanism for the induction of the bean defense genes phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), and chitinase (CHT) by incompatible, compatible, and nonpathogenic bacteria. Here, we show that bean plants infiltrated with isolates of P. s. phaseolicola failed to produce transcripts for PAL, CHS, or CHI up to 120 hr after infiltration and CHT transcript accumulation was significantly delayed when compared to the incompatible P. syringae strains. Infiltration of bean plants with 108 cells per mL of P. s. phaseolicola NPS3121 8 hr prior to infiltration with an equal concentration of incompatible P. s. pv tabaci Pt11528 significantly reduced the typical profile of defense transcript accumulation when compared to plants infiltrated with Pt11528 alone. A corresponding suppression of phytoalexin accumulation was also observed. NPS3121 also suppressed PAL, CHS, CHI, and CHT transcript accumulation and phytoalexin production induced by Escherichia coli DH5[alpha] or the elicitor glutathione. Heat-killed NPS3121 cells or cells treated with protein synthesis inhibitors lost the suppressor activity. Taken together, these experiments suggest that NPS3121 has an active mechanism to suppress the accumulation of defense transcripts and phytoalexin biosynthesis in bean.

Transcripts for phenylalanine ammonia-lyase, chalcone synthase, chalcone isomerase, and chitinase accumulated in common bean after infiltration with the Pseudomonas syringae pv tabaci Hrp- mutant Pt11528::Hrp1, even though a hypersensitive reaction did not occur. The temporal pattern of this transcript accumulation was similar to that seen after infiltration with wild-type P. s. tabaci Pt11528, which resulted in a hypersensitive reaction. Escherichia coli DH5[alpha], P. fluorescens Pf101, heat-killed Pt11528 cells, and Pt11528 cells treated with protein synthesis inhibitors also induced accumulation of defense transcripts but not a hypersensitive reaction. In contrast, these transcripts were not detected in plants infiltrated with water or P.s. pv phaseolicola NPS3121, a compatible pathogen that causes halo blight. Phytoalexins were produced in bean after infiltration with Pt11528, Pt11528::Hrp1, Pt11528 cells treated with neomycin, or Pf101, but not in plants infiltrated with NPS3121 or water. These results suggest that there are unique biochemical events associated with the expression of a hypersensitive reaction which are distinct from other plant defense responses such as phytoalexin biosynthesis. In addition, our results support the hypothesis that there is a general, nonspecific mechanism for the induction of defense transcripts and phytoalexins by pathogenic and saprophytic bacteria that is distinct from the more specific mechanism associated with the induction of the hypersensitive reaction.