Curious bivalves: Systematic utility and unusual properties of anomalodesmatan mitochondrial genomes

Phylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes

Despite significant advances in the phylogenomics of bivalves over the past decade, the higher-level phylogeny of Imparidentia (a superorder of Heterodonta) remains elusive. Here, a total of five new mitochondrial sequences (Chama asperella, Chama limbula, Chama dunkeri, Barnea manilensis and Ctena divergens) was added to provide resolution in nodes that required additional study. Although the monophyly of Lucinida remains less clear, the results revealed the overall backbone of the Imparidentia tree and the monophyly of Imparidentia. Likewise, most relationships among the five major Imparidentia lineages-Lucinida, Cardiida, Adapedonta, Myida and Venerida-were addressed with a well-supported topology. Basal relationships of Imparidentia recovered Lucinidae as the sister group to all remaining imparidentian taxa. Thyasiridae is a sister group to other imparidentian bivalves (except Lucinidae species) which is split into Cardiida, Adapedonta and the divergent clade of Neoheterodontei. Neoheterodontei was comprised of Venerida and Myida, the former of which now also contains Chamidae as the sister group to all the remaining venerid taxa. Moreover, molecular divergence times were inferred by calibrating nine nodes in the Imparidentia tree of life by extinct taxa. The origin of these major clades ranged from Ordovician to Permian with the diversification through the Palaeozoic to Mesozoic. Overall, the results obtained in this study demonstrate a better-resolved Imparidentia phylogeny based on mitochondrial genomes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00178-x.

The complete mitochondrial genome of Cuspidaria undata (Bivalvia, Anomalodesmata, Cuspidariidae)

The mitochondrial genome of Cuspidaria undata (Verrill, 1884) was sequenced in full using Illumina HiSeq 2500. The circular mitochondrial DNA (mtDNA) was 16,266 bp in size, encoded 37 genes, and contained 13 protein-coding genes (PCGs), 2 rRNAs and 22 tRNAs. The gene order of the 13 PCGs in this species exhibited extensive rearrangement and differences in comparison to other Cuspidariidae, indicating that gene order is not conserved within this family. Phylogenetic analysis based on 13 PCGs and 2 rRNAs recovered a monophyletic Cuspidariidae.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae

BACKGROUND

Arcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2-3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear.

RESULTS

We assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes.

CONCLUSIONS

The presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

The Phylogenetic Position of the Enigmatic, Polypodium hydriforme (Cnidaria, Polypodiozoa): Insights from Mitochondrial Genomes

Polypodium hydriforme is an enigmatic parasite that belongs to the phylum Cnidaria. Its taxonomic position has been debated: whereas it was previously suggested to be part of Medusozoa, recent phylogenomic analyses based on nuclear genes support the view that P. hydriforme and Myxozoa form a clade called Endocnidozoa. Medusozoans have linear mitochondrial (mt) chromosomes, whereas myxozoans, as most metazoan species, have circular chromosomes. In this work, we determined the structure of the mt genome of P. hydriforme, using Illumina and Oxford Nanopore Technologies reads, and showed that it is circular. This suggests that P. hydriforme is not nested within Medusozoa, as this would entail linearization followed by recirculation. Instead, our results support the view that P. hydriforme is a sister clade to Myxozoa, and mt linearization in the lineage leading to medusozoans occurred after the divergence of Myxozoa + P. hydriforme. Detailed analyses of the assembled P. hydriforme mt genome show that: (1) it is encoded on a single circular chromosome with an estimated size of ∼93,000 base pairs, making it one of the largest metazoan mt genomes; (2) around 78% of the genome encompasses a noncoding region composed of several repeat types; (3) similar to Myxozoa, no mt tRNAs were identified; (4) the codon TGA is a stop codon and does not encode for tryptophan as in other cnidarians; (5) similar to myxozoan mt genomes, it is extremely fast evolving.

Contrasting modes of mitochondrial genome evolution in sister taxa of wood-eating marine bivalves (Teredinidae and Xylophagaidae)

The bivalve families Teredinidae and Xylophagaidae include voracious consumers of wood in shallow and deep-water marine environments, respectively. The taxa are sister clades whose members consume wood as food with the aid of intracellular cellulolytic endosymbionts housed in their gills. This combination of adaptations is found in no other group of animals and was likely present in the common ancestor of both families. Despite these commonalities, the two families have followed dramatically different evolutionary paths with respect to anatomy, life history and distribution. Here we present 42 new mitochondrial genome sequences from Teredinidae and Xylophagaidae and show that distinct trajectories have also occurred in the evolution and organization of their mitochondrial genomes. Teredinidae display significantly greater rates of amino acid substitution but absolute conservation of protein-coding gene order, whereas Xylophagaidae display significantly less amino acid change but have undergone numerous and diverse changes in genome organization since their divergence from a common ancestor. As with many bivalves, these mitochondrial genomes encode two ribosomal RNAs, 12 protein coding genes, and 22 tRNAs; atp8 was not detected. We further show that their phylogeny, as inferred from amino acid sequences of 12 concatenated mitochondrial protein-coding genes, is largely congruent with those inferred from their nuclear genomes based on 18S and 28S ribosomal RNA sequences. Our results provide a robust phylogenetic framework to explore the tempo and mode of mitochondrial genome evolution and offer directions for future phylogenetic and taxonomic studies of wood-boring bivalves.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

The origins, relationships, evolution and conservation of the weirdest marine bivalves: The watering pot shells. A review

The fossil record shows that the two clavagelloid or watering pot families evolved at different times, the Clavagellidae first in the late Mesozoic (100-66mya), the Penicillidae later in the Cenozoic (33-23mya)-the former originally with, thus, a near-global Tethyan distribution, the latter restricted to the Indo-West Pacific. Representatives of the two clavagelloid families, moreover, have wholly different adventitious tube/crypt structures and, thus, methods of formation suggesting that evolutionary experiments have been undertaken to achieve such radical architectural novelties. This has resulted in one of the most surprising examples of convergent evolution in the Bivalvia. But, what were the ancestors of the Clavagelloidea? The shell and internal morphology of representatives of the three recognized genera of the Lyonsiidae, that is, Lyonsia, Entodesma and Mytilimeria, are described. Species of the latter two genera are highly specialized epibenthic, byssate, nestlers and embedded symbionts of ascidian colonies and sponges, respectively. Species of Lyonsia, however, are mostly shallow endobenthic burrowers. On the basis of these studies, it is concluded that species of Lyonsia can be regarded as representative of the ancestral watering pot (Clavagelloidea) condition. Evidence for this conclusion include the mineralogy, characteristics and ligament structure of the shell and features of the anatomy, importantly the modification of the vestigial pedal retractor muscles to form simple (Clavagellidae) and more complex (Penicillidae) proprioreceptors. Such an anatomy-based conclusion is supported to some extent by DNA analyses of representatives of the Lyonsiidae and the two constituent families of the Clavagelloidea. To some extent because all clavagelloids are exceedingly rare hindering such analyses. Such rarity, however, also argues for the strict conservation of all the species of the Clavagelloidea.

Bivalve molluscs as model systems for studying mitochondrial biology

The class Bivalvia is a highly successful and ancient taxon including ∼25,000 living species. During their long evolutionary history bivalves adapted to a wide range of physicochemical conditions, habitats, biological interactions, and feeding habits. Bivalves can have strikingly different size, and despite their apparently simple body plan, they evolved very different shell shapes, and complex anatomic structures. One of the most striking features of this class of animals is their peculiar mitochondrial biology: some bivalves have facultatively anaerobic mitochondria that allow them to survive prolonged periods of anoxia/hypoxia. Moreover, more than 100 species have now been reported showing the only known evolutionarily stable exception to the strictly maternal inheritance of mitochondria in animals, named doubly uniparental inheritance. Mitochondrial activity is fundamental to eukaryotic life, and thanks to their diversity and uncommon features, bivalves represent a great model system to expand our knowledge about mitochondrial biology, so far limited to a few species. We highlight recent works studying mitochondrial biology in bivalves at either genomic or physiological level. A link between these two approaches is still missing, and we believe that an integrated approach and collaborative relationships are the only possible ways to be successful in such endeavour.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Eight new mitogenomes clarify the phylogenetic relationships of Stromboidea within the caenogastropod phylogenetic framework

Members of the gastropod superfamily Stromboidea (Littorinimorpha) are characterised by their elaborate shell morphologies, distinctive mode of locomotion, and often large and colourful eyes. This iconic group comprises over 130 species, including many large and charismatic species. The family Strombidae is of particular interest, largely due to its commercial importance and wide distribution in tropical and subtropical waters. Although a few strombid mitochondrial genomes have been sequenced, data for the other four Recent families in Stromboidea are lacking. In this study we report seven new stromboid mitogenomes obtained from transcriptomic and genomic data, with taxonomic representation from each Recent stromboid family, including the first mitogenomes for Aporrhaidae, Rostellariidae, Seraphsidae and Struthiolariidae. We also report a new mitogenome for the family Xenophoridae. We use these data, along with published sequences, to investigate the relationships among these and other caenogastropod groups. All analyses undertaken in this study support monophyly of Stromboidea as redefined here to include Xenophoridae, a finding consistent with morphological and behavioural data. Consistent with previous morphological and molecular analyses, including those based on mitogenomes, monophyly of Hypsogastropoda is confirmed but monophyly of Littorinimorpha is again rejected.

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.

Interactions between Schistosoma haematobium group species and their Bulinus spp. intermediate hosts along the Niger River Valley

Background

Urogenital schistosomiasis, caused by infection with Schistosoma haematobium, is endemic in Niger but complicated by the presence of Schistosoma bovis, Schistosoma curassoni and S. haematobium group hybrids along with various Bulinus snail intermediate host species. Establishing the schistosomes and snails involved in transmission aids disease surveillance whilst providing insights into snail-schistosome interactions/compatibilities and biology.

Methods

Infected Bulinus spp. were collected from 16 villages north and south of the Niamey region, Niger, between 2011 and 2015. From each Bulinus spp., 20–52 cercariae shed were analysed using microsatellite markers and a subset identified using the mitochondrial (mt) cox1 and nuclear ITS1 + 2 and 18S DNA regions. Infected Bulinus spp. were identified using both morphological and molecular analysis (partial mt cox1 region).

Results

A total of 87 infected Bulinus from 24 sites were found, 29 were molecularly confirmed as B. truncatus, three as B. forskalii and four as B. globosus. The remaining samples were morphologically identified as B. truncatus (n = 49) and B. forskalii (n = 2). The microsatellite analysis of 1124 cercariae revealed 186 cercarial multilocus genotypes (MLGs). Identical cercarial genotypes were frequently (60%) identified from the same snail (clonal populations from a single miracidia); however, several (40%) of the snails had cercariae of different genotypes (2–10 MLG’s) indicating multiple miracidial infections. Fifty-seven of the B. truncatus and all of the B. forskalii and B. globosus were shedding the Bovid schistosome S. bovis. The other B. truncatus were shedding the human schistosomes, S. haematobium (n = 6) and the S. haematobium group hybrids (n = 13). Two B. truncatus had co-infections with S. haematobium and S. haematobium group hybrids whilst no co-infections with S. bovis were observed.

Conclusions

This study has advanced our understanding of human and bovid schistosomiasis transmission in the Niger River Valley region. Human Schistosoma species/forms (S. haematobium and S. haematobium hybrids) were found transmitted only in five villages whereas those causing veterinary schistosomiasis (S. bovis), were found in most villages. Bulinus truncatus was most abundant, transmitting all Schistosoma species, while the less abundant B. forskalii and B. globosus, only transmitted S. bovis. Our data suggest that species-specific biological traits may exist in relation to co-infections, snail-schistosome compatibility and intramolluscan schistosome development.

Variability of mitochondrial ORFans hints at possible differences in the system of doubly uniparental inheritance of mitochondria among families of freshwater mussels (Bivalvia: Unionida)

BACKGROUND

Supernumerary ORFan genes (i.e., open reading frames without obvious homology to other genes) are present in the mitochondrial genomes of gonochoric freshwater mussels (Bivalvia: Unionida) showing doubly uniparental inheritance (DUI) of mitochondria. DUI is a system in which distinct female-transmitted and male-transmitted mitotypes coexist in a single species. In families Unionidae and Margaritiferidae, the transition from dioecy to hermaphroditism and the loss of DUI appear to be linked, and this event seems to affect the integrity of the ORFan genes. These observations led to the hypothesis that the ORFans have a role in DUI and/or sex determination. Complete mitochondrial genome sequences are however scarce for most families of freshwater mussels, therefore hindering a clear localization of DUI in the various lineages and a comprehensive understanding of the influence of the ORFans on DUI and sexual systems. Therefore, we sequenced and characterized eleven new mitogenomes from poorly sampled freshwater mussel families to gather information on the evolution and variability of the ORFan genes and their protein products.

RESULTS

We obtained ten complete plus one almost complete mitogenome sequence from ten representative species (gonochoric and hermaphroditic) of families Margaritiferidae, Hyriidae, Mulleriidae, and Iridinidae. ORFan genes are present only in DUI species from Margaritiferidae and Hyriidae, while non-DUI species from Hyriidae, Iridinidae, and Mulleriidae lack them completely, independently of their sexual system. Comparisons among the proteins translated from the newly characterized ORFans and already known ones provide evidence of conserved structures, as well as family-specific features.

CONCLUSIONS

The ORFan proteins show a comparable organization of secondary structures among different families of freshwater mussels, which supports a conserved physiological role, but also have distinctive family-specific features. Given this latter observation and the fact that the ORFans can be either highly mutated or completely absent in species that secondarily lost DUI depending on their respective family, we hypothesize that some aspects of the connection among ORFans, sexual systems, and DUI may differ in the various lineages of unionids.

Predatory marine bivalves: A review

Most bivalves are suspension feeders. On the deep sea floor, however, some are predators, typically of meiobenthic crustaceans: copepods, cumaceans and ostracods. Propeamusiid scallops are one such group of predators. The largest numbers of predators, however, belong to the bivalve subclass Anomalodesmata and constitute, as currently recognised, some 500 species belonging principally to the Verticordioidea (120), Poromyoidea (75) and Cuspidarioidea (304) with four, two and four constituent families, respectively. A further family, the Parilimyidae, is considered to be derived from the Pholadomyoidea-the anomalodesmatan ancestor. These, generally small (<60mm shell length), nacreous and thin-shelled predators share many anatomical features that formerly allowed them to be collectively classified as the Septibranchia. Although this name is now rarely used, it refers to their possession of a ctenidially-derived septum in the mantle cavity and functioning in prey capture. Generally, there is a trend, possibly evolutionary, from a typical bivalve ctenidium (Parilimyidae and some Verticordioidea) to a complete septum (other Verticordioidea, Poromyoidea and Cuspidarioidea). In addition, the inhalant siphon, foot, labial palps, mouth and its lips play a role in prey capture, and ingestion. Similarly, the stomach is modified to digest such, typically chitinous, ingested prey. Most septibranchs are either consecutive or simultaneous hermaphrodites with self-fertilisation possibly usual and with some evidence in a few of larval brooding. Notwithstanding, the deep sea septibranch species are poorly studied with virtually nothing being known about their wider distributions, ecology, detailed reproductive strategies and life history traits.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.