Spatial and Temporal Analysis of Populations of the Sudden Oak Death Pathogen in Oregon Forests

Bat white-nose disease fungus diversity in time and space

White-nose disease (WND), caused by the psychrophilic fungus Pseudogymnoascusdestructans, represents one of the greatest threats for North American hibernating bats. Research on molecular data has significantly advanced our knowledge of various aspects of the disease, yet more studies are needed regarding patterns of P.destructans genetic diversity distribution. In the present study, we investigate three sites within the native range of the fungus in detail: two natural hibernacula (karst caves) in Bulgaria, south-eastern Europe and one artificial hibernaculum (disused cellar) in Germany, northern Europe, where we conducted intensive surveys between 2014 and 2019. Using 18 microsatellite and two mating type markers, we describe how P.destructans genetic diversity is distributed between and within sites, the latter including differentiation across years and seasons of sampling; across sampling locations within the site; and between bats and hibernaculum walls. We found significant genetic differentiation between hibernacula, but we could not detect any significant differentiation within hibernacula, based on the variables examined. This indicates that most of the pathogen's movement occurs within sites. Genotypic richness of P.destructans varied between sites within the same order of magnitude, being approximately two times higher in the natural caves (Bulgaria) compared to the disused cellar (Germany). Within all sites, the pathogen's genotypic richness was higher in samples collected from hibernaculum walls than in samples collected from bats, which corresponds with the hypothesis that hibernacula walls represent the environmental reservoir of the fungus. Multiple pathogen genotypes were commonly isolated from a single bat (i.e. from the same swab sample) in all study sites, which might be important to consider when studying disease progression.

Openness and Computational Reproducibility in Plant Pathology: Where We Stand and a Way Forward

Open research practices have been highlighted extensively during the last 10 years in many fields of scientific study as essential standards needed to promote transparency and reproducibility of scientific results. Scientific claims can only be evaluated based on how protocols, materials, equipment, and methods were described; data were collected and prepared; and analyses were conducted. Openly sharing protocols, data, and computational code is central to current scholarly dissemination and communication, but in many fields, including plant pathology, adoption of these practices has been slow. We randomly selected 450 articles published from 2012 to 2021 across 21 journals representative of the plant pathology discipline and assigned them scores reflecting their openness and computational reproducibility. We found that most of the articles did not follow protocols for open science and failed to share data or code in a reproducible way. We propose that use of open-source tools facilitates computationally reproducible work and analyses, benefitting not just readers but the authors as well. Finally, we provide ideas and suggest tools to promote open, reproducible computational research practices among plant pathologists. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Cataloging Phytophthora Species of Agriculture, Forests, Horticulture, and Restoration Outplantings in California, U.S.A.: A Sequence-Based Meta-Analysis

California contains a diverse flora, and knowledge of the pathogens that threaten those plants is essential to managing their long-term health. To better understand threats to California plant health, a meta-analysis of Phytophthora detections within the state was conducted using publicly available sequences as a primary source of data rather than published records. Accessions of internal transcribed spacer (ITS) ribosomal DNA were cataloged from 800 Californian Phytophthora isolates, analyzed, and determined to correspond to 80 taxa, including several phylogenetically distinct provisional species. A number of Phytophthora taxa not previously reported from California were identified, including 20 described species. Pathways of introduction and spread were analyzed by categorizing isolates' origins, grouped by land-use: (i) agriculture, (ii) forests and other natural ecosystems, (iii) horticulture and nurseries, or (iv) restoration outplantings. The pooled Phytophthora metacommunities of the restoration outplantings and horticulture land-use categories were the most similar, whereas the communities pooled from forests and agriculture were least similar. Phytophthora cactorum, P. pini, P. pseudocryptogea, and P. syringae were identified in all four land-use categories, while 13 species were found in three. P. gonapodyides was the most common species by number of ITS accessions and exhibited the greatest diversity of ITS haplotypes. P. cactorum, P. ramorum, and P. nicotianae were associated with the greatest number of host genera. In this analysis, the Phytophthora spp. most prevalent in California differ from those compiled from the scientific literature.

Repeated Emergence of Sudden Oak Death in Oregon: Chronology, Impact, and Management

It has been two decades since the first detection of the sudden oak death pathogen Phytophthora ramorum in Oregon forests. Although the epidemic was managed since its first discovery in 2001, at least three invasions of three separate variants (clonal lineages), NA1, EU1, and NA2, are documented to have occurred to date. Control of this epidemic has cost over US$32 million from 2001 to 2020. This is dwarfed by the predicted cost of the closure to the Coos Bay export terminal, estimated at $58 million per year, if the epidemic was allowed to spread unchecked. Management efforts in Oregon have reduced inoculum and limited the spread of the pathogen. An outreach and citizen scientist program has been piloted to help in early detection efforts and search for disease-resistant tanoak. This feature article documents the repeated emergence, impact, costs, and lessons learned from managing this devastating invasive pathogen.

Local Eradication of Phytophthora ramorum Is Effective on Both NA1 and EU1 Lineages in Oregon Tanoak Forests

Sudden oak death (SOD), caused by the oomycete Phytophthora ramorum, has been actively managed in Oregon since its discovery there in 2001. SOD is a devastating disease affecting an ecologically and culturally important tree species in southwestern Oregon. Initially infested with the NA1 lineage, the more virulent EU1 lineage was discovered in 2015. Management has adapted over time in response to experimental findings and administrative limitations. Current management practices present an opportunity to compare the efficacy of treatment on these lineages by analyzing P. ramorum inoculum at untreated and treated sites. Current treatment includes herbicide treatment on host stems followed by felling, piling, and burning on site. Infested sites were visited between 2018 and 2020 (n = 88), where understory vegetation and soil was collected. Generalized linear modeling demonstrated that treatment had a significant impact on P. ramorum prevalence from vegetation samples, with an average of 33% (± 10%) fewer positive samples at treated sites. Linear mixed-effects modeling of a subpopulation of EU1 sites visited before and after treatment showed a similar effect of treatment, with a 43% (± 15%) reduction in P. ramorum prevalence. Prevalence of P. ramorum in soil was not affected by treatment in either analysis. A third analysis taking into consideration recent wildfire incursion into infested areas revealed that wildfire alone is insufficient to reduce prevalence of P. ramorum. These results strongly suggest that management is successfully reducing P. ramorum inoculum found on understory vegetation, and that treatment remains necessary to reduce the spread of this major forest pathogen.

Pathotype Complexity and Genetic Characterization of Phytophthora sojae Populations in Illinois, Indiana, Kentucky, and Ohio

Phytophthora sojae, the causal agent of Phytophthora root and stem rot of soybean, has been managed with single Rps genes since the 1960s but has subsequently adapted to many of these resistance genes, rendering them ineffective. The objective of this study was to examine the pathotype and genetic diversity of P. sojae from soil samples across Illinois, Indiana, Kentucky, and Ohio by assessing which Rps genes were still effective and identifying possible population clusters. There were 218 pathotypes identified from 473 P. sojae isolates with an average of 6.7 out of 15 differential soybean lines exhibiting a susceptible response for each isolate. Genetic characterization of 103 P. sojae isolates from across Illinois, Indiana, Kentucky, and Ohio with 19 simple sequence repeat markers identified 92 multilocus genotypes. There was a moderate level of population differentiation between these four states, with pairwise FST values ranging from 0.026 to 0.246. There were also moderate to high levels of differentiation between fields, with pairwise FST values ranging from 0.071 to 0.537. Additionally, cluster analysis detected the presence of P. sojae population structure across neighboring states. The level of pathotype and genetic diversity, in addition to the identification of population clusters, supports the hypothesis of occasional outcrossing events that allow an increase in diversity and the potential to select for a loss in avirulence to specific resistance genes within regions. The trend of suspected gene flow among neighboring fields is expected to be an ongoing issue with current agricultural practices.

First report of the NA2 clonal lineage of the sudden oak death pathogen, Phytophthora ramorum, infecting tanoak in Oregon forests

Phytophthora ramorum Werres, de Cock & Man in't Veld, causal agent of sudden oak death (SOD) and ramorum leaf blight, is comprised of four clonal lineages in its invasive ranges of North America and Europe (Grünwald et al. 2012, Van Poucke et al. 2012). Of these, three - the NA1, NA2, and EU1 lineages - are found in U.S. nurseries, but only two, the NA1 and EU1 lineages, have been found infecting trees in North American forests (Grünwald et al. 2012, 2016). In the spring of 2021, tanoak (Notholithocarpus densiflorus Manos, Cannon & Oh) displaying symptoms consistent with SOD were detected north of Port Orford (Curry County, Oregon). Symptoms were canopy dieback and blackened petiole and stem lesions on tanoak sprouts. The pathogen isolated on PAR (CMA plus 200 ml/L ampicillin, 10 mg/L rifamycin, 66.7 mg/L PCNB) selective media was determined to be P. ramorum based on characteristic morphology of hyphae, sporangia, and chlamydospores (Werres et al. 2001). Positive identification as P. ramorum was obtained with a lineage-specific LAMP assay targeting an NA2 orphan gene, indicating the presence of the NA2 lineage. NA2 was confirmed by sequencing a portion of the cellulose binding elicitor lectin (CBEL) gene using CBEL5U and CBEL6L primers (Gagnon et al. 2014). Sequences (GenBank accessions MZ733981 and MZ733982) were aligned against reference sequences for all lineages (Gagnon et al. 2014) confirming the presence of NA2. Lineage determination as NA2 was further confirmed at eleven SSR loci (ILVOPrMS145, PrMS39, PrMS9C3, ILVOPrMS79, KI18, KI64, PrMS45, PrMS6, ILVOPrMS131, KI82ab, and PrMS43) using the methods of Kamvar et al. (2015). We completed Koch's postulates using potted tanoaks, wound-inoculated at the midpoint of 1-year old stems with either hyphal plugs or non-colonized agar (n=4 per treatment). Tanoaks were maintained in a growth chamber (20°C-day / 18°C-night temperatures) with regular watering and an 18-photoperiod using F32T8 fluorescent bulbs (Phillips, Eindhoven, The Netherlands). After 7 days, brown to black lesions 1.2 to 2.9 cm in length were observed on the inoculated stems, from which P. ramorum was subsequently re-isolated; no symptoms were observed on the controls, and no pathogens were recovered when plating the wound sites in PAR. This is the first detection of the NA2 lineage causing disease in forests worldwide. The outbreak was found on private and public lands in forests typical to the SOD outbreak in Oregon (mixed conifer and tanoak), and was 33 km north of the closest known P. ramorum infestation. Follow-up ground surveys on adjacent lands have identified over 100 P. ramorum-positive tanoak trees, from which additional NA2 isolates have been recovered from bole cankers. NA2 is thought to be more aggressive than the NA1 lineage (Elliott et al. 2011), which has been present in Curry County since the mid-1990s (Goheen et al. 2017). Eradication of the NA2 lineage is being pursued to slow its further spread and prevent overlap with existing NA1 and EU1 populations. The repeated introductions of novel lineages into the western United States native plant communities highlights the vulnerability of this region to Phytophthora establishment, justifying continued monitoring for P. ramorum in nurseries and forests. References • Elliott, M, et al. 2011. For. Path. 41:7. https://doi.org/10.1111/j.1439-0329.2009.00627.x • Gagnon, M.-C., et al. 2014. Can. J. Plant Pathol. 36:367. https://doi.org/10.1080/07060661.2014.924999 • Goheen, E.M., et al. 2017. For. Phytophthoras 7:45. https://doi: 10.5399/osu/fp.7.1.4030 • Grünwald, N. J., et al. 2012. Trends Microbiol. 20:131. https://doi.org/10.1016/j.tim.2011.12.006 • Grünwald, N. J., et al. 2016. Plant Dis. 100:1024. https://doi.org/10.1094/PDIS-10-15-1169-PDN • Kamvar, Z.N. et al. 2015. Phytopath. 105:982. https://doi.org/10.1094/PHYTO-12-14-0350-FI • Van Poucke, K., et al. 2012. Fungal Biol. 116:1178. https://doi.org/10.1016/j.funbio.2012.09.003 • Werres, S., et al. 2001. Mycol. Res. 105: 1155. https://doi.org/10.1016/S0953-7562(08)61986-3.

Microsatellite Markers from Peronospora tabacina, the Cause of Blue Mold of Tobacco, Reveal Species Origin, Population Structure, and High Gene Flow

Peronospora tabacina is an obligate parasite that causes blue mold of tobacco. The pathogen reproduces primarily by sporangia, whereas the sexual oospores are rarely observed. A collection of 122 isolates of P. tabacina was genotyped using nine microsatellites to assess the population structure of individuals from subpopulations collected from central, southern, and western Europe; the Middle East; Central America; North America; and Australia. Genetic variations among the six subpopulations accounted for ∼8% of the total variation, including moderate levels of genetic differentiation, high gene flow among these subpopulations, and a positive correlation between geographic and genetic distance (r = 0.225; P < 0.001). Evidence of linkage disequilibrium (P < 0.001) showed that populations contained partially clonal subpopulations but that subpopulations from Australia and Mediterranean Europe did not. High genetic variation and population structure among samples could be explained by continuous gene flow across continents via infected transplant exchange and/or long-distance dispersal of sporangia via wind currents. This study analyzed the most numerous P. tabacina collection and allowed conclusions regarding the migration, mutation, and evolutionary history of this obligate biotrophic oomycete. The evidence pointed to the species origin in Australia and identified intracontinental and intercontinental migration patterns of this important pathogen.[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.

Molecular Characterization of Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) Collections from Nine Countries

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat worldwide. To understand the worldwide distribution of its molecular groups, as well as the diversity, differentiation, and migration of the Pst populations, 567 isolates collected from nine countries (China, Pakistan, Italy, Egypt, Ethiopia, Canada, Mexico, Ecuador, and the U.S.) in 2010–2018 were genotyped using 14 codominant simple sequence repeat markers. A total of 433, including 333 new multi-locus genotypes (MLGs), were identified, which were clustered into ten molecular groups (MGs). The MGs and country-wise populations differed in genetic diversity, heterozygosity, and correlation coefficient between the marker and virulence data. Many isolates from different countries, especially the isolates from Mexico, Ecuador, and the U.S., were found to be identical or closely related MLGs, and some of the MGs were present in all countries, indicating Pst migrations among different countries. The analysis of molecular variance revealed 78% variation among isolates, 12% variation among countries, and 10% variation within countries. Only low levels of differentiation were found by the pairwise comparisons of country populations. Of the 10 MGs, 5 were found to be involved in sexual and/or somatic recombination. Identical and closely related MLGs identified from different countries indicated international migrations. The study provides information on the distributions of various Pst genetic groups in different countries and evidence for the global migrations, which should be useful in understanding the pathogen evolution and in stressing the need for continual monitoring of the disease and pathogen populations at the global scale.

Novel Introductions and Epidemic Dynamics of the Sudden Oak Death Pathogen Phytophthora ramorum in Oregon Forests

Sudden oak death caused by Phytophthora ramorum has been actively managed in Oregon since the early 2000s. To date, this epidemic has been driven mostly by the NA1 clonal lineage of P. ramorum, but an outbreak of the EU1 lineage has recently emerged. Here, we contrast the population dynamics of the NA1 outbreak first reported in 2001 to the outbreak of the EU1 lineage first detected in 2015. We performed tests to determine whether any of the lineages were introduced more than once. Infested regions of the forest were sampled between 2013 and 2018 (n = 903), and strains were genotyped at 15 microsatellite loci. Most genotypes observed were transient, with 272 of 358 unique genotypes emerging during one year and disappearing the next year. The diversity of EU1 was very low and isolates were spatially clustered (less than 8 km apart), suggesting a single EU1 introduction. Some forest isolates are genetically similar to isolates collected from a local nursery in 2012, suggesting the introduction of EU1 from this nursery or simultaneous introduction to both the nursery and latently into the forest. In contrast, the older NA1 populations were more polymorphic and spread more than 30 km². A principal component analysis supported two to four independent NA1 introductions. The NA1 and EU1 epidemics infest the same area but show disparate demographics because of the initial introductions of the lineages spaced 10 years apart. Comparing these epidemics provides novel insight regarding patterns of emergence of clonal pathogens in forest ecosystems.

Hybridization and introgression of native and foreign Sorbus tree species in unique environments of protected mountainous areas

Hybridization and introgression are important processes influencing the genetic diversity and evolution of species. These processes are of particular importance in protected areas, where they can lead to the formation of hybrids between native and foreign species and may ultimately result in the loss of parental species from their natural range. Despite their importance, the contribution of hybridization and introgression to genetic diversity in Sorbus genus remains not fully recognized. We analysed the genetic and morphological variability of several Sorbus species including native (Sorbus aria), foreign (S. intermedia) and potentially hybrid (S. carpatica) individuals from the Polish Carpathian range. Patterns of variation at 13 nuclear microsatellite loci show hybridization between the tested species and confirm the existence of the hybrid form S. carpatica. Biometric analysis on leaves, based of 10 metric features and three parameters, identified several characters for preliminary taxonomic classification; however, none of them could be used as a fully diagnostic marker for faultless annotation of S. intermedia and S. carpatica. The genetic structure analysis indicated complex patterns of population differentiation and its diverse origin. The results allow assessment of genetic variation and identification of parental species participating in hybridization. This knowledge will advance the management of genetic diversity and development of conservation strategies for efficient maintenance of the unique protected ecosystem.

Molecular Phylogenomics and Population Structure of Phytophthora pluvialis

Phytophthora pluvialis is an oomycete that was first isolated from soil, water, and tree foliage in mixed Douglas-fir-tanoak forests of the U.S. Pacific Northwest (PNW). It was then identified as the causal agent of red needle cast of radiata pine (Pinus radiata) in New Zealand (NZ). Genotyping-by-sequencing was used to obtain 1,543 single nucleotide polymorphisms across 145 P. pluvialis isolates to characterize the population structure in the PNW and NZ. We tested the hypothesis that P. pluvialis was introduced to NZ from the PNW using genetic distance measurements and population structure analyses among locations between countries. The low genetic distance, population heterozygosity, and lack of geographic structure in NZ suggest a single colonization event from the United States followed by clonal expansion in NZ. The PNW Coast Range was proposed as a presumptive center of origin of the currently known distribution of P. pluvialis based on its geographic range and position as the central cluster in a minimum spanning network. The Coastal cluster of isolates were located at the root of every U.S. cluster and emerged earlier than all NZ clusters. The Coastal cluster had the highest degree of heterozygosity (Hs = 0.254) and median pairwise genetic distance (0.093) relative to any other cluster. Finally, the rapid host diversification between closely related isolates of P. pluvialis in NZ indicate that this pathogen has the potential to infect a broader range of hosts than is currently recognized.

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.

Genetic analysis reveals long-standing population differentiation and high diversity in the rust pathogen Melampsora lini

A priority for research on infectious disease is to understand how epidemiological and evolutionary processes interact to influence pathogen population dynamics and disease outcomes. However, little is understood about how population adaptation changes across time, how sexual vs. asexual reproduction contribute to the spread of pathogens in wild populations and how diversity measured with neutral and selectively important markers correlates across years. Here, we report results from a long-term study of epidemiological and genetic dynamics within several natural populations of the Linum marginale-Melampsora lini plant-pathogen interaction. Using pathogen isolates collected from three populations of wild flax (L. marginale) spanning 16 annual epidemics, we probe links between pathogen population dynamics, phenotypic variation for infectivity and genomic polymorphism. Pathogen genotyping was performed using 1567 genome-wide SNP loci and sequence data from two infectivity loci (AvrP123, AvrP4). Pathogen isolates were phenotyped for infectivity using a differential set. Patterns of epidemic development were assessed by conducting surveys of infection prevalence in one population (Kiandra) annually. Bayesian clustering analyses revealed host population and ecotype as key predictors of pathogen genetic structure. Despite strong fluctuations in pathogen population size and severe annual bottlenecks, analysis of molecular variance revealed that pathogen population differentiation was relatively stable over time. Annually, varying levels of clonal spread (0–44.8%) contributed to epidemics. However, within populations, temporal genetic composition was dynamic with rapid turnover of pathogen genotypes, despite the dominance of only four infectivity phenotypes across the entire study period. Furthermore, in the presence of strong fluctuations in population size and migration, spatial selection may maintain pathogen populations that, despite being phenotypically stable, are genetically highly dynamic.

Melampsora lini is a rust fungus that infects native flax, Linum marginale in south-eastern Australia where its epidemiology and evolution have been intensively studied since 1987. Over that time, substantial diversity in the pathotypic structure of M. lini has been demonstrated but an understanding of how genetic diversity in pathogen populations is maintained through space and time is lacking. Here we integrated phenotypic, genotypic and epidemiological datasets spanning 16 annual epidemics across three host populations to examine long-term pathogen genetic dynamics. The results show that host ecotype is the dominant selective force in the face of strong bottlenecks and annual patterns of genetic turnover. Results from previous studies indicate that in this geographic region, M. lini lacks the capacity to reproduce sexually–we thus expected to find limited genetic diversity and evidence for strong clonality influencing genetic dynamics within growing seasons. However, the breadth of genomic coverage provided by the SNP markers revealed high levels of genotypic variation within M. lini populations. This discovery contrasts with observed phenotypic dynamics as the epidemics of this pathogen were largely dominated by four pathotypes across the study period. Based on a detailed assessment and comparison of pathotypic and genotypic patterns, our study increases the understanding of how genetic diversity is generated and maintained through space and time within wild pathogen populations. The implications for the management of resistance to pathogens in agricultural or conservation contexts are significant: the appearance of clonality may be hiding high levels of pathogen diversity and recombination. Understanding how this diversity is generated could provide new and unique ways to mitigate or suppress the emergence of infectious strains, allowing to efficiently combat harmful diseases.

Tools for Evolutionary and Genetic Analysis (TEGA): A new platform for the management of molecular and environmental data

Population genetics studies the distributions and changes in population allele frequencies in response to processes, such as mutation, natural selection, gene flow, and genetic drift. Researchers daily manage genetic, biological, and environmental data of the samples, storing them in text files or spreadsheets, which makes it difficult to maintain consistency and traceability. Here we present TEGA, a WEB-based stand-alone software developed for the easy analysis and management of population genetics data. It was designed to: 1) facilitate data management, 2) provide a way to execute the analysis procedures, and 3) supply a means to publish data, procedures, and results. TEGA is distributed under the GNU AGPL v3 license. The documentation, source code, and screenshots are available at https://github.com/darioelias/TEGA. In addition, we present Rabid Fish, the first implementation of TEGA in the Genetics Labortory of the Faculty of Humanities and Sciences at the National University of the Litoral, where research focuses on population genetics studies applied to non-model organisms.

Genetic Diversity and Conservation Status of Helianthus verticillatus, an Endangered Sunflower of the Southern United States

Evaluating species diversity and patterns of population genetic variation is an essential aspect of conservation biology to determine appropriate management strategies and preserve the biodiversity of native plants. Habitat fragmentation and potential habitat loss are often an outcome of a reduction in naturally occurring wildfires and controlled prescribed burning, as seen in Helianthus verticillatus (whorled sunflower). This endangered, wild relative of the common sunflower, Helianthus annuus, is endemic to four locations in Alabama, Georgia, and Tennessee, United States. Despite its endangered status, there is no recovery plan for H. verticillatus, and knowledge related to its basic plant biology and importance in ecosystem services is mostly unknown. In this study, we utilized 14 microsatellite loci to investigate fine-scale population structure and genetic diversity of H. verticillatus individuals found on two sampling sites within the Georgia population. Our results indicated moderate genetic diversity and the presence of two distinct genetic clusters. Analyses of molecular variance indicated that the majority of variance was individually based, thus confirming high genetic differentiation and limited gene flow between H. verticillatus collection sites. The evidence of a population bottleneck in these sites suggests a recent reduction in population size that could be explained by habitat loss and population fragmentation. Also, high levels of linkage disequilibrium were detected, putatively suggesting clonal reproduction among these individuals. Our study provides a better understanding of fine-scale genetic diversity and spatial distribution of H. verticillatus populations in Georgia. Our results can underpin an original recovery plan for H. verticillatus that could be utilized for the conservation of this endangered species and to promote its persistence in the wild.

Genetic Differentiation and Clonal Expansion of Chinese Botrytis cinerea Populations from Tomato and Other Crops in China

Botrytis cinerea is an important pathogen of vegetable and fruit crops but little is known about its population structure and genetics in China. We hypothesized that the geographic populations of B. cinerea in China would be genetically differentiated by host, geographic location, and/or year. In this study, we collected 393 B. cinerea isolates representing 28 populations from tomato, cherry, and nectarine from 2006 to 2014 in China. The isolates were analyzed using 14 microsatellite markers, including six new markers that provided more genotyping power than the eight previously published loci. We also investigated the B. cinerea population structure and inferred its mode of reproduction and dispersal based on genotype data. High genotypic diversity was detected in all populations, and clonal reproduction was dominant. Southern China populations harbored more genotypes than northern populations. Differentiation by host plant was evident. Between 2011 and 2012, genotypes changed only slightly among years for Liaoning populations, but they changed substantially among years for the Shanghai and Fujian populations. Clonal dispersal was detected and the farthest dispersal distance was estimated to be about 1,717 km. Two high-frequency genotypes were widely distributed in more than 10 populations and across several years. Our results provide useful, novel information for plant breeding programs and control of B. cinerea in China.

Population Genomic Analyses Reveal Connectivity via Human-Mediated Transport across Populus Plantations in North America and an Undescribed Subpopulation of Sphaerulina musiva

Domestication of plant species has affected the evolutionary dynamics of plant pathogens in agriculture and forestry. A model system for studying the consequences of plant domestication on the evolution of an emergent plant disease is the fungal pathogen Sphaerulina musiva. This ascomycete causes leaf spot and stem canker disease of Populus spp. and their hybrids. A population genomics approach was used to determine the degree of population structure and evidence for selection on the North American population of S. musiva. In total, 122 samples of the fungus were genotyped identifying 120,016 single-nucleotide polymorphisms after quality filtering. In North America, S. musiva has low to moderate degrees of differentiation among locations. Three main genetic clusters were detected: southeastern United States, midwestern United States and Canada, and a new British Columbia cluster (BC2). Population genomics suggest that BC2 is a novel genetic cluster from central British Columbia, clearly differentiated from previously reported S. musiva from coastal British Columbia, and the product of a single migration event. Phenotypic measurements from greenhouse experiments indicate lower aggressiveness of BC2 on Populus trichocarpa. In summary, S. musiva has geographic structure across broad regions indicative of gene flow among clusters. The interconnectedness of the North American S. musiva populations across large geographic distances further supports the hypothesis of anthropogenic-facilitated transport of the pathogen.

Low Genetic Diversity Suggests the Recent Introduction of Dogwood Powdery Mildew to North America

Cornus florida (flowering dogwood) is a popular understory tree endemic to the eastern hardwood forests of the United States. In 1996, dogwood powdery mildew caused by Erysiphe pulchra, an obligate biotrophic fungus of large bracted dogwoods, reached epidemic levels throughout the C. florida growing region. In the late 1990s, both sexual and asexual stages of E. pulchra were regularly observed; thereafter, the sexual stage was found less frequently. We examined the genetic diversity and population structure of 167 E. pulchra samples on C. florida leaves using 15 microsatellite loci. Samples were organized into two separate collection zone data sets, separated as eight zones and two zones, for the subsequent analysis of microsatellite allele length data. Clone correction analysis reduced the sample size to 90 multilocus haplotypes. Our study indicated low genetic diversity, a lack of definitive population structure, low genetic distance among multilocus haplotypes, and significant linkage disequilibrium among zones. Evidence of a population bottleneck was also detected. The results of our study indicated a high probability that E. pulchra reproduces predominately via asexual conidia and lend support to the hypothesis that E. pulchra is an exotic pathogen to North America.[Formula: see text] Copyright © 2019 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Ecology and Evolution of the Sudden Oak Death Pathogen Phytophthora ramorum

The sudden oak and sudden larch death pathogen Phytophthora ramorum emerged simultaneously in the United States on oak and in Europe on Rhododendron in the 1990s. This pathogen has had a devastating impact on larch plantations in the United Kingdom as well as mixed conifer and oak forests in the Western United States. Since the discovery of this pathogen, a large body of research has provided novel insights into the emergence, epidemiology, and genetics of this pandemic. Genetic and genomic resources developed for P. ramorum have been instrumental in improving our understanding of the epidemiology, evolution, and ecology of this disease. The recent reemergence of EU1 in the United States and EU2 in Europe and the discovery of P. ramorum in Asia provide renewed impetus for research on the sudden oak death pathogen.

Film-Forming Polymers and Surfactants Reduce Infection and Sporulation of Phytophthora ramorum on Rhododendron

Phytophthora ramorum, cause of sudden oak death and ramorum leaf blight, can persist undetected in infested nurseries. Many conventional fungicides are effective in reducing or delaying symptom expression but some may confound visual detection of infected plants. We tested film-forming polymers (FFPs) and surfactants for their ability to reduce infection and sporulation of P. ramorum on rhododendron. FFPs (Anti-Stress, Moisturin, Nature Shield, Nu-Film, and Vapor Gard) and surfactants (Tergitol, Zonix, and an unregistered AGAE product) were screened in detached-leaf assays. Anti-Stress, Nu-Film, Zonix, and a Nu-Film-Zonix mixture were additionally tested for durability, protection against exposure to infested water, and a reduction in sporulation. FFP effectiveness was retained for at least 3 weeks of exposure to overhead irrigation and rain. Relative to controls, foliar treatments protected rhododendron branches exposed to infested water. No treatments prevented symptom development when applied postinfection but leaves treated with Anti-Stress, Zonix, and the Nu-Film-Zonix mixture produced significantly fewer sporangia relative to controls. Application of FFPs and surfactants to quarantined, potentially infected plants offers a management tool for reducing infection and sporulation but not symptom expression, thereby limiting disease spread without interfering with disease detection.

The Genetic Structure of Populations of the Douglas-Fir Swiss Needle Cast Fungus Nothophaeocryptopus gaeumannii in New Zealand

Swiss needle cast is a foliar disease of Douglas-fir (Pseudotsuga menziesii) that results in premature foliage loss and reduced growth. The causal fungus, Nothophaeocryptopus gaeumannii, was first detected in New Zealand in 1959 and spread throughout the North and South Islands over the following decades. The contemporary genetic structure of the N. gaeumannii population in New Zealand was assessed by analyzing 468 multilocus SSR genotypes (MLGs) from 2,085 N. gaeumannii isolates collected from 32 sites in the North and South Islands. Overall diversity was lower than that reported from native N. gaeumannii populations in the northwestern United States, which was expected given that N. gaeumannii is introduced in New Zealand. Linkage disequilibrium was significantly higher than expected under random mating, suggesting that population structure is clonal. Populations of N. gaeumannii in the North and South Islands were weakly differentiated, and the isolates collected from sites within the islands were moderately differentiated. This suggests that gene flow has occurred between the N. gaeumannii populations in the North and South Islands, and between the local N. gaeumannii populations within each island. Eighteen isolates of N. gaeumannii Lineage 2, which has previously been reported only from western Oregon, were recovered from two sites in the North Island and four sites in the South Island. The most likely explanation for the contemporary distribution of N. gaeumannii in New Zealand is that it was introduced on infected live seedlings through the forestry or ornamental nursery trade, as the fungus is neither seed borne nor saprobic, and the observed population structure is not consistent with a stochastic intercontinental dispersal event.

Contrasting fine-scale genetic structure of two sympatric clonal plants in an alpine swampy meadow featured by tussocks

Tussocks are unique vegetation structures in wetlands. Many tussock species mainly reproduce by clonal growth, resulting in genetically identical offspring distributed in various spatial patterns. These fine-scale patterns could influence mating patterns and thus the long-term evolution of wetland plants. Here, we contribute the first genetic and clonal structures of two key species in alpine wetlands on the Qinghai–Tibet Plateau, Kobresia tibetica and Blysmus sinocompressus, using > 5000 SNPs identified by 2b-RAD sequencing. The tussock-building species, K. tibetica, has a phalanx (clumping) growth form, but different genets could co-occur within the tussocks, indicating that it is not proper to treat a tussock as one genetic individual. Phalanx growth does not necessarily lead to increased inbreeding in K. tibetica. B. sinocompressus has a guerilla (spreading) growth form, with the largest detected clone size being 18.32 m, but genets at the local scale tend to be inbred offspring. Our results highlight that the combination of clone expansion and seedling recruitment facilitates the contemporary advantage of B. sinocompressus, but its evolutionary potential is limited by the input genetic load of the original genets. The tussocks of K. tibetica are more diverse and a valuable genetic legacy of former well-developed wet meadows, and they are worthy of conservation attention.

Populations of Phytophthora rubi Show Little Differentiation and High Rates of Migration Among States in the Western United States

Population genetics is a powerful tool to understand patterns and evolutionary processes that are involved in plant-pathogen emergence and adaptation to agricultural ecosystems. We are interested in studying the population dynamics of Phytophthora rubi, the causal agent of Phytophthora root rot in raspberry. P. rubi is found in the western United States, where most of the fresh and processed raspberries are produced. We used genotyping-by-sequencing to characterize genetic diversity in populations of P. rubi sampled in the United States and other countries. Our results confirm that P. rubi is a monophyletic species with complete lineage sorting from its sister taxon P. fragariae. Overall, populations of P. rubi show low genetic diversity across the western United States. Demographic analyses suggest that populations of P. rubi from the western United States are the source of pathogen migration to Europe. We found no evidence for population differentiation at a global or regional (western United States) level. Finally, our results provide evidence of migration from California and Oregon into Washington. This report provides new insights into the evolution and structure of global and western United States populations of the raspberry pathogen P. rubi, indicating that human activity might be involved in moving the pathogen among regions and fields.

Best Practices for Population Genetic Analyses

Population genetic analysis is a powerful tool to understand how pathogens emerge and adapt. However, determining the genetic structure of populations requires complex knowledge on a range of subtle skills that are often not explicitly stated in book chapters or review articles on population genetics. What is a good sampling strategy? How many isolates should I sample? How do I include positive and negative controls in my molecular assays? What marker system should I use? This review will attempt to address many of these practical questions that are often not readily answered from reading books or reviews on the topic, but emerge from discussions with colleagues and from practical experience. A further complication for microbial or pathogen populations is the frequent observation of clonality or partial clonality. Clonality invariably makes analyses of population data difficult because many assumptions underlying the theory from which analysis methods were derived are often violated. This review provides practical guidance on how to navigate through the complex web of data analyses of pathogens that may violate typical population genetics assumptions. We also provide resources and examples for analysis in the R programming environment.

Confirmation of independent introductions of an exotic plant pathogen of Cornus species, Discula destructiva, on the east and west coasts of North America

Cornus florida (flowering dogwood) and C. nuttallii (Pacific dogwood) are North American native tree species that belong to the big-bracted group of dogwoods. Cornus species are highly valued for their ornamental characteristics, and have fruits that contain high fat content for animals. Also, they are an important understory tree in natural forests. Dogwood anthracnose, caused by Discula destructiva, was observed in the late 1970s on the east and west coasts of the United States and by 1991 had quickly spread throughout most of the native ranges of C. florida and C. nuttalli. We investigated the genetic diversity and population structure of 93 D. destructiva isolates using 47 microsatellite loci developed from the sequenced genome of the type strain of D. destructiva. Clone-corrected data indicated low genetic diversity and the presence of four genetic clusters that corresponded to two major geographic areas, the eastern United States and the Pacific Northwest, and to the two collection time periods when the isolates were collected (pre- and post-1993). Linkage disequilibrium was present in five out of six subpopulations, suggesting that the fungus only reproduced asexually. Evidence of population bottlenecks was indicated across four identified genetic clusters, and was probably the result of the limited number of founding individuals on both coasts. These results support the hypothesis that D. destructiva is an exotic pathogen with independent introductions on the east and west coasts of North America. We also tested the cross-amplification of these microsatellite primers to other Discula species. Genomic DNA from 17 isolates of four other Discula species and two isolates of Juglanconis species (formerly Melanconis species) were amplified by 17 of 47 primer pairs. These primers may be useful for investigating the genetic diversity and population structure of these Discula species.

Genome sequences of 12 isolates of the EU1 lineage of Phytophthora ramorum, a fungus-like pathogen that causes extensive damage and mortality to a wide range of trees and other plants

Here we present genome sequences for twelve isolates of the invasive pathogen Phytophthora ramorum EU1. The assembled genome sequences and raw sequence data are available via BioProject accession number PRJNA177509. These data will be useful in developing molecular tools for specific detection and identification of this pathogen.

Ecology and pathology of Phytophthora ITS clade 3 species in forests in western Oregon, USA

Phytophthora species are widespread and diverse in forest ecosystems, but little is known about their ecology. We explore ecological attributes of the closely related clade 3 species that occur sympatrically in western North American forests. We address the population structure, pathology, and epidemiology of P. ilicis, P. nemorosa, P. pluvialis, P. pseudosyringae, and P. psychrophila. Phytophthora species were isolated from plant tissues, rainwater falling through the forest canopy, streams, and soils in forests in western Oregon. Species identifications were based on morphology in culture with molecular confirmation using COX spacer and internal transcribed spacer (ITS) sequences. All five clade 3 Phytophthora species are present in western Oregon forests, although P. ilicis (only 1 forest isolate) and P. psychrophila (only 12 isolates) are apparently rare. P. ilicis is known only from holly in horticultural situations and once from a naturalized seedling in an urban forest. The known distribution of P. nemorosa in forest settings coincides with the ranges of its principle hosts, tanoak and myrtlewood, in Oregon and California. Although it is regularly identified from streams within the tanoak range, it has not been recovered from streams beyond that range. P. pluvialis is primarily associated with Douglas-fir canopies. It was identified from scattered locations throughout western Oregon in rain traps beneath Douglas-fir plantations and from diseased needles. P. pseudosyringae is also isolated from tanoak and myrtlewood in southwest Oregon and California, but its distribution, in streams at least, extends throughout much of western Oregon. P. psychrophila in Oregon is known only from rain traps beneath tanoak trees. Little intraspecific variation was detected in the nuclear rDNA ITS of clade 3 species. Variation in the mitochondrial COX spacer region was more frequent, with 2 to 10 haplotypes identified in the clade 3 species, for which we had multiple isolates.

Genomic Analyses of Dominant U.S. Clonal Lineages of Phytophthora infestans Reveals a Shared Common Ancestry for Clonal Lineages US11 and US18 and a Lack of Recently Shared Ancestry Among All Other U.S. Lineages

Populations of the potato and tomato late-blight pathogen Phytophthora infestans are well known for emerging as novel clonal lineages. These successions of dominant clones have historically been named US1 through US24, in order of appearance, since their first characterization using molecular markers. Hypothetically, these lineages can emerge through divergence from other U.S. lineages, recombination among lineages, or as novel, independent lineages originating outside the United States. We tested for the presence of phylogenetic relationships among U.S. lineages using a population of 31 whole-genome sequences, including dominant U.S. clonal lineages as well as available samples from global populations. We analyzed ancestry of the whole mitochondrial genome and samples of nuclear loci, including supercontigs 1.1 and 1.5 as well as several previously characterized coding regions. We found support for a shared ancestry among lineages US11 and US18 from the mitochondrial genome as well as from one nuclear haplotype on each supercontig analyzed. The other nuclear haplotype from each sample assorted independently, indicating an independent ancestry. We found no support for emergence of any other of the U.S. lineages from a common ancestor shared with the other U.S. lineages. Each of the U.S. clonal lineages fit a model where populations of new clonal lineages emerge via migration from a source population that is sexual in nature and potentially located in central Mexico or elsewhere. This work provides novel insights into patterns of emergence of clonal lineages in plant pathogen genomes.

Developing educational resources for population genetics in R: an open and collaborative approach

The r computing and statistical language community has developed a myriad of resources for conducting population genetic analyses. However, resources for learning how to carry out population genetic analyses in r are scattered and often incomplete, which can make acquiring this skill unnecessarily difficult and time consuming. To address this gap, we developed an online community resource with guidance and working demonstrations for conducting population genetic analyses in r. The resource is freely available at http://popgen.nescent.org and includes material for both novices and advanced users of r for population genetics. To facilitate continued maintenance and growth of this resource, we developed a toolchain, process and conventions designed to (i) minimize financial and labour costs of upkeep; (ii) to provide a low barrier to contribution; and (iii) to ensure strong quality assurance. The toolchain includes automatic integration testing of every change and rebuilding of the website when new vignettes or edits are accepted. The process and conventions largely follow a common, distributed version control-based contribution workflow, which is used to provide and manage open peer review by designated website editors. The online resources include detailed documentation of this process, including video tutorials. We invite the community of population geneticists working in r to contribute to this resource, whether for a new use case of their own, or as one of the vignettes from the 'wish list' we maintain, or by improving existing vignettes.

Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality

To gain a detailed understanding of how plant microbes evolve and adapt to hosts, pesticides, and other factors, knowledge of the population dynamics and evolutionary history of populations is crucial. Plant pathogen populations are often clonal or partially clonal which requires different analytical tools. With the advent of high throughput sequencing technologies, obtaining genome-wide population genetic data has become easier than ever before. We previously contributed the R package poppr specifically addressing issues with analysis of clonal populations. In this paper we provide several significant extensions to poppr with a focus on large, genome-wide SNP data. Specifically, we provide several new functionalities including the new function mlg.filter to define clone boundaries allowing for inspection and definition of what is a clonal lineage, minimum spanning networks with reticulation, a sliding-window analysis of the index of association, modular bootstrapping of any genetic distance, and analyses across any level of hierarchies.