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.

Ancient DNA genomics and the renaissance of herbaria

Herbaria are undergoing a renaissance as valuable sources of genomic data for exploring plant evolution, ecology, and diversity. Ancient DNA retrieved from herbarium specimens can provide unprecedented glimpses into past plant communities, their interactions with biotic and abiotic factors, and the genetic changes that have occurred over time. Here, we highlight recent advances in the field of herbarium genomics and discuss the challenges and opportunities of combining data from modern and time-stamped historical specimens. We also describe how integrating herbarium genomics data with other data types can yield substantial insights into the evolutionary and ecological processes that shape plant communities. Herbarium genomic analysis is a tool for understanding plant life and informing conservation efforts in the face of dire environmental challenges.

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.

Genotypic and Phenotypic Structure of the Population of Phytophthora infestans in Egypt Revealed the Presence of European Genotypes

Late blight disease of potato and tomato, caused by Phytophthora infestans, results in serious losses to Egyptian and global potato and tomato production. To understand the structure and dynamics of the Egyptian population of P. infestans, 205 isolates were collected from potato and tomato plants during three growing seasons in 2010–2012. The characterization was achieved by mating-type assay, metalaxyl sensitivity assay, and virulence pattern. Additionally, genotyping of 85 Egyptian isolates and 15 reference UK isolates was performed using 12 highly informative microsatellite (SSR) markers David E. L. Cooke and five effector (RxLR) genes. Mating-type testing showed that 58% (118 of 205) of the isolates belonged to mating type A1, 35% (71 isolates) to mating type A2, and the rest 8% (16 isolates) were self-fertile. The phenotype of metalaxyl response was represented as 45% resistant, 43% sensitive, and 12% as intermediate. Structure analysis grouped the 85 identified genotypes into two main clonal lineages. The first clonal lineage comprised 21 isolates belonging to A2 mating type and 8 self-fertile isolates. This clonal lineage was identified as Blue_13 or EU_13_A2. The second main clonal lineage comprised 55 isolates and was identified as EU_23_A1. A single isolate with a novel SSR genotype that formed a distinct genetic grouping was also identified. The effector sequencing showed good correspondence with the virulence data and highlighted differences in the presence and absence of loci as well as nucleotide polymorphism that affect gene function. This study indicated a changing population of P. infestans in Egypt and discusses the findings in the context of late blight management.

Variance of allele balance calculated from low coverage sequencing data infers departure from a diploid state

Background

Polyploidy and heterokaryosis are common and consequential genetic phenomena that increase the number of haplotypes in an organism and complicate whole-genome sequence analysis. Allele balance has been used to infer polyploidy and heterokaryosis in diverse organisms using read sets sequenced to greater than 50× whole-genome coverage. However, sequencing to adequate depth is costly if applied to multiple individuals or large genomes.

Results

We developed VCFvariance.pl to utilize the variance of allele balance to infer polyploidy and/or heterokaryosis at low sequence coverage. This analysis requires as little as 10× whole-genome coverage and reduces the allele balance profile down to a single value, which can be used to determine if an individual has two or more haplotypes. This approach was validated using simulated, synthetic, and authentic read sets from the oomycete species Bremia lactucae and Phytophthora infestans, the fungal species Saccharomyces cerevisiae, and the plant species Arabidopsis arenosa. This approach was deployed to determine that nine of 21 genotyped European race-type isolates of Bremia lactucae were inconsistent with diploidy and therefore likely heterokaryotic.

Conclusions

Variance of allele balance is a reliable metric to detect departures from a diploid state, including polyploidy, heterokaryosis, a mixed sample, or chromosomal copy number variation. Deploying this strategy is computationally inexpensive, can reduce the cost of sequencing by up to 80%, and used to test any organism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04685-z.

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.

Quantitative resistance differences between and within natural populations of Solanum chilense against the oomycete pathogen Phytophthora infestans

The wild tomato species Solanum chilense is divided into geographically and genetically distinct populations that show signs of defense gene selection and differential phenotypes when challenged with several phytopathogens, including the oomycete causal agent of late blight Phytophthora infestans. To better understand the phenotypic diversity of this disease resistance in S. chilense and to assess the effect of plant genotype versus pathogen isolate, respectively, we evaluated infection frequency in a systematic approach and with large sample sizes. We studied 85 genetically distinct individuals representing nine geographically separated populations of S. chilense. This showed that differences in quantitative resistance can be observed between but also within populations at the level of individual plants. Our data also did not reveal complete immunity in any of the genotypes. We further evaluated the resistance of a subset of the plants against P. infestans isolates with diverse virulence properties. This confirmed that the relative differences in resistance phenotypes between individuals were mainly determined by the plant genotype under consideration with modest effects of pathogen isolate used in the study. Thus, our report suggests that the observed quantitative resistance against P. infestans in natural populations of a wild tomato species S. chilense is the result of basal defense responses that depend on the host genotype and are pathogen isolate-unspecific.

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.

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.

Mitochondrial genome sequence of Phytophthora sansomeana and comparative analysis of Phytophthora mitochondrial genomes

Phytophthora sansomeana infects soybean and causes root rot. It was recently separated from the species complex P. megasperma sensu lato. In this study, we sequenced and annotated its complete mitochondrial genome and compared it to that of nine other Phytophthora species. The genome was assembled into a circular molecule of 39,618 bp with a 22.03% G+C content. Forty-two protein coding genes, 25 tRNA genes and two rRNA genes were annotated in this genome. The protein coding genes include 14 genes in the respiratory complexes, four ATP synthase genes, 16 ribosomal proteins genes, a tatC translocase gene, six conserved ORFs and a unique orf402. The tRNA genes encode tRNAs for 19 amino acids. Comparison among mitochondrial genomes of 10 Phytophthora species revealed three inversions, each covering multiple genes. These genomes were conserved in gene content with few exceptions. A 3' truncated atp9 gene was found in P. nicotianae. All 10 Phytophthora species, as well as other oomycetes and stramenopiles, lacked tRNA genes for threonine in their mitochondria. Phylogenomic analysis using the mitochondrial genomes supported or enhanced previous findings of the phylogeny of Phytophthora spp.

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.

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.

Using herbaria to study global environmental change

During the last centuries, humans have transformed global ecosystems. With their temporal dimension, herbaria provide the otherwise scarce long-term data crucial for tracking ecological and evolutionary changes over this period of intense global change. The sheer size of herbaria, together with their increasing digitization and the possibility of sequencing DNA from the preserved plant material, makes them invaluable resources for understanding ecological and evolutionary species' responses to global environmental change. Following the chronology of global change, we highlight how herbaria can inform about long-term effects on plants of at least four of the main drivers of global change: pollution, habitat change, climate change and invasive species. We summarize how herbarium specimens so far have been used in global change research, discuss future opportunities and challenges posed by the nature of these data, and advocate for an intensified use of these 'windows into the past' for global change research and beyond.

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.

MtDNA allows the sensitive detection and haplotyping of the crayfish plague disease agent Aphanomyces astaci showing clues about its origin and migration

The oomycete Aphanomyces astaci, the causative agent of crayfish plague, is listed as one of the 100 worst invasive species in the world, destroying the native crayfish populations throughout Eurasia. The aim of this study was to examine the potential of selected mitochondrial (mt) genes to track the diversity of the crayfish plague pathogen A. astaci. Two sets of primers were developed to amplify the mtDNA of ribosomal rnnS and rnnL subunits. We confirmed two main lineages, with four different haplogroups and five haplotypes among 27 studied A. astaci strains. The haplogroups detected were (1) the A-haplogroup with the a-haplotype strains originating from Orconectes sp., Pacifastacus leniusculus and Astacus astacus; (2) the B-haplogroup with the b-haplotype strains originating from the P. leniusculus; (3) the D-haplogroup with the d1 and d2-haplotypes strains originating from Procambarus clarkii; and (4) the E-haplogroup with the e-haplotype strains originating from the Orconectes limosus. The described markers are stable and reliable and the results are easily repeatable in different laboratories. The present method has high applicability as it allows the detection and characterization of the A. astaci haplotype in acute disease outbreaks in the wild, directly from the infected crayfish tissue samples.

Reinforcing plant evolutionary genomics using ancient DNA

Improved understanding of ancient DNA (aDNA) biochemical properties coupled with application of next generation sequencing (NGS) methods enabled sequencing and authenticating genomes of historical samples. This advancement ignited a revolution in plant evolutionary genomics by allowing direct observations of past molecular diversity. Analyses of genomes sequenced from temporally distributed samples of Gossypium sp., Phytophthora infestans and Arabidopsis thaliana improved our understanding of the evolutionary rates and time scales at which genome remodeling takes place. Comparison of historical samples of barley (Hordeum vulgare) and maize (Zea mays ssp. mays) with their present-day counterparts enabled assessment of selection during different stages of domestication. These examples show how aDNA already improved our evolutionary inferences. Increasing quality and amount of sequencing data retrieved from historical plants will further advance our understanding of plant evolution.

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.

Mitochondrial genomes and comparative genomics of Aphanomyces astaci and Aphanomyces invadans

The genus Aphanomyces (Saprolegniales, Oomycetes) includes species with a variety of ecologies from saprotrophs to plant and animal parasites. Two important species in this genus are A. astaci, the cause of crayfish plague and its close relative, A. invadans, which causes the epizootic ulcerative syndrome on fish. In this study, we have assembled and annotated the mitochondrial (mt) genomes of A. astaci and A. invadans from the whole genome shotgun sequence reads (PRJNA187372; PRJNA258292, respectively). The assembly was generated from A. astaci Pc-genotype strain APO3 and A. invadans strain NJM9701. The sizes of the mtDNAs were 49,489 bp and 49,061 bp for A. astaci and A. invadans, respectively. The species shared similar genetic content and organization encoding 35 proteins, two ribosomal RNAs, three putative open reading frames and 33 transfer RNAs of 19 amino acids for peptide synthesis. Both species also had a large inverted repeat region (LIR) of approximately 12 kb, the LIR contained large and small ribosomal RNAs and eight protein coding genes. These annotated mt genomes serve as a valuable genetic backbone for further development of diagnostic methods and phylogenetic and migration studies of the animal parasitic species of Aphanomyces.

Genetic Variation within Clonal Lineages of Phytophthora infestans Revealed through Genotyping-By-Sequencing, and Implications for Late Blight Epidemiology

Genotyping-by-sequencing (GBS) was performed on 257 Phytophthora infestans isolates belonging to four clonal lineages to study within-lineage diversity. The four lineages used in the study were US-8 (n = 28), US-11 (n = 27), US-23 (n = 166), and US-24 (n = 36), with isolates originating from 23 of the United States and Ontario, Canada. The majority of isolates were collected between 2010 and 2014 (94%), with the remaining isolates collected from 1994 to 2009, and 2015. Between 3,774 and 5,070 single-nucleotide polymorphisms (SNPs) were identified within each lineage and were used to investigate relationships among individuals. K-means hierarchical clustering revealed three clusters within lineage US-23, with US-23 isolates clustering more by collection year than by geographic origin. K-means hierarchical clustering did not reveal significant clustering within the smaller US-8, US-11, and US-24 data sets. Neighbor-joining (NJ) trees were also constructed for each lineage. All four NJ trees revealed evidence for pathogen dispersal and overwintering within regions, as well as long-distance pathogen transport across regions. In the US-23 NJ tree, grouping by year was more prominent than grouping by region, which indicates the importance of long-distance pathogen transport as a source of initial late blight inoculum. Our results support previous studies that found significant genetic diversity within clonal lineages of P. infestans and show that GBS offers sufficiently high resolution to detect sub-structuring within clonal populations.

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.

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.

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

As the oomycete pathogen causing potato late blight disease, Phytophthora infestans triggered the famous 19th-century Irish potato famine and remains the leading cause of global commercial potato crop destruction. But the geographic origin of the genotype that caused this devastating initial outbreak remains disputed, as does the New World center of origin of the species itself. Both Mexico and South America have been proposed, generating considerable controversy. Here, we readdress the pathogen's origins using a genomic data set encompassing 71 globally sourced modern and historical samples of P. infestans and the hybrid species P. andina, a close relative known only from the Andean highlands. Previous studies have suggested that the nuclear DNA lineage behind the initial outbreaks in Europe in 1845 is now extinct. Analysis of P. andina's phased haplotypes recovered eight haploid genome sequences, four of which represent a previously unknown basal lineage of P. infestans closely related to the famine-era lineage. Our analyses further reveal that clonal lineages of both P. andina and historical P. infestans diverged earlier than modern Mexican lineages, casting doubt on recent claims of a Mexican center of origin. Finally, we use haplotype phasing to demonstrate that basal branches of the clade comprising Mexican samples are occupied by clonal isolates collected from wild Solanum hosts, suggesting that modern Mexican P. infestans diversified on Solanum tuberosum after a host jump from a wild species and that the origins of P. infestans are more complex than was previously thought.

Mitochondrial genome sequences reveal evolutionary relationships of the Phytophthora 1c clade species

Phytophthora infestans is one of the most destructive plant pathogens of potato and tomato globally. The pathogen is closely related to four other Phytophthora species in the 1c clade including P. phaseoli, P. ipomoeae, P. mirabilis and P. andina that are important pathogens of other wild and domesticated hosts. P. andina is an interspecific hybrid between P. infestans and an unknown Phytophthora species. We have sequenced mitochondrial genomes of the sister species of P. infestans and examined the evolutionary relationships within the clade. Phylogenetic analysis indicates that the P. phaseoli mitochondrial lineage is basal within the clade. P. mirabilis and P. ipomoeae are sister lineages and share a common ancestor with the Ic mitochondrial lineage of P. andina. These lineages in turn are sister to the P. infestans and P. andina Ia mitochondrial lineages. The P. andina Ic lineage diverged much earlier than the P. andina Ia mitochondrial lineage and P. infestans. The presence of two mitochondrial lineages in P. andina supports the hybrid nature of this species. The ancestral state of the P. andina Ic lineage in the tree and its occurrence only in the Andean regions of Ecuador, Colombia and Peru suggests that the origin of this species hybrid in nature may occur there.

Genome-enabled analysis of plant-pathogen migration

Trade in plant and plant products has profoundly affected the global distribution and diversity of plant pathogens. Identification of migration pathways can be used to monitor or manage pathogen movement for proactive disease management or quarantine measures. Genomics-based genetic marker discovery is allowing unprecedented collection of population genetic data for plant pathogens. These data can be used for detailed analysis of the ancestry of population samples and therefore for analysis of migration. Reconstruction of migration histories has confirmed previous hypotheses based on observational data and led to unexpected new findings on the origins of pathogens and source populations for past and recent migration. The choice of software for analysis depends on the type of migration being studied and the reproductive mode of the pathogen. Biased sampling and complex population structures are potential challenges to accurate inference of migration pathways.

The Top 10 oomycete pathogens in molecular plant pathology

Oomycetes form a deep lineage of eukaryotic organisms that includes a large number of plant pathogens which threaten natural and managed ecosystems. We undertook a survey to query the community for their ranking of plant-pathogenic oomycete species based on scientific and economic importance. In total, we received 263 votes from 62 scientists in 15 countries for a total of 33 species. The Top 10 species and their ranking are: (1) Phytophthora infestans; (2, tied) Hyaloperonospora arabidopsidis; (2, tied) Phytophthora ramorum; (4) Phytophthora sojae; (5) Phytophthora capsici; (6) Plasmopara viticola; (7) Phytophthora cinnamomi; (8, tied) Phytophthora parasitica; (8, tied) Pythium ultimum; and (10) Albugo candida. This article provides an introduction to these 10 taxa and a snapshot of current research. We hope that the list will serve as a benchmark for future trends in oomycete research.

Ancient genomics

The past decade has witnessed a revolution in ancient DNA (aDNA) research. Although the field's focus was previously limited to mitochondrial DNA and a few nuclear markers, whole genome sequences from the deep past can now be retrieved. This breakthrough is tightly connected to the massive sequence throughput of next generation sequencing platforms and the ability to target short and degraded DNA molecules. Many ancient specimens previously unsuitable for DNA analyses because of extensive degradation can now successfully be used as source materials. Additionally, the analytical power obtained by increasing the number of sequence reads to billions effectively means that contamination issues that have haunted aDNA research for decades, particularly in human studies, can now be efficiently and confidently quantified. At present, whole genomes have been sequenced from ancient anatomically modern humans, archaic hominins, ancient pathogens and megafaunal species. Those have revealed important functional and phenotypic information, as well as unexpected adaptation, migration and admixture patterns. As such, the field of aDNA has entered the new era of genomics and has provided valuable information when testing specific hypotheses related to the past.

Computational analyses of ancient pathogen DNA from herbarium samples: challenges and prospects

The application of DNA sequencing technology to the study of ancient DNA has enabled the reconstruction of past epidemics from genomes of historically important plant-associated microbes. Recently, the genome sequences of the potato late blight pathogen Phytophthora infestans were analyzed from 19th century herbarium specimens. These herbarium samples originated from infected potatoes collected during and after the Irish potato famine. Herbaria have therefore great potential to help elucidate past epidemics of crops, date the emergence of pathogens, and inform about past pathogen population dynamics. DNA preservation in herbarium samples was unexpectedly good, raising the possibility of a whole new research area in plant and microbial genomics. However, the recovered DNA can be extremely fragmented resulting in specific challenges in reconstructing genome sequences. Here we review some of the challenges in computational analyses of ancient DNA from herbarium samples. We also applied the recently developed linkage method to haplotype reconstruction of diploid or polyploid genomes from fragmented ancient DNA.

The changing face of the molecular evolutionary clock

The molecular clock has played an important role in biological research, both as a description of the evolutionary process and as a tool for inferring evolutionary timescales. Genomic data have provided valuable insights into the molecular clock, allowing the patterns and causes of evolutionary rate variation to be characterized in increasing detail. I explain how genome sequences offer exciting opportunities for estimating the timescale of the Tree of Life. I describe the different approaches that have been used to deal with the computational and statistical challenges encountered in molecular clock analyses of genomic data. Finally, I offer a perspective on the future of molecular clocks, highlighting some of the key limitations and the most promising research directions.