<abstract xmlns="http://www.w3.org/1999/xhtml">

<p><italic>Heterodera carotae</italic>, the carrot cyst nematode, is a significant pest affecting carrot globally. Here we present the draft genome of <italic>H. carotae</italic>, which was generated from short read libraries from Illumina HiSeq technology, and the corresponding genome annotation.</p>
</abstract>

Heterodera carotae, the carrot cyst nematode, is a significant pest affecting carrot globally. Here we present the draft genome of H. carotae, which was generated from short read libraries from Illumina HiSeq technology, and the corresponding genome annotation.

<div xmlns="http://www.w3.org/1999/xhtml">
<p>The carrot cyst nematode, <italic>Heterodera carotae</italic>, is a considerable pest affecting carrot growing regions around the world, including much of Europe, South Africa, Mexico, Cyprus, Chile, Ontario, Canada, and Michigan in the United States (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_002">Berney and Bird, 1992</a>; <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_022">Subbotin <italic>et al</italic>., 2010</a>; <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_026">Yu <italic>et al</italic>., 2017</a>; <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_006">Escobar-Avila <italic>et al</italic>., 2018</a>; <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_016">Madani <italic>et al</italic>., 2018</a>; <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_003">CABI, 2021</a>). Although <italic>H. carotae</italic> has a host range limited to species in the Apiaceae family, mainly wild and cultivated carrot (<italic>Daucus carota</italic>), infestation of commercial production fields by this nematode can be devastating with yield losses between 12% and 80% (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_009">Greco <italic>et al</italic>., 1993</a>; <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_022">Subbotin <italic>et al</italic>., 2010</a>). Despite the significance of this pest, there are no genomic or transcriptomic resources publically available for this nematode. To date, research utilizing molecular tools in <italic>H. carotae</italic> has been limited to the development of a molecular diagnostic for populations of <italic>H. carotae</italic> from Ontario, Canada, and Northern Italy (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_016">Madani <italic>et al</italic>., 2018</a>) and microsatellite genotyping of <italic>H. carotae</italic> to examine the gene flow in French populations (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_007">Esquibet <italic>et al</italic>., 2020</a>). The <italic>H. carotae</italic> genome will provide a valuable resource to researchers tackling pest management of <italic>H. carotae</italic> and to the wider cyst nematode research community.</p>
<p>In 2019, 1 kg of soil was collected from a carrot field in Calama, Antofagasta region, Chile. The soil was placed in a paper bag. At the Nematology Laboratory of Agriculture and Livestock Service (Santiago, Chile) the soil samples were placed on trays and allowed to dry for at least a week. Cysts were then extracted from 250 g of soil using a modified Fenwick can (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_008">Fenwick, 1940</a>). The cysts were picked from the samples and placed in DESS (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_025">Yoder <italic>et al</italic>., 2006</a>). The samples were shipped to the USDA-ARS Horticultural Crops Research Unit in Corvallis, OR for DNA extraction and sequencing. From the picked cysts, 25 were hand selected, broken open in sterilized water using a scalpel, and the eggs were collected using a Pasteur pipette. Genomic DNA was extracted from the eggs using the QIAmp DNA Micro Kit (Qiagen, Hilden, Germany) following manufacturer's instructions. Sequencing and library preparation were done at the Center for Qualitative Life Sciences at Oregon State University (Corvallis, OR). The NEBNext Ultra II DNA Library Prep Kit for Illumina (San Diego, CA) was used to prepare a whole genome shotgun library from egg DNA and sequencing was performed on the Illumina HiSeq 3000 platform. The cysts were also sent to the USDA-ARS Mycology and Nematology Genetic Diversity and Biology Laboratory (Beltsville, MD) to confirm species identification through morphometrics.</p>
<p>To assemble the genome, raw sequencing data was first subjected to adaptor removal and quality control (<italic>Q</italic> = 20) using Trim Galore! v. 0.6.6 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_014">Krueger, 2020</a>) resulting in 33,127,853 150-bp paired-end reads. Using the filtered reads, a de novo assembly was generated using metaSPAdes v. 3.15.3 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_018">Nurk <italic>et al</italic>., 2017</a>). The de novo meta-assembly was visualized using the Blob Tools workflow (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_015">Kumar <italic>et al</italic>., 2013</a>) to identify and remove contaminating contigs. Each contig in the meta-assembly was assigned a phylum ID based on BLAST similarity (<italic>E</italic>-value &lt; 10e<sup>−25</sup>) to sequences in the NCBI “nt” database or other plant-parasitic nematode genomes (<italic>H. schatti</italic> [NCBI BioProject:PRJNA722882], <italic>H. glycines</italic> [NCBI BioProject:PRJNA381081], <italic>Meloidogyne incognita</italic> [ENA Project:PRJEB8714], <italic>Radopholus similis</italic> [NCBI BioProject:PRJNA541590], <italic>Globodera rostochiensis</italic> [ENA Project: PRJEB13504], <italic>Globodera pallida</italic> [ENA Project: PRJEB123]). The read coverage and GC content of each contig were used to visualize the assembly. Reads that belonged to contigs that were identified as Nematoda or had no assigned identity were retained and used to assemble the <italic>H. carotae</italic> genome using de novo assembler SPAdes v. 3.15.3 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_001">Bankevich <italic>et al</italic>., 2012</a>). This second assembly was again subjected to the Blob Tools BLAST-based filtering of reads to remove any remaining contamination before reassembly using SPAdes to achieve the final assembly of the <italic>H. carotae</italic> genome. The final assembly was visualized using each contig's coverage, GC content, and phylum identity (<a ref-type="fig" href="#j_jofnem-2022-0014_fig_001">Supplementary Fig. 1</a>). Using Pilon v. 1.22 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_024">Walker <italic>et al</italic>., 2014</a>) gap-filling, mis-assembly correction, base correction, and scaffolding were performed on the final assembly. The genome statistics of the corrected version of the final assembly were calculated using QUAST (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_010">Gurevich <italic>et al</italic>., 2013</a>).</p>
<p>The goal of this genome sequencing effort was to rapidly provide usable data that was of sufficient quality to the nematology community to expand genomic resources available for cyst nematodes. The scaffolded assembly of the <italic>H. carotae</italic> genome was 95,118,078 bp in 17,839 scaffolds (<a ref-type="table" href="#j_jofnem-2022-0014_tab_001">Table 1</a>). Two other <italic>Heterodera</italic> species have publicly available genomes, <italic>H. glycines</italic> and <italic>H. schachtii</italic>, which are 1.66 and 1.88X larger than the <italic>H. carotae</italic> genome, respectively (<a ref-type="table" href="#j_jofnem-2022-0014_tab_001">Table 1</a>). Unlike <italic>H. carotae</italic>, the genomes for <italic>H. glycines</italic> and <italic>H. schachtii</italic> were assembled in fewer segments and consisted of 9 and 395 scaffolds, respectively. Using the raw data as input and a kmer size of 21, GenomeScope (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_023">Vurture <italic>et al</italic>., 2017</a>) was used to determine the potential genome size and repeat length for the <italic>H. carotae</italic> genome which were estimated at 120 Mb and 81 Mb, respectively. The level of duplication and the number of repetitive regions is unclear in the <italic>H. schachtii</italic> genome. The <italic>H. glycines</italic> genome contains 34% repeated regions and 18.7 Mb of tandem duplicates (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_017">Masonbrink <italic>et al</italic>., 2019</a>), indicating that the fragmented assembly of <italic>H. carotae</italic> may be due in part to the inability of Illumina sequencing to resolve repetitive and low complexity portions of the genome. Additionally, some repetitive regions of the genome may have been lost due to the high stringency standards used when BLAST filtering for contamination. The N50 for <italic>H. glycines</italic> and <italic>H. schachtii</italic> are 17 Mb and 1.2 Mb, respectively, whereas the N50 for <italic>H. carotae</italic> is 13,935 bp (<a ref-type="table" href="#j_jofnem-2022-0014_tab_001">Table 1</a>). The GC content of each genome is relatively similar with 39.39%, 36.66%, and 33.23% for <italic>H. carotae</italic>, <italic>H. glycines</italic>, and <italic>H. schachtii,</italic> respectively (<a ref-type="table" href="#j_jofnem-2022-0014_tab_001">Table 1</a>).</p>
<table-wrap id="j_jofnem-2022-0014_tab_001" position="float"><label>Table 1</label><caption><p>Comparison of the genome assembly statistic of <italic>Heterodera</italic> species.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top"><bold>Assembly statistic</bold></th>
<th align="center" valign="top"><bold><italic>H. glycines</italic> (PRJNA381081)<a ref-type="table-fn" href="#j_jofnem-2022-0014_table-fn_001"><sup>a</sup></a></bold></th>
<th align="center" valign="top"><bold><italic>H. schachtii</italic> (PRJNA722882)<a ref-type="table-fn" href="#j_jofnem-2022-0014_table-fn_001"><sup>a</sup></a></bold></th>
<th align="center" valign="top"><bold><italic>H. carotae</italic> (PRJNA774818)<a ref-type="table-fn" href="#j_jofnem-2022-0014_table-fn_001"><sup>a</sup></a></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">Size (bp)</td>
<td align="center" valign="top">157,978,452</td>
<td align="center" valign="top">179,246,932</td>
<td align="center" valign="top">95,118,078</td>
</tr>
<tr>
<td align="left" valign="top">Number of scaffolds</td>
<td align="center" valign="top">9</td>
<td align="center" valign="top">395</td>
<td align="center" valign="top">17,839</td>
</tr>
<tr>
<td align="left" valign="top">Largest scaffold (bp)</td>
<td align="center" valign="top">23,985,585</td>
<td align="center" valign="top">6,046,013</td>
<td align="center" valign="top">113,425</td>
</tr>
<tr>
<td align="left" valign="top">GC content (%)</td>
<td align="center" valign="top">36.66</td>
<td align="center" valign="top">33.23</td>
<td align="center" valign="top">39.39</td>
</tr>
<tr>
<td align="left" valign="top">N50 value (bp)</td>
<td align="center" valign="top">17,907,690</td>
<td align="center" valign="top">1,273,070</td>
<td align="center" valign="top">13,935</td>
</tr>
<tr>
<td align="left" valign="top">No. contigs >5,000 bp</td>
<td align="center" valign="top">9</td>
<td align="center" valign="top">395</td>
<td align="center" valign="top">5,030</td>
</tr>
<tr>
<td align="left" valign="top">No. contigs >10,000 bp</td>
<td align="center" valign="top">9</td>
<td align="center" valign="top">359</td>
<td align="center" valign="top">2,755</td>
</tr>
<tr>
<td align="left" valign="top">No. contigs >25,000 bp</td>
<td align="center" valign="top">9</td>
<td align="center" valign="top">309</td>
<td align="center" valign="top">699</td>
</tr>
<tr>
<td align="left" valign="top">No. contigs >50,000 bp</td>
<td align="center" valign="top">9</td>
<td align="center" valign="top">269</td>
<td align="center" valign="top">103</td>
</tr>
<tr>
<td align="center" valign="top" colspan="4"><hr/></td>
</tr>
<tr>
<td align="left" valign="top">Complete BUSCOs (%)</td>
<td align="center" valign="top">1,400 (44.7)</td>
<td align="center" valign="top">1,422 (45.4)</td>
<td align="center" valign="top">1,259 (40.2)</td>
</tr>
<tr>
<td align="left" valign="top">Complete and single-copy BUSCOs (%)</td>
<td align="center" valign="top">1,291 (41.2)</td>
<td align="center" valign="top">1,372 (43.8)</td>
<td align="center" valign="top">1,238 (39.5)</td>
</tr>
<tr>
<td align="left" valign="top">Complete and duplicated BUSCOs (%)</td>
<td align="center" valign="top">109 (3.5)</td>
<td align="center" valign="top">50 (1.6)</td>
<td align="center" valign="top">21 (0.7)</td>
</tr>
<tr>
<td align="left" valign="top">Fragmented BUSCOs (%)</td>
<td align="center" valign="top">109 (3.5)</td>
<td align="center" valign="top">120 (3.8)</td>
<td align="center" valign="top">143 (4.6)</td>
</tr>
<tr>
<td align="left" valign="top">Missing BUSCOs (%)</td>
<td align="center" valign="top">1,622 (51.8)</td>
<td align="center" valign="top">1,589 (50.8)</td>
<td align="center" valign="top">1,729 (55.2)</td>
</tr>
<tr>
<td align="center" valign="top" colspan="4"><hr/></td>
</tr>
<tr>
<td align="left" valign="top">Predicted protein coding genes</td>
<td align="center" valign="top">29,679</td>
<td align="center" valign="top">26,768</td>
<td align="center" valign="top">17,212</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn-group>
<fn id="j_jofnem-2022-0014_table-fn_001"><label>a</label><p>NCBI GenBank BioProject accession number.</p></fn>
</fn-group>
</table-wrap-foot>
</table-wrap>
<p>To determine the completeness of the <italic>H. carotae</italic> genome BUSCO v.5.2.2 was used (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_020">Simão <italic>et al</italic>., 2015</a>). BUSCO was also run on the <italic>H. glycines</italic> and <italic>H. schachtii</italic> genomes for comparison (<a ref-type="table" href="#j_jofnem-2022-0014_tab_001">Table 1</a>). The <italic>H. carotae</italic> genome was ~5% less complete than the <italic>H. glycines</italic> and <italic>H. schachtii</italic> genomes. All three genomes shared 927 complete genes and 1,391 missing BUSCO genes. Across all three genomes, there are 1,613 complete BUSCO genes. In <a ref-type="fig" href="#j_jofnem-2022-0014_fig_002">Supplementary Fig. 2</a>, a Venn diagram can be found depicting the overlap of complete BUSCO genes across the three <italic>Heterodera</italic> species. Although all three genomes had complete BUSCO scores below 50% they are in line with other plant-parasitic nematode genomes that range in completeness from 40% to 60% (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_012">Howe <italic>et al</italic>., 2016</a>, <a ref-type="bibr" href="#j_jofnem-2022-0014_ref_013">2017</a>).</p>
<p>In addition to confirming specimen identity using morphometrics, <italic>coxI</italic> and <italic>hsp90</italic> were used to place the specimen in a phylogenetic context. All available sequences of <italic>coxI</italic> for <italic>H. carotae</italic> and <italic>Heterodera cruciferae</italic> were obtained from NCBI, along with a random selection of <italic>coxI</italic> accessions from 31 other <italic>Heterodera</italic> species, and five accessions of <italic>coxI</italic> from <italic>Meloidogyne</italic> species as an outgroup. All NCBI accessions used are denoted in <a ref-type="fig" href="#j_jofnem-2022-0014_fig_003">Supplementary Fig. 3</a>. The <italic>coxI</italic> sequence was retrieved from the <italic>H. carotae</italic> assembly using usearch11 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_005">Edgar, 2010</a>) with <italic>H. carotae</italic> isolate U45 (NCBI Accession KX463301.1) as the search sequence. All sequences were aligned using MUSCLE v.3.8.31 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_004">Edgar, 2004</a>), Newick phylogenetic trees were created using FastTree v.2.1.10 with default parameters (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_019">Price <italic>et al</italic>., 2010</a>), and trees were visualized using TreeGraph 2 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_021">Stöver and Müller, 2010</a>). As shown in <a ref-type="fig" href="#j_jofnem-2022-0014_fig_003">Supplementary Fig. 3</a>, <italic>H. carotae</italic> from this study formed a clade with three other <italic>H. carotae</italic> isolates collected from carrot in the Tepeaca Valley, Puebla, Mexico (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_006">Escobar-Avila <italic>et al</italic>., 2018</a>) supported by a 0.836 bootstrap value.</p>
<p>The same process described above was used to examine the phylogenetics of <italic>hsp90</italic> from the Chilean <italic>H. carotae. H. carotae</italic> clone 2449 (NCBI Accession JQ316187.1) was used as the search sequence to retrieve <italic>hsp90</italic> from the <italic>H. carotae</italic> genome. <italic>Hsp90</italic> sequences were obtained from NCBI for phylogenetic comparison to the <italic>H. carotae</italic> in this study from four <italic>Globodera</italic> species, six <italic>Heterodera</italic> species, all available <italic>H. carotae</italic> accessions, two <italic>Cactodera</italic> species, and six <italic>Meloidogyne</italic> species as an outgroup. NCBI accessions and species are denoted in <a ref-type="fig" href="#j_jofnem-2022-0014_fig_004">Supplementary Fig. 4</a>. Similar to the <italic>coxI</italic> results, the <italic>H. carotae</italic> from Chile fell into a clade separate from other <italic>Heterodera</italic> species but with other <italic>H. carotae</italic> isolates, supported by a bootstrap value of 0.943 (<a ref-type="fig" href="#j_jofnem-2022-0014_fig_004">Supplementary Fig. 4</a>).</p>
<p>The <italic>H. carotae</italic> genome was annotated using BRAKER v. 2.1.5 (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_011">Hoff <italic>et al</italic>., 2016</a>) with protein hints from <italic>H. glycines</italic> (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_017">Masonbrink <italic>et al</italic>., 2019</a>). In order to compare across other <italic>Heterodera</italic> species genomes, <italic>H. schachtii</italic> was also annotated with the same parameters as <italic>H. carotae</italic>. <italic>H. schachtii</italic> and <italic>H. glycines</italic> have 9,556 and 12,467 more genes, respectively, than <italic>H. carotae,</italic> which has 17,212 protein coding genes identified in the genome (<a ref-type="table" href="#j_jofnem-2022-0014_tab_001">Table 1</a>). In the <italic>H. glycines</italic> genome there is a strong colocalization of effector genes in highly duplicated scaffolds (<a ref-type="bibr" href="#j_jofnem-2022-0014_ref_017">Masonbrink <italic>et al</italic>., 2019</a>). Some protein coding genes could be lost in <italic>H. carotae</italic> genome due to the collapse of repeated and duplicated regions in the current assembly. Additionally, <italic>H. glycines</italic> was annotated with the addition of RNAseq data, which is not available for <italic>H. carotae</italic>. The addition of RNAseq data could both improve the predicted proteins calls in the <italic>H. carotae</italic> annotation and increase the number of predicted proteins.</p>
<p>The <italic>H. carotae</italic> genome assembly, annotation, and raw data can be found under the NCBI BioProject PRJNA774818. This version of the <italic>H. carotae</italic> genome provides a starting point for further investigation into the basic biology of <italic>H. carotae</italic> and a resource for the greater cyst nematode genomics community. Further sequencing and improvements to the <italic>H. carotae</italic> genome would increase its usefulness; however, in its current state this genome offers researchers a resource for discovering novel <italic>H. carotae</italic> effectors, developing diagnostic markers, or examining genomic similarities across cyst species.</p>
</div>

Genome Announcement: The Draft Genome of the Carrot Cyst Nematode <italic>Heterodera carotae</italic>

Journal of Nematology

This work is licensed under the Creative Commons Attribution 4.0 International License.

Genome Announcement: The Draft Genome of the Carrot Cyst Nematode 

Heterodera carotae, the carrot cyst nematode, is a significant pest affecting carrot globally. Here we present the draft genome of H. carotae, which was...

Genome Announcement: The Draft Genome of the Carrot Cyst Nematode Heterodera carotae

Article Category: Research Paper

Data publikacji: 29 maj 2022

Zakres stron: -

Otrzymano: 21 sty 2022

DOI: https://doi.org/10.2478/jofnem-2022-0014

Słowa kluczowe
carrot cyst nematode, draft genome, genomics,

© 2022 Catherine L. Wram et al., published by Sciendo

This work is licensed under the Creative Commons Attribution 4.0 International License.

Genome Announcement: The Draft Genome of the Carrot Cyst Nematode Heterodera carotae

Article Category: Research Paper

Data publikacji: 29 maj 2022

Zakres stron: -

Otrzymano: 21 sty 2022

DOI: https://doi.org/10.2478/jofnem-2022-0014

Słowa kluczowecarrot cyst nematode, draft genome, genomics,

© 2022 Catherine L. Wram et al., published by Sciendo

This work is licensed under the Creative Commons Attribution 4.0 International License.

Słowa kluczowe
carrot cyst nematode, draft genome, genomics,