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

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.

Extraordinary Genome Instability and Widespread Chromosome Rearrangements During Vegetative Growth

The haploid genome of the pathogenic fungus Zymoseptoria tritici is contained on "core" and "accessory" chromosomes. While 13 core chromosomes are found in all strains, as many as eight accessory chromosomes show presence/absence variation and rearrangements among field isolates. The factors influencing these presence/absence polymorphisms are so far unknown. We investigated chromosome stability using experimental evolution, karyotyping, and genome sequencing. We report extremely high and variable rates of accessory chromosome loss during mitotic propagation in vitro and in planta Spontaneous chromosome loss was observed in 2 to >50% of cells during 4 weeks of incubation. Similar rates of chromosome loss in the closely related Zymoseptoria ardabiliae suggest that this extreme chromosome dynamic is a conserved phenomenon in the genus. Elevating the incubation temperature greatly increases instability of accessory and even core chromosomes, causing severe rearrangements involving telomere fusion and chromosome breakage. Chromosome losses do not affect the fitness of Zymoseptoria tritici in vitro, but some lead to increased virulence, suggesting an adaptive role of this extraordinary chromosome instability.

Pan-Parastagonospora Comparative Genome Analysis-Effector Prediction and Genome Evolution

We report a fungal pan-genome study involving Parastagonospora spp., including 21 isolates of the wheat (Triticum aestivum) pathogen Parastagonospora nodorum, 10 of the grass-infecting Parastagonospora avenae, and 2 of a closely related undefined sister species. We observed substantial variation in the distribution of polymorphisms across the pan-genome, including repeat-induced point mutations, diversifying selection and gene gains and losses. We also discovered chromosome-scale inter and intraspecific presence/absence variation of some sequences, suggesting the occurrence of one or more accessory chromosomes or regions that may play a role in host-pathogen interactions. The presence of known pathogenicity effector loci SnToxA, SnTox1, and SnTox3 varied substantially among isolates. Three P. nodorum isolates lacked functional versions for all three loci, whereas three P. avenae isolates carried one or both of the SnTox1 and SnTox3 genes, indicating previously unrecognized potential for discovering additional effectors in the P. nodorum-wheat pathosystem. We utilized the pan-genomic comparative analysis to improve the prediction of pathogenicity effector candidates, recovering the three confirmed effectors among our top-ranked candidates. We propose applying this pan-genomic approach to identify the effector repertoire involved in other host-microbe interactions involving necrotrophic pathogens in the Pezizomycotina.

Genetic Structure of Zymoseptoria tritici in Northern France at Region, Field, Plant, and Leaf Layer Scales

Population genetic structure of the worldwide-distributed wheat pathogen Zymoseptoria tritici has been extensively studied at large geographical scales, but to a much less extent at small or local spatial scales. A total of 627 single-conidial fungal isolates were sampled from several locations in northern France (Hauts-de-France Region) to assess fungal genetic structure at region, field, plant, and leaf layer scales, using highly polymorphic microsatellite markers and mating type idiomorphs. Important and overall similar levels of both gene and genotype diversities (gene diversity values of ≥0.44 and haplotype frequencies of ≥94%) were found at all the examined scales. Such rates of diversity are likely due to an active sexual recombination in the investigated areas, as revealed by equal proportions of the two mating types scored in all sampled populations. Interestingly, a rare occurrence of clones among lesions from a same leaf, as well as among leaves from different plant leaf layers (e.g., upper versus lower leaves), was highlighted, indicating that ascospores contribute much more than expected to Z. tritici epidemics, compared with pycnidiospores. Population structure and analyses of molecular variance revealed significant genetic differentiation at the regional scale (G_ST = 0.23) and, as expected, not at the other more local scales (G_ST ≤ 0.01). Further analyses using Bayesian and unweighted neighbor-joining statistical methods detected six genetic clusters within the regional population, overall distributed according to the locations from which the isolates were sampled. Neither clear directional relative migration linked to the geographical distribution of the locations, nor isolation by distance, were observed. Separate evolutionary trajectories caused by selection and adaptations to habitat heterogeneity could be the main forces shaping such structuration. This study provides new insights into the epidemiology and the genetic structure of Z. tritici at small local and, for the first time, at single plant and leaf layer scales. Such findings would be helpful in implementing effective control strategies.

Improved prediction of fungal effector proteins from secretomes with EffectorP 2.0

Plant-pathogenic fungi secrete effector proteins to facilitate infection. We describe extensive improvements to EffectorP, the first machine learning classifier for fungal effector prediction. EffectorP 2.0 is now trained on a larger set of effectors and utilizes a different approach based on an ensemble of classifiers trained on different subsets of negative data, offering different views on classification. EffectorP 2.0 achieves an accuracy of 89%, compared with 82% for EffectorP 1.0 and 59.8% for a small size classifier. Important features for effector prediction appear to be protein size, protein net charge as well as the amino acids serine and cysteine. EffectorP 2.0 decreases the number of predicted effectors in secretomes of fungal plant symbionts and saprophytes by 40% when compared with EffectorP 1.0. However, EffectorP 1.0 retains value, and combining EffectorP 1.0 and 2.0 results in a stringent classifier with a low false positive rate of 9%. EffectorP 2.0 predicts significant enrichments of effectors in 12 of 13 sets of infection-induced proteins from diverse fungal pathogens, whereas a small cysteine-rich classifier detects enrichment in only seven of 13. EffectorP 2.0 will fast track the prioritization of high-confidence effector candidates for functional validation and aid in improving our understanding of effector biology. EffectorP 2.0 is available at http://effectorp.csiro.au.

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.

A fungal avirulence factor encoded in a highly plastic genomic region triggers partial resistance to septoria tritici blotch

Cultivar-strain specificity in the wheat-Zymoseptoria tritici pathosystem determines the infection outcome and is controlled by resistance genes on the host side, many of which have been identified. On the pathogen side, however, the molecular determinants of specificity remain largely unknown. We used genetic mapping, targeted gene disruption and allele swapping to characterise the recognition of the new avirulence factor Avr3D1. We then combined population genetic and comparative genomic analyses to characterise the evolutionary trajectory of Avr3D1. Avr3D1 is specifically recognised by wheat cultivars harbouring the Stb7 resistance gene, triggering a strong defence response without preventing pathogen infection and reproduction. Avr3D1 resides in a cluster of putative effector genes located in a genome region populated by independent transposable element insertions. The gene was present in all 132 investigated strains and is highly polymorphic, with 30 different protein variants identified. We demonstrated that specific amino acid substitutions in Avr3D1 led to evasion of recognition. These results demonstrate that quantitative resistance and gene-for-gene interactions are not mutually exclusive. Localising avirulence genes in highly plastic genomic regions probably facilitates accelerated evolution that enables escape from recognition by resistance proteins.

The MSDIN family in amanitin-producing mushrooms and evolution of the prolyl oligopeptidase genes

The biosynthetic pathway for amanitins and related cyclic peptides in deadly Amanita (Amanitaceae) mushrooms represents the first known ribosomal cyclic peptide pathway in the Fungi. Amanitins are found outside of the genus in distantly related agarics Galerina (Strophariaceae) and Lepiota (Agaricaceae). A long-standing question in the field persists: why is this pathway present in these phylogenetically disjunct agarics? Two deadly mushrooms, A. pallidorosea and A. subjunquillea, were deep sequenced, and sequences of biosynthetic genes encoding MSDINs (cyclic peptide precursor) and prolyl oligopeptidases (POPA and POPB) were obtained. The two Amanita species yielded 29 and 18 MSDINs, respectively. In addition, two MSDIN sequences were cloned from L. brunneoincarnata basidiomes. The toxin MSDIN genes encoding amatoxins or phallotoxins from the three genera were compared, and a phylogenetic tree constructed. Prolyl oligopeptidase B (POPB), a key enzyme in the biosynthetic pathway, was used in phylogenetic reconstruction to infer the evolutionary history of the genes. Phylogenies of POPB and POPA based on both coding and amino acid sequences showed very different results: while POPA genes clearly reflected the phylogeny of the host species, POPB did not; strikingly, it formed a well-supported monophyletic clade, despite that the species belong to different genera in disjunct families. POPA, a known house-keeping gene, was shown to be restricted in a branch containing only Amanita species and the phylogeny resembled that of those Amanita species. Phylogenetic analyses of MSDIN and POPB genes showed tight coordination and disjunct distribution. A POPB gene tree was compared with a corresponding species tree, and distances and substitution rates were compared. The result suggested POPB genes have significant smaller distances and rates than the house-keeping rpb2, discounting massive gene loss. Under this assumption, the incongruency between the gene tree and species tree was shown with strong support. Additionally, k-mer analyses consistently cluster Galerina and Amanita POPB genes, while Lepiota POPB is distinct. Our result suggests that horizontal gene transfer (HGT), at least between Amanita and Galerina, was involved in the acquisition of POPB genes, which may shed light on the evolution of the α-amanitin biosynthetic pathway.

Colletotrichum higginsianum as a Model for Understanding Host⁻Pathogen Interactions: A Review

Colletotrichum higginsianum is a hemibiotrophic ascomycetous fungus that causes economically important anthracnose diseases on numerous monocot and dicot crops worldwide. As a model pathosystem, the Colletotrichum⁻Arabidopsis interaction has the significant advantage that both organisms can be manipulated genetically. The goal of this review is to provide an overview of the system and to point out recent significant studies that update our understanding of the pathogenesis of C. higginsianum and resistance mechanisms of Arabidopsis against this hemibiotrophic fungus. The genome sequence of C. higginsianum has provided insights into how genome structure and pathogen genetic variability has been shaped by transposable elements, and allows systematic approaches to longstanding areas of investigation, including infection structure differentiation and fungal⁻plant interactions. The Arabidopsis-Colletotrichum pathosystem provides an integrated system, with extensive information on the host plant and availability of genomes for both partners, to illustrate many of the important concepts governing fungal⁻plant interactions, and to serve as an excellent starting point for broad perspectives into issues in plant pathology.

Annotation resource of tandem repeat-containing secretory proteins in sixty fungi

Fungal secretory proteins that interact with host plants are regarded as effectors. Because fungal effectors rarely contain conserved sequence features, identification and annotation of fungal effectors from predicted secretory proteins are difficult using outward comparison methods such as BLAST or hidden Markov model. In desire of more sequence features to prioritize research interests of fungal secretory proteins, this study developed a pipeline to identify tandem repeat (TR) domain within putative secretory proteins and tested a hypothesis that at least one type of TR domain in non-orthologous secretory proteins has emerged from convergent evolution for plant pathogenicity. There were 2804 types of TR domains and a total of 2925 TR-containing secretory proteins found from 60 fungi. There was no conserved type of TR domain shared only by plant pathogens, indicating functional divergence for different types of TR domain and TR-containing secretory proteins. The annotation resource of putative fungal TR-containing secretory proteins provides new sequence features that will be useful for the community interested in fungal effector biology.

Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat

BACKGROUND

Fungal plant pathogens pose major threats to crop yield and sustainable food production if they are highly adapted to their host and the local environment. Variation in gene expression contributes to phenotypic diversity within fungal species and affects adaptation. However, very few cases of adaptive regulatory changes have been reported in fungi and the underlying mechanisms remain largely unexplored. Fungal pathogen genomes are highly plastic and harbor numerous insertions of transposable elements, which can potentially contribute to gene expression regulation. In this work, we elucidated how transposable elements contribute to variation in melanin accumulation, a quantitative trait in fungi that affects survival under stressful conditions.

RESULTS

We demonstrated that differential transcriptional regulation of the gene encoding the transcription factor Zmr1, which controls expression of the genes in the melanin biosynthetic gene cluster, is responsible for variation in melanin accumulation in the fungal plant pathogen Zymoseptoria tritici. We show that differences in melanin levels between two strains of Z. tritici are due to two levels of transcriptional regulation: (1) variation in the promoter sequence of Zmr1 and (2) an insertion of transposable elements upstream of the Zmr1 promoter. Remarkably, independent insertions of transposable elements upstream of Zmr1 occurred in 9% of Z. tritici strains from around the world and negatively regulated Zmr1 expression, contributing to variation in melanin accumulation.

CONCLUSIONS

Our studies identified two levels of transcriptional control that regulate the synthesis of melanin. We propose that these regulatory mechanisms evolved to balance the fitness costs associated with melanin production against its positive contribution to survival in stressful environments.

Comparative genomics of Botryosphaeria dothidea and B. kuwatsukai, causal agents of apple ring rot, reveals both species expansion of pathogenicity-related genes and variations in virulence gene content during speciation

Ring rot, one of the most destructive diseases of apple worldwide, is caused primarily by Botryosphaeria dothidea and B. kuwatsukai. Here, we sequenced the genomes of B. dothidea strain PG45 (44.3 Mb with 5.12 % repeat rate) and B. kuwatsukai epitype strain PG2 (48.0 Mb with 13.02 % repeat rate), and conducted a comparative analysis of these two genomes, as well as other sequenced fungal genomes, in order to understand speciation and distinctive patterns of evolution of pathogenicity-related genes. Pair-wise genome alignments revealed that the two species are highly syntenic (96.74 % average sequence identity). Both species encode a significant number of pathogenicity-related genes, e.g. carbohydrate active enzymes (CAZYs), plant cell wall degrading enzymes (PCWDEs), secondary metabolites (SMs) biosynthetic enzymes, cytochrome P450 enzymes (CYPs), and secreted peptidases, in comparison to all additional sequenced fungal species involved in various life-styles. The number of pathogenicity-related genes in B. dothidea and B. kuwatsukai is higher than other genomes of Botryosphaeriaceae pathogens (Macrophomina phaseolina and Neofusicoccum parvum), suggesting a secondary round of Botryosphaeria-lineage expansion in the family. There were, however, also significant differences in the genomes of the two Botryosphaeria species. Botryosphaeria kuwatsukai, which infects only apple and pear, apparently lost a set of SMs genes, CAZYs and PCWDEs, possibly as a result of host specialization. Botryosphaeria kuwatsukai contained significantly more transposable elements and higher value of repeat induced point (RIP) index than B. dothidea. Our results will be instrumental in understanding how both phytopathogens interact with their plant hosts and in designing efficient strategies for disease control and molecular breeding to help ensure global apple production and food security.

SNARE-Encoding Genes VdSec22 and VdSso1 Mediate Protein Secretion Required for Full Virulence in Verticillium dahliae

Proteins that mediate cellular and subcellular membrane fusion are key factors in vesicular trafficking in all eukaryotic cells, including the secretion and transport of plant pathogen virulence factors. In this study, we identified vesicle-fusion components that included 22 soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs), four Sec1/Munc18 (SM) family proteins, and 10 Rab GTPases encoded in the genome of the vascular wilt pathogen Verticillium dahliae Vd991. Targeted deletion of two SNARE-encoding genes in V. dahliae, VdSec22 and VdSso1, significantly reduced virulence of both mutants on cotton, relative to the wild-type Vd991 strain. Comparative analyses of the secreted protein content (exoproteome) revealed that many enzymes involved in carbohydrate hydrolysis were regulated by VdSec22 or VdSso1. Consistent with a role of these enzymes in plant cell-wall degradation, pectin, cellulose, and xylan utilization were reduced in the VdSec22 or VdSso1 mutant strains along with a loss of exoproteome cytotoxic activity on cotton leaves. Comparisons with a pathogenicity-related exoproteome revealed that several known virulence factors were not regulated by VdSec22 or VdSso1, but some of the proteins regulated by VdSec22 or VdSso1 displayed different characteristics, including the lack of a typical signal peptide, suggesting that V. dahliae employs more than one secretory route to transport proteins to extracellular sites during infection.

Meiosis Leads to Pervasive Copy-Number Variation and Distorted Inheritance of Accessory Chromosomes of the Wheat Pathogen Zymoseptoria tritici

Meiosis is one of the most conserved molecular processes in eukaryotes. The fidelity of pairing and segregation of homologous chromosomes has a major impact on the proper transmission of genetic information. Aberrant chromosomal transmission can have major phenotypic consequences, yet the mechanisms are poorly understood. Fungi are excellent models to investigate processes of chromosomal transmission, because many species have highly polymorphic genomes that include accessory chromosomes. Inheritance of accessory chromosomes is often unstable and chromosomal losses have little impact on fitness. We analyzed chromosomal inheritance in 477 progeny coming from two crosses of the fungal wheat pathogen Zymoseptoria tritici. For this, we developed a high-throughput screening method based on restriction site-associated DNA sequencing that generated dense coverage of genetic markers along each chromosome. We identified rare instances of chromosomal duplications (disomy) in core chromosomes. Accessory chromosomes showed high overall frequencies of disomy. Chromosomal rearrangements were found exclusively on accessory chromosomes and were more frequent than disomy. Accessory chromosomes present in only one of the parents in an analyzed cross were inherited at significantly higher rates than the expected 1:1 segregation ratio. Both the chromosome and the parental background had significant impacts on the rates of disomy, losses, rearrangements, and distorted inheritance. We found that chromosomes with higher sequence similarity and lower repeat content were inherited more faithfully. The large number of rearranged progeny chromosomes identified in this species will enable detailed analyses of the mechanisms underlying chromosomal rearrangement.

Multiple independent origins for a subtelomeric locus associated with growth rate in Fusarium circinatum

Fusarium is a diverse assemblage that includes a large number of species of considerable medical and agricultural importance. Not surprisingly, whole genome sequences for many Fusarium species have been published or are in the process of being determined, the availability of which is invaluable for deciphering the genetic basis of key phenotypic traits. Here we investigated the distribution, genic composition, and evolutionary history of a locus potentially determining growth rate in the pitch canker pathogen F. circinatum. We found that the genomic region underlying this locus is highly conserved amongst F. circinatum and its close relatives, except for the presence of a 12 000 base pair insertion in all of the examined isolates of F. circinatum. This insertion encodes for five genes and our phylogenetic analyses revealed that each was most likely acquired through horizontal gene transfer from polyphyletic origins. Our data further showed that this region is located in a region low in G+C content and enriched for repetitive sequences and transposable elements, which is situated near the telomere of Chromosome 3 of F. circinatum. As have been shown for other fungi, these findings thus suggest that the emergence of the unique 12 000 bp region in F. circinatum is linked to the dynamic evolutionary processes associated with subtelomeres that, in turn, have been implicated in the ecological adaptation of fungal pathogens.

Origin and Spread of Spliceosomal Introns: Insights from the Fungal Clade Zymoseptoria

Spliceosomal introns are a key feature of eukaryote genome architecture and have been proposed to originate from selfish group II introns from an endosymbiotic bacterium, that is, the ancestor of mitochondria. However, the mechanisms underlying the wide spread of spliceosomal introns across eukaryotic genomes have been obscure. In this study, we characterize the dynamic evolution of spliceosomal introns in the fungal genus Zymoseptoria at different evolutionary scales, that is, within a genome, among conspecific strains within species, and between different species. Within the genome, spliceosomal introns can proliferate in unrelated genes and intergenic regions. Among conspecific strains, spliceosomal introns undergo rapid turnover (gains and losses) and frequent sequence exchange between geographically distinct strains. Furthermore, spliceosomal introns could undergo introgression between distinct species, which can further promote intron invasion and proliferation. The dynamic invasion and proliferation processes of spliceosomal introns resemble the life cycles of mobile selfish (group I/II) introns, and these intron movements, at least in part, account for the dramatic processes of intron gain and intron loss during eukaryotic evolution.

Proximal Phenotyping and Machine Learning Methods to Identify Septoria Tritici Blotch Disease Symptoms in Wheat

Phenotyping with proximal sensors allow high-precision measurements of plant traits both in the controlled conditions and in the field. In this work, using machine learning, an integrated analysis was done from the data obtained from spectroradiometer, infrared thermometer, and chlorophyll fluorescence measurements to identify most predictive proxy measurements for studying Septoria tritici blotch (STB) disease of wheat. The random forest (RF) models for chlorosis and necrosis identified photosystem II quantum yield (QY) and vegetative indices (VIs) associated with the biochemical composition of leaves as the top predictive variables for identifying disease symptoms. The RF model for chlorosis was validated with a validation set (R2: 0.80) and in an independent test set (R2: 0.55). Based on the results, it can be concluded that the proxy measurements for photosystem II, chlorophyll content, carotenoid, and anthocyanin levels and leaf surface temperature can be successfully used to detect STB. Further validation of these results in the field will enable application of these predictive variables for detection of STB in the field.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

Fungal Genomes and Insights into the Evolution of the Kingdom

The kingdom Fungi comprises species that inhabit nearly all ecosystems. Fungi exist as both free-living and symbiotic unicellular and multicellular organisms with diverse morphologies. The genomes of fungi encode genes that enable them to thrive in diverse environments, invade plant and animal cells, and participate in nutrient cycling in terrestrial and aquatic ecosystems. The continuously expanding databases of fungal genome sequences have been generated by individual and large-scale efforts such as Génolevures, Broad Institute's Fungal Genome Initiative, and the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). These efforts have produced a catalog of fungal genes and genomic organization. The genomic datasets can be utilized to better understand how fungi have adapted to their lifestyles and ecological niches. Large datasets of fungal genomic and transcriptomic data have enabled the use of novel methodologies and improved the study of fungal evolution from a molecular sequence perspective. Combined with microscopes, petri dishes, and woodland forays, genome sequencing supports bioinformatics and comparative genomics approaches as important tools in the study of the biology and evolution of fungi.

Comparative Transcriptomics Reveals How Wheat Responds to Infection by Zymoseptoria tritici

The fungus Zymoseptoria tritici causes septoria tritici blotch (STB) on wheat, an important disease globally and the most damaging wheat disease in Europe. Despite the global significance of STB, the molecular basis of wheat defense against Z. tritici is poorly understood. Here, we use a comparative transcriptomic study to investigate how wheat responds to infection by four distinct strains of Z. tritici. We examined the response of wheat across the entire infection cycle, identifying both shared responses to the four strains and strain-specific responses. We found that the early asymptomatic phase is characterized by strong upregulation of genes encoding receptor-like kinases and pathogenesis-related proteins, indicating the onset of a defense response. We also identified genes that were differentially expressed among the four fungal strains, including genes related to defense. Genes involved in senescence were induced during both the asymptomatic phase and at late stages of infection, suggesting manipulation of senescence processes by both the plant and the pathogen. Our findings illustrate the need, when identifying important genes affecting disease resistance in plants, to include multiple pathogen strains.

Verticillium dahliae transcription factor VdFTF1 regulates the expression of multiple secreted virulence factors and is required for full virulence in cotton

Fungal transcription factors (TFs) implicated in the regulation of virulence gene expression have been identified in a number of plant pathogens. In Verticillium dahliae, despite its agricultural importance, few regulators of transcription have been characterized. In this study, a T-DNA insertion mutant with significantly reduced virulence towards cotton was identified. The T-DNA was traced to VdFTF1, a gene encoding a TF containing a Fungal_trans domain. Transient expression in onion epidermal cells indicated that VdFTF1 is localized to the nucleus. The VdFTF1-deletion strains displayed normal vegetative growth, mycelial pigmentation and conidial morphology, but exhibited significantly reduced virulence on cotton, suggesting that VdFTF1 is required exclusively for pathogenesis. Comparisons of global transcription patterns of wild-type and VdFTF1-deletion strains indicated that VdFTF1 affected the expression of 802 genes, 233 of which were associated with catalytic processes. These genes encoded 69 potentially secreted proteins, 43 of which contained a carbohydrate enzyme domain known to participate in pathogenesis during infection of cotton. Targeted gene deletion of one VdFTF1-regulated gene resulted in significantly impaired vascular colonization, as measured by quantitative polymerase chain reaction, as well as aggressiveness and symptom severity in cotton. In conclusion, VdFTF1, which encodes a TF containing a Fungal_trans domain, regulates the gene expression of plant cell wall degradation enzymes in V. dahliae, which are required for full virulence on cotton.

Distinct Trajectories of Massive Recent Gene Gains and Losses in Populations of a Microbial Eukaryotic Pathogen

Differences in gene content are a significant source of variability within species and have an impact on phenotypic traits. However, little is known about the mechanisms responsible for the most recent gene gains and losses. We screened the genomes of 123 worldwide isolates of the major pathogen of wheat Zymoseptoria tritici for robust evidence of gene copy number variation. Based on orthology relationships in three closely related fungi, we identified 599 gene gains and 1,024 gene losses that have not yet reached fixation within the focal species. Our analyses of gene gains and losses segregating in populations showed that gene copy number variation arose preferentially in subtelomeres and in proximity to transposable elements. Recently lost genes were enriched in virulence factors and secondary metabolite gene clusters. In contrast, recently gained genes encoded mostly secreted protein lacking a conserved domain. We analyzed the frequency spectrum at loci segregating a gene presence–absence polymorphism in four worldwide populations. Recent gene losses showed a significant excess in low-frequency variants compared with genome-wide single nucleotide polymorphism, which is indicative of strong negative selection against gene losses. Recent gene gains were either under weak negative selection or neutral. We found evidence for strong divergent selection among populations at individual loci segregating a gene presence–absence polymorphism. Hence, gene gains and losses likely contributed to local adaptation. Our study shows that microbial eukaryotes harbor extensive copy number variation within populations and that functional differences among recently gained and lost genes led to distinct evolutionary trajectories.

Fine-Scale Recombination Maps of Fungal Plant Pathogens Reveal Dynamic Recombination Landscapes and Intragenic Hotspots

Meiotic recombination is an important driver of evolution. Variability in the intensity of recombination across chromosomes can affect sequence composition, nucleotide variation, and rates of adaptation. In many organisms, recombination events are concentrated within short segments termed recombination hotspots. The variation in recombination rate and positions of recombination hotspot can be studied using population genomics data and statistical methods. In this study, we conducted population genomics analyses to address the evolution of recombination in two closely related fungal plant pathogens: the prominent wheat pathogen Zymoseptoria tritici and a sister species infecting wild grasses Z. ardabiliae We specifically addressed whether recombination landscapes, including hotspot positions, are conserved in the two recently diverged species and if recombination contributes to rapid evolution of pathogenicity traits. We conducted a detailed simulation analysis to assess the performance of methods of recombination rate estimation based on patterns of linkage disequilibrium, in particular in the context of high nucleotide diversity. Our analyses reveal overall high recombination rates, a lack of suppressed recombination in centromeres, and significantly lower recombination rates on chromosomes that are known to be accessory. The comparison of the recombination landscapes of the two species reveals a strong correlation of recombination rate at the megabase scale, but little correlation at smaller scales. The recombination landscapes in both pathogen species are dominated by frequent recombination hotspots across the genome including coding regions, suggesting a strong impact of recombination on gene evolution. A significant but small fraction of these hotspots colocalize between the two species, suggesting that hotspot dynamics contribute to the overall pattern of fast evolving recombination in these species.

Genome-Enhanced Detection and Identification (GEDI) of plant pathogens

Plant diseases caused by fungi and Oomycetes represent worldwide threats to crops and forest ecosystems. Effective prevention and appropriate management of emerging diseases rely on rapid detection and identification of the causal pathogens. The increase in genomic resources makes it possible to generate novel genome-enhanced DNA detection assays that can exploit whole genomes to discover candidate genes for pathogen detection. A pipeline was developed to identify genome regions that discriminate taxa or groups of taxa and can be converted into PCR assays. The modular pipeline is comprised of four components: (1) selection and genome sequencing of phylogenetically related taxa, (2) identification of clusters of orthologous genes, (3) elimination of false positives by filtering, and (4) assay design. This pipeline was applied to some of the most important plant pathogens across three broad taxonomic groups: Phytophthoras (Stramenopiles, Oomycota), Dothideomycetes (Fungi, Ascomycota) and Pucciniales (Fungi, Basidiomycota). Comparison of 73 fungal and Oomycete genomes led the discovery of 5,939 gene clusters that were unique to the targeted taxa and an additional 535 that were common at higher taxonomic levels. Approximately 28% of the 299 tested were converted into qPCR assays that met our set of specificity criteria. This work demonstrates that a genome-wide approach can efficiently identify multiple taxon-specific genome regions that can be converted into highly specific PCR assays. The possibility to easily obtain multiple alternative regions to design highly specific qPCR assays should be of great help in tackling challenging cases for which higher taxon-resolution is needed.

A Near-Complete Haplotype-Phased Genome of the Dikaryotic Wheat Stripe Rust Fungus Puccinia striiformis f. sp. tritici Reveals High Interhaplotype Diversity

A long-standing biological question is how evolution has shaped the genomic architecture of dikaryotic fungi. To answer this, high-quality genomic resources that enable haplotype comparisons are essential. Short-read genome assemblies for dikaryotic fungi are highly fragmented and lack haplotype-specific information due to the high heterozygosity and repeat content of these genomes. Here, we present a diploid-aware assembly of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici based on long reads using the FALCON-Unzip assembler. Transcriptome sequencing data sets were used to infer high-quality gene models and identify virulence genes involved in plant infection referred to as effectors. This represents the most complete Puccinia striiformis f. sp. tritici genome assembly to date (83 Mb, 156 contigs, N50 of 1.5 Mb) and provides phased haplotype information for over 92% of the genome. Comparisons of the phase blocks revealed high interhaplotype diversity of over 6%. More than 25% of all genes lack a clear allelic counterpart. When we investigated genome features that potentially promote the rapid evolution of virulence, we found that candidate effector genes are spatially associated with conserved genes commonly found in basidiomycetes. Yet, candidate effectors that lack an allelic counterpart are more distant from conserved genes than allelic candidate effectors and are less likely to be evolutionarily conserved within the P. striiformis species complex and Pucciniales In summary, this haplotype-phased assembly enabled us to discover novel genome features of a dikaryotic plant-pathogenic fungus previously hidden in collapsed and fragmented genome assemblies.IMPORTANCE Current representations of eukaryotic microbial genomes are haploid, hiding the genomic diversity intrinsic to diploid and polyploid life forms. This hidden diversity contributes to the organism's evolutionary potential and ability to adapt to stress conditions. Yet, it is challenging to provide haplotype-specific information at a whole-genome level. Here, we take advantage of long-read DNA sequencing technology and a tailored-assembly algorithm to disentangle the two haploid genomes of a dikaryotic pathogenic wheat rust fungus. The two genomes display high levels of nucleotide and structural variations, which lead to allelic variation and the presence of genes lacking allelic counterparts. Nonallelic candidate effector genes, which likely encode important pathogenicity factors, display distinct genome localization patterns and are less likely to be evolutionary conserved than those which are present as allelic pairs. This genomic diversity may promote rapid host adaptation and/or be related to the age of the sequenced isolate since last meiosis.

Stress and sexual reproduction affect the dynamics of the wheat pathogen effector AvrStb6 and strobilurin resistance

Host resistance and fungicide treatments are cornerstones of plant-disease control. Here, we show that these treatments allow sex and modulate parenthood in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that the Z. tritici-wheat interaction complies with the gene-for-gene model by identifying the effector AvrStb6, which is recognized by the wheat resistance protein Stb6. Recognition triggers host resistance, thus implying removal of avirulent strains from pathogen populations. However, Z. tritici crosses on wheat show that sex occurs even with an avirulent parent, and avirulence alleles are thereby retained in subsequent populations. Crossing fungicide-sensitive and fungicide-resistant isolates under fungicide pressure results in a rapid increase in resistance-allele frequency. Isolates under selection always act as male donors, and thus disease control modulates parenthood. Modeling these observations for agricultural and natural environments reveals extended durability of host resistance and rapid emergence of fungicide resistance. Therefore, fungal sex has major implications for disease control.

Genome sequencing and comparative genomics reveal the potential pathogenic mechanism of Cercospora sojina Hara on soybean

Frogeye leaf spot, caused by Cercospora sojina Hara, is a common disease of soybean in most soybean-growing countries of the world. In this study, we report a high-quality genome sequence of C. sojina by Single Molecule Real-Time sequencing method. The 40.8-Mb genome encodes 11,655 predicated genes, and 8,474 genes are revealed by RNA sequencing. Cercospora sojina genome contains large numbers of gene clusters that are involved in synthesis of secondary metabolites, including mycotoxins and pigments. However, much less carbohydrate-binding module protein encoding genes are identified in C. sojina genome, when compared with other phytopathogenic fungi. Bioinformatics analysis reveals that C. sojina harbours about 752 secreted proteins, and 233 of them are effectors. During early infection, the genes for metabolite biosynthesis and effectors are significantly enriched, suggesting that they may play essential roles in pathogenicity. We further identify 13 effectors that can inhibit BAX-induced cell death. Taken together, our results provide insights into the infection mechanisms of C. sojina on soybean.

Widespread signatures of selection for secreted peptidases in a fungal plant pathogen

BACKGROUND

Fungal plant pathogens secrete a large arsenal of hydrolytic enzymes during the course of infection, including peptidases. Secreted peptidases have been extensively studied for their role as effectors. In this study, we combined transcriptomics, comparative genomics and evolutionary analyses to investigate all 39 secreted peptidases in the fungal wheat pathogen Zymoseptoria tritici and its close relatives Z. pseudotritici and Z. ardabiliae.

RESULTS

RNA-seq data revealed that a majority of the secreted peptidases displayed differential transcription during the course of Z. tritici infection, indicative of specialization for different stages in the life cycle. Evolutionary analyses detected widespread evidence of adaptive evolution acting on at least 28 of the peptidases. A few peptidases displayed lineage-specific rates of molecular evolution, suggesting altered selection pressure in Z. tritici following host specialization on domesticated wheat. The peptidases belonging to MEROPS families A1 and G1 emerged as a particularly interesting group that may play key roles in host-pathogen co-evolution, host adaptation and pathogenicity. Sister genes in the A1 and G1 families showed accelerated substitution rates after gene duplications.

CONCLUSIONS

These results suggest widespread evolution of secreted peptidases leading to novel gene functions, consistent with predicted models of "escape from adaptive conflict" and "neo-functionalization". Our analyses identified candidate genes worthy of functional analyses that may encode effector functions, for example by suppressing plant defenses during the biotrophic phase of infection.

Pangenome analyses of the wheat pathogen Zymoseptoria tritici reveal the structural basis of a highly plastic eukaryotic genome

Background

Structural variation contributes substantially to polymorphism within species. Chromosomal rearrangements that impact genes can lead to functional variation among individuals and influence the expression of phenotypic traits. Genomes of fungal pathogens show substantial chromosomal polymorphism that can drive virulence evolution on host plants. Assessing the adaptive significance of structural variation is challenging, because most studies rely on inferences based on a single reference genome sequence.

Results

We constructed and analyzed the pangenome of Zymoseptoria tritici, a major pathogen of wheat that evolved host specialization by chromosomal rearrangements and gene deletions. We used single-molecule real-time sequencing and high-density genetic maps to assemble multiple genomes. We annotated the gene space based on transcriptomics data that covered the infection life cycle of each strain. Based on a total of five telomere-to-telomere genomes, we constructed a pangenome for the species and identified a core set of 9149 genes. However, an additional 6600 genes were exclusive to a subset of the isolates. The substantial accessory genome encoded on average fewer expressed genes but a larger fraction of the candidate effector genes that may interact with the host during infection. We expanded our analyses of the pangenome to a worldwide collection of 123 isolates of the same species. We confirmed that accessory genes were indeed more likely to show deletion polymorphisms and loss-of-function mutations compared to core genes.

Conclusions

The pangenome construction of a highly polymorphic eukaryotic pathogen showed that a single reference genome significantly underestimates the gene space of a species. The substantial accessory genome provides a cradle for adaptive evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0457-4) contains supplementary material, which is available to authorized users.

Comparative genomics reveals cotton-specific virulence factors in flexible genomic regions in Verticillium dahliae and evidence of horizontal gene transfer from Fusarium

Verticillium dahliae isolates are most virulent on the host from which they were originally isolated. Mechanisms underlying these dominant host adaptations are currently unknown. We sequenced the genome of V. dahliae Vd991, which is highly virulent on its original host, cotton, and performed comparisons with the reference genomes of JR2 (from tomato) and VdLs.17 (from lettuce). Pathogenicity-related factor prediction, orthology and multigene family classification, transcriptome analyses, phylogenetic analyses, and pathogenicity experiments were performed. The Vd991 genome harbored several exclusive, lineage-specific (LS) genes within LS regions (LSRs). Deletion mutants of the seven genes within one LSR (G-LSR2) in Vd991 were less virulent only on cotton. Integration of G-LSR2 genes individually into JR2 and VdLs.17 resulted in significantly enhanced virulence on cotton but did not affect virulence on tomato or lettuce. Transcription levels of the seven LS genes in Vd991 were higher during the early stages of cotton infection, as compared with other hosts. Phylogenetic analyses suggested that G-LSR2 was acquired from Fusarium oxysporum f. sp. vasinfectum through horizontal gene transfer. Our results provide evidence that horizontal gene transfer from Fusarium to Vd991 contributed significantly to its adaptation to cotton and may represent a significant mechanism in the evolution of an asexual plant pathogen.

A high-resolution genetic map of the cereal crown rot pathogen Fusarium pseudograminearum provides a near-complete genome assembly

Fusarium pseudograminearum is an important pathogen of wheat and barley, particularly in semi-arid environments. Previous genome assemblies for this organism were based entirely on short read data and are highly fragmented. In this work, a genetic map of F. pseudograminearum has been constructed for the first time based on a mapping population of 178 individuals. The genetic map, together with long read scaffolding of a short read-based genome assembly, was used to give a near-complete assembly of the four F. pseudograminearum chromosomes. Large regions of synteny between F. pseudograminearum and F. graminearum, the related pathogen that is the primary causal agent of cereal head blight disease, were previously proposed in the core conserved genome, but the construction of a genetic map to order and orient contigs is critical to the validation of synteny and the placing of species-specific regions. Indeed, our comparative analyses of the genomes of these two related pathogens suggest that rearrangements in the F. pseudograminearum genome have occurred in the chromosome ends. One of these rearrangements includes the transposition of an entire gene cluster involved in the detoxification of the benzoxazolinone (BOA) class of plant phytoalexins. This work provides an important genomic and genetic resource for F. pseudograminearum, which is less well characterized than F. graminearum. In addition, this study provides new insights into a better understanding of the sexual reproduction process in F. pseudograminearum, which informs us of the potential of this pathogen to evolve.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

Cut-and-Paste Transposons in Fungi with Diverse Lifestyles

Transposable elements (TEs) shape genomes via recombination and transposition, lead to chromosomal rearrangements, create new gene neighborhoods, and alter gene expression. They play key roles in adaptation either to symbiosis in Amanita genus or to pathogenicity in Pyrenophora tritici-repentis. Despite growing evidence of their importance, the abundance and distribution of mobile elements replicating in a "cut-and-paste" fashion is barely described so far. In order to improve our knowledge on this old and ubiquitous class of transposable elements, 1,730 fungal genomes were scanned using both de novo and homology-based approaches. DNA TEs have been identified across the whole data set and display uneven distribution from both DNA TE classification and fungal taxonomy perspectives. DNA TE content correlates with genome size, which confirms that many transposon families proliferate simultaneously. In contrast, it is independent from intron density, average gene distance and GC content. TE count is associated with species' lifestyle and tends to be elevated in plant symbionts and decreased in animal parasites. Lastly, we found that fungi with both RIP and RNAi systems have more total DNA TE sequences but less elements retaining a functional transposase, what reflects stringent control over transposition.

Forward Genetics Approach Reveals Host Genotype-Dependent Importance of Accessory Chromosomes in the Fungal Wheat Pathogen Zymoseptoria tritici

The fungal wheat pathogen Zymoseptoria tritici possesses a large complement of accessory chromosomes showing presence/absence polymorphism among isolates. These chromosomes encode hundreds of genes; however, their functional role and why the chromosomes have been maintained over long evolutionary times are so far not known. In this study, we addressed the functional relevance of eight accessory chromosomes in reference isolate IPO323. We induced chromosome losses by inhibiting the β-tubulin assembly during mitosis using carbendazim and generated several independent isogenic strains, each lacking one of the accessory chromosomes. We confirmed chromosome losses by electrophoretic karyotyping and whole-genome sequencing. To assess the importance of the individual chromosomes during host infection, we performed in planta assays comparing disease development results in wild-type and chromosome mutant strains. Loss of the accessory chromosomes 14, 16, 18, 19, and 21 resulted in increased virulence on wheat cultivar Runal but not on cultivars Obelisk, Titlis, and Riband. Moreover, some accessory chromosomes affected the switch from biotrophy to necrotrophy as strains lacking accessory chromosomes 14, 18, 19, and 21 showed a significantly earlier onset of necrosis than the wild type on the Runal cultivar. In general, we observed that the timing of the lifestyle switch affects the fitness of Z. tritici. Taking the results together, this study was the first to use a forward-genetics approach to demonstrate a cultivar-dependent functional relevance of the accessory chromosomes of Z. tritici during host infection.

Zymoseptoria tritici is a prominent fungal pathogen of wheat of worldwide distribution. This fungus shows a remarkable genome organization, with a large number of chromosomes that are present in only some isolates and therefore considered to be “accessory” chromosomes. To date, the function of these accessory chromosomes in Z. tritici has been unknown, although their maintenance in the species over evolutionary times suggests a functional relevance. Here we deleted whole accessory chromosomes to test the effect of these chromosomes on host specificity and virulence of the fungus. We show for the first time that some accessory chromosomes of Z. tritici affect the fitness of the fungus during host infection in a cultivar-dependent manner. These results show that the accessory chromosomes encode host-specific virulence determinants having a negative effect on fitness. Understanding the population dynamic of the accessory chromosomes and the molecular interaction of pathogen and plant traits is crucial to improve wheat-breeding strategies.

Plasticity of the MFS1 Promoter Leads to Multidrug Resistance in the Wheat Pathogen Zymoseptoria tritici

The ascomycete Zymoseptoria tritici is the causal agent of Septoria leaf blotch on wheat. Disease control relies mainly on resistant wheat cultivars and on fungicide applications. The fungus displays a high potential to circumvent both methods. Resistance against all unisite fungicides has been observed over decades. A different type of resistance has emerged among wild populations with multidrug-resistant (MDR) strains. Active fungicide efflux through overexpression of the major facilitator gene MFS1 explains this emerging resistance mechanism. Applying a bulk-progeny sequencing approach, we identified in this study a 519-bp long terminal repeat (LTR) insert in the MFS1 promoter, a relic of a retrotransposon cosegregating with the MDR phenotype. Through gene replacement, we show the insert as a mutation responsible for MFS1 overexpression and the MDR phenotype. Besides this type I insert, we found two different types of promoter inserts in more recent MDR strains. Type I and type II inserts harbor potential transcription factor binding sites, but not the type III insert. Interestingly, all three inserts correspond to repeated elements present at different genomic locations in either IPO323 or other Z. tritici strains. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici and which contribute to its adaptive potential. IMPORTANCE Disease control through fungicides remains an important means to protect crops from fungal diseases and to secure the harvest. Plant-pathogenic fungi, especially Zymoseptoria tritici, have developed resistance against most currently used active ingredients, reducing or abolishing their efficacy. While target site modification is the most common resistance mechanism against single modes of action, active efflux of multiple drugs is an emerging phenomenon in fungal populations reducing additionally fungicides' efficacy in multidrug-resistant strains. We have investigated the mutations responsible for increased drug efflux in Z. tritici field strains. Our study reveals that three different insertions of repeated elements in the same promoter lead to multidrug resistance in Z. tritici. The target gene encodes the membrane transporter MFS1 responsible for drug efflux, with the promoter inserts inducing its overexpression. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici.

The ZtVf1 transcription factor regulates development and virulence in the foliar wheat pathogen Zymoseptoria tritici

The dimorphic fungal pathogen, Zymoseptoria tritici undergoes discrete developmental changes to complete its life cycle on wheat. Molecular mechanisms underlying morphogenesis during infection process of Z. tritici are poorly understood. In this study, we have investigated the role of ZtVf1 gene encoding a transcription factor belonging to C₂-H₂ subfamily. In planta assays revealed that ZtVf1 is required for virulence. Reduced necrotic lesions and low pycnidia density within the lesions resulted in significantly reduced virulence of ZtVf1 mutants. Cytological analysis showed that the impaired virulence of ZtVf1 mutants attributed to reduced penetration and colonization along with hampered pycnidia differentiation. In vitro phenotyping showed that ZtVf1 deletion affects hyphal branching and biomass production suggesting that the reduced tissue colonization by the ZtVf1 mutant might be due to lower hyphal branching and less fungal biomass production. In addition, the majority of infected substomatal cavities by the ZtVf1 mutant filled with compacted mycelia mat that did not differentiate to mature pycnidia indicating that the impaired melanization negatively affected pycnidia formation and maturation. The ZtVf1 might target multiple genes belonging to different cellular processes whose identification is of eminent interest to increase our understanding of this pathosystem. Overall, the data provided in this study indicates that attenuated pathogenicity of ZtVf1 mutant is due to involvement of this gene in the regulation of both early and late stages of infection.

A conserved fungal glycosyltransferase facilitates pathogenesis of plants by enabling hyphal growth on solid surfaces

Pathogenic fungi must extend filamentous hyphae across solid surfaces to cause diseases of plants. However, the full inventory of genes which support this is incomplete and many may be currently concealed due to their essentiality for the hyphal growth form. During a random T-DNA mutagenesis screen performed on the pleomorphic wheat (Triticum aestivum) pathogen Zymoseptoria tritici, we acquired a mutant unable to extend hyphae specifically when on solid surfaces. In contrast "yeast-like" growth, and all other growth forms, were unaffected. The inability to extend surface hyphae resulted in a complete loss of virulence on plants. The affected gene encoded a predicted type 2 glycosyltransferase (ZtGT2). Analysis of >800 genomes from taxonomically diverse fungi highlighted a generally widespread, but discontinuous, distribution of ZtGT2 orthologues, and a complete absence of any similar proteins in non-filamentous ascomycete yeasts. Deletion mutants of the ZtGT2 orthologue in the taxonomically un-related fungus Fusarium graminearum were also severely impaired in hyphal growth and non-pathogenic on wheat ears. ZtGT2 expression increased during filamentous growth and electron microscopy on deletion mutants (ΔZtGT2) suggested the protein functions to maintain the outermost surface of the fungal cell wall. Despite this, adhesion to leaf surfaces was unaffected in ΔZtGT2 mutants and global RNAseq-based gene expression profiling highlighted that surface-sensing and protein secretion was also largely unaffected. However, ΔZtGT2 mutants constitutively overexpressed several transmembrane and secreted proteins, including an important LysM-domain chitin-binding virulence effector, Zt3LysM. ZtGT2 likely functions in the synthesis of a currently unknown, potentially minor but widespread, extracellular or outer cell wall polysaccharide which plays a key role in facilitating many interactions between plants and fungi by enabling hyphal growth on solid matrices.

Overview of craft brewing specificities and potentially associated microbiota

The brewing process differs slightly in craft breweries as compared to industrial breweries, as there are fewer control points. This affects the microbiota of the final product. Beer contains several antimicrobial properties that protect it from pathogens, such as low pH, low oxygen and high carbon dioxide content, and the addition of hops. However, these hurdles have limited power controlling spoilage organisms. Contamination by these organisms can originate in the raw materials, persist in the environment, and be introduced by using flavoring ingredients later in the process. Spoilage is a prominent issue in brewing, and can cause quality degradation resulting in consumer rejection and product waste. For example, lactic acid bacteria are predominately associated with producing a ropy texture and haze, along with producing diacetyl which gives the beer butter flavor notes. Other microorganisms may not affect flavor or aroma, but can retard fermentation by consuming nutrients needed by fermentation yeast. Quality control in craft breweries today relies on culturing methods to detect specific spoilage organisms. Using media can be beneficial for detecting the most common beer spoilers, such as Lactobacillus and Pediococci. However, these methods are time consuming with long incubation periods. Molecular methods such as community profiling or high throughput sequencing are better used for identifying entire populations of beer. These methods allow for detection, differentiation, and identification of taxa.

Proteomic Characterization of Armillaria mellea Reveals Oxidative Stress Response Mechanisms and Altered Secondary Metabolism Profiles

Armillaria mellea is a major plant pathogen. Yet, the strategies the organism uses to infect susceptible species, degrade lignocellulose and other plant material and protect itself against plant defences and its own glycodegradative arsenal are largely unknown. Here, we use a combination of gel and MS-based proteomics to profile A. mellea under conditions of oxidative stress and changes in growth matrix. 2-DE and LC-MS/MS were used to investigate the response of A. mellea to H₂O₂ and menadione/FeCl₃ exposure, respectively. Several proteins were detected with altered abundance in response to H₂O₂, but not menadione/FeCl₃ (i.e., valosin-containing protein), indicating distinct responses to these different forms of oxidative stress. One protein, cobalamin-independent methionine synthase, demonstrated a common response in both conditions, which may be a marker for a more general stress response mechanism. Further changes to the A. mellea proteome were investigated using MS-based proteomics, which identified changes to putative secondary metabolism (SM) enzymes upon growth in agar compared to liquid cultures. Metabolomic analyses revealed distinct profiles, highlighting the effect of growth matrix on SM production. This establishes robust methods by which to utilize comparative proteomics to characterize this important phytopathogen.

Genome sequence of the ectophytic fungus Ramichloridium luteum reveals unique evolutionary adaptations to plant surface niche

BACKGROUND

Ectophytic fungi occupy the waxy plant surface, an extreme environment characterized by prolonged desiccation, nutrient limitation, and exposure to solar radiation. The nature of mechanisms that facilitate adaptation to this environment remains unclear. In this study, we sequenced the complete genome of an ectophytic fungus, Ramichloridium luteum, which colonizes the surface of apple fruit, and carried out comparative genomic and transcriptome analysis.

RESULTS

The R. luteum genome was 28.18 Mb and encoded 9466 genes containing 1.85% repetitive elements. Compared with cell-penetrating pathogens, genes encoding plant cell wall degrading enzymes (PCWDEs), PTH11-like G protein-coupled receptors (GPCRs) and effectors were drastically reduced. In contrast, genes encoding cutinases and secretory lipases were strikingly expanded, and four of nine secretory lipases were probably acquired by horizontal gene transfer from Basidiomycota. Transcriptomic analysis revealed elevated expression of genes involved in cuticle degradation (cutinase, secretory lipase) and stress responses (melanin biosynthesis, aquaporins, lysozymes and HOG pathway).

CONCLUSIONS

Taken together, our results highlight genomic features associated with evolution of surface niche adaptation by the ectophytic fungus R. luteum, namely the contraction of PCWDEs, PTH11-like GPCRs and effectors, and the expansion of cuticle degradation and stress tolerance.

Dose-dependent selection drives lineage replacement during the experimental evolution of SDHI fungicide resistance in Zymoseptoria tritici

Fungicide resistance is a constant threat to agricultural production worldwide. Molecular mechanisms of fungicide resistance have been studied extensively in the wheat pathogen Zymoseptoria tritici. However, less is known about the evolutionary processes driving resistance development. In vitro evolutionary studies give the opportunity to investigate this. Here, we examine the adaptation of Z. tritici to fluxapyroxad, a succinate dehydrogenase (Sdh) inhibitor. Replicate populations of Z. tritici derived from the sensitive isolate IPO323 were exposed to increasing concentrations of fluxapyroxad with or without UV mutagenesis. After ten increases in fungicide concentration, sensitivity had decreased dramatically, with replicate populations showing similar phenotypic trajectories. Sequencing the Sdh subunit B, C, and D encoding genes identified seven mutations associated with resistance to fluxapyroxad. Mutation frequency over time was measured with a pyrosequencing assay, revealing sequential lineage replacement in the UV‐mutagenized populations but not in the untreated populations. Repeating selection from set time‐points with different fungicide concentrations revealed that haplotype replacement of Sdh variants was driven by dose‐dependent selection as fungicide concentration changed, and was not mutation‐limited. These findings suggest that fungicide field applications may select for highly insensitive Sdh variants with higher resistance factors if the fungicide concentration is increased to achieve a better disease control. However, in the absence or presence of lower fungicide concentrations, the spread of these strains might be restricted if the underlying Sdh mutations carry fitness penalties.

A role for random, humidity-dependent epiphytic growth prior to invasion of wheat by Zymoseptoria tritici

Zymoseptoria tritici causes Septoria leaf blotch of wheat. The prevailing paradigm of the Z. tritici-wheat interaction assumes fungal ingress through stomata within 24-48h, followed by days of symptomless infection. This is extrapolated from studies testing the mode of fungal ingress under optimal infection conditions. Here, we explicitly assess the timing of entry, using GFP-tagged Z. tritici. We show that early entry is comparatively rare, and extended epiphytic growth possible. We test the hypotheses that our data diverge from earlier studies due to: i. random ingress of Z. tritici into the leaf, with some early entry events; ii. previous reliance upon fungal stains, combined with poor attachment of Z. tritici to the leaf, leading to increased likelihood of observing internal versus external growth, compared to using GFP; iii. use of exceptionally high humidity to promote entry in previous studies. We combine computer simulation of leaf-surface growth with thousands of in planta observations to demonstrate that while spores germinate rapidly on the leaf, over 95% of fungi remain epiphytic, growing randomly over the leaf for ten days or more. We show that epiphytic fungi are easily detached from leaves by rinsing and that humidity promotes epiphytic growth, increasing infection rates. Together, these results explain why epiphytic growth has been dismissed and early ingress assumed. The prolonged epiphytic phase should inform studies of pathogenicity and virulence mutants, disease control strategies, and interpretation of the observed low in planta growth, metabolic quiescence and evasion of plant defences by Zymoseptoria during symptomless infection.

Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

BACKGROUND

The ascomycete fungus Colletotrichum higginsianum causes anthracnose disease of brassica crops and the model plant Arabidopsis thaliana. Previous versions of the genome sequence were highly fragmented, causing errors in the prediction of protein-coding genes and preventing the analysis of repetitive sequences and genome architecture.

RESULTS

Here, we re-sequenced the genome using single-molecule real-time (SMRT) sequencing technology and, in combination with optical map data, this provided a gapless assembly of all twelve chromosomes except for the ribosomal DNA repeat cluster on chromosome 7. The more accurate gene annotation made possible by this new assembly revealed a large repertoire of secondary metabolism (SM) key genes (89) and putative biosynthetic pathways (77 SM gene clusters). The two mini-chromosomes differed from the ten core chromosomes in being repeat- and AT-rich and gene-poor but were significantly enriched with genes encoding putative secreted effector proteins. Transposable elements (TEs) were found to occupy 7% of the genome by length. Certain TE families showed a statistically significant association with effector genes and SM cluster genes and were transcriptionally active at particular stages of fungal development. All 24 subtelomeres were found to contain one of three highly-conserved repeat elements which, by providing sites for homologous recombination, were probably instrumental in four segmental duplications.

CONCLUSION

The gapless genome of C. higginsianum provides access to repeat-rich regions that were previously poorly assembled, notably the mini-chromosomes and subtelomeres, and allowed prediction of the complete SM gene repertoire. It also provides insights into the potential role of TEs in gene and genome evolution and host adaptation in this asexual pathogen.