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

Genome sequencing and comparative genomics reveal insights into pathogenicity and evolution of Fusarium zanthoxyli, the causal agent of stem canker in prickly ash

BACKGROUND

Fusarium zanthoxyli is a destructive pathogen causing stem canker in prickly ash, an ecologically and economically important forest tree. However, the genome lack of F. zanthoxyli has hindered research on its interaction with prickly ash and the development of precise control strategies for stem canker.

RESULTS

In this study, we sequenced and annotated a relatively high-quality genome of F. zanthoxyli with a size of 43.39 Mb, encoding 11,316 putative genes. Pathogenicity-related factors are predicted, comprising 495 CAZymes, 217 effectors, 156 CYP450s, and 202 enzymes associated with secondary metabolism. Besides, a comparative genomics analysis revealed Fusarium and Colletotrichum diverged from a shared ancestor approximately 141.1 ~ 88.4 million years ago (MYA). Additionally, a phylogenomic investigation of 12 different phytopathogens within Fusarium indicated that F. zanthoxyli originated approximately 34.6 ~ 26.9 MYA, and events of gene expansion and contraction within them were also unveiled. Finally, utilizing conserved domain prediction, the results revealed that among the 59 unique genes, the most enriched domains were PnbA and ULP1. Among the 783 expanded genes, the most enriched domains were PKc_like kinases and those belonging to the APH_ChoK_Like family.

CONCLUSION

This study sheds light on the genetic basis of F. zanthoxyli's pathogenicity and evolution which provides valuable information for future research on its molecular interactions with prickly ash and the development of effective strategies to combat stem canker.

Resolving the Globisporangium ultimum (Pythium ultimum) species complex

The Globisporangium ultimum (formerly Pythium ultimum) species complex was previously composed of two morphological varieties: var. ultimum and var. sporangiiferum. Prior attempts to resolve this morphology-based species complex using molecular techniques have been inconclusive or conflicting. The increased availability of sequenced genomes and isolates identified as G. ultimum var. ultimum and var. sporangiiferum has allowed us to examine these relationships at a higher resolution and with a broader scope than previously possible. Using comparative genomics, we identified highly variable gene regions and designed primers for four new protein-coding genes for phylogenetics. These were then used alongside three known markers to generate a nuclear multigene genealogy of the species complex. From a collection of 163 isolates belonging to the target taxa, a subset of 29 was chosen to be included in this study (verified with nuclear rDNA internal transcribed spacer 1 [ITS1] and mitochondrial cytochrome c oxidase subunit 1 [cox1] sequences). Seventeen isolates of var. ultimum were selected to be representative of variations in genotype, morphology, and geographic collection location. The 12 isolates of var. sporangiiferum included all available specimens identified either morphologically (in previous studies) or through sequence similarity with ITS1 and cox1. Based on the fulfillment of reciprocal monophyly and observed genealogical concordance under the genealogical concordance phylogenetic species recognition, we determined that the Globisporangium ultimum species complex is composed of four genetically distinct species: Globisporangium ultimum, Globisporangium sporangiiferum, Globisporangium solveigiae, and Globisporangium bothae.

Advancing the automation of plant nucleic acid extraction for rapid diagnosis of plant diseases in space

Human space exploration missions will continue the development of sustainable plant cultivation in what are obviously novel habitat settings. Effective pathology mitigation strategies are needed to cope with plant disease outbreaks in any space-based plant growth system. However, few technologies currently exist for space-based diagnosis of plant pathogens. Therefore, we developed a method of extracting plant nucleic acid that will facilitate the rapid diagnosis of plant diseases for future spaceflight applications. The microHomogenizer^™ from Claremont BioSolutions, originally designed for bacterial and animal tissue samples, was evaluated for plant-microbial nucleic acid extractions. The microHomogenizer^™ is an appealing device in that it provides automation and containment capabilities that would be required in spaceflight applications. Three different plant pathosystems were used to assess the versatility of the extraction process. Tomato, lettuce, and pepper plants were respectively inoculated with a fungal plant pathogen, an oomycete pathogen, and a plant viral pathogen. The microHomogenizer^™, along with the developed protocols, proved to be an effective mechanism for producing DNA from all three pathosystems, in that PCR and sequencing of the resulting samples demonstrated clear DNA-based diagnoses. Thus, this investigation advances the efforts to automate nucleic acid extraction for future plant disease diagnosis in space.

Ecological insights into soil health according to the genomic traits and environment-wide associations of bacteria in agricultural soils

Soil microbiomes are sensitive to current and previous soil conditions, and bacterial 'bioindicators' of biological, physical, and chemical soil properties have considerable potential for soil health assessment. However, the lack of ecological or physiological information for most soil microorganisms limits our ability to interpret the associations of bioindicators and, thus, their utility for guiding management. We identified bioindicators of tillage intensity and twelve soil properties used to rate soil health using a 16S rRNA gene-based survey of farmland across North America. We then inferred the genomic traits of bioindicators and evaluated their environment-wide associations (EWAS) with respect to agricultural management practice, disturbance, and plant associations with 89 studies from agroecosystems. Most bioindicators were either positively correlated with biological properties (e.g., organic matter) or negatively correlated with physical and chemical properties. Higher soil health ratings corresponded with smaller genome size and higher coding density, while lower ratings corresponded with larger genomes and higher rrn copy number. Community-weighted genome size explained most variation in health ratings. EWAS linked prominent bioindicators with the impacts of environmental disturbances. Our findings provide ecological insights into bioindicators of soil properties relevant to soil health management, illustrating the tight coupling of microbiome and soil function.

Oncogenic Functions and Clinical Significance of Circular RNAs in Colorectal Cancer

Colorectal cancer (CRC) is ranked as the second most commonly diagnosed disease in females and the third in males worldwide. Therefore, the finding of new more reliable biomarkers for early diagnosis, for prediction of metastasis, and resistance to conventional therapies is an important challenge in overcoming the disease. The current review presents circular RNAs (circRNAs) with their unique features as potential prognostic and diagnostic biomarkers in CRC. The review highlights the mechanism of action and the role of circRNAs with oncogenic functions in the CRC as well as the association between their expression and clinicopathological characteristics of CRC patients. The comprehension of the role of oncogenic circRNAs in CRC pathogenesis is growing rapidly and the next step is using them as suitable new drug targets in the personalized treatment of CRC patients.

Comparative Analysis of Oomycete Genome Evolution Using the Oomycete Gene Order Browser (OGOB)

The oomycetes are a class of microscopic, filamentous eukaryotes within the stramenopiles–alveolates–rhizaria eukaryotic supergroup. They include some of the most destructive pathogens of animals and plants, such as Phytophthora infestans, the causative agent of late potato blight. Despite the threat they pose to worldwide food security and natural ecosystems, there is a lack of tools and databases available to study oomycete genetics and evolution. To this end, we have developed the Oomycete Gene Order Browser (OGOB), a curated database that facilitates comparative genomic and syntenic analyses of oomycete species. OGOB incorporates genomic data for 20 oomycete species including functional annotations and a number of bioinformatics tools. OGOB hosts a robust set of orthologous oomycete genes for evolutionary analyses. Here, we present the structure and function of OGOB as well as a number of comparative genomic analyses we have performed to better understand oomycete genome evolution. We analyze the extent of oomycete gene duplication and identify tandem gene duplication as a driving force of the expansion of secreted oomycete genes. We identify core genes that are present and microsyntenically conserved (termed syntenologs) in oomycete lineages and identify the degree of microsynteny between each pair of the 20 species housed in OGOB. Consistent with previous comparative synteny analyses between a small number of oomycete species, our results reveal an extensive degree of microsyntenic conservation amongst genes with housekeeping functions within the oomycetes. OGOB is available at https://ogob.ie.

High-throughput identification and diagnostics of pathogens and pests: Overview and practical recommendations

High-throughput identification technologies provide efficient tools for understanding the ecology and functioning of microorganisms. Yet, these methods have been only rarely used for monitoring and testing ecological hypotheses in plant pathogens and pests in spite of their immense importance in agriculture, forestry and plant community dynamics. The main objectives of this manuscript are the following: (a) to provide a comprehensive overview about the state-of-the-art high-throughput quantification and molecular identification methods used to address population dynamics, community ecology and host associations of microorganisms, with a specific focus on antagonists such as pathogens, viruses and pests; (b) to compile available information and provide recommendations about specific protocols and workable primers for bacteria, fungi, oomycetes and insect pests; and (c) to provide examples of novel methods used in other microbiological disciplines that are of great potential use for testing specific biological hypotheses related to pathology. Finally, we evaluate the overall perspectives of the state-of-the-art and still evolving methods for diagnostics and population- and community-level ecological research of pathogens and pests.

MM, Prakash Mishra G, Singh D, Ye W, Tyler BM, Tripathy S. EumicrobeDBLite: a lightweight genomic resource and analytic platform for draft oomycete genomes

We have developed EumicrobeDBLite-a lightweight comprehensive genome resource and sequence analysis platform for oomycete organisms. EumicrobeDBLite is a successor of the VBI Microbial Database (VMD) that was built using the Genome Unified Schema (GUS). In this version, GUS has been greatly simplified with the removal of many obsolete modules and the redesign of others to incorporate contemporary data. Several dependences, such as perl object layers used for data loading in VMD, have been replaced with independent lightweight scripts. EumicrobeDBLite now runs on a powerful annotation engine developed at our laboratory, called 'Genome Annotator Lite'. Currently, this database has 26 publicly available genomes and 10 expressed sequence tag (EST) datasets of oomycete organisms. The browser page has dynamic tracks presenting comparative genomics analyses, coding and non-coding data, tRNA genes, repeats and EST alignments. In addition, we have defined 44 777 core conserved proteins from 12 oomycete organisms which form 2974 clusters. Synteny viewing is enabled by the incorporation of the Genome Synteny Viewer (GSV) tool. The user interface has undergone major changes for ease of browsing. Queryable comparative genomics information, conserved orthologous genes and pathways are among the new key features updated in this database. The browser has been upgraded to enable user upload of GFF files for quick view of genome annotation comparisons. The toolkit page integrates the EMBOSS package and has a gene prediction tool. Annotations for the organisms are updated once every 6 months to ensure quality. The database resource is available at www.eumicrobedb.org.

Detecting alternatively spliced transcript isoforms from single-molecule long-read sequences without a reference genome

Alternative splicing (AS) is a major source of transcript and proteome diversity, but examining AS in species without well-annotated reference genomes remains difficult. Research on both human and mouse has demonstrated the advantages of using Iso-Seq™ data for isoform-level transcriptome analysis, including the study of AS and gene fusion. We applied Iso-Seq™ to investigate AS in Amborella trichopoda, a phylogenetically pivotal species that is sister to all other living angiosperms. Our data show that, compared with RNA-Seq data, the Iso-Seq™ platform provides better recovery on large transcripts, new gene locus identification and gene model correction. Reference-based AS detection with Iso-Seq™ data identifies AS within a higher fraction of multi-exonic genes than observed for published RNA-Seq analysis (45.8% vs. 37.5%). These data demonstrate that the Iso-Seq™ approach is useful for detecting AS events. Using the Iso-Seq-defined transcript collection in Amborella as a reference, we further describe a pipeline for detection of AS isoforms from PacBio Iso-Seq™ without using a reference sequence (de novo). Results using this pipeline show a 66%-76% overall success rate in identifying AS events. This de novoAS detection pipeline provides a method to accurately characterize and identify bona fide alternatively spliced transcripts in any nonmodel system that lacks a reference genome sequence. Hence, our pipeline has huge potential applications and benefits to the broader biology community.

Comparative transcriptomic analysis indicates genes associated with local and systemic resistance to Colletotrichum graminicola in maize

The hemibiotrophic fungus Colletotrichum graminicola may cause severe damage to maize, affecting normal development of the plant and decreasing grain yield. In this context, understanding plant defense pathways at the inoculation site and systemically in uninoculated tissues can help in the development of genetic engineering of resistance against this pathogen. Previous work has discussed the molecular basis of maize - C. graminicola interaction. However, many genes involved in defense have not yet been exploited for lack of annotation in public databases. Here, changes in global gene expression were studied in root, male and female inflorescences of maize under local and systemic fungal infection treatments, respectively. RNA-Seq with qPCR was used to indicate genes involved in plant defense. We found that systemic acquired resistance induction in female inflorescences mainly involves accumulation of salicylic acid (SA)-inducible defense genes (ZmNAC, ZmHSF, ZmWRKY, ZmbZIP and PR1) and potential genes involved in chromatin modification. Furthermore, transcripts involved in jasmonic acid (JA) and ethylene (ET) signaling pathways were also accumulated and may participate in plant immunity. Moreover, several genes were functionally re-annotated based on domain signature, indicating novel candidates to be tested in strategies involving gene knockout and overexpression in plants.

Comparative analysis of Cucurbita pepo metabolism throughout fruit development in acorn squash and oilseed pumpkin

Both the fruit mesocarp and the seeds of winter squash can be used for consumption, although the focus of breeding efforts varies by cultivar. Cultivars bred for fruit consumption are selected for fruit mesocarp quality traits such as carotenoid content, percent dry matter, and percent soluble solids, while these traits are essentially ignored in oilseed pumpkins. To compare fruit development in these two types of squash, we sequenced the fruit transcriptome of two cultivars bred for different purposes: an acorn squash, 'Sweet REBA', and an oilseed pumpkin, 'Lady Godiva'. Putative metabolic pathways were developed for carotenoid, starch, and sucrose synthesis in winter squash fruit and squash homologs were identified for each of the structural genes in the pathways. Gene expression, especially of known rate-limiting and branch point genes, corresponded with metabolite accumulation both across development and between the two cultivars. Thus, developmental regulation of metabolite genes is an important factor in winter squash fruit quality.

Prediction of bacterial associations with plants using a supervised machine-learning approach

Recent scenarios of fresh produce contamination by human enteric pathogens have resulted in severe food-borne outbreaks, and a new paradigm has emerged stating that some human-associated bacteria can use plants as secondary hosts. As a consequence, there has been growing concern in the scientific community about these interactions that have not yet been elucidated. Since this is a relatively new area, there is a lack of strategies to address the problem of food-borne illnesses due to the ingestion of fruits and vegetables. In the present study, we performed specific genome annotations to train a supervised machine-learning model that allows for the identification of plant-associated bacteria with a precision of ∼93%. The application of our method to approximately 9500 genomes predicted several unknown interactions between well-known human pathogens and plants, and it also confirmed several cases for which evidence has been reported. We observed that factors involved in adhesion, the deconstruction of the plant cell wall and detoxifying activities were highlighted as the most predictive features. The application of our strategy to sequenced strains that are involved in food poisoning can be used as a primary screening tool to determine the possible causes of contaminations.

Comparative genomics and phylogenetic discordance of cultivated tomato and close wild relatives

Background. Studies of ancestry are difficult in the tomato because it crosses with many wild relatives and species in the tomato clade that have diverged very recently. As a result, the phylogeny in relation to its closest relatives remains uncertain. By using the coding sequence from Solanum lycopersicum, S. galapagense, S. pimpinellifolium, S. corneliomuelleri, and S. tuberosum and the genomic sequence from S. lycopersicum ‘Heinz’, an heirloom line, S. lycopersicum ‘Yellow Pear’, and two of cultivated tomato’s closest relatives, S. galapagense and S. pimpinellifolium, we have aimed to resolve the phylogenies of these closely related species as well as identify phylogenetic discordance in the reference cultivated tomato.

Results. Divergence date estimates suggest that the divergence of S. lycopersicum, S. galapagense, and S. pimpinellifolium happened less than 0.5 MYA. Phylogenies based on 8,857 coding sequences support grouping of S. lycopersicum and S. galapagense, although two secondary trees are also highly represented. A total of 25 genes in our analysis had sites with evidence of positive selection along the S. lycopersicum lineage. Whole genome phylogenies showed that while incongruence is prevalent in genomic comparisons between these genotypes, likely as a result of introgression and incomplete lineage sorting, a primary phylogenetic history was strongly supported.

Conclusions. Based on analysis of these genotypes, S. galapagense appears to be closely related to S. lycopersicum, suggesting they had a common ancestor prior to the arrival of an S. galapagense ancestor to the Galápagos Islands, but after divergence of the sequenced S. pimpinellifolium. Genes showing selection along the S. lycopersicum lineage may be important in domestication or selection occurring post-domestication. Further analysis of intraspecific data in these species will help to establish the evolutionary history of cultivated tomato. The use of an heirloom line is helpful in deducing true phylogenetic information of S. lycopersicum and identifying regions of introgression from wild species.

circBase: a database for circular RNAs

Recently, several laboratories have reported thousands of circular RNAs (circRNAs) in animals. Numerous circRNAs are highly stable and have specific spatiotemporal expression patterns. Even though a function for circRNAs is unknown, these features make circRNAs an interesting class of RNAs as possible biomarkers and for further research. We developed a database and website, "circBase," where merged and unified data sets of circRNAs and the evidence supporting their expression can be accessed, downloaded, and browsed within the genomic context. circBase also provides scripts to identify known and novel circRNAs in sequencing data. The database is freely accessible through the web server at http://www.circbase.org/.

Current impact and future directions of high throughput sequencing in plant virus diagnostics

The ability to provide a fast, inexpensive and reliable diagnostic for any given viral infection is a key parameter in efforts to fight and control these ubiquitous pathogens. The recent developments of high-throughput sequencing (also called Next Generation Sequencing - NGS) technologies and bioinformatics have drastically changed the research on viral pathogens. It is now raising a growing interest for virus diagnostics. This review provides a snapshot vision on the current use and impact of high throughput sequencing approaches in plant virus characterization. More specifically, this review highlights the potential of these new technologies and their interplay with current protocols in the future of molecular diagnostic of plant viruses. The current limitations that will need to be addressed for a wider adoption of high-throughput sequencing in plant virus diagnostics are thoroughly discussed.

Horizontal gene transfer in eukaryotic plant pathogens

Gene transfer has been identified as a prevalent and pervasive phenomenon and an important source of genomic innovation in bacteria. The role of gene transfer in microbial eukaryotes seems to be of a reduced magnitude but in some cases can drive important evolutionary innovations, such as new functions that underpin the colonization of different niches. The aim of this review is to summarize published cases that support the hypothesis that horizontal gene transfer (HGT) has played a role in the evolution of phytopathogenic traits in fungi and oomycetes. Our survey of the literature identifies 46 proposed cases of transfer of genes that have a putative or experimentally demonstrable phytopathogenic function. When considering the life-cycle steps through which a pathogen must progress, the majority of the HGTs identified are associated with invading, degrading, and manipulating the host. Taken together, these data suggest HGT has played a role in shaping how fungi and oomycetes colonize plant hosts.

Identification of novel target genes for safer and more specific control of root-knot nematodes from a pan-genome mining

Root-knot nematodes are globally the most aggressive and damaging plant-parasitic nematodes. Chemical nematicides have so far constituted the most efficient control measures against these agricultural pests. Because of their toxicity for the environment and danger for human health, these nematicides have now been banned from use. Consequently, new and more specific control means, safe for the environment and human health, are urgently needed to avoid worldwide proliferation of these devastating plant-parasites. Mining the genomes of root-knot nematodes through an evolutionary and comparative genomics approach, we identified and analyzed 15,952 nematode genes conserved in genomes of plant-damaging species but absent from non target genomes of chordates, plants, annelids, insect pollinators and mollusks. Functional annotation of the corresponding proteins revealed a relative abundance of putative transcription factors in this parasite-specific set compared to whole proteomes of root-knot nematodes. This may point to important and specific regulators of genes involved in parasitism. Because these nematodes are known to secrete effector proteins in planta, essential for parasitism, we searched and identified 993 such effector-like proteins absent from non-target species. Aiming at identifying novel targets for the development of future control methods, we biologically tested the effect of inactivation of the corresponding genes through RNA interference. A total of 15 novel effector-like proteins and one putative transcription factor compatible with the design of siRNAs were present as non-redundant genes and had transcriptional support in the model root-knot nematode Meloidogyne incognita. Infestation assays with siRNA-treated M. incognita on tomato plants showed significant and reproducible reduction of the infestation for 12 of the 16 tested genes compared to control nematodes. These 12 novel genes, showing efficient reduction of parasitism when silenced, constitute promising targets for the development of more specific and safer control means.

Plant-parasitic nematodes are annually responsible for more than $100 billion crop yield loss worldwide and those considered as causing most of the damages are root-knot nematodes. These nematodes used to be controlled by chemicals that are now banned from use because of their poor specificity and high toxicity for the environment and human health. In the absence of sustainable alternative solutions, new control means, more specifically targeted against these nematodes and safe for the environment are needed. We searched in root-knot nematode genomes, genes conserved in various plant-damaging species while otherwise absent from the genomes of non target species such as those of chordates, plants, annelids, insect pollinators and mollusks. These genes are probably important for plant parasitism and their absence from non-target species make them interesting candidates for the development of more specific and safer control means. Further bioinformatics pruning of this set of genes yielded 16 novel candidates that could be biologically tested. Using RNA interference, we knocked down each of these 16 genes in a root-knot nematode and tested the effect on plant parasitism efficiency. Out of the 16 tested genes, 12 showed a significant and reproducible diminution of infestation when silenced and are thus particularly promising.

The role of prophage in plant-pathogenic bacteria

A diverse set of phage lineages is associated with the bacterial plant-pathogen genomes sequenced to date. Analysis of 37 genomes revealed 5,169 potential genes (approximately 4.3 Mbp) of phage origin, and at least 50% had no function assigned or are nonessential to phage biology. Some phytopathogens have transcriptionally active prophage genes under conditions that mimic plant infection, suggesting an association between plant disease and prophage transcriptional modulation. The role of prophages within genomes for cell biology varies. For pathogens such as Pectobacterium, Pseudomonas, Ralstonia, and Streptomyces, involvement of prophage in disease symptoms has been demonstrated. In Xylella and Xanthomonas, prophage activity is associated with genome rearrangements and strain differentiation. For other pathogens, prophage roles are yet to be established. This review integrates available information in a unique interface ( http://propnav.esalq.usp.br ) that may be assessed to improve research in prophage biology and its association with genome evolution and pathogenicity.

Genome-enabled perspectives on the composition, evolution, and expression of virulence determinants in bacterial plant pathogens

Genome sequence analyses of bacterial plant pathogens are revealing important insights into the molecular determinants of pathogenicity and, through transcript characterization, responses to environmental conditions, evidence for small RNAs, and validation of uncharacterized genes. Genome comparison sheds further light on the processes impacting pathogen evolution and differences in gene repertoire among isolates contributing to niche specialization. Information derived from pathogen genome analysis is providing tools for use in diagnosis and interference with host-pathogen interactions for the purpose of disease control. However, the existing information infrastructure fails to adequately integrate the increasing numbers of sequence data sets, bioinformatic analyses, and experimental characterization, as required for effective systems-level analysis. Enhanced standardization of data formats at the point of publication is proposed as a possible solution.