Keywords
Phytophthora genome, plant pathogen, Phytophthora cinnamomi
This article is included in the Agriculture, Food and Nutrition gateway.
Phytophthora genome, plant pathogen, Phytophthora cinnamomi
Phytophthora cinnamomi is a highly virulent plant pathogen that has a devastating impact of the Australian ecosystem, namely in the south-western areas of Western Australia and much of the south and east coasts of Victoria and New South Wales1. In the south west Botanical Province of Western Australia, alone, over 40% of the 5710 plant species present have been shown to be susceptible to P. cinnamomi2. Significant genetic and phenotypic variation can occur within a signal clonal linage of P. cinnamomi3 and susceptibility of a given host plant species has been shown to vary from site to site4. Furthermore, despite the general lack of crossing during sexual reproduction, P. cinnamomi excels at adapting to new environments and developing virulence to new host species through asexual growth, making it a deadly and difficult-to-control pathogen. Unravelling how P. cinnamomi is able to adapt so quickly, and remain virulent to a wide range of hosts in Australia, is an important research goal.
The currently available genome of P. cinnamomi var. cinnamomi (Joint Genome Institute (JGI); NCBI Accession no. PRJNA68241) is based on the Rans isolate from Sumatra in 1922, which has been in culture for many decades and may not be representative of the current pathogenic strains present in Australia. Here we report and make available two Australian P. cinnamomi genomes, isolated from geographically very separate areas with different available host species. After analyses of genetic differences between these two P. cinnamomi genomes, it may be that key genes or gene families under high evolutionary pressure can be identified; this may aid further studies on more effective control of this pathogen.
Two isolates of P. cinnamomi were selected from areas of infection on either side of the Australian continent: one from the Brisbane Ranges in southeastern Australia (DU054, A2 mating type)5 and the other from southwestern Western Australia (WA94.26, A2 mating type), both Deakin University culture collection. These isolates were selected to represent possible genetic diversity of P. cinnamomi in Australia arising from geographic isolation, and possible variation of selective pressures due to different host species. Isolates were maintained on V8 agar (V8A) [50 ml unclarified V8 ‘Original’ Juice (Campbells, Australia), 0.5 g CaCO3 and 7.5 g biological agar per 500 mL of distilled water] at 25°C in darkness, as per 5. Genomic DNA was isolated from hyphae using a DNeasy Plant Mini Kit (Qiagen), following the manufacturers protocol. Illumina TruSeq Nano library preparation and sequencing on an Illumina HiSeq 2500 platform were performed by the Australian Genome Resources Facility (Walter and Eliza Hall Institute, Parkville, Australia) generating ~154 million (2x 150bp) raw reads per isolate. Raw reads are available in the NCBI Short Read Archive (SRA) under the Bioproject Accession: PRJNA413098
Raw sequencing data for each isolate was first pre-processed using Trimmomatic v0.336 with the following parameters: ILLUMINACLIP:TruSeq3-PE-2.fa:2:30:10:4 AVGQUAL:30 MINLEN:36, to remove Illumina adapters and filter reads based on quality scores (Phred score). Only reads with average Phred > 30 were retained. To ensure only the desired P. cinnamomi genomes were assembled, a second round of pre-processing was conducted to remove potential contaminants. MetaPhlAn v27, was run with default settings and identified the Paenibacillus genus as a likely contaminate. Using BBMap v0.35 (BBMap - Bushnell. B), we mapped the Trimmomatic-filtered reads to the closest species match (Paenibacillus sp., JDR-2, GenBank accession: GCA_000023585.1, with 2.7% and 2.0% of DU054 and WA94.26 reads mapping, respectively; these Paenibacillus reads were subsequently removed. The remaining reads were then mapped using BBMap to the human genome (GRCh38; NCBI accession: GCA_000001405.15), with < 0.5% (~ 430,000 reads from DU054 and ~ 630,000 from WA94.26) being mapped and subsequently removed from the data set. Thus, the final set of reads (DU054, 149 million reads; WA94.26, 151 million reads) used for the assembly contained highly quality paired-end reads not belonging to either human or bacterial contaminants.
De novo contig assembly of the two genomes was conducted independently, using IDBA-UD v1.1.08. IDBA-UD was run using the following parameters: --mink 20 --maxk 100 --step 20 --min_contig 500 --min_support 2 --min_count 3. Briefly, these conducted a multiple K-mer assembly from k = 20 up to k = 100; only assembled contigs above 500 bp and those with a minimum depth coverage ≥ 3 were kept. As heterogeneous data can increase redundancy in genome assemblies, the IDBA-UD assembled contigs were run through the Redundans pipeline v0.12c9 with the following parameters: -threads 4 -min_length 500. Redundans uses paired-end mapping data to reduce assembled sequence redundancy and scaffold contigs into longer less fragmented sequences. The final assembled genome of DU054 was 62.80 Mb in 14,269 scaffolds with an N50 of 9,951; the longest scaffold was 1.54 Mb in length (Table 1). For WA94.26, the final genome was 68.07 Mb in length, in 10,085 scaffolds with the largest being 1.54 Mb and an N50 of 20,813. GC content remained consistent, at ~ 53%, between both isolate genomes across both assemblies and before and after processing with Redundans. The quality, as measured by the above metrics, of the presented genomes is comparable to the previously available P. cinnamomi var. cinnamomi Rans isolate genome (JGI).. The final genome assemblies are available under the NCBI Bioproject Accession: PRJNA413098.
DU054 | WA94.26 | |||
---|---|---|---|---|
IDBA-UD | Redundans | IDBA-UD | Redundans | |
Assembly size (Mb) | 71.29 | 62.80 | 76.95 | 68.07 |
No. scaffolds | 33,475 | 14,268 | 36,333 | 10,084 |
N50 (bp) | 4,085 | 9,951 | 4,075 | 20,813 |
No. predicted genes | NA | 23,414 | NA | 22,573 |
We used the BUSCO (benchmarking universal single-copy orthologs) pipeline v1.2210 with the default e-value cutoff of 0.01, to assess the completeness of the assembled genomes and compared the results to the previously available Rans isolate. Utilizing the set of 429 conserved eukaryotic single-copy orthologs (hereafter BUSCOs), the analysis indicated 82.5% and 81.8% BUSCO completeness for DU054 and WA94.26 genomes, respectively. For DU054, 335 complete BUSCOs (including single-copy and duplicated BUSCOs) and 19 fragmented BUSCOs were identified, and 333 complete and 19 fragmented BUSCOs in WA94.26 (Table 2). Overall, we find comparable levels of BUSCO completeness with the Rans isolate, suggesting our two Australian isolate assemblies are as complete references as that currently available.
In the absence of any available high quality ESTs (expressed sequence tags) or transcriptome (gene expression) data for P. cinnamomi, we conducted a preliminary protein-coding sequence prediction using GeneMark-ES v4.3211, which utilises a self-training algorithm to identify exon, intron and intergenic regions as well as initiation and termination sites. GeneMark-ES was run using the default settings and a database of predicted gene models (i.e., predicted polypeptides) was constructed for DU054 and WA94.26 genomes. An initial 23,414 gene models were identified in DU054 and 22,573 in WA94.26. Of these, 14,735 pairs of predicted gene models appear to be orthologous between the two genomes (reciprocal best-hit Blastp10, e value ≤ 1e-5). As a preliminary verification of these gene model builds, we identified orthologous counterparts to eight available Phytophthora genomes with more complete annotations (P. infestans12, P. kernoviae13, P. lateralis14, P. nicotianae15, P. parasitica (P1569_v1; Broad Institute), P. ramorum16, P. sojae16 and P. cinnamomi var. cinnamomi). Accordingly, we used OrthoFinder v1.1.1017 with default parameters, except we used DIAMOND18 as the alignment program with the diamond_more_sensitive flag. OrthoFinder first identifies ‘orthogroups’ (an extension of orthologues to include groups of genes descended from a single gene in the last common ancestor of a group of species17) and then orthologues between each pair of species in the comparison. OrthoFinder assigned 88.5% (170,769) of the genes found in all the species to 19,089 orthogroups, and of these 50% of all the genes were contained in orthogroups, which had 10 or more genes within them. We found 2,931 orthogroups that contained genes for each of the species, and of these 1,309 orthogroups consisted entirely of single copy genes, see associated data repository19. Using these single copy orthogroups gene trees were first constructed then the species tree was inferred using the distance-based implemented by fastme20. The resultant species tree (see associated data repository19) exhibits strong congruence to the Phytophthora phylogeny recently published by 21, providing more evidence that the genome assembly and preliminary annotation conducted here is valuable.
In summary, we present the genome assembly of two geographically separated isolates of Phytophthora cinnamomi from Australia, representing the first genome assembly of an Australian-isolated strain. These high-quality genomes will act as a valuable resource, particularly for the further identification and analysis of protein-encoding genes, which are expressed during plant infection, such as members of the avirulence gene families22. These gene families are of specific interest in the development of novel and effective pathogen control mechanisms.
Raw reads are available in the NCBI SRA under the Bioproject Accession: PRJNA413098
The final assemblies are available at DDBJ/EMBL/GenBank under the accessions, PDCY00000000 and PDCZ00000000 and under the Bioproject Accession: PRJNA413098.
Supporting data, including OrthoFinder analysis and BUSCO assessment results can be found in the associated data repository: doi, 10.4225/16/59d15a6917a5e19. Data are available under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0).
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the work clearly and accurately presented and does it cite the current literature?
Partly
Is the study design appropriate and is the work technically sound?
Partly
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
References
1. Studholme DJ, McDougal RL, Sambles C, Hansen E, et al.: Genome sequences of six Phytophthora species associated with forests in New Zealand.Genom Data. 2016; 7: 54-6 PubMed Abstract | Publisher Full TextCompeting Interests: No competing interests were disclosed.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | |||
---|---|---|---|
1 | 2 | 3 | |
Version 2 (revision) 28 Feb 18 |
read | ||
Version 1 08 Nov 17 |
read | read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)