Genomic Characterization of a South American Phytophthora Hybrid Mandates Reassessment of the Geographic Origins of Phytophthora infestans

Developing a crop- wild-reservoir pathogen system to understand pathogen evolution and emergence

Crop pathogens reduce yield and contribute to global malnourishment. Surveillance not only detects presence/absence but also reveals genetic diversity, which can inform our understanding of rapid adaptation and control measures. An often neglected aspect is that pathogens may also use crop wild relatives as alternative hosts. This study develops the beet (Beta vulgaris) rust (Uromyces beticola) system to explore how crop pathogens evolve to evade resistance using a wild reservoir. We test predictions that crop selection will drive virulence gene differentiation and affect rates of sex between crop- and wild-host rust populations. We sequenced, assembled, and annotated the 588 Mb beet rust genome, developed a novel leaf peel pathogen DNA extraction protocol, and analysed genetic diversity in 42 wild and crop isolates. We found evidence for two populations: one containing exclusively wild-host isolates; the other containing all crop-host isolates, plus five wild isolates. Effectors showed greater diversity in the exclusively wild population and greater differentiation between populations. Preliminary evidence suggests the rates of sexual reproduction may differ between populations. This study highlights how differences in pathogen populations might be used to identify genes important for survival on crops and how reproduction might impact adaptation. These findings are relevant to all crop-reservoir systems and will remain unnoticed without comparison to wild reservoirs.

A pangenome analysis reveals the center of origin and evolutionary history of Phytophthora infestans and 1c clade species

We examined the evolutionary history of Phytophthora infestans and its close relatives in the 1c clade. We used whole genome sequence data from 69 isolates of Phytophthora species in the 1c clade and conducted a range of genomic analyses including nucleotide diversity evaluation, maximum likelihood trees, network assessment, time to most recent common ancestor and migration analysis. We consistently identified distinct and later divergence of the two Mexican Phytophthora species, P. mirabilis and P. ipomoeae, from P. infestans and other 1c clade species. Phytophthora infestans exhibited more recent divergence from other 1c clade species of Phytophthora from South America, P. andina and P. betacei. Speciation in the 1c clade and evolution of P. infestans occurred in the Andes. P. andina-P. betacei-P. infestans formed a species complex with indistinct species boundaries, hybridizations between the species, and short times to common ancestry. Furthermore, the distinction between modern Mexican and South American P. infestans proved less discrete, suggesting gene flow between populations over time. Admixture analysis indicated a complex relationship among these populations, hinting at potential gene flow across these regions. Historic P. infestans, collected from 1845-1889, were the first to diverge from all other P. infestans populations. Modern South American populations diverged next followed by Mexican populations which showed later ancestry. Both populations were derived from historic P. infestans. Based on the time of divergence of P. infestans from its closest relatives, P. andina and P. betacei in the Andean region, we consider the Andes to be the center of origin of P. infestans, with modern globalization contributing to admixture between P. infestans populations today from Mexico, the Andes and Europe.

Reconstructing the Global Migration History of Phytophthora infestans Toward Colombia

The evolution of new variants of plant pathogens is one of the biggest challenges to controlling and managing plant diseases. Of the forces driving these evolutionary processes, global migration events are particularly important for widely distributed diseases such as potato late blight, caused by the oomycete Phytophthora infestans. However, little is known about its migration routes outside North America and Europe. This work used genotypic data from population studies to elucidate the migration history originating the Colombian P. infestans population. For this purpose, a dataset of 1,706 P. infestans genotypes was recollected, representing North and South America, Europe, and Asia. Descriptive analysis and historical records from North America and Europe were used to propose three global migration hypotheses, differing on the origin of the disease (Mexico or Peru) and the hypothesis that it returned to South America from Europe. These scenarios were tested using approximate Bayesian computation. According to this analysis, the most probable scenario (posterior probability = 0.631) was the one proposing a Peruvian origin for P. infestans, an initial migration toward Colombia and Mexico, and a later event from Mexico to the United States and then to Europe and Asia, with no return to northern South America. In Colombia, the scenario considering a single migration from Peru and posterior migrations within Colombia was the most probable, with a posterior probability of 0.640. The obtained results support the hypothesis of a Peruvian origin for P. infestans followed by rare colonization events worldwide.

Evolution of Phytophthora infestans on its potato host since the Irish potato famine

Phytophthora infestans is a major oomycete plant pathogen, responsible for potato late blight, which led to the Irish Potato Famine from 1845-1852. Since then, potatoes resistant to this disease have been bred and deployed worldwide. Their resistance (R) genes recognize pathogen effectors responsible for virulence and then induce a plant response stopping disease progression. However, most deployed R genes are quickly overcome by the pathogen. We use targeted sequencing of effector and R genes on herbarium specimens to examine the joint evolution in both P. infestans and potato from 1845-1954. Currently relevant effectors are historically present in P. infestans, but with alternative alleles compared to modern reference genomes. The historic FAM-1 lineage has the virulent Avr1 allele and the ability to break the R1 resistance gene before breeders deployed it in potato. The FAM-1 lineage is diploid, but later, triploid US-1 lineages appear. We show that pathogen virulence genes and host resistance genes have undergone significant changes since the Famine, from both natural and artificial selection.

Reconstructing historic and modern potato late blight outbreaks using text analytics

In 1843, a hitherto unknown plant pathogen entered the US and spread to potato fields in the northeast. By 1845, the pathogen had reached Ireland leading to devastating famine. Questions arose immediately about the source of the outbreaks and how the disease should be managed. The pathogen, now known as Phytophthora infestans, still continues to threaten food security globally. A wealth of untapped knowledge exists in both archival and modern documents, but is not readily available because the details are hidden in descriptive text. In this work, we (1) used text analytics of unstructured historical reports (1843-1845) to map US late blight outbreaks; (2) characterized theories on the source of the pathogen and remedies for control; and (3) created modern late blight intensity maps using Twitter feeds. The disease spread from 5 to 17 states and provinces in the US and Canada between 1843 and 1845. Crop losses, Andean sources of the pathogen, possible causes and potential treatments were discussed. Modern disease discussion on Twitter included near-global coverage and local disease observations. Topic modeling revealed general disease information, published research, and outbreak locations. The tools described will help researchers explore and map unstructured text to track and visualize pandemics.

Plant NLR immunity activation and execution: a biochemical perspective

Plants deploy cell-surface and intracellular receptors to detect pathogen attack and trigger innate immune responses. Inside host cells, families of nucleotide-binding/leucine-rich repeat (NLR) proteins serve as pathogen sensors or downstream mediators of immune defence outputs and cell death, which prevent disease. Established genetic underpinnings of NLR-mediated immunity revealed various strategies plants adopt to combat rapidly evolving microbial pathogens. The molecular mechanisms of NLR activation and signal transmission to components controlling immunity execution were less clear. Here, we review recent protein structural and biochemical insights to plant NLR sensor and signalling functions. When put together, the data show how different NLR families, whether sensors or signal transducers, converge on nucleotide-based second messengers and cellular calcium to confer immunity. Although pathogen-activated NLRs in plants engage plant-specific machineries to promote defence, comparisons with mammalian NLR immune receptor counterparts highlight some shared working principles for NLR immunity across kingdoms.

Exploring the Emergence and Evolution of Plant Pathogenic Microbes Using Historical and Paleontological Sources

Biotechnological advances now permit broad exploration of past microbial communities preserved in diverse substrates. Despite biomolecular degradation, high-throughput sequencing of preserved materials can yield invaluable genomic and metagenomic data from the past. This line of research has expanded from its initial human- and animal-centric foci to include plant-associated microbes (viruses, archaea, bacteria, fungi, and oomycetes), for which historical, archaeological, and paleontological data illuminate past epidemics and evolutionary history. Genetic mechanisms underlying the acquisition of microbial pathogenicity, including hybridization, polyploidization, and horizontal gene transfer, can now be reconstructed, as can gene-for-gene coevolution with plant hosts. Epidemiological parameters, such as geographic origin and range expansion, can also be assessed. Building on published case studies with individual phytomicrobial taxa, the stage is now set for broader, community-wide studies of preserved plant microbiomes to strengthen mechanistic understanding of microbial interactions and plant disease emergence.

Comparative analyses of complete Peronosporaceae (Oomycota) mitogenome sequences - insights into structural evolution and phylogeny

Members of the Peronosporaceae (Oomycota, Chromista), which currently consists of 25 genera and approximately 1000 recognised species, are responsible for disease on a wide range of plant hosts. Molecular phylogenetic analyses over the last two decades have improved our understanding of evolutionary relationships within Peronosporaceae. To date, 16 numbered and three named clades have been recognised; it is clear from these studies that the current taxonomy does not reflect evolutionary relationships. Whole organelle genome sequences are an increasingly important source of phylogenetic information, and in this study we present comparative and phylogenetic analyses of mitogenome sequences from 15 of the 19 currently recognized clades of Peronosporaceae, including 44 newly assembled sequences. Our analyses suggest strong conservation of mitogenome size and gene content across Peronosporaceae but, as previous studies have suggested, limited conservation of synteny. Specifically, we identified 28 distinct syntenies amongst the 71 examined isolates. Moreover, 19 of the isolates contained inverted or direct repeats, suggesting repeated sequences may be more common than previously thought. In terms of phylogenetic relationships, our analyses of 34 concatenated mitochondrial gene sequences resulted in a topology that was broadly consistent with previous studies. However, unlike previous studies concatenated mitochondrial sequences provided strong support for higher level relationships within the family.

Global historic pandemics caused by the FAM-1 genotype of Phytophthora infestans on six continents

The FAM-1 genotype of Phytophthora infestans caused late blight in the 1840s in the US and Europe and was responsible for the Irish famine. We sampled 140 herbarium specimens collected between 1845 and 1991 from six continents and used 12-plex microsatellite genotyping (SSR) to identify FAM-1 and the mtDNA lineage (Herb-1/Ia) present in historic samples. FAM-1 was detected in approximately 73% of the historic specimens and was found on six continents. The US-1 genotype was found later than FAM-1 on all continents except Australia/Oceania and in only 27% of the samples. FAM-1 was the first genotype detected in almost all the former British colonies from which samples were available. The data from historic outbreak samples suggest the FAM-1 genotype was widespread, diverse, and spread to Asia and Africa from European sources. The famine lineage spread to six continents over 144 years, remained widespread and likely spread during global colonization from Europe. In contrast, modern lineages of P. infestans are rapidly displaced and sexual recombination occurs in some regions.

The persistent threat of emerging plant disease pandemics to global food security

Plant disease outbreaks are increasing and threaten food security for the vulnerable in many areas of the world. Now a global human pandemic is threatening the health of millions on our planet. A stable, nutritious food supply will be needed to lift people out of poverty and improve health outcomes. Plant diseases, both endemic and recently emerging, are spreading and exacerbated by climate change, transmission with global food trade networks, pathogen spillover, and evolution of new pathogen lineages. In order to tackle these grand challenges, a new set of tools that include disease surveillance and improved detection technologies including pathogen sensors and predictive modeling and data analytics are needed to prevent future outbreaks. Herein, we describe an integrated research agenda that could help mitigate future plant disease pandemics.

The mitogenome of Phytophthora agathidicida: Evidence for a not so recent arrival of the "kauri killing" Phytophthora in New Zealand

Phytophthora agathidicida is associated with a root rot that threatens the long-term survival of the iconic New Zealand kauri. Although it is widely assumed that this pathogen arrived in New Zealand post-1945, this hypothesis has yet to be formally tested. Here we describe evolutionary analyses aimed at evaluating this and two alternative hypotheses. As a basis for our analyses, we assembled complete mitochondrial genome sequences from 16 accessions representing the geographic range of P. agathidicida as well as those of five other members of Phytophthora clade 5. All 21 mitogenome sequences were very similar, differing little in size with all sharing the same gene content and arrangement. We first examined the temporal origins of genetic diversity using a pair of calibration schemes. Both resulted in similar age estimates; specifically, a mean age of 303.0-304.4 years and 95% HPDs of 206.9-414.6 years for the most recent common ancestor of the included isolates. We then used phylogenetic tree building and network analyses to investigate the geographic distribution of the genetic diversity. Four geographically distinct genetic groups were recognised within P. agathidicida. Taken together the inferred age and geographic distribution of the sampled mitogenome diversity suggests that this pathogen diversified following arrival in New Zealand several hundred to several thousand years ago. This conclusion is consistent with the emergence of kauri dieback disease being a consequence of recent changes in the relationship between the pathogen, host, and environment rather than a post-1945 introduction of the causal pathogen into New Zealand.

Population Structure of Phytophthora infestans from a Single Location in Poland Over a Long Period of Time in Context of Weather Conditions

Phytophthora infestans (Mont.) de Bary is a destructive potato pathogen. Changing weather conditions are among the factors that influence the pathogen population structure. In this study, 237 P. infestans isolates were collected from a single unprotected experimental field in an area with high late-blight pressure located in Boguchwała in the southeastern part of Poland during 15 growing seasons (2000-2014). The isolates were assessed for mating type, mitochondrial haplotype, resistance to metalaxyl, virulence, and polymorphism of 14 single-sequence repeat markers (SSRs). The results revealed 89 unique genotypes among the 237 P. infestans isolates. Eighty-seven isolates belonged to genotype 34_A1, which was detected in all the years of research except 2012. Isolates of P. infestans from individual years were very similar to each other, as shown by Nei's genetic identity based on 14 SSR markers. The obtained results on isolate characteristics were analyzed in terms of meteorological data (air temperature and precipitation) and indicated that frost, long winters, and hot, dry summers did not directly affect the P. infestans population structure. We described the variability in metalaxyl resistance and virulence among isolates of the P. infestans genotype 34_A1.

Comparative Performance of Popular Methods for Hybrid Detection using Genomic Data

Interspecific hybridization is an important evolutionary phenomenon that generates genetic variability in a population and fosters species diversity in nature. The availability of large genome scale datasets has revolutionized hybridization studies to shift from the observation of the presence or absence of hybrids to the investigation of the genomic constitution of hybrids and their genome-specific evolutionary dynamics. Although a handful of methods have been proposed in an attempt to identify hybrids, accurate detection of hybridization from genomic data remains a challenging task. In addition to methods that infer phylogenetic networks or that utilize pairwise divergence, site pattern frequency based and population genetic clustering approaches are popularly used in practice, though the performance of these methods under different hybridization scenarios has not been extensively examined. Here, we use simulated data to comparatively evaluate the performance of four tools that are commonly used to infer hybridization events: the site pattern frequency based methods HyDe and the D-statistic (i.e., the ABBA-BABA test) and the population clustering approaches structure and ADMIXTURE. We consider single hybridization scenarios that vary in the time of hybridization and the amount of incomplete lineage sorting (ILS) for different proportions of parental contributions (γ); introgressive hybridization; multiple hybridization scenarios; and a mixture of ancestral and recent hybridization scenarios. We focus on the statistical power to detect hybridization and the false discovery rate (FDR) for comparisons of the D-statistic and HyDe, and the accuracy of the estimates of γ as measured by the mean squared error for HyDe, structure, and ADMIXTURE. Both HyDe and the D-statistic are powerful for detecting hybridization in all scenarios except those with high ILS, although the D-statistic often has an unacceptably high FDR. The estimates of γ in HyDe are impressively robust and accurate whereas structure and ADMIXTURE sometimes fail to identify hybrids, particularly when the proportional parental contributions are asymmetric (i.e., when γ is close to 0). Moreover, the posterior distribution estimated using structure exhibits multimodality in many scenarios, making interpretation difficult. Our results provide guidance in selecting appropriate methods for identifying hybrid populations from genomic data.

Genome Analysis of Two Newly Emerged Potato Late Blight Isolates Sheds Light on Pathogen Adaptation and Provides Tools for Disease Management

Phytophthora infestans, the causal agent of the Irish Potato Famine in the 1840s, is one of the most destructive crop pathogens that threaten global food security. Host resistance (R) genes may help to control the disease, but recognition by through the gene products can be evaded by newly emerging isolates. Such isolates are dangerous as they may cause disease outbreaks under favorable conditions. However, our lack of knowledge about adaptation in these isolates jeopardizes an apt response to resistance breakdown. Here we performed genome and transcriptome sequencing of HB1501 and HN1602, two field isolates from distinct Chinese geographic regions. We found extensive polymorphisms in these isolates, including gene copy number variations, nucleotide polymorphisms, and gene expression changes. Effector encoding genes, which contribute to virulence, show distinct expression landscapes in P. infestans isolates HB1501 and HN1602. In particular, polymorphisms at multiple effectors required for recognition (Avr loci) enabled these isolates to overcome corresponding R gene based resistance. Although the isolates evolved multiple strategies to evade recognition, we experimentally verified that several R genes such as R8, RB, and Rpi-vnt1.1 remain effective against these isolates and are valuable to potato breeding in the future. In summary, rapid characterization of the adaptation in emerging field isolates through genomic tools inform rational agricultural management to prevent potential future epidemics.

Variola virus genome sequenced from an eighteenth-century museum specimen supports the recent origin of smallpox

Smallpox, caused by the variola virus (VARV), was a highly virulent disease with high mortality rates causing a major threat for global human health until its successful eradication in 1980. Despite previously published historic and modern VARV genomes, its past dissemination and diversity remain debated. To understand the evolutionary history of VARV with respect to historic and modern VARV genetic variation in Europe, we sequenced a VARV genome from a well-described eighteenth-century case from England (specimen P328). In our phylogenetic analysis, the new genome falls between the modern strains and another historic strain from Lithuania, supporting previous claims of larger diversity in early modern Europe compared to the twentieth century. Our analyses also resolve a previous controversy regarding the common ancestor between modern and historic strains by confirming a later date around the seventeenth century. Overall, our results point to the benefit of historic genomes for better resolution of past VARV diversity and highlight the value of such historic genomes from around the world to further understand the evolutionary history of smallpox as well as related diseases. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

The Recovery, Interpretation and Use of Ancient Pathogen Genomes

The ability to sequence genomes from ancient biological material has provided a rich source of information for evolutionary biology and engaged considerable public interest. Although most studies of ancient genomes have focused on vertebrates, particularly archaic humans, newer technologies allow the capture of microbial pathogens and microbiomes from ancient and historical human and non-human remains. This coming of age has been made possible by techniques that allow the preferential capture and amplification of discrete genomes from a background of predominantly host and environmental DNA. There are now near-complete ancient genome sequences for three pathogens of considerable historical interest - pre-modern bubonic plague (Yersinia pestis), smallpox (Variola virus) and cholera (Vibrio cholerae) - and for three equally important endemic human disease agents - Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy) and Treponema pallidum pallidum (syphilis). Genomic data from these pathogens have extended earlier work by paleopathologists. There have been efforts to sequence the genomes of additional ancient pathogens, with the potential to broaden our understanding of the infectious disease burden common to past populations from the Bronze Age to the early 20th century. In this review we describe the state-of-the-art of this rapidly developing field, highlight the contributions of ancient pathogen genomics to multidisciplinary endeavors and describe some of the limitations in resolving questions about the emergence and long-term evolution of pathogens.

Genome-Wide Increased Copy Number is Associated with Emergence of Dominant Clones of the Irish Potato Famine Pathogen Phytophthora infestans

The plant pathogen that caused the Irish potato famine, Phytophthora infestans, continues to reemerge globally. These modern epidemics are caused by clonally reproducing lineages. In contrast, a sexual mode of reproduction is observed at its center of origin in Mexico. We conducted a comparative genomic analysis of 47 high-coverage genomes to infer changes in genic copy number. We included samples from sexual populations at the center of origin as well as several dominant clonal lineages sampled worldwide. We conclude that sexual populations at the center of origin are diploid, as was the lineage that caused the famine, while modern clonal lineages showed increased copy number (3×). Copy number variation (CNV) was found genome-wide and did not to adhere to the two-speed genome hypothesis. Although previously reported, tetraploidy was not found in any of the genomes evaluated. We propose a model of dominant clone emergence supported by the epidemiological record (e.g., EU_13_A2, US-11, US-23) whereby a higher copy number provides fitness, leading to replacement of prior clonal lineages.IMPORTANCE The plant pathogen implicated in the Irish potato famine, Phytophthora infestans, continues to reemerge globally. Understanding changes in the genome during emergence can provide insights useful for managing this pathogen. Previous work has relied on studying individuals from the United States, South America, Europe, and China reporting that these can occur as diploids, triploids, or tetraploids and are clonal. We studied variation in sexual populations at the pathogen's center of origin, in Mexico, where it has been reported to reproduce sexually as well as within clonally reproducing, dominant clones from the United States and Europe. Our results newly show that sexual populations at the center of origin are diploid, whereas populations elsewhere are more variable and show genome-wide variation in gene copy number. We propose a model of evolution whereby new pathogen clones emerge predominantly by increasing the gene copy number genome-wide.

Ancient Plant Genomics in Archaeology, Herbaria, and the Environment

The ancient DNA revolution of the past 35 years has driven an explosion in the breadth, nuance, and diversity of questions that are approachable using ancient biomolecules, and plant research has been a constant, indispensable facet of these developments. Using archaeological, paleontological, and herbarium plant tissues, researchers have probed plant domestication and dispersal, plant evolution and ecology, paleoenvironmental composition and dynamics, and other topics across related disciplines. Here, we review the development of the ancient DNA discipline and the role of plant research in its progress and refinement. We summarize our understanding of long-term plant DNA preservation and the characteristics of degraded DNA. In addition, we discuss challenges in ancient DNA recovery and analysis and the laboratory and bioinformatic strategies used to mitigate them. Finally, we review recent applications of ancient plant genomic research.

Analysis of the lineage of Phytophthora infestans isolates using mating type assay, traditional markers, and next generation sequencing technologies

Phytophthora infestans (Mont.) de Bary, a hemibiotrophic oomycete, has caused severe epidemics of late blight in tomato and potato crops around the world since the Irish Potato Famine in the 1840s. Breeding of late blight resistant cultivars is one of the most effective strategies to overcome this disruptive disease. However, P. infestans is able to break down host resistance and acquire resistance to various fungicides, possibly because of the existence of high genetic variability among P. infestans isolates via sexual and asexual reproduction. Therefore, to manage this disease, it is important to understand the genetic divergence of P. infestans isolates. In this study, we analyzed the genomes of P. infestans isolates collected from Egypt and Japan using various molecular approaches including the mating type assay and genotyping simple sequence repeats, mitochondria DNA, and effector genes. We also analyzed genome-wide single nucleotide polymorphisms using double-digest restriction-site associated DNA sequencing and whole genome resequencing (WGRS). The isolates were classified adequately using high-resolution genome-wide approaches. Moreover, these analyses revealed new clusters of P. infestans isolates in the Egyptian population. Monitoring the genetic divergence of P. infestans isolates as well as breeding of resistant cultivars would facilitate the elimination of the late blight disease.

Genetic Structure and Subclonal Variation of Extant and Recent U.S. Lineages of Phytophthora infestans

The oomycete Phytophthora infestans is an important plant pathogen on potato and tomato crops. We examined the genetic structure of extant 20th and 21st century U.S. lineages of P. infestans and compared them with populations from South America and Mexico to examine genetic relationships and potential sources of lineages. US-23, currently the most prevalent lineage detected in the United States, shared genetic similarity primarily with the BR-1 lineage identified in the 1990s from Bolivia and Brazil. Lineages US-8, US-14, and US-24, predominantly virulent on potato, formed a cluster distinct from other U.S. lineages. Many of the other U.S. lineages shared significant genetic similarity with Mexican populations. The US-1 lineage, dominant in the mid-20th century, clustered with US-1 lineages from Peru. A survey of the presence of RXLR effector PiAVR2 revealed that some lineages carried PiAVR2, its resistance-breaking variant PiAVR2-like, or both. Minimum spanning networks developed from simple sequence repeat genotype datasets from USABlight outbreaks clearly showed the expansion of US-23 over a 6-year time period and geographic substructuring of some lineages in the western United States. Many clonal lineages of P. infestans in the United States have come from introductions from Mexico, but the US-23 and US-1 lineages were most likely introduced from other sources.

Pathogen enrichment sequencing (PenSeq) enables population genomic studies in oomycetes

The oomycete pathogens Phytophthora infestans and P. capsici cause significant crop losses world-wide, threatening food security. In each case, pathogenicity factors, called RXLR effectors, contribute to virulence. Some RXLRs are perceived by resistance proteins to trigger host immunity, but our understanding of the demographic processes and adaptive evolution of pathogen virulence remains poor. Here, we describe PenSeq, a highly efficient enrichment sequencing approach for genes encoding pathogenicity determinants which, as shown for the infamous potato blight pathogen Phytophthora infestans, make up < 1% of the entire genome. PenSeq facilitates the characterization of allelic diversity in pathogen effectors, enabling evolutionary and population genomic analyses of Phytophthora species. Furthermore, PenSeq enables the massively parallel identification of presence/absence variations and sequence polymorphisms in key pathogen genes, which is a prerequisite for the efficient deployment of host resistance genes. PenSeq represents a cost-effective alternative to whole-genome sequencing and addresses crucial limitations of current plant pathogen population studies, which are often based on selectively neutral markers and consequently have limited utility in the analysis of adaptive evolution. The approach can be adapted to diverse microbes and pathogens.

Insights into evolving global populations of Phytophthora infestans via new complementary mtDNA haplotype markers and nuclear SSRs

In many parts of the world the damaging potato late blight pathogen, Phytophthora infestans, is spread as a succession of clonal lineages. The discrimination of genetic diversity within such evolving populations provides insights into the processes generating novel lineages and the pathways and drivers of pathogen evolution and dissemination at local and global scales. This knowledge, in turn, helps optimise management practices. Here we combine two key methods for dissecting mitochondrial and nuclear diversity and resolve intra and inter-lineage diversity of over 100 P. infestans isolates representative of key clonal lineages found globally. A novel set of PCR primers that amplify five target regions are provided for mitochondrial DNA sequence analysis. These five loci increased the number of mtDNA haplotypes resolved from four with the PCR RFLP method to 37 (17, 6, 8 and 4 for Ia, Ib, IIa, and IIb haplotypes, respectively, plus 2 Herb-1 haplotypes). As with the PCR RFLP method, two main lineages, I and II were defined. Group I contained 25 mtDNA haplotypes that grouped broadly according to the Ia and Ib types and resolved several sub-clades amongst the global sample. Group II comprised two distinct clusters with four haplotypes corresponding to the RFLP type IIb and eight haplotypes resolved within type IIa. The 12-plex SSR assay revealed 90 multilocus genotypes providing accurate discrimination of dominant clonal lineages and other genetically diverse isolates. Some association of genetic diversity and geographic region of contemporary isolates was observed; US and Mexican isolates formed a loose grouping, distinct from isolates from Europe, South America and other regions. Diversity within clonal lineages was observed that varied according to the age of the clone. In combination, these fine-scale nuclear and maternally inherited mitochondrial markers enabled a greater level of discrimination among isolates than previously available and provided complementary perspectives on evolutionary questions relating to the diversity, phylogeography and the origins and spread of clonal lineages of P. infestans.

Phytophthora betacei, a new species within Phytophthora clade 1c causing late blight on Solanum betaceum in Colombia

Over the past few years, symptoms akin to late blight disease have been reported on a variety of crop plants in South America. Despite the economic importance of these crops, the causal agents of the diseases belonging to the genus Phytophthora have not been completely characterized. In this study, a new Phytophthora species was described in Colombia from tree tomato (Solanum betaceum), a semi-domesticated fruit grown in northern South America. Comprehensive phylogenetic, morphological, population genetic analyses, and infection assays to characterize this new species, were conducted. All data support the description of the new species, Phytophthora betacei sp. nov. Phylogenetic analyses suggest that this new species belongs to clade 1c of the genus Phytophthora and is a close relative of the potato late blight pathogen, P. infestans. Furthermore, it appeared as the sister group of the P. andina strains collected from wild Solanaceae (clonal lineage EC-2). Analyses of morphological and physiological characters as well as host specificity showed high support for the differentiation of these species. Based on these results, a complete description of the new species is provided and the species boundaries within Phytophthora clade 1c in northern South America are discussed.

Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing

Ancient DNA research has been revolutionized following development of next-generation sequencing platforms. Although a number of such platforms have been applied to ancient DNA samples, the Illumina series are the dominant choice today, mainly because of high production capacities and short read production. Recently a potentially attractive alternative platform for palaeogenomic data generation has been developed, the BGISEQ-500, whose sequence output are comparable with the Illumina series. In this study, we modified the standard BGISEQ-500 library preparation specifically for use on degraded DNA, then directly compared the sequencing performance and data quality of the BGISEQ-500 to the Illumina HiSeq2500 platform on DNA extracted from 8 historic and ancient dog and wolf samples. The data generated were largely comparable between sequencing platforms, with no statistically significant difference observed for parameters including level (P = 0.371) and average sequence length (P = 0718) of endogenous nuclear DNA, sequence GC content (P = 0.311), double-stranded DNA damage rate (v. 0.309), and sequence clonality (P = 0.093). Small significant differences were found in single-strand DNA damage rate (δS; slightly lower for the BGISEQ-500, P = 0.011) and the background rate of difference from the reference genome (θ; slightly higher for BGISEQ-500, P = 0.012). This may result from the differences in amplification cycles used to polymerase chain reaction–amplify the libraries. A significant difference was also observed in the mitochondrial DNA percentages recovered (P = 0.018), although we believe this is likely a stochastic effect relating to the extremely low levels of mitochondria that were sequenced from 3 of the samples with overall very low levels of endogenous DNA. Although we acknowledge that our analyses were limited to animal material, our observations suggest that the BGISEQ-500 holds the potential to represent a valid and potentially valuable alternative platform for palaeogenomic data generation that is worthy of future exploration by those interested in the sequencing and analysis of degraded DNA.

Coexistence of Multiple Endemic and Pandemic Lineages of the Rice Blast Pathogen

The rice blast fungus Magnaporthe oryzae (syn., Pyricularia oryzae) is both a threat to global food security and a model for plant pathology. Molecular pathologists need an accurate understanding of the origins and line of descent of M. oryzae populations in order to identify the genetic and functional bases of pathogen adaptation and to guide the development of more effective control strategies. We used a whole-genome sequence analysis of samples from different times and places to infer details about the genetic makeup of M. oryzae from a global collection of isolates. Analyses of population structure identified six lineages within M. oryzae, including two pandemic on japonica and indica rice, respectively, and four lineages with more restricted distributions. Tip-dating calibration indicated that M. oryzae lineages separated about a millennium ago, long after the initial domestication of rice. The major lineage endemic to continental Southeast Asia displayed signatures of sexual recombination and evidence of DNA acquisition from multiple lineages. Tests for weak natural selection revealed that the pandemic spread of clonal lineages entailed an evolutionary “cost,” in terms of the accumulation of deleterious mutations. Our findings reveal the coexistence of multiple endemic and pandemic lineages with contrasting population and genetic characteristics within a widely distributed pathogen.

The rice blast fungus Magnaporthe oryzae (syn., Pyricularia oryzae) is a textbook example of a rapidly adapting pathogen, and it is responsible for one of the most damaging diseases of rice. Improvements in our understanding of Magnaporthe oryzae’s diversity and evolution are required to guide the development of more effective control strategies. We used genome sequencing data for samples from around the world to infer the evolutionary history of M. oryzae. We found that M. oryzae diversified about 1,000 years ago, separating into six main lineages: two pandemic on japonica and indica rice, respectively, and four with more restricted distributions. We also found that a lineage endemic to continental Southeast Asia displayed signatures of sexual recombination and the acquisition of genetic material from multiple lineages. This work provides a population-level genomic framework for defining molecular markers for the control of rice blast and investigations of the molecular basis of differences in pathogenicity between M. oryzae lineages.

Large sub-clonal variation in Phytophthora infestans from recent severe late blight epidemics in India

The population structure of the Phytophthora infestans populations that caused the recent 2013–14 late blight epidemic in eastern India (EI) and northeastern India (NEI) was examined. The data provide new baseline information for populations of P. infestans in India. A migrant European 13_A2 genotype was responsible for the 2013–14 epidemic, replacing the existing populations. Mutations have generated substantial sub-clonal variation with 24 multi-locus genotypes (MLGs) found, of which 19 were unique variants not yet reported elsewhere globally. Samples from West Bengal were the most diverse and grouped alongside MLGs found in Europe, the UK and from neighbouring Bangladesh but were not linked directly to most samples from south India. The pathogen population was broadly more aggressive on potato than on tomato and resistant to the fungicide metalaxyl. Pathogen population diversity was higher in regions around the international borders with Bangladesh and Nepal. Overall, the multiple shared MLGs suggested genetic contributions from UK and Europe in addition to a sub-structure based on the geographical location within India. Our data indicate the need for improved phytosanitary procedures and continuous surveillance to prevent the further introduction of aggressive lineages of P. infestans into the country.

Gene flow analysis method, the D-statistic, is robust in a wide parameter space

Background

We evaluated the sensitivity of the D-statistic, a parsimony-like method widely used to detect gene flow between closely related species. This method has been applied to a variety of taxa with a wide range of divergence times. However, its parameter space and thus its applicability to a wide taxonomic range has not been systematically studied. Divergence time, population size, time of gene flow, distance of outgroup and number of loci were examined in a sensitivity analysis.

Result

The sensitivity study shows that the primary determinant of the D-statistic is the relative population size, i.e. the population size scaled by the number of generations since divergence. This is consistent with the fact that the main confounding factor in gene flow detection is incomplete lineage sorting by diluting the signal. The sensitivity of the D-statistic is also affected by the direction of gene flow, size and number of loci. In addition, we examined the ability of the f-statistics, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\widehat{f}}_G $$\end{document}f^G and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\widehat{f}}_{hom} $$\end{document}f^hom, to estimate the fraction of a genome affected by gene flow; while these statistics are difficult to implement to practical questions in biology due to lack of knowledge of when the gene flow happened, they can be used to compare datasets with identical or similar demographic background.

Conclusions

The D-statistic, as a method to detect gene flow, is robust against a wide range of genetic distances (divergence times) but it is sensitive to population size. The D-statistic should only be applied with critical reservation to taxa where population sizes are large relative to branch lengths in generations.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-2002-4) contains supplementary material, which is available to authorized users.

Signatures of selection and host-adapted gene expression of the Phytophthora infestans RNA silencing suppressor PSR2

Phytophthora infestans is a devastating pathogen in agricultural systems. Recently, an RNA silencing suppressor (PSR2, 'Phytophthora suppressor of RNA silencing 2') has been described in P. infestans. PSR2 has been shown to increase the virulence of Phytophthora pathogens on their hosts. This gene is one of the few effectors present in many economically important Phytophthora species. In this study, we investigated: (i) the evolutionary history of PSR2 within and between species of Phytophthora; and (ii) the interaction between sequence variation, gene expression and virulence. In P. infestans, the highest PiPSR2 expression was correlated with decreased symptom expression. The highest gene expression was observed in the biotrophic phase of the pathogen, suggesting that PSR2 is important during early infection. Protein sequence conservation was negatively correlated with host range, suggesting host range as a driver of PSR2 evolution. Within species, we detected elevated amino acid variation, as observed for other effectors; however, the frequency spectrum of the mutations was inconsistent with strong balancing selection. This evolutionary pattern may be related to the conservation of the host target(s) of PSR2 and the absence of known corresponding R genes. In summary, our study indicates that PSR2 is a conserved effector that acts as a master switch to modify plant gene regulation early during infection for the pathogen's benefit. The conservation of PSR2 and its important role in virulence make it a promising target for pathogen management.

Historic Late Blight Outbreaks Caused by a Widespread Dominant Lineage of Phytophthora infestans (Mont.) de Bary

Phytophthora infestans (Mont.) de Bary, the causal agent of potato late blight, was responsible for the Irish potato famine of the 1840s. Initial disease outbreaks occurred in the US in 1843, two years prior to European outbreaks. We examined the evolutionary relationships and source of the 19th-century outbreaks using herbarium specimens of P. infestans from historic (1846-1970) and more recent isolates (1992-2014) of the pathogen. The same unique SSR multilocus genotype, named here as FAM-1, caused widespread outbreaks in both US and Europe. The FAM-1 lineage shared allelic diversity and grouped with the oldest specimens collected in Colombia and Central America. The FAM-1 lineage of P. infestans formed a genetic group that was distinct from more recent aggressive lineages found in the US. The US-1 lineage formed a second, mid-20th century group. Recent modern US lineages and the oldest Mexican lineages formed a genetic group with recent Mexican lineages, suggesting a Mexican origin of recent US lineages. A survey of mitochondrial haplotypes in a larger set of global herbarium specimens documented the more frequent occurrence of the HERB-1 (type Ia) mitochondrial haplotype in archival collections from 1866-75 and 1906-1915 and the rise of the Ib mitochondrial lineage (US-1) between 1946-1955. The FAM-1 SSR lineage survived for almost 100 years in the US, was geographically widespread, and was displaced first in the mid-20th century by the US-1 lineage and then by distinct new aggressive lineages that migrated from Mexico.

Population Genomics of Fungal and Oomycete Pathogens

We are entering a new era in plant pathology in which whole-genome sequences of many individuals of a pathogen species are becoming readily available. Population genomics aims to discover genetic mechanisms underlying phenotypes associated with adaptive traits such as pathogenicity, virulence, fungicide resistance, and host specialization, as genome sequences or large numbers of single nucleotide polymorphisms become readily available from multiple individuals of the same species. This emerging field encompasses detailed genetic analyses of natural populations, comparative genomic analyses of closely related species, identification of genes under selection, and linkage analyses involving association studies in natural populations or segregating populations resulting from crosses. The era of pathogen population genomics will provide new opportunities and challenges, requiring new computational and analytical tools. This review focuses on conceptual and methodological issues as well as the approaches to answering questions in population genomics. The major steps start with defining relevant biological and evolutionary questions, followed by sampling, genotyping, and phenotyping, and ending in analytical methods and interpretations. We provide examples of recent applications of population genomics to fungal and oomycete plant pathogens.

Detecting hybridization using ancient DNA

It is well established that related species hybridize and that this can have varied but significant effects on speciation and environmental adaptation. It should therefore come as no surprise that hybridization is not limited to species that are alive today. In the last several decades, advances in technologies for recovering and sequencing DNA from fossil remains have enabled the assembly of high-coverage genome sequences for a growing diversity of organisms, including many that are extinct. Thanks to the development of new statistical approaches for detecting and quantifying admixture from genomic data, genomes from extinct populations have proven useful both in revealing previously unknown hybridization events and informing the study of hybridization between living organisms. Here, we review some of the key recent statistical innovations for detecting ancient hybridization using genomewide sequence data and discuss how these innovations have revised our understanding of human evolutionary history.