Genomes and virulence difference between two physiological races of Phytophthora nicotianae

Complete Genome Sequence Analysis of Bacillus subtilis MC4-2 Strain That against Tobacco Black Shank Disease

The MC4-2 bacterium strain was isolated and purified from the Periplaneta americana intestine as a biocontrol agent with good antagonistic effect against the pathogens of a soil-borne disease called tobacco black shank. The MC4-2 strain was found to have good broad-spectrum inhibition by plate stand-off test. Based on 16S rRNA and gyrB genes, ANI analysis, and other comparative genomics methods, it was determined that the MC4-2 strain was Bacillus subtilis. The complete genome sequence showed that the genome size was 4,076,630 bp, the average GC content was 43.78%, and the total number of CDSs was 4,207. Genomic prediction analysis revealed that a total of 145 genes were annotated by the CAZy, containing mainly GH and CE enzymes that break down carbohydrates such as glucose, chitin, starch, and alginate, and a large number of enzymes involved in glycosylation were present. A total of ten secondary metabolite clusters were predicted, six clusters of which were annotated as surfactin, bacillaene, fengycin, bacillibactin, subtilosin A, and bacilysin. The present investigation found the biological control mechanism of B. subtilis MC4-2, which provides a strong theoretical basis for the best use of this strain in biological control methods and provides a reference for the subsequent development of agents of this bacterium.

Characterizing the Dynamics of Virulence and Fungicide Resistance of Phytophthora capsici in Michigan Vegetable Fields Reveals Loci Associated with Virulence

The oomycete Phytophthora capsici is a destructive pathogen infecting more than 50 plant species and is one of the most serious threats to cucurbit production. Phytophthora blight caused by Phytophthora capsici can affect all plant growth stages, and fungicides and cultural controls are used to limit losses. Dissecting pathogen virulence and fungicide resistance can provide insights into pathogenic mechanisms and inform effective management practices to control P. capsici. In this study, we assessed virulence, mefenoxam sensitivity, and genetic diversity of nine P. capsici populations collected from Cucurbitaceae, Solanaceae, and Fabaceae host families in Michigan from 2002 to 2016. We developed 992 simple sequence repeats (SSRs) in the P. capsici genome and identified 60 SSRs located within or close to RXLR-class (Arginine-any amino acid-Leucine-Arginine) effectors and 29 SSRs within or close to effector CRN (CRinkling and Necrosis) family protein, which represent 62 RXLR and 34 putative CRNs. Population structure analysis shows that mefenoxam resistance was not associated with the year of collection, host type, or location, but there were significant differences in virulence among the populations. Using the general linear model and mixed linear model-based association analyses with all effector-related SSR markers, we identified four SSR markers significantly associated with at least one of the virulence-related parameters. Of these, one (Pce_SC18) was in a predicted CRN effector and had high identity with the putative PhCRN37 effector in the pathogen Plasmopara halstedii, which can be further verified for virulence identification in P. capsici.

Deciphering the rhizosphere bacteriome associated with biological control of tobacco black shank disease

Introduction

The black shank disease seriously affects the health of tobacco plants. Conventional control methods have limitations in terms of effectiveness or economic aspects and cause public health concerns. Thus, biological control methods have come into the field, and microorganisms play a key role in suppressing tobacco black shank disease.

Methods

In this study, we examined the impact of soil microbial community on black shank disease basing on the structural difference of bacterial communities in rhizosphere soils. We used Illumina sequencing to compare the bacterial community diversity and structure in different rhizosphere soil samples in terms of healthy tobacco, tobacco showing typical black shank symptoms, and tobacco treated with the biocontrol agent, Bacillus velezensis S719.

Results

We found that Alphaproteobacteria in the biocontrol group, accounted for 27.2% of the ASVs, was the most abundant bacterial class among three groups. Heatmap and LEfSe analyses were done to determine the distinct bacterial genera in the three sample groups. For the healthy group, Pseudomonas was the most significant genus; for the diseased group, Stenotrophomonas exhibited the strongest enrichment trend, and Sphingomonas showed the highest linear discriminant analysis score, and was even more abundant than Bacillus; for the biocontrol group, Bacillus, and Gemmatimonas were the largely distributed genus. In addition, co-occurrence network analysis confirmed the abundance of taxa, and detected a recovery trend in the network topological parameters of the biocontrol group. Further functional prediction also provided a possible explanation for the bacterial community changes with related KEGG annotation terms.

Discussion

These findings will improve our knowledge of plant-microbe interactions and the application of biocontrol agents to improve plant fitness, and may contribute to the selection of biocontrol strains.

The Identification and Characterization of Endopolygalacturonases in a South African Isolate of Phytophthora cinnamomi

Phytophthora cinnamomi is an economically important plant pathogen that has caused devastating losses to the avocado industry worldwide. To facilitate penetration and successful colonization of the host plant, pathogens have been reported to secrete polygalacturonases (PGs). Although a large PG gene family has been reported in P. cinnamomi, in-depth bioinformatics analyses and characterization of these genes is still lacking. In this study we used bioinformatics tools and molecular biology techniques to identify and characterize endopolygalacturonases in the genome of a South African P. cinnamomi isolate, GKB4. We identified 37 PGs, with 19 characteristics of full-length PGs. Although eight PcPGs were induced in planta during infection, only three showed significant up- and down-regulation when compared with in vitro mycelial growth, suggesting their possible roles in infection. The phylogenetic analysis of PcPGs showed both gain and loss of introns in the evolution of PGs in P. cinnamomi. Furthermore, 17 PGs were related to characterized PGs from oomycete species, providing insight on possible function. This study provides new data on endoPGs in P. cinnamomi and the evolution of introns in PcPG genes. We also provide a baseline for future functional characterization of PGs suspected to contribute to P. cinnamomi pathogenicity/virulence in avocado.

Genome Resources for the Ex-Type of Phytophthora citricola, and Well-Authenticated Isolates of P. hibernalis, P. nicotianae, and P. syringae

Phytophthora is one of the most important genera of plant pathogens, with many members causing high economic losses worldwide. To build robust molecular identification systems, it is very important to have information from well-authenticated specimens and, in preference, the ex-type specimens. The reference genomes of well-authenticated specimens form a critical foundation for genetics, biological research, and diagnostic applications. In this study, we describe four draft Phytophthora genome resources for the ex-type of Phytophthora citricola BL34 (P0716 WPC) (118 contigs for 50 Mb), and well-authenticated specimens of P. syringae BL57G (P10330 WPC) (591 contigs for 75 Mb), P. hibernalis BL41G (P3822 WPC) (404 contigs for 84 Mb), and P. nicotianae BL162 (P6303 WPC) (3,984 contigs for 108 Mb) generated with MinION long-read high-throughput sequencing technology (Oxford Nanopore Technologies). Using the quality reads, we assembled high-coverage genomes of P. citricola with 291× coverage and 16,662 annotated genes; P. nicotianae with 205× coverage and 29,271 annotated genes; P. syringae with 76× coverage and 23,331 annotated genes, and P. hibernalis with 42× coverage and 21,762 annotated genes. With the availability of genome sequences and their annotations, we predict that these draft genomes will be accommodating for various basic and applied research, including diagnostics to protect global agriculture.

Impacts of Continued Exposure to a Susceptible Host Genotype on Aggressiveness of Phytophthora nicotianae Isolates Adapted to Multiple Sources of Partial Resistance

Pathogen adaptation can threaten the durability of partial resistance. Mixed plantings of susceptible and partially resistant varieties may prolong the effectiveness of partial resistance, but little is known about how continued exposure to a susceptible genotype can change the aggressiveness of pathogen isolates adapted to a source of partial resistance. The objective of this study was to examine the effects of continued exposure to a highly susceptible tobacco genotype on isolates of Phytophthora nicotianae that had been adapted to partial resistance. Isolates of P. nicotianae previously adapted to two sources of partial resistance were continually exposed to either the original host of adaptation or a susceptible host. After six generations of host exposure, isolates obtained from the partially resistant and the susceptible hosts were compared for their aggressiveness on the resistant host and for differences in expression of genes associated with pathogenicity and aggressiveness. Results suggested that exposure to the susceptible tobacco genotype reduced aggressiveness of isolates adapted to partial resistance in K 326 Wz/- but not of isolates adapted to partial resistance in Fla 301. Quantification of pathogenicity-associated gene expression using qRT-PCR suggested the rapid change in aggressiveness of isolates adapted to Wz-sourced partial resistance may have resulted from modification in gene expression in multiple genes.

Homoeologous chromosome exchange explains the creation of a QTL affecting soil-borne pathogen resistance in tobacco

Crop plant partial resistance to plant pathogens controlled by quantitative trait loci (QTL) is desirable in cultivar development programmes because of its increased durability. Mechanisms underlying such resistance are difficult to study. We performed RNA-seq analyses for tobacco (Nicotiana tabacum) nearly isogenic lines (NILs) with and without favourable allele(s) at Phn7.1, a major QTL influencing partial resistance to the soil-borne pathogens Phytophthora nicotianae and Ralstonia solanacearum. Based upon combined analyses of transcriptome-based sequence variation and gene expression profiles, we concluded that allelic variability at the Phn7.1 locus was likely generated from homoeologous exchange, which led to deletion of low-expressing members of the SAR8.2 gene family and duplication of high-expressing SAR8.2 genes from a different subgenome of allotetraploid tobacco. The high expression of endogenous Phn7.1-associated SAR8.2 genes was correlated with observed resistance to P. nicotianae. Our findings suggest a role for genomic rearrangements in the generation of favourable genetic variability affecting resistance to pathogens in plants.

Genome Analysis of Phytophthora nicotianae JM01 Provides Insights into Its Pathogenicity Mechanisms

Phytophthora nicotianae is a widely distributed plant pathogen that can cause serious disease and cause significant economic losses to various crops, including tomatoes, tobacco, onions, and strawberries. To understand its pathogenic mechanisms and explore strategies for controlling diseases caused by this pathogen, we sequenced and analyzed the whole genome of Ph. nicotianae JM01. The Ph. nicotianae JM01 genome was assembled using a combination of approaches including shotgun sequencing, single-molecule sequencing, and the Hi-C technique. The assembled Ph. nicotianae JM01 genome is about 95.32 Mb, with contig and scaffold N50 54.23 kb and 113.15 kb, respectively. The average GC content of the whole-genome is about 49.02%, encoding 23,275 genes. In addition, we identified 19.15% of interspersed elements and 0.95% of tandem elements in the whole genome. A genome-wide phylogenetic tree indicated that Phytophthora diverged from Pythium approximately 156.32 Ma. Meanwhile, we found that 252 and 285 gene families showed expansion and contraction in Phytophthora when compared to gene families in Pythium. To determine the pathogenic mechanisms Ph. nicotianae JM01, we analyzed a suite of proteins involved in plant-pathogen interactions. The results revealed that gene duplication contributed to the expansion of Cell Wall Degrading Enzymes (CWDEs) such as glycoside hydrolases, and effectors such as Arg-Xaa-Leu-Arg (RXLR) effectors. In addition, transient expression was performed on Nicotiana benthamiana by infiltrating with Agrobacterium tumefaciens cells containing a cysteine-rich (SCR) protein. The results indicated that SCR can cause symptoms of hypersensitive response. Moreover, we also conducted comparative genome analysis among four Ph. nicotianae genomes. The completion of the Ph. nicotianae JM01 genome can not only help us understand its genomic characteristics, but also help us discover genes involved in infection and then help us understand its pathogenic mechanisms.

A Novel Disease of Mung Bean, Phytophthora Stem Rot Caused by a New Forma Specialis of Phytophthora vignae

An emerging soilborne disease resembling Phytophthora stem rot was observed on mung bean plants grown in Anhui, China. To identify the causal agent, diseased plants and soil samples from 13 fields were collected to isolate the pathogen. Twenty-two Phytophthora isolates were recovered from the samples and detailed identification was conducted. Based on morphological and molecular characterizations, all of the isolates were consistently identified as P. vignae. Phylogenetic analysis using eight nuclear loci sequences of the internal transcribed spacer region, rRNA gene large subunit, a partial sequence of the β-tubulin gene, translation elongation factor 1α, 60S ribosomal protein L10, the enolase gene, heat shock protein 90, and triose phosphate isomerase/glyceraldehyde-3-phosphate dehydrogenase and a mitochondrial locus cytochrome c oxidase subunit I revealed that the mung bean isolates grouped in the same clade as P. vignae and its two formae speciales, P. vignae f. sp. adzukicola and P. vignae f. sp. vignae. A host specificity test showed that the mung bean isolates of P. vignae were pathogenic toward mung bean with a much stronger virulence and toward adzuki bean with a relatively weak virulence, but they were nonpathogenic to the other tested legume crops, including soybean, cowpea, pea, common bean, faba bean, and chickpea. The host range of mung bean isolates significantly differs from those of P. vignae f. sp. adzukicola and P. vignae f. sp. vignae based on our results and on previous studies. Thus, the pathogen causing Phytophthora stem rot of mung bean is proposed as a new forma specialis of P. vignae, designated as P. vignae f. sp. mungcola.

Predicting Lifestyle from Positive Selection Data and Genome Properties in Oomycetes

As evidenced in parasitism, host and niche shifts are a source of genomic and phenotypic diversification. Exemplary is a reduction in the core metabolism as parasites adapt to a particular host, while the accessory genome often maintains a high degree of diversification. However, selective pressures acting on the genome of organisms that have undergone recent lifestyle or host changes have not been fully investigated. Here, we developed a comparative genomics approach to study underlying adaptive trends in oomycetes, a eukaryotic phylum with a wide and diverse range of economically important plant and animal parasitic lifestyles. Our analysis reveals converging evolution on biological processes for oomycetes that have similar lifestyles. Moreover, we find that certain functions, in particular carbohydrate metabolism, transport, and signaling, are important for host and environmental adaptation in oomycetes. Given the high correlation between lifestyle and genome properties in our oomycete dataset, together with the known convergent evolution of fungal and oomycete genomes, we developed a model that predicts plant pathogenic lifestyles with high accuracy based on functional annotations. These insights into how selective pressures correlate with lifestyle may be crucial to better understand host/lifestyle shifts and their impact on the genome.

Diverse Trajectories Drive the Expression of a Giant Virus in the Oomycete Plant Pathogen Phytophthora parasitica

Giant viruses of amoebas, recently classified in the class Megaviricetes, are a group of viruses that can infect major eukaryotic lineages. We previously identified a set of giant virus sequences in the genome of Phytophthora parasitica, an oomycete and a devastating major plant pathogen. How viral insertions shape the structure and evolution of the invaded genomes is unclear, but it is known that the unprecedented functional potential of giant viruses is the result of an intense genetic interplay with their hosts. We previously identified a set of giant virus sequences in the genome of P. parasitica, an oomycete and a devastating major plant pathogen. Here, we show that viral pieces are found in a 550-kb locus and are organized in three main clusters. Viral sequences, namely RNA polymerases I and II and a major capsid protein, were identified, along with orphan sequences, as a hallmark of giant viruses insertions. Mining of public databases and phylogenetic reconstructions suggest an ancient association of oomycetes and giant viruses of amoeba, including faustoviruses, African swine fever virus (ASFV) and pandoraviruses, and that a single viral insertion occurred early in the evolutionary history of oomycetes prior to the Phytophthora–Pythium radiation, estimated at ∼80 million years ago. Functional annotation reveals that the viral insertions are located in a gene sparse region of the Phytophthora genome, characterized by a plethora of transposable elements (TEs), effectors and other genes potentially involved in virulence. Transcription of viral genes was investigated through analysis of RNA-Seq data and qPCR experiments. We show that most viral genes are not expressed, and that a variety of mechanisms, including deletions, TEs insertions and RNA interference may contribute to transcriptional repression. However, a gene coding a truncated copy of RNA polymerase II along a set of neighboring sequences have been shown to be expressed in a wide range of physiological conditions, including responses to stress. These results, which describe for the first time the endogenization of a giant virus in an oomycete, contribute to challenge our view of Phytophthora evolution.

Comparative analysis of Phytophthora genomes reveals oomycete pathogenesis in crops

The oomycete genus Phytophthora includes devastating plant pathogens that are found in almost all ecosystems. We sequenced the genomes of two quarantined Phytophthora species–P. fragariae and P. rubi. Comparing these Phytophthora species and related genera allowed reconstruction of the phylogenetic relationships within the genus Phytophthora and revealed Phytophthora genomic features associated with infection and pathogenicity. We found that several hundred Phytophthora genes are putatively inherited from red algae, but Phytophthora does not have vestigial plastids originating from phototrophs. The horizontally-transferred Phytophthora genes are abundant transposons that “transmit” exogenous gene to Phytophthora species thus bring about the gene recombination possibility. Several expansion events of Phytophthora gene families associated with cell wall biogenesis can be used as mutational targets to elucidate gene function in pathogenic interactions with host plants. This work enhanced the understanding of Phytophthora evolution and will also be helpful for the design of phytopathological control strategies.

Comparative transcriptome analysis reveals resistant and susceptible genes in tobacco cultivars in response to infection by Phytophthora nicotianae

Phytophthora nicotianae is highly pathogenic to Solanaceous crops and is a major problem in tobacco production. The tobacco cultivar Beihart1000-1 (BH) is resistant, whereas the Xiaohuangjin 1025 (XHJ) cultivar is susceptible to infection. Here, BH and XHJ were used as models to identify resistant and susceptible genes using RNA sequencing (RNA-seq). Roots were sampled at 0, 6, 12, 24, and 60 h post infection. In total, 23,753 and 25,187 differentially expressed genes (DEGs) were identified in BH and XHJ, respectively. By mapping upregulated DEGs to the KEGG database, changes of the rich factor of "plant pathogen interaction pathway" were corresponded to the infection process. Of all the DEGs in this pathway, 38 were specifically regulated in BH. These genes included 11 disease-resistance proteins, 3 pathogenesis-related proteins, 4 RLP/RLKs, 2 CNGCs, 7 calcium-dependent protein kinases, 4 calcium-binding proteins, 1 mitogen-activated protein kinase kinase, 1 protein EDS1L, 2 WRKY transcription factors, 1 mannosyltransferase, and 1 calmodulin-like protein. By combining the analysis of reported susceptible (S) gene homologs and DEGs in XHJ, 9 S gene homologs were identified, which included 1 calmodulin-binding transcription activator, 1 cyclic nucleotide-gated ion channel, 1 protein trichome birefringence-like protein, 1 plant UBX domain-containing protein, 1 ADP-ribosylation factor GTPase-activating protein, 2 callose synthases, and 2 cellulose synthase A catalytic subunits. qRT-PCR was used to validate the RNA-seq data. The comprehensive transcriptome dataset described here, including candidate resistant and susceptible genes, will provide a valuable resource for breeding tobacco plants resistant to P. nicotianae infections.

Improved Genome Sequence and Gene Annotation Resource for the Potato Late Blight Pathogen Phytophthora infestans

Phytophthora infestans is a devastating pathogen causing potato late blight (Solanum tuberosum). Here we report the sequencing, assembly and genome annotation for two Phytophthora infestans isolates sampled in Republic of Korea. Genome sequencing was carried out using long read (Oxford Nanopore) and short read (Illumina Nextseq) sequencing technologies that significantly improved the contiguity and quality of P. infestans genome assembly. Our resources would help researchers better understand the molecular mechanisms by which P. infestans causes late blight disease in the future.

Characterization of Two Satellite DNA Families in the Genome of the Oomycete Plant Pathogen Phytophthora parasitica

Satellite DNA is a class of repetitive sequences that are organized in long arrays of tandemly repeated units in most eukaryotes. Long considered as selfish DNA, satellite sequences are now proposed to contribute to genome integrity. Despite their potential impact on the architecture and evolution of the genome, satellite DNAs have not been investigated in oomycetes due to the paucity of genomic data and the difficulty of assembling highly conserved satellite arrays. Yet gaining knowledge on the structure and evolution of genomes of oomycete pathogens is crucial to understanding the mechanisms underlying adaptation to their environment and to proposing efficient disease control strategies. A de novo assembly of the genome of Phytophthora parasitica, an important oomycete plant pathogen, led to the identification of several families of tandemly repeated sequences varying in size, copy number, and sequence conservation. Among them, two abundant families, designated as PpSat1 and PpSat2, displayed typical features of satellite DNA and were collectively designated as PpSat. These two satellite families differ by their length, sequence, organization, genomic environment, and evolutionary dynamics. PpSat1, but not PpSat2, presented homologs among oomycetes. This observation, as well as the characterization of transcripts of PpSat families, suggested that these satellite DNA families likely play a conserved role within this important group of pathogens.

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.

Evaluation of High-Resolution Melting for Rapid Differentiation of Phytophthora Hybrids and Their Parental Species

Phytophthora species hybrids have been repeatedly reported as causing damaging diseases to cultivated and wild plants. Two known hybrids, P. andina and P. × pelgrandis, are pathogens of Solanaceae and ornamentals, respectively, although the extent of their host ranges are unknown. P. andina emerged from hybridization of P. infestans and an unidentified related species, whereas P. × pelgrandis emerged from P. nicotianae and P. cactorum. Considering that hybrids and parental species can coexist in the same regions and to distinguish them usually requires cloning or whole genome sequencing, we aimed to develop a rapid tool to distinguish them. Specifically, we used high-resolution melting (HRM) assays to differentiate genotypes based on their amplicon melting profiles. We designed primers for P. × pelgrandis and parental species based on available sequences of P. nicotianae and P. cactorum nuclear genes containing polymorphisms between species. For P. andina, heterozygous sites from Illumina short reads were used for the same purpose. We identified multiple amplicons exhibiting differences in melting curves between parental species and hybrids. We propose HRM as a rapid method for differentiation of P. andina and P. × pelgrandis hybrids from parental species that could be employed to advance research on these pathogens.

Effectors of Phytophthora pathogens are powerful weapons for manipulating host immunity

This article provides an overview of the interactions between Phytophthora effectors and plant immune system components, which form a cross-linked complex network that regulates plant pathogen resistance. Pathogens secrete numerous effector proteins into plants to promote infections. Several Phytophthora species (e.g., P. infestans, P. ramorum, P. sojae, P. capsici, P. cinnamomi, and P. parasitica) are notorious pathogens that are extremely damaging to susceptible plants. Analyses of genomic data revealed that Phytophthora species produce a large group of effector proteins, which are critical for pathogenesis. And, the targets and functions of many identified Phytophthora effectors have been investigated. Phytophthora effectors can affect various aspects of plant immune systems, including plant cell proteases, phytohormones, RNAs, the MAPK pathway, catalase, the ubiquitin proteasome pathway, the endoplasmic reticulum, NB-LRR proteins, and the cell membrane. Clarifying the effector-plant interactions is important for unravelling the functions of Phytophthora effectors during pathogenesis. In this article, we review the effectors identified in recent decades and provide an overview of the effector-directed regulatory network in plants following infections by Phytophthora species.

The Nuclear-Localized RxLR Effector PvAvh74 From Plasmopara viticola Induces Cell Death and Immunity Responses in Nicotiana benthamiana

Downy mildew is one of the most serious diseases of grapevine (Vitis spp). The causal agent of grapevine downy mildew, Plasmopara viticola, is an obligate biotrophic oomycete. Although oomycete pathogens such as P. viticola are known to secrete RxLR effectors to manipulate host immunity, there have been few studies of the associated mechanisms by which these may act. Here, we show that a candidate P. viticola RxLR effector, PvAvh74, induces cell death in Nicotiana benthamiana leaves. Using agroinfiltration, we found that nuclear localization, two putative N-glycosylation sites, and 427 amino acids of the PvAvh74 carboxyl terminus were necessary for cell-death-inducing activity. Using virus-induced gene silencing (VIGS), we found that PvAvh74-induced cell death in N. benthamiana requires EDS1, NDR1, SGT1, RAR1, and HSP90, but not BAK1. The MAPK cascade components MEK2, WIPK, and SIPK were also involved in PvAvh74-induced cell death in N. benthamiana. Transient expression of PvAvh74 could suppress Phytophthora capsici colonization of N. benthamiana, which suggests that PvAvh74 elicits plant immune responses. Suppression of P. capsici colonization also was dependent on nuclear localization of PvAvh74. Additionally, PvAvh74-triggered cell death could be suppressed by another effector, PvAvh8, from the same isolate. This work provides a framework to further investigate the interactions of PvAvh74 and other RxLR effectors with host immunity.

A High-Quality Grapevine Downy Mildew Genome Assembly Reveals Rapidly Evolving and Lineage-Specific Putative Host Adaptation Genes

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high contiguity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant–pathogen interactions. Plasmopara viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.

A High-Quality Grapevine Downy Mildew Genome Assembly Reveals Rapidly Evolving and Lineage-Specific Putative Host Adaptation Genes

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high contiguity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant-pathogen interactions. Plasmopara viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.

Study of the whole genome, methylome and transcriptome of Cordyceps militaris

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

Characterization of the Black Shank Pathogen, Phytophthora nicotianae, Across North Carolina Tobacco Production Areas

Black shank disease of tobacco, caused by the oomycete Phytophthora nicotianae, is a major threat to production in the United States and tobacco-producing areas worldwide. In a statewide survey of North Carolina, the rapid shift from race 0 to race 1 was documented. Collected pathogen isolates were characterized phenotypically for mating type and mefenoxam sensitivity, and genotypically by comparing sequences from three cytoplasmic and two nuclear regions. Both the A1 and A2 mating types were found throughout the state. When both mating types were recovered from the same field, pairings of isolates yielded viable oospores, indicating for the first time the potential for sexual sporulation by P. nicotianae in natural populations. Because the loss of complete resistance required a renewed use of the fungicide mefenoxam, a subset of the survey isolates was screened for sensitivity to the fungicide. All isolates were sensitive, with a mean effective concentration to inhibit 50% of hyphal growth of 0.4 μg/ml that was similar across mating types and races. Molecular characterization of 226 isolates revealed that the pathogen exists as multiple clonal types within the state. Genetic diversity among the pathogen population and the potential for sexual recombination may help explain the ability of the pathogen to rapidly adapt to host resistance genes.

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.

Complete Mitochondrial Genome of Phytophthora nicotianae and Identification of Molecular Markers for the Oomycetes

Phytophthora nicotianae is one of the most destructive plant pathogens affecting a variety of plants, causing black shank of tobacco, among several other devastating diseases. Herein, we assembled the mitochondrial genome of P. nicotianae and analyzed its gene content and genome structure, performed comparative mitochondrial genomics analysis, and assessed phylogenetic relationships among oomycetes species. The circular mitogenome is 37,561 bp long, with 38 protein-coding genes, 25 transfer RNA (tRNA) genes, and 2 ribosomal RNA genes (rrnl and rrns). The mitochondrial genome showed a biased A/T usage versus G/C. The overall gene content and size of the P. nicotianae mitogenome are identical to those of other published Phytophthora mitogenomes. Interestingly, collinearity analysis using an existing ∼10 k inversion region (including 11 genes and 8 tRNAs) revealed that Phytophthora andina, Phytophthora infestans, Phytophthora mirabilis, Phytophthora ipomoeae, and Phytophthora phaseoli differed from Phytophthora nicotianae, Phytophthora sojae, Phytophthora ramorum, and Phytophthora polonica. Moreover, inverted repeat regions were found to be absent among species of the Peronosporales when compared with species from the Pythiales and Saprolegniales. A phylogenomic investigation based on 29 protein-coding genes demonstrated that Phytophthora is monophyletic, and placed P. nicotianae close to the clade including P. mirabilis, P. ipomoeae, P. andina, P. infestans, and P. phaseoli. Furthermore, we discovered six new candidate DNA molecular markers (rpl6, atp8, nad11, rps2, rps3, and rps4) based on these mitogenomes that would be suitable for species identification in the oomycetes, which have the same identification level as the whole mitogenome and ribosomal DNA sequences. These new molecular markers can not only provide a quick preview of the species without mitogenome information, but will also help to gain better understanding of the oomycetes pathogens and developing treatment or monitoring strategies.

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.

Narrow-Leafed Lupin (Lupinus angustifolius) β1- and β6-Conglutin Proteins Exhibit Antifungal Activity, Protecting Plants against Necrotrophic Pathogen Induced Damage from Sclerotinia sclerotiorum and Phytophthora nicotianae

Vicilins (7S globulins) are seed storage proteins and constitute the main protein family in legume seeds, particularly in narrow-leafed lupin (Lupinus angustifolius L.; NLL), where seven vicilin genes, called β1- to β7-conglutin have been identified. Vicilins are involved in germination processes supplying amino acids for seedling growth and plant development, as well as in some cases roles in plant defense and protection against pathogens. The roles of NLL β-conglutins in plant defense are unknown. Here the potential role of five NLL β-conglutin family members in protection against necrotrophic fungal pathogens was investigated and it was demonstrated that recombinant purified 6xHis-tagged β1- and β6-conglutin proteins exhibited the strongest in vitro growth inhibitory activity against a range of necrotrophic fungal pathogens compared to β2, β3, and β4 conglutins. To examine activity in vivo, two representative necrotrophic pathogens, the fungus Sclerotinia sclerotiorum and oomycete Phytophthora nicotianae were used. Transient expression of β1- and β6-conglutin proteins in Nicotiana benthamiana leaves demonstrated in vivo growth suppression of both of these pathogens, resulting in low percentages of hyphal growth and elongation in comparison to control treated leaves. Cellular studies using β1- and β6-GFP fusion proteins showed these conglutins localized to the cell surface including plasmodesmata. Analysis of cellular death following S. sclerotiorum or P. nicotianae revealed both β1- and β6-conglutins suppressed pathogen induced cell death in planta and prevented pathogen induced suppression of the plant oxidative burst as determined by protein oxidation in infected compared to mock-inoculated leaves.