Complete mitochondrial genome of Coelomactra antiquata (Mollusca: Bivalvia): The first representative from the family Mactridae with novel gene order and unusual tandem repeats

SMRT Sequencing of the Full-Length Transcriptome of the Coelomactra antiquata

Coelomactra antiquata is an important aquatic economic shellfish with high medicinal value. However, because C. antiquata has no reference genome, a lot of molecular biology research cannot be carried out, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. antiquata. In the present study, we conducted the first full-length transcriptome analysis of C. antiquata by using PacBio single-molecule real-time (SMRT) sequencing technology. The results identified a total of 39,209 unigenes with an average length of 2,732 bp, 23,338 CDSs, 251 AS events, 9,881 lncRNAs, 20,106 SSRs, and 2,316 TFs. Subsequently, 59.22% (23,220) of the unigenes were successfully annotated, of which 23,164, 18,711, 15,840, 13,534, and 13,474 unigenes could be annotated using NR, Swiss-prot, KOG, GO, and KEGG databases, respectively. This study lays the foundation for the follow-up research of molecular biology and provides a reference for studying the more medicinal value of C. antiquata.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

The first representative of Coelomactra antiquata mitochondrial genome from Liaoning (China) and phylogenetic consideration

Coelomactra antiquata is a famous delicacy and a promising new candidate for aquaculture, which belongs to the family Mactridae (Mollusca: Veneroida). The complete mitochondrial genome of C. antiquata (Liao Ning province, in China, LN) was finished, which is the first representative from this province. The results showed that the total length of LN-mtDNA sequence is 16,797 bp and the content of A + T is 65.01%. It encodes 35 genes, including 12 protein-coding genes, 21 transfer RNA genes and two ribosomal RNA genes. All coding genes are encoded on the heavy strand. Compared with the typical gene content of animal mitochondrial genomes, atp8 and trnSer(UCN) genes are missing in the mitochondrial genome. The complete mitochondrial genome contains 26 non-coding regions (1598 bp), one major non-coding region consists of 1046 bp in which 4.9 tandem repeat sequences (99bp per sequence) was observed. The phylogenetic tree showed that Liaoning population was clustered into one clade with Shandong (Rizhao, Jiaonan and Jimo) and Guangxi (Beihai) populations, meanwhile all of them are far from the Fujian populations (Pingtan, Zhangzhou and Changle). So, Liaoning, Shandong and Guangxi populations have a close relationship. Actually, Fujian is located between Liaoning, Shandong and Guangxi. So, the result challenges the previously assumed relevance between geographic distance and genetic distance. The genetic distance of Liaoning C. antiquata and Fujian (Changle, Zhangzhou and Pingtan) C. antiquata (0.176–0.177) is greater than the genetic distance between Mytilus galloprovincialis and Mytilus trossulus (0.160). The genetic difference of Liaoning population and Fujian populations has reached species level.

Comparative genomic analysis of mitochondrial protein-coding genes in Veneroida clams: Analysis of superfamily-specific genomic and evolutionary features

Veneroida is the largest order of bivalves, and these clams are commercially important in Asian countries. Although numerous studies have focused on the genomic characters of individual species or genera in Veneroida, superfamily-specific genomic characters have not been determined. In this study, we performed a comparative genomic analysis of 12 mitochondrial protein coding genes (PCGs) from 25 clams in six Veneroida superfamilies to determine genomic and evolutionary features of each superfamily. Length and distribution of nucleotides encoding the PCGs were too variable to define superfamily-specific genomic characters. Phylogenetic analysis revealed that PCGs are suitable for classification of species in three superfamilies: Cardioidea, Mactroidea, and Veneroidea. However, one species classified in Tellinoidea, Sinonovacula constricta, was evolutionarily closer to Solenoidea clams than Tellinoidea clams. dN/dS analysis showed that positively selected sites in NADH dehydrogenase subunit, nd4 and subunit of ATP synthase, atp6 were present in Mactroidea. Differences in selected sites in the nd4 and atp6 could be caused by superfamily-level differences in sodium transport or ATP synthesis functions, respectively. These differences in selected sites in NADH may have conferred these animals, which have low motility and do not generally move, with increased flexibility to maintain homeostasis in the face of osmotic pressure. Our study provides insight into evolutionary traits as well as facilitates identification of veneroids.

The first mitochondrial genome of Coelomactra antiquata (Mollusca: Veneroida: Mactridae) from Guangxi (China) and potential molecular markers

The complete mitochondrial genome of Coelomactra antiquata (Guangxi, in China, GX) was determined. It is 16 801 bp in length and is the first representative from this province. The mitochondrial genome encodes 35 genes, including 12 PCGs, two ribosomal RNA, and 21 transfer RNA genes. Atp8 and trnSer(UCN) genes are missing, compared with the typical gene content of animal mitochondrial genomes. Three (cob, nad1, nad4, and nad6) of the 12 PCGs in the mitochondrial genome initiate with the ATA, while other PCGs start with ATG. Two PCGs (atp6 and nad4L) end with incomplete stop codons (T-), and the remaining ones have complete stop codons (TAA or TAG). The largest non-coding region of the C. antiquata (GX) contains one section of tandem repeats (5 × 99 bp). Among all PCGs and rRNAs, the nad5 gene contains the maximum polymorphic sites (430), followed by nad4 (261) and cox2 (240). Two ribosomal RNA genes (srRNA and lrRNA) and cox1 are most conservative. The proportions of polymorphic sites in six genes (nad4, nad2, nad6, nad5, cox2, and nad3) are more than 20% (ranging from 20.25% to 25.21%). These high variable genes can be used as molecular markers in the population genetic analysis of the species.

The complete mitogenome of the giant clam Tridacna squamosa (Heterodonta: Bivalvia: Tridacnidae)

The complete mitochondrial genome of the commercially and ecologically important and internationally vulnerable giant clam Tridacna squamosa was recovered by genome skimming using the MiSeq platform. The T. squamosa mitogenome has 20,930 base pairs (62.35% A+T content) and is made up of 12 protein-coding genes, 2 ribosomal subunit genes, 24 transfer RNAs, and a 2594 bp non-coding AT-rich region. The mitogenome has a relatively large insertion in the atp6 gene. This is the first mitogenome to be sequenced from the genus Tridacna, and the family Tridacnidae and represents a new gene order.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities

BACKGROUND

Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform.

RESULTS

The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations.

CONCLUSIONS

Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.

The complete mitochondrial DNA of Tegillarca granosa and comparative mitogenomic analyses of three Arcidae species

To better understand the characteristics and the evolutionary dynamics of mt genomes in Arcidae, the complete mitochondrial genome of Tegillarca granosa was firstly determined and compared with other two Arcidae species (Scapharca broughtonii and Scapharca kagoshimensis). The complete mitochondrial genome of T. granosa was 31,589 bp in length, including 12 protein-coding genes, 2 rRNA genes and 23 tRNA genes, and a major non-coding region. Three tandem repeat fragments were identified in the major non-coding region and the tandem repeat motifs of these fragments can be folded into stem-loop structures. The mitochondrial genome of the three species has several common features such as the AT content, the arrangement of the protein-coding genes, the codon usage of the protein-coding genes and AT/GC skew. However, a high level of variability is presented in the size of the genome, the number of tRNA genes and the length of non-coding sequences in the three mitogenomes. According to the phylogenetic analyses, these mitogenome-level characters are correlated with their phylogenetic relationships. It is the absence of the duplicated tRNAs and large non-coding sequences that are responsible for the length divergence of mitogenomes between T. granosa and other two Arcidae species. The phylogenetic analyses were conducted based on 12 partitioned protein genes, which support the relationship at the family level: (((Pectinidae+Ostreidae)+Mytilidae)+Arcidae).

The complete mitogenome of the marine bivalve Lutraria rhynchaena Jonas 1844 (Heterodonta: Bivalvia: Mactridae)

The complete mitochondrial genome of the commercially important snout otter clam Lutraria rhynchaena was obtained from low-coverage shotgun sequencing data on the MiSeq platform. The L. rhynchaena mitogenome has 16,927 base pairs (69% A + T content) and made up of 12 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a 953 bp non-coding AT-rich region. This is the first mitogenome to be sequenced from the genus Lutraria, and the seventh to be reported for the family Mactridae.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with conchifera

Background

Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet.

Results

In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene.

Conclusions

We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.