Accessories Make the Outfit: Accessory Chromosomes and Other Dispensable DNA Regions in Plant-Pathogenic Fungi.

Fungal pathogen genomes can often be divided into core and accessory regions. Accessory regions ARs) may be comprised of either ARs (within core chromosomes (CCs) or wholly dispensable (accessory) chromosomes (ACs). Fungal ACs and ARs typically accumulate mutations and structural rearrangements more rapidly over time than CCs and many harbor genes relevant to host-pathogen interactions. These regions are of particular interest in plant pathology and include host-specific virulence factors and secondary metabolite synthesis gene clusters. This review outlines known ACs and ARs in fungal genomes, methods used for their detection, their common properties that differentiate them from the core genome, and what is currently known of their various roles in pathogenicity. Reports on the evolutionary processes generating and shaping AC and AR compartments are discussed, including repeat induced point mutation and breakage fusion bridge cycles. Previously ACs have been studied extensively within key genera, including Fusarium, Zymoseptoria, and Alternaria, but are growing in frequency of observation and perceived importance across a wider range of fungal species. Recent advances in sequencing technologies permit affordable genome assembly and resequencing of populations that will facilitate further discovery and routine screening of ACs.

Genome sequencing and comparative genomics of the broad host-range pathogen Rhizoctonia solani AG8.

Rhizoctonia solani is a soil-borne basidiomycete fungus with a necrotrophic lifestyle which is classified into fourteen reproductively incompatible anastomosis groups (AGs). One of these, AG8, is a devastating pathogen causing bare patch of cereals, brassicas and legumes. R. solani is a multinucleate heterokaryon containing significant heterozygosity within a single cell. This complexity posed significant challenges for the assembly of its genome. We present a high quality genome assembly of R. solani AG8 and a manually curated set of 13,964 genes supported by RNA-seq. The AG8 genome assembly used novel methods to produce a haploid representation of its heterokaryotic state. The whole-genomes of AG8, the rice pathogen AG1-IA and the potato pathogen AG3 were observed to be syntenic and co-linear. Genes and functions putatively relevant to pathogenicity were highlighted by comparing AG8 to known pathogenicity genes, orthology databases spanning 197 phytopathogenic taxa and AG1-IA. We also observed SNP-level "hypermutation" of CpG dinucleotides to TpG between AG8 nuclei, with similarities to repeat-induced point mutation (RIP). Interestingly, gene-coding regions were widely affected along with repetitive DNA, which has not been previously observed for RIP in mononuclear fungi of the Pezizomycotina. The rate of heterozygous SNP mutations within this single isolate of AG8 was observed to be higher than SNP mutation rates observed across populations of most fungal species compared. Comparative analyses were combined to predict biological processes relevant to AG8 and 308 proteins with effector-like characteristics, forming a valuable resource for further study of this pathosystem. Predicted effector-like proteins had elevated levels of non-synonymous point mutations relative to synonymous mutations (dN/dS), suggesting that they may be under diversifying selection pressures. In addition, the distant relationship to sequenced necrotrophs of the Ascomycota suggests the R. solani genome sequence may prove to be a useful resource in future comparative analysis of plant pathogens.

Transcriptome analysis of the fungal pathogen Fusarium oxysporum f. sp. medicaginis during colonisation of resistant and susceptible Medicago truncatula hosts identifies differential pathogenicity profiles and novel candidate effectors.

BACKGROUND: Pathogenic members of the Fusarium oxysporum species complex are responsible for vascular wilt disease on many important crops including legumes, where they can be one of the most destructive disease causing necrotrophic fungi. We previously developed a model legume-infecting pathosystem based on the reference legume Medicago truncatula and a pathogenic F. oxysporum forma specialis (f. sp.) medicaginis (Fom). To dissect the molecular pathogenicity arsenal used by this root-infecting pathogen, we sequenced its transcriptome during infection of a susceptible and resistant host accession. RESULTS: High coverage RNA-Seq of Fom infected root samples harvested from susceptible (DZA315) or resistant (A17) M. truncatula seedlings at early or later stages of infection (2 or 7 days post infection (dpi)) and from vegetative (in vitro) samples facilitated the identification of unique and overlapping sets of in planta differentially expressed genes. This included enrichment, particularly in DZA315 in planta up-regulated datasets, for proteins associated with sugar, protein and plant cell wall metabolism, membrane transport, nutrient uptake and oxidative processes. Genes encoding effector-like proteins were identified, including homologues of the F. oxysporum f. sp. lycopersici Secreted In Xylem (SIX) proteins, and several novel candidate effectors based on predicted secretion, small protein size and high in-planta induced expression. The majority of the effector candidates contain no known protein domains but do share high similarity to predicted proteins predominantly from other F. oxysporum ff. spp. as well as other Fusaria (F. solani, F. fujikori, F. verticilloides, F. graminearum and F. pseudograminearum), and from another wilt pathogen of the same class, a Verticillium species. Overall, this suggests these novel effector candidates may play important roles in Fusaria and wilt pathogen virulence. CONCLUSION: Combining high coverage in planta RNA-Seq with knowledge of fungal pathogenicity protein features facilitated the identification of differentially expressed pathogenicity associated genes and novel effector candidates expressed during infection of a resistant or susceptible M. truncatula host. The knowledge from this first in depth in planta transcriptome sequencing of any F. oxysporum ff. spp. pathogenic on legumes will facilitate the dissection of Fusarium wilt pathogenicity mechanisms on many important legume crops.

Comparative genomics and prediction of conditionally dispensable sequences in legume-infecting Fusarium oxysporum formae speciales facilitates identification of candidate effectors.

BACKGROUND: Soil-borne fungi of the Fusarium oxysporum species complex cause devastating wilt disease on many crops including legumes that supply human dietary protein needs across many parts of the globe. We present and compare draft genome assemblies for three legume-infecting formae speciales (ff. spp.): F. oxysporum f. sp. ciceris (Foc-38-1) and f. sp. pisi (Fop-37622), significant pathogens of chickpea and pea respectively, the world's second and third most important grain legumes, and lastly f. sp. medicaginis (Fom-5190a) for which we developed a model legume pathosystem utilising Medicago truncatula. RESULTS: Focusing on the identification of pathogenicity gene content, we leveraged the reference genomes of Fusarium pathogens F. oxysporum f. sp. lycopersici (tomato-infecting) and F. solani (pea-infecting) and their well-characterised core and dispensable chromosomes to predict genomic organisation in the newly sequenced legume-infecting isolates. Dispensable chromosomes are not essential for growth and in Fusarium species are known to be enriched in host-specificity and pathogenicity-associated genes. Comparative genomics of the publicly available Fusarium species revealed differential patterns of sequence conservation across F. oxysporum formae speciales, with legume-pathogenic formae speciales not exhibiting greater sequence conservation between them relative to non-legume-infecting formae speciales, possibly indicating the lack of a common ancestral source for legume pathogenicity. Combining predicted dispensable gene content with in planta expression in the model legume-infecting isolate, we identified small conserved regions and candidate effectors, four of which shared greatest similarity to proteins from another legume-infecting ff. spp. CONCLUSIONS: We demonstrate that distinction of core and potential dispensable genomic regions of novel F. oxysporum genomes is an effective tool to facilitate effector discovery and the identification of gene content possibly linked to host specificity. While the legume-infecting isolates didn't share large genomic regions of pathogenicity-related content, smaller regions and candidate effector proteins were highly conserved, suggesting that they may play specific roles in inducing disease on legume hosts.

Next-generation re-sequencing as a tool for rapid bioinformatic screening of presence and absence of genes and accessory chromosomes across isolates of Zymoseptoria tritici.

The wheat pathogen Zymoseptoria tritici possesses a large number of accessory chromosomes that may be present or absent in its genome. The genome of the reference isolate IPO323 has been assembled to a very high standard and contains 21 full length chromosome sequences, 8 of which represent accessory chromosomes. The IPO323 reference, when combined with low-cost next-generation sequencing and bioinformatics, can be used as a powerful tool to assess the presence or absence of accessory chromosomes. We present an outline of a range of bioinformatics techniques that can be applied to the analysis of presence-absence variation among accessory chromosomes across 13 novel isolates of Z. tritici.

The genome sequence of the biocontrol fungus Metarhizium anisopliae and comparative genomics of Metarhizium species.

BACKGROUND: Metarhizium anisopliae is an important fungal biocontrol agent of insect pests of agricultural crops. Genomics can aid the successful commercialization of biopesticides by identification of key genes differentiating closely related species, selection of virulent microbial isolates which are amenable to industrial scale production and formulation and through the reduction of phenotypic variability. The genome of Metarhizium isolate ARSEF23 was recently published as a model for M. anisopliae, however phylogenetic analysis has since re-classified this isolate as M. robertsii. We present a new annotated genome sequence of M. anisopliae (isolate Ma69) and whole genome comparison to M. robertsii (ARSEF23) and M. acridum (CQMa 102). RESULTS: Whole genome analysis of M. anisopliae indicates significant macrosynteny with M. robertsii but with some large genomic inversions. In comparison to M. acridum, the genome of M. anisopliae shares lower sequence homology. While alignments overall are co-linear, the genome of M. acridum is not contiguous enough to conclusively observe macrosynteny. Mating type gene analysis revealed both MAT1-1 and MAT1-2 genes present in M. anisopliae suggesting putative homothallism, despite having no known teleomorph, in contrast with the putatively heterothallic M. acridum isolate CQMa 102 (MAT1-2) and M. robertsii isolate ARSEF23 (altered MAT1-1). Repetitive DNA and RIP analysis revealed M. acridum to have twice the repetitive content of the other two species and M. anisopliae to be five times more RIP affected than M. robertsii. We also present an initial bioinformatic survey of candidate pathogenicity genes in M. anisopliae. CONCLUSIONS: The annotated genome of M. anisopliae is an important resource for the identification of virulence genes specific to M. anisopliae and development of species- and strain- specific assays. New insight into the possibility of homothallism and RIP affectedness has important implications for the development of M. anisopliae as a biopesticide as it may indicate the potential for greater inherent diversity in this species than the other species. This could present opportunities to select isolates with unique combinations of pathogenicity factors, or it may point to instability in the species, a negative attribute in a biopesticide.

Reference genome assembly for Australian Ascochyta lentis isolate Al4.

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use, and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

Reference Genome Assembly for Australian Ascochyta rabiei Isolate ArME14.

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

Identification of a Polyketide Synthase Gene Responsible for Ascochitine Biosynthesis in Ascochyta fabae and Its Abrogation in Sister Taxa.

The polyketide-derived secondary metabolite ascochitine is produced by species in the Didymellaceae family, including but not restricted to Ascochyta species pathogens of cool-season food legumes. Ascochitine is structurally similar to the well-known mycotoxin citrinin and exhibits broad-spectrum phytotoxicity and antimicrobial activities. Here, we identified a polyketide synthase (PKS) gene (denoted pksAC) responsible for ascochitine production in the filamentous fungus Ascochyta fabae Deletion of the pksAC prevented production of ascochitine and its derivative ascochital in A. fabae The putative ascochitine biosynthesis gene cluster comprises 11 genes that have undergone rearrangement and gain-and-loss events relative to the citrinin biosynthesis gene cluster in Monascus ruber Interestingly, we also identified pksAC homologs in two recently diverged species, A. lentis and A. lentis var. lathyri, that are sister taxa closely related to ascochitine producers such as A. fabae and A. viciae-villosae However, nonsense mutations have been independently introduced in coding sequences of the pksAC homologs of A. lentis and A. lentis var. lathyri that resulted in loss of ascochitine production. Despite its reported phytotoxicity, ascochitine was not a pathogenicity factor in A. fabae infection and colonization of faba bean (Vicia faba L.). Ascochitine was mainly produced from mature hyphae at the site of pycnidial formation, suggesting a possible protective role of the compound against other microbial competitors in nature. This report highlights the evolution of gene clusters harnessing the structural diversity of polyketides and a mechanism with the potential to alter secondary metabolite profiles via single nucleotide polymorphisms in closely related fungal species.IMPORTANCE Fungi produce a diverse array of secondary metabolites, many of which are of pharmacological importance whereas many others are noted for mycotoxins, such as aflatoxin and citrinin, that can threaten human and animal health. The polyketide-derived compound ascochitine, which is structurally similar to citrinin mycotoxin, has been considered to be important for pathogenicity of legume-associated Ascochyta species. Here, we identified the ascochitine polyketide synthase (PKS) gene in Ascochyta fabae and its neighboring genes that may be involved in ascochitine biosynthesis. Interestingly, the ascochitine PKS genes in other legume-associated Ascochyta species have been mutated, encoding truncated PKSs. This indicated that point mutations may have contributed to genetic diversity for secondary metabolite production in these fungi. We also demonstrated that ascochitine is not a pathogenicity factor in A. fabae The antifungal activities and production of ascochitine during sporulation suggested that it may play a role in competition with other saprobic fungi in nature.

Utilizing Gene Tree Variation to Identify Candidate Effector Genes in Zymoseptoria tritici.

Zymoseptoria tritici is a host-specific, necrotrophic pathogen of wheat. Infection by Z. tritici is characterized by its extended latent period, which typically lasts 2 wks, and is followed by extensive host cell death, and rapid proliferation of fungal biomass. This work characterizes the level of genomic variation in 13 isolates, for which we have measured virulence on 11 wheat cultivars with differential resistance genes. Between the reference isolate, IPO323, and the 13 Australian isolates we identified over 800,000 single nucleotide polymorphisms, of which ∼10% had an effect on the coding regions of the genome. Furthermore, we identified over 1700 probable presence/absence polymorphisms in genes across the Australian isolates using de novo assembly. Finally, we developed a gene tree sorting method that quickly identifies groups of isolates within a single gene alignment whose sequence haplotypes correspond with virulence scores on a single wheat cultivar. Using this method, we have identified < 100 candidate effector genes whose gene sequence correlates with virulence toward a wheat cultivar carrying a major resistance gene.

Comprehensive Annotation of the Parastagonospora nodorum Reference Genome Using Next-Generation Genomics, Transcriptomics and Proteogenomics.

Parastagonospora nodorum, the causal agent of Septoria nodorum blotch (SNB), is an economically important pathogen of wheat (Triticum spp.), and a model for the study of necrotrophic pathology and genome evolution. The reference P. nodorum strain SN15 was the first Dothideomycete with a published genome sequence, and has been used as the basis for comparison within and between species. Here we present an updated reference genome assembly with corrections of SNP and indel errors in the underlying genome assembly from deep resequencing data as well as extensive manual annotation of gene models using transcriptomic and proteomic sources of evidence (https://github.com/robsyme/Parastagonospora_nodorum_SN15). The updated assembly and annotation includes 8,366 genes with modified protein sequence and 866 new genes. This study shows the benefits of using a wide variety of experimental methods allied to expert curation to generate a reliable set of gene models.

Evaluation of Secretion Prediction Highlights Differing Approaches Needed for Oomycete and Fungal Effectors.

The steadily increasing number of sequenced fungal and oomycete genomes has enabled detailed studies of how these eukaryotic microbes infect plants and cause devastating losses in food crops. During infection, fungal and oomycete pathogens secrete effector molecules which manipulate host plant cell processes to the pathogen's advantage. Proteinaceous effectors are synthesized intracellularly and must be externalized to interact with host cells. Computational prediction of secreted proteins from genomic sequences is an important technique to narrow down the candidate effector repertoire for subsequent experimental validation. In this study, we benchmark secretion prediction tools on experimentally validated fungal and oomycete effectors. We observe that for a set of fungal SwissProt protein sequences, SignalP 4 and the neural network predictors of SignalP 3 (D-score) and SignalP 2 perform best. For effector prediction in particular, the use of a sensitive method can be desirable to obtain the most complete candidate effector set. We show that the neural network predictors of SignalP 2 and 3, as well as TargetP were the most sensitive tools for fungal effector secretion prediction, whereas the hidden Markov model predictors of SignalP 2 and 3 were the most sensitive tools for oomycete effectors. Thus, previous versions of SignalP retain value for oomycete effector prediction, as the current version, SignalP 4, was unable to reliably predict the signal peptide of the oomycete Crinkler effectors in the test set. Our assessment of subcellular localization predictors shows that cytoplasmic effectors are often predicted as not extracellular. This limits the reliability of secretion predictions that depend on these tools. We present our assessment with a view to informing future pathogenomics studies and suggest revised pipelines for secretion prediction to obtain optimal effector predictions in fungi and oomycetes.

Phoma medicaginis stimulates the induction of the octadecanoid and phenylpropanoid pathways in Medicago truncatula.

Gene expression changes and metabolite abundances were measured during the interaction of Medicago truncatula with the fungal necrotrophic pathogen Phoma medicaginis in leaf tissue of susceptible and resistant accessions. Over 330 genes were differentially expressed in plants infected with P. medicaginis relative to mock-inoculated plants at 12 h post-inoculation. Of these, 191 were induced in either the resistant or the susceptible accession, with 143 genes repressed. Expression changes were observed in genes involved in the oxidative burst, cell wall strengthening and lipid metabolism, as well as several transcription factors. Genes related to salicylic acid, jasmonate and ethylene responses were up-regulated, as well as genes leading to the production of jasmonic acid. Significant induction of genes in the phenylpropanoid pathway leading to lignin and isoflavonoid biosynthesis occurred. High-pressure liquid chromatography with UV detection (HPLC-UV) identified several phenolic compounds induced by P. medicaginis, as well as constitutively higher levels of phenolic compounds, in the resistant M. truncatula accession. Differentially regulated genes induced in both the resistant and susceptible accessions, but with different kinetics, and constitutively more highly expressed and induced phenolic compounds provide candidates for functional analysis. Taken together, these results highlight the importance of the octadecanoid and phenylpropanoid pathways in defence against this necrotrophic pathogen.

The Medicago truncatula reference accession A17 has an aberrant chromosomal configuration.

Medicago truncatula (barrel medic) has emerged as a model legume and accession A17 is the reference genotype selected for the sequencing of the genome. In the present study we compare the A17 chromosomal configuration with that of other accessions by examining pollen viability and genetic maps of intraspecific hybrids. Hybrids derived from crosses between M. truncatula accessions, representative of the large genetic variation within the germplasm collection, were evaluated for pollen viability using Alexander's stain. Genetic maps were generated for the following crosses: SA27063 x SA3054 (n = 94), SA27063 x A17 (n = 92), A17 x Borung (n = 99) and A17 x A20 (n = 69). All F(1) individuals derived from crosses involving A17 showed 50% pollen viability or less. Examination of the recombination frequencies between markers of chromosomes 4 and 8 revealed an apparent genetic linkage between the lower arms of these chromosomes in genetic maps derived from A17. Semisterility and unexpected linkage relationship are both good indicators of a reciprocal translocation. The implications of the A17 distinctive chromosomal rearrangement on studies of genetic mapping, genome sequencing and synteny between species are discussed.