High-Quality Reference Genome Sequence for the Oomycete Vegetable Pathogen Phytophthora capsici Strain LT1534

The oomycete Phytophthora capsici is a destructive pathogen of a wide range of vegetable hosts, especially peppers and cucurbits. A 94.17-Mb genome assembly was constructed using PacBio and Illumina data and annotated with support from transcriptome sequencing (RNA-Seq) reads. P hytophthora capsici is a highly destructive pathogen of vegetables worldwide, especially those in the Solanaceae and Cucurbitaceae families (Fig. 1). Efforts to understand the role of genes in pathogenesis and host range (1–8), population structure (9, 10), genomic variation (11), and development (12) benefit from an accurate genome assembly. P. capsici belongs to the eukaryotic phylum Oomycota, many members of which have repeat-rich genomes that are difficult to assemble with short reads. Strain LT1534, A2 mating type, is an inbred strain derived from crossing isolates from infected cucurbits in Michigan (cucumber) and Tennessee (pumpkin) (11). An earlier assembly for LT1534, generated with Roche 454 and Sanger sequencing (11), was 64Mb in 917 scaffolds (N50, 706 kb). A single oospore-derived isolate (LT1534) was maintained axenically on PARP (25 ppm pimaricin, 100 ppm ampicillin, 25 ppm rifampicin, and 100 ppm pentachloronitrobenzene) agar plates, and a small weft of mycelium was transferred to clarified V8 juice broth amended with PARP and grown at 22 to 23°C for 7 days prior to harvesting, freeze-drying, disruption, and extraction of genomic DNA (13) using a GeneJet genomic DNA purification kit (Thermo Fisher). Libraries were constructed with a TruSeq DNA kit and sequenced on an Illumina HiSeq X Ten system by Novogene (Shenzen, China), yielding 34.9 million, 2 150-bp read pairs (10.48 Gbp). PacBio libraries were constructed by the National Center for Genome Resources (Santa Fe, NM) following the 20-kb protocol (14). One single-molecule real-time (SMRT) cell was sequenced on a PacBio RS II system using P6 polymerase and C4 chemistry, yielding 964,374 subreads (n= 964,374; N50, 12.5 kb; maximum, 47.5 kb; total, 8.15 Gbp) as processed by the SMRT pipeline v2.3.0.1. PacBio reads were corrected with Illumina by LoRDEC v0.9 (15) using the default parameters. To provide gene model support, mRNA was extracted from 2.5-, 3-, and 4.5-day V8 agar cultures grown at 23°C under 12-h light/dark conditions using a Spectrum plant total RNA kit (Sigma). Strand-specific transcriptome sequencing (RNA-Seq) libraries were constructed and sequenced by Cofactor Genomics (St. Louis, MO) using oligo(dT) priming on an Illumina NextSeq 500 instrument to obtain;25 million single-end reads per library. Read trimming, correction, and genome assembly were performed with MaSuRCA v3.3.1 (16, 17) with Illumina and LoRDEC-corrected PacBio reads (LHE_COVERAGE=25 Citation Stajich JE, Vu AL, Judelson HS, Vogel GM, Gore MA, Carlson MO, Devitt N, Jacobi J, Mudge J, Lamour KH, Smart CD. 2021. Highquality reference genome sequence for the oomycete vegetable pathogen Phytophthora capsici strain LT1534. Microbiol Resour Announc 10:e00295-21. https://doi.org/10 .1128/MRA.00295-21. Editor Antonis Rokas, Vanderbilt University Copyright © 2021 Stajich et al. This is an openaccess article distributed under the terms of the Creative Commons Attribution 4.0 International license. Address correspondence to Jason E. Stajich, jason.stajich@ucr.edu. * Present address: Andrea L. Vu, Department of Molecular Biomedical Sciences, North Carolina State University, Raleigh, North Carolina, USA; Maryn O. Carlson, Committee on Genetics, Genomics, and Systems Biology, University of Chicago, Chicago, Illinois, USA. Received 22 March 2021 Accepted 7 May 2021 Published 27 May 2021 Volume 10 Issue 21 e00295-21 mra.asm.org 1 GENOME SEQUENCES

[1]  J. Stajich,et al.  Funannotate v1.8.1: Eukaryotic genome annotation , 2020 .

[2]  M. Zhang,et al.  A Phytophthora capsici effector suppresses plant immunity via interaction with EDS1 , 2020, Molecular plant pathology.

[3]  Steven L Salzberg,et al.  Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype , 2019, Nature Biotechnology.

[4]  Wenbo Ma,et al.  The RXLR effector PcAvh1 is required for full virulence of Phytophthora capsici. , 2019, Molecular plant-microbe interactions : MPMI.

[5]  D. Shen,et al.  A Phytophthora capsici Effector Targets ACD11 Binding Partners that Regulate ROS-Mediated Defense Response in Arabidopsis. , 2019, Molecular plant.

[6]  The UniProt Consortium,et al.  UniProt: a worldwide hub of protein knowledge , 2018, Nucleic Acids Res..

[7]  Zhenglu Yang,et al.  dbCAN2: a meta server for automated carbohydrate-active enzyme annotation , 2018, Nucleic Acids Res..

[8]  L. Quesada-Ocampo,et al.  Analysis of microsatellites from transcriptome sequences of Phytophthora capsici and applications for population studies , 2018, Scientific Reports.

[9]  Neil D. Rawlings,et al.  The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database , 2017, Nucleic Acids Res..

[10]  Robert M. Waterhouse,et al.  BUSCO Applications from Quality Assessments to Gene Prediction and Phylogenomics , 2017, bioRxiv.

[11]  Silvio C. E. Tosatto,et al.  InterPro in 2017—beyond protein family and domain annotations , 2016, Nucleic Acids Res..

[12]  Luis Pedro Coelho,et al.  Fast Genome-Wide Functional Annotation through Orthology Assignment by eggNOG-Mapper , 2016, bioRxiv.

[13]  S. Salzberg,et al.  Hybrid assembly of the large and highly repetitive genome of Aegilops tauschii, a progenitor of bread wheat, with the mega-reads algorithm , 2016, bioRxiv.

[14]  Evgeny M. Zdobnov,et al.  BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs , 2015, Bioinform..

[15]  D. Shen,et al.  A Virulence Essential CRN Effector of Phytophthora capsici Suppresses Host Defense and Induces Cell Death in Plant Nucleus , 2015, PloS one.

[16]  James K. Hane,et al.  CodingQuarry: highly accurate hidden Markov model gene prediction in fungal genomes using RNA-seq transcripts , 2015, BMC Genomics.

[17]  Chao Xie,et al.  Fast and sensitive protein alignment using DIAMOND , 2014, Nature Methods.

[18]  Leena Salmela,et al.  LoRDEC: accurate and efficient long read error correction , 2014, Bioinform..

[19]  Adam M Phillippy,et al.  Long-read, whole-genome shotgun sequence data for five model organisms , 2014, Scientific Data.

[20]  J. Vega-Arreguín,et al.  Recognition of an Avr3a homologue plays a major role in mediating nonhost resistance to Phytophthora capsici in Nicotiana species. , 2014, Molecular plant-microbe interactions : MPMI.

[21]  Michael Roberts,et al.  The MaSuRCA genome assembler , 2013, Bioinform..

[22]  E. Huitema,et al.  Characterization of cell death inducing Phytophthora capsici CRN effectors suggests diverse activities in the host nucleus , 2013, Front. Plant Sci..

[23]  Xiao-Ren Chen,et al.  RNA-Seq Reveals Infection-Related Gene Expression Changes in Phytophthora capsici , 2013, PloS one.

[24]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[25]  P. Hedley,et al.  Identification and Characterisation CRN Effectors in Phytophthora capsici Shows Modularity and Functional Diversity , 2013, PloS one.

[26]  James R. Knight,et al.  Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable pathogen Phytophthora capsici. , 2012, Molecular plant-microbe interactions : MPMI.

[27]  E. Huitema,et al.  The oomycete broad-host-range pathogen Phytophthora capsici. , 2012, Molecular plant pathology.

[28]  L. Quesada-Ocampo,et al.  Investigating the genetic structure of Phytophthora capsici populations. , 2011, Phytopathology.

[29]  N. Friedman,et al.  Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data , 2011, Nature Biotechnology.

[30]  M. Borodovsky,et al.  Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. , 2008, Genome research.

[31]  David Haussler,et al.  Using native and syntenically mapped cDNA alignments to improve de novo gene finding , 2008, Bioinform..

[32]  Jonathan E. Allen,et al.  Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments , 2007, Genome Biology.

[33]  K. Lamour,et al.  A strategy for recovering high quality genomic DNA from a large number of Phytophthora isolates. , 2006, Mycologia.

[34]  Ewan Birney,et al.  Automated generation of heuristics for biological sequence comparison , 2005, BMC Bioinformatics.

[35]  Steven Salzberg,et al.  TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders , 2004, Bioinform..

[36]  Ian Korf,et al.  Gene finding in novel genomes , 2004, BMC Bioinformatics.