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

Phylogenetic and Comparative Genomics Study of Cephalopina titillator Based on Mitochondrial Genomes

Camel bot fly (Cephalopina titillator) larvae cause myiasis in domesticated and wild camels, resulting in significant economic losses to the camel industry and posing a serious global public health concern. To date, only one mitochondrial genome (mitogenome) of C. titillator isolated from the Alxa Bactrian camel has been reported. Herein, C. titillator was isolated from the Junggar Bactrian camel to assemble a complete circular mitogenome with a length of 16,552 bp encoding 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The mitogenome showed a high A + T content (73.31%), positive AT-skew (0.12), and negative GC-skew (-0.34) base composition patterns. All protein-coding genes (PCGs) employed ATG, ATA, ATT, GTG, or TCG as the start codons and TAA, TAG, or single T as the stop codons. Similar to other parasites in the Oestridae subfamily, the mitogenome was structurally conserved, with genes retaining the same order and direction as those in the ancestral insect mitogenome. The phylogenetic analysis clustered this species with the Oestrinae, showing that the subfamily did not exhibit monophyly. C. titillator isolated from the Junggar Bactrian camel was found to be a sister lineage to that isolated from the Alxa Bactrian camel. Despite the lack of data on the mitogenome of C. titillator isolated from dromedaries in the Middle East, phylogenetic analysis of C. titillator isolated from Xinjiang revealed one distinct lineage of the Xinjiang camel nasal bot fly. In conclusion, this study reports the complete mitogenome of Xinjiang C. titillator for the first time, providing valuable data for future studies on the phylogenetic relationships in this subfamily.

Chemical Synthesis and Insecticidal Activity Research Based on α-Conotoxins

The escalating resistance of agricultural pests to chemical insecticides necessitates the development of novel, efficient, and safe biological insecticides. Conus quercinus, a vermivorous cone snail, yields a crude venom rich in peptides for marine worm predation. This study screened six α-conotoxins with insecticidal potential from a previously constructed transcriptome database of C. quercinus, characterized by two disulfide bonds. These conotoxins were derived via solid-phase peptide synthesis (SPPS) and folded using two-step iodine oxidation for further insecticidal activity validation, such as CCK-8 assay and insect bioassay. The final results confirmed the insecticidal activities of the six α-conotoxins, with Qc1.15 and Qc1.18 exhibiting high insecticidal activity. In addition, structural analysis via homology modeling and functional insights from molecular docking offer a preliminary look into their potential insecticidal mechanisms. In summary, this study provides essential references and foundations for developing novel insecticides.

Novel insights into adaptive evolution based on the unusual AT-skew in Acheilognathus gracilis mitogenome and phylogenetic relationships of bitterling

Acheilognathus gracilis, a bitterling species, distribute in lower reaches of Yangtze River. They are identified as the top-priority bitterling species for conservation as having high evolutionary distinctiveness and are at risk of extinction. In present study, we first sequenced the complete mitogenome of A. gracilis and analyzed its phylogenetic position using 13 PCGs. The A. gracilis mitogenome is 16,774 bp in length, including 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, a control region and the origin of the light strand replication. The overall base composition of A. gracilis in descending order is T 27.9 %, A 27.7 %, C 26.1 % and G 18.3 %, shows a unusual AT-skew with slightly negative. Further investigation revealed A. gracilis uses excess T over A in NADH dehydrogenase 5 (nd5), whereas the most of other bitterlings are biased toward to use A not T, implying there is likely to be unique strategy of adaptive evolution in A. gracilis. We also compared 13 PCGs of 30 bitterling mitogenomes and the results exhibit highly conservative. Phylogenetic trees constructed by 13 PCGs strongly support the monophyly of Acheilognathus and the paraphyly of Rhodeus and Tanakia. Current results will provide valuable information for follow-up research on conservation of species facing with serious population decline and can provide novel insights into the phylogenetic analysis and evolutionary biology research.

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.

Mitogenome Characterization of Four Conus Species and Comparative Analysis

Cone snails, as a type of marine organism, have rich species diversity. Traditionally, classifications of cone snails were based mostly on radula, shell, and anatomical characters. Because of these phenotypic features' high population variability and propensity for local adaptation and convergence, identifying species can be difficult and occasionally inaccurate. In addition, mitochondrial genomes contain high phylogenetic information, so complete mitogenomes have been increasingly employed for inferring molecular phylogeny. To enrich the mitogenomic database of cone snails (Caenogastropoda: Conidae), mitogenomes of four Conus species, i.e., C. imperialis (15,505 bp), C. literatus (15,569 bp), C. virgo (15,594 bp), and C. marmoreus (15,579 bp), were characterized and compared. All 4 of these mitogenomes included 13 protein-coding genes, 2 ribosomal RNA genes, 22 tRNA genes, and non-coding regions. All the Protein Codon Genes (PCGs) of both newly sequenced mitogenomes used TAA or TAG as a terminal codon. Most PCGs used conventional start codon ATG, but an alternative initiation codon GTG was detected in a gene (NADH dehydrogenase subunit 4 (nad4)) of C. imperialis. In addition, the phylogenetic relationships were reconstructed among 20 Conus species on the basis of PCGs, COX1, and the complete mitogenome using both Bayesian Inference (BI) and Maximum Likelihood (ML). The phylogenetic results supported that C. litteratus, C. quercinus, and C. virgo were clustered together as a sister group (PP = 1, BS = 99), but they did not support the phylogenetic relation of C. imperialis and C. tribblei (PP = 0.79, BS = 50). In addition, our study established that PCGs and complete mitogenome are the two useful markers for phylogenetic inference of Conus species. These results enriched the data of the cone snail's mitochondrion in the South China Sea and provided a reliable basis for the interpretation of the phylogenetic relationship of the cone snail based on the mitochondrial genome.

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 genome of an endemic cichlid Etroplus canarensis from Western Ghats, India (Perciformes: Cichlidae) and molecular phylogenetic analysis

BACKGROUND

The Indian endemic cichlid Etroplus canarensis (Canara pearl spot) is an endangered fish and is one among the three Indian cichlids (Etroplinae) that had a restricted distribution in the South Canara region of Karnataka, India. Despite considerable investigations, the phylogeny of Indian Cichlids and its biogeographical origin is still ambiguous and remains a question under discussion which is scrutinized based on whole mitogenomes in the present study.

METHODS AND RESULTS

We report the 16,339 bp complete mitochondrial genome of E. canarensis for the first time using the next-generation sequencing methods. Comparison of gene arrangement and genome characterization was found to commensurate with the previous reports on two Indian cichlid fishes, E. suratensis and E. maculatus. ND6 has been identified as a gene with the highest evolutionary rate and COI and COII is the most conserved gene based on p-genetic distance calculation. Substitution rate (ka/ks) was found to be very low indicating a reduced rate of evolution among subfamily Etroplinae accounting for its subsided species divergence of Indian cichlids. Phylogenetic analysis of Indian cichlids based on a combined dataset of 12 protein-coding genes representing cichlids generated high posterior probability values pillaring paraphyletic nature of Indian Malagasy lineage and monophyletic Indian genus Etroplus.

CONCLUSION

The mitogenome sequence of E. canarensis may provide fundamental molecular data useful for further researches on genetic diversity, endemicity and the conservation of this endangered freshwater fish.

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.

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

The first Conus genome assembly reveals a primary genetic central dogma of conopeptides in C. betulinus

Although there are various Conus species with publicly available transcriptome and proteome data, no genome assembly has been reported yet. Here, using Chinese tubular cone snail (C. betulinus) as a representative, we sequenced and assembled the first Conus genome with original identification of 133 genome-widely distributed conopeptide genes. After integration of our genomics, transcriptomics, and peptidomics data in the same species, we established a primary genetic central dogma of diverse conopeptides, assuming a rough number ratio of ~1:1:1:10s for the total genes: transcripts: proteins: post-translationally modified peptides. This ratio may be special for this worm-hunting Conus species, due to the high diversity of various Conus genomes and the big number ranges of conopeptide genes, transcripts, and peptides in previous reports of diverse Conus species. Only a fraction (45.9%) of the identified conotopeptide genes from our achieved genome assembly are transcribed with transcriptomic evidence, and few genes individually correspond to multiple transcripts possibly due to intraspecies or mutation-based variances. Variable peptide processing at the proteomic level, generating a big diversity of venom conopeptides with alternative cleavage sites, post-translational modifications, and N-/C-terminal truncations, may explain how the 133 genes and ~123 transcripts can generate thousands of conopeptides in the venom of individual C. betulinus. We also predicted many conopeptides with high stereostructural similarities to the putative analgesic ω-MVIIA, addiction therapy AuIB and insecticide ImI, suggesting that our current genome assembly for C. betulinus is a valuable genetic resource for high-throughput prediction and development of potential pharmaceuticals.

Complete Mitogenome of a Leaf-Mining Buprestid Beetle, Trachys auricollis, and Its Phylogenetic Implications

A complete mitogenome of Trachys auricollis is reported, and a mitogenome-based phylogenetic tree of Elateriformia with all protein-coding genes (PCGs), rRNAs, and tRNAs is presented for the first time. The complete mitochondrial genome of T. auricollis is 16,429 bp in size and contains 13 PCGs, two rRNA genes, 22 tRNA genes, and an A + T-rich region. The A + T content of the entire genome is approximately 71.1%, and the AT skew and GC skew are 0.10 and -0.20, respectively. According to the the nonsynonymous substitution rate to synonymous substitution rates (Ka/Ks) of all PCGs, the highest and lowest evolutionary rates were observed for atp8 and cox1, respectively, which is a common finding among animals. The start codons of all PCGs are the typical ATN. Ten PCGs have complete stop codons, but three have incomplete stop codons with T or TA. As calculated based on the relative synonymous codon usage (RSCU) values, UUA(L) is the codon with the highest frequency. Except for trnS1, all 22 tRNA genes exhibit typical cloverleaf structures. The A + T-rich region of T. auricollis is located between rrnS and the trnI-trnG-trnM gene cluster, with six 72-bp tandem repeats. Both maximum likelihood (ML) and Bayesian (BI) trees suggest that Buprestoidea is close to Byrrhoidea and that Buprestoidea and Byrrhoidea are sister groups of Elateroidea, but the position of Psephenidae is undetermined. The inclusion of tRNAs might help to resolve the phylogeny of Coleoptera.

High-Throughput Identification and Analysis of Novel Conotoxins from Three Vermivorous Cone Snails by Transcriptome Sequencing

The venom of each Conus species consists of a diverse array of neurophysiologically active peptides, which are mostly unique to the examined species. In this study, we performed high-throughput transcriptome sequencing to extract and analyze putative conotoxin transcripts from the venom ducts of 3 vermivorous cone snails (C. caracteristicus, C. generalis, and C. quercinus), which are resident in offshore waters of the South China Sea. In total, 118, 61, and 48 putative conotoxins (across 22 superfamilies) were identified from the 3 Conus species, respectively; most of them are novel, and some possess new cysteine patterns. Interestingly, a series of 45 unassigned conotoxins presented with a new framework of C-C-C-C-C-C, and their mature regions were sufficiently distinct from any other known conotoxins, most likely representing a new superfamily. O- and M-superfamily conotoxins were the most abundant in transcript number and transcription level, suggesting their critical roles in the venom functions of these vermivorous cone snails. In addition, we identified numerous functional proteins with potential involvement in the biosynthesis, modification, and delivery process of conotoxins, which may shed light on the fundamental mechanisms for the generation of these important conotoxins within the venom duct of cone snails.

High Throughput Identification of Novel Conotoxins from the Vermivorous Oak Cone Snail (Conus quercinus) by Transcriptome Sequencing

The primary objective of this study was to realize the large-scale discovery of conotoxin sequences from different organs (including the venom duct, venom bulb and salivary gland) of the vermivorous Oak cone snail, Conus quercinus. Using high-throughput transcriptome sequencing, we identified 133 putative conotoxins that belong to 34 known superfamilies, of which nine were previously reported while the remaining 124 were novel conotoxins, with 17 in new and unassigned conotoxin groups. A-, O₁-, M-, and I₂- superfamilies were the most abundant, and the cysteine frameworks XIII and VIII were observed for the first time in the A- and I₂-superfamilies. The transcriptome data from the venom duct, venom bulb and salivary gland showed considerable inter-organizational variations. Each organ had many exclusive conotoxins, and only seven of all the inferred mature peptides were common in the three organs. As expected, most of the identified conotoxins were synthesized in the venom duct at relatively high levels; however, a number of conotoxins were also identified in the venom bulb and the salivary gland with very low transcription levels. Therefore, various organs have different conotoxins with high diversity, suggesting greater contributions from several organs to the high-throughput discovery of new conotoxins for future drug development.

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.