Mutation and selection on the wobble nucleotide in tRNA anticodons in marine bivalve mitochondrial genomes

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

The queen conch mitogenome: intra- and interspecific mitogenomic variability in Strombidae and phylogenetic considerations within the Hypsogastropoda

Aliger gigas is an economically important and vulnerable marine species. We present a new mitogenome of A. gigas from the Mexican Caribbean and use the eight publicly available Strombidae mitogenomes to analyze intra- and interspecific variation. We present the most complete phylogenomic understanding of Hypsogastropoda to date (17 superfamilies, 39 families, 85 genera, 109 species) to revisit the phylogenetic position of the Stromboidea and evaluate divergence times throughout the phylogeny. The A. gigas mitogenome comprises 15,460 bp including 13 PCGs, 22 tRNAs, and two rRNAs. Nucleotide diversity suggested divergence between the Mexican and Colombian lineages of A. gigas. Interspecific divergence showed high differentiation among Strombidae species and demonstrated a close relationship between A. gigas and Strombus pugilis, between Lambis lambis and Harpago chiragra, and among Tridentarius dentatus/Laevistrombus canarium/Ministrombus variabilis. At the intraspecific level, the gene showing the highest differentiation is ATP8 and the lowest is NAD4L, whereas at the interspecific level the NAD genes show the highest variation and the COX genes the lowest. Phylogenomic analyses confirm that Stromboidea belongs in the non-Latrogastropoda clade and includes Xenophoridea. The phylogenomic position of other superfamilies, including those of previously uncertain affiliation, is also discussed. Finally, our data indicated that Stromboidea diverged into two principal clades in the early Cretaceous while Strombidae diversified in the Paleocene, and lineage diversification within A. gigas took place in the Pleistocene.

Mitogenomic architecture of the multivalent endemic black clam (Villorita cyprinoides) and its phylogenetic implications

The Indian black clam Villorita cyprinoides (Family: Cyrenidae), an extractive commercially exploited species with aquaculture importance contributing more than 70% of clam fishery in India, is endemic to the Indian peninsula. Currently, there is very sparse information, especially on the molecular data of Villorita. The present study aims to provide a comprehensive knowledge of mitogenome architecture and assess the phylogenetic status of Cyrenidae. This has resulted in reporting the first complete mitogenome of V. cyprinoides using next-generation sequencing technology. The A+T circular mitogenome was 15,880 bp long, exhibiting 13 protein-coding genes (PCGs) including ATP8 (absent in several bivalves), 22 transfer RNA, and two ribosomal RNA genes residing in the heavy strand in a clockwise orientation and a gene order akin to Corbicula fluminea. The molecular phylogeny inferred from a concatenated multi-gene sequence [14 mitochondrial (12 PCGs, rrnS and rrnL) and two nuclear genes (Histone H3, 18S rRNA)] from 47 representative species of superorder Imparidentia, clustered V. cyprinoides and Cyrenid clams to a single clade supporting the monophyly of Cyrenidae. The subsequent mitochondrial gene order analysis substantiates the close relationship of V. cyprinoides and C. fluminea, analogous to phylogenetic output. The multilocus tree topology calibrated with verified fossil data deciphered the origin and diversification of Cyrenid clams during late Triassic-early Jurassic. The data derived from this study shall contribute remarkably for further insights on cryptic species identification, molecular characterization of bivalve mitogenomes and mitochondrial evolutionary history of genus Villorita. Moreover, complete mitogenome can aid in potential marker development for assessing the genetic health of black clam populations.

New insights into the factors affecting synonymous codon usage in human infecting Plasmodium species

Codon usage bias is due to the non-random usage of synonymous codons for coding amino acids. The synonymous sites are under weak selection, and codon usage bias is maintained by the equilibrium in mutational bias, genetic drift and selection pressure. The differential codon usage choices are also relevant to human infecting Plasmodium species. Recently, P. knowlesi switches its natural host, long-tailed macaques, and starts infecting humans. This review focuses on the comparative analysis of codon usage choices among human infecting P. falciparum and P. vivax along with P. knowlesi species taking their coding sequence data. The variation in GC content, amino acid frequencies, effective number of codons and other factors plays a crucial role in determining synonymous codon choices. Within species codon choices are more similar for P. vivax and P. knowlesi in comparison with P. falciparum species. This study suggests that synonymous codon choice modulates the gene expression level, mRNA stability, ribosome speed, protein folding, translation efficiency and its accuracy in Plasmodium species, and provides a valuable information regarding the codon usage pattern to facilitate gene cloning as well as expression and transfection studies for malaria causing species.

The complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus (Bivalvia, Veneridae)

Despite the large number of animal complete mitochondrial genomes currently available in public databases, knowledge about mitochondrial genomics in invertebrates is uneven. This paper reports, for the first time, the complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus, also known as the European clam. Ruditapes decussatus is morphologically and ecologically similar to the Manila clam Ruditapes philippinarum, which has been recently introduced for aquaculture in the very same habitats of Ruditapes decussatus, and that is replacing the native species. Currently the production of the European clam is almost insignificant, nonetheless it is considered a high value product, and therefore it is an economically important species, especially in Portugal, Spain and Italy. In this work we: (i) assembled Ruditapes decussatus mitochondrial genome from RNA-Seq data, and validated it by Sanger sequencing; (ii) analyzed and characterized the Ruditapes decussatus mitochondrial genome, comparing its features with those of other venerid bivalves; (iii) assessed mitochondrial sequence polymorphism (SP) and copy number variation (CNV) of tandem repeats across 26 samples. Despite using high-throughput approaches we did not find evidence for the presence of two sex-linked mitochondrial genomes, typical of the doubly uniparental inheritance of mitochondria, a phenomenon known in ∼100 bivalve species. According to our analyses, Ruditapes decussatus is more genetically similar to species of the Genus Paphia than to the congeneric Ruditapes philippinarum, a finding that bolsters the already-proposed need of a taxonomic revision. We also found a quite low genetic variability across the examined samples, with few SPs and little variability of the sequences flanking the control region (Largest Unassigned Regions (LURs). Strikingly, although we found low nucleotide variability along the entire mitochondrial genome, we observed high levels of length polymorphism in the LUR due to CNV of tandem repeats, and even a LUR length heteroplasmy in two samples. It is not clear if the lack of genetic variability in the mitochondrial genome of Ruditapes decussatus is a cause or an effect of the ongoing replacement of Ruditapes decussatus with the invasive Ruditapes philippinarum, and more analyses, especially on nuclear sequences, are required to assess this point.

The complete mitochondrial genome of Crassostrea gasar (Bivalvia: Ostreidae)

The complete mitochondrial genome of Crassostrea gasar was sequenced using the Ion Proton technology in combination with 454 Roche GS-FLX plataform data. We assembled a 17,686 bp complete circular mitochondrial genome, containing 13 protein-coding genes, a major non-coding region (MNR), two ribosomal RNA genes and 24 transfer RNA genes. Phylogenetic analysis of concatenated amino acid sequences from mitochondria showed monophyletic clades formed with high bootstrap values. This is the first complete mitochondrial sequence of an oyster from South America. Mitogenome sequence was deposited in GenBank under the accession number KR856227.

Intraspecific Variation in Mitogenomes of Five Crassostrea Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China's coast, which are shaped by unknown mechanisms in a north-south divide.

Comparative mitogenomic analysis reveals cryptic species: A case study in Mactridae (Mollusca: Bivalvia)

The Chinese surf clam Mactra chinensis Philippi, 1846 is a commercially important marine bivalve belonging to the family Mactridae (Mollusca: Bivalvia). In this study, the M. chinensis mitochondrial genomic features are analyzed. The genome has 34 genes on the same strand, lacking atp8 and both trnS (trnS1 and trnS2) as compared with the typical gene content of metazoan mitochondrial genomes. The A+T content of M. chinensis mitochondrial genome is 63.72%, which is slightly lower than that of M. veneriformis (67.59%) and Coelomactra antiquata (64.33% and 64.14% for the samples from Ri Zhao, Shandong Province, and Zhang Zhou, Fujian Province, China, respectively) in the same family. There are 22 NCRs in the M. chinensis mitochondrial genome, accounting for 12.91% of the genome length. The longest NCR (1,075bp in length) is located between trnT and trnQ. A TRS (127bp×8.15) accounts for 96.3% (1,035/1,075) of this NCR. The occurrence of TRS in NCR is shared by the two Mactra mitochondrial genomes, but is not found in the two Coelomactra mitochondrial genomes. A phylogenetic tree constructed based on 12 PCGs of 25 bivalve mitochondrial genomes shows that all seven genera (Mactra, Coelomactra, Paphia, Meretrix, Solen, Mytilus, and Crassostrea) constitute monophyletic groups with very high support values. Pairwise genetic distance analyses indicate that the genetic distance of C. antiquata from the two localities is 0.084, which is greater than values between congeneric species, such as those in Mactra, Mytilus, Meretrix, and Crassostrea. The results show that the C. antiquata from the two localities represent cryptic species.

Mutation and selection cause codon usage and bias in mitochondrial genomes of ribbon worms (Nemertea)

The phenomenon of codon usage bias is known to exist in many genomes and it is mainly determined by mutation and selection. To understand the patterns of codon usage in nemertean mitochondrial genomes, we use bioinformatic approaches to analyze the protein-coding sequences of eight nemertean species. Neutrality analysis did not find a significant correlation between GC12 and GC3. ENc-plot showed a few genes on or close to the expected curve, but the majority of points with low-ENc values are below it. ENc-plot suggested that mutational bias plays a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that GC and AT were not used proportionally and we propose that codons containing A or U at third position are used preferentially in nemertean species, regardless of whether corresponding tRNAs are encoded in the mitochondrial DNA. Context-dependent analysis indicated that the nucleotide at the second codon position slightly affects synonymous codon choices. These results suggested that mutational and selection forces are probably acting to codon usage bias in nemertean mitochondrial genomes.

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.

Base-pairing versatility determines wobble sites in tRNA anticodons of vertebrate mitogenomes

BACKGROUND

Vertebrate mitochondrial genomes typically have one transfer RNA (tRNA) for each synonymous codon family. This limited anticodon repertoire implies that each tRNA anticodon needs to wobble (establish a non-Watson-Crick base pairing between two nucleotides in RNA molecules) to recognize one or more synonymous codons. Different hypotheses have been proposed to explain the factors that determine the nucleotide composition of wobble sites in vertebrate mitochondrial tRNA anticodons. Until now, the two major postulates--the "codon-anticodon adaptation hypothesis" and the "wobble versatility hypothesis"--have not been formally tested in vertebrate mitochondria because both make the same predictions regarding the composition of anticodon wobble sites. The same is true for the more recent "wobble cost hypothesis".

PRINCIPAL FINDINGS

In this study we have analyzed the occurrence of synonymous codons and tRNA anticodon wobble sites in 1553 complete vertebrate mitochondrial genomes, focusing on three fish species with mtDNA codon usage bias reversal (L-strand is GT-rich). These mitogenomes constitute an excellent opportunity to study the evolution of the wobble nucleotide composition of tRNA anticodons because due to the reversal the predictions for the anticodon wobble sites differ between the existing hypotheses. We observed that none of the wobble sites of tRNA anticodons in these unusual mitochondrial genomes coevolved to match the new overall codon usage bias, suggesting that nucleotides at the wobble sites of tRNA anticodons in vertebrate mitochondrial genomes are determined by wobble versatility.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that, at wobble sites of tRNA anticodons in vertebrate mitogenomes, selection favors the most versatile nucleotide in terms of wobble base-pairing stability and that wobble site composition is not influenced by codon usage. These results are in agreement with the "wobble versatility hypothesis".

Complete mitochondrial DNA sequence of the European flat oyster Ostrea edulis confirms Ostreidae classification

BACKGROUND

Because of its typical architecture, inheritance and small size, mitochondrial (mt) DNA is widely used for phylogenetic studies. Gene order is generally conserved in most taxa although some groups show considerable variation. This is particularly true in the phylum Mollusca, especially in the Bivalvia. During the last few years, there have been significant increases in the number of complete mitochondrial sequences available. For bivalves, 35 complete mitochondrial genomes are now available in GenBank, a number that has more than doubled in the last three years, representing 6 families and 23 genera. In the current study, we determined the complete mtDNA sequence of O. edulis, the European flat oyster. We present an analysis of features of its gene content and genome organization in comparison with other Ostrea, Saccostrea and Crassostrea species.

RESULTS

The Ostrea edulis mt genome is 16 320 bp in length and codes for 37 genes (12 protein-coding genes, 2 rRNAs and 23 tRNAs) on the same strand. As in other Ostreidae, O. edulis mt genome contains a split of the rrnL gene and a duplication of trnM. The tRNA gene set of O. edulis, Ostrea denselamellosa and Crassostrea virginica are identical in having 23 tRNA genes, in contrast to Asian oysters, which have 25 tRNA genes (except for C. ariakensis with 24). O. edulis and O. denselamellosa share the same gene order, but differ from other Ostreidae and are closer to Crassostrea than to Saccostrea. Phylogenetic analyses reinforce the taxonomic classification of the 3 families Ostreidae, Mytilidae and Pectinidae. Within the Ostreidae family the results also reveal a closer relationship between Ostrea and Saccostrea than between Ostrea and Crassostrea.

CONCLUSIONS

Ostrea edulis mitogenomic analyses show a high level of conservation within the genus Ostrea, whereas they show a high level of variation within the Ostreidae family. These features provide useful information for further evolutionary analysis of oyster mitogenomes.

Identifying the true oysters (Bivalvia: Ostreidae) with mitochondrial phylogeny and distance-based DNA barcoding

Oysters (family Ostreidae), with high levels of phenotypic plasticity and wide geographic distribution, are a challenging group for taxonomists and phylogenetics. As a useful tool for molecular species identification, DNA barcoding might offer significant potential for oyster identification and taxonomy. This study used two mitochondrial fragments, cytochrome c oxidase I (COI) and the large ribosomal subunit (16S rDNA), to assess whether oyster species could be identified by phylogeny and distance-based DNA barcoding techniques. Relationships among species were estimated by the phylogenetic analyses of both genes, and then pairwise inter- and intraspecific genetic divergences were assessed. Species forming well-differentiated clades in the molecular phylogenies were identical for both genes even when the closely related species were included. Intraspecific variability of 16S rDNA overlapped with interspecific divergence. However, average intra- and interspecific genetic divergences for COI were 0-1.4% (maximum 2.2%) and 2.6-32.2% (minimum 2.2%), respectively, indicating the existence of a barcoding gap. These results confirm the efficacy of species identification in oysters via DNA barcodes and phylogenetic analysis.