identifier	taxonID	type	CVterm	format	language	title	description	additionalInformationURL	UsageTerms	rights	Owner	contributor	creator	bibliographicCitation
03A04416CC6E18507AEC74CDC256FDD2.text	03A04416CC6E18507AEC74CDC256FDD2.taxon	http://purl.org/dc/dcmitype/Text	http://rs.tdwg.org/ontology/voc/SPMInfoItems#GeneralDescription	text/html	en	Naididae Ehrenberg 1831	<div><p>3.1. Mitogenome characterization of Naididae</p><p>The P. hammoniensis mitogenome assembly was 14,881 bp long and was completely annotated. The partial mitogenomes of St. lacustris, C. diaphanus, C. diaphanus sp. B, Sl. appendiculata, and St. fossularis were 14,813; 15,724; 15,694; 15,699, and 15,719 bp long, respectively. All genes are dispersed and transcribed in the heavy chain, like other annelids (Oceguera-Figueroa et al., 2016; Weigert et al., 2016) (Figure 1; Table S2). The mitogenome sequences of P. hammoniensis, St. fossularis, St. lacustris, C. diaphanus, C. diaphanus sp. B, and Sl. appendiculata can be accessed from NCBI under accession numbers OQ654101, PP909790, PP909791, PP909792, PP909793 and PP893275, respectively.</p><p>As seen in other invertebrates, these mitogenomes showed an A+T bias of 67.8– 70.9%). Several adjacent genes overlapped from two (trnY-trnG, trnM-rrnS, trnAtrnS2 for all species) to 29 (trnQ- nd6, nd2-cox1 for C. diaphanus) nucleotides (Table S2). In the mitogenomes, the overall lengths of PCGs, tRNAs, and rRNAs were 11928, 1389, and 1982 bp, respectively. Conserved atp6-atp8 overlap was missing in this mitogenome, and accelerated evolution seems to lead to both differentiation of gene order and accumulation of intergenic regions. The conserved transcriptional sequence WHWGHTW identifiedinEcdysozoawasdetectedinthe P. hammoniensis mitogenome, 14 bp upstream of the trnT gene (Aydemir et al., 2023).</p><p>All PCGs were found to start with the standard ATG (methionine) codon, except for nd6 (ATT) in P. hammoniensis; nd1 (GTG) in St. lacustris; nd2 (CTT), nd1 (GTG), and nd3 (ATT) in St. fossularis; nd2 (TTT), nd4L (ATA), and nd1 (GTG) in Sl. appendiculata; nd3 (ATA) and nd1 (GTG) in C. diaphanus and nd2 (ATC), nd1 (GTG), nd3 (ATA), and atp8 (TTA) in C. diaphanus sp. B. While cox1, cox3, cytB, atp6, nd5, nd4L, and nd4 genes had complete TAA triplets as stop codons, cox2, atp8, nd6, nd1, nd3, and nd2 genes had incomplete T-stop codons that are completed via polyadenylation processes. Like the total mitogenome, all PCGs had negative AT- and negative GC-skew values. From the ORF Finder analysis, the P. hammoniensis mitogenome is capable of coding 780 different ORFs that start with any sense codon and are longer than 30 nucleotides (Supplementary file S1). Interestingly, 445 of these 780 ORFs are potentially encoded in the light chain, which does not code any known mitochondrial gene. Based on BLAST searches, 39 of 780 potential ORFs aligned with protein sequences in the nonredundant protein database (Supplementary file S2). Thirteen of the 39 were formal mitochondrial PCGs (atp6, atp8, cox1, cox2, cox3, cytB, nd1, nd2, nd3, nd4, nd4L, nd5, nd6) and five of the 39 were different framed forms of cox1 and cox2 sequences. Eighteen of the 39 aligned amino acid sequences were potentially coded in the light chain, eight of which seemed to be an alternative form of some nd genes (nd1, nd3, nd4 and nd5) in different organisms (Supplementary file S2). Similarly, 726, 727, 757, 770 and 736 potential ORF sequences were detected in the St. fossularis, St. lacustris, C. diaphanus, C. diaphanus sp. B, and Sl. appendiculata mitogenomes, respectively.</p><p>The average length of tRNAs was 63.13 bp and varied from 56 bp in trnT to 69 bp in trnQ (Table S2). All tRNAs constructed had a conserved clover-leaf structure except trnS1 and trnT, both of which had a DHU-replacement (Figure S1). A reduced secondary structure pattern is common in invertebrates for trnS1, but trnT reduction is unique to P. hammoniensis in Clitellata . trnE, trnP, trnS1, and trnS2 tRNAs had negative AT- and negative GC- skew values. Unlike the total mitogenome and PCGs, tRNA genes tended to have at least one positive skewness value (Table S2).</p><p>The rrnS and rrnL genes were 737 and 1245 bp long, respectively, and they had positive AT- and negative GCskew values (Table S2). Both rRNA genes had conserved putative secondary structures (Figures 2a and 2b). The putative secondary structures of rrnS and rrnL genes comprised three domains with 30 helices and five domains (domain III is absent in invertebrates) with 49 helices, respectively.</p><p>Noncoding regions of the mitogenome of P. hammoniensis were dispersed between 11 consecutive gene clusters and had a total length of 472 bp, varying from 2 (trnH-trnR and rrnS-trnV) to 208 (trnR- nd5) bp long. Intergenic regions did not exhibit significant matches in the BLAST search with low query coverage values (e&gt; 0.01). In Clitellata, the putative control region is generally located in the atp6-nd5 gene cluster (trnW- atp6 -trnHtrnR- nd5 -trnF), and the P. hammoniensis mitogenome had 39-bp-long total intergenic sequences in this region (Table S2). In the nd5 -trnF intergenic region, transcription polarity was dramatically changed. This region probably has the replication origin that can construct a stable secondary structure with significant ΔG value (−13.60 kcal/mol) (Figure S2).</p></div>	https://treatment.plazi.org/id/03A04416CC6E18507AEC74CDC256FDD2	Public Domain	No known copyright restrictions apply. See Agosti, D., Egloff, W., 2009. Taxonomic information exchange and copyright: the Plazi approach. BMC Research Notes 2009, 2:53 for further explanation.		Plazi	Mercan, Deniz;Horenkamp, Austin;Aydemir, Merve Nur;Aydemir, Habeş Bilal;Neubig, Kurt;Arslan, Naime;Anderson, Frank E.	Mercan, Deniz, Horenkamp, Austin, Aydemir, Merve Nur, Aydemir, Habeş Bilal, Neubig, Kurt, Arslan, Naime, Anderson, Frank E. (2025): Comparative mitogenomics of Clitellata reveals mitogenome organization can affect the mode of phylogeny. Turkish Journal of Zoology 49 (4): 197-206, DOI: 10.55730/1300-0179.3225, URL: https://doi.org/10.55730/1300-0179.3225
03A04416CC6C185E7849776CC07DFA04.text	03A04416CC6C185E7849776CC07DFA04.taxon	http://purl.org/dc/dcmitype/Text	http://rs.tdwg.org/ontology/voc/SPMInfoItems#GeneralDescription	text/html	en	Clitellata	<div><p>3.2. Comparative analysis of Clitellata mitogenomes</p><p>The size of complete Clitellata mitogenomes varied from 14,407 bp in Erpobdella octoculata (Hirudinea: Erpobdellidae) to 18,528 bp in Acanthobdella peledina (Hirudinea: Acanthobdellidae), with an average of 15,096 bp in length. Total mitogenomes of Clitellata species were biased to A+T and generally have negative AT- and GCskew values. Exceptionally, 13 species from Hirudinea and two species from Oligochaeta had moderately positive GC-skew values, with most of Oligochaeta (32 of 46) mildly biased to positive AT- skewness (Table S3). Average nucleotide compositions were as follows: T = 35.3%, A = 33.4%, C = 17.5%, and G = 13.8%.</p><p>Total PCG size varied from 10,979 bp ( Haemadipsa tianmushana) to 11,244 bp ( Theromyzon tessulatum), with an average of 11,102.1 bp long. Like the total mitogenome data, total PCGs of Clitellata species were biased to A+T and generally had negative AT- and GC-skew values.</p><p>Amino acid preferences of Clitellata were represented mostly by hydrophobic amino acids as expected with an A+T nucleotide bias: leucine (14.75%), serine (9.80%), isoleucine (8.86%), and methionine (8.37%). Cysteine (1.04%), arginine (1.60%), and glutamine (1.74%) amino acid preferences had the lowest levels (Table S4).</p><p>RSCU values from whole sampling are given in</p><p>Figure 3. All stop codons were removed to avoid the bias of incomplete stop codons. All 4-fold degenerate codons preferred A in the third codon position. All 2-fold degenerate codons were biased to keep A or T nucleotides consistent with the A+T bias of the total mitogenome. On the contrary, anticodons of complementary tRNAs are conserved among Clitellata and are biased to G and T nucleotides in the first position. Anticodons of 9 tRNAs start with G or C nucleotides and therefore some interactions do not exhibit exact matches between the tRNA anticodon and the mRNA codon, as previously documented in the trnP interaction in hymenopteran insects (Aydemir and Korkmaz, 2020). The situation in Clitellata appears more dramatic. This pattern is explained by the fact that the transcription of mitogenomes can be partially restricted by mitochondrial tRNAs (Shtolz and Mishmar, 2023). On the other hand, this is one reason why the most conserved interaction position is the second codon (anticodon) position (Błażej et al., 2018).</p><p>The average rrnS gene length is 763.4 bp and ranges from 688 ( Hirudinaria manillensis and Poecilobdella manillensis) to 801 bp ( Perionyx excavatus). Average nucleotide compositions of rrnS are as follows: T = 29.7%, A = 38.5%, C = 16.1%, and G = 15.7%. The rrnS gene has a positive AT-skew (except for Poecilobdella javanica) but variable GC-skew (32 of 78 have positive values, and 45 of 78 have negative values) (Table S5). The size of the rrnL gene varied from 1027 bp in Enchytraeus irregularis to 12,812 bp in Amynthas deogyusanensis, with an average of 1206.4 bp in length. Average nucleotide compositions of rrnL are calculated as T = 31.0%, A = 39.8%, C = 15.1%, and G = 14.1% (Table S6). As seen with the rrnS gene, the rrnL gene of Clitellata generally has positive AT- and variable GC-skew values.</p><p>The length of tRNA genes varied from 41 (trnR of W. pigra) to 83 (trnS of some Hirudinea organisms) with an average of 63.4 bp in length (Table S7). Average nucleotide compositions of tRNAs are calculated as follows: T = 33.3%, A = 35.9%, C = 14.7%, and G = 16.1%.</p><p>The results of dN/dS calculations indicated that all mitochondrial PCGs are subjected to negative selection (dN/dS &lt;0.5) in Clitellata (Figure 4). Mutations that may occur in the cox1 gene are strongly and rapidly purified and following that, the cox3, cytB, nd1 and cox2 genes are subjected to densified negative selection. It is expected for cox and cytB genes because of their critical roles in OXPHOS (Castellana et al., 2011), but purifying selection on the nd1 gene is unexpected, especially despite the relaxation of negative selection observed in other nd subunits. The most likely reason may be that the nd1 gene product functions as a binding site for the hydrophobic quinone during OXPHOS. This is essential as it interfaces between hydrophobic and hydrophilic subunits, and it is crucial to eliminate mutations in the nd1 gene to stabilize this interaction (Garvin et al., 2015). On the other hand, the atp8, nd6, nd2, and nd4L genes had the highest dN/dS values and seem to have evolved relatively rapidly among mitochondrial PCGs (Figure 4).</p><p>From the CREx analysis, nine different gene order patterns were observed within Clitellata (Figure 5). A minimum of three and a maximum of 20 rearrangement breakpoints were detected between the reference mitogenome and the clitellate mitogenomes (Table S8). Accordingly, the transposition and tandem-duplicationrandom-loss (TDRL) event of trnG-trnY- atp8, the TDRL event of trnM-trnD gene clusters and the transposition of trnC were characterised. The most divergent gene order was determined in Olavius algarvensis and generally four different patterns were determined in Naididae . Except for Naididae (with five representatives) and Lumbriculidae ( Lumbriculus variegatus), the oligochaete species shared a common mitochondrial gene order. Within Hirudinea, two unique gene orders were generally detected in two main groups (Erpobdelliformes + Oceanobdelliformes and Hirudiniformes) with a few exceptions. Interestingly, the mitochondrial gene order of P. hammoniensis was identical to Theromyzon tessulatum (Hirudinea: Glossiphoniidae) (Figure 5).</p><p>3.3. Phylogenetic analysis of Clitellata</p><p>In the phylogenetic analysis of the newly generated mitogenomes combined with all available clitellate mitogenomes on GenBank, St. lacustris, St. fossularis, C. diaphanus, C. diaphanus sp. B, and Sl. appendiculata clustered with the other naidines with high support (Figure 5). These findings are congruent with previous studies, with Naidinae forming a sister clade to Tubificinae, with near maximum support on every branch. In previous studies, within the species used in this study, Naididae is recovered as a sister clade to most other clitellate taxa, but the mitogenome tree recovers Enchytraeidae (represented by E. irregularis) in that position (Erséus et al., 2020). Mitogenomes from representatives of Phreodrilidae and Propappidae could aid in recovering a similar tree to the phylogenomic tree in Erseus et al., 2020. In this analysis, the leeches are recovered as the sister taxon to Lumbriculidae, which matches a previous phylogenomic study (Erséus et al., 2020). The unusual relationships recovered in our mitogenomic tree may reflect phylogenetic reality, but they could also be a result of long-branch attraction between Hirudinida and Naididae, which itself may be due in part to the relatively poor sampling of mitogenomes from Clitellata to date—there are no published mitogenomes available for several major lineages represented in Erséus et al. (2020). Levels of bootstrap support for many taxa, particularly within Megascolecidae, are low, so relationships between members of this family and other families are less trustworthy.</p><p>This is the most complete mitogenome study using all available clitellate mitogenomes to date, with six newly created mitogenomes for a total of 90 species. Previous studies of clitellate evolution examined the evolution of carnivory within the group, and the timing and nature of transitions among different habitat types (i.e. marine, freshwater, and terrestrial) (Erséus et al., 2020; Horn et al., 2019; Mack et al., 2023; Rousset et al., 2008). The mitogenome data here could aid future studies of clitellate phylogeny and evolution, but accurate phylogenetic inferences will likely require sampling of mitogenomes from additional clitellate lineages, including Branchiobdellidae, Capilloventridae, Haplotaxidae, Parvidrillidae, Phreodrillidae and Propappidae . To conduct a more thorough mitogenome study, additional specimens from the following clades should also be obtained: Alluroididae, Capilloventridae, Crassiclitellata, Haplotaxidae, Moniligastridae, Narapidae, Parvidrillidae, and Randiellidae, as well as additional members of Naididae .</p></div>	https://treatment.plazi.org/id/03A04416CC6C185E7849776CC07DFA04	Public Domain	No known copyright restrictions apply. See Agosti, D., Egloff, W., 2009. Taxonomic information exchange and copyright: the Plazi approach. BMC Research Notes 2009, 2:53 for further explanation.		Plazi	Mercan, Deniz;Horenkamp, Austin;Aydemir, Merve Nur;Aydemir, Habeş Bilal;Neubig, Kurt;Arslan, Naime;Anderson, Frank E.	Mercan, Deniz, Horenkamp, Austin, Aydemir, Merve Nur, Aydemir, Habeş Bilal, Neubig, Kurt, Arslan, Naime, Anderson, Frank E. (2025): Comparative mitogenomics of Clitellata reveals mitogenome organization can affect the mode of phylogeny. Turkish Journal of Zoology 49 (4): 197-206, DOI: 10.55730/1300-0179.3225, URL: https://doi.org/10.55730/1300-0179.3225
