Ten new complete mitochondrial genomes of pulmonates (Mollusca: Gastropoda) and their impact on phylogenetic relationships

Complementing aculiferan mitogenomics: comparative characterization of mitochondrial genomes of Solenogastres (Mollusca, Aplacophora)

BACKGROUND

With the advances in high-throughput sequencing and bioinformatic pipelines, mitochondrial genomes have become increasingly popular for phylogenetic analyses across different clades of invertebrates. Despite the vast rise in available mitogenomic datasets of molluscs, one class of aplacophoran molluscs - Solenogastres (or Neomeniomorpha) - is still neglected.

RESULTS

Here, we present six new mitochondrial genomes from five families of Solenogastres (Amphimeniidae, Gymnomeniidae, Proneomeniidae, Pruvotinidae, Simrothiellidae), including the first complete mitogenomes, thereby now representing three of the four traditional orders. Solenogaster mitogenomes are variable in size (ranging from approximately 15,000 bp to over 17,000 bp). The gene order of the 13 protein coding genes and two rRNA genes is conserved in three blocks, but considerable variation occurs in the order of the 22 tRNA genes. Based on phylogenetic analyses and reconstruction of ancestral mitochondrial genomes of Aculifera, the position of (1) trnD gene between atp8 and atp6, (2) trnT and P genes between atp6 and nad5, and (3) trnL1 gene between G and E, resulting in a 'MCYWQGL1E'-block of tRNA genes, are all three considered synapomorphies for Solenogastres. The tRNA gene block 'KARNI' present in Polyplacophora and several conchiferan taxa is dissolved in Solenogastres.

CONCLUSION

Our study shows that mitogenomes are suitable to resolve the phylogenetic relationships among Aculifera and within Solenogastres, thus presenting a cost and time efficient compromise to approach evolutionary history in these clades.

The Characterization of the Mitochondrial Genome of Fulgoraria rupestris and Phylogenetic Considerations within the Neogastropoda

Fulgoraria rupestris is a predatory marine gastropod belonging to Neogastropoda and possessing considerable taxonomic significance. However, research on this species remains limited. We acquired the complete mitochondrial genome of F. rupestris through second-generation sequencing and conducted an analysis of its genome structural features. The mitochondrial genome of F. rupestris spans a total length of 16,223 bp and encompasses 37 genes (13 protein-coding genes (PCGs), 22 transfer RNAs, and 2 ribosomal RNAs). Notably, most tRNAs exhibit the typical cloverleaf structure, but there is an absence of the Dihydrouridine (DHU) arm in the trnS1 and trnS2 genes. The A + T content is 68.67%, indicating a pronounced AT bias. Additionally, we conducted a selection pressure analysis on the mitochondrial genomes of four species within Volutidae, revealing that all PCGs are subjected to purifying selection. In comparison to other species within Neogastropoda, F. rupestris shares an identical gene arrangement. Additionally, based on mitochondrial genome sequences of the 13 PCGs from 50 species within Neogastropoda, we constructed a phylogenetic tree. The phylogenetic tree indicates F. rupestris forms a clade with species within the family Volutidae (Cymbium olla, Neptuneopsis gilchristi, and Melo melo). This study serves as a valuable reference for future research on F. rupestris, offering insights for the upcoming phylogenetic and taxonomic classification within Neogastropoda. Furthermore, the findings provide valuable information for the development of genetic resources in this context.

New mitogenomes of Runcinidae and Facelinidae: two understudied heterobranch families (Mollusca: Gastropoda)

Here, we present the mitochondrial sequences of two sea slugs (Heterobranchia): Runcina aurata and Facelina auriculata, the latter being the type species of the family. The mitochondrial genomes are 14,282 and 14,171bp in length, respectively, with a complete set of 13 PCGs, 2 rRNAs, and 22 tRNAs. None of the mitogenomes show gene reorganization, keeping the standard mitogenomic structure of Heterobranchia. Nucleotide composition differs significantly between them, with R. aurata showing the most AT-rich mitogenome (25.7% GC content) reported to date in Heterobranchia, and F. auriculata showing a rich GC content (35%) compared with other heterobranch mitochondrial genomes.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Mitogenomic phylogeny of Cypraeidae (Gastropoda: Mesogastropoda)

Cowries (Family Cypraeidae) are widely distributed in tropical and subtropical seas, with the highest diversity of cowries in the Indo-Pacific region. However, the classification of Cypraeidae, especially at the lower taxonomic levels, is still controversial. In the present study, we determined the complete mitochondrial genomes of 10 cowries. All the newly sequenced mtDNA encode 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and several non-coding regions. The reconstructed mitogenomic phylogeny of Cypraeidae recognized two monophyletic clades, with the first clade comprised of Erroneinae, Cypraeinae and Luriinae and the second clade formed by the single subfamily Erosarinae. The congeneric genetic distance values fall within 0.118–0.144, lower than those above genus level ranging from 0.163 to 0.271, consistent with the current division of genera within Cypraeidae. The divergence time estimated here indicated that the ancestor of Indo-Pacific cowries diversified during the Paleocene, and the closure of the Tethys Seaway might lead to the speciation events of several Indo-Pacific species. This study suggests that the complete mtDNA is a promising tool to improve the phylogenetic resolution of Cypraeidae, and mtDNA could also provide important information for future species delimitation especially within the cowries that possess different morphological phenotypes.

The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei

The enigmatic scaphopods, or tusk shells, are a small and rare group of molluscs whose phylogenomic position among the Conchifera is undetermined, and the taxonomy within this class also needs revision. Such work is hindered by there only being a very few mitochondrial genomes in this group that are currently available. Here, we present the assembly and annotation of the complete mitochondrial genome from Dentaliida Pictodentalium vernedei, whose mitochondrial genome is 14,519 bp in size, containing 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The nucleotide composition was skewed toward A-T, with a 71.91% proportion of AT content. Due to the mitogenome-based phylogenetic analysis, we defined P. vernedei as a sister to Graptacme eborea in Dentaliida. Although a few re-arrangements occurred, the mitochondrial gene order showed deep conservation within Dentaliida. Yet, such a gene order in Dentaliida largely diverges from Gadilida and other molluscan classes, suggesting that scaphopods have the highest degree of mitogenome arrangement compared to other molluscs.

The complete mitochondrial genome of the terrestrial snail Monachacartusiana (O.F. Müller, 1774) (Gastropoda, Eupulmonata, Hygromiidae)

The mitochondrial genome of Monachacartusiana is the first complete mitochondrial sequence described for the pulmonate snail genus Monacha and for the family Hygromiidae. The identified mitogenome has a length of 13,894 bp and encodes 13 proteins, 22 tRNAs, and two rRNAs. A phylogenetic analysis of available mitogenomes from representatives of helicoid families shows a sister group relationship of Hygromiidae and Geomitridae, which have been recently recognised as separate families.

Phylogenomic resolution of the root of Panpulmonata, a hyperdiverse radiation of gastropods: new insight into the evolution of air breathing

Transitions to terrestriality have been associated with major animal radiations including land snails and slugs in Stylommatophora (>20 000 described species), the most successful lineage of 'pulmonates' (a non-monophyletic assemblage of air-breathing gastropods). However, phylogenomic studies have failed to robustly resolve relationships among traditional pulmonates and affiliated marine lineages that comprise clade Panpulmonata (Mollusca, Gastropoda), especially two key taxa: Sacoglossa, a group including photosynthetic sea slugs, and Siphonarioidea, intertidal limpet-like snails with a non-contractile pneumostome (narrow opening to a vascularized pallial cavity). To clarify the evolutionary history of the panpulmonate radiation, we performed phylogenomic analyses on datasets of up to 1160 nuclear protein-coding genes for 110 gastropods, including 40 new transcriptomes for Sacoglossa and Siphonarioidea. All 18 analyses recovered Sacoglossa as the sister group to a clade we named Pneumopulmonata, within which Siphonarioidea was sister to the remaining lineages in most analyses. Comparative modelling indicated shifts to marginal habitat (estuarine, mangrove and intertidal zones) preceded and accelerated the evolution of a pneumostome, present in the pneumopulmonate ancestor along with a one-sided plicate gill. These findings highlight key intermediate stages in the evolution of air-breathing snails, supporting the hypothesis that adaptation to marginal zones played an important role in major sea-to-land transitions.

Global diversification and evolutionary history of onchidiid slugs (Gastropoda, Pulmonata)

Many marine species are specialized to specific parts of a habitat. In a mangrove forest, for instance, species may be restricted to the mud surface, the roots and trunks of mangrove trees, or rotting logs, which can be regarded as distinct microhabitats. Shifts to new microhabitats may be an important driver of sympatric speciation. However, the evolutionary history of these shifts is still poorly understood in most groups of marine organisms, because it requires a well-supported phylogeny with relatively complete taxon sampling. Onchidiid slugs are an ideal case study for the evolutionary history of habitat and microhabitat shifts because onchidiid species are specialized to different tidal zones and microhabitats in mangrove forests and rocky shores, and the taxonomy of the family in the Indo-West Pacific has been recently revised in a series of monographs. Here, DNA sequences for onchidiid species from the North and East Pacific, the Caribbean, and the Atlantic are used to reconstruct phylogenetic relationships among Onchidella species, and are combined with new data for Indo-West Pacific species to reconstruct a global phylogeny of the family. The phylogenetic relationships of onchidiid slugs are reconstructed based on three mitochondrial markers (COI, 12S, 16S) and three nuclear markers (28S, ITS2, H3) and nearly complete taxon sampling (all 13 genera and 62 of the 67 species). The highly-supported phylogeny presented here suggests that ancestral onchidiids most likely lived in the rocky intertidal, and that a lineage restricted to the tropical Indo-West Pacific colonized new habitats, including mudflats, mangrove forests, and high-elevation rainforests. Many onchidiid species in the Indo-West Pacific diverged during the Miocene, around the same time that a high diversity of mangrove plants appears in the fossil record, while divergence among Onchidella species occurred earlier, likely beginning in the Eocene. It is demonstrated that ecological specialization to microhabitats underlies the divergence between onchidiid genera, as well as the diversification through sympatric speciation in the genera Wallaconchis and Platevindex. The geographic distributions of onchidiid species also indicate that allopatric speciation played a key role in the diversification of several genera, especially Onchidella and Peronia. The evolutionary history of several morphological traits (penial gland, rectal gland, dorsal eyes, intestinal loops) is examined in relation to habitat and microhabitat evolutionary transitions and that the rectal gland of onchidiids is an adaptation to high intertidal and terrestrial habitats.

The mitochondrial genome of the semi-slug Omalonyx unguis (Gastropoda: Succineidae) and the phylogenetic relationships within Stylommatophora

Here we report the first complete mitochondrial genome of the semi-slug Omalonyx unguis (d’Orbigny, 1836) (Gastropoda: Succineidae). Sequencing was performed on a specimen from Argentina. Assembly was performed using Sanger data and Illumina next generation sequencing (NGS). The mitogenome was 13,984 bp in length and encoded the 37 typical Metazoan genes. A potential origin for mitochondrial DNA replication was found in a non-coding intergenic spacer (49 bp) located between cox3 and tRNA-Ile genes, and its secondary structure was characterized. Secondary structure models of the tRNA genes of O. unguis largely agreed with those proposed for other mollusks. Secondary structure models for the two rRNA genes were also obtained. To our knowledge, the 12S-rRNA model derived here is the first complete one available for mollusks. Phylogenetic analyses based on the mitogenomes of O. unguis and 37 other species of Stylommatophora were performed using amino acid sequences from the 13 protein-coding genes. Our results located Succineoidea as a sister group of Helicoidea + Urocoptoidea, similar to previous studies based on mitochondrial genomes. The gene arrangement of O. unguis was identical to that reported for another species of Succineoidea. The unique rearrangements observed for this group within Stylommatophora, may constitute synapomorphies for the superfamily.

Evolutionary genomics of endangered Hawaiian tree snails (Achatinellidae: Achatinellinae) for conservation of adaptive capacity

Phylogenomic studies can provide insights into speciation, adaptation, and extinction, while providing a roadmap for conservation. Hawaiian tree snails are a model system for an adaptive radiation facing an extinction crisis. In the last 5 years, nearly all populations of Hawaiian tree snails across the 30 remaining species in the subfamily Achatinellinae (Achatinellidae) have declined from hundreds or thousands in the wild down to undetectable levels. Nearly 100 species historically occurred across dramatic environmental gradients on five of the Hawaiian Islands, but habitat loss, overcollection, and predation by invasive species have decimated populations. As such, this system offers the opportunity to integrate efforts to conserve evolutionary potential into conservation planning for a rapidly declining subfamily. Here, we used genome-wide, restriction-site associated DNA sequencing (RADseq), along with mitochondrial genome reconstruction, to resolve evolutionary relationships to inform conservation efforts. Phylogenetic analysis of nearly 400k genome-wide SNPs from 59 populations and 25 species across six genera in the family Achatinellidae, was generally concordant with taxonomy, geography, and mtDNA with several notable exceptions; mtDNA was unable to resolve some deeper nodes (e.g., the monophyly of Achatinella), while SNP data did not resolve as many shallow nodes. Both phylogenetic and coalescent analysis revealed deep divergences between populations within Achatinella mustelina that were consistent with species-level differences. Given cryptic species-level divergence within populations that are geographically proximate, they are at higher risk of extirpation from invasive predators and climate change than previously assumed. This study clarifies evolutionary relationships within this model system for adaptive radiation, forming the basis for conservation strategies such as translocation, captive rearing, and hybridization trials to prevent the loss of capacity to adapt to rapidly changing environmental conditions.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Assessment of mitochondrial genomes for heterobranch gastropod phylogenetics

Background

Heterobranchia is a diverse clade of marine, freshwater, and terrestrial gastropod molluscs. It includes such disparate taxa as nudibranchs, sea hares, bubble snails, pulmonate land snails and slugs, and a number of (mostly small-bodied) poorly known snails and slugs collectively referred to as the “lower heterobranchs”. Evolutionary relationships within Heterobranchia have been challenging to resolve and the group has been subject to frequent and significant taxonomic revision. Mitochondrial (mt) genomes can be a useful molecular marker for phylogenetics but, to date, sequences have been available for only a relatively small subset of Heterobranchia.

Results

To assess the utility of mitochondrial genomes for resolving evolutionary relationships within this clade, eleven new mt genomes were sequenced including representatives of several groups of “lower heterobranchs”. Maximum likelihood analyses of concatenated matrices of the thirteen protein coding genes found weak support for most higher-level relationships even after several taxa with extremely high rates of evolution were excluded. Bayesian inference with the CAT + GTR model resulted in a reconstruction that is much more consistent with the current understanding of heterobranch phylogeny. Notably, this analysis recovered Valvatoidea and Orbitestelloidea in a polytomy with a clade including all other heterobranchs, highlighting these taxa as important to understanding early heterobranch evolution. Also, dramatic gene rearrangements were detected within and between multiple clades. However, a single gene order is conserved across the majority of heterobranch clades.

Conclusions

Analysis of mitochondrial genomes in a Bayesian framework with the site heterogeneous CAT + GTR model resulted in a topology largely consistent with the current understanding of heterobranch phylogeny. However, mitochondrial genomes appear to be too variable to serve as good phylogenetic markers for robustly resolving a number of deeper splits within this clade.

Genetic assessment of the variation and distribution of the species of Salinator (Panpulmonata: Amphibolidae) in south-eastern Australia

Amphibolidae is one of the most abundant families of gastropods in estuarine environments of south-eastern Australia. However, the range limits of the species of Salinator, one of the family’s two genera in the region, remain unclear partly due to uncertainty of identifications based solely on shell morphology. Insufficient data have been collected to address questions regarding the genetic variability of any of the species of Salinator. Here, DNA sequences from a segment of the cytochrome c oxidase subunit I and 28S ribosomal RNA genes were collected to investigate the distribution and variation of the three Salinator species in the region, these being S.fragilis, S.rhamphidia and S.tecta. The results demonstrate a large range extension in S.rhamphidia and suggest that S.tecta may have limited distribution in Tasmania. In contrast to previously-studied estuarine Mollusca in the south-eastern coasts of the mainland and Tasmania, S.rhamphidia has regional differentiation. There is evidence of genetic disequilibrium within S.fragilis, suggesting that it may presently comprise contributions from two distinct sets of populations.

Systematic revision of the genus Peronia Fleming, 1822 (Gastropoda, Euthyneura, Pulmonata, Onchidiidae)

The genus Peronia Fleming, 1822 includes all the onchidiid slugs with dorsal gills. Its taxonomy is revised for the first time based on a large collection of fresh material from the entire Indo-West Pacific, from South Africa to Hawaii. Nine species are supported by mitochondrial (COI and 16S) and nuclear (ITS2 and 28S) sequences as well as comparative anatomy. All types available were examined and the nomenclatural status of each existing name in the genus is addressed. Of 31 Peronia species-group names available, 27 are regarded as invalid (twenty-one synonyms, sixteen of which are new, five nomina dubia, and one homonym), and four as valid: Peronia peronii (Cuvier, 1804), Peronia verruculata (Cuvier, 1830), Peronia platei (Hoffmann, 1928), and Peronia madagascariensis (Labbé, 1934a). Five new species names are created: P. griffithsi Dayrat & Goulding, sp. nov., P. okinawensis Dayrat & Goulding, sp. nov., P. setoensis Dayrat & Goulding, sp. nov., P. sydneyensis Dayrat & Goulding, sp. nov., and P. willani Dayrat & Goulding, sp. nov. Peronia species are cryptic externally but can be distinguished using internal characters, with the exception of P. platei and P. setoensis. The anatomy of most species is described in detail here for the first time. All the secondary literature is commented on and historical specimens from museum collections were also examined to better establish species distributions. The genus Peronia includes two species that are widespread across the Indo-West Pacific (P. verruculata and P. peronii) as well as endemic species: P. okinawensis and P. setoensis are endemic to Japan, and P. willani is endemic to Northern Territory, Australia. Many new geographical records are provided, as well as a key to the species using morphological traits.

Sequence comparison of the mitochondrial genomes in two species of the genus Nerita (Gastropoda: Neritimorpha: Neritidae): phylogenetic implications and divergence time estimation for Neritimorpha

Many Nerita species live in warm-water environments, and they are some of the few organisms from the intertidal zone that can live in both freshwater and seawater. Previous comparative studies of the mitogenomes of Nerita species suggest that the genome rearrangements are very conservative. Generally, the species possess a set of similar mitochondrial gene arrangements, but nucleotide sequences can be used to elucidate phylogenetic relationships at various levels of divergence. Here, the mitogenomes of Nerita undata and Nerita balteata were sequenced and found to be 15,583 bp and 15,571 bp, respectively. The mitogenomes of both species contain 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The nucleotides of the two mitogenomes are highly similar, with the same gene composition and genomic organization as those present in other Nerita species. The tRNA secondary structures were different from those of other gastropods: trnS2 is not folded into typical secondary structures, and the dihydrouridine (DHU) arm simply forms a loop. The phylogenetic analysis showed that Neritimorpha is a sister group of Vetigastropoda and Caenogastropoda. Nerita balteata is a sister group of Nerita versicolor and Nerita undata, and all three species belong to Neritimorpha. This study contributes towards the comparative mitogenomic analysis of Neritidae and phylogenetic considerations among Neritimorpha species. The estimation of divergence time revealed that the two Nerita species were differentiated in the late Paleogene of the Cenozoic Era, and their evolution may be related to environmental changes.

From Persian Gulf to Indonesia: interrelated phylogeographic distance and chemistry within the genus Peronia (Onchidiidae, Gastropoda, Mollusca)

The knowledge of relationships between taxa is essential to understand and explain the chemical diversity of the respective groups. Here, twelve individuals of the panpulmonate slug Peronia persiae from two localities in Persian Gulf, and one animal of P. verruculata from Bangka Island, Indonesia, were analyzed in a phylogenetic and chemotaxonomic framework. Based on the ABGD test and haplotype networking using COI gene sequences of Peronia specimens, nine well-supported clades were found. Haplotype network analysis highlighted a considerable distance between the specimens of P. persiae and other clades. Metabolomic analysis of both species using tandem mass spectrometry-based GNPS molecular networking revealed a large chemical diversity within Peronia of different clades and localities. While P. persiae from different localities showed a highly similar metabolome, only few identical chemical features were found across the clades. The main common metabolites in both Peronia species were assigned as polypropionate esters of onchitriols and ilikonapyrones, and osmoprotectant amino acid-betaine compounds. On the other hand, the isoflavonoids genistein and daidzein were exclusively detected in P. persiae, while cholesterol and conjugated chenodeoxycholic acids were only found in P. verruculata. Flavonoids, bile acids, and amino acid-betaine compounds were not reported before from Onchidiidae, some are even new for panpulmonates. Our chemical analyses indicate a close chemotaxonomic relation between phylogeographically distant Peronia species.

The complete mitogenome of Arion vulgaris Moquin-Tandon, 1855 (Gastropoda: Stylommatophora): mitochondrial genome architecture, evolution and phylogenetic considerations within Stylommatophora

Stylommatophora is one of the most speciose orders of Gastropoda, including terrestrial snails and slugs, some of which are economically important as human food, agricultural pests, vectors of parasites or due to invasiveness. Despite their great diversity and relevance, the internal phylogeny of Stylommatophora has been debated. To date, only 34 stylommatophoran mitogenomes were sequenced. Here, the complete mitogenome of an invasive pest slug, Arion vulgaris Moquin-Tandon, 1855 (Stylommatophora: Arionidae), was sequenced using next generation sequencing, analysed and compared with other stylommatophorans. The mitogenome of A. vulgaris measures 14,547 bp and contains 13 protein-coding, two rRNA, 22 tRNA genes, and one control region, with an A + T content of 70.20%. All protein coding genes (PCGs) are initiated with ATN codons except for COX1, ND5 and ATP8 and all are ended with TAR or T-stop codons. All tRNAs were folded into a clover-leaf secondary structure except for trnC and trnS1 (AGN). Phylogenetic analyses confirmed the position of A. vulgaris within the superfamily Arionoidea, recovered a sister group relationship between Arionoidea and Orthalicoidea, and supported monophyly of all currently recognized superfamilies within Stylommatophora except for the superfamily Helicoidea. Initial diversification time of the Stylommatophora was estimated as 138.55 million years ago corresponding to Early Cretaceous. The divergence time of A. vulgaris and Arion rufus (Linnaeus, 1758) was estimated as 15.24 million years ago corresponding to one of Earth's most recent, global warming events, the Mid-Miocene Climatic Optimum. Furthermore, selection analyses were performed to investigate the role of different selective forces shaping stylommatophoran mitogenomes. Although purifying selection is the predominant selective force shaping stylommatophoran mitogenomes, six genes (ATP8, COX1, COX3, ND3, ND4 and ND6) detected by the branch-specific aBSREL approach and three genes (ATP8, CYTB and ND4L) detected by codon-based BEB, FUBAR and MEME approaches were exposed to diversifying selection. The positively selected substitutions at the mitochondrial PCGs of stylommatophoran species seems to be adaptive to environmental conditions and affecting mitochondrial ATP production or protection from reactive oxygen species effects. Comparative analysis of stylommatophoran mitogenome rearrangements using MLGO revealed conservatism in Stylommatophora; exceptions refer to potential apomorphies for several clades including rearranged orders of trnW-trnY and of trnE-trnQ-rrnS-trnM-trnL2-ATP8-trnN-ATP6-trnR clusters for the genus Arion. Generally, tRNA genes tend to be rearranged and tandem duplication random loss, transitions and inversions are the most basic mechanisms shaping stylommatophoran mitogenomes.

The first complete mitochondrial genome of a cyclophorid land snail, with implications for architaenioglossan relationships (Mollusca, Caenogastropoda, Cyclophoroidea)

Herein, we described the complete mitochondrial genome ('mitogenome') of the Chinese land snail Cyclophorus martensianus, which is the first published mitogenome for the caenogastropod family Cyclophoridae. This mitogenome is 15,308 bp long consisting of 37 genes: 13 protein-coding genes (PCGs), 22 tRNA genes and two rRNA genes. The A + T content (69.6%) is distinctly higher than the GC content (30.4%). Most PCGs start with ATN as initiation codons, and all PCGs have the conventional stop codons TAA and TAG. Overall, the biochemical properties and gene order of the cyclophorid mitogenome are very similar to those of other documented caenogastropod mitogenomes. We corroborate the findings of earlier studies that mitochondrial gene order is rather conserved among caenogastropods. Caenogastropoda are the taxonomically, morpho-anatomically and ecologically most diverse group of living gastropods comprising lineages inhabiting marine, freshwater, and terrestrial environments. Traditionally, the three most speciose groups of non-marine caenogastropods have all been placed in a single group, the Architaenioglossa. This group contains two freshwater lineages, the Viviparoidea and Ampullaroidea, and the terrestrial Cyclophoroidea. However, architaenioglossan relationships have remained controversial with several morphology-based on molecular phylogenetic analyses supporting contradicting phylogenetic hypotheses. In order to more robustly resolve the relationships among the Architaenioglossa, we have performed phylogenetic analyses of caenogastropod mitogenomes including the new mitogenome of Cyclophorus martensianus. Our phylogenetic reconstructions are based on the amino acid sequences of all protein-coding genes and consistently recovered Architaenioglossa as non-monophyletic.

A novel gene arrangement among the Stylommatophora by the complete mitochondrial genome of the terrestrial slug Meghimatium bilineatum (Gastropoda, Arionoidea)

Stylommatophora is a main clade of Gastropoda that encompasses approximately 112 gastropod families and may exceed a total of 30,000 species. Twenty-four complete stylommatophoran mitogenomes have been sequenced to date, yet our understanding of mitochondrial evolution in stylommatophorans is still in its infancy. To further expand the set of available mitogenomes, we sequenced the mitogenome of Meghimatium bilineatum (Arionoidea: Philomycidae), a widespread land slug in East Asia. This is the first report on a mitogenome of the superfamily Arionoidea, and indeed on a terrestrial slug. The mitogenome of Meghimatium bilineatum comprises 13,972 bp and exhibits a novel, highly distinctive gene arrangement among the Stylommatophora. Phylogenetic reconstructions based on the sequences of all protein-coding genes consistently recovered Meghimatium bilineatum as sister-group of the Succineidae. A phylogenetic reconstruction based on gene order, however, suggested a highly divergent tree topology, which is less credible when taking into account prior knowledge of stylommatophoran relationships. Our CREx (Common interval Rearrangement Explorer) analysis suggested that three successive events of tandem duplication random loss (TDRL) best explain the evolutionary process of gene order rearrangement in Meghimatium bilineatum from an ancestral stylommatophoran mitogenome. The present example offers new insights into the mechanisms of mitogenome rearrangements in gastropods at large and into the usefulness of mitogenomic gene order as a phylogenetic marker.

Incorporation of deep-sea and small-sized species provides new insights into gastropods phylogeny

The use of phylogeny with uneven or limited taxon sampling may bias our interpretation of organismal evolution, for instance, the origin(s) of the deep-sea animals. The Mollusca is the second most speciose phylum, in which the Gastropoda forms the largest group. However, the currently proposed hypotheses of gastropod phylogeny are mainly based on part of their taxonomic diversity, notably on the large-sized and shallow-water species. In this study, we aimed at correcting this bias by reconstructing the phylogeny with new mitogenomes of deep-sea gastropods including Anatoma sp., Bathysciadiidae sp., Bayerotrochus teramachii, Calliotropis micraulax, Coccocrater sp., Cocculina subcompressa, Lepetodrilus guaymasensis, Peltospira smaragdina, Perotrochus caledonicus, Pseudococculinidae sp., and Shinkailepas briandi. This dataset provided the first reports of the mitogenomes for the Cocculiniformia, three vetigastropod superfamilies: Pleurotomarioidea, Lepetelloidea, and Scissurelloidea, and the neritimorph family Phenacolepadidae. The addition of deep-sea representatives also allowed us to evaluate the evolution of habitat use in gastropods. Our results showed a strongly supported sister-group relationship between the deep-sea lineages Cocculiniformia and Neomphalina. Within the Vetigastropoda, the Pleurotomarioidea was revealed as the sister-group of the remaining vetigastropods. Although this clade was presently restricted to the deep sea, fossil records showed that it has only recently invaded this habitat, thus suggesting that shallow waters was the ancestral habitat for the Vetigastropoda. The deep-sea Lepetelloidea and Lepetodriloidea formed a well-supported clade, with the Scissurelloidea sister to it, suggesting an early transition from shallow water to deep sea in this lineage. In addition, the switch between different chemosynthetic habitats was also observed in deep-sea gastropod lineages, notably in Neomphalina and Lepetelloidea. In both cases, the biogenic substrates appeared as the putative ancestral habitat, confirming the previously proposed hypothesis of a wooden-step to deep-sea vents scenario of evolution of habitat use for these taxa.

The complete mitogenome of Helixpomatia and the basal phylogeny of Helicinae (Gastropoda, Stylommatophora, Helicidae)

A complete mitochondrial genome of the Roman snail Helixpomatia Linnaeus, 1758 has been sequenced. The length and gene order correspond to that of other available helicid mitogenomes. We used the mitogenome sequence to reappraise the relationships among the four presumed principal groups of the helicid subfamily Helicinae. The results support the idea that the subfamily is divided between two western Palaearctic diversification centres: Iberian Peninsula and western Maghreb in the west, and Anatolia, the Aegean and Caucasus in the east. One group, the tribe Helicini, diversified in the east and the remaining three currently recognised tribes in the west. However, the exact relationships among lineages of the non-Helicini tribes could not be resolved.

The Complete Mitochondrial Genome Sequences of the Philomycus bilineatus (Stylommatophora: Philomycidae) and Phylogenetic Analysis

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.

New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks

Long-branch attraction (LBA) is a well-known artifact in phylogenetic reconstruction. Sparse taxon sampling and extreme heterogeneity of evolutionary rates among lineages generate propitious situations for LBA, even defying probabilistic methods of phylogenetic inference. A clear example illustrating LBA challenges is the difficulty of reconstructing the deep gastropod phylogeny, particularly using mitochondrial (mt) genomes. Previous studies consistently obtained unorthodox phylogenetic relationships due to the LBA between the mitogenomes of patellogastropods (true limpets, represented only by Lottia digitalis), heterobranchs, and outgroup taxa. Here, we use the reconstruction of the gastropod mitogenomic phylogeny as a case exercise to test the effect of key methodological approaches proposed to counteract LBA, including the selection of slow-evolving representatives, the use of different outgroups, the application of site-heterogeneous evolutionary models, and the removal of fast-evolving sites. In this regard, we sequenced three new patellogastropod mt genomes, which displayed shorter branches than the one of Lottia as well as gene organizations more similar to that of the hypothetical gastropod ancestor. Phylogenetic analyses incorporating the mt genomes of Patella ferruginea, Patella vulgata, and Cellana radiata allowed eliminating the artificial clustering of Patellogastropoda and Heterobranchia that had prevailed in previous studies. Furthermore, the use of site-heterogeneous models with certain combinations of lineages within the outgroup allowed eliminating also the LBA between Heterobranchia and the outgroup, and recovering Apogastropoda (i.e., Caenogastropoda + Heterobranchia). Hence, for the first time, we were able to obtain a mitogenomic phylogeny of gastropods that is congruent with both morphological and nuclear datasets.

Comparative Characterization of the Complete Mitochondrial Genomes of the Three Apple Snails (Gastropoda: Ampullariidae) and the Phylogenetic Analyses

The apple snails Pomacea canaliculata, Pomacea diffusa and Pomacea maculate (Gastropoda: Caenogastropoda: Ampullariidae) are invasive pests causing massive economic losses and ecological damage. We sequenced and characterized the complete mitochondrial genomes of these snails to conduct phylogenetic analyses based on comparisons with the mitochondrial protein coding sequences of 47 Caenogastropoda species. The gene arrangements, distribution and content were canonically identical and consistent with typical Mollusca except for the tRNA-Gln absent in P. diffusa. An identifiable control region (d-loop) was absent. Bayesian phylogenetic analysis indicated that all the Ampullariidae species clustered on the same branch. The genus Pomacea clustered together and then with the genus Marisa. The orders Architaenioglossa and Sorbeoconcha clustered together and then with the order Hypsogastropoda. Furthermore, the intergenic and interspecific taxonomic positions were defined. Unexpectedly, Ceraesignum maximum, Dendropoma gregarium, Eualetes tulipa and Thylacodes squamigerus, traditionally classified in order Hypsogastropoda, were isolated from the order Hypsogastropoda in the most external branch of the Bayesian inference tree. The divergence times of the Caenogastropoda indicated that their evolutionary process covered four geological epochs that included the Quaternary, Neogene, Paleogene and Cretaceous periods. This study will facilitate further investigation of species identification to aid in the implementation of effective management and control strategies of these invasive species.

The mitochondrial genome of the planorbid snail Planorbella duryi

The complete mitochondrial genome of a freshwater planorbid snail, Planorbella duryi (Mollusca, Gastropoda) was recovered from de novo assembly of genomic sequences generated with the Illumina NextSeq500 platform. The P. duryi mitogenome (14,217 base pairs) is AT rich (72.69%) and comprises 13 protein-coding genes, two ribosomal subunit genes, and 22 transfer RNAs. The gene order is identical to that of Biomphalaria glabrata and other snail species in the family Planorbidae. This is the first full characterization of a mitochondrial genome of the genus Planorbella.

Comparative immunological study of the snail Physella acuta (Hygrophila, Pulmonata) reveals shared and unique aspects of gastropod immunobiology

The freshwater snail Physella acuta was selected to expand the perspective of comparative snail immunology. Analysis of Physella acuta, belonging to the Physidae, taxonomic sister family to Planorbidae, affords family-level comparison of immune features characterized from Biomphalaria glabrata, the model snail often used to interpret general gastropod immunity. To capture constitutive and induced immune sequences, transcriptomes of an individual Physella acuta snail, 12 h post injection with bacteria (Gram -/+) and one sham-exposed snail were recorded with 454 pyrosequencing. Assembly yielded a combined reference transcriptome containing 24,288 transcripts. Additionally, genomic Illumina reads were obtained (∼15-fold coverage). Recovery of transcripts for two macin-like antimicrobial peptides (AMPs), 12 aplysianins, four LBP/BPIs and three physalysins indicated that Physella acuta shares a similar organization of antimicrobial defenses with Biomphalaria glabrata, contrasting a modest AMP arsenal with a diverse set of antimicrobial proteins. The lack of predicted transmembrane domains in all seven Physella acuta PGRP transcripts supports the notion that gastropods do not employ cell-bound PGRP receptors, different from ecdysozoan invertebrates yet similar to mammals (vertebrate deuterostomes). The well-documented sequence diversification by Biomphalaria glabrata FREPs (immune lectins comprising immunoglobulin superfamily domains and fibrinogen domains), resulting from somatic mutations of a large FREP gene family is hypothesized to be unique to Planorbidae; Physella acuta revealed just two bonafide FREP genes and these were not diversified. Furthermore, the flatworm parasite Echinostoma paraensei, confirmed here to infect both snail species, did not evoke from Physella acuta the abundant expression of FREP proteins at 2, 4 and 8 days post exposure that was previously observed from Biomphalaria glabrata. The Physella acuta reference transcriptome also revealed 24 unique transcripts encoding proteins consisting of a single fibrinogen-related domain (FReDs), with a short N-terminal sequence encoding either a signal peptide, transmembrane domain or no predicted features. The Physella acuta FReDs are candidate immune genes based on implication of similar sequences in immunity of bivalve molluscs. Overall, comparative analysis of snails of sister families elucidated the potential for taxon-specific immune features and investigation of strategically selected species will provide a more comprehensive view of gastropod immunity.

A comparison of mitochondrial genomes from five species in three genera suggests polyphyly in the subfamily Achatinellinae (Gastropoda: Pulmonata: Stylommatophora: Achatinellidae)

We compare the complete mitochondrial genomes of Achatinella fulgens, A. mustelina, A. sowerbyana, Partulina redfieldi, and Perdicella helena, five species of Hawaiian tree snails across three genera. Mitogenomes ranged in length from 15,187 to 16,793 base pairs, with a base composition of A (36.4-37.4%); T (42.2-42.7%); C (8.8-9.2%); and G (11.3-11.8%). Similar with other pulmonates, these mitogenomes contain 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes, with the order conserved among genera. Our study suggests polyphyly in the current arrangement of the subfamily Achatinellinae, part of a spectacular radiation in the Hawaiian Islands.

A revision of Peronina Plate, 1893 (Gastropoda : Euthyneura : Onchidiidae) based on mitochondrial and nuclear DNA sequences, morphology and natural history

Peronina Plate, 1893 is a genus of onchidiids that live on the mud in mangrove forests. Peronina can be identified in the field by the lung opening at the margin between the ventral hyponotum and the dorsal notum, and by the distinctive scalloped notum edge. This genus was previously known only from the holotype of the type species, Peronina alta Plate, 1893, from eastern India. Onchidium tenerum Stoliczka, 1869 is moved to Peronina and applies to the same species as Peronina alta. Peronina species are described using an integrative approach (natural history, comparative anatomy and DNA sequences). Mitochondrial COI and 16S sequences and nuclear ITS2 and 28S sequences are used to independently test species boundaries. Mitochondrial sequences yielded three units separated by a large barcode gap, but nuclear sequences yielded two units. Because these two units are congruent with differences in the male copulatory apparatus, they are accepted as species. Explanations for highly divergent COI haplotypes within one species are discussed. Peronina tenera (Stoliczka, 1869) is distributed in the Bay of Bengal and the Strait of Malacca, while P. zulfigari Goulding & Dayrat, sp. nov. is endemic to the Strait of Malacca. The two species differ internally but are cryptic externally.

Complete sequence analysis of the mitochondrial genome of Auriculastra duplicata (Mollusca, Gastropoda, Ellobiidae)

The mitochondrial genome of the gastropod Auriculastra duplicata was completely sequenced. It was 13,920 bp in length and comprised 37 genes; two rrn genes and 22 trn genes. Phylogenetic analyses based on the concatenated protein-coding genes depicted the polyphyly of all species belonging to the family Ellobiidae; however, monophyly was observed among all species belonging to the subfamily Ellobiinae, in which A. duplicata clustered consistently with Auriculinella bidentata.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.

The complete mitogenome of Orcula dolium (Draparnaud, 1801); ultra-deep sequencing from a single long-range PCR using the Ion-Torrent PGM

BACKGROUND

With the increasing capacity of present-day next-generation sequencers the field of mitogenomics is rapidly changing. Enrichment of the mitochondrial fraction, is no longer necessary for obtaining mitogenomic data. Despite the benefits, shotgun sequencing approaches also have disadvantages. They do not guarantee obtaining the complete mitogenome, generally require larger amounts of input DNA and coverage is low compared to sequencing with enrichment strategies. If the mitogenome could be amplified in a single amplification, additional time and costs for sample preparation might outweigh these disadvantages.

RESULTS

A sequence of the complete mitochondrial genome of the pupilloid landsnail Orcula dolium is presented. The mitogenome was amplified in a single long-range (LR) PCR and sequenced on an Ion Torrent PGM (Life Technologies). The length is 14,063 nt and the average depth of coverage is 1112 X. This is the first published mitogenome for a member of the family Orculidae. It has the typical metazoan makeup of 13 protein coding genes (PCGs), 2 ribosomal RNAs (12S and 16S) and 22 transfer RNAs (tRNAs). Orcula is positioned between Pupilla and the Vertiginidae as the sister-group of Gastrocopta and Vertigo, together. An ancestral gene order reconstruction shows that Orthurethra in contrast to other Stylommatophora, have tRNA-H before tRNA-G and that the gene order in the 'non-achatinoid' clade is identical to that of closely related non-stylommatophoran taxa.

CONCLUSIONS

We show it is feasible to ultra-deep sequence a mitogenome from a single LR-PCR. This approach is particularly relevant to studies that have low concentrations of input DNA. It results in a more efficient use of NGS capacity (only the targeted fraction is sequenced) and is an effective selection against nuclear mitochondrial inserts (NUMTS). In contrast to previous studies based in particular on 28S, our results indicate that phylogeny reconstructions based on complete mitogenomes might be more suitable to resolve deep relationships within Stylommatophora. Ancestral gene order reconstructions reveal rearrangements that characterize systematic groups.

Burrowers from the Past: Mitochondrial Signatures of Ordovician Bivalve Infaunalization

Bivalves and gastropods are the two largest classes of extant molluscs. Despite sharing a huge number of features, they do not share a key ecological one: gastropods are essentially epibenthic, although most bivalves are infaunal. However, this is not the ancestral bivalve condition; Cambrian forms were surface crawlers and only during the Ordovician a fundamental infaunalization process took place, leading to bivalves as we currently know them. This major ecological shift is linked to the exposure to a different redox environoments (hypoxic or anoxic) and with the Lower Devonian oxygenation event. We investigated selective signatures on bivalve and gastropod mitochondrial genomes with respect to a time calibrated mitochondrial phylogeny by means of dN/dS ratios. We were able to detect 1) a major signal of directional selection between the Ordovician and the Lower Devonian for bivalve mitochondrial Complex I, and 2) an overall higher directional selective pressure on bivalve Complex V with respect to gastropods. These and other minor dN/dS patterns and timings are discussed, showing that the Ordovician infaunalization event left heavy traces in bivalve mitochondrial genomes.

The complete mitochondrial genomes of two freshwater snails provide new protein-coding gene rearrangement models and phylogenetic implications

BACKGROUND

Mitochondrial (mt) genome sequences are widely used for species identification and to study the phylogenetic relationships among Gastropoda. However, to date, limited data are available as taxon sampling is narrow. In this study we sequenced the complete mt genomes of the freshwater gastropods Radix swinhoei (Lymnaeidae) and Planorbarius corneus (Planorbidae). Based on these sequences, we investigated the gene rearrangement in these two species and the relationships with respect to the ancestral gene order and assessed their phylogenetic relationships.

METHODS

The complete mt genomes of R. swinhoei and P. corneus were sequenced using Illumina-based paired-end sequencing and annotated by comparing the sequence information with that of related gastropod species. Putative models of mitochondrial gene rearrangements were predicted for both R. swinhoei and P. corneus, using Reishia clavigera mtDNA structure as the ancestral gene order. The phylogenetic relationships were inferred using thirteen protein sequences based on Maximum likelihood and Bayesian inference analyses.

RESULTS

The complete circular mt genome sequences of R. swinhoei and P. corneus were 14,241 bp and 13,687 bp in length, respectively. Comparison of the gene order demonstrated complex rearrangement events in Gastropoda, both for tRNA genes and protein-coding genes. The phylogenetic analyses showed that the family Lymnaeidae was more closely related to the family Planorbidae, consistent with previous classification. Nevertheless, due to the position recovered for R. swinhoei, the family Lymnaeidae was not monophyletic.

CONCLUSION

This study provides the complete mt genomes of two freshwater snails, which will aid the development of useful molecular markers for epidemiological, ecological and phylogenetic studies. Additionally, the predicted models for mt gene rearrangement might provide novel insights into mt genome evolution in gastropods.

Beyond Conus: Phylogenetic relationships of Conidae based on complete mitochondrial genomes

Understanding how the extraordinary taxonomic and ecological diversity of cone snails (Caenogastropoda: Conidae) evolved requires a statistically robust phylogenetic framework, which thus far is not available. While recent molecular phylogenies have been able to distinguish several deep lineages within the family Conidae, including the genera Profundiconus, Californiconus, Conasprella, and Conus (and within this one, several subgenera), phylogenetic relationships among these genera remain elusive. Moreover, the possibility that additional deep lineages may exist within the family is open. Here, we reconstructed with probabilistic methods a molecular phylogeny of Conidae using the newly sequenced complete or nearly complete mitochondrial (mt) genomes of the following nine species that represent all main Conidae lineages and potentially new ones: Profundiconus teramachii, Californiconus californicus, Conasprella wakayamaensis, Lilliconus sagei, Pseudolilliconus traillii, Conus (Kalloconus) venulatus, Conus (Lautoconus) ventricosus, Conus (Lautoconus) hybridus, and Conus (Eugeniconus) nobilis. To test the monophyly of the family, we also sequenced the nearly complete mt genomes of the following three species representing closely related conoidean families: Benthomangelia sp. (Mangeliidae), Tomopleura sp. (Borsoniidae), and Glyphostoma sp. (Clathurellidae). All newly sequenced conoidean mt genomes shared a relatively constant gene order with rearrangements limited to tRNA genes. The reconstructed phylogeny recovered with high statistical support the monophyly of Conidae and phylogenetic relationships within the family. The genus Profundiconus was placed as sister to the remaining genera. Within these, a clade including Californiconus and Lilliconus+Pseudolilliconus was the sister group of Conasprella to the exclusion of Conus. The phylogeny included a new lineage whose relative phylogenetic position was unknown (Lilliconus) and uncovered thus far hidden diversity within the family (Pseudolilliconus). Moreover, reconstructed phylogenetic relationships allowed inferring that the peculiar diet of Californiconus based on worms, mollusks, crustaceans and fish is derived, and reinforce the hypothesis that the ancestor of Conidae was a worm hunter. A chronogram was reconstructed under an uncorrelated relaxed molecular clock, which dated the origin of the family shortly after the Cretaceous-Tertiary boundary (about 59million years ago) and the divergence among main lineages during the Paleocene and the Eocene (56-30million years ago).

Two complete mitochondrial genomes from Praticolella mexicana Perez, 2011 (Polygyridae) and gene order evolution in Helicoidea (Mollusca, Gastropoda)

Helicoidea is a diverse group of land snails with a global distribution. While much is known regarding the relationships of helicoid taxa, comparatively little is known about the evolution of the mitochondrial genome in the superfamily. We sequenced two complete mitochondrial genomes from Praticolellamexicana Perez, 2011 representing the first such data from the helicoid family Polygyridae, and used them in an evolutionary analysis of mitogenomic gene order. We found the mitochondrial genome of Praticolellamexicana to be 14,008 bp in size, possessing the typical 37 metazoan genes. Multiple alternate stop codons are used, as are incomplete stop codons. Mitogenome size and nucleotide content is consistent with other helicoid species. Our analysis of gene order suggested that Helicoidea has undergone four mitochondrial rearrangements in the past. Two rearrangements were limited to tRNA genes only, and two involved protein coding genes.

The complete mitochondrial genome of Achatinella sowerbyana (Gastropoda: Pulmonata: Stylommatophora: Achatinellidae)

In this study, we report the complete mitochondrial genome sequence of Achatinella sowerbyana, an endangered Hawaiian tree snail. The mitogenome is 15,374 bp in length and has a base composition of A (36.4%), T (42.7%), C (9.1%), and G (11.8%). Similar to other pulmonates, it contains 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. The gene order is the same as its sister species, A. mustelina. To our knowledge, this is the second mitochondrial genome sequenced within the superfamily Achatinelloidea, and will aid in the examination of the genus Achatinella, part of a spectacular radiation in the Hawaiian Islands.

Positive selection on panpulmonate mitogenomes provide new clues on adaptations to terrestrial life

BACKGROUND

Transitions from marine to intertidal and terrestrial habitats resulted in a significant adaptive radiation within the Panpulmonata (Gastropoda: Heterobranchia). This clade comprises several groups that invaded the land realm independently and in different time periods, e.g., Ellobioidea, Systellomatophora, and Stylommatophora. Thus, mitochondrial genomes of panpulmonate gastropods are promising to screen for adaptive molecular signatures related to land invasions.

RESULTS

We obtained three complete mitochondrial genomes of terrestrial panpulmonates, i.e., the ellobiid Carychium tridentatum, and the stylommatophorans Arion rufus and Helicella itala. Our dataset consisted of 50 mitogenomes comprising almost all major panpulmonate lineages. The phylogenetic tree based on mitochondrial genes supports the monophyly of the clade Panpulmonata. Terrestrial lineages were sampled from Ellobioidea (1 sp.) and Stylommatophora (9 spp.). The branch-site test of positive selection detected significant non-synonymous changes in the terrestrial branches leading to Carychium (Ellobiodea) and Stylommatophora. These convergent changes occurred in the cob and nad5 genes (OXPHOS complex III and I, respectively).

CONCLUSIONS

The convergence of the non-synonymous changes in cob and nad5 suggest possible ancient episodes of positive selection related to adaptations to non-marine habitats. The positively selected sites in our data are in agreement with previous results in vertebrates suggesting a general pattern of adaptation to the new metabolic requirements. The demand for energy due to the colonization of land (for example, to move and sustain the body mass in the new habitat) and the necessity to tolerate new conditions of abiotic stress may have changed the physiological constraints in the early terrestrial panpulmonates and triggered adaptations at the mitochondrial level.

The mitochondrial genome of the land snail Cernuella virgata (Da Costa, 1778): the first complete sequence in the family Hygromiidae (Pulmonata, Stylommatophora)

The land snail Cernuella virgata (da Costa, 1778) is widely considered as a pest to be quarantined in most countries. In this study, the complete mitochondrial genome of Cernuella virgata is published. The mitochondrial genome has a length of 14,147 bp a DNA base composition of 29.07% A, 36.88% T, 15.59% C and 18.46% G, encoding 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes and two ribosomal RNA (rRNA) genes. The complete nucleotide composition was biased toward adenine and thymine, A+T accounting for 69.80%. Nine PCGs and 14 tRNA genes are encoded on the J strand, and the other four PCGs and eight tRNA genes are encoded on the N strand. The genome also includes 16 intergenic spacers. All PCGs start strictly with ATN, and have conventional stop codons (TAA and TAG). All tRNAs fold into the classic cloverleaf structure, except tRNA(Arg) , tRNA(Ser(UCN)) , tRNA(Ser(AGN)) and tRNA(Pro) . The first three lack the dihydrouridine arm while the last lacks the TψC arm. There are 502 bp long noncoding regions and 418bp long gene overlaps in the whole mitochondrial genome, accounting for 3.54% and 2.95% of the total length respectively. Phylogenetic analyses based on the sequences of the protein coding genes revealed a sister group relationship between the Hygromiidae and the Helicidae.

The complete mitochondrial genome of Achatinella mustelina (Gastropoda: Pulmonata: Stylommatophora)

In this study, we report the complete mitochondrial genome sequence of Achatinella mustelina, an endangered Hawaiian tree snail. The mitogenome is 16 323 bp in length and has a base composition of A (34.7%), T (42.6%), C (12.7%) and G (10.0%). Similar to other Pulmonates, it contains 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. To our knowledge, this is the first mitochondrial genome sequenced within the Achatinelloidea superfamily, which contains a high number of endangered species. As such, this mitogenome will be useful in conservation genetics studies.

Complete mitochondrial genome of the Endangered species Ellobium chinense (Pulmonata, Ellobiidae) from Korea

The complete mitochondrial genome of Ellobium chinense (Ellobioidea, Ellobiidae), an Endangered species in South Korea, is reported here for the first time. The mitogenome of E. chinense is 13,979 base pairs in total length and includes 13 PCGs, small and large rRNAs, and 21 tRNAs. Twelve genes are encoded on the light-strand and 24 genes on the heavy-strand. Compared to four other ellobiid species, the PCGs of E. chinense have a conserved gene order except for the positions of ND4L and ND4. These data provide useful molecular information for phylogenetic studies concerning ellobiids and related species.