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

<p><italic>Bradynema listronoti</italic> is an insect-parasitic nematode known to infect the carrot weevil, <italic>Listronotus oregonensis</italic>. We present the first sequence for this species and for any Allantonematidae, produced with a combination of short and long reads. The draft genome of <italic>B. listronoti</italic> is 80.6 Mb in size, assembled in 152 scaffolds.</p></abstract>

Bradynema listronoti is an insect-parasitic nematode known to infect the carrot weevil, Listronotus oregonensis. We present the first sequence for this species and for any Allantonematidae, produced with a combination of short and long reads. The draft genome of B. listronoti is 80.6 Mb in size, assembled in 152 scaffolds.

<div xmlns="http://www.w3.org/1999/xhtml">
<p>The genus <italic>Bradynema</italic> (Nematoda: Allantonematidae), first characterized by Strassen in 1892, currently comprises nine species of insect-parasitic nematodes with a diverse set of hosts. In 2007, a novel species, <italic>Bradynema listronoti</italic>, was isolated from the hemocoel of adult carrot weevils, <italic>Listronotus oregonensis</italic> (Coleoptera: Curculionidae), an important pest of horticultural crops (<a ref-type="bibr" href="#j_jofnem-2022-0047_ref_022">Zeng <italic>et al</italic>., (2007)</a>). It was shown that <italic>B. listronoti</italic> castrates its host by inhibiting the maturation of the reproductive system in females (<a ref-type="bibr" href="#j_jofnem-2022-0047_ref_006">Gagnon <italic>et al</italic>., (2007)</a> 2019), and its potential as a biological control agent is currently investigated. In comparison to entomopathogenic nematodes, this group has been studied very little, despite its importance in the regulation of insect populations. In this study, we combined short- and long-read sequencing data to assemble a high-quality genome for <italic>B. listronoti</italic>.</p>
<p><italic>B. listronoti</italic> specimens (strain AAFC-BrLi-01) used in this study were isolated from a colony of infected carrot weevils at the Saint-Jean-sur-Richelieu Research and Development Centre of Agriculture and Agri-Food Canada. This colony was established following the initial isolation of <italic>B. listronoti</italic> by <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_022">Zeng <italic>et al</italic>. (2007)</a>. Genomic DNA was obtained from an individual juvenile (J3) nematode and a pool of J3. For the single nematode, DNA was extracted and amplified using the REPLI-g Advanced DNA Single Cell Kit (Qiagen, Hilden, Germany). To resolve the hyperbranching generated by this whole-genome amplification (WGA) protocol, 5 μg of WGA gDNA was treated with T7 endonuclease, and short fragments were removed using the Circulomics SRE XS Kit (Circulomics, Baltimore, MD). From this size-selected material, two Nanopore MinION libraries were generated, from 1 μg of genomic material each, using the Ligation Sequencing Kit chemistry (SQK-LSK109, Oxford Nanopore Technologies, Oxford, United Kingdom). Short fragments were eliminated using the L Fragment Buffer provided in the kit prior to sequencing on MinION R9.4.1 flow cells. Using the same WGA gDNA, a short-read shotgun sequencing library was prepared using the NEBNext Ultra II DNA Library Preparation Kit and paired-end-sequenced on the Illumina NovaSeq6000 platform (2 × 150 bp) at the McGill University and Génome Québec Innovation Centre. In a second experiment, high-molecular weight (HMW) gDNA was prepared from a pool of ~15,000 J3 nematodes for downstream assembly scaffolding purposes. Genomic DNA was extracted using the Monarch HMW DNA Extraction Kit for Tissue (New England Biolabs, Ipswich, MA) following the manufacturer’s recommendations. The sequencing library was prepared from 1 μg of genomic material using the Ligation Sequencing Kit (Q20+) chemistry (SQK-LSK112) and sequenced on a MinION R10.4 flow cell. Messenger RNA was extracted and amplified from multiple <italic>B. listronoti</italic> stages (egg, J2, J3, and J4) using the method described in <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_018">Ste-Croix <italic>et al</italic>. (2021)</a> with three biological replicates of each. Sequencing libraries were prepared using the NEBNext Ultra II DNA Library Preparation Kit and paired-end-sequenced on the NovaSeq6000 platform (2 <sub>×</sub> 150 bp). Genome size was estimated using GenomeScope (<a ref-type="bibr" href="#j_jofnem-2022-0047_ref_019">Vurture <italic>et al</italic>., 2017)</a> with k-mer frequency distributions derived from single-nematode whole-genome shotgun sequencing data.</p>
<p>The different sequencing technologies generated 33 Gb of long-read sequences alongside 5.3 Gb of Q20+ long-read sequences derived from HMW gDNA, and 20.7 Gb of paired-end short reads. Additionally, mRNA sequencing generated 706 M reads for an estimated 210 Gb of data. For long- and short-read sequencing data, sequences were trimmed and filtered based on quality using Guppy (v 6.0.6) and Fastp (v 0.20.1 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_005">Chen <italic>et al</italic>., 2018)</a>, respectively. The genome was assembled using sequencing data generated from a single nematode to reduce the level of heterozygosity and facilitate assembly. First, long-read sequences were assembled using Flye (v 2.9-b1768 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_010">Kolmogorov <italic>et al</italic>., 2019)</a> in default setting. Second, short-reads were assembled using HASLR (v 0.8a1 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_007">Haghshenas <italic>et al</italic>., 2020)</a> assisted by long-read sequences. These assemblies were then merged using Quickmerge (v 0.3) prior to scaffolding using LINKS (v 1.8.7 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_021">Warren <italic>et al</italic>., 2015)</a> with Q20+ long reads derived from the pool of J3. This assembly was further scaffolded using information provided by RNA-seq with Rascaf (v 1.0.2 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_015">Song <italic>et al</italic>., 2016)</a>. Finally, this scaffolded assembly was polished and corrected using five cycles of Pilon (v 1.23 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_020">Walker <italic>et al</italic>., 2014)</a> supplemented with single-nematode whole-genome sequencing data. The final corrected draft was analyzed for contaminating scaffolds with BlobToolKit (v 3.1.6 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_004">Challis <italic>et al</italic>., 2020),</a> and the quality and completeness of the genome were evaluated using BUSCO with the built-in Augustus gene predictor option set (v 5.3.2 – eukaryota_odb10 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_013">Manni <italic>et al</italic>., 2021).</a> The genome was then masked using RepeatMasker (v 4.1.2-p1 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_014">Smit <italic>et al</italic>., 2015)</a>, annotated using BRAKER (<a ref-type="bibr" href="#j_jofnem-2022-0047_ref_016">Stanke <italic>et al</italic>., 2006, 2008;</a> <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_011">Li <italic>et al</italic>., 2009;</a> <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_002">Barnett <italic>et al</italic>., 2011;</a> <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_012">Lomsadze <italic>et al</italic>., 2014;</a> <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_008">Hoff <italic>et al</italic>., 2016, 2019;</a> <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_003">Bruna <italic>et al</italic>., 2021)</a>, retrained using the RNA-seq data. The draft genome size of <italic>B. listronoti</italic> was 80,637,718 bp assembled into 152 scaffolds, with an average G + C content of 33.08%, an L50 of 7 contigs, and an N50 of 4,373,361 bp. Completeness analysis using eukaryotic BUSCO genes yielded a score of 75.3% (C: 75.3% [S: 61.6%, D: 13.7%], F: 9.0%, M: 15.7%, n: 255). This is slightly lower than that obtained with the genome of <italic>Deladenus siricidicola</italic> (C: 87.1% [S: 86.3%, D: 0.8%], F: 5.5%, M: 7.4%, n: 255), a phylogenetically related species of insect-parasitic nematode (assembly accession GCA_009724625.1). Also, 13.7% of these BUSCO genes were found to be duplicated, indicating that a low level of haplotypic duplication may remain. This could also explain the small difference between the genome size (80.6 Mb) and its estimation by GenomeScope (69.8 Mb). BRAKER predicted 16,831 protein-coding genes. Functional annotation was accomplished by querying these predictions against the NCBI nr database using BLASTP (v. 2.6.0 – <a ref-type="bibr" href="#j_jofnem-2022-0047_ref_003">Altschul <italic>et al</italic>., 1990)</a>. This is the first genome release for <italic>B. listronoti</italic> and for any Allantonematidae species.</p>
<sec id="j_jofnem-2022-0047_s_001"><div>Data availability and accession numbers</div>
<p>All the raw sequences and genome assembly files are available under the NCBI BioProject PRJNA842945. This whole-genome shotgun project has been deposited at DDBJ/ENA/GenBank under the accession JANLNO000000000.</p>
</sec>
</div>

Draft Genome of the Insect-Parasitic Nematode <italic>Bradynema listronoti</italic>

Journal of Nematology

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

Draft Genome of the Insect-Parasitic Nematode 

Bradynema listronoti is an insect-parasitic nematode known to infect the carrot weevil, Listronotus oregonensis. We present the first sequence for this species...

Draft Genome of the Insect-Parasitic Nematode Bradynema listronoti

Article Category: Research Paper

Published Online: Dec 14, 2022

Page range: -

Received: Jun 02, 2022

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

Keywords
Bradynema listronoti, entomoparasitic nematode, genomics, MinION, NovaSeq, whole-genome amplification

© 2022 His Majesty the King in Right of Canada. published by Sciendo

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

Draft Genome of the Insect-Parasitic Nematode Bradynema listronoti

Article Category: Research Paper

Published Online: Dec 14, 2022

Page range: -

Received: Jun 02, 2022

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

KeywordsBradynema listronoti, entomoparasitic nematode, genomics, MinION, NovaSeq, whole-genome amplification

© 2022 His Majesty the King in Right of Canada. published by Sciendo

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

Keywords
Bradynema listronoti, entomoparasitic nematode, genomics, MinION, NovaSeq, whole-genome amplification