Analysis of loss of pathogenicity mutants reveals that repeat-induced point mutations can occur in the Dothideomycete Leptosphaeria maculans

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

Acidovorax pan-genome reveals specific functional traits for plant beneficial and pathogenic plant-associations

Beta-proteobacteria belonging to the genus Acidovorax have been described from various environments. Many strains can interact with a range of hosts, including humans and plants, forming neutral, beneficial or detrimental associations. In the frame of this study, we investigated the genomic properties of 52 bacterial strains of the genus Acidovorax, isolated from healthy roots of Lotus japonicus, with the intent of identifying traits important for effective plant-growth promotion. Based on single-strain inoculation bioassays with L. japonicus, performed in a gnotobiotic system, we distinguished seven robust plant-growth promoting strains from strains with no significant effects on plant-growth. We showed that the genomes of the two groups differed prominently in protein families linked to sensing and transport of organic acids, production of phytohormones, as well as resistance and production of compounds with antimicrobial properties. In a second step, we compared the genomes of the tested isolates with those of plant pathogens and free-living strains of the genus Acidovorax sourced from public repositories. Our pan-genomics comparison revealed features correlated with commensal and pathogenic lifestyle. We showed that commensals and pathogens differ mostly in their ability to use plant-derived lipids and in the type of secretion-systems being present. Most free-living Acidovorax strains did not harbour any secretion-systems. Overall, our data indicate that Acidovorax strains undergo extensive adaptations to their particular lifestyle by horizontal uptake of novel genetic information and loss of unnecessary genes.

Accessory Chromosome-Acquired Secondary Metabolism in Plant Pathogenic Fungi: The Evolution of Biotrophs Into Host-Specific Pathogens

Accessory chromosomes are strain- or pathotype-specific chromosomes that exist in addition to the core chromosomes of a species and are generally not considered essential to the survival of the organism. Among pathogenic fungal species, accessory chromosomes harbor pathogenicity or virulence factor genes, several of which are known to encode for secondary metabolites that are involved in plant tissue invasion. Accessory chromosomes are of particular interest due to their capacity for horizontal transfer between strains and their dynamic "crosstalk" with core chromosomes. This review focuses exclusively on secondary metabolism (including mycotoxin biosynthesis) associated with accessory chromosomes in filamentous fungi and the role accessory chromosomes play in the evolution of secondary metabolite gene clusters. Untargeted metabolomics profiling in conjunction with genome sequencing provides an effective means of linking secondary metabolite products with their respective biosynthetic gene clusters that reside on accessory chromosomes. While the majority of literature describing accessory chromosome-associated toxin biosynthesis comes from studies of Alternaria pathotypes, the recent discovery of accessory chromosome-associated biosynthetic genes in Fusarium species offer fresh insights into the evolution of biosynthetic enzymes such as non-ribosomal peptide synthetases (NRPSs), polyketide synthases (PKSs) and regulatory mechanisms governing their expression.

Large-scale transcriptomics to dissect 2 years of the life of a fungal phytopathogen interacting with its host plant

BACKGROUND

The fungus Leptosphaeria maculans has an exceptionally long and complex relationship with its host plant, Brassica napus, during which it switches between different lifestyles, including asymptomatic, biotrophic, necrotrophic, and saprotrophic stages. The fungus is also exemplary of "two-speed" genome organisms in the genome of which gene-rich and repeat-rich regions alternate. Except for a few stages of plant infection under controlled conditions, nothing is known about the genes mobilized by the fungus throughout its life cycle, which may last several years in the field.

RESULTS

We performed RNA-seq on samples corresponding to all stages of the interaction of L. maculans with its host plant, either alive or dead (stem residues after harvest) in controlled conditions or in field experiments under natural inoculum pressure, over periods of time ranging from a few days to months or years. A total of 102 biological samples corresponding to 37 sets of conditions were analyzed. We show here that about 9% of the genes of this fungus are highly expressed during its interactions with its host plant. These genes are distributed into eight well-defined expression clusters, corresponding to specific infection lifestyles or to tissue-specific genes. All expression clusters are enriched in effector genes, and one cluster is specific to the saprophytic lifestyle on plant residues. One cluster, including genes known to be involved in the first phase of asymptomatic fungal growth in leaves, is re-used at each asymptomatic growth stage, regardless of the type of organ infected. The expression of the genes of this cluster is repeatedly turned on and off during infection. Whatever their expression profile, the genes of these clusters are enriched in heterochromatin regions associated with H3K9me3 or H3K27me3 repressive marks. These findings provide support for the hypothesis that part of the fungal genes involved in niche adaptation is located in heterochromatic regions of the genome, conferring an extreme plasticity of expression.

CONCLUSION

This work opens up new avenues for plant disease control, by identifying stage-specific effectors that could be used as targets for the identification of novel durable disease resistance genes, or for the in-depth analysis of chromatin remodeling during plant infection, which could be manipulated to interfere with the global expression of effector genes at crucial stages of plant infection.

Limitations of transcriptome-based prediction of pathogenicity genes in the plant pathogen Leptosphaeria maculans

Identification of pathogenicity determinants in Leptosphaeria maculans, a major cause of disease of oilseed crops, has been a focus of research for many years. A wealth of gene expression information from RNA sequencing promises to illuminate the mechanisms by which the fungus is able to cause blackleg disease. However, to date, no studies have tested the hypothesis that high gene transcript levels during infection correlate with importance to disease progression. In this study, we use CRISPR-Cas9 to disrupt 11 genes that are highly expressed during the early stages of disease and show that none of these genes are crucial for fungal pathogenicity on Brassica napus. This finding suggests that in order to understand the pathogenicity of this fungus more sophisticated techniques than simple expression analysis will need to be employed.

Sterol Demethylation Inhibitor Fungicide Resistance in Leptosphaeria maculans is Caused by Modifications in the Regulatory Region of ERG11

Blackleg is a worldwide disease of canola (Brassica napus), caused by a complex of fungal species in the genus Leptosphaeria, that impacts canola production and seed quality. Demethylation inhibitor (DMI) fungicides that target sterol 14α-demethylase are an integral part of disease control. Here, we report six DMI-resistant isolates of Leptosphaeria maculans and two different types of genetic modification related to the resistance. Analysis of the regulatory region of the DMI target gene ERG11 (also known as CYP51) revealed a 275-bp insertion in two of the isolates and three long terminal repeat retrotransposons (5,263, 5,267, and 5,248 bp) inserted in the promoter region of three resistant isolates. Genetic approaches confirmed that these elements are responsible for DMI resistance in L. maculans and crosses show segregation consistent with a single locus. Reverse-transcription quantitative PCR assays demonstrated that the 275-bp insertion increases ERG11 gene expression, conferring DMI fungicide resistance both in vitro and in planta. Moreover, transformation of a susceptible isolate of L. maculans with ERG11 driven by a promoter containing the 275-bp insertion increased resistance to tebuconazole. A minimal shift of the values of concentration whereby 50% of the mycelial growth is inhibited in vitro was observed in resistant isolates containing long terminal repeat retrotransposons; nevertheless, these isolates were able to develop significant lesions on cotyledons from fungicide-treated seedlings. This is the first report of genetic modifications in L. maculans relating to DMI fungicide resistance.

Enhancing peptaibols production in the biocontrol fungus Trichoderma longibrachiatum SMF2 by elimination of a putative glucose sensor

Trichoderma spp. are main producers of peptide antibiotics known as peptaibols. While peptaibols have been shown to possess a range of biological activities, molecular understanding of the regulation of their production is largely unclear, which hampers the production improvement through genetic engineering. Here, we demonstrated that the orthologue of glucose sensors in the outstanding biocontrol fungus Trichoderma longibrachiatum SMF2, TlSTP1, participates in the regulation of peptaibols production. Deletion of Tlstp1 markedly impaired hyphal growth and conidiation, but significantly increased peptaibols yield by 5-fold for Trichokonins A and 2.6-fold for Trichokonins B. Quantitative real-time polymerase chain reaction analyses showed that the increased peptaibols production occurs at the transcriptional levels of the two nonribosomal peptide synthetase encoding genes, tlx1 and tlx2. Transcriptome analyses of the wild type and the Tlstp1 mutant strains indicated that TlSTP1 exerts a regulatory effect on a set of genes that are involved in a number of metabolic and cellular processes, including synthesis of several other secondary metabolites. These results suggest an important role of TlSTP1 in the regulation of vegetative growth and peptaibols production in T. longibrachiatum SMF2 and provide insights into construction of peptaibol-hyperproducing strains through genetic engineering.

Analysis of Repeat Induced Point (RIP) Mutations in Leptosphaeria maculans Indicates Variability in the RIP Process Between Fungal Species

Gene duplication contributes to evolutionary potential, yet many duplications in a genome arise from the activity of "selfish" genetic elements such as transposable elements. Fungi have a number of mechanisms by which they limit the expansion of transposons, including Repeat Induced Point mutation (RIP). RIP has been best characterized in the Sordariomycete Neurospora crassa, wherein duplicated DNA regions are recognized after cell fusion, but before nuclear fusion during the sexual cycle, and then mutated. While "signatures" of RIP appear in the genome sequences of many fungi, the species most distant from N. crassa in which the process has been experimentally demonstrated to occur is the Dothideomycete Leptosphaeria maculans In the current study, we show that similar to N. crassa, nonlinked duplications can trigger RIP; however, the frequency of the generated RIP mutations is extremely low in L maculans (< 0.1%) and requires a large duplication to initiate RIP, and that multiple premeiotic mitoses are involved in the RIP process. However, a single sexual cycle leads to the generation of progeny with unique haplotypes, despite progeny pairs being generated from mitosis. We hypothesize that these different haplotypes may be the result of the deamination process occurring post karyogamy, leading to unique mutations within each of the progeny pairs. These findings indicate that the RIP process, while common to many fungi, differs between fungi and that this impacts on the fate of duplicated DNA.

Genomic and Genetic Insights Into a Cosmopolitan Fungus, Paecilomyces variotii (Eurotiales)

Species in the genus Paecilomyces, a member of the fungal order Eurotiales, are ubiquitous in nature and impact a variety of human endeavors. Here, the biology of one common species, Paecilomyces variotii, was explored using genomics and functional genetics. Sequencing the genome of two isolates revealed key genome and gene features in this species. A striking feature of the genome was the two-part nature, featuring large stretches of DNA with normal GC content separated by AT-rich regions, a hallmark of many plant-pathogenic fungal genomes. These AT-rich regions appeared to have been mutated by repeat-induced point (RIP) mutations. We developed methods for genetic transformation of P. variotii, including forward and reverse genetics as well as crossing techniques. Using transformation and crossing, RIP activity was identified, demonstrating for the first time that RIP is an active process within the order Eurotiales. A consequence of RIP is likely reflected by a reduction in numbers of genes within gene families, such as in cell wall degradation, and reflected by growth limitations on P. variotii on diverse carbon sources. Furthermore, using these transformation tools we characterized a conserved protein containing a domain of unknown function (DUF1212) and discovered it is involved in pigmentation.

New insights into Phakopsora pachyrhizi infection based on transcriptome analysis in planta

Asian soybean rust (ASR) is one of the most destructive diseases affecting soybeans. The causative agent of ASR, the fungus Phakopsora pachyrhizi, presents characteristics that make it difficult to study in vitro, limiting our knowledge of plant-pathogen dynamics. Therefore, this work used leaf lesion laser microdissection associated with deep sequencing to determine the pathogen transcriptome during compatible and incompatible interactions with soybean. The 36,350 generated unisequences provided an overview of the main genes and biological pathways that were active in the fungus during the infection cycle. We also identified the most expressed transcripts, including sequences similar to other fungal virulence and signaling proteins. Enriched P. pachyrhizi transcripts in the resistant (PI561356) soybean genotype were related to extracellular matrix organization and metabolic signaling pathways and, among infection structures, in amino acid metabolism and intracellular transport. Unisequences were further grouped into gene families along predicted sequences from 15 other fungi and oomycetes, including rust fungi, allowing the identification of conserved multigenic families, as well as being specific to P. pachyrhizi. The results revealed important biological processes observed in P. pachyrhizi, contributing with information related to fungal biology and, consequently, a better understanding of ASR.

Effectiveness of Rlm7 resistance against Leptosphaeria maculans (phoma stem canker) in UK winter oilseed rape cultivars

The Rlm7 gene in Brassica napus is an important source of resistance for control of phoma stem canker on oilseed rape caused by the fungus Leptosphaeria maculans. This study shows the first report of L. maculans isolates virulent against Rlm7 in the UK. Leptosphaeria maculans isolates virulent against Rlm7 represented 3% of the pathogen population when cultivars with the Rlm7 gene represented 5% of the UK oilseed rape area in 2012/13. However, the Rlm7 gene has been widely used since then, representing >15% of the UK oilseed rape area in 2015/16. Winter oilseed rape field experiments included cultivars with the Rlm7 gene, with the Rlm4 gene or without Rlm genes and took place at five sites in the UK over four cropping seasons. An increase in phoma leaf spotting severity on Rlm7 cultivars in successive seasons was observed. Major resistance genes played a role in preventing severe phoma leaf spotting at the beginning of the cropping season and, in addition, quantitative resistance (QR) in the cultivars examined made an important contribution to control of phoma stem canker development at the end of the cropping season. Deployment of the Rlm7 resistance gene against L. maculans in cultivars with QR in combination with sustainable disease management practices will prolong the use of this gene for effective control of phoma stem canker epidemics.

Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

BACKGROUND

Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species.

RESULTS

To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates.

CONCLUSION

Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.

Reference Assembly and Annotation of the Pyrenophora teres f. teres Isolate 0-1

Pyrenophora teres f. teres, the causal agent of net form net blotch (NFNB) of barley, is a destructive pathogen in barley-growing regions throughout the world. Typical yield losses due to NFNB range from 10 to 40%; however, complete loss has been observed on highly susceptible barley lines where environmental conditions favor the pathogen. Currently, genomic resources for this economically important pathogen are limited to a fragmented draft genome assembly and annotation, with limited RNA support of the P. teres f. teres isolate 0-1. This research presents an updated 0-1 reference assembly facilitated by long-read sequencing and scaffolding with the assistance of genetic linkage maps. Additionally, genome annotation was mediated by RNAseq analysis using three infection time points and a pure culture sample, resulting in 11,541 high-confidence gene models. The 0-1 genome assembly and annotation presented here now contains the majority of the repetitive content of the genome. Analysis of the 0-1 genome revealed classic characteristics of a “two-speed” genome, being compartmentalized into GC-equilibrated and AT-rich compartments. The assembly of repetitive AT-rich regions will be important for future investigation of genes known as effectors, which often reside in close proximity to repetitive regions. These effectors are responsible for manipulation of the host defense during infection. This updated P. teres f. teres isolate 0-1 reference genome assembly and annotation provides a robust resource for the examination of the barley–P. teres f. teres host–pathogen coevolution.

A silver bullet in a golden age of functional genomics: the impact of Agrobacterium-mediated transformation of fungi

The implementation of Agrobacterium tumefaciens as a transformation tool revolutionized approaches to discover and understand gene functions in a large number of fungal species. A. tumefaciens mediated transformation (AtMT) is one of the most transformative technologies for research on fungi developed in the last 20 years, a development arguably only surpassed by the impact of genomics. AtMT has been widely applied in forward genetics, whereby generation of strain libraries using random T-DNA insertional mutagenesis, combined with phenotypic screening, has enabled the genetic basis of many processes to be elucidated. Alternatively, AtMT has been fundamental for reverse genetics, where mutant isolates are generated with targeted gene deletions or disruptions, enabling gene functional roles to be determined. When combined with concomitant advances in genomics, both forward and reverse approaches using AtMT have enabled complex fungal phenotypes to be dissected at the molecular and genetic level. Additionally, in several cases AtMT has paved the way for the development of new species to act as models for specific areas of fungal biology, particularly in plant pathogenic ascomycetes and in a number of basidiomycete species. Despite its impact, the implementation of AtMT has been uneven in the fungi. This review provides insight into the dynamics of expansion of new research tools into a large research community and across multiple organisms. As such, AtMT in the fungi, beyond the demonstrated and continuing power for gene discovery and as a facile transformation tool, provides a model to understand how other technologies that are just being pioneered, e.g. CRISPR/Cas, may play roles in fungi and other eukaryotic species.

Sit4-Associated Protein is Required for Pathogenicity of Leptosphaeria maculans on Brassica napus

An insertional mutant with reduced pathogenicity on Brassica napus was identified in the plant pathogenic fungus Leptosphaeria maculans. The transfer-DNA molecule from Agrobacterium tumefaciens inserted into a gene encoding a protein with similarity to Sit4-associated proteins (SAPs). In contrast to Saccharomyces cerevisiae which has four members of the SAP family, there is a single copy of the gene in L. maculans. The mutant had normal spore production and spore germination, but altered hyphal branching, suggesting that nutrient signaling is impaired in the strain. This is the first time that a SAP gene has been mutated in a filamentous fungus and links the function of SAP proteins to plant pathogenesis and hyphal branching.

The Evolution of Orphan Regions in Genomes of a Fungal Pathogen of Wheat

Fungal plant pathogens rapidly evolve virulence on resistant hosts through mutations in genes encoding proteins that modulate the host immune responses. The mutational spectrum likely includes chromosomal rearrangements responsible for gains or losses of entire genes. However, the mechanisms creating adaptive structural variation in fungal pathogen populations are poorly understood. We used complete genome assemblies to quantify structural variants segregating in the highly polymorphic fungal wheat pathogen Zymoseptoria tritici The genetic basis of virulence in Z. tritici is complex, and populations harbor significant genetic variation for virulence; hence, we aimed to identify whether structural variation led to functional differences. We combined single-molecule real-time sequencing, genetic maps, and transcriptomics data to generate a fully assembled and annotated genome of the highly virulent field isolate 3D7. Comparative genomics analyses against the complete reference genome IPO323 identified large chromosomal inversions and the complete gain or loss of transposable-element clusters, explaining the extensive chromosomal-length polymorphisms found in this species. Both the 3D7 and IPO323 genomes harbored long tracts of sequences exclusive to one of the two genomes. These orphan regions contained 296 genes unique to the 3D7 genome and not previously known for this species. These orphan genes tended to be organized in clusters and showed evidence of mutational decay. Moreover, the orphan genes were enriched in genes encoding putative effectors and included a gene that is one of the most upregulated putative effector genes during wheat infection. Our study showed that this pathogen species harbored extensive chromosomal structure polymorphism that may drive the evolution of virulence.

IMPORTANCE

Pathogen outbreak populations often harbor previously unknown genes conferring virulence. Hence, a key puzzle of rapid pathogen evolution is the origin of such evolutionary novelty in genomes. Chromosomal rearrangements and structural variation in pathogen populations likely play a key role. However, identifying such polymorphism is challenging, as most genome-sequencing approaches only yield information about point mutations. We combined long-read technology and genetic maps to assemble the complete genome of a strain of a highly polymorphic fungal pathogen of wheat. Comparisons against the reference genome of the species showed substantial variation in the chromosome structure and revealed large regions unique to each assembled genome. These regions were enriched in genes encoding likely effector proteins, which are important components of pathogenicity. Our study showed that pathogen populations harbor extensive polymorphism at the chromosome level and that this polymorphism can be a source of adaptive genetic variation in pathogen evolution.

OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as "OcculterCut" (https://sourceforge.net/projects/occultercut, last accessed April 30, 2016), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighboring repeat-rich regions. Over time RIP perpetuates "two speeds" of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-891) contains supplementary material, which is available to authorized users.

Flippase (FLP) recombinase-mediated marker recycling in the dermatophyte Arthroderma vanbreuseghemii

Biological processes can be elucidated by investigating complex networks of relevant factors and genes. However, this is not possible in species for which dominant selectable markers for genetic studies are unavailable. To overcome the limitation in selectable markers for the dermatophyte Arthroderma vanbreuseghemii (anamorph: Trichophyton mentagrophytes), we adapted the flippase (FLP) recombinase-recombination target (FRT) site-specific recombination system from the yeast Saccharomyces cerevisiae as a selectable marker recycling system for this fungus. Taking into account practical applicability, we designed FLP/FRT modules carrying two FRT sequences as well as the flp gene adapted to the pathogenic yeast Candida albicans (caflp) or a synthetic codon-optimized flp (avflp) gene with neomycin resistance (nptII) cassette for one-step marker excision. Both flp genes were under control of the Trichophyton rubrum copper-repressible promoter (PCTR4). Molecular analyses of resultant transformants showed that only the avflp-harbouring module was functional in A. vanbreuseghemii. Applying this system, we successfully produced the Ku80 recessive mutant strain devoid of any selectable markers. This strain was subsequently used as the recipient for sequential multiple disruptions of secreted metalloprotease (fungalysin) (MEP) or serine protease (SUB) genes, producing mutant strains with double MEP or triple SUB gene deletions. These results confirmed the feasibility of this system for broad-scale genetic manipulation of dermatophytes, advancing our understanding of functions and networks of individual genes in these fungi.

The LmSNF1 gene is required for pathogenicity in the canola blackleg pathogen Leptosphaeria maculans

Leptosphaeria maculans is a fungal pathogen causing blackleg in canola. Its virulence has been attributed, among other factors, to the activity of hydrolytic cell wall degrading enzymes (CWDEs). Studies on the pathogenicity function of CWDEs in plant pathogenic fungi have been difficult due to gene redundancy. In microorganisms many CWDE genes are repressed by glucose and derepressed by the function of the sucrose non-fermenting protein kinase 1 gene (SNF1). To address the molecular function of SNF1 in L. maculans, the ortholog of SNF1 (LmSNF1) was cloned and functionally characterized using a gene knockout strategy. Growth of the LmSNF1 knockout strains was severely disrupted, as was sporulation, spore germination and the ability to attach on the plant surface. When inoculated on canola cotyledons, the LmSNF1 knockout strains could not cause any symptoms, indicating the loss of pathogenicity. The expression of 11 selected CWDE genes and a pathogenicity gene (LopB) was significantly down-regulated in the LmSNF1 knockout strains. In conclusion, knockout of LmSNF1 prevents L. maculans from properly derepressing the production of CWDEs, compromises the utilization of certain carbon sources, and impairs fungal pathogenicity on canola.

The transcriptional response to nonself in the fungus Podospora anserina

In fungi, heterokaryon incompatibility is a nonself recognition process occurring when filaments of different isolates of the same species fuse. Compatibility is controlled by so-called het loci and fusion of strains of unlike het genotype triggers a complex incompatibility reaction that leads to the death of the fusion cell. Herein, we analyze the transcriptional changes during the incompatibility reaction in Podospora anserina. The incompatibility response was found to be associated with a massive transcriptional reprogramming: 2231 genes were up-regulated by a factor 2 or more during incompatibility. In turn, 2441 genes were down-regulated. HET, NACHT, and HeLo domains previously found to be involved in the control of heterokaryon incompatibility were enriched in the up-regulated gene set. In addition, incompatibility was characterized by an up-regulation of proteolytic and other hydrolytic activities, of secondary metabolism clusters and toxins and effector-like proteins. The up-regulated set was found to be enriched for proteins lacking orthologs in other species and chromosomal distribution of the up-regulated genes was uneven with up-regulated genes residing preferentially in genomic islands and on chromosomes IV and V. There was a significant overlap between regulated genes during incompatibility in P. anserina and Neurospora crassa, indicating similarities in the incompatibility responses in these two species. Globally, this study illustrates that the expression changes occurring during cell fusion incompatibility in P. anserina are in several aspects reminiscent of those described in host-pathogen or symbiotic interactions in other fungal species.

The dispensable chromosome of Leptosphaeria maculans shelters an effector gene conferring avirulence towards Brassica rapa

Phytopathogenic fungi frequently contain dispensable chromosomes, some of which contribute to host range or pathogenicity. In Leptosphaeria maculans, the stem canker agent of oilseed rape (Brassica napus), the minichromosome was previously suggested to be dispensable, without evidence for any role in pathogenicity. Using genetic and genomic approaches, we investigated the inheritance and molecular determinant of an L. maculans-Brassica rapa incompatible interaction. Single gene control of the resistance was found, while all markers located on the L. maculans minichromosome, absent in the virulent parental isolate, co-segregated with the avirulent phenotype. Only one candidate avirulence gene was identified on the minichromosome, validated by complementation experiments and termed AvrLm11. The minichromosome was frequently lost following meiosis, but the frequency of isolates lacking it remained stable in field populations sampled at a 10-yr time interval, despite a yearly sexual stage in the L. maculans life cycle. This work led to the cloning of a new 'lost in the middle of nowhere' avirulence gene of L. maculans, interacting with a B. rapa resistance gene termed Rlm11 and introgressed into B. napus. It demonstrated the dispensability of the L. maculans minichromosome and suggested that its loss generates a fitness deficit.

Sex in cheese: evidence for sexuality in the fungus Penicillium roqueforti

Although most eukaryotes reproduce sexually at some moment of their life cycle, as much as a fifth of fungal species were thought to reproduce exclusively asexually. Nevertheless, recent studies have revealed the occurrence of sex in some of these supposedly asexual species. For industrially relevant fungi, for which inoculums are produced by clonal-subcultures since decades, the potentiality for sex is of great interest for strain improvement strategies. Here, we investigated the sexual capability of the fungus Penicillium roqueforti, used as starter for blue cheese production. We present indirect evidence suggesting that recombination could be occurring in this species. The screening of a large sample of strains isolated from diverse substrates throughout the world revealed the existence of individuals of both mating types, even in the very same cheese. The MAT genes, involved in fungal sexual compatibility, appeared to evolve under purifying selection, suggesting that they are still functional. The examination of the recently sequenced genome of the FM 164 cheese strain enabled the identification of the most important genes known to be involved in meiosis, which were found to be highly conserved. Linkage disequilibria were not significant among three of the six marker pairs and 11 out of the 16 possible allelic combinations were found in the dataset. Finally, the detection of signatures of repeat induced point mutations (RIP) in repeated sequences and transposable elements reinforces the conclusion that P. roqueforti underwent more or less recent sex events. In this species of high industrial importance, the induction of a sexual cycle would open the possibility of generating new genotypes that would be extremely useful to diversify cheese products.

Patterns of repeat-induced point mutation in transposable elements of basidiomycete fungi

Transposable elements (TEs) are ubiquitous genomic parasites that have prompted the evolution of genome defense systems that restrict their activity. Repeat-induced point mutation (RIP) is a homology-dependent genome defense that introduces C-to-T transition mutations in duplicated DNA sequences and is thought to control the proliferation of selfish repetitive DNA. Here, we determine the taxonomic distribution of hypermutation patterns indicative of RIP among basidiomycetes. We quantify C-to-T transition mutations in particular di- and trinucleotide target sites for TE-like sequences from nine fungal genomes. We find evidence of RIP-like patterns of hypermutation at TpCpG trinucleotide sites in repetitive sequences from all species of the Pucciniomycotina subphylum of the Basidiomycota, Microbotryum lychnidis-dioicae, Puccinia graminis, Melampsora laricis-populina, and Rhodotorula graminis. In contrast, we do not find evidence for RIP-like hypermutation in four species of the Agaricomycotina and Ustilaginomycotina subphyla of the Basidiomycota. Our results suggest that a RIP-like process and the specific nucleotide context for mutations are conserved within the Pucciniomycotina subphylum. These findings imply that coevolutionary interactions between TEs and a hypermutating genome defense are stable over long evolutionary timescales.

The [Het-s] prion of Podospora anserina and its role in heterokaryon incompatibility

[Het-s] is a prion from the filamentous fungus Podospora anserina and corresponds to a self-perpetuating amyloid aggregate of the HET-s protein. This prion protein is involved in a fungal self/non-self discrimination process termed heterokaryon incompatibility corresponding to a cell death reaction occurring upon fusion of genetically unlike strains. Two antagonistic allelic variants of this protein exist: HET-s, the prion form of which corresponds to [Het-s] and HET-S, incapable of prion formation. Fusion of a [Het-s] and HET-S strain triggers the incompatibility reaction, so that interaction of HET-S with the [Het-s] prion leads to cell death. HET-s and HET-S are highly homologous two domain proteins with a N-terminal globular domain termed HeLo and a C-terminal unstructured prion forming domain (PFD). The structure of the prion form of the HET-s PFD has been solved by solid state NMR and corresponds to a very well ordered β-solenoid fold with a triangular hydrophobic core. The ability to form this β-solenoid fold is retained in a distant homolog of HET-s from another fungal species. A model for the mechanism of [Het-s]/HET-S incompatibility has been proposed. It is believe that when interacting with the [Het-s] prion seed, the HET-S C-terminal region adopts the β-solenoid fold. This would act as a conformational switch to induce refolding and activation of the HeLo domain which then would exert its toxicity by a yet unknown mechanism.

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.

In silico reversal of repeat-induced point mutation (RIP) identifies the origins of repeat families and uncovers obscured duplicated genes

Background

Repeat-induced point mutation (RIP) is a fungal genome defence mechanism guarding against transposon invasion. RIP mutates the sequence of repeated DNA and over time renders the affected regions unrecognisable by similarity search tools such as BLAST.

Results

DeRIP is a new software tool developed to predict the original sequence of a RIP-mutated region prior to the occurrence of RIP. In this study, we apply deRIP to the genome of the wheat pathogen Stagonospora nodorum SN15 and predict the origin of several previously uncharacterised classes of repetitive DNA.

Conclusions

Five new classes of transposon repeats and four classes of endogenous gene repeats were identified after deRIP. The deRIP process is a new tool for fungal genomics that facilitates the identification and understanding of the role and origin of fungal repetitive DNA. DeRIP is open-source and is available as part of the RIPCAL suite at http://www.sourceforge.net/projects/ripcal.

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.

De novo assembly of a 40 Mb eukaryotic genome from short sequence reads: Sordaria macrospora, a model organism for fungal morphogenesis

Filamentous fungi are of great importance in ecology, agriculture, medicine, and biotechnology. Thus, it is not surprising that genomes for more than 100 filamentous fungi have been sequenced, most of them by Sanger sequencing. While next-generation sequencing techniques have revolutionized genome resequencing, e.g. for strain comparisons, genetic mapping, or transcriptome and ChIP analyses, de novo assembly of eukaryotic genomes still presents significant hurdles, because of their large size and stretches of repetitive sequences. Filamentous fungi contain few repetitive regions in their 30-90 Mb genomes and thus are suitable candidates to test de novo genome assembly from short sequence reads. Here, we present a high-quality draft sequence of the Sordaria macrospora genome that was obtained by a combination of Illumina/Solexa and Roche/454 sequencing. Paired-end Solexa sequencing of genomic DNA to 85-fold coverage and an additional 10-fold coverage by single-end 454 sequencing resulted in approximately 4 Gb of DNA sequence. Reads were assembled to a 40 Mb draft version (N50 of 117 kb) with the Velvet assembler. Comparative analysis with Neurospora genomes increased the N50 to 498 kb. The S. macrospora genome contains even fewer repeat regions than its closest sequenced relative, Neurospora crassa. Comparison with genomes of other fungi showed that S. macrospora, a model organism for morphogenesis and meiosis, harbors duplications of several genes involved in self/nonself-recognition. Furthermore, S. macrospora contains more polyketide biosynthesis genes than N. crassa. Phylogenetic analyses suggest that some of these genes may have been acquired by horizontal gene transfer from a distantly related ascomycete group. Our study shows that, for typical filamentous fungi, de novo assembly of genomes from short sequence reads alone is feasible, that a mixture of Solexa and 454 sequencing substantially improves the assembly, and that the resulting data can be used for comparative studies to address basic questions of fungal biology.

Repeat-induced point mutation (RIP) as an alternative mechanism of evolution toward virulence in Leptosphaeria maculans

Three avirulence genes, AvrLm1, AvrLm6, and AvrLm4-7, were recently identified in Leptosphaeria maculans and found to be localized as solo genes within large noncoding, heterochromatin-like regions mainly composed of retrotransposons, truncated and degenerated by repeat-induced point mutation (RIP). The Rlm6 resistance gene has been overcome within 3 years in outdoor experiments in France and, here, we investigate the molecular basis of evolution toward virulence at the AvrLm6 locus. A region of 235 kb was sequenced in a virulent isolate and showed the deletion of AvrLm6 and three divergent mosaics of retrotransposons. AvrLm6 was found to be absent from 66% of 70 virulent isolates, with multiple events of deletion. The sequencing of virulent alleles in 24 isolates revealed a few cases of point mutations that had created stop codons in the sequence. The most frequent mutation events, however, were RIP, leading to the modification of 4 to 9% of the bases compared with the avirulent allele and generating 2 to 4 stop codons. Thus, RIP is described for the first time as an efficient mechanism leading to virulence and the multiple patterns of mutation observed suggest that multiple RIP events could occur independently in a single field population during 1 year.

RIPCAL: a tool for alignment-based analysis of repeat-induced point mutations in fungal genomic sequences

Background

Repeat-induced point mutation (RIP) is a fungal-specific genome defence mechanism that alters the sequences of repetitive DNA, thereby inactivating coding genes. Repeated DNA sequences align between mating and meiosis and both sequences undergo C:G to T:A transitions. In most fungi these transitions preferentially affect CpA di-nucleotides thus altering the frequency of certain di-nucleotides in the affected sequences. The majority of previously published in silico analyses were limited to the comparison of ratios of pre- and post-RIP di-nucleotides in putatively RIP-affected sequences – so-called RIP indices. The analysis of RIP is significantly more informative when comparing sequence alignments of repeated sequences. There is, however, a dearth of bioinformatics tools available to the fungal research community for alignment-based RIP analysis of repeat families.

Results

We present RIPCAL , a software tool for the automated analysis of RIP in fungal genomic DNA repeats, which performs both RIP index and alignment-based analyses. We demonstrate the ability of RIPCAL to detect RIP within known RIP-affected sequences of Neurospora crassa and other fungi. We also predict and delineate the presence of RIP in the genome of Stagonospora nodorum – a Dothideomycete pathogen of wheat. We show that RIP has affected different members of the S. nodorum rDNA tandem repeat to different extents depending on their genomic contexts.

Conclusion

The RIPCAL alignment-based method has considerable advantages over RIP indices for the analysis of whole genomes. We demonstrate its application to the recently published genome assembly of S. nodorum.

The Lmpma1 gene of Leptosphaeria maculans encodes a plasma membrane H+-ATPase isoform essential for pathogenicity towards oilseed rape

Following Agrobacterium tumefaciens-mediated mutagenesis in Leptosphaeria maculans, we identified the mutant 210, displaying total loss of pathogenicity towards its host plant (Brassica napus). Microscopic observations showed that m210 is unable to germinate on the host leaf surface and is thus blocked at the pre-penetration stage. The pathogenicity phenotype is linked with a single T-DNA insertion into the promoter region of a typical plasma membrane H(+)-ATPase-encoding gene, termed Lmpma1, thus leading to a twofold reduction in Lmpma1 expression. Since LmPMA1 is involved in intracellular pH homeostasis, we postulate that reduction in LmPMA1 activity disturbs the electrochemical transmembrane gradient in m210, thus leading to conidia defective in turgor pressure generation on leaf surface. Whole genome survey showed that L. maculans possesses a second plasma membrane H(+)-ATPase-encoding gene, termed Lmpma2. Silencing experiments, expression analyses and phylogenetic studies allowed us to highlight the essential role assumed by the Lmpma1 isoform in L.maculans pathogenicity.

The Lmgpi15 gene, encoding a component of the glycosylphosphatidylinositol anchor biosynthesis pathway, is required for morphogenesis and pathogenicity in Leptosphaeria maculans

Random insertional mutagenesis was used to investigate pathogenicity determinants in Leptosphaeria maculans. One tagged nonpathogenic mutant, termed m20, was analysed in detail here. The mutant phenotype was investigated by microscopic analyses of infected plant tissues and in vitro growth assays. Complementation and silencing experiments were used to identify the altered gene. Its function was determined by bioinformatics analyses, cell biology experiments and functional studies. The mutant was blocked at the invasive growth phase after an unaffected initial penetration stage, and displayed a reduced growth rate and an aberrant hyphal morphology in vitro. The T-DNA insertion occurred in the intergenic region between two head-to-tail genes, leading to a complex deregulation of their expression. The unique gene accounting for the mutant phenotype was suggested to be the orthologue of the poorly conserved Saccharomyces cerevisiae gpi15, which encodes for one component of the glycosylphosphatidylinositol (GPI) anchor biosynthesis pathway. Consistent with this predicted function, a functional translational fusion with the green fluorescent protein (GFP) was targeted to the endoplasmic reticulum. Moreover, the mutant exhibited an altered cell wall and addition of glucosamine relieved growth defects. It is concluded that the GPI anchor biosynthetic pathway is required for morphogenesis, cell wall integrity and pathogenicity in Leptosphaeria maculans.

A cytosine methyltransferase homologue is essential for sexual development in Aspergillus nidulans

BACKGROUND

The genome defense processes RIP (repeat-induced point mutation) in the filamentous fungus Neurospora crassa, and MIP (methylation induced premeiotically) in the fungus Ascobolus immersus depend on proteins with DNA methyltransferase (DMT) domains. Nevertheless, these proteins, RID and Masc1, respectively, have not been demonstrated to have DMT activity. We discovered a close homologue in Aspergillus nidulans, a fungus thought to have no methylation and no genome defense system comparable to RIP or MIP.

PRINCIPAL FINDINGS

We report the cloning and characterization of the DNA methyltransferase homologue A (dmtA) gene from Aspergillus nidulans. We found that the dmtA locus encodes both a sense (dmtA) and an anti-sense transcript (tmdA). Both transcripts are expressed in vegetative, conidial and sexual tissues. We determined that dmtA, but not tmdA, is required for early sexual development and formation of viable ascospores. We also tested if DNA methylation accumulated in any of the dmtA/tmdA mutants we constructed, and found that in both asexual and sexual tissues, these mutants, just like wild-type strains, appear devoid of DNA methylation.

CONCLUSIONS/SIGNIFICANCE

Our results demonstrate that a DMT homologue closely related to proteins implicated in RIP and MIP has an essential developmental function in a fungus that appears to lack both DNA methylation and RIP or MIP. It remains formally possible that DmtA is a bona fide DMT, responsible for trace, undetected DNA methylation that is restricted to a few cells or transient but our work supports the idea that the DMT domain present in the RID/Masc1/DmtA family has a previously undescribed function.

Repeat induced point mutation in two asexual fungi, Aspergillus niger and Penicillium chrysogenum

Repeat induced point mutation (RIP) is a gene silencing mechanism present in fungal genomes. During RIP, duplicated sequences are efficiently and irreversibly mutated by transitions from C:G to T:A. For the first time, we have identified traces of RIP in transposable elements of Aspergillus niger and Penicillium chrysogenum, two biotechnologically relevant fungi. We found that RIP in P. chrysogenum has affected a large set of sequences, which also contain other mutations. On the other hand, RIP in A. niger is limited to only few sequences, but literally all mutations are RIP-like. Surprisingly, RIP occurred only in transposon sequences that have disrupted open reading frames in A. niger, a phenomenon not yet reported for other fungi. In both fungal species, we identified two sequences with strong sequence similarity to Neurospora crassa RID. RID is a putative DNA methyltransferase and the only known enzyme involved in the RIP process. Our findings suggest that both A. niger and P. chrysogenum either had a sexual past or have a sexual potential. These findings have important implications for future strain development of these fungi.

Dothideomycete plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum

Stagonospora nodorum is a major necrotrophic fungal pathogen of wheat (Triticum aestivum) and a member of the Dothideomycetes, a large fungal taxon that includes many important plant pathogens affecting all major crop plant families. Here, we report the acquisition and initial analysis of a draft genome sequence for this fungus. The assembly comprises 37,164,227 bp of nuclear DNA contained in 107 scaffolds. The circular mitochondrial genome comprises 49,761 bp encoding 46 genes, including four that are intron encoded. The nuclear genome assembly contains 26 classes of repetitive DNA, comprising 4.5% of the genome. Some of the repeats show evidence of repeat-induced point mutations consistent with a frequent sexual cycle. ESTs and gene prediction models support a minimum of 10,762 nuclear genes. Extensive orthology was found between the polyketide synthase family in S. nodorum and Cochliobolus heterostrophus, suggesting an ancient origin and conserved functions for these genes. A striking feature of the gene catalog was the large number of genes predicted to encode secreted proteins; the majority has no meaningful similarity to any other known genes. It is likely that genes for host-specific toxins, in addition to ToxA, will be found among this group. ESTs obtained from axenic mycelium grown on oleate (chosen to mimic early infection) and late-stage lesions sporulating on wheat leaves were obtained. Statistical analysis shows that transcripts encoding proteins involved in protein synthesis and in the production of extracellular proteases, cellulases, and xylanases predominate in the infection library. This suggests that the fungus is dependant on the degradation of wheat macromolecular constituents to provide the carbon skeletons and energy for the synthesis of proteins and other components destined for the developing pycnidiospores.

A critical assessment of Agrobacterium tumefaciens-mediated transformation as a tool for pathogenicity gene discovery in the phytopathogenic fungus Leptosphaeria maculans

We evaluated the usefulness and robustness of Agrobacterium tumefaciens-mediated transformation (ATMT) as a high-throughput transformation tool for pathogenicity gene discovery in the filamentous phytopathogen Leptosphaeria maculans. Thermal asymmetric interlaced polymerase chain reaction allowed us to amplify the left border (LB) flanking sequence in 135 of 400 transformants analysed, and indicated a high level of preservation of the T-DNA LB. In addition, T-DNA preferentially integrated as a single copy in gene-rich regions of the fungal genome, with a probable bias towards intergenic and/or regulatory regions. A total of 53 transformants out of 1388 (3.8%) showed reproducible pathogenicity defects when inoculated on cotyledons of Brassica napus, with diverse altered phenotypes. Co-segregation of the altered phenotype with the T-DNA integration was observed for 6 of 12 transformants crossed. If extrapolated to the whole collection, this indicates that 1.9% of the collection actually corresponds to tagged pathogenicity mutants. The preferential insertion into gene-rich regions along with the high ratio of tagged mutants renders ATMT a tool of choice for large-scale gene discovery in 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.

Molecular cloning and genetic analysis of a symbiosis-expressed gene cluster for lolitrem biosynthesis from a mutualistic endophyte of perennial ryegrass

Lolitrems are potent tremorgenic mycotoxins that are synthesised by clavicipitaceous fungal endophytes of the Epichloë/Neotyphodium group in association with grasses. These indole-diterpenes confer major ecological benefits on the grass-endophyte symbiotum. A molecular signature for diterpene biosynthesis is the presence of two geranylgeranyl diphosphate (GGPP) synthases. Using degenerate primers for conserved domains of fungal GGPP synthases, we cloned two such genes, ltmG and ggsA, from Neotyphodium lolii. Adjacent to ltmG are two genes, ltmM and ltmK, that are predicted to encode an FAD-dependent monooxygenase and a cytochrome P450 monooxygenase, respectively. The cluster of ltm genes is flanked by AT-rich retrotransposon DNA that appears to have undergone extensive repeat induced point (RIP) mutation. Epichloë festucae, the sexual ancestor of N. lolii, contains an identical ltm gene cluster, but lacks the retrotransposon "platform'' on the right flank. Associations established between perennial ryegrass and an E. festucae mutant deleted for ltmM lack detectable levels of lolitrems. A wild-type copy of ltmM complemented this phenotype, as did paxM from Penicillium paxilli. Northern hybridization and RT-PCR analysis showed that all three genes are weakly expressed in culture but strongly induced in planta. The relative endophyte biomass in these associations was estimated by real-time PCR to be between 0.3 and 1.9%. Taking this difference into account, the steady-state levels of the ltm transcripts are about 100-fold greater than the levels of the endogenous ryegrass beta-tubulin (beta -Tub1) and actin (Act1) RNAs. Based on these results we propose that ltmG, ltmM and ltmK are members of a set of genes required for lolitrem biosynthesis in E. festucae and N. lolii.

The stem canker (blackleg) fungus, Leptosphaeria maculans, enters the genomic era

SUMMARY Leptosphaeria maculans is the most ubiquitous pathogen of Brassica crops, and mainly oilseed brassicas (oilseed rape, canola), causing the devastating 'stem canker' or 'blackleg'. This review summarizes our current knowledge on the pathogen, from taxonomic issues to specific life traits. It mainly illustrates the importance of formal genetics approaches on the pathogen side to dissect the interaction with the host plants. In addition, this review presents the main current research topics on L. maculans and focuses on the L. maculans genome initiative recently begun, including its main research issues.

TAXONOMY

Leptosphaeria maculans (Desm.) Ces. & de Not. (anamorph Phoma lingam Tode ex Fr.). Kingdom Fungi, Phylum Ascomycota, Class Dothideomycetes (Loculoascomycetes), Order Pleosporales, Genus Leptosphaeria, Species maculans.

HOST RANGE

cultivated Brassicas such as Brassica napus (oilseed rape, canola), B. rapa, B. juncea, B. oleracea, etc., along with numerous wild crucifers species. Arabidopsis thaliana was recently reported to be a potential host for L. maculans. Primary disease symptoms are greyish-green collapse of cotyledon or leaf tissue, without a visible margin, bearing tiny black spots (pycnidia). The fungus then develops an endophytic symptomless growth for many months. Secondary symptoms, at the end of the growing season, are dry necroses of the crown tissues with occasional blackening (stem canker or blackleg) causing lodging of the plants. Pseudothecia differentiate on leftover residues. Seedling damping-off and premature ripening are also reported under certain environmental conditions.

USEFUL WEBSITES

Leptosphaeria maculans sequencing project at Genoscope: http://www.genoscope.cns.fr/externe/English/Projets/Projet_DM/organisme_DM.html; the SECURE site: http://www.secure.rothamsted.ac.uk/ the 'Blackleg' group at the University of Melbourne: http://www.botany.unimelb.edu.au/blackleg/overview.htm.

The sirodesmin biosynthetic gene cluster of the plant pathogenic fungus Leptosphaeria maculans

Sirodesmin PL is a phytotoxin produced by the fungus Leptosphaeria maculans, which causes blackleg disease of canola (Brassica napus). This phytotoxin belongs to the epipolythiodioxopiperazine (ETP) class of toxins produced by fungi including mammalian and plant pathogens. We report the cloning of a cluster of genes with predicted roles in the biosynthesis of sirodesmin PL and show via gene disruption that one of these genes (encoding a two-module non-ribosomal peptide synthetase) is essential for sirodesmin PL biosynthesis. Of the nine genes in the cluster tested, all are co-regulated with the production of sirodesmin PL in culture. A similar cluster is present in the genome of the opportunistic human pathogen Aspergillus fumigatus and is most likely responsible for the production of gliotoxin, which is also an ETP. Homologues of the genes in the cluster were also identified in expressed sequence tags of the ETP producing fungus Chaetomium globosum. Two other fungi with publicly available genome sequences, Magnaporthe grisea and Fusarium graminearum, had similar gene clusters. A comparative analysis of all four clusters is presented. This is the first report of the genes responsible for the biosynthesis of an ETP.