# A draft genome for a species of Halicephalobus
(Panagrolaimidae)

Erik J Ragsdale, Georgios Koutsovoulos, Joseph F Biddle

► To cite this version:

Erik J Ragsdale, Georgios Koutsovoulos, Joseph F Biddle. A draft genome for a species of Halicephalobus (Panagrolaimidae). Journal of Nematology, 2019, 51 (68), pp.1-4. 10.21307/jofnem-2019-068. hal-02529176

HAL Id: hal-02529176

https://hal.inrae.fr/hal-02529176

Submitted on 2 Apr 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A draft genome for a species of *Halicephalobus*
(Panagrolaimidae)

Erik J. Ragsdale,¹*
Georgios Koutsovoulos² and
Joseph F. Biddle¹

¹Department of Biology, Indiana
University, 915 E. 3rd St.,
Bloomington, IN, 47405.

²Institute national de la recherche
agronomique, Sophia-Antipolis,
PACA, France.

*E-mail: ragsdale@indiana.edu

This paper was edited by
Dee Denver.

Received for publication June 19,
2019.

## Abstract

*Halicephalobus* is a clade of small, exclusively parthenogenic nematodes that have sometimes colonized remarkable habitats. Given their phylogenetic closeness to other parthenogenic panagrolaimid species with which they likely share a sexually reproducing ancestor, *Halicephalobus* species provide a point of comparison for parallelisms in the evolution of asexuality. Here, we present a draft genome of a putatively new species of *Halicephalobus* isolated from termites in Japan.

## Keywords

Comparative genomics, Free-living, Panagrolaimorpha, Parthenogenesis.

Panagrolaimidae are a nematode family represented by lifestyles under the most extreme environments known to animals, including desiccation (Ricci and Pagani 1997; Treonis and Wall, 2005), freezing (Wharton and Barclay, 1993), and an unusual range of pH conditions (Peters, 1928). Phylogenetically close families that include important parasites of insects and vertebrates, Steinernematidae and Strongyloididae, respectively, have enjoyed considerable attention from sequencing efforts. However, free-living species of Panagrolaimidae have lagged behind, as have free-living species of Clade IV (sensu Blaxter et al., 1998) more generally. A major contribution to the genomics of Panagrolaimidae was the genome sequence of *Panagrellus redivivus*, a species championed for decades as a laboratory model for genetics and development, especially as a satellite to *Caenorhabditis elegans* (Srinivasan et al., 2013). More recently, several species of *Panagrolaimus* have been sequenced to study the evolutionary signatures of parthenogenesis and tolerance of extreme environments (Schiffer et al. under review). To provide resources for another member of this nematode family, we present a draft genome for a strain of *Halicephalobus*.

*Halicephalobus* species are minute nematodes that often inhabit soil-like environments and are particularly well-known from organic-rich substrates (Andrássy, 1984; Anderson et al., 1998; Steel et al., 2012). Also,

like other panagrolaimids, these nematodes have been reported from demanding, unusual habitats. For example, *H. gingivalis*, typically a resident of soil and humus, is an opportunistic, blood-inhabiting pathogen of mammals, especially horses and occasionally humans (Anderson and Bemrick, 1965; Hoogstraten and Young, 1975; Blunden et al., 1987). Another species of this genus, *H. mephisto*, was discovered on microbial biofilms at depths of over three kilometers within the earth's crust (Borgonie et al., 2011). The ability of *Halicephalobus* species to colonize such habitats suggests that comparative genomics of the genus may yield clues into the mechanisms enabling unusual lifestyles. It is possible that the mode of reproduction in all known *Halicephalobus* species, obligate parthenogenesis, also aids in the colonization of new, sometimes extreme habitats, by requiring only a single individual to found a population. Because the common ancestor of *Panagrolaimus* and *Halicephalobus* was most likely gonochoristic (dioecious) (Lewis et al., 2009), parthenogenesis may have evolved twice independently in these two genera. Consequently, the genome for a *Halicephalobus* species may reveal illuminating parallelisms in how asexuality influences genomic evolution.

Here, we have sequenced the genome of a *Halicephalobus* strain (NKZ332) isolated from termites in Japan. Because of the paucity of reliable diagnostic
morphology for the genus (Anderson et al., 1998), combined with the inability to test species boundaries through crossing experiments, assignment of a valid species binomial to this or other species of the genus is not infallible. Indeed, taxonomic revision of this genus will rely on molecular data (Nadler et al., 2003), interpreted in terms of evolutionary or phylogenetic species concepts (e.g., Borgonie et al., 2011). Based on bionomic evidence, specifically the nematodes' colonization of rotting wood, the sequenced isolate herein would be best assigned to *H. similigaster*, a morphospecies known from Europe (Andrássy, 1952; Köhler, 2011). However, the present strain shows a sharp distinction in rRNA markers, specifically <95% similarity, from all other sequenced *Halicephalobus* isolates, particularly those from similar habitats in other locations (Foley et al., 2018). This similarity is much less than that between biological species in other nematode groups (Kanzaki et al., 2013; Félix et al., 2014). Given its molecular and geographic distinctness from other known *Halicephalobus* strains and morphospecies, the present strain likely has an evolutionary fate (sensu Wiley, 1978) separate from other nominal species. Therefore, we anticipate that molecular comparisons of *Halicephalobus* strains from both Europe and East Asia will support and delimit *Halicephalobus* sp. NKZ332 as a new species.

To sequence the genome of *Halicephalobus* sp. NKZ332, which we collected from xenic bacterial culture, we used the Illumina NextSeq platform to generate 75 million, 300-bp, paired-end reads of 150-bp and 500-bp insert sizes. Reads were trimmed with Cutadapt (Martin, 2011) and error-corrected with Reckoner (Długosz and Deorowicz, 2017). We first used Minia (Drezen et al., 2014) to perform a preliminary single-end assembly, with which we validated insert size and screened for possible non-nematode contaminants using Blobtools (Kumar et al., 2013). We identified two contaminants, *Stenotophomonas maltophilia* at high coverage and *Pseudomonas* sp. at low coverage, both with a higher GC content than the predicted nematode scaffolds' GC content (38%). We then assembled the genome using the SPAdes assembler (Bankevich et al., 2012) and removed scaffolds that were identified as contamination, at low coverage (<10×), or less than 500 bp in length. We made a first-pass annotation with the MAKER pipeline using evidence drawn from the Swiss-Prot database, the hidden Markov model (HMM) profile computed from GeneMark (Lukashin and Borodovsky, 1998), and the Augustus (Stanke and Waack, 2003) profile generated with BUSCO (Waterhouse et al., 2018). We then used the predicted genes from the MAKER output to train and run Augustus to produce the final annotation.

Our assembly of *Halicephalobus* sp. NKZ332 spans 47 Mb in 5,085 contigs with an N50 of 60 kb. This assembly size is smaller than the genomes of most sequenced nematodes (International Helminth Genomes Consortium, 2018). Likewise, the assembly is much smaller than those of other *Panagrolaimus* species, whose genomes are ~140 to ~180 Mb in gonochoristic species and even larger (~230–270 Mb) in parthenogenic strains, which likely originated from an allopolyploidy event (Schiffer et al. under review). In our annotation, we predicted a total of 11,023 genes with a mean of 1,511 bp in length. We assessed the completeness of the genome assembly using BUSCO (Simão et al. 2015), with which our annotation shared 83.2% (79.8% complete genes) of the nematode ortholog set. Given such representation of conserved orthologs in the annotation, the genome of *Halicephalobus* sp. NKZ332 is strikingly minimal. Finally, using the MITOS Web Server (Bernt et al. 2013), we assembled the mitochondrial genome sequence, which has a length of 13,886 bp and GC content of 20%.

In summary, we present a compact nematode genome representing a genus of small, parthenogenic, and ecologically opportunistic nematodes. The nuclear and mitochondrial genome sequences have been deposited in NCBI GenBank under accession numbers VOSG00000000 and MN207311, respectively. Raw reads have been deposited in the NCBI Sequence Read Archive (SAMN12324525).

## Acknowledgments

The authors thank the Indiana University Center for Genomics and Bioinformatics for genomic library preparation and sequencing. The authors also thank Natsumi Kanzaki, who collected the isolate sequenced by the authors, for providing the strain.

## References

Anderson, R. C., Linder, K. E. and Peregrine, A. S. 1998. *Halicephalobus gingivalis* (Stefanski, 1954) from a fatal infection in a horse in Ontario, Canada with comments on the validity of *H. deletrix* and a review of the genus. Parasite 5:255–61.

Anderson, R. V. and Bemrick, W. J. 1965. *Micronema deletrix*, a saprophagous nematode inhabiting a nasal tumor of a horse. Proceedings of the Helminthological Society of Washington 32:74–5.

Andrássy, I. 1952. Freilebende Nematoden aus dem Bükk-Gebirge. Annales historico-naturales Musei Nationalis Hungarici 2:13–65.
Andrássy, I. 1984. Klasse Nematoda, Gustav Fischer Verlag, Stuttgart.

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., Lesin, V. M., Nikolenko, S. I., Pham, S., Prjibelski, A. D., Pyshkin, A. V., Sirotkin, A. V., Vyahhi, N., Tesler, G., Alekseyev, M. A. and Pevzner, P. A. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. Journal of Computational Biology 19:455-77.

Bernt, M., Donath, A., Jühling, F., Externbrink, F., Florentz, C., Fritzsch, G., Pütz, J., Middendorf, M. and Stadler, P. F. 2013. MITOS: improved *de novo* metazoan mitochondrial genome annotation. Molecular Phylogenetics and Evolution 69:313-19.

Blaxter, M., De Ley, P., Garey, J. R., Liu, L. X., Scheldeman, P., Vierstraete, A., Vanfleteren, J. R., Mackey, L. Y., Dorris, M., Frisse, L. M., Vida, J. T. and Thomas, W. K. 1998. A molecular evolutionary framework for the phylum Nematoda. Nature 392:71-5.

Blunden, A. S., Khalil, L. F. and Webbon, P. M. 1987. *Halicephalobus deletrix* infection in a horse. Equine Veterinary Journal 19:255-60.

Borgonie, G., Garcia-Moyano, A., Litthauer, D., Bert, W., Bester, A., van Heerden, E., Moller, C., Erasmus, M. and Onstott, T. C. 2011. Nematoda from the terrestrial deep subsurface of South Africa. Nature 474:79-82.

Długosz, M. and Deorowicz, S. 2017. RECKONER: read error corrector based on KMC. Bioinformatics 33:1086-9.

Drezen, E., Rizk, G., Chikhi, R., Deltel, C., Lemaitre, C., Peterlongo, P. and Lavenier, D. 2014. GATB: genome assembly & analysis tool box. Bioinformatics 30:2959-61.

Félix, M. A., Braendle, C. and Cutter, A. D. 2014. A streamlined system for species diagnosis in *Caenorhabditis* (Nematoda: Rhabditidae) with name designations for 15 distinct biological species. PLoS One 9:e0118327.

Foley, J. R., Chouvenc, T., Giblin-Davis, R. M., Su, N. Y. and Kanzaki, N. 2018. Phoresy and within-colony transmission of nematodes associated with alates of subterranean termites (Isoptera: Rhinotermitidae). Environmental Entomology 47:1107-16.

Hoogstraten, J. and Young, W. G. 1975. Meningoencephalomyelitis due to the saprophagous nematode, *Micronema deletrix*. The Canadian Journal of Neurological Sciences, 2:121-6.

International Helminth Genomes Consortium. 2018. Comparative genomics of the major parasitic worms. Nature Genetics. 51:163-74.

Kanzaki, N., Ragsdale, E. J., Herrmann, M., Susoy, V. and Sommer, R. J. 2013. Two androdioecious and one dioecious new species of *Pristionchus* (Nematoda: Diplogastridae): new reference points for the evolution of reproductive mode. Journal of Nematology 45:172-94.

Köhler, A. 2011. Description of *Halicephalobus similigaster* (Andrássy, 1952) and a comparison with

its junior synonyms *H. minutus* (Körner, 1954) syn. n. and *H. parvus* (Körner, 1954) syn. n. (Nematoda: Panagrolaimidae). Nematology 13:287-93.

Kumar, S., Jones, M., Koutsovoulos, G., Clarke, M. and Blaxter, M. 2013. Blobology: exploring raw genome data for contaminants, symbionts and parasites using taxon-annotated GC-coverage plots. Frontiers in Genetics 4:237.

Lewis, S. C., Dyal, L. A., Hilburn, C. F., Weitz, S., Liau, W. S., Lamunyon, C. W. and Denver, D. R. 2009. Molecular evolution in *Panagrolaimus* nematodes: origins of parthenogenesis, hermaphroditism and the Antarctic species *P. davidi*. BMC Evolutionary Biology 9:15.

Lukashin, A. V. and Borodovsky, M. 1998. GeneMark.hmm: new solutions for gene finding. Nucleic Acids Research 26:1107-15.

Martin, M. 2011. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17:10-2.

Nadler, S. A., Carreno, R. A., Adams, B. J., Kinde, H., Baldwin, J. G. and Mundo-Ocampo, M. 2003. Molecular phylogenetics and diagnosis of soil and clinical isolates of *Halicephalobus gingivalis* (Nematoda: Cephalobina: Panagrolaimoidea), an opportunistic pathogen of horses. International Journal for Parasitology 33:1115-25.

Peters, B. G. 1928. On the bionomics of the vinegar eelworm. Journal of Helminthology 6:1-38.

Ricci, C. and Pagani, M. 1997. Desiccation of *Panagrolaimus rigidus* (Nematoda): survival, reproduction and the influence on the internal clock. Hydrobiologia 347:1-13.

Schiffer, P. H., Danchin, E. G. J., Burnell, A. M., Schiffer, A.-M., Creevey, C. J., Wong, S., Dix, I., O'Mahony, G., Culleton, B. A., Rencurel, C., Stier, G., Martínez-Salazar, A., Marconi, A., Trivedi, U., Kroiher, M., Thorne, M. A. S., Schierenberg, E., Wiehe, T. and Blaxter, M. Under review. Signatures of the evolution of parthenogenesis and cryptobiosis in panagrolaimid nematodes. bioRxiv 159152.

Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. and Zdobnov, E. M. 2015. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31(19):3210-2.

Srinivasan, J., Dillman, A. R., Macchietto, M. G., Heikkinen, L., Lakso, M., Fracchia, K. M., Antoshechkin, I., Mortazavi, A., Wong, G. and Sternberg, P. W. 2013. The draft genome and transcriptome of *Panagrellus redivivus* are shaped by the harsh demands of a free-living lifestyle. Genetics 193:1279-95.

Stanke, M. and Waack, S. 2003. Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics 19:ii215-ii225.

Steel, H., Vandecasteele, B., Willekens, K., Sabbe, K., Moens, T. and Bert, W. 2012. Nematode
communities and macronutrients in composts and
compost-amended soils as affected by feedstock
composition. Applied Soil Ecology 61:100–12.

Treonis, A. M. and Wall, D. H. 2005. Soil nematodes and desiccation survival in the extreme arid environment of the Antarctic Dry Valleys. Integrative and Comparative Biology 45:741–50.

Waterhouse, R. M., Seppey, M., Simão, F. A., Manni, M., Ioannidis, P., Klioutchnikov, G., Kriventseva, E. V. and

Zdobnov, E. M. 2018. BUSCO applications from quality
assessments to gene prediction and phylogenomics.
Molecular Biology and Evolution 35:543–48.

Wharton, D. A. and Barclay, S. 1993. Anhydrobiosis in the free-living antarctic nematode *Panagrolaimus davidi* (Nematoda: Rhabditida). Fundamental and Applied Nematology 16:17–22.

Wiley, E. O. 1978. The evolutionary species concept reconsidered. Systematic Zoology 27:17–26.