Breakage-fusion-bridge cycles and large insertions contribute to the rapid evolution of accessory chromosomes in a fungal pathogen

Identification of a novel genetic locus conferring virulence in the wheat tan spot pathogen Pyrenophora tritici-repentis

The ascomycete Pyrenophora tritici-repentis (Ptr) is the causal agent of tan spot, a common and economically important wheat disease worldwide. Three necrotrophic effectors (NEs), known as Ptr ToxA, Ptr ToxB, and Ptr ToxC, have been identified from the fungal pathogen as major virulence factors. The race 2 isolate 86-124 which produces Ptr ToxA is capable of causing disease on wheat lines that is insensitive to Ptr ToxA, suggesting the use of additional NEs during the infection. To identify new NE gene(s) from 86-124, we developed a biparental fungal population from a cross between this isolate and the race 5 isolate DW5 using genetically modified heterothallic strains. The fungal population was genotyped with SNP and SSR markers as well as the ToxA gene, the mating type genes, and six ToxB loci. Each progeny was phenotyped onto the hard red spring wheat line CDC-Osler, which is insensitive to both Ptr ToxA and Ptr ToxB, but is highly susceptible to 86-124. The constructed genetic map consisted of 11 linkage groups that corresponded to the 11 chromosomes (chr) of the Ptr reference genome. ToxA and mating type genes mapped to the expected positions. Five of the six ToxB copies were tightly linked, residing at the distal end of chr 11, while the sixth copy was localized to the distal end of chr 5. Composite interval mapping revealed a major QTL on the distal end of chr 2 conferring virulence toward CDC-Osler by 86-124. This locus was designated as VirOsler1. Genomic sequence alignment at the locus showed a region of approximately 900 kb at the end of chr 2 absent in DW5. The identification of VirOsler1 locus provides clear evidence that the wheat tan spot pathogen uses additional virulence factors that interact with unidentified host factors for disease susceptibility.

Candida albicans isolates contain frequent heterozygous structural variants and transposable elements within genes and centromeres

The human fungal pathogen Candida albicans poses a significant burden on global health, causing high rates of mortality and antifungal drug resistance. C. albicans is a heterozygous diploid organism that reproduces asexually. Structural variants (SVs) are an important source of genomic rearrangement, particularly in species that lack sexual recombination. To comprehensively investigate SVs across clinical isolates of C. albicans, we conducted long-read sequencing and genome-wide SV analysis in three distantly related clinical isolates. Our work includes a new, comprehensive analysis of transposable element (TE) composition, location, and diversity. SVs and TEs are frequently close to coding sequences and many SVs are heterozygous, suggesting that SVs might impact gene and allele-specific expression. Most SVs are uniquely present in only one clinical isolate, indicating that SVs represent a significant source of intraspecies genetic variation. We identify multiple, distinct SVs at the centromeres of Chromosome 4 and Chromosome 5, including inversions and transposon polymorphisms. These two chromosomes are often aneuploid in drug-resistant clinical isolates and can form isochromosome structures with breakpoints near the centromere. Further screening of 100 clinical isolates confirms the widespread presence of centromeric SVs in C. albicans, often appearing in a heterozygous state, indicating that SVs are contributing to centromere evolution in C. albicans Together, these findings highlight that SVs and TEs are common across diverse clinical isolates of C. albicans and that the centromeres of this organism are important sites of genome rearrangement.

From Lesions to Lessons: Two Decades of Filamentous Plant Pathogen Genomics

Many filamentous microorganisms, such as fungi and oomycetes, have evolved the ability to colonize plants and cause devastating crop diseases. Coevolutionary conflicts with their hosts have shaped the genomes of these plant pathogens. Over the past 20 years, genomics and genomics-enabled technologies have revealed remarkable diversity in genome size, architecture, and gene regulatory mechanisms. Technical and conceptual advances continue to provide novel insights into evolutionary dynamics, diversification of distinct genomic compartments, and facilitated molecular disease diagnostics. In this review, we discuss how genomics has advanced our understanding of genome organization and plant-pathogen coevolution and provide a perspective on future developments in the field. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

The role of gene copy number variation in antimicrobial resistance in human fungal pathogens

Faced with the burden of increasing resistance to antifungals in many fungal pathogens and the constant emergence of new drug-resistant strains, it is essential to assess the importance of various resistance mechanisms. Fungi have relatively plastic genomes and can tolerate genomic copy number variation (CNV) caused by aneuploidy and gene amplification or deletion. In many cases, these genomic changes lead to adaptation to stressful conditions, including those caused by antifungal drugs. Here, we specifically examine the contribution of CNVs to antifungal resistance. We undertook a thorough literature search, collecting reports of antifungal resistance caused by a CNV, and classifying the examples of CNV-conferred resistance into four main mechanisms. We find that in human fungal pathogens, there is little evidence that gene copy number plays a major role in the emergence of antifungal resistance compared to other types of mutations. We discuss why we might be underestimating their importance and new approaches being used to study them.

Repeat-induced point mutations driving Parastagonospora nodorum genomic diversity are balanced by selection against non-synonymous mutations

Parastagonospora nodorum is necrotrophic fungal pathogen of wheat with significant genomic resources. Population-level pangenome data for 173 isolates, of which 156 were from Western Australia (WA) and 17 were international, were examined for overall genomic diversity and effector gene content. A heterothallic core population occurred across all regions of WA, with asexually-reproducing clonal clusters in dryer northern regions. High potential for SNP diversity in the form of repeat-induced point mutation (RIP)-like transitions, was observed across the genome, suggesting widespread 'RIP-leakage' from transposon-rich repetitive sequences into non-repetitive regions. The strong potential for RIP-like mutations was balanced by negative selection against non-synonymous SNPs, that was observed within protein-coding regions. Protein isoform profiles of known effector loci (SnToxA, SnTox1, SnTox3, SnTox267, and SnTox5) indicated low-levels of non-synonymous and high-levels of silent RIP-like mutations. Effector predictions identified 186 candidate secreted predicted effector proteins (CSEPs), 69 of which had functional annotations and included confirmed effectors. Pangenome-based effector isoform profiles across WA were distinct from global isolates and were conserved relative to population structure, and may enable new approaches for monitoring crop disease pathotypes.

Using recurrent neural networks to detect supernumerary chromosomes in fungal strains causing blast diseases

The genomes of the fungus Magnaporthe oryzae that causes blast diseases on diverse grass species, including major crops, have indispensable core-chromosomes and may contain supernumerary chromosomes, also known as mini-chromosomes. These mini-chromosomes are speculated to provide effector gene mobility, and may transfer between strains. To understand the biology of mini-chromosomes, it is valuable to be able to detect whether a M. oryzae strain possesses a mini-chromosome. Here, we applied recurrent neural network models for classifying DNA sequences as arising from core- or mini-chromosomes. The models were trained with sequences from available core- and mini-chromosome assemblies, and then used to predict the presence of mini-chromosomes in a global collection of M. oryzae isolates using short-read DNA sequences. The model predicted that mini-chromosomes were prevalent in M. oryzae isolates. Interestingly, at least one mini-chromosome was present in all recent wheat isolates, but no mini-chromosomes were found in early isolates collected before 1991, indicating a preferential selection for strains carrying mini-chromosomes in recent years. The model was also used to identify assembled contigs derived from mini-chromosomes. In summary, our study has developed a reliable method for categorizing DNA sequences and showcases an application of recurrent neural networks in predictive genomics.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

Retracing the evolution of Pneumocystis species, with a focus on the human pathogen Pneumocystis jirovecii

SUMMARYEvery human being is presumed to be infected by the fungus Pneumocystis jirovecii at least once in his or her lifetime. This fungus belongs to a large group of species that appear to exclusively infect mammals, with P. jirovecii being the only one known to cause disease in humans. The mystery of P. jirovecii origin and speciation is just beginning to unravel. Here, we provide a review of the major steps of P. jirovecii evolution. The Pneumocystis genus likely originated from soil or plant-associated organisms during the period of Cretaceous ~165 million years ago and successfully shifted to mammals. The transition coincided with a substantial loss of genes, many of which are related to the synthesis of nutrients that can be scavenged from hosts or cell wall components that could be targeted by the mammalian immune system. Following the transition, the Pneumocystis genus cospeciated with mammals. Each species specialized at infecting its own host. Host specialization is presumably built at least partially upon surface glycoproteins, whose protogene was acquired prior to the genus formation. P. jirovecii appeared at ~65 million years ago, overlapping with the emergence of the first primates. P. jirovecii and its sister species P. macacae, which infects macaques nowadays, may have had overlapping host ranges in the distant past. Clues from molecular clocks suggest that P. jirovecii did not cospeciate with humans. Molecular evidence suggests that Pneumocystis speciation involved chromosomal rearrangements and the mounting of genetic barriers that inhibit gene flow among species.

Segmental duplications drive the evolution of accessory regions in a major crop pathogen

Many pathogens evolved compartmentalized genomes with conserved core and variable accessory regions (ARs) that carry effector genes mediating virulence. The fungal plant pathogen Fusarium oxysporum has such ARs, often spanning entire chromosomes. The presence of specific ARs influences the host range, and horizontal transfer of ARs can modify the pathogenicity of the receiving strain. However, how these ARs evolve in strains that infect the same host remains largely unknown. We defined the pan-genome of 69 diverse F. oxysporum strains that cause Fusarium wilt of banana, a significant constraint to global banana production, and analyzed the diversity and evolution of the ARs. Accessory regions in F. oxysporum strains infecting the same banana cultivar are highly diverse, and we could not identify any shared genomic regions and in planta-induced effectors. We demonstrate that segmental duplications drive the evolution of ARs. Furthermore, we show that recent segmental duplications specifically in accessory chromosomes cause the expansion of ARs in F. oxysporum. Taken together, we conclude that extensive recent duplications drive the evolution of ARs in F. oxysporum, which contribute to the evolution of virulence.

Frequent horizontal chromosome transfer between asexual fungal insect pathogens

Entire chromosomes are typically only transmitted vertically from one generation to the next. The horizontal transfer of such chromosomes has long been considered improbable, yet gained recent support in several pathogenic fungi where it may affect the fitness or host specificity. To date, it is unknown how these transfers occur, how common they are, and whether they can occur between different species. In this study, we show multiple independent instances of horizontal transfers of the same accessory chromosome between two distinct strains of the asexual entomopathogenic fungus Metarhizium robertsii during experimental co-infection of its insect host, the Argentine ant. Notably, only the one chromosome-but no other-was transferred from the donor to the recipient strain. The recipient strain, now harboring the accessory chromosome, exhibited a competitive advantage under certain host conditions. By phylogenetic analysis, we further demonstrate that the same accessory chromosome was horizontally transferred in a natural environment between M. robertsii and another congeneric insect pathogen, Metarhizium guizhouense. Hence, horizontal chromosome transfer is not limited to the observed frequent events within species during experimental infections but also occurs naturally across species. The accessory chromosome that was transferred contains genes that may be involved in its preferential horizontal transfer or support its establishment. These genes encode putative histones and histone-modifying enzymes, as well as putative virulence factors. Our study reveals that both intra- and interspecies horizontal transfer of entire chromosomes is more frequent than previously assumed, likely representing a not uncommon mechanism for gene exchange.

Rapid mini-chromosome divergence among fungal isolates causing wheat blast outbreaks in Bangladesh and Zambia

The fungal pathogen, Magnaporthe oryzae Triticum pathotype, causing wheat blast disease was first identified in South America and recently spread across continents to South Asia and Africa. Here, we studied the genetic relationship among isolates found on the three continents. Magnaporthe oryzae strains closely related to a South American field isolate B71 were found to have caused the wheat blast outbreaks in South Asia and Africa. Genomic variation among isolates from the three continents was examined using an improved B71 reference genome and whole-genome sequences. We found strong evidence to support that the outbreaks in Bangladesh and Zambia were caused by the introductions of genetically separated isolates, although they were all close to B71 and, therefore, collectively referred to as the B71 branch. In addition, B71 branch strains carried at least one supernumerary mini-chromosome. Genome assembly of a Zambian strain revealed that its mini-chromosome was similar to the B71 mini-chromosome but with a high level of structural variation. Our findings show that while core genomes of the multiple introductions are highly similar, the mini-chromosomes have undergone marked diversification. The maintenance of the mini-chromosome and rapid genomic changes suggest the mini-chromosomes may serve important virulence or niche adaptation roles under diverse environmental conditions.

Shared infection strategy of a fungal pathogen across diverse lineages of land plants, the Fusarium example

Plants engage with a wide variety of microorganisms either in parasitic or mutualistic relationships, which have helped them to adapt to terrestrial ecosystems. Microbial interactions have driven plant evolution and led to the emergence of complex interaction outcomes via suppression of host defenses by evolving pathogens. The evolution of plant-microbe interactions is shaped by conserved host and pathogen gene modules and fast-paced lineage-specific adaptability which determines the interaction outcome. Recent findings from different microbes ranging from bacteria, oomycetes, and fungi suggest recurrent concepts in establishing interactions with evolutionarily distant plant hosts, but also clade-specific adaptation that ultimately contributes to pathogenicity. Here, we revisit some of the latest features that illustrate shared colonization strategies of the fungal pathogen Fusarium oxysporum on distant plant lineages and lineage-specific adaptability of mini-chromosomal units encoding effectors, for shaping host-specific pathogenicity in angiosperms.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Quantitative trait locus mapping of osmotic stress response in the fungal wheat pathogen Zymoseptoria tritici

Osmotic stress is a ubiquitous and potent stress for all living organisms, but few studies have investigated the genetic basis of salt tolerance in filamentous fungi. The main aim of this study was to identify regions of the genome associated with tolerance to potassium chloride (KCl) in the wheat pathogen Zymoseptoria tritici. A secondary aim was to identify candidate genes affecting salt tolerance within the most promising chromosomal regions. We achieved these aims with a quantitative trait locus (QTL) mapping study using offspring from 2 crosses grown in vitro in the presence or absence of osmotic stress imposed by 0.75 M KCl. We identified significant QTL for most of the traits in both crosses. Several QTLs overlapped with QTL identified in earlier studies for other traits, and some QTL explained trait variation in both the control and salt stress environments. A significant QTL on chromosome 3 explained variation in colony radius at 8-day postinoculation (dpi) in the KCl environment as well as colony radius KCl tolerance at 8 dpi. The QTL peak had a high logarithm of the odds ratio (LOD) and encompassed an interval containing only 36 genes. Six of these genes present promising candidates for functional analyses. A gene ontology (GO) enrichment analysis of QTL unique to the KCl environment found evidence for the enrichment of functions involved in osmotic stress responses.

Virulent strains of Zymoseptoria tritici suppress the host immune response and facilitate the success of avirulent strains in mixed infections

Plants interact with a plethora of pathogenic microorganisms in nature. Pathogen-plant interaction experiments focus mainly on single-strain infections, typically ignoring the complexity of multi-strain infections even though mixed infections are common and critical for the infection outcome. The wheat pathogen Zymoseptoria tritici forms highly diverse fungal populations in which several pathogen strains often colonize the same leaf. Despite the importance of mixed infections, the mechanisms governing interactions between a mixture of pathogen strains within a plant host remain largely unexplored. Here we demonstrate that avirulent pathogen strains benefit from being in mixed infections with virulent strains. We show that virulent strains suppress the wheat immune response, allowing avirulent strains to colonize the apoplast and to reproduce. Our experiments indicate that virulent strains in mixed infections can suppress the plant immune system, probably facilitating the persistence of avirulent pathogen strains in fields planted with resistant host plants.

Disentangling the complex gene interaction networks between rice and the blast fungus identifies a new pathogen effector

Studies focused solely on single organisms can fail to identify the networks underlying host-pathogen gene-for-gene interactions. Here, we integrate genetic analyses of rice (Oryza sativa, host) and rice blast fungus (Magnaporthe oryzae, pathogen) and uncover a new pathogen recognition specificity of the rice nucleotide-binding domain and leucine-rich repeat protein (NLR) immune receptor Pik, which mediates resistance to M. oryzae expressing the avirulence effector gene AVR-Pik. Rice Piks-1, encoded by an allele of Pik-1, recognizes a previously unidentified effector encoded by the M. oryzae avirulence gene AVR-Mgk1, which is found on a mini-chromosome. AVR-Mgk1 has no sequence similarity to known AVR-Pik effectors and is prone to deletion from the mini-chromosome mediated by repeated Inago2 retrotransposon sequences. AVR-Mgk1 is detected by Piks-1 and by other Pik-1 alleles known to recognize AVR-Pik effectors; recognition is mediated by AVR-Mgk1 binding to the integrated heavy metal-associated (HMA) domain of Piks-1 and other Pik-1 alleles. Our findings highlight how complex gene-for-gene interaction networks can be disentangled by applying forward genetics approaches simultaneously to the host and pathogen. We demonstrate dynamic coevolution between an NLR integrated domain and multiple families of effector proteins.

Heterochromatin in the fungal plant pathogen, Zymoseptoria tritici: Control of transposable elements, genome plasticity and virulence

Heterochromatin is a repressive chromatin state that plays key roles in the functional organisation of eukaryotic genomes. In fungal plant pathogens, effector genes that are required for host colonization tend to be associated with heterochromatic regions of the genome that are enriched with transposable elements. It has been proposed that the heterochromatin environment silences effector genes in the absence of host and dynamic chromatin remodelling facilitates their expression during infection. Here we discuss this model in the context of the key wheat pathogen, Zymoseptoria tritici. We cover progress in understanding the deposition and recognition of heterochromatic histone post translational modifications in Z. tritici and the role that heterochromatin plays in control of genome plasticity and virulence.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Zymoseptoria tritici white-collar complex integrates light, temperature and plant cues to initiate dimorphism and pathogenesis

Transitioning from spores to hyphae is pivotal to host invasion by the plant pathogenic fungus Zymoseptoria tritici. This dimorphic switch can be initiated by high temperature in vitro (~27 °C); however, such a condition may induce cellular heat stress, questioning its relevance to field infections. Here, we study the regulation of the dimorphic switch by temperature and other factors. Climate data from wheat-growing areas indicate that the pathogen sporadically experiences high temperatures such as 27 °C during summer months. However, using a fluorescent dimorphic switch reporter (FDR1) in four wild-type strains, we show that dimorphic switching already initiates at 15-18 °C, and is enhanced by wheat leaf surface compounds. Transcriptomics reveals 1261 genes that are up- or down-regulated in hyphae of all strains. These pan-strain core dimorphism genes (PCDGs) encode known effectors, dimorphism and transcription factors, and light-responsive proteins (velvet factors, opsins, putative blue light receptors). An FDR1-based genetic screen reveals a crucial role for the white-collar complex (WCC) in dimorphism and virulence, mediated by control of PCDG expression. Thus, WCC integrates light with biotic and abiotic cues to orchestrate Z. tritici infection.

Asexual reproductive potential trumps virulence as a predictor of competitive ability in mixed infections

Natural infections frequently involve several co-infecting pathogen strains. These mixed infections can affect the extent of the infection, the transmission success of the pathogen and the eventual epidemic outcome. To date, few studies have investigated how mixed infections affect transmission between hosts. Zymoseptoria tritici is a highly diverse wheat pathogen in which multiple strains often coexist in the same lesion. Here we demonstrate that the most competitive strains often exclude their competitors during serial passages of mixed infections. The outcome of the competition depended on both the host genotype and the genotypes of the competing pathogen strains. Differences in virulence among the strains were not associated with competitive advantages during transmission, while differences in reproductive potential had a strong effect on strain competitive ability. Overall, our findings suggest that host specialization is determined mainly by the ability to successfully transmit offspring to new hosts during mixed infections.

How large and diverse are field populations of fungal plant pathogens? The case of Zymoseptoria tritici

Pathogen populations differ in the amount of genetic diversity they contain. Populations carrying higher genetic diversity are thought to have a greater evolutionary potential than populations carrying less diversity. We used published studies to estimate the range of values associated with two critical components of genetic diversity, the number of unique pathogen genotypes and the number of spores produced during an epidemic, for the septoria tritici blotch pathogen Zymoseptoria tritici. We found that wheat fields experiencing typical levels of infection are likely to carry between 3.1 and 14.0 million pathogen genotypes per hectare and produce at least 2.1-9.9 trillion pycnidiospores per hectare. Given the experimentally derived mutation rate of 3 × 10-10 substitutions per site per cell division, we estimate that between 27 and 126 million pathogen spores carrying adaptive mutations to counteract fungicides and resistant cultivars will be produced per hectare during a growing season. This suggests that most of the adaptive mutations that have been observed in Z. tritici populations can emerge through local selection from standing genetic variation that already exists within each field. The consequences of these findings for disease management strategies are discussed.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

Chromosomal polymorphism of the Ceratocystis fimbriata species complex in Brazil

Ceratocystis fimbriata is an important pathogen that causes wilt in several plant species. Despite the importance of this pathogen, knowledge about its karyotypic polymorphism and genomic architecture is limited. The main objective of this study was to investigate the karyotype of isolates of the C. fimbriata species complex from different host plants and geographical origins in Brazil. First, the identity of the isolates was confirmed conducting multilocus sequence analysis (MLSA) phylogeny using β-tubulin (TUBB), translation elongation factor 1α (TEF-1α) and mating-type (MAT1 and MAT2) gene sequences. To investigate the chromosomal polymorphism, two conditions of pulsed-field gel electrophoresis (PFGE) were used and the karyotypes of the isolates obtained. The retrotransposon-microsatellite amplified polymorphism (REMAP) molecular marker was utilized to assess the genetic variability among isolates. In the MLSA utilizing the concatenated gene sequences, Ceratocystis cacaofunesta and C. fimbriata formed separate clades, but considerable variation among C. fimbriata isolates was observed. Polymorphism in chromosome number and size was found, indicating the existence of genomic differences among isolates and occurrence of chromosomal rearrangements in the species complex. The number of chromosomes varied from seven to nine and the estimated minimum chromosome sizes were estimated to be between 2.7 to 6.0 Mbp. Small polymorphic chromosomes ranging from 1.2 to 1.8 Mbp were observed in all isolates, raising the hypothesis that they could be supernumerary chromosomes. REMAP analysis revealed a high genetic variability and that isolates from the same host tend to group together in a same cluster. Our results bring new insights into the chromosomal diversity and genome organization of the C. fimbriata complex.

Non-Mendelian transmission of accessory chromosomes in fungi

Non-Mendelian transmission has been reported for various genetic elements, ranging from small transposons to entire chromosomes. One prime example of such a transmission pattern are B chromosomes in plants and animals. Accessory chromosomes in fungi are similar to B chromosomes in showing presence/absence polymorphism and being non-essential. How these chromosomes are transmitted during meiosis is however poorly understood—despite their often high impact on the fitness of the host. For several fungal organisms, a non-Mendelian transmission or a mechanistically unique meiotic drive of accessory chromosomes have been reported. In this review, we provide an overview of the possible mechanisms that can cause the non-Mendelian transmission or meiotic drives of fungal accessory chromosomes. We compare processes responsible for the non-Mendelian transmission of accessory chromosomes for different fungal eukaryotes and discuss the structural traits of fungal accessory chromosomes affecting their meiotic transmission. We conclude that research on fungal accessory chromosomes, due to their small size, ease of sequencing, and epigenetic profiling, can complement the study of B chromosomes in deciphering factors that influence and regulate the non-Mendelian transmission of entire chromosomes.

Determinants of endophytic and pathogenic lifestyle in root colonizing fungi

Plant-fungal interactions in the soil crucially impact crop productivity and can range from highly beneficial to detrimental. Accumulating evidence suggests that some root-colonizing fungi shift between endophytic and pathogenic behaviour depending on the host species and that combinations of effector proteins collectively shape the fungal lifestyle on a given plant. In this review we discuss recent advances in our understanding of how fungal infection strategies on roots can lead to contrasting outcomes for the host. We highlight functional similarities and differences in compatibility determinants that control the colonization of specific-cell layers within plant roots, ultimately shaping the continuum between endophytic and pathogenic lifestyle.

A Conserved Hypothetical Gene Is Required but Not Sufficient for Ptr ToxC Production in Pyrenophora tritici-repentis

The fungus Pyrenophora tritici-repentis causes tan spot, an important foliar disease of wheat worldwide. The fungal pathogen produces three necrotrophic effectors, namely Ptr ToxA, Ptr ToxB, and Ptr ToxC to induce necrosis or chlorosis in wheat. Both Ptr ToxA and Ptr ToxB are proteins, and their encoding genes have been cloned. Ptr ToxC was characterized as a low-molecular weight molecule 20 years ago but the one or more genes controlling its production in P. tritici-repentis are unknown. Here, we report the genetic mapping, molecular cloning, and functional analysis of a fungal gene that is required for Ptr ToxC production. The genetic locus controlling the production of Ptr ToxC, termed ToxC, was mapped to a subtelomeric region using segregating biparental populations, genome sequencing, and association analysis. Additional marker analysis further delimited ToxC to a 173-kb region. The predicted genes in the region were examined for presence/absence polymorphism in different races and isolates leading to the identification of a single candidate gene. Functional validation showed that this gene was required but not sufficient for Ptr ToxC production, thus it is designated as ToxC1. ToxC1 encoded a conserved hypothetical protein likely located on the vacuole membrane. The gene was highly expressed during infection, and only one haplotype was identified among 120 isolates sequenced. Our work suggests that Ptr ToxC is not a protein and is likely produced through a cascade of biosynthetic pathway. The identification of ToxC1 is a major step toward revealing the Ptr ToxC biosynthetic pathway and studying its molecular interactions with host factors.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Untargeted Metabolomic Profiling of Fungal Species Populations

This chapter describes protocols for the development of consensus chemical phenotypes or "metabolomes" of fungal populations using ultra-high pressure liquid chromatography coupled to high resolution mass spectrometry (UPLC-HRMS). Isolates are cultured using multiple media conditions to elicit the expression of diverse secondary metabolite biosynthetic gene clusters. The mycelium and spent culture media are extracted using organic solvents and profiled by ultra-high pressure chromatography coupled with a high resolution Thermo Orbitrap XL mass spectrometer with the ability to trap and fragment ions to general MS² spectra. MS data preprocessing is explained and illustrated using the freely available software MZMine 2. Through data processing, binary matrices of mass features can be generated and then combined into a consensus secondary metabolite phenotype of all isolates grown in all media conditions. The production of consensus chemical phenotypes is useful for screening large fungal populations (both inter and intra-species populations) for isolates potentially expressing novel secondary metabolites or analogs of known secondary metabolites.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen

Genome evolution is driven by the activity of transposable elements (TEs). The spread of TEs can have deleterious effects including the destabilization of genome integrity and expansions. However, the precise triggers of genome expansions remain poorly understood because genome size evolution is typically investigated only among deeply divergent lineages. Here, we use a large population genomics dataset of 284 individuals from populations across the globe of Zymoseptoria tritici, a major fungal wheat pathogen. We built a robust map of genome-wide TE insertions and deletions to track a total of 2456 polymorphic loci within the species. We show that purifying selection substantially depressed TE frequencies in most populations, but some rare TEs have recently risen in frequency and likely confer benefits. We found that specific TE families have undergone a substantial genome-wide expansion from the pathogen's center of origin to more recently founded populations. The most dramatic increase in TE insertions occurred between a pair of North American populations collected in the same field at an interval of 25 years. We find that both genome-wide counts of TE insertions and genome size have increased with colonization bottlenecks. Hence, the demographic history likely played a major role in shaping genome evolution within the species. We show that both the activation of specific TEs and relaxed purifying selection underpin this incipient expansion of the genome. Our study establishes a model to recapitulate TE-driven genome evolution over deeper evolutionary timescales.

Genome and transcriptome of a pathogenic yeast, Candida nivariensis

We present a highly contiguous genome and transcriptome of the pathogenic yeast, Candida nivariensis. We sequenced both the DNA and RNA of this species using both the Oxford Nanopore Technologies and Illumina platforms. We assembled the genome into an 11.8 Mb draft composed of 16 contigs with an N50 of 886 Kb, including a circular mitochondrial sequence of 28 Kb. Using direct RNA nanopore sequencing and Illumina cDNA sequencing, we constructed an annotation of our new assembly, supplemented by lifting over genes from Saccharomyces cerevisiae and Candida glabrata.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

A genome-wide genetic linkage map and reference quality genome sequence for a new race in the wheat pathogen Pyrenophora tritici-repentis

Pyrenophora tritici-repentis is an ascomycete fungus that causes tan spot of wheat. The disease has a worldwide distribution and can cause significant yield and quality losses in wheat production. The fungal pathogen is homothallic in nature, which means it can undergo sexual reproduction by selfing to produce pseudothecia on wheat stubble for seasonal survival. Since homothallism precludes the development of bi-parental fungal populations, no genetic linkage map has been developed for P. tritici-repentis for mapping and map-based cloning of fungal virulence genes. In this work, we created two heterothallic strains by deleting one of the mating type genes in each of two parental isolates 86-124 (race 2) and AR CrossB10 (a new race) and developed a bi-parental fungal population between them. The draft genome sequences of the two parental isolates were aligned to the Pt-1C-BFP reference sequence to mine single nucleotide polymorphisms (SNPs). A total of 225 SNP markers were developed for genotyping the entire population. Additionally, 75 simple sequence repeat, and two gene markers were also developed and used in the genotyping. The resulting linkage map consisted of 13 linkage groups spanning 5,075.83 cM in genetic distance. Because the parental isolate AR CrossB10 is a new race and produces Ptr ToxC, it was sequenced using long-read sequencing platforms and de novo assembled into contigs. The majority of the contigs were further anchored into chromosomes with the aid of the linkage maps. The whole genome comparison of AR CrossB10 to the reference genome of M4 revealed a few chromosomal rearrangements. The genetic linkage map and the new AR CrossB10 genome sequence are valuable tools for gene cloning in P. tritici-repentis.

Accessory Chromosome-Acquired Secondary Metabolism in Plant Pathogenic Fungi: The Evolution of Biotrophs Into Host-Specific Pathogens

Accessory chromosomes are strain- or pathotype-specific chromosomes that exist in addition to the core chromosomes of a species and are generally not considered essential to the survival of the organism. Among pathogenic fungal species, accessory chromosomes harbor pathogenicity or virulence factor genes, several of which are known to encode for secondary metabolites that are involved in plant tissue invasion. Accessory chromosomes are of particular interest due to their capacity for horizontal transfer between strains and their dynamic "crosstalk" with core chromosomes. This review focuses exclusively on secondary metabolism (including mycotoxin biosynthesis) associated with accessory chromosomes in filamentous fungi and the role accessory chromosomes play in the evolution of secondary metabolite gene clusters. Untargeted metabolomics profiling in conjunction with genome sequencing provides an effective means of linking secondary metabolite products with their respective biosynthetic gene clusters that reside on accessory chromosomes. While the majority of literature describing accessory chromosome-associated toxin biosynthesis comes from studies of Alternaria pathotypes, the recent discovery of accessory chromosome-associated biosynthetic genes in Fusarium species offer fresh insights into the evolution of biosynthetic enzymes such as non-ribosomal peptide synthetases (NRPSs), polyketide synthases (PKSs) and regulatory mechanisms governing their expression.

Maintenance of variation in virulence and reproduction in populations of an agricultural plant pathogen

Genetic diversity within pathogen populations is critically important for predicting pathogen evolution, disease outcomes and prevalence. However, we lack a good understanding of the processes maintaining genetic variation and constraints on pathogen life-history evolution. Here, we analysed interactions between 12 wheat host genotypes and 145 strains of Zymoseptoria tritici from five global populations to investigate the evolution and maintenance of variation in pathogen virulence and reproduction. We found a strong positive correlation between virulence (amount of leaf necrosis) and reproduction (pycnidia density within lesions), with substantial variation in both traits maintained within populations. On average, highly virulent isolates exhibited higher reproduction, which might increase transmission potential in agricultural fields planted to homogeneous hosts at a high density. We further showed that pathogen strains with a narrow host range (i.e. specialists) for reproduction were on average less virulent, and those with a broader host range (i.e. generalists) were on average less fecund on a given specific host. These costs associated with adaptation to different host genotypes might constrain the emergence of generalists by disrupting the directional evolution of virulence and fecundity. We conclude that selection favouring pathogen strains that are virulent across diverse hosts, coupled with selection that maximizes fecundity on specific hosts, may explain the maintenance of these pathogenicity traits within and among populations.

Partial pathogenicity chromosomes in Fusarium oxysporum are sufficient to cause disease and can be horizontally transferred

In Fusarium oxysporum f.sp. lycopersici, all effector genes reported so far - also called SIX genes - are located on a single accessory chromosome which is required for pathogenicity and can also be horizontally transferred to another strain. To narrow down the minimal region required for virulence, we selected partial pathogenicity chromosome deletion strains by fluorescence-assisted cell sorting of a strain in which the two arms of the pathogenicity chromosome were labelled with GFP and RFP respectively. By testing the virulence of these deletion mutants, we show that the complete long arm and part of the short arm of the pathogenicity chromosome are not required for virulence. In addition, we demonstrate that smaller versions of the pathogenicity chromosome can also be transferred to a non-pathogenic strain and they are sufficient to turn the non-pathogen into a pathogen. Surprisingly, originally non-pathogenic strains that had received a smaller version of the pathogenicity chromosome were much more aggressive than recipients with a complete pathogenicity chromosome. Whole genome sequencing analysis revealed that partial deletions of the pathogenicity chromosome occurred mainly close to repeats, and that spontaneous duplication of sequences in accessory regions is frequent both in chromosome deletion strains and in horizontal transfer strains.

Accessory Chromosomes in Fusarium oxysporum

Most genomes within the species complex of Fusarium oxysporum are organized into two compartments: the core chromosomes (CCs) and accessory chromosomes (ACs). As opposed to CCs, which are conserved and vertically transmitted to carry out essential housekeeping functions, lineage- or strain-specific ACs are believed to be initially horizontally acquired through unclear mechanisms. These two genomic compartments are different in terms of gene density, the distribution of transposable elements, and epigenetic markers. Although common in eukaryotes, the functional importance of ACs is uniquely emphasized among fungal species, specifically in relationship to fungal pathogenicity and their adaptation to diverse hosts. With a focus on the cross-kingdom fungal pathogen F. oxysporum, this review provides a summary of the differences between CCs and ACs based on current knowledge of gene functions, genome structures, and epigenetic signatures, and explores the transcriptional crosstalk between the core and accessory genomes.

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

Threats to global food security from emerging fungal and oomycete crop pathogens

Emerging fungal and oomycete pathogens infect staple calorie crops and economically important commodity crops, thereby posing a significant risk to global food security. Our current agricultural systems - with emphasis on intensive monoculture practices - and globalized markets drive the emergence and spread of new pathogens and problematic traits, such as fungicide resistance. Climate change further promotes the emergence of pathogens on new crops and in new places. Here we review the factors affecting the introduction and spread of pathogens and current disease control strategies, illustrating these with the historic example of the Irish potato famine and contemporary examples of soybean rust, wheat blast and blotch, banana wilt and cassava root rot. Our Review looks to the future, summarizing what we see as the main challenges and knowledge gaps, and highlighting the direction that research must take to face the challenge of emerging crop pathogens.

The Many Questions about Mini Chromosomes in Colletotrichum spp

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

Integrative Omics Analysis Reveals a Limited Transcriptional Shock After Yeast Interspecies Hybridization

The formation of interspecific hybrids results in the coexistence of two diverged genomes within the same nucleus. It has been hypothesized that negative epistatic interactions and regulatory interferences between the two sub-genomes may elicit a so-called genomic shock involving, among other alterations, broad transcriptional changes. To assess the magnitude of this shock in hybrid yeasts, we investigated the transcriptomic differences between a newly formed Saccharomyces cerevisiae × Saccharomyces uvarum diploid hybrid and its diploid parentals, which diverged ∼20 mya. RNA sequencing (RNA-Seq) based allele-specific expression (ASE) analysis indicated that gene expression changes in the hybrid genome are limited, with only ∼1-2% of genes significantly altering their expression with respect to a non-hybrid context. In comparison, a thermal shock altered six times more genes. Furthermore, differences in the expression between orthologous genes in the two parental species tended to be diminished for the corresponding homeologous genes in the hybrid. Finally, and consistent with the RNA-Seq results, we show a limited impact of hybridization on chromatin accessibility patterns, as assessed with assay for transposase-accessible chromatin using sequencing (ATAC-Seq). Overall, our results suggest a limited genomic shock in a newly formed yeast hybrid, which may explain the high frequency of successful hybridization in these organisms.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

Repeat-Induced Point Mutations Drive Divergence between Fusarium circinatum and Its Close Relatives

The Repeat-Induced Point (RIP) mutation pathway is a fungal-specific genome defense mechanism that counteracts the deleterious effects of transposable elements. This pathway permanently mutates its target sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence of RIP in the pitch canker pathogen, Fusarium circinatum, and its close relatives in the Fusarium fujikuroi species complex (FFSC). Our results showed that the examined fungi all exhibited hallmarks of RIP, but that they differed in terms of the extent to which their genomes were affected by this pathway. RIP mutations constituted a large proportion of all the FFSC genomes, including both core and dispensable chromosomes, although the latter were generally more extensively affected by RIP. Large RIP-affected genomic regions were also much more gene sparse than the rest of the genome. Our data further showed that RIP-directed sequence diversification increased the variability between homologous regions of related species, and that RIP-affected regions can interfere with homologous recombination during meiosis, thereby contributing to post-mating segregation distortion. Taken together, these findings suggest that RIP can drive the independent divergence of chromosomes, alter chromosome architecture, and contribute to the divergence among F. circinatum and other members of this economically important group of fungi.

Interspecific Gene Exchange Introduces High Genetic Variability in Crop Pathogen

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

Genomic features and evolution of the conditionally dispensable chromosome in the tangerine pathotype of Alternaria alternata

The tangerine pathotype of the ascomycete fungus Alternaria alternata is the causal agent of citrus brown spot, which can result in significant losses of both yield and marketability for tangerines worldwide. A conditionally dispensable chromosome (CDC), which harbours the host-selective ACT toxin gene cluster, is required for tangerine pathogenicity of A. alternata. To understand the genetic makeup and evolution of the tangerine pathotype CDC, we isolated and sequenced the CDCs of the A. alternata Z7 strain and analysed the function and evolution of their genes. The A. alternata Z7 strain has two CDCs (~1.1 and ~0.8 Mb, respectively), and the longer Z7 CDC contains all but one contig of the shorter one. Z7 CDCs contain 254 predicted protein-coding genes, which are enriched in functional categories associated with 'metabolic process' (55 genes, P = 0.037). Relatively few of the CDC genes can be classified as carbohydrate-active enzymes (CAZymes) (4) and transporters (19) and none as kinases. Evolutionary analysis of the 254 CDC proteins showed that their evolutionary conservation tends to be restricted within the genus Alternaria and that the CDC genes evolve faster than genes in the essential chromosomes, likely due to fewer selective constraints. Interestingly, phylogenetic analysis suggested that four of the 25 genes responsible for the ACT toxin production were likely transferred from Colletotrichum (Sordariomycetes). Functional experiments showed that two of them are essential for the virulence of the tangerine pathotype of A. alternata. These results provide new insights into the function and evolution of CDC genes in Alternaria.

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

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

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

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