Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

Candida albicans natural diversity: a resource to dissect fungal commensalism and pathogenesis

Candida albicans is a ubiquitous fungus of humans. It is not only a component of the oral and intestinal microbiota of most healthy adults but also a major cause of mucosal disorders and life-threatening disseminated infections. Until recently, research on the biology and pathogenesis of the fungus was largely based on a single clinical isolate. We review investigations that have started to dissect a diverse set of C. albicans strains. Using different approaches to leverage the species' phenotypic and/or genetic diversity, these studies illuminate the wide range of interactions between fungus and host. While connecting genetic variants to phenotypes of interest remains challenging, research on C. albicans' natural diversity is central to understand fungal commensalism and pathogenesis.

Genetic Diversity and Classification of Colletotrichum sublineola Pathotypes Using a Standard Set of Sorghum Differentials

Anthracnose, incited by Colletotrichum sublineola, is the most destructive foliar disease of sorghum and, under severe conditions, yield losses can exceed 80% on susceptible cultivars. The hyper-variable nature of the pathogen makes its management challenging despite the occurrence of several resistant sources. In this study, the genetic variability and pathogenicity of 140 isolates of C. sublineola, which were sequenced using restriction site-associated sequencing (RAD-Seq), resulted in 1244 quality SNPs. The genetic relationship based on the SNP data showed low to high genetic diversity based on isolates' origin. Isolates from Georgia and North Carolina were grouped into multiple clusters with some level of genetic relationships to each other. Even though some isolates from Texas formed a cluster, others clustered with isolates from Puerto Rico. The isolates from Puerto Rico showed scattered distribution, indicating the diverse nature of these isolates. A population structure and cluster analysis revealed that the genetic variation was stratified into eight populations and one admixture group. The virulence pattern of 30 sequenced isolates on 18 sorghum differential lines revealed 27 new pathotypes. SC748-5, SC112-14, and Brandes were resistant to all the tested isolates, while BTx623 was susceptible to all. Line TAM428 was susceptible to all the pathotypes, except for pathotype 26. Future use of the 18 differentials employed in this study, which contains cultivars/lines which have been used in the Americas, Asia, and Africa, could allow for better characterization of C. sublineola pathotypes at a global level, thus accelerating the development of sorghum lines with stable resistance to the anthracnose pathogen.

Phenotypic Variants of Azole-Resistant Aspergillus Fumigatus that Co-exist in Human Respiratory Samples are Genetically Highly Related

Respiratory specimens obtained from patients with chronic forms of aspergillosis contain phenotypic variants of azole-resistant Aspergillus fumigatus (ARAF) that co-exist in the airway. Here we aimed to study whether phenotypic variants of ARAF that co-exist in clinical specimens were genetically distinct. A panel of six phenotypic variants of ARAF cultured from two sputum samples collected from two patients with chronic aspergillosis were included. Preliminary identification of all isolates was obtained using MALDI–ToF mass spectrometry and confirmed by AsperGenius^® real-time PCR assay. Antifungal susceptibility testing was determined using EUCAST E.Def 9.3 microbroth dilution. Genomic DNA libraries were constructed with the Illumina TruSeq Nano kit. Prepared whole-genome libraries were sequenced on an Illumina HiSeq 2500. Whole genome data were converted into presence/absence of a SNP with respect to the Af293 reference genome. Colonies of ARAF that co-existed in one respiratory sample demonstrated marked phenotypic diversity. Two cyp51A polymorphisms were found among azole-resistant isolates: TR₃₄/L98H/T289A/I364V/G448S was consistently present in four variants with a pan-azole resistant phenotype and TR₃₄/L98H was detected in two variants (itraconazole MIC > 16 mg/L). WGS typing showed that despite marked phenotypic variation, each sample contained a population of highly genetically related azole-resistant A. fumigatus variants. Our SNP analysis suggest that mechanisms additional to genetic-based variation are responsible for phenotypic diversity. Our data demonstrate that the phenotypic variants of ARAF that co-exist in clinical specimens are highly clonal and strongly suggest their origination from a single common ancestor.

Molecular and Morphological Characterization of Exserohilum turcicum (Passerini) Leonard and Suggs Causing Northern Corn Leaf Blight of Maize in Bihar

Maize is considered the third most important cereal crop in Asia after rice and wheat. Many diseases affect this crop due to the cultivation of various hybrids. This research aimed to characterize the causative agent of northern corn leaf blight disease in Bihar, India, caused by Exserohilum turcicum (Passerini) Leonard and Suggs. Leaf samples were collected from infected fields in five maize growing districts of Bihar in 2020-2022. A total of 45 fungal isolates from 135 samples were examined for cultural, morphological, and molecular characteristics and were identified as E. turcicum. The isolates were grouped into four groups based on colony color, i.e., olivaceous brown, blackish brown, whitish black, and grayish, and into two groups based on regular and irregular margins. The conidial shapes were observed to be elongated and spindle-shaped with protruding hilum, with conidial septa ranging from 2-12. Similarly, conidial length varied from 52.94 μm to 144.12 μm. β-tubulin gene sequences analysis made it possible to verify the identities of fungal strains and the phylogenetic relationships of all isolates, which were clustered in the same clade. The β-tubulin gene sequences of all the isolates showed a high level of similarity (100%) with reference isolates from GenBank accession numbers KU670342.1, KU670344.1, KU670343.1, KU670341.1, and KU670340.1. The findings of this study will serve as a baseline for future studies and will help to minimize yield losses.

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.

Evolution of the human pathogenic lifestyle in fungi

Fungal pathogens cause more than a billion human infections every year, resulting in more than 1.6 million deaths annually. Understanding the natural history and evolutionary ecology of fungi is helping us understand how disease-relevant traits have repeatedly evolved. Different types and mechanisms of genetic variation have contributed to the evolution of fungal pathogenicity and specific genetic differences distinguish pathogens from non-pathogens. Insights into the traits, genetic elements, and genetic and ecological mechanisms that contribute to the evolution of fungal pathogenicity are crucial for developing strategies to both predict emergence of fungal pathogens and develop drugs to combat them.

Genomic landscape of a relict fir-associated fungus reveals rapid convergent adaptation towards endophytism

Comparative and pan-genomic analyses of the endophytic fungus Pezicula neosporulosa (Helotiales, Ascomycota) from needles of the relict fir, Abies beshanzuensis, showed expansions of carbohydrate metabolism and secondary metabolite biosynthetic genes characteristic for unrelated plant-beneficial helotialean, such as dark septate endophytes and ericoid mycorrhizal fungi. The current species within the relatively young Pliocene genus Pezicula are predominantly saprotrophic, while P. neosporulosa lacks such features. To understand the genomic background of this putatively convergent evolution, we performed population analyses of 77 P. neosporulosa isolates. This revealed a mosaic structure of a dozen non-recombining and highly genetically polymorphic subpopulations with a unique mating system structure. We found that one idiomorph of a probably duplicated mat1-2 gene was found in putatively heterothallic isolates, while the other co-occurred with mat1-1 locus suggesting homothallic reproduction for these strains. Moreover, 24 and 81 genes implicated in plant cell-wall degradation and secondary metabolite biosynthesis, respectively, showed signatures of the balancing selection. These findings highlight the evolutionary pattern of the two gene families for allowing the fungus a rapid adaptation towards endophytism and facilitating diverse symbiotic interactions.

Comparative Genomics and Gene Pool Analysis Reveal the Decrease of Genome Diversity and Gene Number in Rice Blast Fungi by Stable Adaption with Rice

Magnaporthe oryzae caused huge losses in rice and wheat production worldwide. Comparing to long-term co-evolution history with rice, wheat-infecting isolates were new-emerging. To reveal the genetic differences between rice and wheat blast on global genomic scale, 109 whole-genome sequences of M. oryzae from rice, wheat, and other hosts were reanalyzed in this study. We found that the rice lineage had gone through stronger selective sweep and fewer conserved genes than those of Triticum and Lolium lineages, which indicated that rice blast fungi adapted to rice by gene loss and rapid evolution of specific loci. Furthermore, 228 genes associated with host adaptation of M. oryzae were found by presence/absence variation (PAV) analyses. The functional annotation of these genes found that the fine turning of genes gain/loss involved with transport and transcription factor, thiol metabolism, and nucleotide metabolism respectively are major mechanisms for rice adaption. This result implies that genetic base of specific host plant may lead to gene gain/loss variation of pathogens, so as to enhance their adaptability to host. Further characterization of these specific loci and their roles in adaption and evaluation of the fungi may eventually lead to understanding of interaction mechanism and develop new strategies of the disease management.

Fungi as mediators linking organisms and ecosystems

Fungi form a major and diverse component of most ecosystems on Earth. They are both micro and macroorganisms with high and varying functional diversity as well as great variation in dispersal modes. With our growing knowledge of microbial biogeography, it has become increasingly clear that fungal assembly patterns and processes differ from other microorganisms such as bacteria, but also from macroorganisms such as plants. The success of fungi as organisms and their influence on the environment lies in their ability to span multiple dimensions of time, space, and biological interactions, that is not rivalled by other organism groups. There is also growing evidence that fungi mediate links between different organisms and ecosystems, with the potential to affect the macroecology and evolution of those organisms. This suggests that fungal interactions are an ecological driving force, interconnecting different levels of biological and ecological organisation of their hosts, competitors, and antagonists with the environment and ecosystem functioning. Here we review these emerging lines of evidence by focusing on the dynamics of fungal interactions with other organism groups across various ecosystems. We conclude that the mediating role of fungi through their complex and dynamic ecological interactions underlie their importance and ubiquity across Earth's ecosystems.

Neurospora from Natural Populations: Population Genomics Insights into the Life History of a Model Microbial Eukaryote

The ascomycete filamentous fungus Neurospora crassa played a historic role in experimental biology and became a model system for genetic research. Stimulated by a systematic effort to collect wild strains initiated by Stanford geneticist David Perkins, the genus Neurospora has also become a basic model for the study of evolutionary processes, speciation, and population biology. In this chapter, we will first trace the history that brought Neurospora into the era of population genomics. We will then cover the major contributions of population genomic investigations using Neurospora to our understanding of microbial biogeography and speciation, and review recent work using population genomics and genome-wide association mapping that illustrates the unique potential of Neurospora as a model for identifying the genetic basis of (potentially adaptive) phenotypes in filamentous fungi. The advent of population genomics has contributed to firmly establish Neurospora as a complete model system and we hope our review will entice biologists to include Neurospora in their research.

One Clonal Lineage of Calonectria pseudonaviculata Is Primarily Responsible for the Boxwood Blight Epidemic in the United States

Boxwood blight caused by Calonectria pseudonaviculata and C. henricotiae is destroying cultivated and native boxwood worldwide, with profound negative economic impacts on the horticulture industry. First documented in the United States in 2011, the disease has now occurred in 30 states. Previous research showed that global C. pseudonaviculata populations prior to 2014 had a clonal structure, and only the MAT1-2 idiomorph was observed. In this study, we examined C. pseudonaviculata genetic diversity and population structure in the United States after 2014, following the expansion of the disease across the country over the past 5 years. Two hundred eighteen isolates from 21 states were genotyped by sequencing 11 simple sequence repeat (SSR) loci and by MAT1 idiomorph typing. All isolates presented C. pseudonaviculata-specific alleles, indicating that C. henricotiae is still absent in the U.S. states sampled. The presence of only the MAT1-2 idiomorph and gametic linkage disequilibrium suggests the prevalence of asexual reproduction. The contemporary C. pseudonaviculata population is characterized by a clonal structure and composed of 13 multilocus genotypes (SSR-MLGs) unevenly distributed across the United States. These SSR-MLGs grouped into two clonal lineages (CLs). The predominant lineage CL2 (93% of isolates) is the primary contributor to U.S. disease expansion. The contemporary U.S. C. pseudonaviculata population is not geographically subdivided and not genetically differentiated from the U.S. population prior to 2014, but is significantly differentiated from the main European population, which is largely composed of CL1. Our findings provide insights into the boxwood blight epidemic that are critical for disease management and breeding of resistant boxwood cultivars.

The Rise of Fungi: A Report on the CIFAR Program Fungal Kingdom: Threats & Opportunities Inaugural Meeting

The first meeting of the CIFAR Fungal Kingdom: Threats & Opportunities research program saw the congregation of experts on fungal biology to address the most pressing threats fungi pose to global health, agriculture, and biodiversity. This report covers the research discussed during the meeting and the advancements made toward mitigating the devastating impact of fungi on plants, animals, and humans.

Identification of CSP Types and Genotypic Variability of Clinical and Environmental Isolates of Aspergillus fumigatus from Different Geographic Origins

The CSP (cell surface protein) microsatellite marker is useful for typing Aspergillus fumigatus isolates and determining relationships at the subpopulation level because it has shown high discriminatory power. In the present study, 90 A. fumigatus isolates from Mexico (MX), Argentina (AR), France (FR), and Peru (PE) were identified through a phylogenetic analysis using the benA gene fragment and were typed with the CSP microsatellite, and the types were identified using the nomenclature recommended in the literature. Genetic variability was analyzed through haplotype diversity, nucleotide diversity, polymorphic sites, and nucleotide differences between pairs of sequences. The population structure was evaluated using the Tajima's D statistic. No new CSP types were recorded in the MX, FR, and PE isolates, while in the AR isolates, two new CSP types were identified (t25 and t26). The most common CSP types in the studied populations were t01, t02, t03, and t04A; these results are consistent with findings in other countries. In addition, the genetic diversity parameters we obtained revealed that the greatest genetic diversity was found in the MX population, followed by AR and FR. No population structure was identified among the isolates studied.

Colletotrichum: species complexes, lifestyle, and peculiarities of some sources of genetic variability

The genus Colletotrichum comprises species with different lifestyles but is mainly known for phytopathogenic species that infect crops of agronomic relevance causing considerable losses. The fungi of the genus Colletotrichum are distributed in species complexes and within each complex some species have particularities regarding their lifestyle. The most commonly found and described lifestyles in Colletotrichum are endophytic and hemibiotrophic phytopathogenic. Several of these phytopathogenic species show wide genetic variability, which makes long-term maintenance of resistance in plants difficult. Different mechanisms may play an important role in the emergence of genetic variants but are not yet fully understood in this genus. These mechanisms include heterokaryosis, a parasexual cycle, sexual cycle, transposable element activity, and repeat-induced point mutations. This review provides an overview of the genus Colletotrichum, the species complexes described so far and the most common lifestyles in the genus, with a special emphasis on the mechanisms that may be responsible, at least in part, for the emergence of new genotypes under field conditions.

Fungal community assembly in drought-stressed sorghum shows stochasticity, selection, and universal ecological dynamics

Community assembly of crop-associated fungi is thought to be strongly influenced by deterministic selection exerted by the plant host, rather than stochastic processes. Here we use a simple, sorghum system with abundant sampling to show that stochastic forces (drift or stochastic dispersal) act on fungal community assembly in leaves and roots early in host development and when sorghum is drought stressed, conditions when mycobiomes are small. Unexpectedly, we find no signal for stochasticity when drought stress is relieved, likely due to renewed selection by the host. In our experimental system, the host compartment exerts the strongest effects on mycobiome assembly, followed by the timing of plant development and lastly by plant genotype. Using a dissimilarity-overlap approach, we find a universality in the forces of community assembly of the mycobiomes of the different sorghum compartments and in functional guilds of fungi.

Fungal community assembly on crop plants is thought to be driven by deterministic selection exerted by the host. Here Gao et al. use a sorghum system to show that stochastic forces act on fungal community assembly in leaves and roots early in host development and when sorghum is drought stressed.

Population Structure and Genetic Diversity among Isolates of Coccidioides posadasii in Venezuela and Surrounding Regions

Coccidioides posadasii is a pathogenic fungus that causes coccidioidomycosis in many arid regions of the Americas. One of these regions is bordered by the Caribbean Sea, and the surrounding landscape may play an important role in the dispersion of C. posadasii across South America through southeastern Mexico, Honduras, Guatemala, and Venezuela. Comparative phylogenomic analyses of C. posadasii reveal that clinical strains from Venezuela are genetically distinct from the North American populations found in (i) Arizona and (ii) Texas, Mexico, and the rest of South America (TX/MX/SA). We find evidence for admixture between the Venezuela and the North American populations of C. posadasii in Central America. Additionally, the proportion of Venezuelan alleles in the admixed population decreases as latitude (and distance from Venezuela) increases. Our results indicate that the population in Venezuela may have been subjected to a recent bottleneck and shows a strong population structure. This analysis provides insight into potential for Coccidioides spp. to invade new regions.IMPORTANCE Valley Fever is a fungal disease caused by two species of fungi: Coccidioides immitis and C. posadasii These fungi are found throughout the arid regions of North and South America; however, our understanding of genetic diversity and disease in South America is limited. In this report, we analyze 10 new genomes of Coccidioides posadasii from regions bordering the Caribbean Sea. We show that these populations are distinct and that isolates from Venezuela are likely a result of a recent bottleneck. These data point to patterns that might be observed when investigating recently established populations.

Fungal clones win the battle, but recombination wins the war

Clonal reproduction is common in fungi and fungal-like organisms during epidemics and invasion events. The success of clonal fungi shaped systems for their classification and some pathogens are tacitly treated as asexual. We argue that genetic recombination driven by sexual reproduction must be a starting hypothesis when dealing with fungi for two reasons: (1) Clones eventually crash because they lack adaptability; and (2) fungi find a way to exchange genetic material through recombination, whether sexual, parasexual, or hybridisation. Successful clones may prevail over space and time, but they are the product of recombination and the next successful clone will inevitably appear. Fungal pathogen populations are dynamic rather than static, and they need genetic recombination to adapt to a changing environment.

Correlating Genotype and Phenotype in the Asexual Yeast Candida orthopsilosis Implicates ZCF29 in Sensitivity to Caffeine

Candida orthopsilosis is diploid asexual yeast that causes human disease. Most C. orthopsilosis isolates arose from at least four separate hybridizations between related, but not identical, parents. Here, we used population genomics data to correlate genotypic and phenotypic variation in 28 C. orthopsilosis isolates. We used cosine similarity scores to identify 65 variants with potential high-impact (deleterious effects) that correlated with specific phenotypes. Of these, 19 were Single Nucleotide Polymorphisms (SNPs) that changed stop or start codons, or splice sites. One variant resulted in a premature stop codon in both alleles of the gene ZCF29 in C. orthopsilosis isolate 185, which correlated with sensitivity to nystatin and caffeine. We used CRISPR-Cas9 editing to introduce this polymorphism into two resistant C. orthopsilosis isolates. Introducing the stop codon resulted in sensitivity to caffeine and to ketoconazole, but not to nystatin. Our analysis shows that it is possible to associate genomic variants with phenotype in asexual Candida species, but that only a small amount of genomic variation can be easily explored.

Understanding Adaptation, Coevolution, Host Specialization, and Mating System in Castrating Anther-Smut Fungi by Combining Population and Comparative Genomics

Anther-smut fungi provide a powerful system to study host-pathogen specialization and coevolution, with hundreds of Microbotryum species specialized on diverse Caryophyllaceae plants, castrating their hosts through manipulation of the hosts' reproductive organs to facilitate disease transmission. Microbotryum fungi have exceptional genomic characteristics, including dimorphic mating-type chromosomes, that make this genus anexcellent model for studying the evolution of mating systems and their influence on population genetics structure and adaptive potential. Important insights into adaptation, coevolution, host specialization, and mating system evolution have been gained using anther-smut fungi, with new insights made possible by the recent advent of genomic approaches. We illustrate with Microbotryum case studies how using a combination of comparative genomics, population genomics, and transcriptomics approaches enables the integration of different evolutionary perspectives across different timescales. We also highlight current challenges and suggest future studies that will contribute to advancing our understanding of the mechanisms underlying adaptive processes in populations of fungal pathogens.

First microsatellite markers developed and applied for the genetic diversity study and population structure of Didymella pisi associated with ascochyta blight of dry pea in Montana

Didymella pisi is the predominant causal pathogen of ascochyta blight of dry pea causing yield losses in Montana, where 415 000 acres were planted to dry pea in 2018. Thirty-three microsatellite markers were developed for dry pea pathogenic fungus, Didymella pisi, these markers were used to analyze genetic diversity and population structure of 205 isolates from four different geographical regions of Montana. These loci produced a total of 216 alleles with an average of 1.63 alleles per microsatellite marker. The polymorphic information content values ranged from 0.020 to 0.990 with an average of 0.323. The average observed heterozygosity across all loci varied from 0.000 to 0.018. The gene diversity among the loci ranged from 0.003 to 0.461. Unweighted Neighbor-joining and population structure analysis grouped these 205 isolates into two major sub-groups. The clusters did not match the geographic origin of the isolates. Analysis of molecular variance showed 85 % of the total variation within populations and only 15 % among populations. There was moderate genetic variation in the total populations (PhiPT = 0.153). Information obtained from this study could be useful as a base to design strategies for improved management such as breeding for resistance to ascochyta blight of dry pea in Montana.

A gene graveyard in the genome of the fungus Podospora comata

Mechanisms involved in fine adaptation of fungi to their environment include differential gene regulation associated with single nucleotide polymorphisms and indels (including transposons), horizontal gene transfer, gene copy amplification, as well as pseudogenization and gene loss. The two Podospora genome sequences examined here emphasize the role of pseudogenization and gene loss, which have rarely been documented in fungi. Podospora comata is a species closely related to Podospora anserina, a fungus used as model in several laboratories. Comparison of the genome of P. comata with that of P. anserina, whose genome is available for over 10 years, should yield interesting data related to the modalities of genome evolution between these two closely related fungal species that thrive in the same types of biotopes, i.e., herbivore dung. Here, we present the genome sequence of the mat + isolate of the P. comata reference strain T. Comparison with the genome of the mat + isolate of P. anserina strain S confirms that P. anserina and P. comata are likely two different species that rarely interbreed in nature. Despite having a 94-99% of nucleotide identity in the syntenic regions of their genomes, the two species differ by nearly 10% of their gene contents. Comparison of the species-specific gene sets uncovered genes that could be responsible for the known physiological differences between the two species. Finally, we identified 428 and 811 pseudogenes (3.8 and 7.2% of the genes) in P. anserina and P. comata, respectively. Presence of high numbers of pseudogenes supports the notion that difference in gene contents is due to gene loss rather than horizontal gene transfers. We propose that the high frequency of pseudogenization leading to gene loss in P. anserina and P. comata accompanies specialization of these two fungi. Gene loss may be more prevalent during the evolution of other fungi than usually thought.

Gene Presence-Absence Polymorphism in Castrating Anther-Smut Fungi: Recent Gene Gains and Phylogeographic Structure

Gene presence-absence polymorphisms segregating within species are a significant source of genetic variation but have been little investigated to date in natural populations. In plant pathogens, the gain or loss of genes encoding proteins interacting directly with the host, such as secreted proteins, probably plays an important role in coevolution and local adaptation. We investigated gene presence-absence polymorphism in populations of two closely related species of castrating anther-smut fungi, Microbotryum lychnidis-dioicae (MvSl) and M. silenes-dioicae (MvSd), from across Europe, on the basis of Illumina genome sequencing data and high-quality genome references. We observed presence-absence polymorphism for 186 autosomal genes (2% of all genes) in MvSl, and only 51 autosomal genes in MvSd. Distinct genes displayed presence-absence polymorphism in the two species. Genes displaying presence-absence polymorphism were frequently located in subtelomeric and centromeric regions and close to repetitive elements, and comparison with outgroups indicated that most were present in a single species, being recently acquired through duplications in multiple-gene families. Gene presence-absence polymorphism in MvSl showed a phylogeographic structure corresponding to clusters detected based on SNPs. In addition, gene absence alleles were rare within species and skewed toward low-frequency variants. These findings are consistent with a deleterious or neutral effect for most gene presence-absence polymorphism. Some of the observed gene loss and gain events may however be adaptive, as suggested by the putative functions of the corresponding encoded proteins (e.g., secreted proteins) or their localization within previously identified selective sweeps. The adaptive roles in plant and anther-smut fungi interactions of candidate genes however need to be experimentally tested in future studies.