Transcriptome-Derived Amplicon Sequencing Markers Elucidate the U.S. Podosphaera macularis Population Structure Across Feral and Commercial Plantings of Humulus lupulus

Population Genetic Structure of the Rubber Tree Powdery Mildew Pathogen (Erysiphe quercicola) from China

In order to manage agricultural pathogens, it is crucial to understand the population structure underlying epidemics. Rubber tree powdery mildew, caused by Erysiphe quercicola, is a serious threat to rubber plantations worldwide, especially in subtropical environments including all rubber tree-growing regions in China. However, the population structure of the pathogen is uncertain. In this study, 16 polymorphic microsatellite markers were used to genotype powdery mildew samples from the main rubber tree-growing regions including Yunnan (YN), Hainan (HN), western Guangdong (WG), and eastern Guangdong (EG). YN had higher genotypic diversity (Simpson's indices), genotypic evenness, Nei's gene diversity, allelic richness, and private allelic richness than the other regions. Cluster analysis, discriminant analysis of principal components, pairwise divergence, and shared multilocus genotype analyses all showed that YN differed significantly from the other regions. The genetic differentiation was small among the other three regions (HN, WG, and EG). Analysis of molecular variance indicated that the variability among regions accounted for 22.37% of the total variability. Genetic differentiation was significantly positively correlated (Rxy = 0.772, P = 0.001) with geographic distance. Linkage equilibrium analysis suggested possible occurrence of sexual recombination although asexual reproduction predominates in E. quercicola. The results suggested that although significant genetic differentiation of E. quercicola occurred between YN and the other regions, pathogen populations from the other three regions lacked genetic differentiation.

Identification of quantitative trait loci associated with R1-mediated resistance to powdery mildew and sex determination in hop (Humulus lupulus L.)

KEY MESSAGE

Two QTL were identified using linkage mapping approaches, one on hop linkage group 3 (qHl_Chr3.PMR1) associated with powdery mildew resistance and a second on linkage group 10 (cqHl_ChrX.SDR1) associated with sex determination. Hop (Humulus lupulus L.) is a dioecious species cultivated for use in beer. Hop powdery mildew, caused by Podosphaera macularis, is a constraint in many growing regions. Thus, identifying markers associated with powdery mildew resistance and sex provides the opportunity to pyramid R-genes and select female plants as seedlings, respectively. Our objectives were to characterize the genetic basis of R1-mediated resistance in the cultivar Zenith which provides resistance to pathogen races in the US, identify quantitative trait loci (QTL) associated with R1 and sex, and develop markers for molecular breeding-based approaches. Phenotypic evaluation of the population indicated that R1-based resistance and sex are inherited monogenically. We constructed a genetic map using 1339 single nucleotide polymorphisms (SNPs) based upon genotype-by-sequencing of 128 F₁ progeny derived from a Zenith × USDA 21058M biparental population. SNPs were assigned to 10 linkage groups comprising a map length of 1204.97 cM with an average density of 0.94 cM/marker. Quantitative trait locus mapping identified qHl_Chr3.PMR1, associated with R1 on linkage group 3 (LOD = 23.57, R² = 57.2%), and cqHl_ChrX.SDR1, associated with sex on linkage group 10 (LOD = 5.42, R² = 25.0%). Kompetitive allele-specific PCR (KASP) assays were developed for both QTL and assessed against diverse germplasm. Our results indicate that KASP markers associated with R1 may be limited to materials that are pedigree-related to Zenith, whereas markers associated with sex may be transferable across populations. The high-density map, QTL, and associated KASP markers will enable selecting for sex and R1-mediated resistance in hop.

Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities

Many species of fungal plant pathogens coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. The rice blast fungus Pyricularia oryzae is a widespread model plant pathogen displaying population subdivision. However, most studies of natural variation in P. oryzae have been limited in genomic or geographic resolution, and host adaptation is the only factor that has been investigated extensively as a contributor to population subdivision. In an effort to complement previous studies, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus covering a broad geographical range. Using single-nucleotide polymorphism genotyping data for 886 isolates sampled from 152 sites in 51 countries, we showed that population subdivision of P. oryzae in one recombining and three clonal lineages with broad distributions persisted with deeper sampling. We also extended previous findings by showing further population subdivision of the recombining lineage into one international and three Asian clusters, and by providing evidence that the three clonal lineages of P. oryzae were found in areas with different prevailing environmental conditions, indicating niche separation. Pathogenicity tests and bioinformatic analyses using an extended set of isolates and rice varieties indicated that partial specialization to rice subgroups contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong additional barriers to gene flow between clonal lineages. Our results demonstrate that the spread of a fungal pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct, widely distributed, lineages.

A Comprehensive Characterization of Ecological and Epidemiological Factors Driving Perennation of Podosphaera macularis Chasmothecia on Hop (Humulus lupulus)

Hop powdery mildew, caused by the ascomycete fungus Podosphaera macularis, is a consistent threat to sustainable hop production. The pathogen utilizes two reproductive strategies for overwintering and perennation: (i) asexual vegetative hyphae on dormant buds that emerge the following season as infected shoots; and (ii) sexual ascocarps (chasmothecia), which are discharged during spring rain events. We demonstrate that P. macularis chasmothecia, in the absence of any asexual P. macularis growth forms, are a viable overwintering source capable of causing early season infection two to three orders of magnitude greater than that reported for perennation via asexual growth. Two epidemiological models were defined that describe (i) temperature-driven maturation of P. macularis chasmothecia; and (ii) ascosporic discharge in response to duration of leaf wetness and prevailing temperatures. P. macularis ascospores were confirmed to be infectious at temperatures ranging from 5 to 20°C. The organism's chasmothecia were also found to adhere tightly to the host tissue on which they formed, suggesting that these structures likely overwinter wherever hop tissue senesces within a hop yard. These observations suggest that existing early season disease management practices are especially crucial to controlling hop powdery mildew in the presence of P. macularis chasmothecia. Furthermore, these insights provide a baseline for the validation of weather-driven models describing maturation and release of P. macularis ascospores, models that can eventually be incorporated into hop disease management programs.

High Genetic Diversity in Predominantly Clonal Populations of the Powdery Mildew Fungus Podosphaera leucotricha from U.S. Apple Orchards

Apple powdery mildew (APM), caused by Podosphaera leucotricha, is a constant threat to apple production worldwide. Very little is known about the biology and population structure of this pathogen in the United States and other growing regions, which affects APM management. A total of 253 P. leucotricha isolates, sampled from 10 apple orchards in Washington, New York, and Virginia, were genetically characterized with novel single sequence repeat and mating type markers. Eighty-three multilocus genotypes (MLGs) were identified, most of which were unique to a given orchard. Each isolate carried either a MAT1-1 or a MAT1-2 idiomorph at the mating type locus, indicating that P. leucotricha is heterothallic. Virulence tests on detached apple leaves showed that the 10 most frequent P. leucotricha MLGs were avirulent on a line containing a major resistance gene. Analysis of molecular variance showed significant differentiation (P < 0.001) among populations, a result supported by principal coordinate analysis revealing three genetic groups, each represented by nonoverlapping MLGs from Washington, New York, and Virginia. A Bayesian cluster analysis showed genetic heterogeneity between Washington populations, and a relative migration analysis indicated substantial gene flow among neighboring orchards. Random mating tests indicated that APM epidemics during the active cycle were dominated by clonal reproduction. However, the presence of sexual structures in orchards, the likelihood that five repeated MLGs resulted from sexual reproduction, and high genotypic diversity observed in some populations suggest that sexual spores play some role in APM epidemics. IMPORTANCE Understanding the population biology and epidemiology of plant pathogens is essential to develop effective strategies for controlling plant diseases. Herein, we gathered insights into the population biology of P. leucotricha populations from conventional and organic apple orchards in the United States. We showed genetic heterogeneity between P. leucotricha populations in Washington and structure between populations from different U.S. regions, suggesting that short-distance spore dispersal plays an important role in the disease's epidemiology. We presented evidence that P. leucotricha is heterothallic and that populations likely result from a mixed (i.e., sexual and asexual) reproductive system, revealing that the sexual stage contributes to apple powdery mildew epidemics. We showed that the major resistance gene Pl-1 is valuable for apple breeding because virulent isolates have most likely not emerged yet in U.S. commercial orchards. These results will be important to achieve sustainability of disease management strategies and maintenance of plant health in apple orchards.

Population Genetic Analyses of Botrytis cinerea Isolates From Michigan Vineyards Using a High-Throughput Marker System Approach

As sequencing costs continue to decrease, new tools are being developed for assessing pathogen diversity and population structure. Traditional marker types, such as microsatellites, are often more cost effective than single-nucleotide polymorphism (SNP) panels when working with small numbers of individuals, but may not allow for fine scale evaluation of low or moderate structure in populations. Botrytis cinerea is a necrotrophic plant pathogen with high genetic variability that can infect more than 200 plant species worldwide. A panel of 52 amplicons were sequenced for 82 isolates collected from four Michigan vineyards representing 2 years of collection and varying fungicide resistance. A panel of nine microsatellite markers previously described was also tested across 74 isolates from the same population. A microsatellite and SNP marker analysis of B. cinerea populations was performed to assess the genetic diversity and population structure of Michigan vineyards, and the results from both marker types were compared. Both methods were able to detect population structure associated with resistance to the individual fungicides thiabendazole and boscalid, and multiple fungicide resistance (MFR). Microsatellites were also able to differentiate population structure associated with another fungicide, fluopyram, while SNPs were able to additionally differentiate structure based on year. For both methods, AMOVA results were similar, with microsatellite results explaining a smaller portion of the variation compared with the SNP results. The SNP-based markers presented here were able to successfully differentiate population structure similar to microsatellite results. These SNP markers represent new tools to discriminate B. cinerea isolates within closely related populations using multiple targeted sequences.

Development of a Diagnostic Assay for Race Differentiation of Podosphaera macularis

Hop powdery mildew (caused by Podosphaera macularis) was confirmed in the Pacific Northwest in 1996. Before 2012, the most common race of P. macularis was able to infect plants that possessed powdery mildew resistance based on the R-genes Rb, R3, and R5. After 2012, two additional races of P. macularis were discovered that can overcome the resistance gene R6 and the partial resistance found in the cultivar Cascade. These three races now occur throughout the region, which can complicate management and research efforts because of uncertainty on which race(s) may be present in the region and able to infect susceptible hop genotypes. Current methods for determining the races of P. macularis are labor intensive, costly, and typically require more than 14 days to obtain results. We sought to develop a molecular assay to differentiate races of the fungus possessing virulence on plants with R6, referred to as V6-virulent, from other races. The transcriptomes of 46 isolates of P. macularis were sequenced to identify loci and variants unique to V6 isolates. Fourteen primer pairs were designed for 10 candidate loci that contained single nucleotide polymorphisms (SNP) and short insertion-deletion polymorphisms. Two differentially labeled locked nucleic acid probes were designed for a contig that contained a conserved SNP associated with V6-virulence. The resulting duplexed real-time PCR assay was validated against 46 V6 and 54 non-V6 P. macularis isolates collected from the United States and Europe. The assay had perfect discrimination of V6-virulence among isolates of P. macularis originating from the western U.S. but failed to predict V6-virulence in three isolates collected from Europe. The specificity of the assay was tested with different species of powdery mildew fungi and other microorganisms associated with hop. Weak nonspecific amplification occurred with powdery mildew fungi collected from Vitis vinifera, Fragaria sp., and Zinnia sp.; however, nonspecification amplification is not a concern when differentiating pathogen race from colonies on hop. The assay has practical applications in hop breeding, epidemiological studies, and other settings where rapid confirmation of pathogen race is needed.

Population Genomic- and Phylogenomic-Enabled Advances to Increase Insight Into Pathogen Biology and Epidemiology

Population genetics has been a key discipline in phytopathology for many years. The recent rise in cost-effective, high-throughput DNA sequencing technologies, allows sequencing of dozens, if not hundreds of specimens, turning population genetics into population genomics and opening up new, exciting opportunities as described in this Focus Issue. Without the limitations of genetic markers and the availability of whole or near whole-genome data, population genomics can give new insights into the biology, evolution and adaptation, and dissemination patterns of plant-associated microbes.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

Anthropocene marks the era when human activity is making a significant impact on earth, its ecological and biogeographical systems. The domestication and intensification of agricultural and forest production systems have had a large impact on plant and tree health. Some pathogens benefitted from these human activities and have evolved and adapted in response to the expansion of crop and forest systems, resulting in global outbreaks. Global pathogen genomics data including population genomics and high-quality reference assemblies are crucial for understanding the evolution and adaptation of pathogens. Crops and forest trees have remarkably different characteristics, such as reproductive time and the level of domestication. They also have different production systems for disease management with more intensive management in crops than forest trees. By comparing and contrasting results from pathogen population genomic studies done on widely different agricultural and forest production systems, we can improve our understanding of pathogen evolution and adaptation to different selection pressures. We find that in spite of these differences, similar processes such as hybridization, host jumps, selection, specialization, and clonal expansion are shaping the pathogen populations in both crops and forest trees. We propose some solutions to reduce these impacts and lower the probability of global pathogen outbreaks so that we can envision better management strategies to sustain global food production as well as ecosystem services.