The identification of members of the Onchidiidae is based on morphological characters; this is often time-consuming and can be inconclusive. In order to explore the species diversity of onchidiids in China, we provide a phylogeny constructed using partial sequences of two mitochondrial genes (16S rRNA and COI) and one nuclear ribosomal RNA gene (28S rRNA) from 32 samples comprising five genera. The topology, using both Bayesian and Maximum Likelihood inference methods, showed that the taxa clustered in two main groups of six species, one of which included Platevindex mortoni, Platevindex sp. and Onchidium ‘struma’; the other included Paraoncidium reevesii, Onchidella sp. and Peronia verruculata. It is clear that COI will be useful in discriminating onchidiid species-group taxa.

The Onchidiidae family is ideal for studying the evolution of marine invertebrate species from sea to wetland environments. However, comparative studies of Onchidiidae species are rare. A total of 40 samples were collected from four species (10 specimens per onchidiid), and their histological and molecular differences were systematically evaluated to elucidate the morphological foundations underlying the adaptations of these species. A histological analysis was performed to compare the structures of respiratory organs (gill, lung sac, dorsal skin) among onchidiids, and transcriptome sequencing of four representative onchidiids was performed to investigate the molecular mechanisms associated with their respective habitats. Twenty-six SNP markers of Onchidium reevesii revealed some DNA polymorphisms determining visible traits. Non-muscle myosin heavy chain II (NMHC II) and myosin heavy chain (MyHC), which play essential roles in amphibian developmental processes, were found to be differentially expressed in different onchidiids and tissues. The species with higher terrestrial ability and increased integrated expression of Os-MHC (NMHC II gene) and the MyHC gene, illustrating that the expression levels of these genes were associated with the evolutionary degree. This study provides a comprehensive analysis of the adaptions of a diverse and widespread group of invertebrates, the Onchidiidae. Some onchidiids can breathe well through gills and skin when under seawater, and some can breathe well through lung sacs and skin when in wetlands. A histological comparison of respiratory organs and the relative expression levels of two genes provided insights into the adaptions of onchidiids that allowed their transition from shallow seas to wetlands. This work provides a valuable reference and might encourage further study.

The mitochondrial genome (mitogenome) can provide information for phylogenetic analyses and evolutionary biology. We first sequenced, annotated, and characterized the mitogenome of Philomycus bilineatus in this study. The complete mitogenome was 14,347 bp in length, containing 13 protein-coding genes (PCGs), 23 transfer RNA genes, two ribosomal RNA genes, and two non-coding regions (A + T-rich region). There were 15 overlap locations and 18 intergenic spacer regions found throughout the mitogenome of P. bilineatus. The A + T content in the mitogenome was 72.11%. All PCGs used a standard ATN as a start codon, with the exception of cytochrome c oxidase 1 (cox1) and ATP synthase F0 subunit 8 (atp8) with TTG and GTG. Additionally, TAA or TAG was identified as the typical stop codon. All transfer RNA (tRNA) genes had a typical clover-leaf structure, except for trnS1 (AGC), trnS2 (TCA), and trnK (TTT). A phylogenetic analysis with another 37 species of gastropods was performed using Bayesian inference, based on the amino acid sequences of 13 mitochondrial PCGs. The results indicated that P. bilineatus shares a close ancestry with Meghimatium bilineatum. It seems more appropriate to reclassify it as Arionoidea rather than Limacoidea, as previously thought. Our research may provide a new meaningful insight into the evolution of P. bilineatus.