The mitochondrial genome of Conus textile, coxI-coxII intergenic sequences and Conoidean evolution

Comparative Mitogenome Analysis of Two Native Apple Snail Species (Ampullariidae, Pomacea) from Peruvian Amazon

Apple snails of the genus Pomacea Perry, 1810 (Mollusca: Caenogastropoda: Ampullariidae) are native to the Neotropics and exhibit high species diversity, holding cultural and ecological significance as an important protein source in Peru. However, most genetic studies in Pomacea have focused mostly on invasive species, especially in Southeast Asia, where they are considered important pests. In this study, we assembled and annotated the mitochondrial genomes of two Pomacea species native to the Peruvian Amazon: Pomacea reevei Ampuero & Ramírez, 2023 and Pomacea aulanieri (Deville & Hupé, 1850). The mitogenomes of P. reevei and P. aulanieri comprise 15,660 and 16,096 bp, respectively, and contain the typical 37 genes of the animal mitochondria with a large control region of 292 bp in P. reevei and 524 bp in P. aulanieri-which fall within the range of what is currently known in Pomacea. Comparisons with previously published mitogenomes in Pomacea revealed differences in the overlapping of adjacent genes, the size of certain protein-coding genes (PCGs) and the secondary structure of some tRNAs that are consistent with the phylogenetic relationships between these species. These findings provide valuable insights into the systematics and genomics of the genus Pomacea.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Utilisation of a mitochondrial intergenic region for species differentiation of fruit flies (Diptera: Tephritidae) in South Africa

Background

Fruit flies (Diptera: Tephritidae) comprise species of agricultural and economic importance. Five such fruit fly species are known to affect commercial fruit production and export in South Africa: Ceratitis capitata, Ceratitis cosyra, Ceratitis rosa, Ceratitis quilicii, and Bactrocera dorsalis. Management practices for these pests include monitoring, application of pest control products, post-harvest disinfestation measures and inspection of consignments both prior to shipment and at ports of entry. In activities relating to monitoring and inspection, accurate identification of these pests to species level is required. While morphological keys for adult stages of these fruit fly species have been well developed, morphological keys for earlier life stages remain problematic. In instances where closely related species cannot be reliably distinguished morphologically, there is a need for molecular tools to assist in identifying these five fruit fly species during surveillance practices, where sequencing-based approaches would be beneficial.

Results

Two complete mitochondrial genomes were assembled for each fruit fly species investigated using high throughput sequencing data generated in this study. A single primer set was designed to amplify a region between tRNAile and tRNAmet. The amplicon consists of a partial segment of tRNAile, intergenic region I (tRNAile - tRNAgln), the complete sequence of tRNAgln, intergenic region II (tRNAgln - tRNAmet), and a partial segment of tRNAmet. PCR amplicons were generated for 20 specimens of each species, five of which were colony adult males, five colony larvae, and 10 wild, trap-collected specimens. Upon analysis of the amplicon, intergenic region I was identified as the most informative region, allowing for unambiguous identification of the five fruit fly species. The similarity in intergenic region II was too high between C. rosa and C. quilicii for accurate differentiation of these species.

Conclusion

The identity of all five fruit flies investigated in this study can be determined through sequence analysis of the mitochondrial intergenic regions. Within the target amplicon, intergenic region I (tRNAile - tRNAgln) shows interspecific variation sufficient for species differentiation based on multiple sequence alignment. The variation in the length of intergenic region I is proposed as a potential tool for accurately identifying these five fruit flies in South Africa.

Complete mitochondrial genome of Conus lischkeanus Weinkauff, 1875 (Neogastropoda, Conidae) and phylogenetic implications of the evolutionary diversification of dietary types of Conus species

The family Conidae, commonly known as cone snails, is one of the most intriguing gastropod groups owing to their diverse array of feeding behaviors (diets) and toxin peptides (conotoxins). Conuslischkeanus Weinkauff, 1875 is a worm-hunting species widely distributed from Africa to the Northwest Pacific. In this study, we report the mitochondrial genome sequence of C.lischkeanus and inferred its phylogenetic relationship with other Conus species. Its mitochondrial genome is a circular DNA molecule (16,120 bp in size) composed of 37 genes: 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes. Phylogenetic analyses of concatenated nucleotide sequences of 13 PCGs and two ribosomal RNA genes showed that C.lischkeanus belongs to the subgenus Lividoconus group, which is grouped with species of the subgenus Virgiconus, and a member of the largest assemblage of worm-hunting (vermivorous) species at the most basal position in this group. Mitochondrial genome phylogeny supports the previous hypothesis that the ancestral diet of cone snails was worm-hunting, and that other dietary types (molluscivous or piscivorous) have secondarily evolved multiple times from different origins. This new, complete mitochondrial genome information provides valuable insights into the mitochondrial genome diversity and molecular phylogeny of Conus species.

The complete mitochondrial DNA genome of a cone snail, Conus betulinus (Neogastropoda: Conidae), from the South China sea

The complete mitochondrial genome of the tubular cone snail Conus betulinus is presented in this study. The C. betulinus mitochondrial genome was 16,240 bp with 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a non-coding AT-rich region (D-loop). The overall base composition was estimated to be 25.67% for A, 38.26% for T, 21.38% for G, and 14.69% for C, with a high A + T content of 63.93%. Phylogenetic analyses based on 13 PCGs showed the close relationship of vermivorous C. betulinus with the common ancestor of molluscivorous Conus textile and Conus gloriamaris, providing a basis for further studies on the phylogenetics of cone snails according to their dietary type.

Chemical Synthesis and NMR Solution Structure of Conotoxin GXIA from Conus geographus

Conotoxins are disulfide-rich peptides found in the venom of cone snails. Due to their exquisite potency and high selectivity for a wide range of voltage and ligand gated ion channels they are attractive drug leads in neuropharmacology. Recently, cone snails were found to have the capability to rapidly switch between venom types with different proteome profiles in response to predatory or defensive stimuli. A novel conotoxin, GXIA (original name G117), belonging to the I₃-subfamily was identified as the major component of the predatory venom of piscivorous Conus geographus. Using 2D solution NMR spectroscopy techniques, we resolved the 3D structure for GXIA, the first structure reported for the I₃-subfamily and framework XI family. The 32 amino acid peptide is comprised of eight cysteine residues with the resultant disulfide connectivity forming an ICK+1 motif. With a triple stranded β-sheet, the GXIA backbone shows striking similarity to several tarantula toxins targeting the voltage sensor of voltage gated potassium and sodium channels. Supported by an amphipathic surface, the structural evidence suggests that GXIA is able to embed in the membrane and bind to the voltage sensor domain of a putative ion channel target.

The complete mitochondrial genome of a cold seep gastropod Phymorhynchus buccinoides (Neogastropoda: Conoidea: Raphitomidae)

Phymorhynchus is a genus of deep-sea snails that are most distributed in hydrothermal vent or cold seep environments. In this study, we presented the complete mitochondrial genome of P. buccinoides, a cold seep snail from the South China Sea. It is the first mitochondrial genome of a cold seep member of the superfamily Conoidea. The mitochondrial genome is 15,764 bp in length, and contains 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes. These genes are encoded on the positive strand, except for 8 tRNA genes that are encoded on the negative strand. The start codon ATG and 3 types of stop codons, TAA, TAG and the truncated termination codon T, are used in the 13 PCGs. All 13 PCGs in the 26 species of Conoidea share the same gene order, while several tRNA genes have been translocated. Phylogenetic analysis revealed that P. buccinoides clustered with Typhlosyrinx sp., Eubela sp., and Phymorhynchus sp., forming the Raphitomidae clade, with high support values. Positive selection analysis showed that a residue located in atp6 (18 S) was identified as the positively selected site with high posterior probabilities, suggesting potential adaption to the cold seep environment. Overall, our data will provide a useful resource on the evolutionary adaptation of cold seep snails for future studies.

Prevention of scar hyperplasia in the skin by conotoxin: A prospective review

Scars are often considered to be skin problems that affect beauty. The tension acting on the edge of the wound is the main factor causing the scar hyperplasia. At present, the clinical use of botulinum toxin A (BTX-A) around the wound to cause transient muscle paralysis reduce the muscle movement around the wound and wound tension to prevent scar hyperplasia during wound healing. But the use of BTX-A to prevent scarring requires the use of a syringe. The syringe can cause trauma and pain when it pricks the skin for BTX-A injection. The conotoxin which is secreted by the poison glands on the inside of the venom tube and capsule of the snail provides a simple and effective way to prevent skin scar hyperplasia. We reviewed the classification of conotoxin, the conotoxin's mechanism of preventing scar hyperplasia, and the research direction of conotoxin in the future and provided reference for promoting the application of conotoxin in preventing skin scar hyperplasia.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Mitochondrial genome sequencing of a vermivorous cone snail Conus quercinus supports the correlative analysis between phylogenetic relationships and dietary types of Conus species

Complete mitochondrial genome (mitogenome) sequence of a worm-hunting cone snail, Conus quercinus, was reported in this study. Its mitogenome, the longest one (16,460 bp) among reported Conus specie, is composed of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one D-loop region. The mitochondrial gene arrangement is highly-conserved and identical to other reported. However, the D-loop region of C. quercinus is the longest (943 bp) with the higher A+T content (71.3%) and a long AT tandem repeat stretch (68 bp). Subsequent phylogenetic analysis demonstrated that three different dietary types (vermivorous, molluscivorous and piscivorous) of cone snails are clustered separately, suggesting that the phylogenetics of cone snails is related to their dietary types. In conclusion, our current work improves our understanding of the mitogenomic structure and evolutionary status of the vermivorous C. quercinus, which support the putative hypothesis that the Conus ancestor was vermivorous.

Cone Snails: A Big Store of Conotoxins for Novel Drug Discovery

Marine drugs have developed rapidly in recent decades. Cone snails, a group of more than 700 species, have always been one of the focuses for new drug discovery. These venomous snails capture prey using a diverse array of unique bioactive neurotoxins, usually named as conotoxins or conopeptides. These conotoxins have proven to be valuable pharmacological probes and potential drugs due to their high specificity and affinity to ion channels, receptors, and transporters in the nervous systems of target prey and humans. Several research groups, including ours, have examined the venom gland of cone snails using a combination of transcriptomic and proteomic sequencing, and revealed the existence of hundreds of conotoxin transcripts and thousands of conopeptides in each Conus species. Over 2000 nucleotide and 8000 peptide sequences of conotoxins have been published, and the number is still increasing quickly. However, more than 98% of these sequences still lack 3D structural and functional information. With the rapid development of genomics and bioinformatics in recent years, functional predictions and investigations on conotoxins are making great progress in promoting the discovery of novel drugs. For example, ω-MVIIA was approved by the U.S. Food and Drug Administration in 2004 to treat chronic pain, and nine more conotoxins are at various stages of preclinical or clinical evaluation. In short, the genus Conus, the big family of cone snails, has become an important genetic resource for conotoxin identification and drug development.

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).

Characterization of the complete mitochondrial genome of Conus tribblei Walls, 1977

The genus Conus sensu lato consists of 500-700 species. However, the mitochondrial genomes of only few species have been fully sequenced and reported so far. In this study, the complete mitochondrial genome of Conus tribblei, a member of the poorly known subgenus Splinoconus is sequenced with the mean coverage of 604×. The mitochondrial genome is 15 570 bp long and consists of genes encoding for 13 respiratory chain proteins, 22 tRNA and 2 rRNA. The gene organization is highly conserved among the Conus species. The longest intergenic region between tRNA-Phe and cytochrome c oxidase subunit III (cox3), which in C. tribblei is 169 bp long and contains a 112 bp long segment of inverted repeat, represents the putative control region. The control regions of Conus species exhibited variability in the length and position of the inverted repeats. Therefore, this region may have the potential to be used as a genetic marker for species discrimination.

One, four or 100 genera? A new classification of the cone snails

We present a new classification for the genus Conus sensu lato (family Conidae), based on molecular phylogenetic analyses of 329 species. This classification departs from both the traditional classification in only one genus and from a recently proposed shell- and radula-based classification scheme that separates members of this group into five families and 115 genera. Roughly 140 genus-group names are available for Recent cone snails. We propose to place all cone snails within a single family (Conidae) containing four genera-Conus, Conasprella, Profundiconus and Californiconus (with Conus alone encompassing about 85% of known species)-based on the clear separation of cone snails into four distinct and well-supported groups/lineages in molecular phylogenetic analyses. Within Conus and Conasprella, we recognize 57 and 11 subgenera, respectively, that represent well-supported subgroupings within these genera, which we interpret as evidence of intrageneric distinctiveness. We allocate the 803 Recent species of Conidae listed as valid in the World Register of Marine Species into these four genera and 71 subgenera, with an estimate of the confidence for placement of species in these taxonomic categories based on whether molecular or radula and/or shell data were used in these determinations. Our proposed classification effectively departs from previous schemes by (1) limiting the number of accepted genera, (2) retaining the majority of species within the genus Conus and (3) assigning members of these genera to species groups/subgenera to enable the effective communication of these groups, all of which we hope will encourage acceptance of this scheme.

Molecular phylogeny and evolution of the cone snails (Gastropoda, Conoidea)

We present a large-scale molecular phylogeny that includes 320 of the 761 recognized valid species of the cone snails (Conus), one of the most diverse groups of marine molluscs, based on three mitochondrial genes (COI, 16S rDNA and 12S rDNA). This is the first phylogeny of the taxon to employ concatenated sequences of several genes, and it includes more than twice as many species as the last published molecular phylogeny of the entire group nearly a decade ago. Most of the numerous molecular phylogenies published during the last 15years are limited to rather small fractions of its species diversity. Bayesian and maximum likelihood analyses are mostly congruent and confirm the presence of three previously reported highly divergent lineages among cone snails, and one identified here using molecular data. About 85% of the species cluster in the single Large Major Clade; the others are divided between the Small Major Clade (∼12%), the Conus californicus lineage (one species), and a newly defined clade (∼3%). We also define several subclades within the Large and Small major clades, but most of their relationships remain poorly supported. To illustrate the usefulness of molecular phylogenies in addressing specific evolutionary questions, we analyse the evolution of the diet, the biogeography and the toxins of cone snails. All cone snails whose feeding biology is known inject venom into large prey animals and swallow them whole. Predation on polychaete worms is inferred as the ancestral state, and diet shifts to molluscs and fishes occurred rarely. The ancestor of cone snails probably originated from the Indo-Pacific; rather few colonisations of other biogeographic provinces have probably occurred. A new classification of the Conidae, based on the molecular phylogeny, is published in an accompanying paper.

Physella acuta: atypical mitochondrial gene order among panpulmonates (Gastropoda)

Mitochondrial (mt) sequences are frequently used for phylogenetic reconstruction and for identification of species of molluscs. This study expands the phylogenetic range of Hygrophila (Panpulmonata) for which such sequence data are available by characterizing the full mt genome of the invasive freshwater snail Physella acuta (Physidae). The mt genome sequences of two P. acuta isolates from Stubblefield Lake, New Mexico, USA, differed in length (14,490 vs 14,314 bp) and showed 11.49% sequence divergence, whereas ITS1 and ITS2 sequences from the nuclear genome differed by 1.75%. The mt gene order of P. acuta (cox1, P, nad6, nad5, nad1, D, F, cox2, Y, W, nad4L, C, Q, atp6, R, E, rrnS, M, T, cox3, I, nad2, K, V, rrnL, L1, A, cytb, G, H, L2, atp8, N, nad2, S1, S2, nad4) differs considerably from the relatively conserved gene order within Panpulmonata. Phylogenetic trees show that the 13 protein-encoding mt gene sequences (equivalent codons) of P. acuta group according to gastropod phylogeny, yet branch lengths and dN/dS ratios for P. acuta indicate elevated amino acid substitutions relative to other gastropods. This study indicates that mt sequences of P. acuta are phylogenetically informative despite a considerable intraspecific divergence and the atypical gene order in its mt genome.

The mitochondrial genome of the venomous cone snail Conus consors

Cone snails are venomous predatory marine neogastropods that belong to the species-rich superfamily of the Conoidea. So far, the mitochondrial genomes of two cone snail species (Conus textile and Conus borgesi) have been described, and these feed on snails and worms, respectively. Here, we report the mitochondrial genome sequence of the fish-hunting cone snail Conus consors and describe a novel putative control region (CR) which seems to be absent in the mitochondrial DNA (mtDNA) of other cone snail species. This possible CR spans about 700 base pairs (bp) and is located between the genes encoding the transfer RNA for phenylalanine (tRNA-Phe, trnF) and cytochrome c oxidase subunit III (cox3). The novel putative CR contains several sequence motifs that suggest a role in mitochondrial replication and transcription.

A very short, functionally constrained sequence diagnoses cone snails in several Conasprella clades

The traditional taxonomy of ca. 700 cone snails assigns all species to a single genus, Conus Linnaeus 1758. However, an increasing body of evidence suggests that some belong to a phylogenetically distinct clade that is sometimes referred to as Conasprella. Previous work (Kraus et al., 2011) showed that a short (259 bp) conserved intronic sequence (CIS) of the γ-glutamyl carboxylase gene (intron 9) can be used to delineate deep phylogenetic relationships among some groups of Conus. The work described here uses intron 9 (338 bp) to resolve problematic relationships among the conasprellans and to distinguish them from Conus proper. Synapomorphic mutations at just 39 sites can resolve several groups within Conasprella because the informative region of intron 9 is so well conserved that the phylogenetic signal is not obscured by homoplasies at conflicting sites. Intron 9 also unambiguously distinguishes Conasprella as a whole from Conus because the conserved regions that are so well conserved within each group are not alignable and clearly not homologous between them. This pattern suggests that expression of the γ-glutamyl carboxylase gene may have undergone a functionally significant change in Conus or Conasprella shortly after they diverged.

A comparative study of nemertean complete mitochondrial genomes, including two new ones for Nectonemertes cf. mirabilis and Zygeupolia rubens, may elucidate the fundamental pattern for the phylum Nemertea

BACKGROUND

The mitochondrial genome is important for studying genome evolution as well as reconstructing the phylogeny of organisms. Complete mitochondrial genome sequences have been reported for more than 2200 metazoans, mainly vertebrates and arthropods. To date, from a total of about 1275 described nemertean species, only three complete and two partial mitochondrial DNA sequences from nemerteans have been published. Here, we report the entire mitochondrial genomes for two more nemertean species: Nectonemertes cf. mirabilis and Zygeupolia rubens.

RESULTS

The sizes of the entire mitochondrial genomes are 15365 bp for N. cf. mirabilis and 15513 bp for Z. rubens. Each circular genome contains 37 genes and an AT-rich non-coding region, and overall nucleotide composition is AT-rich. In both species, there is significant strand asymmetry in the distribution of nucleotides, with the coding strand being richer in T than A and in G than C. The AT-rich non-coding regions of the two genomes have some repeat sequences and stem-loop structures, both of which may be associated with the initiation of replication or transcription. The 22 tRNAs show variable substitution patterns in nemerteans, with higher sequence conservation in genes located on the H strand. Gene arrangement of N. cf. mirabilis is identical to that of Paranemertes cf. peregrina, both of which are Hoplonemertea, while that of Z. rubens is the same as in Lineus viridis, both of which are Heteronemertea. Comparison of the gene arrangements and phylogenomic analysis based on concatenated nucleotide sequences of the 12 mitochondrial protein-coding genes revealed that species with closer relationships share more identical gene blocks.

CONCLUSION

The two new mitochondrial genomes share many features, including gene contents, with other known nemertean mitochondrial genomes. The tRNA families display a composite substitution pathway. Gene order comparison to the proposed ground pattern of Bilateria and some lophotrochozoans suggests that the nemertean ancestral mitochondrial gene order most closely resembles the heteronemertean type. Phylogenetic analysis proposes a sister-group relationship between Hetero- and Hoplonemertea, which supports one of two recent alternative hypotheses of nemertean phylogeny.

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

BACKGROUND

Reconstructing the higher relationships of pulmonate gastropods has been difficult. The use of morphology is problematic due to high homoplasy. Molecular studies have suffered from low taxon sampling. Forty-eight complete mitochondrial genomes are available for gastropods, ten of which are pulmonates. Here are presented the new complete mitochondrial genomes of the ten following species of pulmonates: Salinator rhamphidia (Amphiboloidea); Auriculinella bidentata, Myosotella myosotis, Ovatella vulcani, and Pedipes pedipes (Ellobiidae); Peronia peronii (Onchidiidae); Siphonaria gigas (Siphonariidae); Succinea putris (Stylommatophora); Trimusculus reticulatus (Trimusculidae); and Rhopalocaulis grandidieri (Veronicellidae). Also, 94 new pulmonate-specific primers across the entire mitochondrial genome are provided, which were designed for amplifying entire mitochondrial genomes through short reactions and closing gaps after shotgun sequencing.

RESULTS

The structural features of the 10 new mitochondrial genomes are provided. All genomes share similar gene orders. Phylogenetic analyses were performed including the 10 new genomes and 17 genomes from Genbank (outgroups, opisthobranchs, and other pulmonates). Bayesian Inference and Maximum Likelihood analyses, based on the concatenated amino-acid sequences of the 13 protein-coding genes, produced the same topology. The pulmonates are paraphyletic and basal to the opisthobranchs that are monophyletic at the tip of the tree. Siphonaria, traditionally regarded as a basal pulmonate, is nested within opisthobranchs. Pyramidella, traditionally regarded as a basal (non-euthyneuran) heterobranch, is nested within pulmonates. Several hypotheses are rejected, such as the Systellommatophora, Geophila, and Eupulmonata. The Ellobiidae is polyphyletic, but the false limpet Trimusculus reticulatus is closely related to some ellobiids.

CONCLUSIONS

Despite recent efforts for increasing the taxon sampling in euthyneuran (opisthobranchs and pulmonates) molecular phylogenies, several of the deeper nodes are still uncertain, because of low support values as well as some incongruence between analyses based on complete mitochondrial genomes and those based on individual genes (18S, 28S, 16S, CO1). Additional complete genomes are needed for pulmonates (especially for Williamia, Otina, and Smeagol), as well as basal heterobranchs closely related to euthyneurans. Increasing the number of markers for gastropod (and more broadly mollusk) phylogenetics also is necessary in order to resolve some of the deeper nodes -although clearly not an easy task. Step by step, however, new relationships are being unveiled, such as the close relationships between the false limpet Trimusculus and ellobiids, the nesting of pyramidelloids within pulmonates, and the close relationships of Siphonaria to sacoglossan opisthobranchs. The additional genomes presented here show that some species share an identical mitochondrial gene order due to convergence.

Characterization of the Conus bullatus genome and its venom-duct transcriptome

BACKGROUND

The venomous marine gastropods, cone snails (genus Conus), inject prey with a lethal cocktail of conopeptides, small cysteine-rich peptides, each with a high affinity for its molecular target, generally an ion channel, receptor or transporter. Over the last decade, conopeptides have proven indispensable reagents for the study of vertebrate neurotransmission. Conus bullatus belongs to a clade of Conus species called Textilia, whose pharmacology is still poorly characterized. Thus the genomics analyses presented here provide the first step toward a better understanding the enigmatic Textilia clade.

RESULTS

We have carried out a sequencing survey of the Conus bullatus genome and venom-duct transcriptome. We find that conopeptides are highly expressed within the venom-duct, and describe an in silico pipeline for their discovery and characterization using RNA-seq data. We have also carried out low-coverage shotgun sequencing of the genome, and have used these data to determine its size, genome-wide base composition, simple repeat, and mobile element densities.

CONCLUSIONS

Our results provide the first global view of venom-duct transcription in any cone snail. A notable feature of Conus bullatus venoms is the breadth of A-superfamily peptides expressed in the venom duct, which are unprecedented in their structural diversity. We also find SNP rates within conopeptides are higher compared to the remainder of C. bullatus transcriptome, consistent with the hypothesis that conopeptides are under diversifying selection.

Sessile snails, dynamic genomes: gene rearrangements within the mitochondrial genome of a family of caenogastropod molluscs

BACKGROUND

Widespread sampling of vertebrates, which comprise the majority of published animal mitochondrial genomes, has led to the view that mitochondrial gene rearrangements are relatively rare, and that gene orders are typically stable across major taxonomic groups. In contrast, more limited sampling within the Phylum Mollusca has revealed an unusually high number of gene order arrangements. Here we provide evidence that the lability of the molluscan mitochondrial genome extends to the family level by describing extensive gene order changes that have occurred within the Vermetidae, a family of sessile marine gastropods that radiated from a basal caenogastropod stock during the Cenozoic Era.

RESULTS

Major mitochondrial gene rearrangements have occurred within this family at a scale unexpected for such an evolutionarily young group and unprecedented for any caenogastropod examined to date. We determined the complete mitochondrial genomes of four species (Dendropoma maximum, D. gregarium, Eualetes tulipa, and Thylacodes squamigerus) and the partial mitochondrial genomes of two others (Vermetus erectus and Thylaeodus sp.). Each of the six vermetid gastropods assayed possessed a unique gene order. In addition to the typical mitochondrial genome complement of 37 genes, additional tRNA genes were evident in D. gregarium (trnK) and Thylacodes squamigerus (trnV, trnLUUR). Three pseudogenes and additional tRNAs found within the genome of Thylacodes squamigerus provide evidence of a past duplication event in this taxon. Likewise, high sequence similarities between isoaccepting leucine tRNAs in Thylacodes, Eualetes, and Thylaeodus suggest that tRNA remolding has been rife within this family. While vermetids exhibit gene arrangements diagnostic of this family, they also share arrangements with littorinimorph caenogastropods, with which they have been linked based on sperm morphology and primary sequence-based phylogenies.

CONCLUSIONS

We have uncovered major changes in gene order within a family of caenogastropod molluscs that are indicative of a highly dynamic mitochondrial genome. Studies of mitochondrial genomes at such low taxonomic levels should help to illuminate the dynamics of gene order change, since the telltale vestiges of gene duplication, translocation, and remolding have not yet been erased entirely. Likewise, gene order characters may improve phylogenetic hypotheses at finer taxonomic levels than once anticipated and aid in investigating the conditions under which sequence-based phylogenies lack resolution or prove misleading.

Index-free de novo assembly and deconvolution of mixed mitochondrial genomes

Second-generation sequencing technology has allowed a very large increase in sequencing throughput. In order to make use of this high throughput, we have developed a pipeline for sequencing and de novo assembly of multiple mitochondrial genomes without the costs of indexing. Simulation studies on a mixture of diverse animal mitochondrial genomes showed that mitochondrial genomes could be reassembled from a high coverage of short (35 nt) reads, such as those generated by a second-generation Illumina Genome Analyzer. We then assessed this experimentally with long-range polymerase chain reaction products from mitochondria of a human, a rat, a bird, a frog, an insect, and a mollusc. Comparison with reference genomes was used for deconvolution of the assembled contigs rather than for mapping of sequence reads. As proof of concept, we report the complete mollusc mitochondrial genome of an olive shell (Amalda northlandica). It has a very unusual putative control region, which contains a structure that would probably only be detectable by next-generation sequencing. The general approach has considerable potential, especially when combined with indexed sequencing of different groups of genomes.

Conservation and variation in mitochondrial genomes of gastropods Oncomelania hupensis and Tricula hortensis, intermediate host snails of Schistosoma in China

The complete mitochondrial genomes of intermediate host snails for Schistosoma in China were sequenced, including the sub-species Oncomelania hupensis hupensis in two types, and O. hupensis robertsoni, intermediate hosts for S. japonicum, and Tricula hortensis, the intermediate host of S. sinensium. Four genomes have completely the same gene order as in other caenogastropods, containing 13 protein-coding genes and 22 transfer RNAs. The gene size, start codon and termination codon are mostly the same for all protein-coding genes. However, pairwise sequence alignments revealed quite different degrees of variation. The ribbed-shelled O. hupensis hupensis and the smooth-shelled but with varix O. hupensis hupensis had a lower level of genetic distance (3.1% for protein-coding genes), but the coden usages differed obviously in the mitochondrial genomes of these two types of snails, implying that their genetic difference may be larger than previously recognized. The mean genetic distance between O. hupensis hupensis and O. hupensis robertsoni was 12% for protein-coding genes, indicating a higher degree of genetic difference. In consideration of the difference in morphology and distribution, we considered that O. hupensis hupensis and O. hupensis robertsoni can be considered as separate species. The ribbed-shelled O. hupensishupensis and smooth-shelled O. hupensis robertsoni were phylogenetically clustered together within a same clade, which was then clustered with T. hortensis, confirming their close relationship. However, species or sub-species in the Oncomelania from southeastern Asian countries should be included in future study in order to resolve the phylogenetic relationship and origination of all snails in the genus.

Defining a Clade by Morphological, Molecular and Toxinological Criteria: Distinctive Forms related to Conus praecellens A. Adams, 1854

We carried out a definition of the species group to which Conus praecellens A. Adams 1854 belongs using a combination of comparative morphological data, molecular phylogeny based on standard genetic markers and toxinological markers. Prior to this work, Conus praecellens was generally postulated to belong to a clade of similarly high-spired, smaller Conusspecies such as Conus pagodus Kiener, 1845, Conus memiae (Habe & Kosuge, 1970) and Conus arcuatus Broderip & Sowerby, 1829. The molecular phylogeny and toxinological data demonstrate that these prior hypotheses are incorrect, and that instead, Conus praecellens is in a branch of Conus that includes Conus stupa (Kuroda, 1956), Conus stupella (Kuroda, 1956), Conus acutangulus Lamark, 1810 and surprisingly, some species that are morphologically strikingly different, Conus mitratus Sowerby, 1870 and Conus cylindraceus Broderip & Sowerby, 1830. A more careful analysis of the morphologically diverse forms assigned to Conus praecellens suggests that from the Philippine material alone, there are at least three additional undescribed species, Conus andremenezi, Conus miniexcelsus and Conus rizali. A reevaluation of protoconch/early teleoconch morphology also strongly suggests that Conus excelsus Sowerby III, 1908 is related to these species. Together, the different data suggest a clade including the 10 species above that we designate, the Turriconus (Shikama and Habe, 1968) (clade; there are additional distinctive forms within the clade that may be separable at the species level. The phylogenetic definition using the multidisciplinary approach described herein provides a framework for comprehensively investigating biodiverse lineages of animals, such as the cone snails.

Evolution of Conus peptide toxins: analysis of Conus californicus Reeve, 1844

Conus species are characterized by their hyperdiverse toxins, encoded by a few gene superfamilies. Our phylogenies of the genus, based on mitochondrial genes, confirm previous results that C. californicus is highly divergent from all other species. Genetic and biochemical analysis of their venom peptides comprise the fifteen most abundant conopeptides and over 50 mature cDNA transcripts from the venom duct. Although C. californicus venom retains many of the general properties of other Conus species, they share only half of the toxin gene superfamilies found in other Conus species. Thus, in these two lineages, approximately half of the rapidly diversifying gene superfamilies originated after an early Tertiary split. Such results demonstrate that, unlike endogenously acting gene families, these genes are likely to be significantly more restricted in their phylogenetic distribution. In concordance with the evolutionary distance of C. californicus from other species, there are aspects of prey-capture behavior and prey preferences of this species that diverges significantly from all other Conus.

A molecular phylogenetic framework for the Muricidae, a diverse family of carnivorous gastropods

With over 1600 extant described species, the Muricidae are one of the most species-rich and morphologically diverse families of molluscs. As predators of molluscs, polychaetes, anthozoans barnacles and other invertebrates, they form an important component of many benthic communities. Traditionally, the classification of muricids at specific and generic levels has been based primarily on shells, while subfamilies have been defined largely by radular morphology, although the composition and relationships of suprageneric groups have never been studied exhaustively. Here we present the phylogenetic relationships of 77 muricid species belonging to nine of the ten currently recognized subfamilies, based on Bayesian inference and Maximum Likelihood analyses of partial sequences of three mitochondrial (12S, 16S and COI) and one nuclear (28S) genes. The resulting topologies are discussed with respect to traditional subfamilial arrangements, and previous anatomical and molecular findings. We confirm monophyly of each of the subfamilies Ergalataxinae, Rapaninae, Coralliophilinae, Haustrinae, Ocenebrinae and Typhinae as previously defined, but earlier concepts of Muricinae, Trophoninae and Muricopsinae are shown to be polyphyletic. Based on our phylogenetic hypothesis, a new arrangement of these subfamilies is proposed.

Neogastropod phylogenetic relationships based on entire mitochondrial genomes

Background

The Neogastropoda is a highly diversified group of predatory marine snails (Gastropoda: Caenogastropoda). Traditionally, its monophyly has been widely accepted based on several morphological synapomorphies mostly related with the digestive system. However, recent molecular phylogenetic studies challenged the monophyly of Neogastropoda due to the inclusion of representatives of other caenogastropod lineages (e.g. Littorinimorpha) within the group. Neogastropoda has been classified into up to six superfamilies including Buccinoidea, Muricoidea, Olivoidea, Pseudolivoidea, Conoidea, and Cancellarioidea. Phylogenetic relationships among neogastropod superfamilies remain unresolved.

Results

The complete mitochondrial (mt) genomes of seven Neogastropoda (Bolinus brandaris, Cancellaria cancellata, Conus borgesi, Cymbium olla, Fusiturris similis, Nassarius reticulatus, and Terebra dimidiata) and of the tonnoidean Cymatium parthenopeum (Littorinimorpha), a putative sister group to Neogastropoda, were sequenced. In addition, the partial sequence of the mitochondrial genome of the calyptraeoidean Calyptraea chinensis (Littorinimorpha) was also determined. All sequenced neogastropod mt genomes shared a highly conserved gene order with only two instances of tRNA gene translocation. Phylogenetic relationships of Neogastropoda were inferred based on the 13 mt protein coding genes (both at the amino acid and nucleotide level) of all available caenogastropod mitochondrial genomes. Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic analyses failed to recover the monophyly of Neogastropoda due to the inclusion of the tonnoidean Cymatium parthenopeum within the group. At the superfamily level, all phylogenetic analyses questioned the taxonomic validity of Muricoidea, whereas the monophyly of Conoidea was supported by most phylogenetic analyses, albeit weakly. All analyzed families were recovered as monophyletic except Turridae due to the inclusion of Terebridae. Further phylogenetic analyses based on either a four mt gene data set including two additional Littorinimorpha or combining mt and nuclear sequence data also rejected the monophyly of Neogastropoda but rendered rather unresolved topologies. The phylogenetic performance of each mt gene was evaluated under ML. The total number of resolved internal branches of the reference (whole-mt genome) topology was not recovered in any of the individual gene phylogenetic analysis. The cox2 gene recovered the highest number of congruent internal branches with the reference topology, whereas the combined tRNA genes, cox1, and atp8 showed the lowest phylogenetic performance.

Conclusion

Phylogenetic analyses based on complete mt genome data resolved a higher number of internal branches of the caenogastropod tree than individual mt genes. All performed phylogenetic analyses agreed in rejecting the monophyly of the Neogastropoda due to the inclusion of Littorinimorpha lineages within the group. This result challenges morphological evidence, and prompts for further re-evaluation of neogastropod morphological synapomorphies. The important increase in number of analyzed positions with respect to previous studies was not enough to achieve conclusive results regarding phylogenetic relationships within Neogastropoda. In this regard, sequencing of complete mtDNAs from all closely related caenogastropod lineages is needed. Nevertheless, the rapid radiation at the origin of Neogastropoda may not allow full resolution of this phylogeny based only on mt data, and in parallel more nuclear sequence data will also need to be incorporated into the phylogenetic analyses.

Phylogeny and mitochondrial gene order variation in Lophotrochozoa in the light of new mitogenomic data from Nemertea

BACKGROUND

The new animal phylogeny established several taxa which were not identified by morphological analyses, most prominently the Ecdysozoa (arthropods, roundworms, priapulids and others) and Lophotrochozoa (molluscs, annelids, brachiopods and others). Lophotrochozoan interrelationships are under discussion, e.g. regarding the position of Nemertea (ribbon worms), which were discussed to be sister group to e.g. Mollusca, Brachiozoa or Platyhelminthes. Mitochondrial genomes contributed well with sequence data and gene order characters to the deep metazoan phylogeny debate.

RESULTS

In this study we present the first complete mitochondrial genome record for a member of the Nemertea, Lineus viridis. Except two trnP and trnT, all genes are located on the same strand. While gene order is most similar to that of the brachiopod Terebratulina retusa, sequence based analyses of mitochondrial genes place nemerteans close to molluscs, phoronids and entoprocts without clear preference for one of these taxa as sister group.

CONCLUSION

Almost all recent analyses with large datasets show good support for a taxon comprising Annelida, Mollusca, Brachiopoda, Phoronida and Nemertea. But the relationships among these taxa vary between different studies. The analysis of gene order differences gives evidence for a multiple independent occurrence of a large inversion in the mitochondrial genome of Lophotrochozoa and a re-inversion of the same part in gastropods. We hypothesize that some regions of the genome have a higher chance for intramolecular recombination than others and gene order data have to be analysed carefully to detect convergent rearrangement events.

Remarkable inter- and intra-species complexity of conotoxins revealed by LC/MS

Cone snails have evolved an assortment of venom peptides as an evolutionary strategy for rapid prey immobilization and defence. Earlier studies estimated approximately 100 conopeptides per species. In this study we optimized liquid chromatography and electrospray ionization mass spectrometry for the detection of conopeptides in crude venom to characterize conopeptides present in the venom of individual specimens of Conus textile, C. imperialis and C. marmoreus. Using this approach, we have expanded the predicted number of venom peptides 10-fold to an estimate of 1000-1900 conopeptides per species. Our investigation has also revealed a surprisingly high level of intra-species variation that distinguishes cone snails from other venomous species including spiders and scorpions. Given this inherent diversity and variability, more sensitive bioassays and sequencing techniques will be required to fully explore conotoxin bioactivity.

Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida

Background

Mitochondrial genomes are a valuable source of data for analysing phylogenetic relationships. Besides sequence information, mitochondrial gene order may add phylogenetically useful information, too. Sipuncula are unsegmented marine worms, traditionally placed in their own phylum. Recent molecular and morphological findings suggest a close affinity to the segmented Annelida.

Results

The first complete mitochondrial genome of a member of Sipuncula, Sipunculus nudus, is presented. All 37 genes characteristic for metazoan mtDNA were detected and are encoded on the same strand. The mitochondrial gene order (protein-coding and ribosomal RNA genes) resembles that of annelids, but shows several derivations so far found only in Sipuncula. Sequence based phylogenetic analysis of mitochondrial protein-coding genes results in significant bootstrap support for Annelida sensu lato, combining Annelida together with Sipuncula, Echiura, Pogonophora and Myzostomida.

Conclusion

The mitochondrial sequence data support a close relationship of Annelida and Sipuncula. Also the most parsimonious explanation of changes in gene order favours a derivation from the annelid gene order. These results complement findings from recent phylogenetic analyses of nuclear encoded genes as well as a report of a segmental neural patterning in Sipuncula.