Draft genome sequences of Phytophthora kernoviae and Phytophthora ramorum lineage EU2 from Scotland

Krisp: A Python package to aid in the design of CRISPR and amplification-based diagnostic assays from whole genome sequencing data

Recent pandemics like COVID-19 highlighted the importance of rapidly developing diagnostics to detect evolving pathogens. CRISPR-Cas technology has recently been used to develop diagnostic assays for sequence-specific recognition of DNA or RNA. These assays have similar sensitivity to the gold standard qPCR but can be deployed as easy to use and inexpensive test strips. However, the discovery of diagnostic regions of a genome flanked by conserved regions where primers can be designed requires extensive bioinformatic analyses of genome sequences. We developed the Python package krisp to aid in the discovery of primers and diagnostic sequences that differentiate groups of samples from each other, using either unaligned genome sequences or a variant call format (VCF) file as input. Krisp has been optimized to handle large datasets by using efficient algorithms that run in near linear time, use minimal RAM, and leverage parallel processing when available. The validity of krisp results has been demonstrated in the laboratory with the successful design of a CRISPR diagnostic assay to distinguish the sudden oak death pathogen Phytophthora ramorum from closely related Phytophthora species. Krisp is released open source under a permissive license with all the documentation needed to quickly design CRISPR-Cas diagnostic assays.

Draft genome assemblies for tree pathogens Phytophthora pseudosyringae and Phytophthora boehmeriae

Species of Phytophthora, plant pathogenic eukaryotic microbes, can cause disease on many tree species. Genome sequencing of species from this genus has helped to determine components of their pathogenicity arsenal. Here, we sequenced genomes for two widely distributed species, Phytophthora pseudosyringae and Phytophthora boehmeriae, yielding genome assemblies of 49 and 40 Mb, respectively. We identified more than 270 candidate disease promoting RXLR effector coding genes for each species, and hundreds of genes encoding candidate plant cell wall degrading carbohydrate active enzymes (CAZymes). These data boost genome sequence representation across the Phytophthora genus, and form resources for further study of Phytophthora pathogenesis.

Haustorium formation and a distinct biotrophic transcriptome characterize infection of Nicotiana benthamiana by the tree pathogen Phytophthora kernoviae

Phytophthora species cause some of the most serious diseases of trees and threaten forests in many parts of the world. Despite the generation of genome sequence assemblies for over 10 tree-pathogenic Phytophthora species and improved detection methods, there are many gaps in our knowledge of how these pathogens interact with their hosts. To facilitate cell biology studies of the infection cycle we examined whether the tree pathogen Phytophthora kernoviae could infect the model plant Nicotiana benthamiana. We transformed P. kernoviae to express green fluorescent protein (GFP) and demonstrated that it forms haustoria within infected N. benthamiana cells. Haustoria were also formed in infected cells of natural hosts, Rhododendron ponticum and European beech (Fagus sylvatica). We analysed the transcriptome of P. kernoviae in cultured mycelia, spores, and during infection of N. benthamiana, and detected 12,559 transcripts. Of these, 1,052 were predicted to encode secreted proteins, some of which may function as effectors to facilitate disease development. From these, we identified 87 expressed candidate RXLR (Arg-any amino acid-Leu-Arg) effectors. We transiently expressed 12 of these as GFP fusions in N. benthamiana leaves and demonstrated that nine significantly enhanced P. kernoviae disease progression and diversely localized to the cytoplasm, nucleus, nucleolus, and plasma membrane. Our results show that N. benthamiana can be used as a model host plant for studying this tree pathogen, and that the interaction likely involves suppression of host immune responses by RXLR effectors. These results establish a platform to expand the understanding of Phytophthora tree diseases.

Comparative analyses of saprotrophy in Salisapilia sapeloensis and diverse plant pathogenic oomycetes reveal lifestyle-specific gene expression

Oomycetes include many devastating plant pathogens. Across oomycete diversity, plant-infecting lineages are interspersed by non-pathogenic ones. Unfortunately, our understanding of the evolution of lifestyle switches is hampered by a scarcity of data on the molecular biology of saprotrophic oomycetes, ecologically important primary colonizers of dead tissue that can serve as informative reference points for understanding the evolution of pathogens. Here, we established Salisapilia sapeloensis as a tractable system for the study of saprotrophic oomycetes. We generated multiple transcriptomes from S. sapeloensis and compared them with (i) 22 oomycete genomes and (ii) the transcriptomes of eight pathogenic oomycetes grown under 13 conditions. We obtained a global perspective on gene expression signatures of oomycete lifestyles. Our data reveal that oomycete saprotrophs and pathogens use similar molecular mechanisms for colonization but exhibit distinct expression patterns. We identify a S. sapeloensis-specific array and expression of carbohydrate-active enzymes and putative regulatory differences, highlighted by distinct expression levels of transcription factors. Salisapilia sapeloensis expresses only a small repertoire of candidates for virulence-associated genes. Our analyses suggest lifestyle-specific gene regulatory signatures and that, in addition to variation in gene content, shifts in gene regulatory networks underpin the evolution of oomycete lifestyles.

A Global Survey of Carbohydrate Esterase Families 1 and 10 in Oomycetes

Carbohydrate-active enzymes (CAZymes) are a cornerstone in the phytopathogenicity of filamentous microbes. CAZymes are required for every step of a successful infection cycle-from penetration, to nutrient acquisition (during colonization), to exit and dispersal. Yet, CAZymes are not a unique feature of filamentous pathogens. They are found across eukaryotic genomes and including, for example, saprotrophic relatives of major pathogens. Comparative genomics and functional analyses revealed that CAZyme content is shaped by a multitude of factors, including utilized substrate, lifestyle, and host preference. Yet, family size alone says little about usage. Indeed, in a previous study, we found that genes putatively coding for the CAZyme families of carbohydrate esterase (CE)1 and CE10, while not specifically enriched in number, were suggested to have lifestyle-specific gene expression patterns. Here, we used comparative genomics and a clustering approach to understand how the repertoire of the CE1- and CE10-encoding gene families is shaped across oomycete evolution. These data are combined with comparative transcriptomic analyses across homologous clusters within the gene families. We find that CE1 and CE10 have been reduced in number in biotrophic oomycetes independent of the phylogenetic relationship of the biotrophs to each other. The reduction in CE1 is different from that observed for CE10: While in CE10 specific clusters of homologous sequences show convergent reduction, CE1 reduction is caused by species-specific losses. Comparative transcriptomics revealed that some clusters of CE1 or CE10 sequences have a higher expression than others, independent of the species composition within them. Further, we find that CE1- and CE10-encoding genes are mainly induced in plant pathogens and that some homologous genes show lifestyle-specific gene expression levels during infection, with hemibiotrophs showing the highest expression levels.

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.

Genome sequencing of oomycete isolates from Chile supports the New Zealand origin of Phytophthora kernoviae and makes available the first Nothophytophthora sp. genome

Genome sequences were generated for six oomycete isolates collected from forests in Valdivia, Chile. Three of the isolates were identified morphologically as Phytophthora kernoviae, whereas two were similar to other clade 10 Phytophthora species. One isolate was tentatively identified as Nothophytophthora valdiviana based on nucleotide sequence similarity in the cytochrome oxidase 1 gene. This is the first genome sequence for this recently described genus. The genome assembly was more fragmented and contained many duplicated genes when compared with the other Phytophthora sequences. Comparative analyses were performed with genomic sequences of the P. kernoviae isolates from the UK and New Zealand. Although the potential New Zealand origin of P. kernoviae has been suggested, new isolations from Chile had cast doubt on this hypothesis. We present evidence supporting P. kernoviae as having originated in New Zealand. However, investigation of the diversity of oomycete species in Chile has been limited and warrants further exploration. We demonstrate the expediency of genomic analyses in determining phylogenetic relationships between isolates within new and often scantly represented taxonomic groups, such as Phytophthora clade 10 and Nothophytophthora. Data are available on GenBank via BioProject accession number PRJNA352331.

Genomic, Network, and Phylogenetic Analysis of the Oomycete Effector Arsenal

The oomycetes are a class of microscopic, filamentous eukaryotes and include ecologically significant animal and plant pathogens. Oomycetes secrete large arsenals of effector proteins that degrade host cell components, manipulate host immune responses, and induce necrosis, enabling parasitic colonization. In this study, we catalogued the number and evolution of effectors in 37 oomycete species whose genomes have been completely sequenced. Large expansions of effector protein families in Phytophthora species, including glycoside hydrolases, pectinases, and necrosis-inducing proteins, were observed. Species-specific expansions were detected, including chitinases in Aphanomyces astaci and Pythium oligandrum. Novel effectors which may be involved in suppressing animal immune responses were identified in Ap. astaci and Py. oligandrum. Type 2 necrosis-inducing proteins with an unusual phylogenetic history were also located. This work represents an up-to-date in silico catalogue of the effector arsenal of the oomycetes based on the 37 genomes currently available.

The oomycetes are a class of microscopic, filamentous eukaryotes within the stramenopiles-alveolate-Rhizaria (SAR) supergroup and include ecologically significant animal and plant pathogens. Oomycetes secrete large arsenals of effector proteins that degrade host cell components, manipulate host immune responses, and induce necrosis, enabling parasitic colonization. This study investigated the expansion and evolution of effectors in 37 oomycete species in 4 oomycete orders, including Albuginales, Peronosporales, Pythiales, and Saprolegniales species. Our results highlight the large expansions of effector protein families, including glycoside hydrolases, pectinases, and necrosis-inducing proteins, in Phytophthora species. Species-specific expansions, including expansions of chitinases in Aphanomyces astaci and Pythium oligandrum, were detected. Novel effectors which may be involved in suppressing animal immune responses in Ap. astaci and Py. insidiosum were also identified. Type 2 necrosis-inducing proteins with an unusual phylogenetic history were also located in a number of oomycete species. We also investigated the "RxLR" effector complement of all 37 species and, as expected, observed large expansions in Phytophthora species numbers. Our results provide in-depth sequence information on all putative RxLR effectors from all 37 species. This work represents an up-to-date in silico catalogue of the effector arsenal of the oomycetes based on the 37 genomes currently available.

IMPORTANCE The oomycetes are a class of microscopic, filamentous eukaryotes and include ecologically significant animal and plant pathogens. Oomycetes secrete large arsenals of effector proteins that degrade host cell components, manipulate host immune responses, and induce necrosis, enabling parasitic colonization. In this study, we catalogued the number and evolution of effectors in 37 oomycete species whose genomes have been completely sequenced. Large expansions of effector protein families in Phytophthora species, including glycoside hydrolases, pectinases, and necrosis-inducing proteins, were observed. Species-specific expansions were detected, including chitinases in Aphanomyces astaci and Pythium oligandrum. Novel effectors which may be involved in suppressing animal immune responses were identified in Ap. astaci and Py. oligandrum. Type 2 necrosis-inducing proteins with an unusual phylogenetic history were also located. This work represents an up-to-date in silico catalogue of the effector arsenal of the oomycetes based on the 37 genomes currently available.

Draft genomes of two Australian strains of the plant pathogen, Phytophthora cinnamomi

Background: The oomycete plant pathogen, Phytophthora cinnamomi, is responsible for the destruction of thousands of species of native Australian plants, as well as several crops, such as avocado and macadamia, and has one of the widest host-plant ranges of the Phytophthora genus. The current reference genome of P. cinnamomi is based on an atypical strain and has large gaps in its assembly. To further studies of the pathogenicity of this species, especially in Australia, robust genome assemblies of more typical strains are required. Here we report the genome sequencing, draft assembly, and preliminary annotation of two geographically separated Australian strains of P. cinnamomi. Findings:  Some 308 million raw reads were generated for the two strains, DU054 and WA94.26. Independent genome assembly produced final genome sequences of 62.8 Mb (in 14,268 scaffolds) and 68.1 Mb (in 10,084 scaffolds), which are comparable in size and contiguity to other Phytophthora genomes. Gene prediction yielded > 22,000 predicted protein-encoding genes within each genome, while BUSCO assessment showed 94.4% and 91.5% of the stramenopile single-copy orthologs to be present in the assembled genomes, respectively. Conclusions: The assembled genomes of two geographically distant isolates of Phytophthora cinnamomi will provide a valuable resource for further comparative analyses and evolutionary studies of this destructive pathogen, and further annotation of the presented genomes may yield possible targets for novel pathogen control methods.

Draft genomes of two Australian strains of the plant pathogen, Phytophthora cinnamomi

Background: The oomycete plant pathogen, Phytophthora cinnamomi, is responsible for the destruction of thousands of species of native Australian plants, as well as several crops, such as avocado and macadamia, and has one of the widest host-plant ranges of the Phytophthora genus. The currently available genome of P. cinnamomi is based on an atypical strain and has large gaps in its assembly. To further studies of the pathogenicity of this species, especially in Australia, more robust assemblies of the genomes of more typical strains are required. Here we report the genome sequencing, draft assembly, and preliminary annotation of two geographically separated Australian strains of P. cinnamomi. Findings:  Some 308 million raw reads were generated for the two strains. Independent genome assembly produced final genomes of 62.8 Mb (in 14,268 scaffolds) and 68.1 Mb (in 10,084 scaffolds), which are comparable in size and contiguity to other Phytophthora genomes. Gene prediction yielded > 22,000 predicted protein-encoding genes within each genome, while BUSCO assessment showed 82.5% and 81.8% of the eukaryote universal single-copy orthologs to be present in the assembled genomes, respectively. Conclusions: The assembled genomes of two geographically distant isolates of Phytophthora cinnamomi will provide a valuable resource for further comparative analysis and evolutionary studies of this destructive pathogen, and further annotation of the presented genomes may yield possible targets for novel pathogen control methods.

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

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

Draft genome of the oomycete pathogen Phytophthora cactorum strain LV007 isolated from European beech (Fagus sylvatica)

Phytophthora cactorum is a broad host range phytopathogenic oomycete. P. cactorum strain LV007 was isolated from a diseased European Beech (Fagus sylvatica) in Malmö, Sweden in 2016. The draft genome of P. cactorum strain LV007 is 67.81 Mb. It contains 15,567 contigs and 21,876 predicted protein-coding genes. As reported for other phytopathogenic Phytophthora species, cytoplasmic effector proteins including RxLR and CRN families were identified. The genome sequence has been deposited at DDBJ/ENA/GenBank under the accession NBIJ00000000. The version described in this paper is version NBIJ01000000.

Development and Validation of Polymorphic Microsatellite Loci for the NA2 Lineage of Phytophthora ramorum from Whole Genome Sequence Data

Phytophthora ramorum is the causal agent of sudden oak death and sudden larch death, and is also responsible for causing ramorum blight on woody ornamental plants. Many microsatellite markers are available to characterize the genetic diversity and population structure of P. ramorum. However, only two markers are polymorphic in the NA2 lineage, which is predominant in Canadian nurseries. Microsatellite motifs were mined from whole-genome sequence data of six P. ramorum NA2 isolates. Of the 43 microsatellite primer pairs selected, 13 loci displayed different allele sizes among the four P. ramorum lineages, 10 loci displayed intralineage variation in the EU1, EU2, and/or NA1 lineages, and 12 microsatellites displayed polymorphism in the NA2 lineage. Genotyping of 272 P. ramorum NA2 isolates collected in nurseries in British Columbia, Canada, from 2004 to 2013 revealed 12 multilocus genotypes (MLGs). One MLG was dominant when examined over time and across sampling locations, and only a few mutations separated the 12 MLGs. The NA2 population observed in Canadian nurseries also showed no signs of sexual recombination, similar to what has been observed in previous studies. The markers developed in this study can be used to assess P. ramorum inter- and intralineage genetic diversity and generate a better understanding of the population structure and migration patterns of this important plant pathogen, especially for the lesser-characterized NA2 lineage.

Phylogenomic Reconstruction of the Oomycete Phylogeny Derived from 37 Genomes

The oomycetes are a class of eukaryotes and include ecologically significant animal and plant pathogens. Single-gene and multigene phylogenetic studies of individual oomycete genera and of members of the larger classes have resulted in conflicting conclusions concerning interspecies relationships among these species, particularly for the Phytophthora genus. The onset of next-generation sequencing techniques now means that a wealth of oomycete genomic data is available. For the first time, we have used genome-scale phylogenetic methods to resolve oomycete phylogenetic relationships. We used supertree methods to generate single-gene and multigene species phylogenies. Overall, our supertree analyses utilized phylogenetic data from 8,355 oomycete gene families. We have also complemented our analyses with superalignment phylogenies derived from 131 single-copy ubiquitous gene families. Our results show that a genome-scale approach to oomycete phylogeny resolves oomycete classes and clades. Our analysis represents an important first step in large-scale phylogenomic analysis of the oomycetes.

The oomycetes are a class of microscopic, filamentous eukaryotes within the Stramenopiles-Alveolata-Rhizaria (SAR) supergroup which includes ecologically significant animal and plant pathogens, most infamously the causative agent of potato blight Phytophthora infestans. Single-gene and concatenated phylogenetic studies both of individual oomycete genera and of members of the larger class have resulted in conflicting conclusions concerning species phylogenies within the oomycetes, particularly for the large Phytophthora genus. Genome-scale phylogenetic studies have successfully resolved many eukaryotic relationships by using supertree methods, which combine large numbers of potentially disparate trees to determine evolutionary relationships that cannot be inferred from individual phylogenies alone. With a sufficient amount of genomic data now available, we have undertaken the first whole-genome phylogenetic analysis of the oomycetes using data from 37 oomycete species and 6 SAR species. In our analysis, we used established supertree methods to generate phylogenies from 8,355 homologous oomycete and SAR gene families and have complemented those analyses with both phylogenomic network and concatenated supermatrix analyses. Our results show that a genome-scale approach to oomycete phylogeny resolves oomycete classes and individual clades within the problematic Phytophthora genus. Support for the resolution of the inferred relationships between individual Phytophthora clades varies depending on the methodology used. Our analysis represents an important first step in large-scale phylogenomic analysis of the oomycetes.

IMPORTANCE The oomycetes are a class of eukaryotes and include ecologically significant animal and plant pathogens. Single-gene and multigene phylogenetic studies of individual oomycete genera and of members of the larger classes have resulted in conflicting conclusions concerning interspecies relationships among these species, particularly for the Phytophthora genus. The onset of next-generation sequencing techniques now means that a wealth of oomycete genomic data is available. For the first time, we have used genome-scale phylogenetic methods to resolve oomycete phylogenetic relationships. We used supertree methods to generate single-gene and multigene species phylogenies. Overall, our supertree analyses utilized phylogenetic data from 8,355 oomycete gene families. We have also complemented our analyses with superalignment phylogenies derived from 131 single-copy ubiquitous gene families. Our results show that a genome-scale approach to oomycete phylogeny resolves oomycete classes and clades. Our analysis represents an important first step in large-scale phylogenomic analysis of the oomycetes.

Sequencing of the Litchi Downy Blight Pathogen Reveals It Is a Phytophthora Species With Downy Mildew-Like Characteristics

On the basis of its downy mildew-like morphology, the litchi downy blight pathogen was previously named Peronophythora litchii. Recently, however, it was proposed to transfer this pathogen to Phytophthora clade 4. To better characterize this unusual oomycete species and important fruit pathogen, we obtained the genome sequence of Phytophthora litchii and compared it to those from other oomycete species. P. litchii has a small genome with tightly spaced genes. On the basis of a multilocus phylogenetic analysis, the placement of P. litchii in the genus Phytophthora is strongly supported. Effector proteins predicted included 245 RxLR, 30 necrosis-and-ethylene-inducing protein-like, and 14 crinkler proteins. The typical motifs, phylogenies, and activities of these effectors were typical for a Phytophthora species. However, like the genome features of the analyzed downy mildews, P. litchii exhibited a streamlined genome with a relatively small number of genes in both core and species-specific protein families. The low GC content and slight codon preferences of P. litchii sequences were similar to those of the analyzed downy mildews and a subset of Phytophthora species. Taken together, these observations suggest that P. litchii is a Phytophthora pathogen that is in the process of acquiring downy mildew-like genomic and morphological features. Thus P. litchii may provide a novel model for investigating morphological development and genomic adaptation in oomycete pathogens.