Epidemics caused by fungal plant pathogens pose a major threat to agro-ecosystems and impact global food security. High-throughput sequencing enabled major advances in understanding how pathogens cause disease on crops. Hundreds of fungal genomes are now available and analyzing these genomes highlighted the key role of effector genes in disease. Effectors are small secreted proteins that enhance infection by manipulating host metabolism. Fungal genomes carry 100s of putative effector genes, but the lack of homology among effector genes, even for closely related species, challenges evolutionary and functional analyses. Furthermore, effector genes are often found in rapidly evolving chromosome compartments which are difficult to assemble. We review how population and comparative genomics toolsets can be combined to address these challenges. We highlight studies that associated genome-scale polymorphisms with pathogen lifestyles and adaptation to different environments. We show how genome-wide association studies can be used to identify effectors and other pathogenicity-related genes underlying rapid adaptation. We also discuss how the compartmentalization of fungal genomes into core and accessory regions shapes the evolution of effector genes. We argue that an understanding of genome evolution provides important insight into the trajectory of host-pathogen co-evolution.

Evolution of fungicide resistance is a major threat to food production in agricultural ecosystems. Fungal pathogens rapidly evolved resistance to all classes of fungicides applied to the field. Resistance to the commonly used azole fungicides is thought to be driven mainly by mutations in a gene (CYP51) encoding a protein of the ergosterol biosynthesis pathway. However, some fungi gained azole resistance independently of CYP51 mutations and the mechanisms leading to CYP51-independent resistance are poorly understood. We used whole-genome sequencing and genome-wide association studies (GWAS) to perform an unbiased screen of azole resistance loci in Rhynchosporium commune, the causal agent of the barley scald disease. We assayed cyproconazole resistance in 120 isolates collected from nine populations worldwide. We found that mutations in highly conserved genes encoding the vacuolar cation channel YVC1, a transcription activator, and a saccharopine dehydrogenase made significant contributions to fungicide resistance. These three genes were not previously known to confer resistance in plant pathogens. However, YVC1 is involved in a conserved stress response pathway known to respond to azoles in human pathogenic fungi. We also performed GWAS to identify genetic polymorphism linked to fungal growth rates. We found that loci conferring increased fungicide resistance were negatively impacting growth rates, suggesting that fungicide resistance evolution imposed costs. Analyses of population structure showed that resistance mutations were likely introduced into local populations through gene flow. Multilocus resistance evolution to fungicides shows how pathogen populations can evolve a complex genetic architecture for an important phenotypic trait within a short time span.

Coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect specific host genotypes. Hence, spatial heterogeneity in host genotypes coupled with abiotic factors could lead to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. Rhynchosporium commune, the pathogen causing barley scald disease, interacts at least partially in a gene-for-gene manner with its host. We analyzed global field populations of 125 R. commune isolates to identify candidate genes for local adaptation. Whole genome sequencing data showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded protein functions related to biotic and abiotic stress responses. Selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Summary Plant pathogens secrete effector proteins to manipulate the host and facilitate infection. Cognate hosts trigger strong defence responses upon detection of these effectors. Consequently, pathogens and hosts undergo rapid coevolutionary arms races driven by adaptive evolution of effectors and receptors. Because of their high rate of turnover, most effectors are thought to be species‐specific and the evolutionary trajectories are poorly understood. Here, we investigate the necrosis‐inducing protein 1 (NIP1) effector in the multihost pathogen genus Rhynchosporium. We retraced the evolutionary history of the NIP1 locus using whole‐genome assemblies of 146 strains covering four closely related species. NIP1 orthologues were present in all species but the locus consistently segregated presence–absence polymorphisms suggesting long‐term balancing selection. We also identified previously unknown paralogues of NIP1 that were shared among multiple species and showed substantial copy‐number variation within R. commune. The NIP1A paralogue was under significant positive selection suggesting that NIP1A is the dominant effector variant coevolving with host immune receptors. Consistent with this prediction, we found that copy number variation at NIP1A had a stronger effect on virulence than NIP1B. Our analyses unravelled the origins and diversification mechanisms of a pathogen effector family shedding light on how pathogens gain adaptive genetic variation.

Transcription factors (TFs) form the major class of regulatory genes and play key roles in multiple plant stress responses. In most eukaryotic plants, transcription factor (TF) families (WRKY, MADS-box and MYB) activate unique cellular-level abiotic and biotic stress-responsive strategies, which are considered as key determinants for defense and developmental processes. Arabidopsis and rice are two important representative model systems for dicot and monocot plants, respectively. A comprehensive comparative study on 101 OsWRKY, 34 OsMADS box and 122 OsMYB genes (rice genome) and, 71 AtWRKY, 66 AtMADS box and 144 AtMYB genes (Arabidopsis genome) showed various relationships among TFs across species. The phylogenetic analysis clustered WRKY, MADS-box and MYB TF family members into 10, 7 and 14 clades, respectively. All clades in WRKY and MYB TF families and almost half of the total number of clades in the MADS-box TF family are shared between both species. Chromosomal and gene structure analysis showed that the Arabidopsis-rice orthologous TF gene pairs were unevenly localized within their chromosomes whilst the distribution of exon–intron gene structure and motif conservation indicated plausible functional similarity in both species. The abiotic and biotic stress-responsive cis-regulatory element type and distribution patterns in the promoter regions of Arabidopsis and rice WRKY, MADS-box and MYB orthologous gene pairs provide better knowledge on their role as conserved regulators in both species. Co-expression network analysis showed the correlation between WRKY, MADs-box and MYB genes in each independent rice and Arabidopsis network indicating their role in stress responsiveness and developmental processes.

The coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect a specific host. Hence, heterogeneity in the host genotypes and abiotic factors leads to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. We analyzed global field populations of Rhynchosporium commune, the pathogen causing barley scald disease, to identify candidate genes for local adaptation. Whole genome sequencing data generated for 125 isolates showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded proteins with functions related to both biotic and abiotic stresses. We found that selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Microbial pathogens can rapidly adapt to changing environments such as the application of pesticides or host resistance. Copy number variations (CNV) are a major source of adaptive genetic variation for recent adaptation. Here, we analyze how a major fungal pathogen of barley, Rhynchosporium commune, has adapted to host environment, fungicide and temperature challenges.We screen the genomes of 126 isolates sampled across a worldwide set of populations and identify a total of 7'879 gene duplications and 116 gene deletions. Most gene duplications result from segmental chromosomal duplications. We find that genes showing recent gains or losses are enriched in functions related to host exploitation (i.e. effectors and cell wall degrading enzymes). We perform a phylogeny-informed genome-wide association study (GWAS) and identify 191 copy-number variants associated with different pathogenesis and temperature related traits, including a large segmental duplication of CYP51A that has contributed to the emergence of azole resistance. Additionally, we use a genome-wide SNP dataset to replicate the GWAS and contrast it with the CNV-focused analysis.We find that frequencies of adaptive CNV alleles show high variation among populations for traits under strong selection such as fungicide resistance. In contrast, adaptive CNV alleles underpinning temperature adaptation tend to be near fixation. Finally, we show that transposable elements are important drivers of recent gene copy-number variation. Loci showing signatures of recent positive selection are enriched in miniature inverted repeat transposons. Our findings show how extensive segmental duplications create the raw material for recent adaptation in global populations of a fungal pathogen.

Glucosinolates (GSLs) and cyanogenic glycosides (CGs) fulfil functions in plant defence and have been reported to be anticancer agents. Generally, GSL-containing plants do not produce CG, and vice versa, CG-containing plants do not synthesise GSLs. However, the production of both GSL and CG compounds was observed in Carica papaya. Additionally, several studies found both GSL glucotropaeolin and CG prunasin in papaya leaves. The advancement of genome technologies can be explored to elucidate the gene functions and other molecular discoveries in plants that might relate to GSLs and CGs. This review aims to discuss the complex interplay of the rare events whereby these two compounds (GSL and CG) co-occur in a bifurcation pathway in papaya. To our knowledge, this is the first review that highlights novel GSL and CG genes in papaya. Furthermore, species-specific pathways in papaya are also discussed and comprehensively described. The transcription factors involved in regulating GSL and CG biosynthesis pathways are also discussed, accompanied by relevant bioinformatic approaches that can help discover potential regulatory genes that control the production of prunasin and glucotropaeolin in papaya.