Evolution of virulence in fungal plant pathogens: exploiting fungal genomics to control plant disease.

The propensity of a fungal pathogen to evolve virulence depends on features of its biology (e.g. mode of reproduction) and of its genome (e.g. amount of repetitive DNA). Populations of Leptosphaeria maculans, a pathogen of Brassica napus (canola), can evolve and overcome disease resistance bred into canola within three years of commercial release of a cultivar. Avirulence effector genes are key fungal genes that are complementary to resistance genes. In L. maculans these genes are embedded within inactivated transposable elements in genomic regions where they are readily mutated or deleted. The risk of resistance breakdown in the field can be minimised by monitoring disease severity of canola cultivars and virulence of fungal populations using high throughput molecular assays and by sowing canola cultivars with different resistance genes in subsequent years. This strategy has been exploited to avert yield losses due to blackleg disease in Australia.

Transposable element-assisted evolution and adaptation to host plant within the Leptosphaeria maculans-Leptosphaeria biglobosa species complex of fungal pathogens.

BACKGROUND: Many plant-pathogenic fungi have a tendency towards genome size expansion, mostly driven by increasing content of transposable elements (TEs). Through comparative and evolutionary genomics, five members of the Leptosphaeria maculans-Leptosphaeria biglobosa species complex (class Dothideomycetes, order Pleosporales), having different host ranges and pathogenic abilities towards cruciferous plants, were studied to infer the role of TEs on genome shaping, speciation, and on the rise of better adapted pathogens. RESULTS: L. maculans 'brassicae', the most damaging species on oilseed rape, is the only member of the species complex to have a TE-invaded genome (32.5%) compared to the other members genomes (<4%). These TEs had an impact at the structural level by creating large TE-rich regions and are suspected to have been instrumental in chromosomal rearrangements possibly leading to speciation. TEs, associated with species-specific genes involved in disease process, also possibly had an incidence on evolution of pathogenicity by promoting translocations of effector genes to highly dynamic regions and thus tuning the regulation of effector gene expression in planta. CONCLUSIONS: Invasion of L. maculans 'brassicae' genome by TEs followed by bursts of TE activity allowed this species to evolve and to better adapt to its host, making this genome species a peculiarity within its own species complex as well as in the Pleosporales lineage.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea.

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

A gene cluster responsible for biosynthesis of phomenoic acid in the plant pathogenic fungus, Leptosphaeria maculans.

Phomenoic acid, a long chain aliphatic carboxylic acid is a major metabolite produced by Leptosphaeria maculans, a fungal pathogen of Brassica napus (canola). This fungus has 15 predicted polyketide synthases (PKS) and seven of them have the appropriate domains for the biosynthesis of phomenoic acid. The most highly expressed PKS gene after 7 days growth in 10% V8 juice, PKS2, was silenced and the resultant mutant produced very low levels of phomenoic acid, indicating that this PKS is involved in phomenoic acid biosynthesis. This gene is part of a co-regulated cluster of genes. Reduced expression of an adjacent gene encoding the transcriptional regulator C6TF, led to reduced expression of genes for PKS2, P450, a cytochrome P450 monoxygenase, YogA, an alcohol dehydrogenase/quinone reductase, RTA1, a lipid transport exporter superfamily member and MFS, a Major Facilitator Superfamily transporter, as well as a marked reduction in phomenoic acid production. Phomenoic acid is toxic towards another canola pathogen Leptosphaeria biglobosa 'canadensis', but not towards L. maculans and only moderately toxic towards the wheat pathogen Stagonospora nodorum. This molecule is detected in infected stems and stubble of B. napus, but biosynthesis of it does not appear to be essential for pathogenicity of L. maculans. Phomenoic acid may play a role in allowing L. maculans to outcompete other fungi in its environmental niche.

Multigene sequence data reveal morphologically cryptic phylogenetic species within the genus Laccaria in southern Australia.

Laccaria (Hydnangiaceae, Agaricales, Basidiomycota) is one of the more intensively studied ectomycorrhizal genera; however, species boundaries within Laccaria and the closely related Hydnangium and Podohydnangium in Australia have not yet been examined with molecular sequence data. Based on morphological characters, eight native species of Laccaria are currently recognized in Australia, as well as three Hydnangium species and the monotypic Podohydnangium australe. Sequences of the internal transcribed spacer region of nuclear rDNA (ITS), RNA polymerase beta subunit II (rpb2) and translation elongation factor 1 alpha (tef-1α) were generated from 77 collections of Laccaria, Hydnangium and Podohydnangium from Australia. Ten phylogenetic species and a further 11 potential species (represented by singletons) of Laccaria in Australia are delimited from sequence analyses. Most of the morphological species contained cryptic phylogenetic species, but these species were always nested entirely within a given morphological species, although not always as sister taxa. The rpb2 locus performed best as a species barcode with pairwise and patristic distance measures. The ITS sequence region returned the least resolved gene tree of the three regions examined and was the least useful as a barcode region. Based on the phylogenetic topology, there appears to have been multiple gains and/or losses of the ectomycorrhizal association of some species with the myrtle beech, Nothofagus cunninghamii as well as of sequestrate basidiocarps and two-spored basidia.

Evolution of linked avirulence effectors in Leptosphaeria maculans is affected by genomic environment and exposure to resistance genes in host plants.

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a 'gene for gene' manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans.

Spot Form of Net Blotch Resistance in a Diverse Set of Barley Lines in Australia and Canada.

The responses of 95 barley lines and cultivars to spot form of net blotch (SFNB) caused by Pyrenophora teres f. maculata were analyzed as seedlings and adults in Australia and Canada. Cluster analyses revealed complex reaction responses. Only 2 lines (Esperance Orge 289 and TR3189) were resistant to all isolates at the seedling stage, whereas 15 lines and cultivars (81-82/033, Arimont, BYDV-018, CBSS97M00855T-B2-M1-Y1-M2-Y-1M-0Y, CI9776, Keel, Sloop, Torrens, TR326, VB0111, Yarra, VB0229, WI-2477, WI2553, and Wisconsin Pedigree) were resistant toward the two Canadian isolates and mixture of Australian isolates at the adult stages. In Australian field experiments, the effectiveness of SFNB resistance in three barley cultivars (Barque, Cowabbie, and Schooner) and one breeding line (VB9104) with a different source of resistance was tested. Barque, which possessed a resistance gene that provided complete resistance to SFNB, was the most effective and showed no effect on grain yield or quality in the presence of inoculum. Generally, cultivars with seedling or adult resistance had less disease and better grain quality than the susceptible control, Dash, but they were not as effective as Barque. A preliminary differential set of 19 barley lines and cultivars for P. teres f. maculata is proposed.

Identification of cryptic products of the gliotoxin gene cluster using NMR-based comparative metabolomics and a model for gliotoxin biosynthesis.

Gliotoxin, a major product of the gli non-ribosomal peptide synthetase gene cluster, is strongly associated with virulence of the opportunistic human pathogen Aspergillus fumigatus. Despite identification of the gli cluster, the pathway of gliotoxin biosynthesis has remained elusive, in part because few potential intermediates have been identified. In addition, previous studies suggest that knowledge of gli-dependent metabolites is incomplete. Here we use differential analysis by 2D NMR spectroscopy (DANS) of metabolite extracts derived from gli knock-out and wild-type (WT) strains to obtain a detailed inventory of gli-dependent metabolites. DANS-based comparison of the WT metabolome with that of ΔgliZ, a knock-out strain devoid of the gene encoding the transcriptional regulator of the gli cluster, revealed nine novel gliZ-dependent metabolites including unexpected structural motifs. Their identification provides insight into gliotoxin biosynthesis and may benefit studies of the role of the gli cluster in A. fumigatus virulence. Our study demonstrates the utility of DANS for correlating gene expression and metabolite biosynthesis in microorganisms.

The cross-pathway control system regulates production of the secondary metabolite toxin, sirodesmin PL, in the ascomycete, Leptosphaeria maculans.

BACKGROUND: Sirodesmin PL is a secondary metabolite toxin made by the ascomycetous plant pathogen, Leptosphaeria maculans. The sirodesmin biosynthetic genes are clustered in the genome. The key genes are a non-ribosomal peptide synthetase, sirP, and a pathway-specific transcription factor, sirZ. Little is known about regulation of sirodesmin production. RESULTS: Genes involved in regulation of sirodesmin PL in L. maculans have been identified. Two hundred random insertional T-DNA mutants were screened with an antibacterial assay for ones producing low levels of sirodesmin PL. Three such mutants were isolated and each transcribed sirZ at very low levels. One of the affected genes had high sequence similarity to Aspergillus fumigatus cpcA, which regulates the cross-pathway control system in response to amino acid availability. This gene was silenced in L. maculans and the resultant mutant characterised. When amino acid starvation was artificially-induced by addition of 3-aminotriazole for 5 h, transcript levels of sirP and sirZ did not change in the wild type. In contrast, levels of sirP and sirZ transcripts increased in the silenced cpcA mutant. After prolonged amino acid starvation the silenced cpcA mutant produced much higher amounts of sirodesmin PL than the wild type. CONCLUSIONS: Production of sirodesmin PL in L. maculans is regulated by the cross pathway control gene, cpcA, either directly or indirectly via the pathway-specific transcription factor, sirZ.

Fungal pathogenesis: gene clusters unveiled as secrets within the Ustilago maydis code.

The genome sequence of a second plant pathogenic fungus is now available, revealing unique gene clusters encoding secretory proteins that are induced during infection and regulate pathogenesis. Gene clusters play important roles in pathogenic fungi, yet their evolution and maintenance remain a mystery.

Mutations to LmIFRD affect cell wall integrity, development and pathogenicity of the ascomycete Leptosphaeria maculans.

Maintaining cell wall integrity is essential for fungal growth and development. We describe two mutants with altered expression of a gene, LmIFRD, from the ascomycete Leptosphaeria maculans. Truncation of the LmIFRD transcript in a T-DNA insertional mutant led to slower germination, less sporulation and loss-of-pathogenicity towards Brassica napus, whereas silencing of the LmIFRD transcript led to increased germination, sporulation and earlier infection. The increased tolerance to cell wall lysing enzymes and cell wall-disrupting compounds of the T-DNA mutant contrasts with decreased tolerance of the silenced mutant and suggests altered cell wall integrity and accessibility to 1,3-linked glucan and chitin. Lectin binding experiments and monosaccharide analysis revealed altered polysaccharide content and structure within the cell wall of the LmIFRD mutants, notably increased 1,3-linked galactose and chitin within the cell wall of the T-DNA mutant. This is the first analysis of monosaccharide linkage composition of cell walls of spores and mycelia for any dothideomycete.

Cloning, purification and characterisation of brassinin glucosyltransferase, a phytoalexin-detoxifying enzyme from the plant pathogen Sclerotinia sclerotiorum.

The plant-pathogenic fungus Sclerotinia sclerotiorum can detoxify cruciferous phytoalexins such as brassinin via glucosylation. Here we describe a multifaceted approach including genome mining, transcriptional induction, phytoalexin quantification, protein expression and enzyme purification that led to identification of a S. sclerotiorum glucosyltransferase that detoxifies brassinin. Transcription of this gene, denoted as brassinin glucosyltransferase 1 (SsBGT1), was induced significantly in response to the cruciferous phytoalexins camalexin, cyclobrassinin, brassilexin, brassinin and 3-phenylindole, a camalexin analogue. This gene was also up-regulated during infection of Brassica napus leaves. Levels of brassinin decreased significantly between 48 and 72h post-inoculation, with a concomitant increase in levels of 1-beta-d-glucopyranosylbrassinin, the product of the reaction catalysed by SsBGT1. These findings strongly implicate the involvement of this gene during infection of B. napus. This gene was cloned and expressed in Saccharomyces cerevisiae. The purified recombinant enzyme was able to glucosylate brassinin and two other phytoalexins, albeit much less effectively. This is the first report of a fungal gene involved in detoxification of plant defence molecules via glucosylation.

Hunting down fungal secretomes using liquid-phase IEF prior to high resolution 2-DE.

The secreted proteins (secretome) of fungi play a key role in interactions of pathogenic and symbiotic fungi with plants. Using the plant pathogenic fungus Leptosphaeria maculans and symbiont Laccaria bicolor grown in culture, we have established a proteomic protocol for extraction, concentration and resolution of the fungal secretome. As no proteomic data were available on mycelium tissues from both L. maculans and L. bicolor, mycelial proteins were studied; they also helped verifying the purity of secretome samples. The quality of protein extracts was initially assessed by both 1-DE and 2-DE using first a broad pH range for IEF, and then narrower acidic and basic pH ranges, prior to 2-DE. Compared with the previously published protocols for which only dozens of 2-D spots were recovered from fungal secretome samples, up to approximately 2000 2-D spots were resolved by our method. MS identification of proteins along several pH gradients confirmed this high resolution, as well as the presence of major secretome markers such as endopolygalacturonases, beta-glucanosyltransferases, pectate lyases and endoglucanases. Shotgun proteomic experiments evidenced the enrichment of secreted protein within the liquid medium. This is the first description of the proteome of L. maculans and L. bicolor, and the first application of liquid-phase IEF to any fungal extracts.

Secondary metabolite toxins and nutrition of plant pathogenic fungi.

Fungal pathogens derive nutrition from the plants they invade. Some fungi can subvert plant defence responses such as programmed cell death to provide nutrition for their growth and colonisation. Secondary metabolite toxins produced by fungi often play a role in triggering these responses. Knowledge of the biosynthesis of these toxins, and the availability of fungal genome sequences and gene disruption techniques, allows the development of tools for experiments aimed at discovering the role of such toxins in triggering plant cell death and plant disease.

Fungi have three tetraspanin families with distinct functions.

BACKGROUND: Tetraspanins are small membrane proteins that belong to a superfamily encompassing 33 members in human and mouse. These proteins act as organizers of membrane-signalling complexes. So far only two tetraspanin families have been identified in fungi. These are Pls1, which is required for pathogenicity of the plant pathogenic ascomycetes, Magnaporthe grisea, Botrytis cinerea and Colletotrichum lindemuthianum, and Tsp2, whose function is unknown. In this report, we describe a third family of tetraspanins (Tsp3) and a new family of tetraspanin-like proteins (Tpl1) in fungi. We also describe expression of some of these genes in M. grisea and a basidiomycete, Laccaria bicolor, and also their functional analysis in M. grisea. RESULTS: The exhaustive search for tetraspanins in fungal genomes reveals that higher fungi (basidiomycetes and ascomycetes) contain three families of tetraspanins (Pls1, Tsp2 and Tsp3) with different distribution amongst phyla. Pls1 is found in ascomycetes and basidiomycetes, whereas Tsp2 is restricted to basidiomycetes and Tsp3 to ascomycetes. A unique copy of each of PLS1 and TSP3 was found in ascomycetes in contrast to TSP2, which has several paralogs in the basidiomycetes, Coprinus cinereus and Laccaria bicolor. A tetraspanin-like family (Tpl1) was also identified in ascomycetes. Transcriptional analyses in various tissues of L. bicolor and M. grisea showed that PLS1 and TSP2 are expressed in all tissues in L. bicolor and that TSP3 and TPL1 are overexpressed in the sexual fruiting bodies (perithecia) and mycelia of M. grisea, suggesting that these genes are not pseudogenes. Phenotypic analysis of gene replacementmutants Deltatsp3 and Deltatpl1 of M. grisea revealed a reduction of the pathogenicity only on rice, in contrast to Deltapls1 mutants, which are completely non-pathogenic on barley and rice. CONCLUSION: A new tetraspanin family (Tsp3) and a tetraspanin-like protein family (Tpl1) have been identified in fungi. Functional analysis by gene replacement showed that these proteins, as well as Pls1, are involved in the infection process of the plant pathogenic fungus M. grisea. The next challenge will be to decipher the role(s) of tetraspanins in a range of symbiotic, saprophytic and human pathogenic fungi.

A Zn(II)2Cys6 DNA binding protein regulates the sirodesmin PL biosynthetic gene cluster in Leptosphaeria maculans.

A gene, sirZ, encoding a Zn(II)(2)Cys(6) DNA binding protein is present in a cluster of genes responsible for the biosynthesis of the epipolythiodioxopiperazine (ETP) toxin, sirodesmin PL in the ascomycete plant pathogen, Leptosphaeria maculans. RNA-mediated silencing of sirZ gives rise to transformants that produce only residual amounts of sirodesmin PL and display a decrease in the transcription of several sirodesmin PL biosynthetic genes. This indicates that SirZ is a major regulator of this gene cluster. Proteins similar to SirZ are encoded in the gliotoxin biosynthetic gene cluster of Aspergillus fumigatus (gliZ) and in an ETP-like cluster in Penicillium lilacinoechinulatum (PlgliZ). Despite its high level of sequence similarity to gliZ, PlgliZ is unable to complement the gliotoxin-deficiency of a mutant of gliZ in A. fumigatus. Putative binding sites for these regulatory proteins in the promoters of genes in these clusters were predicted using bioinformatic analysis. These sites are similar to those commonly bound by other proteins with Zn(II)(2)Cys(6) DNA binding domains.

Origin and distribution of epipolythiodioxopiperazine (ETP) gene clusters in filamentous ascomycetes.

BACKGROUND: Genes responsible for biosynthesis of fungal secondary metabolites are usually tightly clustered in the genome and co-regulated with metabolite production. Epipolythiodioxopiperazines (ETPs) are a class of secondary metabolite toxins produced by disparate ascomycete fungi and implicated in several animal and plant diseases. Gene clusters responsible for their production have previously been defined in only two fungi. Fungal genome sequence data have been surveyed for the presence of putative ETP clusters and cluster data have been generated from several fungal taxa where genome sequences are not available. Phylogenetic analysis of cluster genes has been used to investigate the assembly and heredity of these gene clusters. RESULTS: Putative ETP gene clusters are present in 14 ascomycete taxa, but absent in numerous other ascomycetes examined. These clusters are discontinuously distributed in ascomycete lineages. Gene content is not absolutely fixed, however, common genes are identified and phylogenies of six of these are separately inferred. In each phylogeny almost all cluster genes form monophyletic clades with non-cluster fungal paralogues being the nearest outgroups. This relatedness of cluster genes suggests that a progenitor ETP gene cluster assembled within an ancestral taxon. Within each of the cluster clades, the cluster genes group together in consistent subclades, however, these relationships do not always reflect the phylogeny of ascomycetes. Micro-synteny of several of the genes within the clusters provides further support for these subclades. CONCLUSION: ETP gene clusters appear to have a single origin and have been inherited relatively intact rather than assembling independently in the different ascomycete lineages. This progenitor cluster has given rise to a small number of distinct phylogenetic classes of clusters that are represented in a discontinuous pattern throughout ascomycetes. The disjunct heredity of these clusters is discussed with consideration to multiple instances of independent cluster loss and lateral transfer of gene clusters between lineages.

Microsatellite and Minisatellite Analysis of Leptosphaeria maculans in Australia Reveals Regional Genetic Differentiation.

ABSTRACT The population genetic structure of the fungal pathogen Leptosphaeria maculans was determined in Australia using six microsatellite and two minisatellite markers. Ascospores were sampled from Brassica napus stubble in disease nurseries and commercial fields in different sites over 2 years. The 13 subpopulations of L. maculans exhibited high gene (H = 0.393 to 0.563) and genotypic diversity, with 357 haplotypes identified among 513 isolates. Although the majority of genetic variation was distributed within subpopulations (85%), 10% occurred between the regions of eastern and Western Australia, and 5% within regions. F(ST) analysis of subpopulation pairs also showed the east-west genetic differentiation, whereas factorial correspondence analysis separated Western Australian subpopulations from eastern ones. Bayesian model-based population structure analyses of multilocus haplotypes inferred three distinct populations, one in Western Australia and an admixture of two in eastern Australia. These two regions are separated by 1,200 km of arid desert that may act as a natural barrier to gene flow, resulting in differentiation by random genetic drift. The genetic differentiation of L. maculans isolates between eastern and Western Australia means that these regions can be treated as different management units, and reinforces the need for widespread disease nurseries in each region to screen breeding lines against a range of genetic and pathogenic populations of L. maculans.

Overexpression of a 3-ketoacyl-CoA thiolase in Leptosphaeria maculans causes reduced pathogenicity on Brassica napus.

Agrobacterium tumefaciens-mediated random mutagenesis was used to generate insertional mutants of the fungus Leptosphaeria maculans. Of 91 transformants screened, only one (A3) produced lesions of reduced size on cotyledons of canola (Brassica napus). Genes flanking the T-DNA insertion had the best matches to an alcohol dehydrogenase class 4 (ADH4)-like gene (Adh4L) and a 3-ketoacyl-CoA thiolase gene (Thiol) and were expressed in mutant A3 in vitro and in planta at significantly higher levels than in the wild type. This is the first report of a T-DNA insertion in fungi causing increased gene expression. Transformants of the wild-type isolate expressing both Adh4L and Thiol under the control of a heterologous promoter had similar pathogenicity to mutant A3. Ectopic expression of only thiolase resulted in loss of pathogenicity, suggesting that thiolase overexpression was primarily responsible for the reduced pathogenicity of the A3 isolate. The thiolase gene encoded a functional protein, as shown by assays in which a nontoxic substrate (2, 4 dichlorophenoxybutyric acid) was converted to a toxic product. The use of a translational fusion with a reporter gene showed thiolase expressed in organelles that are most likely peroxisomes.