The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

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.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Phylogeny and evolutionary radiation of the marine mussels (Bivalvia: Mytilidae) based on mitochondrial and nuclear genes

The marine mussels (Mytilidae) are distributed in the oceans worldwide and occupy various habitats with diverse life styles. However, their taxonomy and phylogeny remain unclear from genus to family level due to equivocal morphological and anatomical characters among some taxa. In this study, we inferred the deep phylogenetic relationships among 42 mytiloid species, 19 genera, and five subfamilies of the extant marine mussels by using two mitochondrial (COI and 16S rRNA) and three nuclear (18S and 28S rRNA, and histone H3) genes. Phylogeny was reconstructed with a combination of five genes using Bayesian inference and maximum likelihood method, and divergence time was estimated for the major nodes using a relaxed clock model with three fossil calibrations. Phylogenetic trees revealed two major clades (Clades 1 and 2). In Clade 1, the deep-sea mussels (subfamily Bathymodiolinae) were sister to subfamily Modiolinae (represented by Modiolus), and then was clustered with Leiosolenus (subfamily Lithophaginae). Clade 2 comprised Lithophaga (Lithophaginae) and subfamily Mytilinae. Additionally, a Modiolus species and Musculus senhousia (subfamily Crenellinae) were positioned within the subfamily Mytilinae. The phylogenetic results strongly indicated monophyly of Mytilidae and Bathymodiolinae, polyphyly of Modiolinae and Lithophaginae, and paraphyly of Mytilinae. Divergence time estimation showed an ancient and gradual divergence in most mussel groups, whereas the deep-sea mussels originated recently and diverged rapidly during the Paleogene. The present study provides new insight into the evolutionary history of the marine mussels, and supports taxonomic revision for this important bivalve group.

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

Divergent chemosymbiosis-related characters in Thyasira cf. gouldi (Bivalvia: Thyasiridae)

Within the marine bivalve family Thyasiridae, some species have bacterial chemosymbionts associated with gill epithelial cells while other species are asymbiotic. Although the abundance of symbionts in a particular thyasirid species may vary, the structure of their gills (i.e., their frontal-abfrontal thickening) does not. We examined gill structure in a species tentatively identified as Thyasira gouldi from a Northwest Atlantic fjord (Bonne Bay, Newfoundland) and found remarkable differences among specimens. Some individuals had thickened gill filaments with abundant symbionts, while others had thin filaments and lacked symbionts. We could differentiate symbiotic and asymbiotic specimens based on the size and outline of their shell as well as 18S rRNA, 28S rRNA and CO1 sequences. The wide morphological, genetic and symbiosis-related disparity described herein suggests that chemosymbiosis may influence host divergence, and that Thyasira gouldi forms a cryptic species complex.

First molecular phylogeny of the circumtropical bivalve family Pinnidae (Mollusca, Bivalvia): evidence for high levels of cryptic species diversity

The family Pinnidae Leach, 1819, includes approximately 50 species of large subtidal and coastal marine bivalves. These commercially important species occur in tropical and temperate waters around the world and are most frequently found in seagrass meadows. The taxonomy of the family has been revised a number of times since the early 20th Century, the most recent revision recognizing 55 species distributed in three genera: Pinna, Atrina and Streptopinna, the latter being monotypic. However, to date no phylogenetic analysis of the family has been conducted using morphological or molecular data. The present study analyzed 306 pinnid specimens from around the world, comprising the three described genera and ca. 25 morphospecies. We sequenced the mitochondrial genes 16S rRNA and cytochrome c oxidase subunit I, and the nuclear ribosomal genes 18S rRNA and 28S rRNA. Phylogenetic analysis of the data revealed monophyly of the genus Atrina but also that the genus Streptopinna is nested within Pinna. Based on the strong support for this relationship we propose a new status for Streptopinna Martens, 1880 and treat it as a subgenus (status nov.) of Pinna Linnaeus, 1758. The phylogeny and the species delimitation analyses suggest the presence of cryptic species in many morphospecies displaying a wide Indo-Pacific distribution, including Pinna muricata, Atrina assimilis, A. exusta and P. (Streptopinna) saccata but also in the Atlantic species A. rigida. Altogether our results highlight the challenges associated with morphological identifications in Pinnidae due to the presence of both phenotypic plasticity and morphological stasis and reveal that many pinnid species are not as widely distributed as previously thought.

Comparative morphology among representatives of main taxa of Scaphopoda and basal protobranch Bivalvia (Mollusca)

This study deals with detailed morphology and anatomy of 4 species of Scaphopoda and 5 species of protobranch Bivalvia. Both classes are traditionally grouped in the taxon Diasoma, which has been questioned by different methodologies, such as molecular and developmental. This study is developed under a phylogenetic methodology with the main concern in performing it in an intelligible and testable methodology. The analyzed Scaphopoda species came from the Brazilian coast and belong to the family Dentaliidae [(1) Coccodentalium carduus; (2) Paradentalium disparile] and Gadiliidae; [(3) Polyschides noronhensis, n. sp. from Fernando de Noronha Archipelago; (4) Gadila braziliensis]. These species represent the main branches of the class Scaphopoda. From protobranch bivalves, representatives of the families Solemyidae [(5) Solemya occidentalis, from Florida; S. notialis, n. sp. from S.E. Brazil], Nuculanidae [(6) Propeleda carpentieri from Florida], and Nuculidae [(7) Ennucula puelcha, from south Brazil] are included. These species represent the main branches of the basal Bivalvia. The descriptions on the anatomy of S. occidentalis and of P. carpentieri are published elsewhere. The remaining are included here, for which a complete taxonomical treatment is performed. Beyond these species, representatives of other taxa are operationally included as part of the ingroup (indices are then shared with them), as a procedure to test the morphological monophyly of Diasoma. These taxa are: two lamellibranch bivalves [(8) Barbatia - Arcidae; (9) Serratina - Tellinidae; both published elsewhere;, and Propilidium (10) Patellogastropoda, and (11) Nautilus, basal Cephalopoda, based on basal taxa. The effective outgroups are (12) Neopilina (Monoplacophora) and (13) Hanleya (Polyplacophora). The phylogenetic analysis based on morphology revealed that the taxon Diasoma is supported by 14 synapomorphies, and is separated from Cyrtosoma (Gastropoda + Cephalopoda). Although they are not the main goal of this paper, the taxa Scaphopoda and Bivalvia are supported by 8 and by 7 synapomorphies respectively. The taxon Protobranchia resulted paraphyletic. Both scaphopod orders resulted monophyletic. The obtained cladogram is: ((((Coccodentalium carduus - Paradentalium disparile) (Polyschides noronhensis - Gadila brasiliensis)) ((Solemya occidentalis - S. notialis) (Propeleda carpenteri (Ennucula puelcha (Barbatia cancellaria - Serratina capsoides))))) (Propilidium curumim - Nautilus pompilius - Lolliguncula brevis)).

Comparative morphology among representatives of main taxa of Scaphopoda and basal protobranch Bivalvia (Mollusca)

This study deals with detailed morphology and anatomy of 4 species of Scaphopoda and 5 species of protobranch Bivalvia. Both classes are traditionally grouped in the taxon Diasoma, which has been questioned by different methodologies, such as molecular and developmental. This study is developed under a phylogenetic methodology with the main concern in performing it in an intelligible and testable methodology. The analyzed Scaphopoda species came from the Brazilian coast and belong to the family Dentaliidae [(1) Coccodentalium carduus; (2) Paradentalium disparile] and Gadiliidae; [(3) Polyschides noronhensis, n. sp. from Fernando de Noronha Archipelago; (4) Gadila braziliensis]. These species represent the main branches of the class Scaphopoda. From protobranch bivalves, representatives of the families Solemyidae [(5) Solemya occidentalis, from Florida; S. notialis, n. sp. from S.E. Brazil], Nuculanidae [(6) Propeleda carpentieri from Florida], and Nuculidae [(7) Ennucula puelcha, from south Brazil] are included. These species represent the main branches of the basal Bivalvia. The descriptions on the anatomy of S. occidentalis and of P. carpentieri are published elsewhere. The remaining are included here, for which a complete taxonomical treatment is performed. Beyond these species, representatives of other taxa are operationally included as part of the ingroup (indices are then shared with them), as a procedure to test the morphological monophyly of Diasoma. These taxa are: two lamellibranch bivalves [(8) Barbatia - Arcidae; (9) Serratina - Tellinidae; both published elsewhere;, and Propilidium (10) Patellogastropoda, and (11) Nautilus, basal Cephalopoda, based on basal taxa. The effective outgroups are (12) Neopilina (Monoplacophora) and (13) Hanleya (Polyplacophora). The phylogenetic analysis based on morphology revealed that the taxon Diasoma is supported by 14 synapomorphies, and is separated from Cyrtosoma (Gastropoda + Cephalopoda). Although they are not the main goal of this paper, the taxa Scaphopoda and Bivalvia are supported by 8 and by 7 synapomorphies respectively. The taxon Protobranchia resulted paraphyletic. Both scaphopod orders resulted monophyletic. The obtained cladogram is: ((((Coccodentalium carduus - Paradentalium disparile) (Polyschides noronhensis - Gadila brasiliensis)) ((Solemya occidentalis - S. notialis) (Propeleda carpenteri (Ennucula puelcha (Barbatia cancellaria - Serratina capsoides))))) (Propilidium curumim - Nautilus pompilius - Lolliguncula brevis)).

Complete sequences of the highly rearranged molluscan mitochondrial genomes of the Scaphopod Graptacme eborea and the bivalve Mytilus edulis

We have determined the complete sequence of the mitochondrial genome of the scaphopod mollusk Graptacme eborea (14,492 nts) and completed the sequence of the mitochondrial genome of the bivalve mollusk Mytilus edulis (16,740 nts). (The name Graptacme eborea is a revision of the species formerly known as Dentalium eboreum.) G. eborea mtDNA contains the 37 genes that are typically found and has the genes divided about evenly between the two strands, but M. edulis contains an extra trnM and is missing atp8, and it has all genes on the same strand. Each has a highly rearranged gene order relative to each other and to all other studied mtDNAs. G. eborea mtDNA has almost no strand skew, but the coding strand of M. edulis mtDNA is very rich in G and T. This is reflected in differential codon usage patterns and even in amino acid compositions. G. eborea mtDNA has fewer noncoding nucleotides than any other mtDNA studied to date, with the largest noncoding region only 24 nt long. Phylogenetic analysis using 2,420 aligned amino acid positions of concatenated proteins weakly supports an association of the scaphopod with gastropods to the exclusion of Bivalvia, Cephalopoda, and Polyplacophora, but it is generally unable to convincingly resolve the relationships among major groups of the Lophotrochozoa, in contrast to the good resolution seen for several other major metazoan groups.