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

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

The Impact of Recombination Hotspots on Genome Evolution of a Fungal Plant Pathogen.

Recombination has an impact on genome evolution by maintaining chromosomal integrity, affecting the efficacy of selection, and increasing genetic variability in populations. Recombination rates are a key determinant of the coevolutionary dynamics between hosts and their pathogens. Historic recombination events created devastating new pathogens, but the impact of ongoing recombination in sexual pathogens is poorly understood. Many fungal pathogens of plants undergo regular sexual cycles, and sex is considered to be a major factor contributing to virulence. We generated a recombination map at kilobase-scale resolution for the haploid plant pathogenic fungus Zymoseptoria tritici. To account for intraspecific variation in recombination rates, we constructed genetic maps from two independent crosses. We localized a total of 10,287 crossover events in 441 progeny and found that recombination rates were highly heterogeneous within and among chromosomes. Recombination rates on large chromosomes were inversely correlated with chromosome length. Short accessory chromosomes often lacked evidence for crossovers between parental chromosomes. Recombination was concentrated in narrow hotspots that were preferentially located close to telomeres. Hotspots were only partially conserved between the two crosses, suggesting that hotspots are short-lived and may vary according to genomic background. Genes located in hotspot regions were enriched in genes encoding secreted proteins. Population resequencing showed that chromosomal regions with high recombination rates were strongly correlated with regions of low linkage disequilibrium. Hence, genes in pathogen recombination hotspots are likely to evolve faster in natural populations and may represent a greater threat to the host.

QTL mapping of fungicide sensitivity reveals novel genes and pleiotropy with melanization in the pathogen Zymoseptoria tritici.

A major problem associated with the intensification of agriculture is the emergence of fungicide resistance. Azoles are ergosterol biosynthesis inhibitors that have been widely used in agriculture and medicine since the 1970s, leading to emergence of increasingly resistant fungal populations. The known genetic mechanisms underlying lower azole sensitivity include mutations affecting the CYP51 gene that encodes the target protein, but in many cases azole resistance is a more complex trait with an unknown genetic basis. We used quantitative trait locus (QTL) mapping to identify genes affecting azole sensitivity in two crosses of Zymoseptoria tritici, the most damaging wheat pathogen in Europe. Restriction site associated DNA sequencing (RADseq) was used to genotype 263 (cross 1) and 261 (cross 2) progeny at ∼ 8500 single nucleotide polymorphisms (SNP) and construct two dense linkage maps. Azole sensitivity was assessed using high-throughput digital image analysis of colonies growing on Petri dishes with or without the fungicide propiconazole. We identified three QTLs for azole sensitivity, including two that contained novel fungicide sensitivity genes. One of these two QTLs contained only 16 candidate genes, among which four most likely candidates were identified. The third QTL contained ERG6, encoding another protein involved in ergosterol biosynthesis. Known genes in QTLs affecting colony growth included CYP51 and PKS1, a gene affecting melanization in Z. tritici. PKS1 showed compelling evidence for pleiotropy, with a rare segregating allele that increased melanization while decreasing growth rate and propiconazole sensitivity. This study resolved the genetic architecture of an important agricultural trait and led to identification of novel genes that are likely to affect azole sensitivity in Z. tritici. It also provided insight into fitness costs associated with lowered azole sensitivity and suggests a novel fungicide mixture strategy.

Quantitative trait locus mapping of melanization in the plant pathogenic fungus Zymoseptoria tritici.

Melanin plays an important role in virulence and antimicrobial resistance in several fungal pathogens. The wheat pathogen Zymoseptoria tritici is important worldwide, but little is known about the genetic architecture of pathogenicity, including the production of melanin. Because melanin production can exhibit complex inheritance, we used quantitative trait locus (QTL) mapping in two crosses to identify the underlying genes. Restriction site-associated DNA sequencing was used to genotype 263 (cross 1) and 261 (cross 2) progeny at ~8500 single-nucleotide polymorphisms and construct two dense linkage maps. We measured gray values, representing degrees of melanization, for single-spore colonies growing on Petri dishes by using a novel image-processing approach that enabled high-throughput phenotyping. Because melanin production can be affected by stress, each offspring was grown in two stressful environments and one control environment. We detected six significant QTL in cross 1 and nine in cross 2, with three QTL shared between the crosses. Different QTL were identified in different environments and at different colony ages. By obtaining complete genome sequences for the four parents and analyzing sequence variation in the QTL confidence intervals, we identified 16 candidate genes likely to affect melanization. One of these candidates was PKS1, a polyketide synthase gene known to play a role in the synthesis of dihydroxynaphthalene melanin. Three candidate quantitative trait nucleotides were identified in PKS1. Many of the other candidate genes were not previously associated with melanization.