Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis

Septoria Leaf Blotch and Reduced Nitrogen Availability Alter WRKY Transcription Factor Expression in a Codependent Manner

A major cause of yield loss in wheat worldwide is the fungal pathogen Zymoseptoria tritici, a hemibiotrophic fungus which causes Septoria leaf blotch, the most destructive wheat disease in Europe. Resistance in commercial wheat varieties is poor, however, a link between reduced nitrogen availability and increased Septoria tolerance has been observed. We have shown that Septoria load is not affected by nitrogen, whilst the fungus is in its first, symptomless stage of growth. This suggests that a link between nitrogen and Septoria is only present during the necrotrophic phase of Septoria infection. Quantitative real-time PCR data demonstrated that WRKYs, a superfamily of plant-specific transcription factors, are differentially expressed in response to both reduced nitrogen and Septoria. WRKY39 was downregulated over 30-fold in response to necrotrophic stage Septoria, whilst changes in the expression of WRKY68a during the late biotrophic phase were dependent on the concentration of nitrogen under which wheat is grown. WRKY68a may therefore mediate a link between nitrogen and Septoria. The potential remains to identify key regulators in the link between nitrogen and Septoria, and as such, elucidate molecular markers for wheat breeding, or targets for molecular-based breeding approaches.

Polycomb Repression without Bristles: Facultative Heterochromatin and Genome Stability in Fungi

Genome integrity is essential to maintain cellular function and viability. Consequently, genome instability is frequently associated with dysfunction in cells and associated with plant, animal, and human diseases. One consequence of relaxed genome maintenance that may be less appreciated is an increased potential for rapid adaptation to changing environments in all organisms. Here, we discuss evidence for the control and function of facultative heterochromatin, which is delineated by methylation of histone H3 lysine 27 (H3K27me) in many fungi. Aside from its relatively well understood role in transcriptional repression, accumulating evidence suggests that H3K27 methylation has an important role in controlling the balance between maintenance and generation of novelty in fungal genomes. We present a working model for a minimal repressive network mediated by H3K27 methylation in fungi and outline challenges for future research.

Wheat pathogen Zymoseptoria tritici N-myristoyltransferase inhibitors: on-target antifungal activity and an unusual metabolic defense mechanism

Zymoseptoria tritici is the causative agent of Septoria tritici blotch (STB), which costs billions of dollars annually to major wheat-producing countries in terms of both fungicide use and crop loss. Agricultural pathogenic fungi have acquired resistance to most commercially available fungicide classes, and the rate of discovery and development of new fungicides has stalled, demanding new approaches and insights. Here we investigate a potential mechanism of targeting an important wheat pathogen Z. tritici via inhibition of N-myristoyltransferase (NMT). We characterize Z. tritici NMT biochemically for the first time, profile the in vivo Z. tritici myristoylated proteome and identify and validate the first Z. tritici NMT inhibitors. Proteomic investigation of the downstream effects of NMT inhibition identified an unusual and novel mechanism of defense against chemical toxicity in Z. tritici through the application of comparative bioinformatics to deconvolute function from the previously largely unannotated Z. tritici proteome. Research into novel fungicidal modes-of-action is essential to satisfy an urgent unmet need for novel fungicide targets, and we anticipate that this study will serve as a useful proteomics and bioinformatics resource for researchers studying Z. tritici.

Genetic diversity and population structure of Zymoseptoria tritici in Ethiopia as revealed by microsatellite markers

Septoria tritici blotch (STB), caused by Zymoseptoria tritici (formerly: Mycosphaerella graminicola or Septoria tritici), is one of the most devastating diseases of wheat globally. Understanding genetic diversity of the pathogen has supreme importance in developing best management strategies. However, there is dearth of information on the genetic structure of Z. tritici populations in Ethiopia. Therefore, the present study was targeted to uncover the genetic diversity and population structure of Z. tritici populations from the major wheat-growing areas of Ethiopia. Totally, 182 Z. tritici isolates representing eight populations were analyzed with 14 microsatellite markers. All the microsatellite loci were polymorphic and highly informative, and hence useful genetic tools to depict the genetic diversity and population structure of the pathogen. A wide range of diversity indices including number of observed alleles, effective number of alleles, Shannon's diversity index, number of private alleles, Nei's gene diversity and percentage of polymorphic loci (PPL) were computed to determine genetic variation within populations. A high within-populations genetic diversity was confirmed with gene diversity index and PPL values ranging from 0.34 - 0.58 and 79-100% with overall mean of 0.45 and 94%, respectively. Analysis of molecular variance (AMOVA) revealed a moderate genetic differentiation where 92% of the total genetic variation resides within populations, leaving only 8% among populations. Cluster (UPGMA), PCoA and STRUCTURE analyses did not group the populations into sharply genetically distinct clusters according to their geographical origins, likely due to high gene flow (Nm = 5.66) and reproductive biology of the pathogen. All individual samples shared alleles from two subgroups (K = 2) evidencing high potential of genetic admixture. In conclusion, the microsatellite markers used in the present study were highly informative and thus, helped to dissect the genetic structures of Z. tritici populations in Ethiopia. Among the studied populations, those of East Shewa, Arsi, South West Shewa and Bale showed a high genetic diversity, and hence these areas can be considered as hot spots for investigations planned on the pathogen and host-pathogen interactions. Therefore, the present study not only enriches missing information in Ethiopia but also provides new insights into the epidemiology and genetic structure of Z. tritici in Africa where the agro-climatic conditions and the wheat cropping systems are different from other parts of the world. Such baseline information is useful for designing and implementing durable and effective management strategies.

The Many Questions about Mini Chromosomes in Colletotrichum spp

Many fungal pathogens carry accessory regions in their genome, which are not required for vegetative fitness. Often, although not always, these regions occur as relatively small chromosomes in different species. Such mini chromosomes appear to be a typical feature of many filamentous plant pathogens. Since these regions often carry genes coding for effectors or toxin-producing enzymes, they may be directly related to virulence of the respective pathogen. In this review, we outline the situation of small accessory chromosomes in the genus Colletotrichum, which accounts for ecologically important plant diseases. We summarize which species carry accessory chromosomes, their gene content, and chromosomal makeup. We discuss the large variation in size and number even between different isolates of the same species, their potential roles in host range, and possible mechanisms for intra- and interspecies exchange of these interesting genetic elements.

EffHunter: A Tool for Prediction of Effector Protein Candidates in Fungal Proteomic Databases

Pathogens are able to deliver small-secreted, cysteine-rich proteins into plant cells to enable infection. The computational prediction of effector proteins remains one of the most challenging areas in the study of plant fungi interactions. At present, there are several bioinformatic programs that can help in the identification of these proteins; however, in most cases, these programs are managed independently. Here, we present EffHunter, an easy and fast bioinformatics tool for the identification of effectors. This predictor was used to identify putative effectors in 88 proteomes using characteristics such as size, cysteine residue content, secretion signal and transmembrane domains.

A chromosome-scale assembly of the smallest Dothideomycete genome reveals a unique genome compaction mechanism in filamentous fungi

Background

The wide variation in the size of fungal genomes is well known, but the reasons for this size variation are less certain. Here, we present a chromosome-scale assembly of ectophytic Peltaster fructicola, a surface-dwelling extremophile, based on long-read DNA sequencing technology, to assess possible mechanisms associated with genome compaction.

Results

At 18.99 million bases (Mb), P. fructicola possesses one of the smallest known genomes sequence among filamentous fungi. The genome is highly compact relative to other fungi, with substantial reductions in repeat content, ribosomal DNA copies, tRNA gene quantity, and intron sizes, as well as intergenic lengths and the size of gene families. Transposons take up just 0.05% of the entire genome, and no full-length transposon was found. We concluded that reduced genome sizes in filamentous fungi such as P. fructicola, Taphrina deformans and Pneumocystis jirovecii occurred through reduction in ribosomal DNA copy number and reduced intron sizes. These dual mechanisms contrast with genome reduction in the yeast fungus Saccharomyces cerevisiae, whose small and compact genome is associated solely with intron loss.

Conclusions

Our results reveal a unique genomic compaction architecture of filamentous fungi inhabiting plant surfaces, and broaden the understanding of the mechanisms associated with compaction of fungal genomes.

Prediction of the secretomes of endophytic and nonendophytic fungi reveals similarities in host plant infection and colonization strategies

Endophytic fungi are microorganisms that inhabit internal plant tissues without causing apparent damage. During the infection process, both endophytic and phytopathogenic fungi secrete proteins to resist or supplant the plant's defense mechanisms. This study analyzed the predicted secretomes of six species of endophytic fungi and compared them with predicted secretomes of eight fungal species with different lifestyles: saprophytic, necrotrophic, hemibiotrophic, and biotrophic. The sizes of the predicted secretomes varied from 260 to 1640 proteins, and the predicted secretomes have a wide diversity of CAZymes, proteases, and conserved domains. Regarding the CAZymes in the secretomes of the analyzed fungi, the most abundant CAZyme families were glycosyl hydrolase and serine proteases. Several predicted proteins have characteristics similar to those found in small, secreted proteins with effector characteristics (SSPEC). The most abundant conserved domains, besides those found in the SSPEC, have oxidation activities, indicating that these proteins can protect the fungus against oxidative stress, against domains with protease activity, which may be involved in the mechanisms of nutrition, or against lytic enzymes secreted by the host plant. This study demonstrates that secretomes of endophytic and nonendophytic fungi share an arsenal of proteins important in the process of infection and colonization of host plants.

Genomic Instability in Fungal Plant Pathogens

Fungi and fungal-like organisms (oomycetes) that cause diseases in plants have impacted human communities for centuries and probably from the dawn of agriculture. In modern agriculture, there is a constant race between new strategies to manage fungal plant pathogens and their ability to adapt. An important component in this race is fungal genetic diversity. Mechanisms such as sexual and parasexual recombination that contribute to the creation of novel allele combinations in fungal plant pathogens are briefly discussed in the first part of this review. Advances in genomics have enabled the investigation of chromosomal aberrations of agriculturally important fungal isolates at the nucleotide level. Some of these cases are summarized in the second part of this review; it is claimed that the effect of chromosomal aberrations on pathogenicity should be studied mechanistically. More data on the effect of gene copy number variations on phenotypes that are relevant to agriculture are especially needed. Genome rearrangements through translocations have shaped the genome of fungal plant pathogens by creating lineage-specific chromosome territories encoding for genes participating in plant diseases. Pathogenicity chromosomes are unique cases of such lineage-specific genetic elements, interestingly these chromosomes can be transferred horizontally and thus transforming a non-pathogenic strain to a pathogenic one. The third part of this review describes our attempts to reveal mutators in fungal plant pathogens by identifying fungi that lack important DNA repair genes or respond to DNA damage in an unconventional way. We found that a group of fungal plant pathogens lack conserved genes that are needed for an important Holliday junction resolution pathway. In addition, in Fusarium oxysporum, the rate-limiting step in dNTP production is not induced under DNA replication stress. This is very different from organisms from bacteria to humans. It remains to be seen if these mechanisms promote genetic instability in fungal plant pathogens.

Interspecific Gene Exchange Introduces High Genetic Variability in Crop Pathogen

Genome analyses have revealed a profound role of hybridization and introgression in the evolution of many eukaryote lineages, including fungi. The impact of recurrent introgression on fungal evolution however remains elusive. Here, we analyzed signatures of introgression along the genome of the fungal wheat pathogen Zymoseptoria tritici. We applied a comparative population genomics approach, including genome data from five Zymoseptoria species, to characterize the distribution and composition of introgressed regions representing segments with an exceptional haplotype pattern. These regions are found throughout the genome, comprising 5% of the total genome and overlapping with > 1,000 predicted genes. We performed window-based phylogenetic analyses along the genome to distinguish regions which have a monophyletic or nonmonophyletic origin with Z. tritici sequences. A majority of nonmonophyletic windows overlap with the highly variable regions suggesting that these originate from introgression. We verified that incongruent gene genealogies do not result from incomplete lineage sorting by comparing the observed and expected length distribution of haplotype blocks resulting from incomplete lineage sorting. Although protein-coding genes are not enriched in these regions, we identify 18 that encode putative virulence determinants. Moreover, we find an enrichment of transposable elements in these regions implying that hybridization may contribute to the horizontal spread of transposable elements. We detected a similar pattern in the closely related species Zymoseptoria ardabiliae, suggesting that hybridization is widespread among these closely related grass pathogens. Overall, our results demonstrate a significant impact of recurrent hybridization on overall genome evolution of this important wheat pathogen.

Karyotypic instability of endoprophase and mitotic cells of Amoeba sp. strain Cont from the "proteus-type" group (Amoebozoa, Euamoebida, Amoebidae)

We performed karyotyping of Amoeba sp. strain Cont. Based on the results of a cytological analysis, we concluded that the chromosome number of Amoeba sp. strain Cont in mitosis was unstable. In all cases they appeared to be hypergaploid (the basic chromosome number is 30), with monosomy of all chromosomes except four shortest ones. The presence of "extrachromosomes" in the nucleus could prolong until the beginning of the anaphase. It was only then that they were ejected from the nucleus and the euploidy (haploidy) was restored. The stage of endoprophase nucleus was revealed in the cell cycle of Amoeba sp. strain Cont. This stage has not yet been found in other amoebae from the "proteus-type" group that had been previously studied (A. proteus strain B and A. borokensis). The maximum number of endoreplication rounds in the strain Cont amoebae nuclear cycle was 4 or 5. The regular extrusion of chromosomes from the nucleus into the cytoplasm occurred in each of the endoreplication rounds. Comparative cytological analysis of A. proteus strain B, A. borokensis and Amoeba sp. strain Cont karyotypes indicated that strain Cont, though rather close to the former two amoebae, is actually a distinct species.

Chromosomal Differentiation in Genetically Isolated Populations of the Marsh-Specialist Crocidura suaveolens (Mammalia: Soricidae)

The genus Crocidura represents a remarkable model for the study of chromosome evolution. This is the case of the lesser white-toothed shrew (Crocidura suaveolens), a representative of the Palearctic group. Although continuously distributed from Siberia to Central Europe, C. suaveolens is a rare, habitat-specialist species in the southwesternmost limit of its distributional range, in the Gulf of Cádiz (Iberian Peninsula). In this area, C. suaveolens is restricted to genetically isolated populations associated to the tidal marches of five rivers (Guadiana, Piedras, Odiel, Tinto and Guadalquivir). This particular distributional range provides a unique opportunity to investigate whether genetic differentiation and habitat specialization was accompanied by chromosomal variation. In this context, the main objective of this study was to determinate the chromosomal characteristics of the habitat-specialist C. suaveolens in Southwestern Iberia, as a way to understand the evolutionary history of this species in the Iberian Peninsula. A total of 41 individuals from six different populations across the Gulf of Cádiz were collected and cytogenetically characterized. We detected four different karyotypes, with diploid numbers (2n) ranging from 2n = 40 to 2n = 43. Two of them (2n = 41 and 2n = 43) were characterized by the presence of B-chromosomes. The analysis of karyotype distribution across lineages and populations revealed an association between mtDNA population divergence and chromosomal differentiation. C. suaveolens populations in the Gulf of Cádiz provide a rare example of true karyotypic polymorphism potentially associated to genetic isolation and habitat specialization in which to investigate the evolutionary significance of chromosomal variation in mammals and their contribution to phenotypic and ecological divergence.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.

RESULTS

Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.

CONCLUSIONS

Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.

RESULTS

Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.

CONCLUSIONS

Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

Stress-Driven Transposable Element De-repression Dynamics and Virulence Evolution in a Fungal Pathogen

Transposable elements (TEs) are drivers of genome evolution and affect the expression landscape of the host genome. Stress is a major factor inducing TE activity; however, the regulatory mechanisms underlying de-repression are poorly understood. Plant pathogens are excellent models to dissect the impact of stress on TEs. The process of plant infection induces stress for the pathogen, and virulence factors (i.e., effectors) located in TE-rich regions become expressed. To dissect TE de-repression dynamics and contributions to virulence, we analyzed the TE expression landscape of four strains of the major wheat pathogen Zymoseptoria tritici. We experimentally exposed strains to nutrient starvation and host infection stress. Contrary to expectations, we show that the two distinct conditions induce the expression of different sets of TEs. In particular, the most highly expressed TEs, including miniature inverted-repeat transposable element and long terminal repeat-Gypsy element, show highly distinct de-repression across stress conditions. Both the genomic context of TEs and the genetic background stress (i.e., different strains harboring the same TEs) were major predictors of de-repression under stress. Gene expression profiles under stress varied significantly depending on the proximity to the closest TEs and genomic defenses against TEs were largely ineffective to prevent de-repression. Next, we analyzed the locus encoding the Avr3D1 effector. We show that the insertion and subsequent silencing of TEs in close proximity likely contributed to reduced expression and virulence on a specific wheat cultivar. The complexity of TE responsiveness to stress across genetic backgrounds and genomic locations demonstrates substantial intraspecific genetic variation to control TEs with consequences for virulence.

Does the Latent Period of Leaf Fungal Pathogens Reflect Their Trophic Type? A Meta-Analysis of Biotrophs, Hemibiotrophs, and Necrotrophs

We performed a meta-analysis to search for a relation between the trophic type and latent period of fungal pathogens. The pathogen incubation period and the level of resistance of the hosts were also investigated. This ecological knowledge would help us to more efficiently regulate crop epidemics for different types of pathogens. We gathered latent period data from 103 studies dealing with 51 fungal pathogens of the three major trophic types (25 biotrophs, 15 hemibiotrophs, and 11 necrotrophs), representing 2,542 mean latent periods. We show that these three trophic types display significantly different latent periods. Necrotrophs exhibited the shortest latent periods (<100 degree-days [DD]), biotrophs had intermediate ones (between 100 and 200 DD), and hemibiotrophs had the longest latent periods (>200 DD). We argue that this relation between trophic type and latent period points to two opposing host exploitation strategies: necrotrophs mount a rapid destructive attack on the host tissue, whereas biotrophs and hemibiotrophs avoid or delay the damaging phase. We query the definition of hemibiotrophic pathogens and discuss whether the length of the latent period is determined by the physiological limits inherent to each trophic type or by the adaptation of pathogens of different trophic types to the contrasting conditions experienced in their interaction with the host.

Genomic Sequencing of Phyllosticta citriasiana Provides Insight Into Its Conservation and Diversification With Two Closely Related Phyllosticta Species Associated With Citrus

Phyllosticta capitalensis, Phyllosticta citricarpa, and Phyllosticta citriasiana are three very important Phyllosticta species associated with citrus. P. capitalensis is an endophyte fungus of citrus while P. citricarpa can cause black spot of citrus (e.g., oranges and mandarins). P. citriasiana was identified recently which is the causal agent of the pomelo tan spot. Here, we present the ∼34 Mb genome of P. citriasiana. The genome is organized in 92 contigs, encompassing 9202 predicted genes. Comparative genomic analyses with two other Phyllosticta species (P. citricarpa and P. capitalensis) associated with citrus was conducted to understand their evolutionary conservation and diversification. Pair-wise genome alignments revealed that these species are highly syntenic. All species encode similar numbers of CAZymes and secreted proteins. However, the molecular functions of the secretome showed that each species contains some enzymes with distinct activities. The three Phyllosticta species investigated shared a core set of 7261 protein families. P. capitalensis had the largest set of orphan genes (1991), in complete contrast to that of P. citriasiana (364) and P. citricarpa (262). Most of the orphan genes are functionally unknown, but they contain a certain number of species-specific secreted proteins. A total of 23 secondary metabolites biosynthesis clusters were identified in the three Phyllosticta species, 21 of them being highly conserved among these species while the remaining two showed whole cluster gain and loss polymorphisms or gene content polymorphisms. Taken together, our study reveals insights into the genetic mechanisms of host adaptation of three species of Phyllosticta associated with citrus and paves the way to identify effectors that function in infection of citrus plants.

Small RNA Bidirectional Crosstalk During the Interaction Between Wheat and Zymoseptoria tritici

Cross-kingdom RNA interference (RNAi) has been shown to play important roles during plant-pathogen interactions, and both plants and pathogens can use small RNAs (sRNAs) to silence genes in each other. This bidirectional cross-kingdom RNAi was still unexplored in the wheat-Zymoseptoria tritici pathosystem. Here, we performed a detailed analysis of the sRNA bidirectional crosstalk between wheat and Z. tritici. Using a combination of small RNA sequencing (sRNA-seq) and microRNA sequencing (mRNA-seq), we were able to identify known and novel sRNAs and study their expression and their action on putative targets in both wheat and Z. tritici. We predicted the target genes of all the sRNAs in either wheat or Z. tritici transcriptome and used degradome analysis to validate the cleavage of these gene transcripts. We could not find any clear evidence of a cross-kingdom RNAi acting by mRNA cleavage in this pathosystem. We also found that the fungal sRNA enrichment was lower in planta than during in vitro growth, probably due to the lower expression of the only Dicer gene of the fungus during plant infection. Our results support the recent finding that Z. tritici sRNAs cannot play important roles during wheat infection. However, we also found that the fungal infection induced wheat sRNAs regulating the expression of specific wheat genes, including auxin-related genes, as an immune response. These results indicate a role of sRNAs in the regulation of wheat defenses during Z. tritici infection. Our findings contribute to improve our understanding of the interactions between wheat and Z. tritici.

Taxonomically Restricted Wheat Genes Interact With Small Secreted Fungal Proteins and Enhance Resistance to Septoria Tritici Blotch Disease

Understanding the nuances of host/pathogen interactions are paramount if we wish to effectively control cereal diseases. In the case of the wheat/Zymoseptoria tritici interaction that leads to Septoria tritici blotch (STB) disease, a 10,000-year-old conflict has led to considerable armaments being developed on both sides which are not reflected in conventional model systems. Taxonomically restricted genes (TRGs) have evolved in wheat to better allow it to cope with stress caused by fungal pathogens, and Z. tritici has evolved specialized effectors which allow it to manipulate its' host. A microarray focused on the latent phase response of a resistant wheat cultivar (cv. Stigg) and susceptible wheat cultivar (cv. Gallant) to Z. tritici infection was mined for TRGs within the Poaceae. From this analysis, we identified two TRGs that were significantly upregulated in response to Z. tritici infection, Septoria-responsive TRG6 and 7 (TaSRTRG6 and TaSRTRG7). Virus induced silencing of these genes resulted in an increased susceptibility to STB disease in cvs. Gallant and Stigg, and significantly so in the latter (2.5-fold increase in STB disease). In silico and localization studies categorized TaSRTRG6 as a secreted protein and TaSRTRG7 as an intracellular protein. Yeast two-hybrid analysis and biofluorescent complementation studies demonstrated that both TaSRTRG6 and TaSRTRG7 can interact with small proteins secreted by Z. tritici (potential effector candidates). Thus we conclude that TRGs are an important part of the wheat-Z. tritici co-evolution story and potential candidates for modulating STB resistance.

The Zymoseptoria tritici ORFeome: A Functional Genomics Community Resource

Libraries of protein-encoding sequences can be generated by identification of open reading frames (ORFs) from a genome of choice that are then assembled into collections of plasmids termed ORFeome libraries. These represent powerful resources to facilitate functional genomic characterization of genes and their encoded products. Here, we report the generation of an ORFeome for Zymoseptoria tritici, which causes the most serious disease of wheat in temperate regions of the world. We screened the genome of strain IP0323 for high confidence gene models, identifying 4,075 candidates from 10,933 predicted genes. These were amplified from genomic DNA, were cloned into the Gateway entry vector pDONR207, and were sequenced, providing a total of 3,022 quality-controlled plasmids. The ORFeome includes genes predicted to encode effectors (n = 410) and secondary metabolite biosynthetic proteins (n = 171) in addition to genes residing at dispensable chromosomes (n = 122) or those that are preferentially expressed during plant infection (n = 527). The ORFeome plasmid library is compatible with our previously developed suite of Gateway destination vectors, which have various combinations of promoters, selection markers, and epitope tags. The Z. tritici ORFeome constitutes a powerful resource for functional genomics and offers unparalleled opportunities to understand the biology of Z. tritici.[Formula: see text] Copyright © 2019 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Comparative Genomics and Transcriptomics To Analyze Fruiting Body Development in Filamentous Ascomycetes

Many filamentous ascomycetes develop three-dimensional fruiting bodies for production and dispersal of sexual spores. Fruiting bodies are among the most complex structures differentiated by ascomycetes; however, the molecular mechanisms underlying this process are insufficiently understood. Previous comparative transcriptomics analyses of fruiting body development in different ascomycetes suggested that there might be a core set of genes that are transcriptionally regulated in a similar manner across species. Conserved patterns of gene expression can be indicative of functional relevance, and therefore such a set of genes might constitute promising candidates for functional analyses. In this study, we have sequenced the genome of the Pezizomycete Ascodesmis nigricans, and performed comparative transcriptomics of developing fruiting bodies of this fungus, the Pezizomycete Pyronema confluens, and the Sordariomycete Sordaria macrospora With only 27 Mb, the A. nigricans genome is the smallest Pezizomycete genome sequenced to date. Comparative transcriptomics indicated that gene expression patterns in developing fruiting bodies of the three species are more similar to each other than to nonsexual hyphae of the same species. An analysis of 83 genes that are upregulated only during fruiting body development in all three species revealed 23 genes encoding proteins with predicted roles in vesicle transport, the endomembrane system, or transport across membranes, and 13 genes encoding proteins with predicted roles in chromatin organization or the regulation of gene expression. Among four genes chosen for functional analysis by deletion in S. macrospora, three were shown to be involved in fruiting body formation, including two predicted chromatin modifier genes.

Genome-scale data resolve ancestral rock-inhabiting lifestyle in Dothideomycetes (Ascomycota)

Dothideomycetes is the most diverse fungal class in Ascomycota and includes species with a wide range of lifestyles. Previous multilocus studies have investigated the taxonomic and evolutionary relationships of these taxa but often failed to resolve early diverging nodes and frequently generated inconsistent placements of some clades. Here, we use a phylogenomic approach to resolve relationships in Dothideomycetes, focusing on two genera of melanized, extremotolerant rock-inhabiting fungi, Lichenothelia and Saxomyces, that have been suggested to be early diverging lineages. We assembled phylogenomic datasets from newly sequenced (4) and previously available genomes (238) of 242 taxa. We explored the influence of tree inference methods, supermatrix vs. coalescent-based species tree, and the impact of varying amounts of genomic data. Overall, our phylogenetic reconstructions provide consistent and well-supported topologies for Dothideomycetes, recovering Lichenothelia and Saxomyces among the earliest diverging lineages in the class. In addition, many of the major lineages within Dothideomycetes are recovered as monophyletic, and the phylogenomic approach implemented strongly supports their relationships. Ancestral character state reconstruction suggest that the rock-inhabiting lifestyle is ancestral within the class.

Pfcyp51 exclusively determines reduced sensitivity to 14α-demethylase inhibitor fungicides in the banana black Sigatoka pathogen Pseudocercospora fijiensis

The haploid fungus Pseudocercospora fijiensis causes black Sigatoka in banana and is chiefly controlled by extensive fungicide applications, threatening occupational health and the environment. The 14α-Demethylase Inhibitors (DMIs) are important disease control fungicides, but they lose sensitivity in a rather gradual fashion, suggesting an underlying polygenic genetic mechanism. In spite of this, evidence found thus far suggests that P. fijiensis cyp51 gene mutations are the main responsible factor for sensitivity loss in the field. To better understand the mechanisms involved in DMI resistance, in this study we constructed a genetic map using DArTseq markers on two F1 populations generated by crossing two different DMI resistant strains with a sensitive strain. Analysis of the inheritance of DMI resistance in the F₁ populations revealed two major and discrete DMI-sensitivity groups. This is an indicative of a single major responsible gene. Using the DMI-sensitivity scorings of both F1 populations and the generation of genetic linkage maps, the sensitivity causal factor was located in a single genetic region. Full agreement was found for genetic markers in either population, underlining the robustness of the approach. The two maps indicated a similar genetic region where the Pfcyp51 gene is found. Sequence analyses of the Pfcyp51 gene of the F₁ populations also revealed a matching bimodal distribution with the DMI resistant. Amino acid substitutions in P. fijiensis CYP51 enzyme of the resistant progeny were previously correlated with the loss of DMI sensitivity. In addition, the resistant progeny inherited a Pfcyp51 gene promoter insertion, composed of a repeat element with a palindromic core, also previously correlated with increased gene expression. This genetic approach confirms that Pfcyp51 is the single explanatory gene for reduced sensitivity to DMI fungicides in the analysed P. fijiensis strains. Our study is the first genetic analysis to map the underlying genetic factors for reduced DMI efficacy.

Local adaptation drives the diversification of effectors in the fungal wheat pathogen Parastagonospora nodorum in the United States

Filamentous fungi rapidly evolve in response to environmental selection pressures in part due to their genomic plasticity. Parastagonospora nodorum, a fungal pathogen of wheat and causal agent of septoria nodorum blotch, responds to selection pressure exerted by its host, influencing the gain, loss, or functional diversification of virulence determinants, known as effector genes. Whole genome resequencing of 197 P. nodorum isolates collected from spring, durum, and winter wheat production regions of the United States enabled the examination of effector diversity and genomic regions under selection specific to geographically discrete populations. 1,026,859 SNPs/InDels were used to identify novel loci, as well as SnToxA and SnTox3 as factors in disease. Genes displaying presence/absence variation, predicted effector genes, and genes localized on an accessory chromosome had significantly higher pN/pS ratios, indicating a higher rate of sequence evolution. Population structure analyses indicated two P. nodorum populations corresponding to the Upper Midwest (Population 1) and Southern/Eastern United States (Population 2). Prevalence of SnToxA varied greatly between the two populations which correlated with presence of the host sensitivity gene Tsn1 in the most prevalent cultivars in the corresponding regions. Additionally, 12 and 5 candidate effector genes were observed to be under diversifying selection among isolates from Population 1 and 2, respectively, but under purifying selection or neutrally evolving in the opposite population. Selective sweep analysis revealed 10 and 19 regions that had recently undergone positive selection in Population 1 and 2, respectively, involving 92 genes in total. When comparing genes with and without presence/absence variation, those genes exhibiting this variation were significantly closer to transposable elements. Taken together, these results indicate that P. nodorum is rapidly adapting to distinct selection pressures unique to spring and winter wheat production regions by rapid adaptive evolution and various routes of genomic diversification, potentially facilitated through transposable element activity.

The ZtvelB Gene Is Required for Vegetative Growth and Sporulation in the Wheat Pathogen Zymoseptoria tritici

The ascomycete fungus Zymoseptoria tritici is the causal agent of Septoria Tritici Blotch (STB), a major disease of wheat across Europe. Current understanding of the genetic components and the environmental cues which influence development and pathogenicity of this fungus is limited. The velvet B gene, velB, has conserved roles in development, secondary metabolism, and pathogenicity across fungi. The function of this gene is best characterised in the model ascomycete fungus Aspergillus nidulans, where it is involved in co-ordinating the light response with downstream processes. There is limited knowledge of the role of light in Z. tritici, and of the molecular mechanisms underpinning the light response. We show that Z. tritici is able to detect light, and that the vegetative morphology of this fungus is influenced by light conditions. We also identify and characterise the Z. tritici velB gene, ZtvelB, by gene disruption. The ΔztvelB deletion mutants were fixed in a filamentous growth pattern and are unable to form yeast-like vegetative cells. Their morphology was similar under light and dark conditions, showing an impairment in light-responsive growth. In addition, the ΔztvelB mutants produced abnormal pycnidia that were impaired in macropycnidiospore production but could still produce viable infectious micropycnidiospores. Our results show that ZtvelB is required for yeast-like growth and asexual sporulation in Z. tritici, and we provide evidence for a role of ZtvelB in integrating light perception and developmental regulation in this important plant pathogenic fungus.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

The emerging blast disease on wheat is proving even harder to control than the ancient, still-problematic rice blast disease. Potential wheat resistance identified using strains isolated soon after disease emergence are no longer effective in controlling recent aggressive field isolates from wheat in South America and South Asia. We construct a high-quality assembly of an aggressive, recently-isolated wheat blast fungal strain and the first assembled mini-chromosome genome sequence of wheat and rice blast pathogens. We report that recent wheat pathogens can contain one or two highly-variable dispensable mini-chromosomes, each with an amalgamation of fungal effector genes and other sequences that are duplicated or absent from indispensable core chromosome ends. Well-studied effectors found on different core chromosomes in rice pathogens appear side-by-side in wheat pathogen mini-chromosomes. The rice pathogen often overcomes deployed resistance genes by deleting triggering effector genes. We propose that the fast-evolving effector-rich compartment of the wheat blast fungus is a combination of core chromosome ends and mobile mini-chromosomes that are easily lost from individual strains. Localization of effectors on mini-chromosomes would therefore accelerate pathogen adaptation in the field.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

A novel starch-binding laccase from the wheat pathogen Zymoseptoria tritici highlights the functional diversity of ascomycete laccases

BACKGROUND

Laccases are multicopper oxidases, which are assigned into auxiliary activity family 1 (AA1) in the CAZy database. These enzymes, catalyzing the oxidation of phenolic and nonphenolic substrates coupled to reduction of O2 to H2O, are increasingly attractive as eco-friendly oxidation biocatalysts. Basidiomycota laccases are well characterized due to their potential in de-lignification of lignocellulose. By contrast, insight into the biochemical diversity of Ascomycota counterparts from saprophytes and plant pathogens is scarce.

RESULTS

Here, we report the properties of the laccase from the major wheat pathogen Zymoseptoria tritici (ZtrLac1A), distinguished from common plant fungal pathogens by an apoplastic infection strategy. We demonstrate that ZtrLac1A is appended to a functional starch-binding module and displays an activity signature disfavoring relatively apolar phenolic redox mediators as compared to the related biochemically characterized laccases. By contrast, the redox potential of ZtrLac1A (370 mV vs. SHE) is similar to ascomycetes counterparts. The atypical specificity is consistent with distinctive sequence substitutions and insertions in loops flanking the T1 site and the enzyme C-terminus compared to characterized laccases.

CONCLUSIONS

ZtrLac1A is the first reported modular laccase appended to a functional starch-specific carbohydrate binding module of family 20 (CBM20). The distinct specificity profile of ZtrLac1A correlates to structural differences in the active site region compared to previously described ascomycetes homologues. These differences are also highlighted by the clustering of the sequence of ZtrLac1A in a distinct clade populated predominantly by plant pathogens in the phylogenetic tree of AA1 laccases. The possible role of these laccases in vivo merits further investigations. These findings expand our toolbox of laccases for green oxidation and highlight the binding functionality of CBM-appended laccases as versatile immobilization tags.

sRNA Profiling Combined With Gene Function Analysis Reveals a Lack of Evidence for Cross-Kingdom RNAi in the Wheat - Zymoseptoria tritici Pathosystem

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonization and disease development. Here we characterized RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assessed their contribution to pathogenesis. Computational analysis of fungal sRNA and host mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict their mRNA targets in wheat. 389 in planta-induced sRNA loci were identified. sRNAs generated from some of these loci were predicted to target wheat mRNAs including those potentially involved in pathogen defense. However, molecular approaches failed to validate targeting of selected wheat mRNAs by fungal sRNAs. Mutant strains of Z. tritici carrying deletions of genes encoding key components of RNAi such as Dicer-like (DCL) and Argonaute (AGO) proteins were generated, and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results did suggest the existence of non-canonical DCL-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA targeting essential fungal genes applied in vitro or generated from an RNA virus vector in planta in a procedure known as HIGS (Host-Induced Gene Silencing) was ineffective in preventing Z. tritici growth or disease. We also demonstrated that Z. tritici is incapable of dsRNA uptake. Collectively, our data suggest that RNAi approaches for gene function analyses in this fungal species and potentially also as a control measure may not be as effective as has been demonstrated for some other plant pathogenic fungi.

Acquisition and Loss of Secondary Metabolites Shaped the Evolutionary Path of Three Emerging Phytopathogens of Wheat

White grain disorder is a recently emerged wheat disease in Australia, caused by Eutiarosporella darliae, E. pseudodarliae, and E. tritici-australis. The disease cycle of these pathogens and the molecular basis of their interaction with wheat are poorly understood. To address this knowledge gap, we undertook a comparative genomics analysis focused on the secondary metabolite gene repertoire among these three species. This analysis revealed a diverse array of secondary metabolite gene clusters in these pathogens, including modular polyketide synthase genes. These genes have only been previously associated with bacteria and this is the first report of such genes in fungi. Subsequent phylogenetic analyses provided strong evidence that the modular PKS genes were horizontally acquired from a bacterial or a protist species. We also uncovered a secondary metabolite gene cluster with three polyketide/nonribosomal peptide synthase genes (Hybrid-1, -2, and -3) in E. darliae and E. pseudodarliae. In contrast, only remnant and partial genes homologous to this cluster were identified in E. tritici-australis, suggesting loss of this cluster. Homologues of Hybrid-2 in other fungi have been proposed to facilitate disease in woody plants, suggesting a possible alternative host range for E. darliae and E. pseudodarliae. Subsequent assays confirmed that E. darliae and E. pseudodarliae were both pathogenic on woody plants, but E. tritici-australis was not, implicating woody plants as potential host reservoirs for the fungi. Combined, these data have advanced our understanding of the lifestyle and potential host-range of these recently emerged wheat pathogens and shed new light on fungal secondary metabolism.

Comparative genomics of chytrid fungi reveal insights into the obligate biotrophic and pathogenic lifestyle of Synchytrium endobioticum

Synchytrium endobioticum is an obligate biotrophic soilborne Chytridiomycota (chytrid) species that causes potato wart disease, and represents the most basal lineage among the fungal plant pathogens. We have chosen a functional genomics approach exploiting knowledge acquired from other fungal taxa and compared this to several saprobic and pathogenic chytrid species. Observations linked to obligate biotrophy, genome plasticity and pathogenicity are reported. Essential purine pathway genes were found uniquely absent in S. endobioticum, suggesting that it relies on scavenging guanine from its host for survival. The small gene-dense and intron-rich chytrid genomes were not protected for genome duplications by repeat-induced point mutation. Both pathogenic chytrids Batrachochytrium dendrobatidis and S. endobioticum contained the largest amounts of repeats, and we identified S. endobioticum specific candidate effectors that are associated with repeat-rich regions. These candidate effectors share a highly conserved motif, and show isolate specific duplications. A reduced set of cell wall degrading enzymes, and LysM protein expansions were found in S. endobioticum, which may prevent triggering plant defense responses. Our study underlines the high diversity in chytrids compared to the well-studied Ascomycota and Basidiomycota, reflects characteristic biological differences between the phyla, and shows commonalities in genomic features among pathogenic fungi.

Progress and Challenges: Development and Implementation of CRISPR/Cas9 Technology in Filamentous Fungi

Widely distributed in various environmental niches, filamentous fungi play an important role in industry, drug development, and plant/animal health. Manipulation of the genome and the coding sequences are essential for a better understanding of the function of genes and their regulation, but traditional genetic approaches in some filamentous fungi are either inefficient or nonfunctional. The rapid development and wide implementation of CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats /(CRISPR)-associated protein-9 nuclease) technology for various model and non-model organisms has provided the initial framework to adapt this gene editing technology for filamentous fungi. In this review, an overview of the CRISPR/Cas9 tools and strategies that have been developed for different filamentous fungi is presented, including integration of the CAS9 gene into the genome, transient expression of Cas9/sgRNA, the AMA1-based plasmid approach, and the Cas9 RNP method. The various applications of CRISPR/Cas9 technology in filamentous fungi that have been implemented are explored, with particular emphasis on gene disruption/deletion and precise genome modification through gene tagging and alteration in gene regulation. Potential challenges that are confronted when developing a CRISPR/Cas9 system for filamentous fungi are also discussed such as the nuclear localization sequence for the CAS9 gene, potential off-target effects, and highly efficient transformation methods. Overcoming these obstacles may further facilitate wide application of this technology. As a simple, economical, and powerful tool, CRISPR/Cas9 systems have the potential for future implementation into many molecular aspects of filamentous fungi.

The genomic determinants of adaptive evolution in a fungal pathogen

Unravelling the strength, frequency, and distribution of selective variants along the genome as well as the underlying factors shaping this distribution are fundamental goals of evolutionary biology. Antagonistic host-pathogen coevolution is thought to be a major driver of genome evolution between interacting species. While rapid evolution of pathogens has been documented in several model organisms, the genetic mechanisms of their adaptation are still poorly understood and debated, particularly the role of sexual reproduction. Here, we apply a population genomic approach to infer genome-wide patterns of selection among 13 isolates of Zymoseptoria tritici, a fungal pathogen characterized by extremely high genetic diversity, gene density, and recombination rates. We report that the genome of Z. tritici undergoes a high rate of adaptive substitutions, with 44% of nonsynonymous substitutions being adaptive on average. This fraction reaches 68% in so-called effector genes encoding determinants of pathogenicity, and the distribution of fitness effects differs in this class of genes as they undergo adaptive mutations with stronger positive fitness effects, but also more slightly deleterious mutations. Besides the globally high rate of adaptive substitutions, we report a negative relationship between pN/pS and the fine-scale recombination rate and a strong positive correlation between the rate of adaptive nonsynonymous substitutions (ωa) and recombination rate. This result suggests a pervasive role of both background selection and Hill-Robertson interference even in a species with an exceptionally high recombination rate (60 cM/Mb on average). While transposable elements (TEs) have been suggested to contribute to adaptation by creating compartments of fast-evolving genomic regions, we do not find a significant effect of TEs on the rate of adaptive mutations. Overall our study suggests that sexual recombination is a significant driver of genome evolution, even in rapidly evolving organisms subject to recurrent mutations with large positive effects.

Global population genomics of the forest pathogen Dothistroma septosporum reveal chromosome duplications in high dothistromin-producing strains

Dothistroma needle blight is one of the most devastating pine tree diseases worldwide. New and emerging epidemics have been frequent over the last 25 years, particularly in the Northern Hemisphere, where they are in part associated with changing weather patterns. One of the main Dothistroma needle blight pathogens, Dothistroma septosporum, has a global distribution but most molecular plant pathology research has been confined to Southern Hemisphere populations that have limited genetic diversity. Extensive genomic and transcriptomic data are available for a D. septosporum reference strain from New Zealand, where an introduced clonal population of the pathogen predominates. Due to the global importance of this pathogen, we determined whether the genome of this reference strain is representative of the species worldwide by sequencing the genomes of 18 strains sampled globally from different pine hosts. Genomic polymorphism shows substantial variation within the species, clustered into two distinct groups of strains with centres of diversity in Central and South America. A reciprocal chromosome translocation uniquely identifies the New Zealand strains. Globally, strains differ in their production of the virulence factor dothistromin, with extremely high production levels in strain ALP3 from Germany. Comparisons with the New Zealand reference revealed that several strains are aneuploids; for example, ALP3 has duplications of three chromosomes. Increased gene copy numbers therefore appear to contribute to increased production of dothistromin, emphasizing that studies of population structure are a necessary adjunct to functional analyses of genetic polymorphisms to identify the molecular basis of virulence in this important forest pathogen.

De novo genome assembly of the cichlid fish Astatotilapia latifasciata reveals a higher level of genomic polymorphism and genes related to B chromosomes

Supernumerary B chromosomes (Bs) are accessory elements to the regular chromosome set (As) and have been observed in a huge diversity of eukaryotic species. Although extensively investigated, the biological significance of Bs remains enigmatic. Here, we present de novo genome assemblies for the cichlid fish Astatotilapia latifasciata, a well-known model to study Bs. High coverage data with Illumina sequencing was obtained for males and females with 0B (B-), 1B, and 2B (B+) chromosomes to provide information regarding the diversity among these genomes. The draft assemblies comprised 771 Mb for the B- genome and 781 Mb for the B+ genome. Comparative analysis of the B+ and B- assemblies reveals syntenic discontinuity, duplicated blocks and several insertions, deletions, and inversions indicative of rearrangements in the B+ genome. Hundreds of transposable elements and 1546 protein coding sequences were annotated in the duplicated B+ regions. Our work contributes a list of thousands of genes harbored on the B chromosome, with functions in several biological processes, including the cell cycle.

Genomic analysis of the aggressive tree pathogen Ceratocystis albifundus

The overall goal of this study was to determine whether the genome of an important plant pathogen in Africa, Ceratocystis albifundus, is structured into subgenomic compartments, and if so, to establish how these compartments are distributed across the genome. For this purpose, the publicly available genome of C. albifundus was complemented with the genome sequences for four additional isolates using the Illumina HiSeq platform. In addition, a reference genome for one of the individuals was assembled using both PacBio and Illumina HiSeq technologies. Our results showed a high degree of synteny between the five genomes, although several regions lacked detectable long-range synteny. These regions were associated with the presence of accessory genes, lower genetic similarity, variation in read-map depth, as well as transposable elements and genes associated with host-pathogen interactions (e.g. effectors and CAZymes). Such patterns are regarded as hallmarks of accelerated evolution, particularly of accessory subgenomic compartments in fungal pathogens. Our findings thus showed that the genome of C. albifundus is made-up of core and accessory subgenomic compartments, which is an important step towards characterizing its pangenome. This study also highlights the value of comparative genomics for understanding mechanisms that may underly and influence the biology and evolution of pathogens.

Destabilization of chromosome structure by histone H3 lysine 27 methylation

Chromosome and genome stability are important for normal cell function as instability often correlates with disease and dysfunction of DNA repair mechanisms. Many organisms maintain supernumerary or accessory chromosomes that deviate from standard chromosomes. The pathogenic fungus Zymoseptoria tritici has as many as eight accessory chromosomes, which are highly unstable during meiosis and mitosis, transcriptionally repressed, show enrichment of repetitive elements, and enrichment with heterochromatic histone methylation marks, e.g., trimethylation of H3 lysine 9 or lysine 27 (H3K9me3, H3K27me3). To elucidate the role of heterochromatin on genome stability in Z. tritici, we deleted the genes encoding the methyltransferases responsible for H3K9me3 and H3K27me3, kmt1 and kmt6, respectively, and generated a double mutant. We combined experimental evolution and genomic analyses to determine the impact of these deletions on chromosome and genome stability, both in vitro and in planta. We used whole genome sequencing, ChIP-seq, and RNA-seq to compare changes in genome and chromatin structure, and differences in gene expression between mutant and wildtype strains. Analyses of genome and ChIP-seq data in H3K9me3-deficient strains revealed dramatic chromatin reorganization, where H3K27me3 is mostly relocalized into regions that are enriched with H3K9me3 in wild type. Many genome rearrangements and formation of new chromosomes were found in the absence of H3K9me3, accompanied by activation of transposable elements. In stark contrast, loss of H3K27me3 actually increased the stability of accessory chromosomes under normal growth conditions in vitro, even without large scale changes in gene activity. We conclude that H3K9me3 is important for the maintenance of genome stability because it disallows H3K27me3 in regions considered constitutive heterochromatin. In this system, H3K27me3 reduces the overall stability of accessory chromosomes, generating a "metastable" state for these quasi-essential regions of the genome.

Broad-specificity GH131 β-glucanases are a hallmark of fungi and oomycetes that colonize plants

Plant-tissue-colonizing fungi fine-tune the deconstruction of plant-cell walls (PCW) using different sets of enzymes according to their lifestyle. However, some of these enzymes are conserved among fungi with dissimilar lifestyles. We identified genes from Glycoside Hydrolase family GH131 as commonly expressed during plant-tissue colonization by saprobic, pathogenic and symbiotic fungi. By searching all the publicly available genomes, we found that GH131-coding genes were widely distributed in the Dikarya subkingdom, except in Taphrinomycotina and Saccharomycotina, and in phytopathogenic Oomycetes, but neither other eukaryotes nor prokaryotes. The presence of GH131 in a species was correlated with its association with plants as symbiont, pathogen or saprobe. We propose that GH131-family expansions and horizontal-gene transfers contributed to this adaptation. We analysed the biochemical activities of GH131 enzymes whose genes were upregulated during plant-tissue colonization in a saprobe (Pycnoporus sanguineus), a plant symbiont (Laccaria bicolor) and three hemibiotrophic-plant pathogens (Colletotrichum higginsianum, C. graminicola, Zymoseptoria tritici). These enzymes were all active on substrates with β-1,4, β-1,3 and mixed β-1,4/1,3 glucosidic linkages. Combined with a cellobiohydrolase, GH131 enzymes enhanced cellulose degradation. We propose that secreted GH131 enzymes unlock the PCW barrier and allow further deconstruction by other enzymes during plant tissue colonization by symbionts, pathogens and saprobes.

From pathogen to endophyte: an endophytic population of Verticillium dahliae evolved from a sympatric pathogenic population

The fungus Verticillium dahliae causes wilts of several hundred plant species, including potato and mint. Verticillium spp. also colonize sympatric hosts such as mustards and grasses as endophytes. The evolutionary history of and interactions between pathogenic and endophytic of this fungus are unknown. Verticillium dahliae isolates recovered from sympatric potato, mint, mustard and grasses were characterized genotypically with microsatellite markers and phenotypically for pathogenicity. The evolutionary history of pathogenic and endophytic populations was reconstructed and gene flow between populations quantified. Verticillium dahliae was recovered from all hosts. Endophytic populations were genetically and genotypically similar to but marginally differentiated from the potato population, from which they evolved. Bidirectional migration was detected between these populations and endophytic isolates were pathogenic to potato and behaved as endophytes in mustard and barley. Verticillium dahliae colonizes plants as both endophytes and pathogens. A historical host-range expansion together with endophytic and pathogenic capabilities are likely to have enabled infection of and gene flow between asymptomatic and symptomatic host populations despite minor differentiation. The ability of hosts to harbor asymptomatic infections and the stability of asymptomatic infections over time warrants investigation to elucidate the mechanisms involved in the maintenance of endophytism and pathogenesis.

Decoupling a novel Trichormus variabilis-Synechocystis sp. interaction to boost phycoremediation

To conserve freshwater resources, domestic and industrial wastewater is recycled. Algal systems have emerged as an efficient, low-cost option for treatment (phycoremediation) of nutrient-rich wastewater and environmental protection. However, industrial wastewater may contain growth inhibitory compounds precluding algal use in phycoremediation. Therefore, extremophyte strains, which thrive in hostile environments, are sought-after. Here, we isolated such an alga - a strain of Synechocystis sp. we found to be capable of switching from commensal exploitation of the nitrogen-fixing Trichormus variabilis, for survival in nitrogen-deficient environments, to free-living growth in nitrate abundance. In nitrogen depletion, the cells are tethered to polysaccharide capsules of T. variabilis using nanotubular structures, presumably for nitrate acquisition. The composite culture failed to establish in industrial/domestic waste effluent. However, gradual exposure to increasing wastewater strength over time untethered Synechocystis cells and killed off T. variabilis. This switched the culture to a stress-acclimated monoculture of Synechocystis sp., which rapidly grew and flourished in wastewater, with ammonium and phosphate removal efficiencies of 99.4% and 97.5%, respectively. Therefore, this strain of Synechocystis sp. shows great promise for use in phycoremediation, with potential to rapidly generate biomass that can find use as a green feedstock for valuable bio-products in industrial applications.

Sequence Expression of Supernumerary B Chromosomes: Function or Fluff?

B chromosomes are enigmatic heritable elements found in the genomes of numerous plant and animal species. Contrary to their broad distribution, most B chromosomes are non-essential. For this reason, they are regarded as genome parasites. In order to be stably transmitted through generations, many B chromosomes exhibit the ability to "drive", i.e., they transmit themselves at super-Mendelian frequencies to progeny through directed interactions with the cell division apparatus. To date, very little is understood mechanistically about how B chromosomes drive, although a likely scenario is that expression of B chromosome sequences plays a role. Here, we highlight a handful of previously identified B chromosome sequences, many of which are repetitive and non-coding in nature, that have been shown to be expressed at the transcriptional level. We speculate on how each type of expressed sequence could participate in B chromosome drive based on known functions of RNA in general chromatin- and chromosome-related processes. We also raise some challenges to functionally testing these possible roles, a goal that will be required to more fully understand whether and how B chromosomes interact with components of the cell for drive and transmission.

Highly flexible infection programs in a specialized wheat pathogen

Many filamentous plant pathogens exhibit high levels of genomic variability, yet the impact of this variation on host-pathogen interactions is largely unknown. We have addressed host specialization in the wheat pathogen Zymoseptoria tritici. Our study builds on comparative analyses of infection and gene expression phenotypes of three isolates and reveals the extent to which genomic variation translates into phenotypic variation. The isolates exhibit genetic and genomic variation but are similarly virulent. By combining confocal microscopy, disease monitoring, staining of ROS, and comparative transcriptome analyses, we conducted a detailed comparison of the infection processes of these isolates in a susceptible wheat cultivar. We characterized four core infection stages: establishment, biotrophic growth, lifestyle transition, and necrotrophic growth and asexual reproduction that are shared by the three isolates. However, we demonstrate differentiated temporal and spatial infection development and significant differences in the expression profiles of the three isolates during the infection stages. More than 20% of the genes were differentially expressed and these genes were located significantly closer to transposable elements, suggesting an impact of epigenetic regulation. Further, differentially expressed genes were enriched in effector candidates suggesting that isolate-specific strategies for manipulating host defenses are present in Z. tritici. We demonstrate that individuals of a host-specialized pathogen have highly differentiated infection programs characterized by flexible infection development and functional redundancy. This illustrates how high genetic diversity in pathogen populations results in highly differentiated infection phenotypes, which fact needs to be acknowledged to understand host-pathogen interactions and pathogen evolution.

Meiotic drive of female-inherited supernumerary chromosomes in a pathogenic fungus

Meiosis is a key cellular process of sexual reproduction that includes pairing of homologous sequences. In many species however, meiosis can also involve the segregation of supernumerary chromosomes, which can lack a homolog. How these unpaired chromosomes undergo meiosis is largely unknown. In this study we investigated chromosome segregation during meiosis in the haploid fungus Zymoseptoria tritici that possesses a large complement of supernumerary chromosomes. We used isogenic whole chromosome deletion strains to compare meiotic transmission of chromosomes when paired and unpaired. Unpaired chromosomes inherited from the male parent as well as paired supernumerary chromosomes in general showed Mendelian inheritance. In contrast, unpaired chromosomes inherited from the female parent showed non-Mendelian inheritance but were amplified and transmitted to all meiotic products. We concluded that the supernumerary chromosomes of Z. tritici show a meiotic drive and propose an additional feedback mechanism during meiosis, which initiates amplification of unpaired female-inherited chromosomes.

Gene exchange between two divergent species of the fungal human pathogen, Coccidioides

The fungal genus Coccidioides is composed of two species, Coccidioides immitis and Coccidioides posadasii. These two species are the causal agents of coccidioidomycosis, a pulmonary disease also known as valley fever. The two species are thought to have shared genetic material due to gene exchange in spite of their long divergence. To quantify the magnitude of shared ancestry between them, we analyzed the genomes of a population sample from each species. Next, we inferred what is the expected size of shared haplotypes that might be inherited from the last common ancestor of the two species and find a cutoff to find what haplotypes have conclusively been exchanged between species. Finally, we precisely identified the breakpoints of the haplotypes that have crossed the species boundary and measure the allele frequency of each introgression in this sample. We find that introgressions are not uniformly distributed across the genome. Most, but not all, of the introgressions segregate at low frequency. Our results show that divergent species can share alleles, that species boundaries can be porous, and highlight the need for a systematic exploration of gene exchange in fungal species.

High intraspecific genome diversity in the model arbuscular mycorrhizal symbiont Rhizophagus irregularis

Arbuscular mycorrhizal fungi (AMF) are known to improve plant fitness through the establishment of mycorrhizal symbioses. Genetic and phenotypic variations among closely related AMF isolates can significantly affect plant growth, but the genomic changes underlying this variability are unclear. To address this issue, we improved the genome assembly and gene annotation of the model strain Rhizophagus irregularis DAOM197198, and compared its gene content with five isolates of R. irregularis sampled in the same field. All isolates harbor striking genome variations, with large numbers of isolate-specific genes, gene family expansions, and evidence of interisolate genetic exchange. The observed variability affects all gene ontology terms and PFAM protein domains, as well as putative mycorrhiza-induced small secreted effector-like proteins and other symbiosis differentially expressed genes. High variability is also found in active transposable elements. Overall, these findings indicate a substantial divergence in the functioning capacity of isolates harvested from the same field, and thus their genetic potential for adaptation to biotic and abiotic changes. Our data also provide a first glimpse into the genome diversity that resides within natural populations of these symbionts, and open avenues for future analyses of plant-AMF interactions that link AMF genome variation with plant phenotype and fitness.

A New Reference Genome Shows the One-Speed Genome Structure of the Barley Pathogen Ramularia collo-cygni

Ramularia leaf spot has recently emerged as a major threat to barley production world-wide, causing 25% yield loss in many barley growing regions. Here, we provide a new reference genome of the causal agent, the Dothideomycete Ramularia collo-cygni. The assembly of 32 Mb consists of 78 scaffolds. We used RNA-seq to identify 11,622 genes of which 1,303 and 282 are coding for predicted secreted proteins and putative effectors respectively.

The pathogen separated from its nearest sequenced relative, Zymoseptoria tritici ∼27 Ma. We calculated the divergence of the two species on protein level and see remarkably high synonymous and nonsynonymous divergence. Unlike in many other plant pathogens, the comparisons of transposable elements and gene distributions, show a very homogeneous genome for R. collo-cygni. We see no evidence for higher selective pressure on putative effectors or other secreted proteins and repetitive sequences are spread evenly across the scaffolds. These findings could be associated to the predominantly endophytic life-style of the pathogen. We hypothesize that R. collo-cygni only recently became pathogenic and that therefore its genome does not yet show the typical pathogen characteristics. Because of its high scaffold length and improved CDS annotations, our new reference sequence provides a valuable resource for the community for future comparative genomics and population genetics studies.

Phosphopantetheinyl transferase (Ppt)-mediated biosynthesis of lysine, but not siderophores or DHN melanin, is required for virulence of Zymoseptoria tritici on wheat

Zymoseptoria tritici is the causal agent of Septoria tritici blotch (STB) disease of wheat. Z. tritici is an apoplastic fungal pathogen, which does not penetrate plant cells at any stage of infection, and has a long initial period of symptomless leaf colonisation. During this phase it is unclear to what extent the fungus can access host plant nutrients or communicate with plant cells. Several important primary and secondary metabolite pathways in fungi are regulated by the post-translational activator phosphopantetheinyl transferase (Ppt) which provides an essential co-factor for lysine biosynthesis and the activities of non-ribosomal peptide synthases (NRPS) and polyketide synthases (PKS). To investigate the relative importance of lysine biosynthesis, NRPS-based siderophore production and PKS-based DHN melanin biosynthesis, we generated deletion mutants of ZtPpt. The ∆ZtPpt strains were auxotrophic for lysine and iron, non-melanised and non-pathogenic on wheat. Deletion of the three target genes likely affected by ZtPpt loss of function (Aar- lysine; Nrps1-siderophore and Pks1- melanin), highlighted that lysine auxotrophy was the main contributing factor for loss of virulence, with no reduction caused by loss of siderophore production or melanisation. This reveals Ppt, and the lysine biosynthesis pathway, as potential targets for fungicides effective against Z. tritici.

Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism.