PhytoPipe: a phytosanitary pipeline for plant pathogen detection and diagnosis using RNA-seq data.

BACKGROUND: Detection of exotic plant pathogens and preventing their entry and establishment are critical for the protection of agricultural systems while securing the global trading of agricultural commodities. High-throughput sequencing (HTS) has been applied successfully for plant pathogen discovery, leading to its current application in routine pathogen detection. However, the analysis of massive amounts of HTS data has become one of the major challenges for the use of HTS more broadly as a rapid diagnostics tool. Several bioinformatics pipelines have been developed to handle HTS data with a focus on plant virus and viroid detection. However, there is a need for an integrative tool that can simultaneously detect a wider range of other plant pathogens in HTS data, such as bacteria (including phytoplasmas), fungi, and oomycetes, and this tool should also be capable of generating a comprehensive report on the phytosanitary status of the diagnosed specimen. RESULTS: We have developed an open-source bioinformatics pipeline called PhytoPipe (Phytosanitary Pipeline) to provide the plant pathology diagnostician community with a user-friendly tool that integrates analysis and visualization of HTS RNA-seq data. PhytoPipe includes quality control of reads, read classification, assembly-based annotation, and reference-based mapping. The final product of the analysis is a comprehensive report for easy interpretation of not only viruses and viroids but also bacteria (including phytoplasma), fungi, and oomycetes. PhytoPipe is implemented in Snakemake workflow with Python 3 and bash scripts in a Linux environment. The source code for PhytoPipe is freely available and distributed under a BSD-3 license. CONCLUSIONS: PhytoPipe provides an integrative bioinformatics pipeline that can be used for the analysis of HTS RNA-seq data. PhytoPipe is easily installed on a Linux or Mac system and can be conveniently used with a Docker image, which includes all dependent packages and software related to analyses. It is publicly available on GitHub at https://github.com/healthyPlant/PhytoPipe and on Docker Hub at https://hub.docker.com/r/healthyplant/phytopipe .

Novel Echarate Virus Variant Isolated from Patient with Febrile Illness, Chanchamayo, Peru.

A new phlebovirus variant was isolated from an acute febrile patient in Chanchamayo, Peru. Genome characterization and p-distance analyses based on complete open reading frames revealed that the virus is probably a natural reassortant of the Echarate virus (large and small segments) with a yet-unidentified phlebovirus (M segment).

Complete genome sequence of zoysia mosaic virus, a novel member of the genus Poacevirus.

Here, we report the detection and characterization of the genome of a novel poacevirus isolated from Zoysia matrella (Merrill) imported into the United States from Japan. The novel virus, tentatively named "zoysia mosaic virus" (ZoMV), is a single-stranded RNA virus with a genome of 9,728 nucleotides (nt) in length, encoding a large putative polyprotein of 3,119 amino acids (aa). The ZoMV genome is closely related to the triticum mosaic virus (TriMV; FJ263671) genome, with 57.18% nt and 51.74% aa sequence identity in the polyprotein region. Moreover, phylogenetic analysis showed that ZoMV is closely related to all other members of the genus Poacevirus. A survey of imported grasses showed that ZoMV was detected only in zoysiagrass. This is the first report of the complete genome sequence of a novel viral pathogen of zoysiagrass of the genus Poacevirus, for which we propose the binomial species name "Poacevirus zoisiae".

Telomere-to-telomere genome assembly of the aflatoxin biocontrol agent Aspergillus flavus isolate La3279 isolated from maize in Nigeria.

Here, we report the complete genome of the non-aflatoxigenic Aspergillus flavus isolate La3279, which is an active ingredient of the aflatoxin biocontrol product Aflasafe. The chromosome-scale assembly clarifies the deletion pattern in the aflatoxin biosynthesis gene cluster and corrects a misidentified assembly previously published for this isolate.

Virome profiling of fig wasps (Ceratosolen spp.) reveals virus diversity spanning four realms.

We investigated the virome of agaonid fig wasps (Ceratosolen spp.) inside syconia ("fruits") of various Ficus trees fed upon by frugivores such as pteropodid bats in Sub-Saharan Africa. This virome includes representatives of viral families spanning four realms and includes near-complete genome sequences of three novel viruses and fragments of five additional potentially novel viruses evolutionarily associated with insects, fungi, plants, and vertebrates. Our study provides evidence that frugivorous animals are exposed to a plethora of viruses by coincidental consumption of fig wasps, which are obligate pollinators of figs worldwide.

A Recombinant Chimeric Cedar Virus-Based Surrogate Neutralization Assay Platform for Pathogenic Henipaviruses.

The henipaviruses, Nipah virus (NiV), and Hendra virus (HeV) can cause fatal diseases in humans and animals, whereas Cedar virus is a nonpathogenic henipavirus. Here, using a recombinant Cedar virus (rCedV) reverse genetics platform, the fusion (F) and attachment (G) glycoprotein genes of rCedV were replaced with those of NiV-Bangladesh (NiV-B) or HeV, generating replication-competent chimeric viruses (rCedV-NiV-B and rCedV-HeV), both with and without green fluorescent protein (GFP) or luciferase protein genes. The rCedV chimeras induced a Type I interferon response and utilized only ephrin-B2 and ephrin-B3 as entry receptors compared to rCedV. The neutralizing potencies of well-characterized cross-reactive NiV/HeV F and G specific monoclonal antibodies against rCedV-NiV-B-GFP and rCedV-HeV-GFP highly correlated with measurements obtained using authentic NiV-B and HeV when tested in parallel by plaque reduction neutralization tests (PRNT). A rapid, high-throughput, and quantitative fluorescence reduction neutralization test (FRNT) using the GFP-encoding chimeras was established, and monoclonal antibody neutralization data derived by FRNT highly correlated with data derived by PRNT. The FRNT assay could also measure serum neutralization titers from henipavirus G glycoprotein immunized animals. These rCedV chimeras are an authentic henipavirus-based surrogate neutralization assay that is rapid, cost-effective, and can be utilized outside high containment.

Genomic and virologic characterization of samples from a shipboard outbreak of COVID-19 reveals distinct variants within limited temporospatial parameters.

Early in the pandemic, in March of 2020, an outbreak of COVID-19 occurred aboard the aircraft carrier USS Theodore Roosevelt (CVN-71), during deployment in the Western Pacific. Out of the crew of 4,779 personnel, 1,331 service members were suspected or confirmed to be infected with SARS-CoV-2. The demographic, epidemiologic, and laboratory findings of service members from subsequent investigations have characterized the outbreak as widespread transmission of virus with relatively mild symptoms and asymptomatic infection among mostly young healthy adults. At the time, there was no available vaccination against COVID-19 and there was very limited knowledge regarding SARS-CoV-2 mutation, dispersal, and transmission patterns among service members in a shipboard environment. Since that time, other shipboard outbreaks from which data can be extracted have occurred, but these later shipboard outbreaks have occurred largely in settings where the majority of the crew were vaccinated, thereby limiting spread of the virus, shortening duration of the outbreaks, and minimizing evolution of the virus within those close quarters settings. On the other hand, since the outbreak on the CVN-71 occurred prior to widespread vaccination, it continued over the course of roughly two months, infecting more than 25% of the crew. In order to better understand genetic variability and potential transmission dynamics of COVID-19 in a shipboard environment of immunologically naïve, healthy individuals, we performed whole-genome sequencing and virus culture from eighteen COVID-19-positive swabs collected over the course of one week. Using the unique variants identified in those genomes, we detected seven discrete groups of individuals within the population aboard CVN-71 infected with viruses of distinct genomic signature. This is in stark contrast to a recent outbreak aboard another U.S. Navy ship with >98% vaccinated crew after a port visit in Reykjavik, Iceland, where the outbreak lasted only approximately 2 weeks and the virus was clonal. Taken together, these results demonstrate the utility of sequencing from complex clinical samples for molecular epidemiology and they also suggest that a high rate of vaccination among a population in close communities may greatly reduce spread, thereby restricting evolution of the virus.

Conservation and Loss of a Putative Iron Utilization Gene Cluster among Genotypes of Aspergillus flavus.

Iron is an essential component for growth and development. Despite relative abundance in the environment, bioavailability of iron is limited due to oxidation by atmospheric oxygen into insoluble ferric iron. Filamentous fungi have developed diverse pathways to uptake and use iron. In the current study, a putative iron utilization gene cluster (IUC) in Aspergillus flavus was identified and characterized. Gene analyses indicate A. flavus may use reductive as well as siderophore-mediated iron uptake and utilization pathways. The ferroxidation and iron permeation process, in which iron transport depends on the coupling of these two activities, mediates the reductive pathway. The IUC identified in this work includes six genes and is located in a highly polymorphic region of the genome. Diversity among A. flavus genotypes is manifested in the structure of the IUC, which ranged from complete deletion to a region disabled by multiple indels. Molecular profiling of A. flavus populations suggests lineage-specific loss of IUC. The observed variation among A. flavus genotypes in iron utilization and the lineage-specific loss of the iron utilization genes in several A. flavus clonal lineages provide insight on evolution of iron acquisition and utilization within Aspergillus section Flavi. The potential divergence in capacity to acquire iron should be taken into account when selecting A. flavus active ingredients for biocontrol in niches where climate change may alter iron availability.

Immunological and Genetic Investigation of SARS-CoV-2 Reinfection in an Otherwise Healthy, Young Marine Recruit.

We used epidemiologic and viral genetic information to identify a case of likely reinfection in an otherwise healthy, young Marine recruit enrolled in the prospective, longitudinal COVID-19 Health Action Response for Marines (CHARM) study, and we paired these findings with serological studies. This participant had a positive RT-PCR to SARS-CoV-2 upon routine sampling on study day 7, although he was asymptomatic at that time. He cleared the infection within seven days. On study day 46, he had developed symptoms consistent with COVID-19 and tested positive by RT-PCR for SARS-CoV-2 again. Viral whole genome sequencing was conducted from nares swabs at multiple time points. The day 7 sample was determined to be lineage B.1.340, whereas both the day 46 and day 49 samples were B.1.1. The first positive result for anti-SARS-CoV-2 IgM serology was collected on day 49 and for IgG on day 91. This case appears most consistent with a reinfection event. Our investigation into this case is unique in that we compared sequence data from more than just paired specimens, and we also assayed for immune response after both the initial infection and the later reinfection. These data demonstrate that individuals who have experienced an infection with SARS-CoV-2 may fail to generate effective or long-lasting immunity, similar to endemic human beta coronaviruses.

Effect of slow desiccation and freezing on gene transcription and stress survival of an Antarctic nematode.

Nematodes are the dominant soil animals of the Antarctic Dry Valleys and are capable of surviving desiccation and freezing in an anhydrobiotic state. Genes induced by desiccation stress have been successfully enumerated in nematodes; however, little is known about gene regulation by Antarctic nematodes that can survive multiple types and incidences of environmental stress. In order to reveal the molecular response of anhydrobiotic survival, we investigated the genetic response of an Antarctic nematode, Plectus murrayi, which can survive desiccation and freezing. Using molecular approaches, we provide insight into the regulation of desiccation-induced transcripts during different stages of stress survival under conditions characteristic of the Antarctic Dry Valley environment. Our study revealed that exposure to slow desiccation and freezing plays an important role in the transcription of stress, metabolism and signal transduction-related genes and improves desiccation and freezing survival compared with nematodes exposed to fast desiccation and freezing. Temporal analyses of gene expression showed that pre-exposure to mild stress promotes survival of harsher stress. Our results further show that exposure to slow dehydration not only improves extreme desiccation survival but also promotes enhanced cold tolerance. We also provide evidence that slow dehydration can enhance freeze tolerance in an Antarctic nematode. Expression profiling of P. murrayi transcripts is an important step in understanding the genome-level response of this nematode to different environmental stressors.

Desiccation survival in an Antarctic nematode: molecular analysis using expressed sequenced tags.

BACKGROUND: Nematodes are the dominant soil animals in Antarctic Dry Valleys and are capable of surviving desiccation and freezing in an anhydrobiotic state. Genes induced by desiccation stress have been successfully enumerated in nematodes; however we have little knowledge of gene regulation by Antarctic nematodes which can survive multiple environmental stresses. To address this problem we investigated the genetic responses of a nematode species, Plectus murrayi, that is capable of tolerating Antarctic environmental extremes, in particular desiccation and freezing. In this study, we provide the first insight into the desiccation induced transcriptome of an Antarctic nematode through cDNA library construction and suppressive subtractive hybridization. RESULTS: We obtained 2,486 expressed sequence tags (ESTs) from 2,586 clones derived from the cDNA library of desiccated P. murrayi. The 2,486 ESTs formed 1,387 putative unique transcripts of which 523 (38%) had matches in the model-nematode Caenorhabditis elegans, 107 (7%) in nematodes other than C. elegans, 153 (11%) in non-nematode organisms and 605 (44%) had no significant match to any sequences in the current databases. The 1,387 unique transcripts were functionally classified by using Gene Ontology (GO) hierarchy and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The results indicate that the transcriptome contains a group of transcripts from diverse functional areas. The subtractive library of desiccated nematodes showed 80 transcripts differentially expressed during desiccation stress, of which 28% were metabolism related, 19% were involved in environmental information processing, 28% involved in genetic information processing and 21% were novel transcripts. Expression profiling of 14 selected genes by quantitative Real-time PCR showed 9 genes significantly up-regulated, 3 down-regulated and 2 continuously expressed in response to desiccation. CONCLUSION: The establishment of a desiccation EST collection for Plectus murrayi, a useful model in assessing the structural, physiological, biochemical and genetic aspects of multiple stress tolerance, is an important step in understanding the genome level response of this nematode to desiccation stress. The type of transcript analysis performed in this study sets the foundation for more detailed functional and genome level analyses of the genes involved in desiccation tolerance in nematodes.

Transcriptional profiling of trait deterioration in the insect pathogenic nematode Heterorhabditis bacteriophora.

BACKGROUND: The success of a biological control agent depends on key traits, particularly reproductive potential, environmental tolerance, and ability to be cultured. These traits can deteriorate rapidly when the biological control agent is reared in culture. Trait deterioration under laboratory conditions has been widely documented in the entomopathogenic nematode (EPN) Heterorhabditis bacteriophora (Hb) but the specific mechanisms behind these genetic processes remain unclear. This research investigates the molecular mechanisms of trait deterioration of two experimental lines of Hb, an inbred line (L5M) and its original parental line (OHB). We generated transcriptional profiles of two experimental lines of Hb, identified the differentially expressed genes (DEGs) and validated their differential expression in the deteriorated line. RESULTS: An expression profiling study was performed between experimental lines L5M and OHB of Hb with probes for 15,220 ESTs from the Hb transcriptome. Microarray analysis showed 1,185 DEGs comprising of 469 down- and 716 up-regulated genes in trait deteriorated nematodes. Analysis of the DEGs showed that trait deterioration involves massive changes of the transcripts encoding enzymes involved in metabolism, signal transduction, virulence and longevity. We observed a pattern of reduced expression of enzymes related to primary metabolic processes and induced secondary metabolism. Expression of sixteen DEGs in trait deteriorated nematodes was validated by quantitative reverse transcription-PCR (qRT-PCR) which revealed similar expression kinetics for all the genes tested as shown by microarray. CONCLUSION: As the most closely related major entomopathogen to C. elegans, Hb provides an attractive near-term application for using a model organism to better understand interspecies interactions and to enhance our understanding of the mechanisms underlying trait deterioration in biological control agents. This information could also be used to improve the beneficial traits of biological control agents and better understand fundamental aspects of nematode parasitism and mutualism.

Monitoring Aspergillus flavus Genotypes in a Multi-Genotype Aflatoxin Biocontrol Product With Quantitative Pyrosequencing.

Aflatoxins pose significant food security and public health risks, decrease productivity and profitability of animal industries, and hamper trade. To minimize aflatoxin contamination in several crops, a biocontrol technology based on atoxigenic strains of Aspergillus flavus is commercially used in the United States and some African countries. Significant efforts are underway to popularize the use of biocontrol in Africa by various means including incentives. The purpose of this study was to develop quantitative pyrosequencing assays for rapid, simultaneous quantification of proportions of four A. flavus biocontrol genotypes within complex populations of A. flavus associated with maize crops in Nigeria to facilitate payment of farmer incentives for Aflasafe (a biocontrol product) use. Protocols were developed to confirm use of Aflasafe by small scale farmers in Nigeria. Nested PCR amplifications followed by sequence by synthesis pyrosequencing assays were required to quantify frequencies of the active ingredients and, in so doing, confirm successful use of biocontrol by participating farmers. The entire verification process could be completed in 3-4 days proving a savings over other monitoring methods in both time and costs and providing data in a time frame that could work with the commercial agriculture scheme. Quantitative pyrosequencing assays represent a reliable tool for rapid detection, quantification, and monitoring of multiple A. flavus genotypes within complex fungal communities, satisfying the requirements of the regulatory community and crop end-users that wish to determine which purchased crops were treated with the biocontrol product. Techniques developed in the current study can be modified for monitoring other crop-associated fungi.

Degeneration of aflatoxin gene clusters in Aspergillus flavus from Africa and North America.

Aspergillus flavus is the most common causal agent of aflatoxin contamination of food and feed. However, aflatoxin-producing potential varies widely among A. flavus genotypes with many producing no aflatoxins. Some non-aflatoxigenic genotypes are used as biocontrol agents to prevent contamination. Aflatoxin biosynthesis genes are tightly clustered in a highly conserved order. Gene deletions and presence of single nucleotide polymorphisms (SNPs) in aflatoxin biosynthesis genes are often associated with A. flavus inability to produce aflatoxins. In order to identify mechanisms of non-aflatoxigenicity in non-aflatoxigenic genotypes of value in aflatoxin biocontrol, complete cluster sequences of 35 A. flavus genotypes from Africa and North America were analyzed. Inability of some genotypes to produce aflatoxin resulted from deletion of biosynthesis genes. In other genotypes, non-aflatoxigenicity originated from SNP formation. The process of degeneration differed across the gene cluster; genes involved in early biosynthesis stages were more likely to be deleted while genes involved in later stages displayed high frequencies of SNPs. Comparative analyses of aflatoxin gene clusters provides insight into the diversity of mechanisms of non-aflatoxigenicity in A. flavus genotypes used as biological control agents. The sequences provide resources for both diagnosis of non-aflatoxigenicity and monitoring of biocontrol genotypes during biopesticide manufacture and in the environment.

Comparative genomics reveals insight into virulence strategies of plant pathogenic oomycetes.

The kingdom Stramenopile includes diatoms, brown algae, and oomycetes. Plant pathogenic oomycetes, including Phytophthora, Pythium and downy mildew species, cause devastating diseases on a wide range of host species and have a significant impact on agriculture. Here, we report comparative analyses on the genomes of thirteen straminipilous species, including eleven plant pathogenic oomycetes, to explore common features linked to their pathogenic lifestyle. We report the sequencing, assembly, and annotation of six Pythium genomes and comparison with other stramenopiles including photosynthetic diatoms, and other plant pathogenic oomycetes such as Phytophthora species, Hyaloperonospora arabidopsidis, and Pythium ultimum var. ultimum. Novel features of the oomycete genomes include an expansion of genes encoding secreted effectors and plant cell wall degrading enzymes in Phytophthora species and an over-representation of genes involved in proteolytic degradation and signal transduction in Pythium species. A complete lack of classical RxLR effectors was observed in the seven surveyed Pythium genomes along with an overall reduction of pathogenesis-related gene families in H. arabidopsidis. Comparative analyses revealed fewer genes encoding enzymes involved in carbohydrate metabolism in Pythium species and H. arabidopsidis as compared to Phytophthora species, suggesting variation in virulence mechanisms within plant pathogenic oomycete species. Shared features between the oomycetes and diatoms revealed common mechanisms of intracellular signaling and transportation. Our analyses demonstrate the value of comparative genome analyses for exploring the evolution of pathogenesis and survival mechanisms in the oomycetes. The comparative analyses of seven Pythium species with the closely related oomycetes, Phytophthora species and H. arabidopsidis, and distantly related diatoms provide insight into genes that underlie virulence.

Carbohydrate-active enzymes in pythium and their role in plant cell wall and storage polysaccharide degradation.

Carbohydrate-active enzymes (CAZymes) are involved in the metabolism of glycoconjugates, oligosaccharides, and polysaccharides and, in the case of plant pathogens, in the degradation of the host cell wall and storage compounds. We performed an in silico analysis of CAZymes predicted from the genomes of seven Pythium species (Py. aphanidermatum, Py. arrhenomanes, Py. irregulare, Py. iwayamai, Py. ultimum var. ultimum, Py. ultimum var. sporangiiferum and Py. vexans) using the "CAZymes Analysis Toolkit" and "Database for Automated Carbohydrate-active Enzyme Annotation" and compared them to previously published oomycete genomes. Growth of Pythium spp. was assessed in a minimal medium containing selected carbon sources that are usually present in plants. The in silico analyses, coupled with our in vitro growth assays, suggest that most of the predicted CAZymes are involved in the metabolism of the oomycete cell wall with starch and sucrose serving as the main carbohydrate sources for growth of these plant pathogens. The genomes of Pythium spp. also encode pectinases and cellulases that facilitate degradation of the plant cell wall and are important in hyphal penetration; however, the species examined in this study lack the requisite genes for the complete saccharification of these carbohydrates for use as a carbon source. Genes encoding for xylan, xyloglucan, (galacto)(gluco)mannan and cutin degradation were absent or infrequent in Pythium spp.. Comparative analyses of predicted CAZymes in oomycetes indicated distinct evolutionary histories. Furthermore, CAZyme gene families among Pythium spp. were not uniformly distributed in the genomes, suggesting independent gene loss events, reflective of the polyphyletic relationships among some of the species.

Alternative splicing of a multi-drug transporter from Pseudoperonospora cubensis generates an RXLR effector protein that elicits a rapid cell death.

Pseudoperonospora cubensis, an obligate oomycete pathogen, is the causal agent of cucurbit downy mildew, a foliar disease of global economic importance. Similar to other oomycete plant pathogens, Ps. cubensis has a suite of RXLR and RXLR-like effector proteins, which likely function as virulence or avirulence determinants during the course of host infection. Using in silico analyses, we identified 271 candidate effector proteins within the Ps. cubensis genome with variable RXLR motifs. In extending this analysis, we present the functional characterization of one Ps. cubensis effector protein, RXLR protein 1 (PscRXLR1), and its closest Phytophthora infestans ortholog, PITG_17484, a member of the Drug/Metabolite Transporter (DMT) superfamily. To assess if such effector-non-effector pairs are common among oomycete plant pathogens, we examined the relationship(s) among putative ortholog pairs in Ps. cubensis and P. infestans. Of 271 predicted Ps. cubensis effector proteins, only 109 (41%) had a putative ortholog in P. infestans and evolutionary rate analysis of these orthologs shows that they are evolving significantly faster than most other genes. We found that PscRXLR1 was up-regulated during the early stages of infection of plants, and, moreover, that heterologous expression of PscRXLR1 in Nicotiana benthamiana elicits a rapid necrosis. More interestingly, we also demonstrate that PscRXLR1 arises as a product of alternative splicing, making this the first example of an alternative splicing event in plant pathogenic oomycetes transforming a non-effector gene to a functional effector protein. Taken together, these data suggest a role for PscRXLR1 in pathogenicity, and, in total, our data provide a basis for comparative analysis of candidate effector proteins and their non-effector orthologs as a means of understanding function and evolutionary history of pathogen effectors.

Expression profiling of Cucumis sativus in response to infection by Pseudoperonospora cubensis.

The oomycete pathogen, Pseudoperonospora cubensis, is the causal agent of downy mildew on cucurbits, and at present, no effective resistance to this pathogen is available in cultivated cucumber (Cucumis sativus). To better understand the host response to a virulent pathogen, we performed expression profiling throughout a time course of a compatible interaction using whole transcriptome sequencing. As described herein, we were able to detect the expression of 15,286 cucumber genes, of which 14,476 were expressed throughout the infection process from 1 day post-inoculation (dpi) to 8 dpi. A large number of genes, 1,612 to 3,286, were differentially expressed in pair-wise comparisons between time points. We observed the rapid induction of key defense related genes, including catalases, chitinases, lipoxygenases, peroxidases, and protease inhibitors within 1 dpi, suggesting detection of the pathogen by the host. Co-expression network analyses revealed transcriptional networks with distinct patterns of expression including down-regulation at 2 dpi of known defense response genes suggesting coordinated suppression of host responses by the pathogen. Comparative analyses of cucumber gene expression patterns with that of orthologous Arabidopsis thaliana genes following challenge with Hyaloperonospora arabidopsidis revealed correlated expression patterns of single copy orthologs suggesting that these two dicot hosts have similar transcriptional responses to related pathogens. In total, the work described herein presents an in-depth analysis of the interplay between host susceptibility and pathogen virulence in an agriculturally important pathosystem.

mRNA-Seq analysis of the Pseudoperonospora cubensis transcriptome during cucumber (Cucumis sativus L.) infection.

Pseudoperonospora cubensis, an oomycete, is the causal agent of cucurbit downy mildew, and is responsible for significant losses on cucurbit crops worldwide. While other oomycete plant pathogens have been extensively studied at the molecular level, Ps. cubensis and the molecular basis of its interaction with cucurbit hosts has not been well examined. Here, we present the first large-scale global gene expression analysis of Ps. cubensis infection of a susceptible Cucumis sativus cultivar, 'Vlaspik', and identification of genes with putative roles in infection, growth, and pathogenicity. Using high throughput whole transcriptome sequencing, we captured differential expression of 2383 Ps. cubensis genes in sporangia and at 1, 2, 3, 4, 6, and 8 days post-inoculation (dpi). Additionally, comparison of Ps. cubensis expression profiles with expression profiles from an infection time course of the oomycete pathogen Phytophthora infestans on Solanum tuberosum revealed similarities in expression patterns of 1,576-6,806 orthologous genes suggesting a substantial degree of overlap in molecular events in virulence between the biotrophic Ps. cubensis and the hemi-biotrophic P. infestans. Co-expression analyses identified distinct modules of Ps. cubensis genes that were representative of early, intermediate, and late infection stages. Collectively, these expression data have advanced our understanding of key molecular and genetic events in the virulence of Ps. cubensis and thus, provides a foundation for identifying mechanism(s) by which to engineer or effect resistance in the host.

The Comprehensive Phytopathogen Genomics Resource: a web-based resource for data-mining plant pathogen genomes.

The Comprehensive Phytopathogen Genomics Resource (CPGR) provides a web-based portal for plant pathologists and diagnosticians to view the genome and trancriptome sequence status of 806 bacterial, fungal, oomycete, nematode, viral and viroid plant pathogens. Tools are available to search and analyze annotated genome sequences of 74 bacterial, fungal and oomycete pathogens. Oomycete and fungal genomes are obtained directly from GenBank, whereas bacterial genome sequences are downloaded from the A Systematic Annotation Package (ASAP) database that provides curation of genomes using comparative approaches. Curated lists of bacterial genes relevant to pathogenicity and avirulence are also provided. The Plant Pathogen Transcript Assemblies Database provides annotated assemblies of the transcribed regions of 82 eukaryotic genomes from publicly available single pass Expressed Sequence Tags. Data-mining tools are provided along with tools to create candidate diagnostic markers, an emerging use for genomic sequence data in plant pathology. The Plant Pathogen Ribosomal DNA (rDNA) database is a resource for pathogens that lack genome or transcriptome data sets and contains 131 755 rDNA sequences from GenBank for 17 613 species identified as plant pathogens and related genera. Database URL: http://cpgr.plantbiology.msu.edu.