Cause and Effectors: Whole-Genome Comparisons Reveal Shared but Rapidly Evolving Effector Sets among Host-Specific Plant-Castrating Fungi

Characterization of Microbotryum lychnidis-dioicae Secreted Effector Proteins, Their Potential Host Targets, and Localization in a Heterologous Host Plant

Microbotryum lychnidis-dioicae is an obligate fungal species colonizing the plant host, Silene latifolia. The fungus synthesizes and secretes effector proteins into the plant host during infection to manipulate the host for completion of the fungal lifecycle. The goal of this study was to continue functional characterization of such M. lychnidis-dioicae effectors. Here, we identified three putative effectors and their putative host-plant target proteins. MVLG_02245 is highly upregulated in M. lychnidis-dioicae during infection; yeast two-hybrid analysis suggests it targets a tubulin α-1 chain protein ortholog in the host, Silene latifolia. A potential plant protein interacting with MVLG_06175 was identified as CASP-like protein 2C1 (CASPL2C1), which facilitates the polymerization of the Casparian strip at the endodermal cells. Proteins interacting with MVLG_05122 were identified as CSN5a or 5b, involved in protein turnover. Fluorescently labelled MVLG_06175 and MVLG_05122 were expressed in the heterologous plant, Arabidopsis thaliana. MVLG_06175 formed clustered granules at the tips of trichomes on leaves and in root caps, while MVLG_05122 formed a band structure at the base of leaf trichomes. Plants expressing MVLG_05122 alone were more resistant to infection with Fusarium oxysporum. These results indicate that the fungus might affect the formation of the Casparian strip in the roots and the development of trichomes during infection as well as alter plant innate immunity.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

Using machine learning to predict protein-protein interactions between a zombie ant fungus and its carpenter ant host

Parasitic fungi produce proteins that modulate virulence, alter host physiology, and trigger host responses. These proteins, classified as a type of "effector," often act via protein-protein interactions (PPIs). The fungal parasite Ophiocordyceps camponoti-floridani (zombie ant fungus) manipulates Camponotus floridanus (carpenter ant) behavior to promote transmission. The most striking aspect of this behavioral change is a summit disease phenotype where infected hosts ascend and attach to an elevated position. Plausibly, interspecific PPIs drive aspects of Ophiocordyceps infection and host manipulation. Machine learning PPI predictions offer high-throughput methods to produce mechanistic hypotheses on how this behavioral manipulation occurs. Using D-SCRIPT to predict host-parasite PPIs, we found ca. 6000 interactions involving 2083 host proteins and 129 parasite proteins, which are encoded by genes upregulated during manipulated behavior. We identified multiple overrepresentations of functional annotations among these proteins. The strongest signals in the host highlighted neuromodulatory G-protein coupled receptors and oxidation-reduction processes. We also detected Camponotus structural and gene-regulatory proteins. In the parasite, we found enrichment of Ophiocordyceps proteases and frequent involvement of novel small secreted proteins with unknown functions. From these results, we provide new hypotheses on potential parasite effectors and host targets underlying zombie ant behavioral manipulation.

Comparative transcriptome analysis of juniper branches infected by Gymnosporangium spp. highlights their different infection strategies associated with cytokinins

BACKGROUND

Gymnosporangium asiaticum and G. yamadae can share Juniperus chinensis as the telial host, but the symptoms are completely different. The infection of G. yamadae causes the enlargement of the phloem and cortex of young branches as a gall, but not for G. asiaticum, suggesting that different molecular interaction mechanisms exist the two Gymnosporangium species with junipers.

RESULTS

Comparative transcriptome analysis was performed to investigate genes regulation of juniper in responses to the infections of G. asiaticum and G. yamadae at different stages. Functional enrichment analysis showed that genes related to transport, catabolism and transcription pathways were up-regulated, while genes related to energy metabolism and photosynthesis were down-regulated in juniper branch tissues after infection with G. asiaticum and G. yamadae. The transcript profiling of G. yamadae-induced gall tissues revealed that more genes involved in photosynthesis, sugar metabolism, plant hormones and defense-related pathways were up-regulated in the vigorous development stage of gall compared to the initial stage, and were eventually repressed overall. Furthermore, the concentration of cytokinins (CKs) in the galls tissue and the telia of G. yamadae was significantly higher than in healthy branch tissues of juniper. As well, tRNA-isopentenyltransferase (tRNA-IPT) was identified in G. yamadae with highly expression levels during the gall development stages.

CONCLUSIONS

In general, our study provided new insights into the host-specific mechanisms by which G. asiaticum and G. yamadae differentially utilize CKs and specific adaptations on juniper during their co-evolution.

The use of evolutionary analyses to predict functionally relevant traits in filamentous plant pathogens

Identifying traits involved in plant-pathogen interactions is one of the major objectives in molecular plant pathology. Evolutionary analyses may assist in the identification of genes encoding traits that are involved in virulence and local adaptation, including adaptation to agricultural intervention strategies. In the past decades, the number of available genome sequences of fungal plant pathogens has rapidly increased, providing a rich source for the discovery of functionally important genes as well as inference of species histories. Positive selection in the form of diversifying or directional selection leaves particular signatures in genome alignments and can be identified with statistical genetics methods. This review summarises the concepts and approaches used in evolutionary genomics and lists major discoveries related to plant-pathogen adaptative evolution. We underline the significant contribution of evolutionary genomics in discovering virulence-related traits and the study of plant-pathogen ecology and adaptive evolution.

A Bacterial Type Three Secretion-Based Delivery System for Functional Characterization of Sporisorium scitamineum Plant Immune Suppressing Effector Proteins

The facultative biotrophic basidiomycete Sporisorium scitamineum causes smut disease in sugarcane. This study applied an assay to identify S. scitamineum candidate effectors (CEs) with plant immunity suppression activities by delivering them into Nicotiana benthamiana cells via the type-three secretion system of Pseudomonas fluorescens EtHAn. Six CEs were individually cloned into the pEDV6 vector and expressed by P. fluorescens EtHAn for translocation into the plant cells. Three CEs (g1052, g3890, and g5159) could suppress pattern-triggered immunity (PTI) responses with high reproducibility in different coinfiltration experiments with P. syringae pv. tomato DC3000. In addition, three CEs (g1052, g4549, and g5159) were also found to be AvrB-induced suppressors of effector-triggered immunity (ETI), demonstrating for the first time that S. scitamineum can defeat both PTI and ETI responses. A transcriptomic analysis at different stages of infection by the smut fungus of three sugarcane cultivars with contrasting responses to the pathogen revealed that suppressors g1052, g3890, g4549, and g5159 were induced at the early stage of infection. By contrast, the two CEs (g2666 and g6610) that did not exhibit suppression activities expressed only at the late stage of infection. Moreover, genomic structures of the CEs and searches for orthologs in other smut species suggested duplication events and further divergence in CEs evolution of S. scitamineum. Thus, the transient assay applied here demonstrated the potential of pEDV6 and P. fluorescens EtHAn as biological tools for identifying plant immune suppressors from S. scitamineum.

Action Mechanisms of Effectors in Plant-Pathogen Interaction

Plant pathogens are one of the main factors hindering the breeding of cash crops. Pathogens, including oomycetes, fungus, and bacteria, secrete effectors as invasion weapons to successfully invade and propagate in host plants. Here, we review recent advances made in the field of plant-pathogen interaction models and the action mechanisms of phytopathogenic effectors. The review illustrates how effectors from different species use similar and distinct strategies to infect host plants. We classify the main action mechanisms of effectors in plant-pathogen interactions according to the infestation process: targeting physical barriers for disruption, creating conditions conducive to infestation, protecting or masking themselves, interfering with host cell physiological activity, and manipulating plant downstream immune responses. The investigation of the functioning of plant pathogen effectors contributes to improved understanding of the molecular mechanisms of plant-pathogen interactions. This understanding has important theoretical value and is of practical significance in plant pathology and disease resistance genetics and breeding.

Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens

Biotrophic plant pathogenic fungi are widely distributed and are among the most damaging pathogenic organisms of agriculturally important crops responsible for significant losses in quality and yield. However, the pathogenesis of obligate parasitic pathogenic microorganisms is still under investigation because they cannot reproduce and complete their life cycle on an artificial medium. The successful lifestyle of biotrophic fungal pathogens depends on their ability to secrete effector proteins to manipulate or evade plant defense response. By integrating genomics, transcriptomics, and effectoromics, insights into how the adaptation of biotrophic plant fungal pathogens adapt to their host populations can be gained. Efficient tools to decipher the precise molecular mechanisms of rust–plant interactions, and standardized routines in genomics and functional pipelines have been established and will pave the way for comparative studies. Deciphering fungal pathogenesis not only allows us to better understand how fungal pathogens infect host plants but also provides valuable information for plant diseases control, including new strategies to prevent, delay, or inhibit fungal development. Our review provides a comprehensive overview of the efforts that have been made to decipher the effector proteins of biotrophic fungal pathogens and demonstrates how rapidly research in the field of obligate biotrophy has progressed.

Lifestyle Transitions in Fusarioid Fungi are Frequent and Lack Clear Genomic Signatures

The fungal genus Fusarium (Ascomycota) includes well-known plant pathogens that are implicated in diseases worldwide, and many of which have been genome sequenced. The genus also encompasses other diverse lifestyles, including species found ubiquitously as asymptomatic-plant inhabitants (endophytes). Here, we produced structurally annotated genome assemblies for five endophytic Fusarium strains, including the first whole-genome data for Fusarium chuoi. Phylogenomic reconstruction of Fusarium and closely related genera revealed multiple and frequent lifestyle transitions, the major exception being a monophyletic clade of mutualist insect symbionts. Differential codon usage bias and increased codon optimisation separated Fusarium sensu stricto from allied genera. We performed computational prediction of candidate secreted effector proteins (CSEPs) and carbohydrate-active enzymes (CAZymes)—both likely to be involved in the host–fungal interaction—and sought evidence that their frequencies could predict lifestyle. However, phylogenetic distance described gene variance better than lifestyle did. There was no significant difference in CSEP, CAZyme, or gene repertoires between phytopathogenic and endophytic strains, although we did find some evidence that gene copy number variation may be contributing to pathogenicity. Large numbers of accessory CSEPs (i.e., present in more than one taxon but not all) and a comparatively low number of strain-specific CSEPs suggested there is a limited specialisation among plant associated Fusarium species. We also found half of the core genes to be under positive selection and identified specific CSEPs and CAZymes predicted to be positively selected on certain lineages. Our results depict fusarioid fungi as prolific generalists and highlight the difficulty in predicting pathogenic potential in the group.

Using structural variants to understand the ecological and evolutionary dynamics of fungal plant pathogens

Deletions, duplications, insertions, inversions and translocations are commonly referred to as structural variants (SVs). Fungal plant pathogens have compact genomes, facilitating the generation of accurate maps of SVs for these species in recent studies. Structural variants have been found to constitute a significant proportion of the standing genetic variation in fungal plant pathogen populations, potentially leading to the generation of accessory genes, regions or chromosomes enriched in pathogenicity factors. Structural variants are involved in the rapid adaptation and ecological traits of pathogens, including host specialization and mating. Long-read sequencing techniques coupled with theoretical and experimental approaches have considerable potential for elucidating the phenotypic effects of SVs and deciphering the evolutionary and genomic mechanisms underlying the formation of SVs in fungal plant pathogens.

The population genetics of adaptation through copy‐number variation in a fungal plant pathogen

Microbial pathogens can adapt rapidly to changing environments such as the application of pesticides or host resistance. Copy number variations (CNVs) are a major source of adaptive genetic variation for recent adaptation. Here, we analyse how a major fungal pathogen of barley, Rhynchosporium commune, has adapted to the host environment and fungicide applications. We screen the genomes of 125 isolates sampled across a worldwide set of populations and identify a total of 7,879 gene duplications and 116 gene deletions. Most gene duplications result from segmental chromosomal duplications. Although CNVs are generally under negative selection, we find that genes affected by CNVs are enriched in functions related to host exploitation (i.e., effectors and cell-wall-degrading enzymes). We perform genome-wide association studies (GWAS) and identify a large segmental duplication of CYP51A that has contributed to the emergence of azole resistance and a duplication encompassing an effector gene affecting virulence. We show that the adaptive CNVs were probably created by recently active transposable element families. Moreover, we find that specific transposable element families are important drivers of recent gene CNV. Finally, we use a genome-wide single nucleotide polymorphism data set to replicate the GWAS and contrast it with the CNV-focused analysis. Together, our findings show how extensive segmental duplications create the raw material for recent adaptation in global populations of a fungal pathogen.

Genome-wide analysis of Claviceps paspali: insights into the secretome of the main species causing ergot disease in Paspalum spp

BACKGROUND

The phytopatogen Claviceps paspali is the causal agent of Ergot disease in Paspalum spp., which includes highly productive forage grasses such as P. dilatatum. This disease impacts dairy and beef production by affecting seed quality and producing mycotoxins that can affect performance in feeding animals. The molecular basis of pathogenicity of C. paspali remains unknown, which makes it more difficult to find solutions for this problem. Secreted proteins are related to fungi virulence and can manipulate plant immunity acting on different subcellular localizations. Therefore, identifying and characterizing secreted proteins in phytopathogenic fungi will provide a better understanding of how they overcome host defense and cause disease. The aim of this work is to analyze the whole genome sequences of three C. paspali isolates to obtain a comparative genome characterization based on possible secreted proteins and pathogenicity factors present in their genome. In planta RNA-seq analysis at an early stage of the interaction of C. paspali with P. dilatatum stigmas was also conducted in order to determine possible secreted proteins expressed in the infection process.

RESULTS

C. paspali isolates had compact genomes and secretome which accounted for 4.6-4.9% of the predicted proteomes. More than 50% of the predicted secretome had no homology to known proteins. RNA-Seq revealed that three protein-coding genes predicted as secreted have mayor expression changes during 1 dpi vs 4 dpi. Also, three of the first 10 highly expressed genes in both time points were predicted as effector-like. CAZyme-like proteins were found in the predicted secretome and the most abundant family could be associated to pectine degradation. Based on this, pectine could be a main component affected by the cell wall degrading enzymes of C. paspali.

CONCLUSIONS

Based on predictions from DNA sequence and RNA-seq, unique probable secreted proteins and probable pathogenicity factors were identified in C. paspali isolates. This information opens new avenues in the study of the biology of this fungus and how it modulates the interaction with its host. Knowledge of the diversity of the secretome and putative pathogenicity genes should facilitate future research in disease management of Claviceps spp.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

Cross-species analysis between the maize smut fungi Ustilago maydis and Sporisorium reilianum highlights the role of transcriptional change of effector orthologs for virulence and disease

The constitution and regulation of effector repertoires shape host-microbe interactions. Ustilago maydis and Sporisorium reilianum are two closely related smut fungi, which both infect maize but cause distinct disease symptoms. Understanding how effector orthologs are regulated in these two pathogens can therefore provide insights into the evolution of different infection strategies. We tracked the infection progress of U. maydis and S. reilianum in maize leaves and used two distinct infection stages for cross-species RNA-sequencing analyses. We identified 207 of 335 one-to-one effector orthologs as differentially regulated during host colonization, which might reflect the distinct disease development strategies. Using CRISPR-Cas9-mediated gene conversion, we identified two differentially expressed effector orthologs with conserved function between two pathogens. Thus, differential expression of functionally conserved genes might contribute to species-specific adaptation and symptom development. Interestingly, another differentially expressed orthogroup (UMAG_05318/Sr10075) showed divergent protein function, providing a possible case for neofunctionalization. Collectively, we demonstrated that the diversification of effector genes in related pathogens can be caused both by alteration on the transcriptional level and through functional diversification of the encoded effector proteins.

Comparative genomics of two inbred lines of the potato cyst nematode Globodera rostochiensis reveals disparate effector family-specific diversification patterns

BACKGROUND

Potato cyst nematodes belong to the most harmful pathogens in potato, and durable management of these parasites largely depends on host-plant resistances. These resistances are pathotype specific. The current Globodera rostochiensis pathotype scheme that defines five pathotypes (Ro1 - Ro5) is both fundamentally and practically of limited value. Hence, resistant potato varieties are used worldwide in a poorly informed manner.

RESULTS

We generated two novel reference genomes of G. rostochiensis inbred lines derived from a Ro1 and a Ro5 population. These genome sequences comprise 173 and 189 scaffolds respectively, marking a ≈ 24-fold reduction in fragmentation as compared to the current reference genome. We provide copy number variations for 19 effector families. Four dorsal gland effector families were investigated in more detail. SPRYSECs, known to be implicated in plant defence suppression, constitute by far the most diversified family studied herein with 60 and 99 variants in Ro1 and Ro5 distributed over 18 and 26 scaffolds. In contrast, CLEs, effectors involved in feeding site induction, show strong physical clustering. The 10 and 16 variants cluster on respectively 2 and 1 scaffolds. Given that pathotypes are defined by their effectoromes, we pinpoint the disparate nature of the contributing effector families in terms of sequence diversification and loss and gain of variants.

CONCLUSIONS

Two novel reference genomes allow for nearly complete inventories of effector diversification and physical organisation within and between pathotypes. Combined with insights we provide on effector family-specific diversification patterns, this constitutes a basis for an effectorome-based virulence scheme for this notorious pathogen.

Positive Selection of Transcription Factors Is a Prominent Feature of the Evolution of a Plant Pathogenic Genus Originating in the Miocene

Tests based on the d_N/d_S statistic are used to identify positive selection of nonsynonymous polymorphisms. Using these tests on alignments of all orthologs from related species can provide insights into which gene categories have been most frequently positively selected. However, longer alignments have more power to detect positive selection, creating a detection bias that could create misleading results from functional enrichment tests. Most studies of positive selection in plant pathogens focus on genes with specific virulence functions, with little emphasis on broader molecular processes. Furthermore, no studies in plant pathogens have accounted for detection bias due to alignment length when performing functional enrichment tests. To address these research gaps, we analyze 12 genomes of the phytopathogenic fungal genus Botrytis, including two sequenced in this study. To establish a temporal context, we estimated fossil-calibrated divergence times for the genus. We find that Botrytis likely originated 16–18 Ma in the Miocene and underwent continuous radiation ending in the Pliocene. An untargeted scan of Botrytis single-copy orthologs for positive selection with three different statistical tests uncovered evidence for positive selection among proteases, signaling proteins, CAZymes, and secreted proteins. There was also a strong overrepresentation of transcription factors among positively selected genes. This overrepresentation was still apparent after two complementary controls for detection bias due to sequence length. Positively selected sites were depleted within DNA-binding domains, suggesting changes in transcriptional responses to internal and external cues or protein–protein interactions have undergone positive selection more frequently than changes in promoter fidelity.

Rapid sequence evolution driven by transposable elements at a virulence locus in a fungal wheat pathogen

BACKGROUND

Plant pathogens cause substantial crop losses in agriculture production and threaten food security. Plants evolved the ability to recognize virulence factors and pathogens have repeatedly escaped recognition due rapid evolutionary change at pathogen virulence loci (i.e. effector genes). The presence of transposable elements (TEs) in close physical proximity of effector genes can have important consequences for gene regulation and sequence evolution. Species-wide investigations of effector gene loci remain rare hindering our ability to predict pathogen evolvability.

RESULTS

Here, we performed genome-wide association studies (GWAS) on a highly polymorphic mapping population of 120 isolates of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe. We identified a major locus underlying significant variation in reproductive success of the pathogen and damage caused on the wheat cultivar Claro. The most strongly associated locus is intergenic and flanked by genes encoding a predicted effector and a serine-type endopeptidase. The center of the locus contained a highly dynamic region consisting of multiple families of TEs. Based on a large global collection of assembled genomes, we show that the virulence locus has undergone substantial recent sequence evolution. Large insertion and deletion events generated length variation between the flanking genes by a factor of seven (5-35 kb). The locus showed also strong signatures of genomic defenses against TEs (i.e. RIP) contributing to the rapid diversification of the locus.

CONCLUSIONS

In conjunction, our work highlights the power of combining GWAS and population-scale genome analyses to investigate major effect loci in pathogens.

Effectors with Different Gears: Divergence of Ustilago maydis Effector Genes Is Associated with Their Temporal Expression Pattern during Plant Infection

Plant pathogens secrete a variety of effector proteins that enable host colonization but are also typical pathogen detection targets for the host immune system. Consequently, effector genes encounter high selection pressures, which typically makes them fast evolving. The corn smut pathogen Ustilago maydis has an effector gene repertoire with a dynamic expression across the different disease stages. We determined the amino acid divergence of U. maydis effector candidates with Sporisorium reilianum orthologs, a close relative of U. maydis. Intriguingly, there are two distinct groups of effector candidates, ones with a respective conserved and diverged protein evolution. Conservatively evolving effector genes especially have their peak expression during the (pre-)penetration stages of the disease cycle. In contrast, expression of divergently evolving effector genes generally peaks during fungal proliferation within the host. To test if this interspecific effector diversity corresponds to intraspecific diversity, we sampled and sequenced a diverse collection of U. maydis strains from the most important maize breeding and production regions in China. Effector candidates with a diverged interspecific evolution had more intraspecific amino acid variation than candidates with a conserved evolution. In conclusion, we highlight diversity in evolution within the U. maydis effector repertoire with dynamically and conservatively evolving members.

Fungal effectors, the double edge sword of phytopathogens

Phyto-pathogenic fungi can cause huge damage to crop production. During millions of years of coexistence, fungi have evolved diverse life-style to obtain nutrients from the host and to colonize upon them. They deploy various proteinaceous as well as non-proteinaceous secreted molecules commonly referred as effectors to sabotage host machinery during the infection process. The effectors are important virulence determinants of pathogenic fungi and play important role in successful pathogenesis, predominantly by avoiding host-surveillance system. However, besides being important for pathogenesis, the fungal effectors end-up being recognized by the resistant cultivars of the host, which mount a strong immune response to ward-off pathogens. Various recent studies involving different pathosystem have revealed the virulence/avirulence functions of fungal effectors and their involvement in governing the outcome of host-pathogen interactions. However, the effectors and their cognate resistance gene in the host remain elusive for several economically important fungal pathogens. In this review, using examples from some of the biotrophic, hemi-biotrophic and necrotrophic pathogens, we elaborate the double-edged functions of fungal effectors. We emphasize that knowledge of effector functions can be helpful in effective management of fungal diseases in crop plants.

Meiotic recombination in the offspring of Microbotryum hybrids and its impact on pathogenicity

BACKGROUND

Hybridization is a central mechanism in evolution, producing new species or introducing important genetic variation into existing species. In plant-pathogenic fungi, adaptation and specialization to exploit a host species are key determinants of evolutionary success. Here, we performed experimental crosses between the two pathogenic Microbotryum species, M. lychnidis-dioicae and M. silenes-acaulis that are specialized to different hosts. The resulting offspring were analyzed on phenotypic and genomic levels to describe genomic characteristics of hybrid offspring and genetic factors likely involved in host-specialization.

RESULTS

Genomic analyses of interspecific fungal hybrids revealed that individuals were most viable if the majority of loci were inherited from one species. Interestingly, species-specific loci were strictly controlled by the species' origin of the mating type locus. Moreover we detected signs of crossing over and chromosome duplications in the genomes of the analyzed hybrids. In Microbotryum, mitochondrial DNA was found to be uniparentally inherited from the a2 mating type. Genome comparison revealed that most gene families are shared and the majority of genes are conserved between the two species, indicating very similar biological features, including infection and pathogenicity processes. Moreover, we detected 211 candidate genes that were retained under host-driven selection of backcrossed lines. These genes and might therefore either play a crucial role in host specialization or be linked to genes that are essential for specialization.

CONCLUSION

The combination of genome analyses with experimental selection and hybridization is a promising way to investigate host-pathogen interactions. This study manifests genetic factors of host specialization that are required for successful biotrophic infection of the post-zygotic stage, but also demonstrates the strong influence of intra-genomic conflicts or instabilities on the viability of hybrids in the haploid host-independent stage.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.