Next Article in Journal
Comprehensive Analysis of the Catalase (CAT) Gene Family and Expression Patterns in Rubber Tree (Hevea brasiliensis) under Various Abiotic Stresses and Multiple Hormone Treatments
Next Article in Special Issue
Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)
Previous Article in Journal
Force-Induced Alignment of Nanofibrillated Bacterial Cellulose for the Enhancement of Cellulose Composite Macrofibers
Previous Article in Special Issue
The Complete Mitochondrial Genomes of Two Rock Scallops (Bivalvia: Spondylidae) Indicate Extensive Gene Rearrangements and Adaptive Evolution Compared with Pectinidae
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 1
10.3390/ijms25010067
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

The Complete Mitochondrial Genome and Phylogenetic Analysis of the Freshwater Shellfish Novaculina chinensis (Bivalvia: Pharidae)

by
Ziquan Zhou
,
Yuxin Song
,
Zewen Zheng
,
Yunguang Liu
,
Haiyan Yao
,
Xiaozhen Rao
* and
Gang Lin
*
Fujian Key Laboratory of Special Marine Bioresource Sustainable Utilization, College of Life Sciences, Fujian Normal University, Fuzhou 350117, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(1), 67; https://doi.org/10.3390/ijms25010067
Submission received: 15 November 2023 / Revised: 13 December 2023 / Accepted: 18 December 2023 / Published: 20 December 2023
(This article belongs to the Special Issue Mitochondrial Genome of Aquatic Animals: Analysis of Structure, Evolution and Diversity)
Browse Figures
Versions Notes

Abstract

:
Razor clams, belonging to the Pharidae and Solenidae families, are ecologically and economically important; however, very little research has been conducted on the Pharidae family. The genus Novaculina is a marine-derived freshwater lineage, and Novaculina chinensis is a rare freshwater species of the Pharidae family. In order to understand the phylogenetic relationships of N. chinensis, we sequenced the mitochondrial genome of the genus Novaculina, which is 16,262 bp in length and consists of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The phylogenetic relationships of 69 Imparidentian mitochondrial genomes (mitogenomes) indicated that N. chineisis is closely related to Sinonovacula constricta of the order Adapedonta. Our study also found that the Ka/Ks ratios of 12 protein-coding genes in the Pharidae family are lower than one, indicating the occurrence of negative purification selection. Morphological observations of the siphons of N. chinensis, Novaculina myanmarensis, and Novaculina gangetica indicate that N. chinensis may be the ancestral clade of the genus Novaculina, which has not been proposed in previous studies. Our study provides useful molecular information on the phylogenetic and evolutionary relationships of Pharidae and also contributes to the conservation and management of the germplasm resources of N. chinensis.

1. Introduction

Bivalves have a wide range of taxa in marine, brackish, and freshwater. Razor clams are important deep-burrowing bivalves, most of which are found in the shallow waters of the tropical, subtropical, and temperate seas, and belong to two families (Pharidae and Solenidae) in the superorder Imparidentia [1,2]. Although Pharidae is mainly composed of marine species, it includes one genus that represents four freshwater species that are geographically separated from each other [3]. These four species are distributed throughout various freshwater drainages in Asia, from the Ganges River located in India to the Yangtze River situated in China (Figure 1) [4].
Novaculina chinensis was first reported in Taihu and Gaoyou Lake, Jiangsu Province, China [5]. It is mainly distributed in the middle and lower reaches of the Yangtze River in China and the Minjiang River in Fuzhou, China [6]. Due to overfishing, water pollution, and habitat changes, the population of N. chinensis has declined significantly in all regions and its germplasm resources are in danger of being depleted [7]. Previous studies of N. chinensis have mainly focused on its reproductive cycle, morphology, nutritive compositions, and the ultrastructure of its sex gonads [7,8,9,10,11]. The research on the genus Novaculina located in Southeast Asia has mainly focused on its biodiversity and biogeography [3,4]. Prior to this work, there was little research at the molecular level for this genus. Only the cox1, 16S rRNA (rrnL), and 28S rRNA (rrnS) genes of N. gangetica and N. myanmarensis were sequenced.
The superorder Imparidentia is a newly defined branch of bivalves encompassing diverse clades of marine, brackish, and freshwater bivalve mollusks [12]. It includes five major lineages: Lucinida, Cardiida, Adapedonta, Myida, and Venerida [2]. However, the phylogenetic relationship of N. chinensis in Imparidentia is not clear due to it not being analyzed before. The mitochondrial genome (mitogenome) of classic metazoans is a small and compact circular molecule that typically includes 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs), with the exception of bivalve shellfish in which loss of the atp8 gene is observed in some species [13,14]. The complete mitochondrial genome has been widely used for studying the phylogenetic reconstruction and adaptative evolution due to its maternal inheritance, low intermolecular recombination, high copy number, and high replacement rate [2,15,16].
In this study, we sequenced and described the complete mitochondrial genome of N. chinensis, with the aim of analyzing the genomic features of its mitogenome, including the genome structure, nucleotide composition, and codon usage, as well as the selection pressure in the Pharidae family. Moreover, we constructed a phylogenetic tree to infer the phylogenetic position of N. chinensis in the superorder Imparidentia. Overall, our study provides useful molecular data to better understand the phylogenetic relationship and evolutionary journey of the genus Novaculina, which is important for the conservation and management of the shellfish germplasm resources.

2. Results and Discussion

2.1. Genome Features

The complete mitogenome of N. chinensis is 16,262 bp in length (Figure 2). It has 12 PCGs (except atp8), 22 tRNA genes, and 2 rRNA genes. All of the genes were encoded on the heavy strand (Table 1). The base composition was as follows: A (28.15%), T (43.45%), G (18.86%), and C (9.54%). The A + T content (71.60%) was higher than the G + C content (28.40%), indicating a significant AT bias (Table 2). The AT skew value was −0.214, while the GC skew value was 0.328. Four overlaps were detected in the mitochondrial genome, and among these, the largest overlap was found between trnE and trnS2. The length of N. chinensis was smaller than other species in the Pharidae family, but the AT content was higher than other species in the Pharidae family (Table S1).
In the present study, the mitochondrial genome structure of N. chinensis was found to be largely consistent with that of other published species in the Pharidae family, showing a high degree of conservation. The mitogenomes of Siliqua minima, Cultellus attenuatus, S. constricta, and N. chinensis had the same genome organization, containing 12 PCGs (except atp8), 22 tRNAs, and 2 rRNAs, respectively; they have the same gene order without gene rearrangement and four mitogenomes varied in size from 16,262 bp to 17,225 bp [15,17,18]. The mitogenome of N. chinensis has a different mitochondrial genome size and nucleotide composition than other species in the Pharidae family. The small genome size and high AT content is a distinctive feature of N. chinensis. Wu et al. [19] and Annam et al. [20] sequenced the mitogenomes of five freshwater mussels in the Unionidae family; they had the same features of mitogenome including 13 PCGs, 22 tRNAs, 2 rRNAs, and one female specific gene (FORF). Among the 38 mitochondrial genes, 11 genes were encoded on the heavy chain and the remaining 27 genes were encoded on the light chain; however, all the mitochondrial genes of N. chinensis were encoded on the heavy strand, yet this phenomenon generally occurs in the mitochondrial genomes of marine bivalves [13,21]. Recently, multi-locus phylogenetic analyses have supported the genus Novaculina as a relict marine freshwater lineage [3,22]. Therefore, our findings have indicated that the mitogenome of N. chinensis retains some of the characteristics of marine bivalves.

2.2. Protein-Coding Genes, Transfer RNA, and Ribosomal RNA Genes

The mitogenome of N. chinensis has 12 PCGs and lacks the atp8 gene. The total length of all 12 PCGs of N. chinensis is 11,502 bp, accounting for 70.73% of the complete length of the mitogenome (Table 2). For all 12 PCGs identified in the N. chinensis mitogenome, three genes (cox1, nad5, and nad6) were initiated with the start codon ATT, and the remaining nine genes had the start codon ATG. The cox2, nad4, nad3, nad6, and cox3 genes carried the termination codon TAG (Table 1). Moreover, the most common termination codon TAA was detected in seven PCGs.
The absence of the atp8 gene has been suggested for most bivalves, especially marine species with the exception the venerid Venerupis philippinarum, the hiatellid Hiatella arctica, and unionid species [15,23]. The absence of atp8 is also observed in the mitogenomes of three species in the Pharidae family: S. minima, C. attenuatus, and S. constricta [15,17,18]. ATP synthase is the final enzymatic complex in the respiratory chain, directly producing ATP as it couples with the electrochemical gradient of the inner mitochondrial membrane [24,25]. Sun et al. [24] suggest that the atp8 genes might be relaxed from selective constraints because of the change in locomotive ability and the reduced energy requirements that emerged in marine bivalve evolution. The lack of the atp8 gene in deep-burrowing bivalves of the Pharidae family may be due to its weak motility and low energy requirements.
Most invertebrates utilized two to four codons to encode its amino acids [26]. The nucleotide relative synonymous codon usages (RSCUs) of N. chinensis are presented (Figure 3, Table 3). UUA (Leu2), GCU (Ala), CCU (Pro), and UCU (Ser) are the most frequently used codons, whereas CUC, CUG (Leu1), AUC (Ile), and GUC (Val) are relatively scarce. As per the RSCU values, codons ending with an A or T were preferred. The two most frequent amino acids in the PCGs of N. chinensis were Ser and Leu.
Similar to most bivalves, the mitochondrial genome of N. chinensis contained 22 transfer RNA genes and 2 ribosomal RNA genes. The size of the 22 tRNA genes varied from 63 to 68 bp, and all of these genes can be classified into typical secondary structures (except two trnS genes). The two ribosomal RNA genes included rrnL and rrnS, with the former having a length of 1232 bp (between nad6 and atp6), and the latter standing at 849 bp (between trnM and cox3).

2.3. KaKs Analysis

The calculation of non-synonymous substitutions (Ka) and synonymous substitutions (Ks) is crucial for constructing a phylogenetic tree and understanding the evolutionary dynamics of protein-coding genes (PCGs) in closely related species [27,28]. The Ka/Ks ratio is used to determine whether selective pressure acted on PCGs during evolution: Ka/Ks > 1, positive selection; Ka/Ks = 1, neutral selection; and Ka/Ks < 1, negative selection [29,30].
To analyze the selection pressure on mitochondrial PCGs of the species in the Pharidae family, the Ka/Ks ratio of 12 PCGs were calculated. The Ka/Ks ratios of 12 PCGs are less than one (ranging from 0.0640 to 0.2311) (Table S2), indicating that those mitochondrial PCGs were under strong negative or purifying selection. The mitochondrial cox1 gene had the smallest Ka/Ks value and bears the largest purifying selection pressure, while the highest Ka/Ks value is nad6. In most metazoans, the mitochondrial genes coding for the cytochrome c oxidase and cytochrome b are more conserved than the genes coding for NADH dehydrogenase [31]. Sun et al. [24] suggested that less motile bivalves survive and reproduce with lower metabolic efficiency, which may have accumulated more non-synonymous mutations in their mitochondrial genomes. The species in the Pharidae family have lower energy demands due to the fact that they are deep-burrowing bivalves. Therefore, we suggest that the species in the Pharidae family have also accumulated more non-synonymous mutations. However, there is relatively little useful molecular information in the Pharidae family and more research is required to understand the evolution of mitochondrial genes.

2.4. Analysis of Phylogenetic

To explore the phylogenetic implications of the N. chinensis mitogenome in Imparidentia, we reconstructed a phylogenetic tree. Using Bayesian inference (BI) and maximum likelihood (ML) analyses of twelve protein-coding genes (except atp8) from sixty-nine species, we obtained nearly identical topologies with both methods (Figure 4), with high support in most nodes. Through the prism of phylogenetic analysis, we can observe the orders of Imparidentia, comprised Venerida, Cardiida, Adapedonta, and Lucinida, which is consistent with the phylogenetic trees constructed by Wang et al. [2] and Lemer et al. [32]; the latter of which utilized transcriptional data and morpho-anatomical features. Additionally, this result is largely consistent with the mitogenome-based phylogenetic tree constructed by Feng et al. [15]. Furthermore, the BI and ML analyses support the monophyly of the families Pharidae, Solenidae, and Hiatellidae, with each forming a distinct assemblage. Both analyses confirmed the sister group relationship of Pharidae and Solenidae, whereas Hiatellidae was identified as the sister to Pharidae and Solenidae. Notably, our findings indicate that N. chinensis and S. constricta in the Pharidae family constitute a new branch that is closely related to C. ttenuates, E. leei, and S. minima within the same family. The phylogenetic relationships between 11 species in the order Adapedonta are (Panopea globose + (Panopea abrupta + Panopea generosa)) + Hiatella arctica + (Solen grandis + Solen strictus) + (S. minima + E. leei + C. attenuatus + (S. constricta + N. chinensis)).
Siphons play a crucial role in bivalve classification [33,34,35]. According to recent studies [3], we find that N. chinensis differs significantly in terms of morphology from N. myanmarensis and N. gangetica, which are distributed in Southeast Asia. Our observations indicate that N. chinensis has long, fused siphons, while the N. myanmarensis and N. gangetica siphons are long and separate (Figure 5 and Figure 6). We did not make comparisons with N. siamensis due to insufficient morphological information. Siphon size is one of the main factors determining the burying depth of benthic bivalves and thus plays a critical role in their survival [36]. Novaculina chinensis, N. myanmarensis, and N. gangetica all have long siphons, the only difference being whether the siphons are separated or not. Bivalves with long, separate siphons that remain active and extensive are better suited for survival than long, fused siphons because they can extend from the bivalve far above the surrounding surface to pick up food [37]. Long, fused siphons are more primitive from an evolutionary standpoint. Therefore, we hypothesize that N. chinensis may be the ancestral lineage of N. myanmarensis. and N. gangetica. However, to understand this process of evolution, additional molecular information on the genus Novaculina is necessary.

3. Materials and Methods

3.1. Sample Collection and DNA Extraction

The samples were collected from Taojiang River, Fuzhou city, Fujian Province. Morphological traits of the specimen were identified based on the following available sources from the literature [5,9]. Samples were photographed using a dissecting microscope. Total genomic DNA was extracted from the adductor muscle tissue using the DNeasy tissue kit (Qiagen, Beijing, China) and the procedure was conducted in accordance with the manufacturer’s protocols. The animal study protocol was approved by the Animal Ethics and Welfare Committee of Fujian Normal University (Approval No. IACUC-20230030).

3.2. Sequencing and Genome Assembly

The purified DNA was fragmented to ~500 bp using the Covaris M220 system, used to construct short-insert libraries according to the manufacturer’s instructions (TruSeqTM Nano DNA Sample Prep Kit, Illumina, Shanghai, China), and was then sequenced on an Illumina NovaSeq 6000 platform (BIOZERON Co., Ltd., Shanghai, China) with a 150 bp paired-end reads length, where raw reads were filtered using Trimmomatic 0.39 [38]. The mitogenome was generated via de novo assembly using GetOrganelle v1.7.5 as a reference for the mitochondrial genomes of closely related species (GeneBank No. EU880278.1 for S. constricta and MW653805.1 for C. attenuatus).

3.3. Annotation and Sequence Analysis of Mitogenome

The mitochondrial genes were annotated using the online MITOS tool [39]. Default parameters were applied to predict protein-coding genes (PCGs), transfer RNA (tRNA) genes, and ribosomal RNA (rRNA) genes. The position of each coding gene was determined using BLAST searches against reference mitochondrial genes. Manual corrections of the start/stop codons of the genes were performed in SnapGene Viewer by referencing the reference mitogenomes. The circular mitogenome map of N. chinensis was drawn using Proksee (CG view) [40]. Relative synonymous codon usage (RSCU) of 12 PCGs was calculated using Phylosuite v1.2.3 [41]. The skew values were calculated according to the following formulas: AT skew = (A − T)/(A + T) and GC skew = (G − C)/(G + C) [42]. We used DnaSP [43] to calculate Ka/Ks ratios for 12 PCGs in the Pharidae family.

3.4. Phylogenetic Analysis

To investigate the phylogenetic relationship of N. chinensis., mitogenome sequences of 69 species in Imparidentia were downloaded from the NCBI (Table 1). Mitogenomes of Mimachlamys nobilis and Azumapecten farreri were used as outgroups. Since most bivalve species lack the gene atp8, the phylogenetic analysis was conducted using amino acid sequences of 12 PCGs. Statistics for basic characteristics and the extraction of sequences were executed using PhyloSuite v1.2.3 [41]. A total of 12 PCGs were aligned using MAFFT v7.490 [44]. Ambiguously aligned fragments of 12 alignments were removed in batches using Gblocks 0.91b [45]. Then, the aligned nucleotide sequences were concatenated. ModelFinder v2.2.0 [46] was used to select the best-fit partition model (edge-linked) using a BIC criterion.
Phylogenies were deduced using Bayesian inference (BI) and maximum likelihood (ML) methods. Maximum likelihood phylogenies were inferred using IQ-TREE v2.2.0 [47] within the Edge-linked partition model for 10,000 ultrafast bootstraps. Bayesian Inference phylogenies were inferred using MrBayes v3.2.7a [48] under a partition model (2 parallel runs, 2,000,000 generations), in which the initial 25% of the sampled data were discarded as burn-in. When the value of the average standard deviation of split frequencies drops below 0.01, it serves as an indication that the BI run has converged. Phylogenetic trees were visualized and annotated using the Interactive Tree of Life (ITOL) (https://itol.embl.de/itol.cgi (accessed on 16 August 2023)) [49].

4. Conclusions

In this study, we sequenced the complete mitochondrial genome of N. chinensis from the genus Novaculina. Its mitogenome contains 12 PCGs (except atp8), 22 tRNAs, and 2 rRNAs. Compared to other marine species in the Pharidae family, it exhibits a high AT content, small mitogenome size, and all genes encoded on the heavy strand. Novaculina chinensis diverges from species found in the freshwater basins of South and Southeast Asia, mainly in terms of the siphon morphology. More mitochondrial genomes from the species of the Pharidae family are necessary to explain this phenomenon and understand the family’s evolutionary progress. These findings also aid in conserving N. chinensis germplasm resources and studying its biological characteristics, providing vital molecular information for researching razor clams, which possess significant ecological and economic importance.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25010067/s1.

Author Contributions

Conceptualization, G.L. and X.R.; methodology, Z.Z. (Ziquan Zhou); software, Z.Z. (Ziquan Zhou); validation, Z.Z. (Ziquan Zhou); formal analysis, Z.Z. (Ziquan Zhou); investigation, X.R. and G.L.; resources, Y.S., Y.L., Z.Z. (Zewen Zheng), and H.Y.; data curation, Y.S.; writing—original draft preparation, Z.Z. (Ziquan Zhou); writing—review and editing, Z.Z. (Ziquan Zhou) and G.L.; visualization, Z.Z. (Ziquan Zhou); supervision, X.R. and G.L.; project administration, X.R. and G.L.; funding acquisition, X.R. and G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Fujian Province (2023J011575) and supported by Water Conservancy Bureau of Minhou Government, Fuzhou City (DH-1911).

Institutional Review Board Statement

The animal study protocol was approved by the Animal Ethics and Welfare Committee of Fujian Normal University (Approval No. IACUC-20230030).

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequenced mitogenome in this study has been deposited in GenBank (NC_077598.1).

Acknowledgments

We are grateful to Fuzhou Institute of Oceanography for providing financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saeedi, H.; Costello, M.J. The Biology, Ecology, and Societal Importance of Razor Clams. In Encyclopedia of the World’s Biomes; Goldstein, M.I., DellaSala, D.A., Eds.; Elsevier: Oxford, UK, 2020; pp. 494–498. [Google Scholar]
  2. Wang, Y.; Yang, Y.; Kong, L.; Sasaki, T.; Li, Q. Phylogenomic Resolution of Imparidentia (Mollusca: Bivalvia) Diversification through Mitochondrial Genomes. Mar. Life Sci. Technol. 2023, 5, 326–336. [Google Scholar] [CrossRef] [PubMed]
  3. Bolotov, I.N.; Vikhrev, I.V.; Lopes-Lima, M.; Lunn, Z.; Chan, N.; Win, T.; Aksenova, O.V.; Gofarov, M.Y.; Kondakov, A.V.; Konopleva, E.S.; et al. Discovery of Novaculina myanmarensis Sp. Nov. (Bivalvia: Pharidae: Pharellinae) Closes the Freshwater Razor Clams Range Disjunction in Southeast Asia. Sci. Rep. 2018, 8, 16325. [Google Scholar] [CrossRef] [PubMed]
  4. Win, T.; Bolotov, I.N.; Vikhrev, I.V.; Lunn, Z.; Chan, N.; Konopleva, E.S.; Gofarov, M.Y.; Tomilova, A.A.; Kondakov, A.V. Phylogeography and Distribution of the Freshwater Razor Clams Novaculina myanmarensis and N. Gangetica in Myanmar, with Notes on Two Doubtful Nominal Taxa Described as Novaculina Members (Bivalvia: Pharidae). Ecol. Monten. 2021, 40, 59–67. [Google Scholar] [CrossRef]
  5. Liu, Y.Y.; Zhang, W.Z. A New Species of Freshwater Razor Clam, Novaculina chinensis, from Jiangsu Province, China. Acta Zootaxonomica Sin. 1979, 4, 356–357. [Google Scholar]
  6. Shu, F.; Zhu, Q.; Zhang, N.; Zhang, X.; Xie, S. Novaculina chinensis Found in Lake Weishan, Shandong Province. Chin. J. Zool. 2013, 48, 278–280. [Google Scholar] [CrossRef]
  7. Rao, X.; Xu, Y.; Chen, Y.; Lin, G. A Study on Induced Spawning of Novaculina chinensis. J. Fujian Teach. Univ. Nat. Sci. Ed. 2003, 78–81+93. [Google Scholar]
  8. Jiang, J.; Chen, Y.; Xu, Y.; Huang, J.; Yan, G. A Study on the Seasonal Changes of Novaculina chinensis in Taojiang. J. Fujian Teach. Univ. Nat. Sci. Ed. 1998, 14, 86–91. [Google Scholar]
  9. Rao, X.; Chen, Y.; Xu, Y. Study on the Morphology of Novaculina Chinensis. J. -Fujian Teach. Univ. Nat. Sci. Ed. 1998, 14, 71–75. [Google Scholar]
  10. Xiaozhen, R.; Yinshan, C.; Wenlie, C.; Youqin, X.; Lianyun, C. Ultrastructural Studies on Spermatogenesis of Novaculina chinensis. Chin. J. Zool. 2000, 35, 2–5. [Google Scholar]
  11. Rao, X. Ultrastructures of Novaculina chinensis during Oogenesis. Chin. J. Appl. Environ. Biol. 2001, 7, 248–253. [Google Scholar]
  12. Bieler, R.; Mikkelsen, P.M.; Collins, T.M.; Glover, E.A.; González, V.L.; Graf, D.L.; Harper, E.M.; Healy, J.; Kawauchi, G.Y.; Sharma, P.P. Investigating the Bivalve Tree of Life–an Exemplar-Based Approach Combining Molecular and Novel Morphological Characters. Invertebr. Syst. 2014, 28, 32–115. [Google Scholar] [CrossRef]
  13. Ghiselli, F.; Gomes-dos-Santos, A.; Adema, C.M.; Lopes-Lima, M.; Sharbrough, J.; Boore, J.L. Molluscan Mitochondrial Genomes Break the Rules. Philos. Trans. R. Soc. B 2021, 376, 20200159. [Google Scholar] [CrossRef]
  14. Lavrov, D.V.; Pett, W. Animal Mitochondrial DNA as We Do Not Know It: Mt-Genome Organization and Evolution in Nonbilaterian Lineages. Genome Biol. Evol. 2016, 8, 2896–2913. [Google Scholar] [CrossRef] [PubMed]
  15. Feng, J.; Guo, Y.; Yan, C.; Ye, Y.; Yan, X.; Li, J.; Xu, K.; Guo, B.; Lü, Z. Novel Gene Rearrangement in the Mitochondrial Genome of Siliqua minima (Bivalvia, Adapedonta) and Phylogenetic Implications for Imparidentia. PLoS ONE 2021, 16, e0249446. [Google Scholar] [CrossRef] [PubMed]
  16. Luo, L.; Xu, Y.; Wang, S.; Zhang, R.; Guo, K.; Xu, W.; Zhao, Z. Complete Mitochondrial Genome Sequence and Phylogenetic Analysis of Procambarus clarkii and Cambaroides dauricus from China. Int. J. Mol. Sci. 2023, 24, 11282. [Google Scholar] [CrossRef] [PubMed]
  17. Li, H.; Yu, R.; Ma, P.; Li, C. Complete Mitochondrial Genome of Cultellus attenuatus and Its Phylogenetic Implications. Mol. Biol. Rep. 2022, 49, 8163–8168. [Google Scholar] [CrossRef] [PubMed]
  18. Zheng, R.; Li, J.; Niu, D. The Complete DNA Sequence of the Mitochondrial Genome of Sinonovacula constricta (Bivalvia: Solecurtidae). Acta Oceanol. Sin. 2010, 29, 88–92. [Google Scholar] [CrossRef]
  19. Wu, R.-W.; Liu, X.-J.; Wang, S.; Roe, K.J.; Ouyang, S.; Wu, X.-P. Analysis of Mitochondrial Genomes Resolves the Phylogenetic Position of Chinese Freshwater Mussels (Bivalvia, Unionidae). ZooKeys 2019, 812, 23–46. [Google Scholar] [CrossRef]
  20. Pavan-Kumar, A.; Varshney, S.; Suman, S.; Das, R.; Chaudhari, A.; Krishna, G. Complete Mitochondrial Genome of Freshwater Pearl Mussel Lamellidens marginalis (Lamarck, 1819) and Its Phylogenetic Relation within Unionidae Family. Mol. Biol. Rep. 2022, 49, 9593–9603. [Google Scholar] [CrossRef]
  21. Xu, M.; Gu, Z.; Huang, J.; Guo, B.; Jiang, L.; Xu, K.; Ye, Y.; Li, J. The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae. Genes 2023, 14, 910. [Google Scholar] [CrossRef]
  22. Bolotov, I.N.; Aksenova, O.V.; Bakken, T.; Glasby, C.J.; Gofarov, M.Y.; Kondakov, A.V.; Konopleva, E.S.; Lopes-Lima, M.; Lyubas, A.A.; Wang, Y.; et al. Discovery of a Silicate Rock-Boring Organism and Macrobioerosion in Fresh Water. Nat. Commun. 2018, 9, 2882. [Google Scholar] [CrossRef] [PubMed]
  23. Breton, S.; Stewart, D.T.; Hoeh, W.R. Characterization of a Mitochondrial ORF from the Gender-Associated mtDNAs of Mytilus Spp. (Bivalvia: Mytilidae): Identification of the “Missing” ATPase 8 Gene. Mar. Genom. 2010, 3, 11–18. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, S.; Li, Q.; Kong, L.; Yu, H. Limited Locomotive Ability Relaxed Selective Constraints on Molluscs Mitochondrial Genomes. Sci. Rep. 2017, 7, 10628. [Google Scholar] [CrossRef] [PubMed]
  25. Sperl, W.; Ješina, P.; Zeman, J.; Mayr, J.A.; Demeirleir, L.; VanCoster, R.; Pickova, A.; Hansikova, H.; Houšt’ková, H.; Krejčík, Z. Deficiency of Mitochondrial ATP Synthase of Nuclear Genetic Origin. Neuromuscul. Disord. 2006, 16, 821–829. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, H.; Zhang, J.; Xia, J.; Yang, J.; Guo, J.; Deng, Z.; Luo, M. Comparative Characterization of the Complete Mitochondrial Genomes of the Three Apple Snails (Gastropoda: Ampullariidae) and the Phylogenetic Analyses. Int. J. Mol. Sci. 2018, 19, 3646. [Google Scholar] [CrossRef] [PubMed]
  27. Mayrose, I.; Doron-Faigenboim, A.; Bacharach, E.; Pupko, T. Towards Realistic Codon Models: Among Site Variability and Dependency of Synonymous and Non-Synonymous Rates. Bioinformatics 2007, 23, i319–i327. [Google Scholar] [CrossRef] [PubMed]
  28. Jin, X.; Guo, X.; Chen, J.; Li, J.; Zhang, S.; Zheng, S.; Wang, Y.; Peng, Y.; Zhang, K.; Liu, Y. The Complete Mitochondrial Genome of Hemigrapsus sinensis (Brachyura, Grapsoidea, Varunidae) and Its Phylogenetic Position within Grapsoidea. Genes Genom. 2023, 45, 377–391. [Google Scholar] [CrossRef]
  29. Yang, Z.; Bielawski, J.P. Statistical Methods for Detecting Molecular Adaptation. Trends Ecol. Evol. 2000, 15, 496–503. [Google Scholar] [CrossRef]
  30. Zhang, Z.; Li, J.; Zhao, X.-Q.; Wang, J.; Wong, G.K.-S.; Yu, J. KaKs_Calculator: Calculating Ka and Ks Through Model Selection and Model Averaging. Genom. Proteom. Bioinform. 2006, 4, 259–263. [Google Scholar] [CrossRef]
  31. Saccone, C.; De Giorgi, C.; Gissi, C.; Pesole, G.; Reyes, A. Evolutionary Genomics in Metazoa: The Mitochondrial DNA as a Model System. Gene 1999, 238, 195–209. [Google Scholar] [CrossRef]
  32. Lemer, S.; Bieler, R.; Giribet, G. Resolving the Relationships of Clams and Cockles: Dense Transcriptome Sampling Drastically Improves the Bivalve Tree of Life. Proc. R. Soc. B 2019, 286, 20182684. [Google Scholar] [CrossRef] [PubMed]
  33. Audino, J.A.; Serb, J.M.; Marian, J.E.A. Phylogeny and Anatomy of Marine Mussels (Bivalvia: Mytilidae) Reveal Convergent Evolution of Siphon Traits. Zool. J. Linn. Soc. 2020, 190, 592–612. [Google Scholar] [CrossRef]
  34. Sartori, A.F.; Printrakoon, C.; Mikkelsen, P.M.; Bieler, R. Siphonal Structure in the Veneridae (Bivalvia: Heterodonta) with an Assessment of Its Phylogenetic Application and a Review of Venerids of the Gulf of Thailand. Raffles Bull. Zool. 2008, 18, 103–125. [Google Scholar]
  35. An Analytical Approach to a Classification of the Bivalvia. Philos. Trans. R. Soc. Lond. B 1978, 284, 425–436. [CrossRef]
  36. Zwarts, L.; Wanink, J. Siphon Size and Burying Depth in Deposit- and Suspension-Feeding Benthic Bivalves. Mar. Biol. 1989, 100, 227–240. [Google Scholar] [CrossRef]
  37. McLusky, D.S.; Elliott, M. The Feeding and Survival Strategies of Estuarine Molluscs. Feed. Surviv. Srategies Estuar. Org. 1981, 15, 109–121. [Google Scholar]
  38. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  39. Bernt, M.; Donath, A.; Jühling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Pütz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de Novo Metazoan Mitochondrial Genome Annotation. Mol. Phylogenetics Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef]
  40. Grant, J.R.; Enns, E.; Marinier, E.; Mandal, A.; Herman, E.K.; Chen, C.; Graham, M.; Van Domselaar, G.; Stothard, P. Proksee: In-Depth Characterization and Visualization of Bacterial Genomes. Nucleic Acids Res. 2023, 51, W484–W492. [Google Scholar] [CrossRef]
  41. Xiang, C.; Gao, F.; Jakovlić, I.; Lei, H.; Hu, Y.; Zhang, H.; Zou, H.; Wang, G.; Zhang, D. Using PhyloSuite for Molecular Phylogeny and Tree-based Analyses. iMeta 2023, 2, e87. [Google Scholar] [CrossRef]
  42. Hassanin, A.; Leger, N.; Deutsch, J. Evidence for Multiple Reversals of Asymmetric Mutational Constraints during the Evolution of the Mitochondrial Genome of Metazoa, and Consequences for Phylogenetic Inferences. Syst. Biol. 2005, 54, 277–298. [Google Scholar] [CrossRef] [PubMed]
  43. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef] [PubMed]
  44. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  45. Talavera, G.; Castresana, J. Improvement of Phylogenies after Removing Divergent and Ambiguously Aligned Blocks from Protein Sequence Alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [PubMed]
  46. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef]
  47. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; Von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef]
  48. Ronquist, F.; Teslenko, M.; Van Der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice across a Large Model Space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef]
  49. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An Online Tool for Phylogenetic Tree Display and Annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
Ijms 25 00067 g001
Figure 1. Distribution range of the genus Novaculina based on those references [3,4,5,6].
Figure 1. Distribution range of the genus Novaculina based on those references [3,4,5,6].
Ijms 25 00067 g001
Ijms 25 00067 g002
Figure 2. Mitogenome map of N. chinensis.
Figure 2. Mitogenome map of N. chinensis.
Ijms 25 00067 g002
Ijms 25 00067 g003
Figure 3. Relative synonymous codon usage (RSCU) of 12 PCGs in the mitogenome of N. chinensis.
Figure 3. Relative synonymous codon usage (RSCU) of 12 PCGs in the mitogenome of N. chinensis.
Ijms 25 00067 g003
Ijms 25 00067 g004
Figure 4. The phylogenetic tree for N. chinensis based on 12 PCGs. Below the nodes, the left is the Bayesian posterior probability (PP) value, and the right is the bootstrap proportion (BP) value. Nodes with no labels were maximally supported (PP/BP = 1/100). The red font is N. chinensis.
Figure 4. The phylogenetic tree for N. chinensis based on 12 PCGs. Below the nodes, the left is the Bayesian posterior probability (PP) value, and the right is the bootstrap proportion (BP) value. Nodes with no labels were maximally supported (PP/BP = 1/100). The red font is N. chinensis.
Ijms 25 00067 g004
Ijms 25 00067 g005
Figure 5. Three species of the genus Novaculina. (A) N. chinensis, (B) N. gangetica, and [3] (C) N. myanmarensis [3].
Figure 5. Three species of the genus Novaculina. (A) N. chinensis, (B) N. gangetica, and [3] (C) N. myanmarensis [3].
Ijms 25 00067 g005
Ijms 25 00067 g006
Figure 6. Soft body morphology of three species of the genus Novaculina. (A) N. chinensis, (B) N. gangetica, and [3] (C) N. myanmarensis [3].
Figure 6. Soft body morphology of three species of the genus Novaculina. (A) N. chinensis, (B) N. gangetica, and [3] (C) N. myanmarensis [3].
Ijms 25 00067 g006
Table 1. Organization of the mitogenome of N. chinensis.
Table 1. Organization of the mitogenome of N. chinensis.
GeneStrandPosition (Start–End)Length(bp)Intergenic_SpacerStart_CodonStop_Codon
trnK-aaa+1–6767---
cox2+113–92281045ATGTAG
trnY-tac+1086–114964163--
nad4l+1157–14442887ATGTAA
trnG-gga+1446–1511661--
trnP-cca+1518–1582656--
nad4+1696–30571362113ATTTAG
trnH-cac+3060–3123642--
trnW-tga+3126–3192672--
trnR-cga+3196–3260653--
trnE-gaa+3279–33456718--
trnS2-tca+3340–340263−6--
nad3+3403–37683660ATGTAG
trnT-aca+3768–383366−1--
trnI-atc+3837–3902663--
trnD-gac+3918–39836615--
trnQ-caa+3984–4051680--
trnC-tgc+4062–41266510--
trnA-gca+4129–4193652--
trnF-ttc+4199–4262645--
cox1+4421–61001680158ATTTAA
trnL2-tta+6128–61926527--
nad1+6193–71199270ATGTAA
trnL1-cta+7124–7191684--
trnV-gta+7192–7255640--
trnN-aac+7256–7321660--
nad5+7322–908217610ATTTAA
cob+9119–10264114636ATGTAA
nad6+10277–1080753112ATGTAG
rrnL+10808–1203912320--
atp6+12040–127386990ATGTAA
trnM-atg+12751–128166612--
rrnS+12894–1374284977--
cox3+13741–14529789−2ATGTAG
trnS1-aga+14529–1459567−1--
nad2+14596–1564810530ATGTAA
Table 2. The nucleotide composition and skewness of the mitogenome of N. chinensis.
Table 2. The nucleotide composition and skewness of the mitogenome of N. chinensis.
RegionSize (bp)A%T%G%C%AT%GC%AT-SkewGC-Skew
Mitogenome16,26228.1543.4518.869.5471.6028.40−0.2140.328
PCGs11,50225.3846.3818.819.4271.7628.24−0.2930.333
tRNAs144432.0636.2218.9112.8168.2831.72−0.0610.192
rRNAs221233.9536.5318.6310.9070.4829.52−0.040.26
Table 3. Codon number and RSCU of 12 PCGs in the mitogenome of N. chinensis. The asterisk (*) in the table indicates the stop codon.
Table 3. Codon number and RSCU of 12 PCGs in the mitogenome of N. chinensis. The asterisk (*) in the table indicates the stop codon.
CodonCountRSCUCodonCountRSCUCodonCountRSCUCodonCountRSCU
UUU (F)4011.94UCU (S)1462.89UAU (Y)1481.91UGU (C)731.78
UUC (F)120.06UCC (S)20.04UAC (Y)70.09UGC (C)90.22
UUA (L)3604.08UCA (S)460.91UAA (*)71.17UGA (W)581.2
UUG (L)850.96UCG (S)90.18UAG (*)50.83UGG (W)390.8
CUU (L)710.8CCU (P)903CAU (H)681.92CGU (R)321.94
CUC (L)10.01CCC (P)30.1CAC (H)30.08CGC (R)20.12
CUA (L)120.14CCA (P)240.8CAA (Q)331.43CGA (R)221.33
CUG (L)10.01CCG (P)30.1CAG (Q)130.57CGG (R)100.61
AUU (I)2551.98ACU (T)802.52AAU (N)961.96AGU (S)561.11
AUC (I)30.02ACC (T)20.06AAC (N)20.04AGC (S)50.1
AUA (M)1611.56ACA (T)391.23AAA (K)731.47AGA (S)1102.18
AUG (M)460.44ACG (T)60.19AAG (K)260.53AGG (S)300.59
GUU (V)2232.56GCU (A)1413.34GAU (D)911.96GGU (G)1652.06
GUC (V)20.02GCC (A)20.05GAC (D)20.04GGC (G)110.14
GUA (V)1071.23GCA (A)240.57GAA (E)571.28GGA (G)881.1
GUG (V)160.18GCG (A)20.05GAG (E)320.72GGG (G)560.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, Z.; Song, Y.; Zheng, Z.; Liu, Y.; Yao, H.; Rao, X.; Lin, G. The Complete Mitochondrial Genome and Phylogenetic Analysis of the Freshwater Shellfish Novaculina chinensis (Bivalvia: Pharidae). Int. J. Mol. Sci. 2024, 25, 67. https://doi.org/10.3390/ijms25010067

AMA Style

Zhou Z, Song Y, Zheng Z, Liu Y, Yao H, Rao X, Lin G. The Complete Mitochondrial Genome and Phylogenetic Analysis of the Freshwater Shellfish Novaculina chinensis (Bivalvia: Pharidae). International Journal of Molecular Sciences. 2024; 25(1):67. https://doi.org/10.3390/ijms25010067

Chicago/Turabian Style

Zhou, Ziquan, Yuxin Song, Zewen Zheng, Yunguang Liu, Haiyan Yao, Xiaozhen Rao, and Gang Lin. 2024. "The Complete Mitochondrial Genome and Phylogenetic Analysis of the Freshwater Shellfish Novaculina chinensis (Bivalvia: Pharidae)" International Journal of Molecular Sciences 25, no. 1: 67. https://doi.org/10.3390/ijms25010067

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop