Sequential breakdown of the Cf-9 leaf mould resistance locus in tomato by Fulvia fulva.

Leaf mould, caused by Fulvia fulva, is a devastating disease of tomato plants. In many commercial tomato cultivars, resistance to this disease is governed by the Cf-9 locus, which encodes five paralogous receptor-like proteins. Two of these proteins confer resistance: Cf-9C recognises the previously identified F. fulva effector Avr9 and provides resistance during all plant growth stages, while Cf-9B recognises the yet-unidentified F. fulva effector Avr9B and provides mature plant resistance only. In recent years, F. fulva strains have emerged that can overcome the Cf-9 locus, with Cf-9C circumvented through Avr9 deletion. To understand how Cf-9B is circumvented, we set out to identify Avr9B. Comparative genomics, transient expression assays and gene complementation experiments were used to identify Avr9B, while gene sequencing was used to assess Avr9B allelic variation across a world-wide strain collection. A strict correlation between Avr9 deletion and resistance-breaking mutations in Avr9B was observed in strains recently collected from Cf-9 cultivars, whereas Avr9 deletion but no mutations in Avr9B were observed in older strains. This research showcases how F. fulva has evolved to sequentially break down the Cf-9 locus and stresses the urgent need for commercial tomato cultivars that carry novel, stacked resistance genes active against this pathogen.

Elucidation of cladofulvin biosynthesis reveals a cytochrome P450 monooxygenase required for anthraquinone dimerization.

Anthraquinones are a large family of secondary metabolites (SMs) that are extensively studied for their diverse biological activities. These activities are determined by functional group decorations and the formation of dimers from anthraquinone monomers. Despite their numerous medicinal qualities, very few anthraquinone biosynthetic pathways have been elucidated so far, including the enzymatic dimerization steps. In this study, we report the elucidation of the biosynthesis of cladofulvin, an asymmetrical homodimer of nataloe-emodin produced by the fungus Cladosporium fulvum A gene cluster of 10 genes controls cladofulvin biosynthesis, which begins with the production of atrochrysone carboxylic acid by the polyketide synthase ClaG and the ß-lactamase ClaF. This compound is decarboxylated by ClaH to yield emodin, which is then converted to chrysophanol hydroquinone by the reductase ClaC and the dehydratase ClaB. We show that the predicted cytochrome P450 ClaM catalyzes the dimerization of nataloe-emodin to cladofulvin. Remarkably, such dimerization dramatically increases nataloe-emodin cytotoxicity against mammalian cell lines. These findings shed light on the enzymatic mechanisms involved in anthraquinone dimerization. Future characterization of the ClaM enzyme should facilitate engineering the biosynthesis of novel, potent, dimeric anthraquinones and structurally related compound families.

Combating a Global Threat to a Clonal Crop: Banana Black Sigatoka Pathogen Pseudocercospora fijiensis (Synonym Mycosphaerella fijiensis) Genomes Reveal Clues for Disease Control.

Black Sigatoka or black leaf streak disease, caused by the Dothideomycete fungus Pseudocercospora fijiensis (previously: Mycosphaerella fijiensis), is the most significant foliar disease of banana worldwide. Due to the lack of effective host resistance, management of this disease requires frequent fungicide applications, which greatly increase the economic and environmental costs to produce banana. Weekly applications in most banana plantations lead to rapid evolution of fungicide-resistant strains within populations causing disease-control failures throughout the world. Given its extremely high economic importance, two strains of P. fijiensis were sequenced and assembled with the aid of a new genetic linkage map. The 74-Mb genome of P. fijiensis is massively expanded by LTR retrotransposons, making it the largest genome within the Dothideomycetes. Melting-curve assays suggest that the genomes of two closely related members of the Sigatoka disease complex, P. eumusae and P. musae, also are expanded. Electrophoretic karyotyping and analyses of molecular markers in P. fijiensis field populations showed chromosome-length polymorphisms and high genetic diversity. Genetic differentiation was also detected using neutral markers, suggesting strong selection with limited gene flow at the studied geographic scale. Frequencies of fungicide resistance in fungicide-treated plantations were much higher than those in untreated wild-type P. fijiensis populations. A homologue of the Cladosporium fulvum Avr4 effector, PfAvr4, was identified in the P. fijiensis genome. Infiltration of the purified PfAVR4 protein into leaves of the resistant banana variety Calcutta 4 resulted in a hypersensitive-like response. This result suggests that Calcutta 4 could carry an unknown resistance gene recognizing PfAVR4. Besides adding to our understanding of the overall Dothideomycete genome structures, the P. fijiensis genome will aid in developing fungicide treatment schedules to combat this pathogen and in improving the efficiency of banana breeding programs.

Specific Hypersensitive Response-Associated Recognition of New Apoplastic Effectors from Cladosporium fulvum in Wild Tomato.

Tomato leaf mold disease is caused by the biotrophic fungus Cladosporium fulvum. During infection, C. fulvum produces extracellular small secreted protein (SSP) effectors that function to promote colonization of the leaf apoplast. Resistance to the disease is governed by Cf immune receptor genes that encode receptor-like proteins (RLPs). These RLPs recognize specific SSP effectors to initiate a hypersensitive response (HR) that renders the pathogen avirulent. C. fulvum strains capable of overcoming one or more of all cloned Cf genes have now emerged. To combat these strains, new Cf genes are required. An effectoromics approach was employed to identify wild tomato accessions carrying new Cf genes. Proteomics and transcriptome sequencing were first used to identify 70 apoplastic in planta-induced C. fulvum SSPs. Based on sequence homology, 61 of these SSPs were novel or lacked known functional domains. Seven, however, had predicted structural homology to antimicrobial proteins, suggesting a possible role in mediating antagonistic microbe-microbe interactions in planta. Wild tomato accessions were then screened for HR-associated recognition of 41 SSPs, using the Potato virus X-based transient expression system. Nine SSPs were recognized by one or more accessions, suggesting that these plants carry new Cf genes available for incorporation into cultivated tomato.

SnTox1, a Parastagonospora nodorum necrotrophic effector, is a dual-function protein that facilitates infection while protecting from wheat-produced chitinases.

SnTox1 induces programmed cell death and the up-regulation of pathogenesis-related genes including chitinases. Additionally, SnTox1 has structural homology to several plant chitin-binding proteins. Therefore, we evaluated SnTox1 for chitin binding and localization. We transformed an avirulent strain of Parastagonospora nodorum as well as three nonpathogens of wheat (Triticum aestivum), including a necrotrophic pathogen of barley, a hemibiotrophic pathogen of sugar beet and a saprotroph, to evaluate the role of SnTox1 in infection and in protection from wheat chitinases. SnTox1 bound chitin and an SnTox1-green fluorescent fusion protein localized to the mycelial cell wall. Purified SnTox1 induced necrosis in the absence of the pathogen when sprayed on the leaf surface and appeared to remain on the leaf surface while inducing both epidermal and mesophyll cell death. SnTox1 protected the different fungi from chitinase degradation. SnTox1 was sufficient to change the host range of a necrotrophic pathogen but not a hemibiotroph or saprotroph. Collectively, this work shows that SnTox1 probably interacts with a receptor on the outside of the cell to induce cell death to acquire nutrients, but SnTox1 accomplishes a second role in that it protects against one aspect of the defense response, namely the effects of wheat chitinases.

Synergistic Action of a Metalloprotease and a Serine Protease from Fusarium oxysporum f. sp. lycopersici Cleaves Chitin-Binding Tomato Chitinases, Reduces Their Antifungal Activity, and Enhances Fungal Virulence.

As part of their defense strategy against fungal pathogens, plants secrete chitinases that degrade chitin, the major structural component of fungal cell walls. Some fungi are not sensitive to plant chitinases because they secrete chitin-binding effector proteins that protect their cell wall against these enzymes. However, it is not known how fungal pathogens that lack chitin-binding effectors overcome this plant defense barrier. Here, we investigated the ability of fungal tomato pathogens to cleave chitin-binding domain (CBD)-containing chitinases and its effect on fungal virulence. Four tomato CBD chitinases were produced in Pichia pastoris and were incubated with secreted proteins isolated from seven fungal tomato pathogens. Of these, Fusarium oxysporum f. sp. lycopersici, Verticillium dahliae, and Botrytis cinerea were able to cleave the extracellular tomato chitinases SlChi1 and SlChi13. Cleavage by F. oxysporum removed the CBD from the N-terminus, shown by mass spectrometry, and significantly reduced the chitinase and antifungal activity of both chitinases. Both secreted metalloprotease FoMep1 and serine protease FoSep1 were responsible for this cleavage. Double deletion mutants of FoMep1 and FoSep1 of F. oxysporum lacked chitinase cleavage activity on SlChi1 and SlChi13 and showed reduced virulence on tomato. These results demonstrate the importance of plant chitinase cleavage in fungal virulence.

Functional analysis of the conserved transcriptional regulator CfWor1 in Cladosporium fulvum reveals diverse roles in the virulence of plant pathogenic fungi.

Fungal Wor1-like proteins are conserved transcriptional regulators that are reported to regulate the virulence of several plant pathogenic fungi by affecting the expression of virulence genes. Here, we report the functional analysis of CfWor1, the homologue of Wor1 in Cladosporium fulvum. Δcfwor1 mutants produce sclerotium-like structures and rough hyphae, which are covered with a black extracellular matrix. These mutants do not sporulate and are no longer virulent on tomato. A CE.CfWor1 transformant that constitutively expresses CfWor1 produces fewer spores with altered morphology and is also reduced in virulence. RNA-seq and RT-qrtPCR analyses suggest that reduced virulence of Δcfwor1 mutants is due to global downregulation of transcription, translation and mitochondrial respiratory chain. The reduced virulence of the CE.CfWor1 transformant is likely due to downregulation of effector genes. Complementation of a non-virulent Δfosge1 (Wor1-homologue) mutant of Fusarium oxysporum f. sp. lycopersici with CfWor1 restored expression of the SIX effector genes in this fungus, but not its virulence. Chimeric proteins of CfWor1/FoSge1 also only partially restored defects of the Δfosge1 mutant, suggesting that these transcriptional regulators have functionally diverged. Altogether, our results suggest that CfWor1 primarily regulates development of C. fulvum, which indirectly affects the expression of a subset of virulence genes.

Regulation of secondary metabolite production in the fungal tomato pathogen Cladosporium fulvum.

Cladosporium fulvum is a non-obligate biotrophic fungal tomato pathogen for which fifteen secondary metabolite (SM) gene clusters were previously identified in its genome. However, most of these SM biosynthetic pathways remain cryptic during growth in planta and in different in vitro conditions. The sole SM produced in vitro is the pigment cladofulvin. In this study, we attempted to activate cryptic pathways in order to identify new compounds produced by C. fulvum. For this purpose, we manipulated orthologues of the global regulators VeA, LaeA and HdaA known to regulate SM biosynthesis in other fungal species. In C. fulvum, deleting or over-expressing these regulators yielded no new detectable SMs. Yet, quantification of cladofulvin revealed that CfHdaA is an activator whilst CfVeA and CfLaeA seemed to act as repressors of cladofulvin production. In the wild type strain, cladofulvin biosynthesis was affected by the carbon source, with highest production under carbon limitation and traces only in presence of saccharose. Repression of cladofulvin production by saccharose was dependent on both CfVeA and CfLaeA. Deletion of CfVeA or CfLaeA caused production of sterile mycelia, whilst Δcfhdaa deletion mutants sporulated, suggesting that cladofulvin production is not linked to asexual reproduction. Profiling the transcription of these regulators showed that CfHdaA-mediated regulation of cladofulvin production is independent of both CfVeA and CfLaeA. Our data suggest CfLaeA directly affects cladofulvin production whilst the effect of CfVeA is indirect, suggesting a role for CfLaeA outside of the Velvet complex. In conclusion, our results showed that regulation of SM production in C. fulvum is different from other fungi and indicate that manipulation of global regulators is not a universal tool to discover new fungal natural products.

Flexible gateway constructs for functional analyses of genes in plant pathogenic fungi.

Genetic manipulation of fungi requires quick, low-cost, efficient, high-throughput and molecular tools. In this paper, we report 22 entry constructs as new molecular tools based on the Gateway technology facilitating rapid construction of binary vectors that can be used for functional analysis of genes in fungi. The entry vectors for single, double or triple gene-deletion mutants were developed using hygromycin, geneticin and nourseothricin resistance genes as selection markers. Furthermore, entry vectors containing green fluorescent (GFP) or red fluorescent (RFP) in combination with hygromycin, geneticin or nourseothricin selection markers were generated. The latter vectors provide the possibility of gene deletion and simultaneous labelling of the fungal transformants with GFP or RFP reporter genes. The applicability of a number of entry vectors was validated in Zymoseptoria tritici, an important fungal wheat pathogen.

Proteome catalog of Zymoseptoria tritici captured during pathogenesis in wheat.

Zymoseptoria tritici is an economically important pathogen of wheat. However, the molecular basis of pathogenicity on wheat is still poorly understood. Here, we present a global survey of the proteins secreted by this fungus in the apoplast of resistant (cv. Shafir) and susceptible (cv. Obelisk) wheat cultivars after inoculation with reference Z. tritici strain IPO323. The fungal proteins present in apoplastic fluids were analyzed by gel electrophoresis and by data-independent acquisition liquid chromatography/mass spectrometry (LC/MS(E)) combined with data-dependent acquisition LC-MS/MS. Subsequent mapping mass spectrometry-derived peptide sequence data against the genome sequence of strain IPO323 identified 665 peptides in the MS(E) and 93 in the LC-MS/MS mode that matched to 85 proteins. The identified fungal proteins, including cell-wall degrading enzymes and proteases, might function in pathogenicity, but the functions of many remain unknown. Most fungal proteins accumulated in cv. Obelisk at the onset of necrotrophy. This inventory provides an excellent basis for future detailed studies on the role of these genes and their encoded proteins during pathogenesis in wheat.

FPLC and liquid-chromatography mass spectrometry identify candidate necrosis-inducing proteins from culture filtrates of the fungal wheat pathogen Zymoseptoria tritici.

Culture filtrates (CFs) of the fungal wheat pathogen Zymoseptoria tritici were assayed for necrosis-inducing activity after infiltration in leaves of various wheat cultivars. Active fractions were partially purified and characterized. The necrosis-inducing factors in CFs are proteinaceous, heat stable and their necrosis-inducing activity is temperature and light dependent. The in planta activity of CFs was tested by a time series of proteinase K (PK) co-infiltrations, which was unable to affect activity 30min after CF infiltrations. This suggests that the necrosis inducing proteins (NIPs) are either absent from the apoplast and likely actively transported into mesophyll cells or protected from the protease by association with a receptor. Alternatively, plant cell death signaling pathways might be fully engaged during the first 30min and cannot be reversed even after PK treatment. Further fractionation of the CFs with the highest necrosis-inducing activity involved fast performance liquid chromatography, SDS-PAGE and mass spectrometry. This revealed that most of the proteins present in the fractions have not been described before. The two most prominent ZtNIP encoding candidates were heterologously expressed in Pichia pastoris and subsequent infiltration assays showed their differential activity in a range of wheat cultivars.

Dual disease resistance mediated by the immune receptor Cf-2 in tomato requires a common virulence target of a fungus and a nematode.

Plants lack the seemingly unlimited receptor diversity of a somatic adaptive immune system as found in vertebrates and rely on only a relatively small set of innate immune receptors to resist a myriad of pathogens. Here, we show that disease-resistant tomato plants use an efficient mechanism to leverage the limited nonself recognition capacity of their innate immune system. We found that the extracellular plant immune receptor protein Cf-2 of the red currant tomato (Solanum pimpinellifolium) has acquired dual resistance specificity by sensing perturbations in a common virulence target of two independently evolved effectors of a fungus and a nematode. The Cf-2 protein, originally identified as a monospecific immune receptor for the leaf mold fungus Cladosporium fulvum, also mediates disease resistance to the root parasitic nematode Globodera rostochiensis pathotype Ro1-Mierenbos. The Cf-2-mediated dual resistance is triggered by effector-induced perturbations of the apoplastic Rcr3(pim) protein of S. pimpinellifolium. Binding of the venom allergen-like effector protein Gr-VAP1 of G. rostochiensis to Rcr3(pim) perturbs the active site of this papain-like cysteine protease. In the absence of the Cf-2 receptor, Rcr3(pim) increases the susceptibility of tomato plants to G. rostochiensis, thus showing its role as a virulence target of these nematodes. Furthermore, both nematode infection and transient expression of Gr-VAP1 in tomato plants harboring Cf-2 and Rcr3(pim) trigger a defense-related programmed cell death in plant cells. Our data demonstrate that monitoring host proteins targeted by multiple pathogens broadens the spectrum of disease resistances mediated by single plant immune receptors.

The genomes of the fungal plant pathogens Cladosporium fulvum and Dothistroma septosporum reveal adaptation to different hosts and lifestyles but also signatures of common ancestry.

We sequenced and compared the genomes of the Dothideomycete fungal plant pathogens Cladosporium fulvum (Cfu) (syn. Passalora fulva) and Dothistroma septosporum (Dse) that are closely related phylogenetically, but have different lifestyles and hosts. Although both fungi grow extracellularly in close contact with host mesophyll cells, Cfu is a biotroph infecting tomato, while Dse is a hemibiotroph infecting pine. The genomes of these fungi have a similar set of genes (70% of gene content in both genomes are homologs), but differ significantly in size (Cfu >61.1-Mb; Dse 31.2-Mb), which is mainly due to the difference in repeat content (47.2% in Cfu versus 3.2% in Dse). Recent adaptation to different lifestyles and hosts is suggested by diverged sets of genes. Cfu contains an α-tomatinase gene that we predict might be required for detoxification of tomatine, while this gene is absent in Dse. Many genes encoding secreted proteins are unique to each species and the repeat-rich areas in Cfu are enriched for these species-specific genes. In contrast, conserved genes suggest common host ancestry. Homologs of Cfu effector genes, including Ecp2 and Avr4, are present in Dse and induce a Cf-Ecp2- and Cf-4-mediated hypersensitive response, respectively. Strikingly, genes involved in production of the toxin dothistromin, a likely virulence factor for Dse, are conserved in Cfu, but their expression differs markedly with essentially no expression by Cfu in planta. Likewise, Cfu has a carbohydrate-degrading enzyme catalog that is more similar to that of necrotrophs or hemibiotrophs and a larger pectinolytic gene arsenal than Dse, but many of these genes are not expressed in planta or are pseudogenized. Overall, comparison of their genomes suggests that these closely related plant pathogens had a common ancestral host but since adapted to different hosts and lifestyles by a combination of differentiated gene content, pseudogenization, and gene regulation.

A Survey on Occurrence of Cladosporium fulvum Identifies Race 0 and Race 2 in Tomato-Growing Areas of Argentina.

The presence of Cladosporium fulvum (syn. Passalora fulva), causal agent of tomato leaf mold, was confirmed in the two main greenhouse-production areas for tomato in Argentina. Using both morphological characters and internal transcribed spacer sequencing, we confirmed the presence of physiological races of this pathogen. A diagnostic multiplex polymerase chain reaction (PCR) was also developed, using primers derived from C. fulvum avirulence (Avr) genes. In all, 20 isolates of Cladosporium spp. were obtained as monospore cultures and 12 were identified as C. fulvum. By this method, we showed that, of these 12 isolates, 5 were race 0 (carrying functional Avr2, Avr4, Avr4E, and Avr9 genes) and 7 were race 2 (lacking the Avr2 gene). Race identity was confirmed by testing their virulence on a set of tomato differentials carrying different Cf resistance genes. All Avr genes could be amplified in single or multiplex PCR using DNA isolated from in vitro grown monospore cultures but only three Avr could be amplified when genomic DNA was isolated from C. fulvum-infected necrotic leaf tissue.

Automated alignment-based curation of gene models in filamentous fungi.

BACKGROUND: Automated gene-calling is still an error-prone process, particularly for the highly plastic genomes of fungal species. Improvement through quality control and manual curation of gene models is a time-consuming process that requires skilled biologists and is only marginally performed. The wealth of available fungal genomes has not yet been exploited by an automated method that applies quality control of gene models in order to obtain more accurate genome annotations. RESULTS: We provide a novel method named alignment-based fungal gene prediction (ABFGP) that is particularly suitable for plastic genomes like those of fungi. It can assess gene models on a gene-by-gene basis making use of informant gene loci. Its performance was benchmarked on 6,965 gene models confirmed by full-length unigenes from ten different fungi. 79.4% of all gene models were correctly predicted by ABFGP. It improves the output of ab initio gene prediction software due to a higher sensitivity and precision for all gene model components. Applicability of the method was shown by revisiting the annotations of six different fungi, using gene loci from up to 29 fungal genomes as informants. Between 7,231 and 8,337 genes were assessed by ABFGP and for each genome between 1,724 and 3,505 gene model revisions were proposed. The reliability of the proposed gene models is assessed by an a posteriori introspection procedure of each intron and exon in the multiple gene model alignment. The total number and type of proposed gene model revisions in the six fungal genomes is correlated to the quality of the genome assembly, and to sequencing strategies used in the sequencing centre, highlighting different types of errors in different annotation pipelines. The ABFGP method is particularly successful in discovering sequence errors and/or disruptive mutations causing truncated and erroneous gene models. CONCLUSIONS: The ABFGP method is an accurate and fully automated quality control method for fungal gene catalogues that can be easily implemented into existing annotation pipelines. With the exponential release of new genomes, the ABFGP method will help decreasing the number of gene models that require additional manual curation.

Transcriptome sequencing uncovers the Avr5 avirulence gene of the tomato leaf mold pathogen Cladosporium fulvum.

The Cf-5 gene of tomato confers resistance to strains of the fungal pathogen Cladosporium fulvum carrying the avirulence gene Avr5. Although Cf-5 has been cloned, Avr5 has remained elusive. We report the cloning of Avr5 using a combined bioinformatic and transcriptome sequencing approach. RNA-Seq was performed on the sequenced race 0 strain (0WU; carrying Avr5), as well as a race 5 strain (IPO 1979; lacking a functional Avr5 gene) during infection of susceptible tomato. Forty-four in planta-induced C. fulvum candidate effector (CfCE) genes of 0WU were identified that putatively encode a secreted, small cysteine-rich protein. An expressed transcript sequence comparison between strains revealed two polymorphic CfCE genes in IPO 1979. One of these conferred avirulence to IPO 1979 on Cf-5 tomato following complementation with the corresponding 0WU allele, confirming identification of Avr5. Complementation also led to increased fungal biomass during infection of susceptible tomato, signifying a role for Avr5 in virulence. Seven of eight race 5 strains investigated escape Cf-5-mediated resistance through deletion of the Avr5 gene. Avr5 is heavily flanked by repetitive elements, suggesting that repeat instability, in combination with Cf-5-mediated selection pressure, has led to the emergence of race 5 strains deleted for the Avr5 gene.

Relocation of genes generates non-conserved chromosomal segments in Fusarium graminearum that show distinct and co-regulated gene expression patterns.

BACKGROUND: Genome comparisons between closely related species often show non-conserved regions across chromosomes. Some of them are located in specific regions of chromosomes and some are even confined to one or more entire chromosomes. The origin and biological relevance of these non-conserved regions are still largely unknown. Here we used the genome of Fusarium graminearum to elucidate the significance of non-conserved regions. RESULTS: The genome of F. graminearum harbours thirteen non-conserved regions dispersed over all of the four chromosomes. Using RNA-Seq data from the mycelium of F. graminearum, we found weakly expressed regions on all of the four chromosomes that exactly matched with non-conserved regions. Comparison of gene expression between two different developmental stages (conidia and mycelium) showed that the expression of genes in conserved regions is stable, while gene expression in non-conserved regions is much more influenced by developmental stage. In addition, genes involved in the production of secondary metabolites and secreted proteins are enriched in non-conserved regions, suggesting that these regions could also be important for adaptations to new environments, including adaptation to new hosts. Finally, we found evidence that non-conserved regions are generated by sequestration of genes from multiple locations. Gene relocations may lead to clustering of genes with similar expression patterns or similar biological functions, which was clearly exemplified by the PKS2 gene cluster. CONCLUSIONS: Our results showed that chromosomes can be functionally divided into conserved and non-conserved regions, and both could have specific and distinct roles in genome evolution and regulation of gene expression.

Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi.

The class Dothideomycetes is one of the largest groups of fungi with a high level of ecological diversity including many plant pathogens infecting a broad range of hosts. Here, we compare genome features of 18 members of this class, including 6 necrotrophs, 9 (hemi)biotrophs and 3 saprotrophs, to analyze genome structure, evolution, and the diverse strategies of pathogenesis. The Dothideomycetes most likely evolved from a common ancestor more than 280 million years ago. The 18 genome sequences differ dramatically in size due to variation in repetitive content, but show much less variation in number of (core) genes. Gene order appears to have been rearranged mostly within chromosomal boundaries by multiple inversions, in extant genomes frequently demarcated by adjacent simple repeats. Several Dothideomycetes contain one or more gene-poor, transposable element (TE)-rich putatively dispensable chromosomes of unknown function. The 18 Dothideomycetes offer an extensive catalogue of genes involved in cellulose degradation, proteolysis, secondary metabolism, and cysteine-rich small secreted proteins. Ancestors of the two major orders of plant pathogens in the Dothideomycetes, the Capnodiales and Pleosporales, may have had different modes of pathogenesis, with the former having fewer of these genes than the latter. Many of these genes are enriched in proximity to transposable elements, suggesting faster evolution because of the effects of repeat induced point (RIP) mutations. A syntenic block of genes, including oxidoreductases, is conserved in most Dothideomycetes and upregulated during infection in L. maculans, suggesting a possible function in response to oxidative stress.

RNA-Seq analysis reveals new gene models and alternative splicing in the fungal pathogen Fusarium graminearum.

BACKGROUND: The genome of Fusarium graminearum has been sequenced and annotated previously, but correct gene annotation remains a challenge. In addition, posttranscriptional regulations, such as alternative splicing and RNA editing, are poorly understood in F. graminearum. Here we took advantage of RNA-Seq to improve gene annotations and to identify alternative splicing and RNA editing in F. graminearum. RESULTS: We identified and revised 655 incorrectly predicted gene models, including revisions of intron predictions, intron splice sites and prediction of novel introns. 231 genes were identified with two or more alternative splice variants, mostly due to intron retention. Interestingly, the expression ratios between different transcript isoforms appeared to be developmentally regulated. Surprisingly, no RNA editing was identified in F. graminearum. Moreover, 2459 novel transcriptionally active regions (nTARs) were identified and our analysis indicates that many of these could be missed genes. Finally, we identified the 5' UTR and/or 3' UTR sequences of 7666 genes. A number of representative novel gene models and alternatively spliced genes were validated by reverse transcription polymerase chain reaction and sequencing of the generated amplicons. CONCLUSIONS: We have developed novel and efficient strategies to identify alternatively spliced genes and incorrect gene models based on RNA-Seq data. Our study identified hundreds of alternatively spliced genes in F. graminearum and for the first time indicated that alternative splicing is developmentally regulated in filamentous fungi. In addition, hundreds of incorrect predicted gene models were identified and revised and thousands of nTARs were discovered in our study, which will be helpful for the future genomic and transcriptomic studies in F. graminearum.

In silico characterization and molecular evolutionary analysis of a novel superfamily of fungal effector proteins.

Most fungal plant pathogens secrete effector proteins during pathogenesis to manipulate their host's defense and promote disease. These are so highly diverse in sequence and distribution, they are essentially considered as species-specific. However, we have recently shown the presence of homologous effectors in fungal species of the Dothideomycetes class. One such example is Ecp2, an effector originally described in the tomato pathogen Cladosporium fulvum but later detected in the plant pathogenic fungi Mycosphaerella fijiensis and Mycosphaerella graminicola as well. Here, using in silico sequence-similarity searches against a database of 135 fungal genomes and GenBank, we extend our queries for homologs of Ecp2 to the fungal kingdom and beyond, and further study their history of diversification. Our analyses show that Ecp2 homologs are members of an ancient and widely distributed superfamily of putative fungal effectors, which we term Hce2 for Homologs of C. fulvum Ecp2. Molecular evolutionary analyses show that the superfamily originated and diversified within the fungal kingdom, experiencing multiple lineage-specific expansions and losses that are consistent with the birth-and-death model of gene family evolution. Newly formed paralogs appear to be subject to diversification early after gene duplication events, whereas at later stages purifying selection acts to preserve diversity and the newly evolved putative functions. Some members of the Hce2 superfamily are fused to fungal Glycoside Hydrolase family 18 chitinases that show high similarity to the Zymocin killer toxin from the dairy yeast Kluyveromyces lactis, suggesting an analogous role in antagonistic interactions. The observed high rates of gene duplication and loss in the Hce2 superfamily, combined with diversification in both sequence and possibly functions within and between species, suggest that Hce2s are involved in adaptation to stresses and new ecological niches. Such findings address the need to rationalize effector biology and evolution beyond the perspective of solely host-microbe interactions.