Phylogenetic studies on three Helicotylenchus species based on 28S rDNA and mtCOI sequence data

Abstract To facilitate the process of spiral nematode species delineation, populations of Helicotylenchus canadensis, H. pseudorobustus, and H. varicaudatus deriving from various locations and diverse natural and anthropogenic environments from Poland were investigated and characterized. For the first time, 28S rDNA sequences are reported for H. canadensis and H. varicaudatus, whereas new mtCOI sequences were acquired for all three analyzed species. A Bayesian phylogenetic analysis of the 28S rDNA fragments revealed that H. canadensis and H. varicaudatus are members of a clade that is a sister group to all other Helicotylenchus species; however, the closest known sister group to H. canadensis is H. vulgaris type A. Both 28S rDNA- and mtCOI-based phylogenetic results suggest that this clade excludes H. pseudorobustus, whose most recent common ancestor with the former species was inferred to be the ancestor of all Helicotylenchus species. Moreover, within the mtCOI sequences obtained from H. pseudorobustus, unlike from the other two, a simultaneous presence of TAG and TAA codons was identified. This may indicate mitochondrial genetic code alterations or other genomic rearrangements in H. pseudorobustus.

Spiral nematodes, which belong to the Helicotylenchus spp. (Hoplolaimidae), are representatives of Tylenchomorpha, which is considered the most successful lineage of plant parasitic nematodes. They evolved by gradual transitions from fungal feeders through facultative plant parasites toward obligatory plant parasites (Holterman et al., 2017). Furthermore, the phylogenetic tree of this group reveals main evolutionary pathways, leading from ectoparasitism or migratory endoparasitism to the most specialized, sedentary endoparasitic lifestyles, exhibited by the most harmful plant-damaging nematode species. Although the majority of research focuses on the latter pests, the investigation of species that developed less sophisticated mechanisms of parasitism should not be overlooked. Analyses of molecular information hidden in their genomes bring us closer to an explanation why some groups of nematodes became more successfully adapted to plant parasitism than others.
To determine particular Helicotylenchus species, a precise evaluation of morphological and morphometrical features is used, usually after the preliminary recognition of the characteristic coiled body shape observed in the relaxed state or after death. However, the available keys do not always allow for proper species identification due to high intra-specific and minor inter-specific variability within the genus (Fortuner, 1984;Fortuner et al., 1984). The proposed diagnostic characters and features can match and overlap between closely related species. Therefore, to facilitate the species identification process as well as to delineate more phylogenetically distant species, which share the same morphology (cryptic species), the support of a molecular approach is needed.
To date, only about 20% of Helicotylenchus species has been molecularly characterized using mostly ribosomal DNA fragments (18S, ITS, 28S rDNA; GenBank resources). Mitochondrial cytochrome c oxidase subunit I (mtCOI) sequences were reported for one species, the recently described H. oleae (Palomares-Rius et al., 2018), and a cytochrome c oxidase subunit II sequence was reported solely for H. dihystera (Riepsamen et al., 2011). A combination of both the nuclear genome-derived large subunit ribosomal DNA gene (28S rDNA) and more variable genes such as mtCOI has been suggested to comprise valuable markers for subsequent phylogenetic analyses in this group Palomares-Rius et al., 2018).
The objectives of this study were as follows: to deliver molecular characteristics of three Hecicotylenchus species commonly occurring in Poland, namely, H. canadensis, H. pseudorobustus and H. varicaudatus, by the use of nuclear 28S rDNA and mtCOI data, and to further evaluate the potential of mtCOI sequences for species identification and phylogenetic study of the genus Helicotylenchus.

Materials and methods
Sampling, nematode extraction, and conservation Each sample (about 1 kg of soil) was taken to a depth of 30 cm from the root zone using a soil sampler. Nematode extraction was performed by the decantation and sieving method, followed by the centrifugal flotation method (van Bezooijen, 2006). Nematodes were killed with hot water at 60°C. The parts of the samples of selected Hoplolaimidae were fixed in TAF and designated for morphological analysis. The rest of the samples were fixed in DESS and given to molecular studies.

Morphological identification
Morphological observations and morphometrical analyses were performed using Leica light microscope with Nomarski differential interference contrast. Morphological identification was performed using identification keys and descriptions by Waseem (1961), Yuen (1964), Sher (1966), Brzeski (1998), and Andrássy (2007. Nematodes fixed in DESS were subjected to morphological vouchering and DNA amplification procedure of Yoder et al. (2006). Nematodes were identified on the temporary slides, and subsequently multifocal images were made for every specimen.

DNA extraction
After morphological identification, analyzed nematode individuals were marked with specific codes and assigned for further molecular studies. Genomic DNA from single nematode specimens was extracted using either GenElute™ Mammalian Genomic DNA Miniprep Kit (Sigma-Aldrich) or QIAamp DNA Micro Kit (Qiagen) according to the manufacturers' instructions. For each sample, DNA was eluted in 30 μ l H 2 O. Extracted DNA was stored at −20°C.

Primers, DNA amplification, and sequencing
The mitochondrial DNA fragments of the mtCOI as well as the genomic DNA fragments of the large subunit rDNA (28S rDNA) were amplified from the collected Hoplolaimidae species. Initially, during the ongoing survey project, the publicly available JB3 and JB4 or JB5 mtCOI primers (Hu et al., 2002;Derycke et al., 2010) were tested for amplification of mtCOI fragments from numerous plant parasitic species. However, from these primers, PCR amplification failed. Therefore, a new set of primers was developed according to the slightly modified primer design methodology, as in Rybarczyk-Mydłowska et al. (2012). An alignment comprising of mtCOI publicly available sequences from different nematode taxa served as a starting point for identification of the most conserved regions, which were subsequently used for design of various variants of forward and reverse primers. GenBank accession numbers of nematode mtCOI sequences used in the alignment are listed in the Supplementary Table 1 (Table S1). Primer combinations that worked best for genus Helicotylenchus and resulted in a successful amplification of partial mtCOI sequences are listed in Table 1. It is worth mentioning that, in their recent work, Palomares-Rius et al. (2018) used a primer pair originally developed by Kanzaki and Futai (2002), which allowed them to amplify approximately 660 bp, covering a broader region of this gene from H. oleae. Remarkably, the use of the hereby presented M3.5F and M8aR primer combination targets a mtCOI region of comparable length of 670, slightly shifted to the 3′ end.
28S rDNA fragments were acquired using MCB1F and MCB1R primers (Dobosz et al., 2013). All primers used in this study were supplied either by Oligo.pl (Warsaw, PL) or by Sigma-Aldrich (USA). The amplification of the 28S rDNA and mtCOI fragments was performed in reactions containing 12.5 μ l JumpStart Taq Ready Mix or RedTaq Ready Mix (Sigma-Aldrich), 1 μ l of corresponding primer (5 μ M), 1 to 3 μ l DNA template and H 2 O with a total volume of 25 μ l. Annealing temperatures used for the amplification of 28S rDNA and mtCOI fragments were 58 to 60°C and 40 to 43°C, respectively. All PCR reactions were performed in Veriti 96-Well Thermal Cycler (Applied Biosystems, Foster City, CA, USA). Amplicons were visualized by UV illumination after Midori Green (Nippon Genetics Europe, Duren, Germany) gel staining and gel electrophoresis. Excess dNTPs and unincorporated primers were removed from the PCR product using the Clean-Up Purification Kit (A&A Biotechnology, Gdynia, Poland). As a final step, the purified DNA Note: i = inosine. a This primer is a modified version of the COIR primer proposed by He et al. (2005). The primer was elongated by 4 nucleotides and two nucleotide modifications were incorporated; b This primer is a shortened and modified version of the COI-F1 primer proposed by Kanzaki and Futai (2002). The primer was shortened from 29 nt to 19 and two nucleotide modification were introduced. , 2 μ l BigDye sequencing buffer (ThermoFisher Scientific), 1.6 μ l (5 μ M) forward or reverse primer, and H 2 O with 10 μ l total volume. The thermal profile for sequencing reactions comprised of an initial denaturation step at 96°C for 1 min, followed by 25 cycles at 96°C for 10 s, 50°C for 5 s and 60°C for 105 s. 28S rDNA and mitochondrial sequences were sequenced with an ABI 3500xL genetic analyzer (Applied Biosystems).

Phylogenetic analyses
Multiple sequence alignments and subsequent comparison of the newly acquired (Tables 1 and 2) and GenBank-derived Helicotylenchus spp. sequences were performed using ClustalW algorithm as implemented in the BioEdit program v. 7.2.5 (Hall, 1999). The final multiple-sequence alignments were 595 and 588 nucleotide sites long in the case of 28S rDNA and mtCOI data, respectively. For the 28S rDNA data set, two Rotylenchus and two Hoplolaimus species were selected as outgroup: R. uniformis (MG653537, acquired in this study), R. magnus (EU280789), H. seinhorsti (DQ328752) and H. galeatus (EU626787). Heterodera glycines (HM462017) and R. uniformis (MG663121-2; this study) sequences were used as outgroups in mtCOI data set. Substitution models were tested using the "FindModel," an online implementation of MODELTEST (Posada and Crandall, 1998). A GTR + I + G substitution model was chosen for both 28S rDNA and mtCOI DNA data sets. The Bayesian phylogenetic trees were constructed with the program MrBayes v. 3.1 (Ronquist and Huelsenbeck, 2003). In both cases, two independent runs were performed with four Markov chains per run. In case of the 28S rDNA, the program was run for 3,000,000 generations and in case of mtCOI for 1,000,000 generations, with a sampling frequency of 100 generations. Burn-in trees of 50,000 and 25,000 generations, respectively, were discarded. For the mtCOI data, the partition by codon position was set. The sampled trees from each run were combined in a single 50% majority-rule tree. The stabilization of the likelihood and parameters was checked with the program Tracer v. 1.6 (Rambaut et al., 2014). Note: Individuals from whom both 28S rDNA and the mtCOI sequences were obtained are marked in grey.

Results
Helicotylenchus canadensis (Waseem, 1961).  Table 3). The morphological characteristics of H. canadensis females (n = 5) are as follows: the body is spiral when relaxed. The lateral field exists with four lines. Head is conical, truncated, bearing four to five annuli, continuous with body contour. Pharyngeal glands overlap intestine on the ventral side. The median oesophagus bulb is oval. Spermatheca is empty. Phasmids' position is three to nine annuli anterior to anus. The tail is slightly curved dorsally, bearing 10 to 16 annuli, a bit ventrally flattened with rounded, striated terminus.
Remarks: The morphological characteristics of H. canadensis under study agree with the first description of the species (Waseem, 1961) as well as following ones given for populations from England (Yuen, 1964), New Zealand (Yeates and Wouts, 1992), Poland (Brzeski, 1998), and Hungary (Andrássy, 2007). Although the ranges of some parameters (V, pharyngeal length, body diameter, tail length, and phasmids' po-sition) slightly differed from those given for the paratypes, they corresponded well to the ranges given by Brzeski (1998) for H. candensis occurring in European regions of the temperate climate.
Locality and habitat: See Table 2.
( Fig. 1B and E; Table 4). The morphological features of the H. pseudorobustus females (n = 13) are as follows: the body is spiral when relaxed. The lateral field exists with four lines. The lip region is hemispherical, continuous with body contour, with four to five annuli. Stylet is present with anteriorly flattened knobs. Pharyngeal glands overlap intestine on ventral side; the median pharyngeal bulb is oval. Spermatheca is empty. Phasmids' position is seven to ten annules anterior to anus. The tail is curved dorsally, bearing 6 to 13 annuli, with small, but distinct rounded ventral projection. Remarks: The morphological features of the investigated H. pseudorobustus individuals were congruent with the descriptions of the topotypes (Sher, 1966;Fortuner, 1984). However, the range of most of the parameters was larger. H. pseudorobustus is a cosmopolitan species and there are many descriptions available of this species from different regions. Consequently, the ranges of many parameters extended, compared to the original ones. Most morphometric parameters of H. pseudorobustus analyzed in our study were consistent with those given by Yeates and Wouts (1992), Brzeski (1998), , and Shokoohi et al. (2018). However, lower ranges of c′, tail length and anal body diameter were slightly smaller (Table 4). This finding confirms the claims of Fortuner et al. (1984) and Subbotin et al. (2015) concerning the large morphological variability of H. pseudorobustus.

Phylogenetic studies on three Helicotylenchus species based on 28S rDNA and mtCOI sequence data
Helicotylenchus varicaudatus (Yuen, 1964). (Fig. 1C and F; Fig. 2A to C; Table 5). The morphological features of the H. varicaudatus females are as follows: (n = 5). Habitus is loose spiral, the lateral field exists with four lines, head is conical, with four to five annuli, continuous with body contour. The length of stylet is 29 to 31.3 μ m, esophageal glands overlap intestine on ventral side. Spermatheca is empty. The tail of specimens is dorsally convex with a nearly straight ventral contour. The connection of dorsal contour with ventral one is almost at right angle. The tail is irregular at ventral terminus, with four to seven annuli and phasmids are located near anus (at anus level up to five annuli anterior to anus).
Remarks: Two H. varicaudatus populations were analyzed: one with and one without males (from Nowy Duninów/ Stary Duninów and Turew, respectively). Although, morphologically, our individuals corresponded well with the first description of H. varicaudatus made by Yuen (1964), ranges of most of the parameters differed from the original ones. H. varicaudatus was recorded in many countries and more descriptions of this species are available (Yeates and Wouts, 1992;Brzeski, 1998;Schreck Reis et al., 2010). Most parameters obtained in the current study were congruent with these descriptions. Nevertheless, tails of some females analyzed in our study were a little longer than tails of H. varicaudatus females described by other authors. It substantiates the claims of many authors (Yuen, 1964;Loof (1984); Brzeski, 1998) on the large variability of tail shape in this species and relates to the origin of the species name. Strikingly, female spermathecae of the bisexual population analyzed in our study were empty. Analyzing Dutch, bisexual population of H. varicaudatus, Loof (1984) also observed that some females had empty spermatheca. Of 86 analyzed females, only 47 had spermatheca filled with sperm. The morphological features of the H. varicaudatus males (n = 4) are as follows: arcuate ventrally, in many respects (number of lines in lateral field, ventral overlapping of oesophageal glands, shape of head, number of head annuli) similar to females. The body is usually shorter than females (612.8 to 715.2 vs 623.8 to 876.55), with shorter stylet than females (25.2 to 27.1 vs 29.0 to 31.3). Bursa is extended to the end of tail. Spicule is slightly arcuate 26.5 to 30.7 µm, gubernaculum 7.9 to 10.2 µm.
Remarks: Although females of H. varicaudatus have been recorded in many countries, males of this species occur extremely rarely. They were noted only in some locations in the Netherlands (Loof, 1984;Bongers, 1994), in Poland (Brzeski, 1998 and in coastal sand dunes in Portugal (Schreck Reis et al., 2010). Therefore, the descriptions of H. varicaudatus males are very scarce. In our study, we found only four males of H. varicaudatus in the rhizosphere of ash-leaved maple. Although these specimens corresponded to the description of Loof (1984) in respect of morphology, they differed in some morphometric parameters in relation to the Portuguese population. The males analyzed in our study had longer spicules than in the Portuguese population, 26.5 to 30.7 µm vs 20 to 25 µm. However, the spicule length of H. varicaudatus males herein did not significantly differ from the spicule length (23 to 28 µm) of males from Poland, as previously described by Brzeski (1998). The gubernaculum of H. varicaudatus males obtained in this study was much longer than those of males from Portugal as well as those from Poland (7.9 to 10.2 µm vs 5 to 7 µm).
Locality and habitat: See Table 2.

Sequence analysis
All sequences reported in this study have been deposited in GenBank and their accession numbers are listed in Table 2. By the use of the newly designed primers, we obtained mitochondrial sequences from 4 representatives of H. canadensis, 11 of H. pseudorobustus and 8 of H. varicaudatus ( Table 2). The particular amplicon and a derived size of the newly obtained sequences varied depending on the primer combination used (Table 1). The designed primers, in most combinations, targeted the region relatively closer to the 5′ end of the mtCOI gene sequence. The longest PCR amplicons, which were approximately 700 bp, were obtained using the M3.5 F and M8aR primer combination, which also  amplified a part of the mitochondrial sequence closer to the 3′ end of the mtCOI gene. However, the most successful primer combination was the M2F and RH1R one, which allowed for the amplification of PCR product of approximately 600 bp from all analyzed Helicotylenchus species. There was no intra-specific nucleotide variation concerning the four mtCOI sequences derived from H. canadensis. The potential minor nucleotide differences concerned few degenerate nucleotides, namely, K (G or T), Y (C or T), and R (A or G). The highest amplification success was obtained for the H. pseudorobustus and H. varicaudatus. By the use of three of the four presented primer combinations, we were able to obtain amplicons and subsequent mitochondrial sequences from 11 and 8 nematode individuals, respectively. The variability of the sequence amplified from H. pseudorobustus involved only some minor ambiguous nucleotides, possibly due to sequencing error. The only significant nucleotide substitution (A/T) was localized in the middle of the sequence at site 418 (MG663108), which had no influence on the potential amino acid sequence translation. Importantly, the more insightful analysis Holotype, Rothamsted, England, acc. (Yuen, 1964) Paratypes, Rothamsted, England, acc. (Yuen, 1964) Populations from New Zealand acc. (Yeates and Wouts, 1992) Populations from temperate Europe acc. (Brzeski, 1998)

17-20
Anal body diam. 14.8 ± 2.8 of the H. pseudorobustus mitochondrial nucleotide sequences revealed that all of them contain multiple TAA and TAG codons, which in the regular invertebrate mitochondrial genetic code would be translated as termination codons. A comparison of these sequences with the publicly available sequences derived from other Hoplolaimidae species (Fig. 3) suggests that Tyrosine (Y) and Phenylalanine (F) codons are the most vulnerable to such substitutions in this family, and that H. oleae-derived sequences also contain TAA codons van den Berg et al., 2016;Palomares-Rius et al., 2018). Both TAA and TAG triplets were also detected in sequences derived from two Scutellonema: S. brachyurus and S. truncatum, reported from study of van den Berg et al. (2013Berg et al. ( , 2017. No amino acid sequence ambiguities were predicted for the two other Helicotylenchus species analyzed here. Furthermore, no intra-specific variation on the nucleotide level was observed in the case of mtCOI sequences obtained from H. varicaudatus individuals. 28S rDNA fragments were successfully acquired from four H. canadensis, nine H. pseudorobustus, and four H. varicaudatus individuals ( Table 2). All fragments were approximately 550 nucleotides long and they were characterized by a lack of intra-specific sequence variation.
The highest values of inter-specific variations were observed between H. pseudorobustus and two other species: 32 to 39% in case of mitochondrial sequences and 12 to 13% in case of the 28S rDNA. The lowest inter-specific variation was between H. canadensis and H. varicaudatus, which was on the level of 16 to 20% in the case of mitochondrial sequences and 5% in the case of 28S rDNA.

28S rDNA and mtCOI-based phylogenetic relationships of the investigated Hoplolaimidae species
Bayesian phylogenetic analysis of the 28S rDNA sequences derived from various Helicotylenchus species resulted in a tree topology (Fig. 4) that agrees with the phylogenetic data presented by Subbotin et al. (2015), Palomares-Rius et al. (2018), and Shokoohi et al. (2018). In total, 12 Helicotylenchus clades (I to XI, as defined in Subbotin et al. 2015;number XII was   Congruent with results of 28S rRNA analyses above, the topology of the phylogenetic tree inferred from mtCOI showed that H. canadensis and H. varicaudatus share the closest pairwise relationship of the three Helicotylenchus species analyzed (Fig. 5) and they were inferred to be in Clade IX of the 28S rDNAbased tree. The species H. pseudorobustus (representative of Clade I) represented a separate lineage, whose most recent common ancestor with the former two species was that of all Helicotylenchus species. However, the relationships between the three distant lineages (Clades I, IX and XII) have not been resolved on the basis of the available mitochondrial data, most probably due to insufficient number of analyzed taxa.

Discussion
Morphological and molecular data in this study included three Helicotylenchus species: H. canadensis, H. pseudorobustus and H. varicaudatus, extracted from various locations and vegetation types in Poland ( Table 2). The morphological evaluation of the H. pseudorobustus nematode individuals and their positioning in the Clade I of the Helicotylenchus 28S rDNA tree agreed with the type characterization of this species, as described in Subbotin et al. (2015). Newly generated 28S rDNA sequences from H. canadensis and H. varicaudatus that are positioned in Clade IX are reported for the first time for those species. Morphological and morphometrical studies on these species were congruent with previous descriptions. The inferred, close phylogenetic relationship of H. canadensis to Helicotylenchus vulgaris (Yuen, 1964) supports morphological comparisons of these two species as being very similar. The most clearly observable differences between these two species were the length of the tail and the width of the annuli at the tail end: the tail of H. vulgaris is shorter than the tail of H. canadensis, and the middle annuli at the end of the H. vulgaris tail are narrower than other annules, whereas the annuli of the H. canadensis tail are of equal width. The mtCOI-based phylogenetic analysis included four newly sequenced Helicotylenchus spp. (Fig. 5) (van Megen et al., 2009;Schreck Reis et al., 2010). The 28S rDNA region is often more variable than the small subunit ribosomal DNA (Hillis and Dixon, 1991;Holterman, 2008) and it has been Figure 4: Phylogeny of the genus Helicotylenchus, as inferred by Bayesian analysis of 28S rDNA. The numbers near nodes indicate posterior probabilities. Roman numerals indicate major clades following Subbotin et al. (2011Subbotin et al. ( , 2015. The original 28S rDNA sequences are in boldface font.  Figure 5: Phylogeny of the genus Helicotylenchus, as inferred by Baysian analysis of partial mtCOI sequences. Roman numerals indicate major clades according to Subbotin et al. (2011Subbotin et al. ( , 2015. The original mtCOI sequences are in boldface font.  Subbotin et al., 2011Subbotin et al., , 2015Palomares-Rius et al., 2018). Accordingly, the majority of the Helicotylenchus sequences deposited in GenBank comprise 28S rDNA data. Although Gen-Bank already included 18S rDNA and 28S rDNA sequences for H. pseudorobustus, this study provides the first mtCOI fragments for this species. Characterized by a very high sequence variability, the mitochondrial data of Helicotylenchus species as well as other nematode genera are slowly increasing, being particularly useful in supporting the process of delineation of the most dubiously designated species. The mitochondrial DNA region encoding for mtCOI gene has been widely explored for barcoding application purposes in many animal organisms (Rach et al., 2017). It was also used in Hoplolaimidae and other nematode groups (Kanzaki and Futai, 2002;Derycke et al., 2010;Riepsamen et al., 2011;Gutiérrez-Gutiérrez et al., 2012;Kumari and Subbotin, 2012;van den Berg et al., 2013;Holguin et al., 2015;van den Berg et al., 2017;Sanchez-Monge et al., 2017;Skwiercz et al., 2017). In this study, to enhance the process of delineation and species identification within genus Helicotylenchus, which was underrepresented in terms of available mtCOI sequences, a new set of mtCOI primers was developed and tested. These newly designed primers offer an alternative to the most frequently used mtCOI primers combinations, especially in a situation when the original ones fail in the amplification of a desired PCR product. Most of the primer combinations presented in this publication allowed for the amplification of the fragments closer to the 5′ end of the assumed cytochrome oxidase c nucleotide sequence. The primer M2F in combination with the RH1R primer was the most successful in our study and allowed for amplification of a mtCOI sequence (600 bp) from all Helicotylenchus species analyzed here. The longest mtCOI sequence, including more nucleotides at its 3′ end, was acquired using the M3.5F and M8aR primer combination on H. pseudorbustus (and R. uniformis that was used as an outgroup species).
The inspection of mitochondrial sequences obtained here shows that all 11 mtCOI sequences from the nematode representatives of H. pseudorobustus were almost identical and included TAA and TAG codons at the same implied amino acid positions. Strikingly, such sequences were obtained each time for the H. pseudorobustus representatives, using three different mtCOI-targeted primer combinations (Table  1), whereas none of the sequences acquired from the other two Helicotylenchus species analyzed here revealed such nucleotide substitutions. In contrast, the homologous triplets in the mtCOI sequences of H. candadensis and H. varicaudatus instead contain triplets translated into tyrosine (TAT) and the phenylalanine (TTT) according to the standard invertebrate mitochondrial genetic code. This finding suggests alterations to the mitochondrial genetic code in H. pseudorobustus, if sequences are indeed translated through putative stop codons in this species.
Genetic code alterations in both nuclear and mitochondrial genomes have been reported in various organisms (Bessho et al., 1992;Telford et al. 2000;Lavrov et al., 2012;Lavrov, 2014;Bezerra et al., 2015;Pánek et al., 2017). In their studies, Jacob et al. (2009) showed that such alterations can also be found in nematodes. In the mitochondrial genome of nematode Radopholus similis (Pratylenchidae), the UAA (TAA) codon is reassigned from translation termination to tyrosine. Consequently, the usage of the TAA codon can be observed in mtCOI sequences acquired from various Hoplolaimidae species, including