Mitotic Recombination and Rapid Genome Evolution in the Invasive Forest Pathogen Phytophthora ramorum

Recombination landscape and karyotypic variations revealed by linkage mapping in the grapevine downy mildew pathogen Plasmopara viticola

Plasmopara viticola, the causal agent of grapevine downy mildew, is a biotrophic oomycete engaged in a tight coevolutionary relationship with its host. Rapid adaptation of the pathogen is favored by annual sexual reproduction that generates genotypic diversity. With the aim of studying the recombination landscape across the P. viticola genome, we generated 2 half-sibling F1 progenies (N = 189 and 162). Using targeted SNP sequencing, between 1,405 and 1,894 markers were included in parental linkage maps, and a consensus map was obtained by integrating 4,509 markers. The reference genome could be assembled into 17 pseudochromosomes, anchoring 88% of its physical length. We observed a strong collinearity between parental genomes and extensive synteny with the downy mildew Peronospora effusa. In the consensus map, the median recombination rate was 13.8 cM/Mb. The local recombination rate was highly variable along chromosomes, and recombination was suppressed in putative centromeric regions. Recombination rate was found negatively correlated with repeats' coverage and positively correlated with gene coverage. However, genes encoding secreted proteins and putative effectors were underrepresented in highly recombining regions. In both progenies, about 5% of the individuals presented karyotypic anomalies. Aneuploidies and triploidies almost exclusively originated from the male-transmitted chromosomes. Triploids resulted from fertilization by diploid gametes, but also from dispermy. Obligatory sexual reproduction each year may explain the lower level of karyotypic variation in P. viticola compared to other oomycetes. The linkage maps will be useful to guide future de novo chromosome-scale assemblies of P. viticola genomes and to perform forward genetics.

Genomic and transcriptomic analyses of Phytophthora cinnamomi reveal complex genome architecture, expansion of pathogenicity factors, and host-dependent gene expression profiles

Phytophthora cinnamomi is a hemibiotrophic oomycete causing Phytophthora root rot in over 5,000 plant species, threatening natural ecosystems, forestry, and agriculture. Genomic studies of P. cinnamomi are limited compared to other Phytophthora spp. despite the importance of this destructive and highly invasive pathogen. The genome of two genetically and phenotypically distinct P. cinnamomi isolates collected from avocado orchards in California were sequenced using PacBio and Illumina sequencing. Genome sizes were estimated by flow cytometry and assembled de novo to 140-141 Mb genomes with 21,111-21,402 gene models. Genome analyses revealed that both isolates exhibited complex heterozygous genomes fitting the two-speed genome model. The more virulent isolate encodes a larger secretome and more RXLR effectors when compared to the less virulent isolate. Transcriptome analysis after P. cinnamomi infection in Arabidopsis thaliana, Nicotiana benthamiana, and Persea americana de Mill (avocado) showed that this pathogen deploys common gene repertoires in all hosts and host-specific subsets, especially among effectors. Overall, our results suggested that clonal P. cinnamomi isolates employ similar strategies as other Phytophthora spp. to increase phenotypic diversity (e.g., polyploidization, gene duplications, and a bipartite genome architecture) to cope with environmental changes. Our study also provides insights into common and host-specific P. cinnamomi infection strategies and may serve as a method for narrowing and selecting key candidate effectors for functional studies to determine their contributions to plant resistance or susceptibility.

Genome sequencing and comparative genome analysis of Rhizoctonia solani AG-3

Rhizoctonia solani AG-3 is a plant pathogenic fungus that belongs to the group of multinucleate Rhizoctonia. According to its internal transcribed spacer (ITS) cluster analysis and host range, it is divided into TB, PT, and TM subgroups. AG-3 TB mainly causes tobacco target spots, AG-3 PT mainly causes potato black scurf, and AG-3 TM mainly causes tomato leaf blight. In our previous study, we found that all 36 tobacco target spot strains isolated from Yunnan (Southwest China) were classified into AG-3 TB subgroup, while only two of the six tobacco target spot strains isolated from Liaoning (Northeast China) were classified into AG-3 TB subgroup, and the remaining four strains were classified into AG-3 TM subgroup, which had a unique taxonomic status, and there was no previous report on the whole genome information of AG-3 TM subgroup. In this study, the whole genomes of R. solani AG-3 strains 3T-1 (AG-3 TM isolated from Liaoning) and MJ-102 (AG-3 TB isolated from Yunnan) isolated from tobacco target spot in Liaoning and Yunnan were sequenced by IIumina and PacBio sequencing platforms. Comparative genomic analysis was performed with the previously reported AG-3 PT strain Rhs1AP, revealing their differences in genomes and virulence factors. The results indicated that the genome size of 3T-1 was 42,103,597 bp with 11,290 coding genes and 49.74% GC content, and the genome size of MJ-102 was 41,908,281 bp with 10,592 coding genes and 48.91% GC content. Through comparative genomic analysis with the previously reported strain Rhs1AP (AG-3 PT), it was found that the GC content between the genomes was similar, but the strains 3T-1 and MJ-102 contained more repetitive sequences. Similarly, there are similarities between their virulence factors, but there are also some differences. In addition, the results of collinearity analysis showed that 3T-1 and MJ-102 had lower similarity and longer evolutionary distance with Rhs1AP, but the genetic relationship between 3T-1 and MJ-102 was closer. This study can lay a foundation for studying the molecular pathogenesis and virulence factors of R. solani AG-3, and revealing its genomic composition will also help to develop more effective disease control strategies.

Genome-Enabled Insights into Downy Mildew Biology and Evolution

Oomycetes that cause downy mildew diseases are highly specialized, obligately biotrophic phytopathogens that can have major impacts on agriculture and natural ecosystems. Deciphering the genome sequence of these organisms provides foundational tools to study and deploy control strategies against downy mildew pathogens (DMPs). The recent telomere-to-telomere genome assembly of the DMP Peronospora effusa revealed high levels of synteny with distantly related DMPs, higher than expected repeat content, and previously undescribed architectures. This provides a road map for generating similar high-quality genome assemblies for other oomycetes. This review discusses biological insights made using this and other assemblies, including ancestral chromosome architecture, modes of sexual and asexual variation, the occurrence of heterokaryosis, candidate gene identification, functional validation, and population dynamics. We also discuss future avenues of research likely to be fruitful in studies of DMPs and highlight resources necessary for advancing our understanding and ability to forecast and control disease outbreaks.

Whole genome resequencing reveal patterns of genetic variation within Colletotrichum acutatum species complex from rubber trees in China

The Colletotrichum acutatum species complex possesses a diverse number of important traits, such as a wide host range and host preference, different modes of reproduction, and different strategies of host infection. Research using comparative genomics has attempted to find correlations between these traits. Here, we used multi-locus techniques and gene genealogical concordance analysis to investigate the phylogenetic relationships and taxonomic status of the Colletotrichum acutatum species complex using field isolates obtained from rubber trees. The results revealed that the dominant species was C. australisinense, followed by C. bannaense, while strain YNJH17109 was identified as C. laticiphilum. The taxonomic status of strains YNLC510 and YNLC511 was undetermined. Using whole-genome single nucleotide polymorphism data to analyze population structure, 18 strains of C. australisinense were subsequently divided into four populations, one of which was derived from an admixture of two populations. In addition, the strains LD1687, GD1628, and YNLC516, did not belong to any populations, and were considered to be admixtures of two or more populations. A split decomposition network analysis also provided evidence for genetic recombination within Colletotrichum acutatum species complex from rubber trees in China. Overall, a weak phylogeographic sub-structure was observed. Analysis also revealed significant differences in morphological characters and levels of virulence between populations.

Forest health in the Anthropocene: the emergence of a novel tree disease is associated with poplar cultivation

Plant domestication and movement are large contributors to the success of new diseases. The introduction of new host species can result in accelerated evolutionary changes in pathogens, affecting long-established coevolutionary dynamics. This has been observed in poplars where severe epidemics of pathogens that were innocuous in their natural pathosystems occurred following host domestication. The North American fungus Sphaerulina musiva is responsible for endemic leaf spots on Populus deltoides. We show that the expansion of poplar cultivation resulted in the emergence of a new lineage of this pathogen that causes stem infections on a new host, P. balsamifera. This suggests a host shift since this is not a known host. Genome analysis of this emerging lineage reveals a mosaic pattern with islands of diversity separated by fixed genome regions, which is consistent with a homoploid hybridization event between two individuals that produced a hybrid swarm. Genome regions of extreme divergence and low diversity are enriched in genes involved in host-pathogen interactions. The specialization of this emerging lineage to a new host and its clonal propagation represents a serious threat to poplars and could affect both natural and planted forests. This work provides a clear example of the changes created by the intensification of tree cultivation that facilitate the emergence of specialized pathogens, jeopardizing the natural equilibrium between hosts and pathogens. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

Chromosome-level assembly of the Phytophthora agathidicida genome reveals adaptation in effector gene families

Phytophthora species are notorious plant pathogens, with some causing devastating tree diseases that threaten the survival of their host species. One such example is Phytophthora agathidicida, the causal agent of kauri dieback - a root and trunk rot disease that kills the ancient, iconic and culturally significant tree species, Agathis australis (New Zealand kauri). A deeper understanding of how Phytophthora pathogens infect their hosts and cause disease is critical for the development of effective treatments. Such an understanding can be gained by interrogating pathogen genomes for effector genes, which are involved in virulence or pathogenicity. Although genome sequencing has become more affordable, the complete assembly of Phytophthora genomes has been problematic, particularly for those with a high abundance of repetitive sequences. Therefore, effector genes located in repetitive regions could be truncated or missed in a fragmented genome assembly. Using a combination of long-read PacBio sequences, chromatin conformation capture (Hi-C) and Illumina short reads, we assembled the P. agathidicida genome into ten complete chromosomes, with a genome size of 57 Mb including 34% repeats. This is the first Phytophthora genome assembled to chromosome level and it reveals a high level of syntenic conservation with the complete genome of Peronospora effusa, the only other completely assembled genome sequence of an oomycete. All P. agathidicida chromosomes have clearly defined centromeres and contain candidate effector genes such as RXLRs and CRNs, but in different proportions, reflecting the presence of gene family clusters. Candidate effector genes are predominantly found in gene-poor, repeat-rich regions of the genome, and in some cases showed a high degree of duplication. Analysis of candidate RXLR effector genes that occur in multicopy gene families indicated half of them were not expressed in planta. Candidate CRN effector gene families showed evidence of transposon-mediated recombination leading to new combinations of protein domains, both within and between chromosomes. Further analysis of this complete genome assembly will help inform new methods of disease control against P. agathidicida and other Phytophthora species, ultimately helping decipher how Phytophthora pathogens have evolved to shape their effector repertoires and how they might adapt in the future.

Genetic variation along an altitudinal gradient in the Phytophthora infestans effector gene Pi02860

Effector genes, together with climatic and other environmental factors, play multifaceted roles in the development of plant diseases. Understanding the role of environmental factors, particularly climate conditions affecting the evolution of effector genes, is important for predicting the long-term value of the genes in controlling agricultural diseases. Here, we collected Phytophthora infestans populations from five locations along a mountainous hill in China and sequenced the effector gene Pi02860 from >300 isolates. To minimize the influence of other ecological factors, isolates were sampled from the same potato cultivar on the same day. We also expressed the gene to visualise its cellular location, assayed its pathogenicity and evaluated its response to experimental temperatures. We found that Pi02860 exhibited moderate genetic variation at the nucleotide level which was mainly generated by point mutation. The mutations did not change the cellular location of the effector gene but significantly modified the fitness of P. infestans. Genetic variation and pathogenicity of the effector gene were positively associated with the altitude of sample sites, possibly due to increased mutation rate induced by the vertical distribution of environmental factors such as UV radiation and temperature. We further found that Pi02860 expression was regulated by experimental temperature with reduced expression as experimental temperature increased. Together, these results indicate that UV radiation and temperature are important environmental factors regulating the evolution of effector genes and provide us with considerable insight as to their future sustainable action under climate and other environmental change.

Population Genetic Analyses of Phytophthora cinnamomi Reveals Three Lineages and Movement Between Natural Vegetation and Avocado Orchards in South Africa

Phytophthora cinnamomi is the causal agent of root rot, canker, and dieback of thousands of plant species around the globe. This oomycete not only causes severe economic losses to forestry and agricultural industries, but also threatens the health of various plants in natural ecosystems. In this study, 380 isolates of P. cinnamomi from four avocado production areas and two regions of natural vegetation in South Africa were investigated using 15 microsatellite markers. These populations were found to have a low level of genetic diversity and consisted of isolates from three lineages. Shared genotypes were detected between isolates from avocado orchards and natural vegetation, indicating the movement of isolates between these areas. The population from the Western Cape natural vegetation had the highest level of genotypic diversity and number of unique alleles, indicating this could be the point of introduction of P. cinnamomi to South Africa. Index of association analysis suggested that five of six populations were under linkage disequilibrium, suggesting a clonal mode of reproduction, whereas genotypes sampled from a recently established avocado orchard in the Western Cape were derived from a randomly recombining population. This study provided novel insights on the genetic diversity and spread of P. cinnamomi in South Africa. It also reported on the predominance of triploidy in natural occurring populations and provided evidence for recombination of P. cinnamomi for the first time. The presence of two dominant genotypes in all avocado production areas in South Africa highlight the importance of considering them in disease management and resistance breeding programs.

Phytophthora: an ancient, historic, biologically and structurally cohesive and evolutionarily successful generic concept in need of preservation

The considerable economic and social impact of the oomycete genus Phytophthora is well known. In response to evidence that all downy mildews (DMs) reside phylogenetically within Phytophthora, rendering Phytophthora paraphyletic, a proposal has been made to split the genus into multiple new genera. We have reviewed the status of the genus and its relationship to the DMs. Despite a substantial increase in the number of described species and improvements in molecular phylogeny the Phytophthora clade structure has remained stable since first demonstrated in 2000. Currently some 200 species are distributed across twelve major clades in a relatively tight monophyletic cluster. In our assessment of 196 species for twenty morphological and behavioural criteria the clades show good biological cohesion. Saprotrophy, necrotrophy and hemi-biotrophy of woody and non-woody roots, stems and foliage occurs across the clades. Phylogenetically less related clades often show strong phenotypic and behavioural similarities and no one clade or group of clades shows the synapomorphies that might justify a unique generic status. We propose the clades arose from the migration and worldwide radiation ~ 140 Mya (million years ago) of an ancestral Gondwanan Phytophthora population, resulting in geographic isolation and clade divergence through drift on the diverging continents combined with adaptation to local hosts, climatic zones and habitats. The extraordinary flexibility of the genus may account for its global 'success'. The 20 genera of the obligately biotrophic, angiosperm-foliage specialised DMs evolved from Phytophthora at least twice via convergent evolution, making the DMs as a group polyphyletic and Phytophthora paraphyletic in cladistic terms. The long phylogenetic branches of the DMs indicate this occurred rather rapidly, via paraphyletic evolutionary 'jumps'. Such paraphyly is common in successful organisms. The proposal to divide Phytophthora appears more a device to address the issue of the convergent evolution of the DMs than the structure of Phytophthora per se. We consider it non-Darwinian, putting the emphasis on the emergent groups (the DMs) rather than the progenitor (Phytophthora) and ignoring the evolutionary processes that gave rise to the divergence. Further, the generic concept currently applied to the DMs is narrower than that between some closely related Phytophthora species. Considering the biological and structural cohesion of Phytophthora, its historic and social impacts and its importance in scientific communication and biosecurity protocol, we recommend that the current broad generic concept is retained by the scientific community.

Genomic biosurveillance detects a sexual hybrid in the sudden oak death pathogen

Invasive exotic pathogens pose a threat to trees and forest ecosystems worldwide, hampering the provision of essential ecosystem services such as carbon sequestration and water purification. Hybridization is a major evolutionary force that can drive the emergence of pathogens. Phytophthora ramorum, an emergent pathogen that causes the sudden oak and larch death, spreads as reproductively isolated divergent clonal lineages. We use a genomic biosurveillance approach by sequencing genomes of P. ramorum from survey and inspection samples and report the discovery of variants of P. ramorum that are the result of hybridization via sexual recombination between North American and European lineages. We show that these hybrids are viable, can infect a host and produce spores for long-term survival and propagation. Genome sequencing revealed genotypic combinations at 54,515 single nucleotide polymorphism loci not present in parental lineages. More than 6,000 of those genotypes are predicted to have a functional impact in genes associated with host infection, including effectors, carbohydrate-active enzymes and proteases. We also observed post-meiotic mitotic recombination that could generate additional genotypic and phenotypic variation and contribute to homoploid hybrid speciation. Our study highlights the importance of plant pathogen biosurveillance to detect variants, including hybrids, and inform management and control.

Sexual reproduction contributes to the evolution of resistance-breaking isolates of the spinach pathogen Peronospora effusa

Peronospora effusa causes downy mildew, the economically most important disease of cultivated spinach worldwide. To date, 19 P. effusa races have been denominated based on their capacity to break spinach resistances, but their genetic diversity and the evolutionary processes that contribute to race emergence are unknown. Here, we performed the first systematic analysis of P. effusa races showing that those emerge by both asexual and sexual reproduction. Specifically, we studied the diversity of 26 P. effusa isolates from 16 denominated races based on mitochondrial and nuclear comparative genomics. Mitochondrial genomes based on long-read sequencing coupled with diversity assessment based on short-read sequencing uncovered two mitochondrial haplogroups, each with distinct genome organization. Nuclear genome-wide comparisons of the 26 isolates revealed that 10 isolates from six races could clearly be divided into three asexually evolving groups, in concordance with their mitochondrial phylogeny. The remaining isolates showed signals of reticulated evolution and discordance between nuclear and mitochondrial phylogenies, suggesting that these evolved through sexual reproduction. Increased understanding of this pathogen's reproductive modes will provide the framework for future studies into the molecular mechanisms underlying race emergence and into the P. effusa-spinach interaction, thus assisting in sustainable production of spinach through knowledge-driven resistance breeding.

Comparative Genomic Analysis Reveals Genetic Variation and Adaptive Evolution in the Pathogenicity-Related Genes of Phytophthora capsici

Phytophthora capsici is an oomycete pathogen responsible for damping off, root rot, fruit rot, and foliar blight in popular vegetable and legume crops. The existence of distinct aggressiveness levels and physiological races among the P. capsici population is a major constraint to developing resistant varieties of host crops. In the present study, we compared the genomes of three P. capsici isolates with different aggressiveness levels to reveal their genomic differences. We obtained genome sequences using short-read and long-read technologies, which yielded an average genome size of 76 Mbp comprising 514 contigs and 15,076 predicted genes. A comparative genomic analysis uncovered the signatures of accelerated evolution, gene family expansions in the pathogenicity-related genes among the three isolates. Resequencing two additional P. capsici isolates enabled the identification of average 1,023,437 SNPs, revealing the frequent accumulation of non-synonymous substitutions in pathogenicity-related gene families. Furthermore, pathogenicity-related gene families, cytoplasmic effectors and ATP binding cassette (ABC) transporters, showed expansion signals in the more aggressive isolates, with a greater number of non-synonymous SNPs. This genomic information explains the plasticity, difference in aggressiveness levels, and genome structural variation among the P. capsici isolates, providing insight into the genomic features related to the evolution and pathogenicity of this oomycete pathogen.

Genomic Approaches to Plant-Pathogen Epidemiology and Diagnostics

Diseases have a significant cost to agriculture. Findings from analyses of whole-genome sequences show great promise for informing strategies to mitigate risks from diseases caused by phytopathogens. Genomic approaches can be used to dramatically shorten response times to outbreaks and inform disease management in novel ways. However, the use of these approaches requires expertise in working with big, complex data sets and an understanding of their pitfalls and limitations to infer well-supported conclusions. We suggest using an evolutionary framework to guide the use of genomic approaches in epidemiology and diagnostics of plant pathogens. We also describe steps that are necessary for realizing these as standard approaches in disease surveillance.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

How Do Pathogens Evolve Novel Virulence Activities?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.We consider the state of knowledge on pathogen evolution of novel virulence activities, broadly defined as anything that increases pathogen fitness with the consequence of causing disease in either the qualitative or quantitative senses, including adaptation of pathogens to host immunity and physiology, host species, genotypes, or tissues, or the environment. The evolution of novel virulence activities as an adaptive trait is based on the selection exerted by hosts on variants that have been generated de novo or arrived from elsewhere. In addition, the biotic and abiotic environment a pathogen experiences beyond the host may influence pathogen virulence activities. We consider host-pathogen evolution, host range expansion, and external factors that can mediate pathogen evolution. We then discuss the mechanisms by which pathogens generate and recombine the genetic variation that leads to novel virulence activities, including DNA point mutation, transposable element activity, gene duplication and neofunctionalization, and genetic exchange. In summary, if there is an (epi)genetic mechanism that can create variation in the genome, it will be used by pathogens to evolve virulence factors. Our knowledge of virulence evolution has been biased by pathogen evolution in response to major gene resistance, leaving other virulence activities underexplored. Understanding the key driving forces that give rise to novel virulence activities and the integration of evolutionary concepts and methods with mechanistic research on plant-microbe interactions can help inform crop protection.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome of the destructive oomycete Phytophthora cinnamomi provides insights into its pathogenicity and adaptive potential

BACKGROUND

Phytophthora cinnamomi is an oomycete pathogen of global relevance. It is considered as one of the most invasive species, which has caused irreversible damage to natural ecosystems and horticultural crops. There is currently a lack of a high-quality reference genome for this species despite several attempts that have been made towards sequencing its genome. The lack of a good quality genome sequence has been a setback for various genetic and genomic research to be done on this species. As a consequence, little is known regarding its genome characteristics and how these contribute to its pathogenicity and invasiveness.

RESULTS

In this work we generated a high-quality genome sequence and annotation for P. cinnamomi using a combination of Oxford Nanopore and Illumina sequencing technologies. The annotation was done using RNA-Seq data as supporting gene evidence. The final assembly consisted of 133 scaffolds, with an estimated genome size of 109.7 Mb, N50 of 1.18 Mb, and BUSCO completeness score of 97.5%. Genome partitioning analysis revealed that P. cinnamomi has a two-speed genome characteristic, similar to that of other oomycetes and fungal plant pathogens. In planta gene expression analysis revealed up-regulation of pathogenicity-related genes, suggesting their important roles during infection and host degradation.

CONCLUSION

This study has provided a high-quality reference genome and annotation for P. cinnamomi. This is among the best assembled genomes for any Phytophthora species assembled to date and thus resulted in improved identification and characterization of pathogenicity-related genes, some of which were undetected in previous versions of genome assemblies. Phytophthora cinnamomi harbours a large number of effector genes which are located in the gene-poor regions of the genome. This unique genomic partitioning provides P. cinnamomi with a high level of adaptability and could contribute to its success as a highly invasive species. Finally, the genome sequence, its annotation and the pathogenicity effectors identified in this study will serve as an important resource that will enable future studies to better understand and mitigate the impact of this important pathogen.

The Destructive Tree Pathogen Phytophthora ramorum Originates from the Laurosilva Forests of East Asia

As global plant trade expands, tree disease epidemics caused by pathogen introductions are increasing. Since ca 2000, the introduced oomycete Phytophthora ramorum has caused devastating epidemics in Europe and North America, spreading as four ancient clonal lineages, each of a single mating type, suggesting different geographical origins. We surveyed laurosilva forests for P. ramorum around Fansipan mountain on the Vietnam-China border and on Shikoku and Kyushu islands, southwest Japan. The surveys yielded 71 P. ramorum isolates which we assigned to eight new lineages, IC1 to IC5 from Vietnam and NP1 to NP3 from Japan, based on differences in colony characteristics, gene x environment responses and multigene phylogeny. Molecular phylogenetic trees and networks revealed the eight Asian lineages were dispersed across the topology of the introduced European and North American lineages. The deepest node within P. ramorum, the divergence of lineages NP1 and NP2, was estimated at 0.5 to 1.6 Myr. The Asian lineages were each of a single mating type, and at some locations, lineages of "opposite" mating type were present, suggesting opportunities for inter-lineage recombination. Based on the high level of phenotypic and phylogenetic diversity in the sample populations, the coalescence results and the absence of overt host symptoms, we conclude that P. ramorum comprises many anciently divergent lineages native to the laurosilva forests between eastern Indochina and Japan.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

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

Phylogeographic Approaches to Characterize the Emergence of Plant Pathogens

Phylogeography combines geographic information with phylogenetic and population genomic approaches to infer the evolutionary history of a species or population in a geographic context. This approach has been instrumental in understanding the emergence, spread, and evolution of a range of plant pathogens. In particular, phylogeography can address questions about where a pathogen originated, whether it is native or introduced, and when and how often introductions occurred. We review the theory, methods, and approaches underpinning phylogeographic inference and highlight applications providing novel insights into the emergence and spread of select pathogens. We hope that this review will be useful in assessing the power, pitfalls, and opportunities presented by various phylogeographic approaches.

Population Genomics of Filamentous Plant Pathogens-A Brief Overview of Research Questions, Approaches, and Pitfalls

With ever-decreasing sequencing costs, research on the population biology of plant pathogens is transitioning from population genetics-using dozens of genetic markers or polymorphism data of several genes-to population genomics-using several hundred to tens of thousands of markers or whole-genome sequence data. The field of population genomics is characterized by rapid theoretical and methodological advances and by numerous steps and pitfalls in its technical and analytical workflow. In this article, we aim to provide a brief overview of topics relevant to the study of population genomics of filamentous plant pathogens and direct readers to more extensive reviews for in-depth understanding. We briefly discuss different types of population genomics-inspired research questions and give insights into the sampling strategies that can be used to answer such questions. We then consider different sequencing strategies, the various options available for data processing, and some of the currently available tools for population genomic data analysis. We conclude by highlighting some of the hurdles along the population genomic workflow, providing cautionary warnings relative to assumptions and technical challenges, and presenting our own future perspectives of the field of population genomics for filamentous plant pathogens.

Patterns of Genetic Diversification in the Invasive Hybrid Plant Pathogen Phytophthora × alni and Its Parental Species P. uniformis

In pathogenic fungi and oomycetes, interspecific hybridization may lead to the formation of new species having a greater impact on natural ecosystems than the parental species. From the early 1990s, a severe alder (Alnus spp.) decline due to an unknown Phytophthora species was observed in several European countries. Genetic analyses revealed that the disease was caused by the triploid hybrid P. × alni, which originated in Europe from the hybridization of P. uniformis and P. × multiformis. Here, we investigated the population structure of P. × alni (158 isolates) and P. uniformis (85 isolates) in several European countries using microsatellite markers. Our analyses confirmed the genetic structure previously observed in other European populations, with P. uniformis populations consisting of at most two multilocus genotypes (MLGs) and P. × alni populations dominated by MLG Pxa-1. The genetic structure of P. × alni populations in the Czech Republic, Hungary and Sweden seemed to reflect the physical isolation of river systems. Most rare P. × alni MLGs showed a loss of heterozygosity (LOH) at one or a few microsatellite loci compared with other MLGs. This LOH may allow a stabilization within the P. × alni genome or a rapid adaptation to stress situations. Alternatively, alleles may be lost because of random genetic drift in small, isolated populations, with no effect on fitness of P. × alni. Additional studies would be necessary to confirm these patterns of population diversification and to better understand the factors driving it.

Genome-wide microsatellites and species specific markers in genus Phytophthora revealed through whole genome analysis

Genome wide microsatellite maps shall support Phytophthora systematics through the development of reliable markers, enabling species discrimination and variability analyses. Whole genome sequences of 17 Phytophthora accessions belonging to 14 species were retrieved from GenBank and the genome-wide microsatellites in each species were mined. A total of 51,200 microsatellites, including dinucleotide to decanucleotide motifs, have been identified across all the species and each one was characterized for uniqueness and repeat number. The P. infestans T30-4 genome had the highest (6873) and P. multivora 3378 had the lowest number of microsatellites (1802). Dinucleotide motifs (63.6%) followed by trinucleotide motifs (30.1%) were most abundant in all the genome. From 14 species, 250 microsatellites which are unique for the respective genomes are detailed along with their primer combinations and product sizes. P. sojae had the highest number of unique microsatellite motifs. Genome wide microsatellite maps for all the 14 Phytophthora species including the chromosome, position, motif, repeat number, forward and reverse primer sequences and expected PCR product size, for every microsatellite are presented. Markers based on the unique microsatellites could be used to identify each species, whereas the ones common to all species could be used to identify the genetic variability. Furthermore, to confirm the results, pure cultures of P. capsici, P. nicotianae and P. palmivora were procured from the Phytophthora Repository, DNA was isolated and the unique markers were screened across the species. The characteristic markers developed have confirmed the genome analysis results.

Loss of Heterozygosity and Base Mutation Rates Vary Among Saccharomyces cerevisiae Hybrid Strains

A growing body of evidence suggests that mutation rates exhibit intra-species specific variation. We estimated genome-wide loss of heterozygosity (LOH), gross chromosomal changes, and single nucleotide mutation rates to determine intra-species specific differences in hybrid and homozygous strains of Saccharomyces cerevisiae. The mutation accumulation lines of the S. cerevisiae hybrid backgrounds - S288c/YJM789 (S/Y) and S288c/RM11-1a (S/R) were analyzed along with the homozygous diploids RM11, S288c, and YJM145. LOH was extensive in both S/Y and S/R hybrid backgrounds. The S/Y background also showed longer LOH tracts, gross chromosomal changes, and aneuploidy. Short copy number aberrations were observed in the S/R background. LOH data from the S/Y and S/R hybrids were used to construct a LOH map for S288c to identify hotspots. Further, we observe up to a sixfold difference in single nucleotide mutation rates among the S. cerevisiae S/Y and S/R genetic backgrounds. Our results demonstrate LOH is common during mitotic divisions in S. cerevisiae hybrids and also highlight genome-wide differences in LOH patterns and rates of single nucleotide mutations between commonly used S. cerevisiae hybrid genetic backgrounds.

Assessing Niche Shifts and Conservatism by Comparing the Native and Post-Invasion Niches of Major Forest Invasive Species

Invasive species experience biotic and abiotic conditions that may (or may not) resemble their native environment. We explored the methodology of determining climatic niches and compared the native and post-invasion niches of four invasive forest pests to determine if these species experienced shifts or changes in their new climatic niches. We used environmental principle components analysis (PCA-env) method to quantify climatic niche shifts, expansions, and temporal changes. Furthermore, we assessed the effect of variable selection in the delineation and comparison of niche space. We found that variable selection influenced the delineation and overlap of each niche, whereas the subset of climatic variables selected from the first two PCA-env axes explained more variance in environmental conditions than the complete set of climatic variables for all four species. Most focal species showed climatic niche shifts in their invasive range and had not yet fully occupied the available niche within the invaded range. Our species varied the proportion of niche overlap between the native and invasive ranges. By comparing native and invasive niches, we can help predict a species’ potential range expansion and invasion potential. Our results can guide monitoring and help inform management of these and other invasive species.

Independent Whole-Genome Duplications Define the Architecture of the Genomes of the Devastating West African Cacao Black Pod Pathogen Phytophthora megakarya and Its Close Relative Phytophthora palmivora

Phytophthora megakarya and P. palmivora are oomycete pathogens that cause black pod rot of cacao (Theobroma cacao), the most economically important disease on cacao globally. While P. palmivora is a cosmopolitan pathogen, P. megakarya, which is more aggressive on cacao than P. palmivora, has been reported only in West and Central Africa where it has been spreading and devastating cacao farms since the 1950s. In this study, we reconstructed the complete diploid genomes of multiple isolates of both species using single-molecule real-time sequencing. Thirty-one additional genotypes were sequenced to analyze inter- and intra-species genomic diversity. The P. megakarya genome is exceptionally large (222 Mbp) and nearly twice the size of P. palmivora (135 Mbp) and most known Phytophthora species (∼100 Mbp on average). Previous reports pointed toward a whole-genome duplication (WGD) in P. palmivora In this study, we demonstrate that both species underwent independent and relatively recent WGD events. In P. megakarya we identified a unique combination of WGD and large-scale transposable element driven genome expansion, which places this genome in the upper range of Phytophthora genome sizes, as well as effector pools with 1,382 predicted RxLR effectors. Finally, this study provides evidence of adaptive evolution of effectors like RxLRs and Crinklers, and discusses the implications of effector expansion and diversification.

Asexual Evolution and Forest Conditions Drive Genetic Parallelism in Phytophthora ramorum

It is commonly assumed that asexual lineages are short-lived evolutionarily, yet many asexual organisms can generate genetic and phenotypic variation, providing an avenue for further evolution. Previous work on the asexual plant pathogen Phytophthora ramorum NA1 revealed considerable genetic variation in the form of Structural Variants (SVs). To better understand how SVs arise and their significance to the California NA1 population, we studied the evolutionary histories of SVs and the forest conditions associated with their emergence. Ancestral state reconstruction suggests that SVs arose by somatic mutations among multiple independent lineages, rather than by recombination. We asked if this unusual phenomenon of parallel evolution between isolated populations is transmitted to extant lineages and found that SVs persist longer in a population if their genetic background had a lower mutation load. Genetic parallelism was also found in geographically distant demes where forest conditions such as host density, solar radiation, and temperature, were similar. Parallel SVs overlap with genes involved in pathogenicity such as RXLRs and have the potential to change the course of an epidemic. By combining genomics and environmental data, we identified an unexpected pattern of repeated evolution in an asexual population and identified environmental factors potentially driving this phenomenon.

Recent advances in oomycete genomics

The oomycetes are a class of ubiquitous, filamentous microorganisms that include some of the biggest threats to global food security and natural ecosystems. Within the oomycete class are highly diverse species that infect a broad range of animals and plants. Some of the most destructive plant pathogens are oomycetes, such as Phytophthora infestans, the agent of potato late blight and the cause of the Irish famine. Recent years have seen a dramatic increase in the number of sequenced oomycete genomes. Here we review the latest developments in oomycete genomics and some of the important insights that have been gained. Coupled with proteomic and transcriptomic analyses, oomycete genome sequences have revealed tremendous insights into oomycete biology, evolution, genome organization, mechanisms of infection, and metabolism. We also present an updated phylogeny of the oomycete class using a phylogenomic approach based on the 65 oomycete genomes that are currently available.

Editing of an effector gene promoter sequence impacts plant-Phytophthora interaction

Pathogen avirulence (Avr) effectors interplay with corresponding plant resistance (R) proteins and activate robust plant immune responses. Although the expression pattern of Avr genes has been tied to their functions for a long time, it is still not clear how Avr gene expression patterns impact plant-microbe interactions. Here, we selected PsAvr3b, which shows a typical effector gene expression pattern from a soybean root pathogen Phytophthora sojae. To modulate gene expression, we engineered PsAvr3b promoter sequences by in situ substitution with promoter sequences from Actin (constitutive expression), PsXEG1 (early expression), and PsNLP1 (later expression) using the CRISPR/Cas9. PsAvr3b driven by different promoters resulted in distinct expression levels across all the tested infection time points. Importantly, those mutants with low PsAvr3b expression successfully colonized soybean plants carrying the cognate R gene Rps3b. To dissect the difference in plant responses to the PsAvr3b expression level, we conducted RNA-sequencing of different infection samples at 24 h postinfection and found soybean immune genes, including a few previously unknown genes that are associated with resistance. Our study highlights that fine-tuning in Avr gene expression impacts the compatibility of plant disease and provides clues to improve crop resistance in disease control management.

Molecular assays to detect the presence and viability of Phytophthora ramorum and Grosmannia clavigera

Wood and wood products can harbor microorganisms that can raise phytosanitary concerns in countries importing or exporting these products. To evaluate the efficacy of wood treatment on the survival of microorganisms of phytosanitary concern the method of choice is to grow microbes in petri dishes for subsequent identification. However, some plant pathogens are difficult or impossible to grow in axenic cultures. A molecular methodology capable of detecting living fungi and fungus-like organisms in situ can provide a solution. RNA represents the transcription of genes and can become rapidly unstable after cell death, providing a proxy measure of viability. We designed and used RNA-based molecular diagnostic assays targeting genes essential to vital processes and assessed their presence in wood colonized by fungi and oomycetes through reverse transcription and real-time polymerase chain reaction (PCR). A stability analysis was conducted by comparing the ratio of mRNA to gDNA over time following heat treatment of mycelial cultures of the Oomycete Phytophthora ramorum and the fungus Grosmannia clavigera. The real-time PCR results indicated that the DNA remained stable over a period of 10 days post treatment in heat-treated samples, whereas mRNA could not be detected after 24 hours for P. ramorum or 96 hours for G. clavigera. Therefore, this method provides a reliable way to evaluate the viability of these pathogens and offers a potential way to assess the effectiveness of existing and emerging wood treatments. This can have important phytosanitary impacts on assessing both timber and non-timber forest products of commercial value in international wood trade.

Genomic biosurveillance of forest invasive alien enemies: A story written in code

The world's forests face unprecedented threats from invasive insects and pathogens that can cause large irreversible damage to the ecosystems. This threatens the world's capacity to provide long-term fiber supply and ecosystem services that range from carbon storage, nutrient cycling, and water and air purification, to soil preservation and maintenance of wildlife habitat. Reducing the threat of forest invasive alien species requires vigilant biosurveillance, the process of gathering, integrating, interpreting, and communicating essential information about pest and pathogen threats to achieve early detection and warning and to enable better decision-making. This process is challenging due to the diversity of invasive pests and pathogens that need to be identified, the diverse pathways of introduction, and the difficulty in assessing the risk of establishment. Genomics can provide powerful new solutions to biosurveillance. The process of invasion is a story written in four chapters: transport, introduction, establishment, and spread. The series of processes that lead to a successful invasion can leave behind a DNA signature that tells the story of an invasion. This signature can help us understand the dynamic, multistep process of invasion and inform management of current and future introductions. This review describes current and future application of genomic tools and pipelines that will provide accurate identification of pests and pathogens, assign outbreak or survey samples to putative sources to identify pathways of spread, and assess risk based on traits that impact the outbreak outcome.

Variation in Susceptibility of Tanoak to the NA1 and EU1 Lineages of Phytophthora ramorum, the Cause of Sudden Oak Death

Phytophthora ramorum, the cause of sudden oak death (SOD), kills tanoak (Notholithocarpus densiflorus) trees in southwestern Oregon and California. Two lineages of P. ramorum are now found in wildland forests of Oregon (NA1 and EU1). In addition to the management of SOD in forest ecosystems, disease resistance could be used as a way to mitigate the impact of P. ramorum. The objectives of this study were to (i) characterize the variability in resistance of N. densiflorus among families using lesion length; (ii) determine whether lineage, isolate, family, or their interactions significantly affect variation in lesion length; and (iii) determine whether there are differences among isolates and among families in terms of lesion length. The parameters isolate nested within lineage (isolate[lineage]) and family × isolate(lineage) interaction explained the majority of the variation in lesion length. There was no significant difference between the NA1 and EU1 lineages in terms of mean lesion length; however, there were differences among the six isolates. Lesions on seedlings collected from surviving trees at infested sites were smaller, on average, than lesions of seedlings collected from trees at noninfested sites (P = 0.0064). The results indicate that there is potential to establish a breeding program for tanoak resistance to SOD and that several isolates of P. ramorum should be used in an artificial inoculation assay.

Expansion of LINEs and species-specific DNA repeats drives genome expansion in Asian Gypsy Moths

Two subspecies of Asian gypsy moth (AGM), Lymantria dispar asiatica and L. dispar japonica, pose a serious alien invasive threat to North American forests. Despite decades of research on the ecology and biology of this pest, limited AGM-specific genomic resources are currently available. Here, we report on the genome sequences and functional content of these AGM subspecies. The genomes of L.d. asiatica and L.d. japonica are the largest lepidopteran genomes sequenced to date, totaling 921 and 999 megabases, respectively. Large genome size in these subspecies is driven by the accumulation of specific classes of repeats. Genome-wide metabolic pathway reconstructions suggest strong genomic signatures of energy-related pathways in both subspecies, dominated by metabolic functions related to thermogenesis. The genome sequences reported here will provide tools for probing the molecular mechanisms underlying phenotypic traits that are thought to enhance AGM invasiveness.

A streamlined and predominantly diploid genome in the tiny marine green alga Chloropicon primus

Tiny marine green algae issued from two deep branches of the Chlorophyta, the Mamiellophyceae and Chloropicophyceae, dominate different regions of the oceans and play key roles in planktonic communities. Considering that the Mamiellophyceae is the sole lineage of prasinophyte algae that has been intensively investigated, the extent to which these two algal groups differ in their metabolic capacities and cellular processes is currently unknown. To address this gap of knowledge, we investigate here the nuclear genome sequence of a member of the Chloropicophyceae, Chloropicon primus. Among the main biological insights that emerge from this 17.4 Mb genome, we find an unexpected diploid structure for most chromosomes and a propionate detoxification pathway in green algae. Our results support the notion that separate events of genome minimization, which entailed differential losses of genes/pathways, have occurred in the Chloropicophyceae and Mamiellophyceae, suggesting different strategies of adaptation to oceanic environments.

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.