Mitochondrial genomic analyses provide new insights into the "missing" atp8 and adaptive evolution of Mytilidae

Comparative mitogenomic analyses of the infraclass Pteriomorphia (Mollusca: Bivalvia) provides novel insights into gene rearrangement and phylogeny

The infraclass Pteriomorphia within the phylum Mollusca, class Bivalvia, encompasses a diverse group of bivalve mollusks characterized by their unique morphological and ecological traits. However, the molecular data of Pteriomorphia remains limited, constraining classification and phylogenetic analysis within the infraclass Pteriomorphia. In this study, we sequenced the mitogenomes of 15 pteriomorphians using high-throughput sequencing technology and conducted comparative genomic analysis with 58 sequences available on GenBank. The newly sequenced mitogenomes all exhibited a circular double-stranded structure, ranging in size from 12,844 to 28,105 bp. Mitogenome features, including nucleotide composition, codon usage, and amino acid content, exhibited significant deviations. The mitochondrial gene orders (MGOs) in pteriomorphians also exhibited significant diversity, with 44 distinct patterns identified (excluding the tRNA gene). The gene rearrangement analysis inferred that the ancestral MGO evolution pathway supported the hypothesis that the common ancestor MGO of Mytilida is the closest to that of the entire infraclass Pteriomorphia. The reconstructed phylogeny also indicated the order Mytilida as the basal group of Pteriomorphia. Furthermore, the analysis revealed that Ostreida and Pectinida were sister taxa, while Ostreoidea and the Pinnoidea were sister taxa within Ostreida. Gene rearrangement analysis further suggested that the MGOs of Arcida and Pectinida likely evolved from the ancestral MGO of Ostreida. Our study provided new insights into the phylogenetic relationships and evolution of gene rearrangements within the subclass Pteriomorphia.

Mitonuclear compatibility is maintained despite relaxed selection on male mitochondrial DNA in bivalves with doubly uniparental inheritance

Mitonuclear coevolution is common in eukaryotes, but bivalve lineages that have doubly uniparental inheritance (DUI) of mitochondria may be an interesting example. In this system, females transmit mtDNA (F mtDNA) to all offspring, while males transmit a different mtDNA (M mtDNA) solely to their sons. Molecular evolution and functional data suggest oxidative phosphorylation (OXPHOS) genes encoded in M mtDNA evolve under relaxed selection due to their function being limited to sperm only (vs. all other tissues for F mtDNA). This has led to the hypothesis that mitonuclear coevolution is less important for M mtDNA. Here, we use comparative phylogenetics, transcriptomics, and proteomics to understand mitonuclear interactions in DUI bivalves. We found nuclear OXPHOS proteins coevolve and maintain compatibility similarly with both F and M mtDNA OXPHOS proteins. Mitochondrial recombination did not influence mitonuclear compatibility and nuclear-encoded OXPHOS genes were not upregulated in tissues with M mtDNA to offset dysfunction. Our results support that selection maintains mitonuclear compatibility with F and M mtDNA despite relaxed selection on M mtDNA. Strict sperm transmission, lower effective population size, and higher mutation rates may explain the evolution of M mtDNA. Our study highlights that mitonuclear coevolution and compatibility may be broad features of eukaryotes.

A chromosome-level genome assembly of the disco clam, Ctenoides ales

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales (C. ales), which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and reported a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein-coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, is a valuable resource for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

Structure, Evolution, and Mitochondrial Genome Analysis of Mussel Species (Bivalvia, Mytilidae)

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

Different Evolutionary Trends of Galloanseres: Mitogenomics Analysis

The two existing clades of Galloanseres, orders Galliformes (landfowl) and Anseriformes (waterfowl), exhibit dramatically different evolutionary trends. Mitochondria serve as primary sites for energy production in organisms, and numerous studies have revealed their role in biological evolution and ecological adaptation. We assembled the complete mitogenome sequences of two species of the genus Aythya within Anseriformes: Aythya baeri and Aythya marila. A phylogenetic tree was constructed for 142 species within Galloanseres, and their divergence times were inferred. The divergence between Galliformes and Anseriformes occurred ~79.62 million years ago (Mya), followed by rapid evolution and diversification after the Middle Miocene (~13.82 Mya). The analysis of selective pressure indicated that the mitochondrial protein-coding genes (PCGs) of Galloanseres species have predominantly undergone purifying selection. The free-ratio model revealed that the evolutionary rates of COX1 and COX3 were lower than those of the other PCGs, whereas ND2 and ND6 had faster evolutionary rates. The CmC model also indicated that most PCGs in Anseriformes exhibited stronger selective constraints. Our study suggests that the distinct evolutionary trends and energy requirements of Galliformes and Anseriformes drive different evolutionary patterns in the mitogenome.

Mitochondrial Genomic Evidence of Selective Constraints in Small-Bodied Terrestrial Cetartiodactyla

Body size may drive the molecular evolution of mitochondrial genes in response to changes in energy requirements across species of different sizes. In this study, we perform selection pressure analysis and phylogenetic independent contrasts (PIC) to investigate the association between molecular evolution of mitochondrial genome protein-coding genes (mtDNA PCGs) and body size in terrestrial Cetartiodactyla. Employing selection pressure analysis, we observe that the average non-synonymous/synonymous substitution rate ratio (ω) of mtDNA PCGs is significantly reduced in small-bodied species relative to their medium and large counterparts. PIC analysis further confirms that ω values are positively correlated with body size (R2 = 0.162, p = 0.0016). Our results suggest that mtDNA PCGs of small-bodied species experience much stronger purifying selection as they need to maintain a heightened metabolic rate. On the other hand, larger-bodied species may face less stringent selective pressures on their mtDNA PCGs, potentially due to reduced relative energy expenditure per unit mass. Furthermore, we identify several genes that undergo positive selection, possibly linked to species adaptation to specific environments. Therefore, despite purifying selection being the predominant force in the evolution of mtDNA PCGs, positive selection can also occur during the process of adaptive evolution.

The branched mitochondrial respiratory chain from the jellyfish Stomolophus sp2 as a probable adaptive response to environmental changes

During their long evolutionary history, jellyfish have faced changes in multiple environmental factors, to which they may selectively fix adaptations, allowing some species to survive and inhabit diverse environments. Previous findings have confirmed the jellyfish's ability to synthesize large ATP amounts, mainly produced by mitochondria, in response to environmental challenges. This study characterized the respiratory chain from the mitochondria of the jellyfish Stomolophus sp2 (previously misidentified as Stomolophus meleagris). The in-gel activity from isolated jellyfish mitochondria confirmed that the mitochondrial respiratory chain contains the four canonical complexes I to IV and F0F1-ATP synthase. Specific additional activity bands, immunodetection, and mass spectrometry identification confirmed the occurrence of four alternative enzymes integrated into a branched mitochondrial respiratory chain of Stomolophus sp2: an alternative oxidase and three dehydrogenases (two NADH type II enzymes and a mitochondrial glycerol-3-phosphate dehydrogenase). The analysis of each transcript sequence, their phylogenetic relationships, and each protein's predicted models confirmed the mitochondrial alternative enzymes' identity and specific characteristics. Although no statistical differences were found among the mean values of transcript abundance of each enzyme in the transcriptomes of jellyfish exposed to three different temperatures, it was confirmed that each gene was expressed at all tested conditions. These first-time reported enzymes in cnidarians suggest the adaptative ability of jellyfish's mitochondria to display rapid metabolic responses, as previously described, to maintain energetic homeostasis and face temperature variations due to climate change.

Comparative mitogenomics of Brachiopods reveals conservatism in articulate species and unusualness in inarticulate species

BACKGROUND

Brachiopods are a phylum of marine invertebrates with over 10,000 fossil species. Today, there are fewer than 500 extant species assigned to the class Articulata or Inarticulata and for which knowledge of evolutionary genetics and genomics is still poor. Until now, complete mitogenome sequences of two inarticulate species and four articulate species were available.

METHODS AND RESULTS

The complete mitogenome of the inarticulate brachiopod species Lingula reevii (20,778 bp) was obtained by using next generation sequencing. It contains 12 protein-coding genes (the annotation of atp8 is unsure), two ribosomal RNA genes, 26 transfer RNA genes, and one supernumerary ORF that is also conserved in the inarticulate species Lingula anatina. It is hypothesized that this ORF could represent a Lingula-specific mtORFan gene (without obvious homology to other genes). Comparative mitogenomics indicate the mitochondrial gene order of L. reevii is unique among brachiopods, and that compared to articulate species, inarticulate species exhibit massive mitogenome rearrangements, deviant ATP8 protein sequences and supernumerary ORFs, possibly representing species- or lineage-specific mtORFan genes.

CONCLUSION

The results of this study enrich genetics knowledge of extant brachiopods, which may eventually help to test hypotheses about their decline.

Taxonomic Status and Phylogenetic Relationship of the Charadriidae Family Based on Complete Mitogenomes

Background

The Charadriiformes provide a good source for researching evolution owing to their diverse distribution, behavior, morphology, and ecology. However, in the Charadrii, family-level relationships remain understudied, and the monophyly of Charadriidae is also a subject of controversy.

Methods

In the present study, we generated complete mitogenomes for two species, Charadrius leschenaultii and Charadrius mongolus, which were found to be 16,905 bp and 16,844 bp in length, respectively. Among the 13 protein codon genes, we observed variation in the rate of non-synonymous substitution rates, with the slowest rate found in COI and the fastest rate observed in ATP8. The Ka/Ks ratio for all Charadriidae species was significantly lower than one, which inferred that the protein-coding genes underwent purifying selection.

Results

Phylogenetic analysis based on the genes of Cyt b, 12S and ND2 revealed that the genus Pluvialis is the sister group of three families (Haematopodidae, Ibidorhynchidae, Recurvirostridae). However, the phylogenetic analysis based on complete mitogenomes indicated that the genus Pluvialis is within the Charadriidae family.

Conclusion

This study highlights the importance of carefully selecting the number of genes used to obtain accurate estimates of the species tree. It also suggests that relying on partial mtDNA genes with fast-evolving rates may lead to misleading results when resolving the Pluvialis sister group. Future research should focus on sequencing more mitogenomes at different taxonomic levels to gain a better understanding of the features and phylogenetic relationships within the Charadriiformes order.

Comparative analysis of the mitochondrial genomes of the family Mactridae (Mollusca: Venerida) and their phylogenetic implications

Taxonomy and phylogenetic relationships within the family Mactridae have remained debatable because of the plasticity of morphological characteristics and the lack of accurate molecular data, thereby resulting in abundant synonyms and taxa rearrangements. Mitochondrial genomes (mitogenomes) have been widely used as powerful tools to reconstruct phylogenies of various groups of mollusks; however, they have not been used for studying the phylogeny of mactrids specifically. In the present study, mitogenomes of seven Mactridae species, namely Mactra chinensis, Mactra cygnus, Mactra quadrangularis, Mactra cumingii, Mactrinula dolabrata, Raeta pulchella, and Raeta sp., were sequenced by Illumina high-throughput sequencing, and a comparative mitochondrial genomic analysis was conducted. The newly sequenced mitogenomes were double-stranded circular molecules, with all functional genes encoded on the heavy strand. All the new mactrid mitogenomes had two rRNA genes (12S and 16S), 13 protein-coding genes (PCGs) (atp6, cox1, cox2, cox3, cytb, nad1, nad2, nad3, nad4, nad4l, nad5, nad6, and atp8), and 22 tRNAs. The mitogenomes showed considerable variation in AT content, GC skew, and AT skew. The results of the phylogenetic analysis confirmed monophyly of the family Mactridae and suggested that genera Mactrinula, Spisula, Rangia, and Mulinia should not be placed under subfamily Mactrinae. Our results supported that potential cryptic species existed in Mactra antiquata. We also proposed subfamily Kymatoxinae should belong to the family Mactridae rather than Anatinellidae and Mactra alta in China should be Mactra cygnus. Additionally, conservation in functional gene arrangement was found in genera Mactra, Raeta, and Lutraria. But gene orders in S. sachalinensis and S. solida were quite different, questioning their congeneric relationship. Our results further suggested that the taxonomy within the family Mactridae requires an integrative revision.

The Complete Mitochondrial Genomes of Two Rock Scallops (Bivalvia: Spondylidae) Indicate Extensive Gene Rearrangements and Adaptive Evolution Compared with Pectinidae

Different from the diverse family Pectinidae, the Spondylidae is a small group with a single genus that shares the sedentary life habit of cementing themselves to the substrate. However, little information related to the genetic diversity of Spondylidae has been reported. In the present study, the complete mitochondrial genomes of Spondylus versicolor and S. spinosus were sequenced and compared with those of pectinids. The mtDNA of S. versicolor and S. spinosus show similar patterns with respect to genome size, AT content, AT skew, GC skew, and codon usage, and their mitogenomic sizes are longer than most pectinid species. The mtDNA of S. spinosus is 27,566 bp in length, encoding 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes, while an additional tRNA-Met was found in the mtDNA of S. versicolor, which is 28,600 bp in length. The monophylies of Spondylidae and Pectinidae were well supported, but the internal relationships within Pectinidae remain unresolved due to the paraphyly of the genus Mimachlamy and the controversial position of the tribe Aequipectinini. The gene orders of S. versicolor and S. spinosus are almost identical but differ greatly from species of the Pectinidae, indicating extensive gene rearrangements compared with Pectinidae. Positive selection analysis revealed evidence of adaptive evolution in the branch of Spondylidae. The present study could provide important information with which to understand the evolutionary progress of the diverse and economically significant marine bivalve Pectinoidea.

Comprehensive analysis of the complete mitochondrial genomes of three Coptis species (C. chinensis, C. deltoidea and C. omeiensis): the important medicinal plants in China

Coptis plants (Ranunculaceae) contain high levels of isoquinoline alkaloids and have a long history of medicinal use. Coptis species are of great value in pharmaceutical industries and scientific research. Mitochondria are considered as one of the central units for receiving stress signals and arranging immediate responses. Comprehensive characterizations of plant mitogenomes are imperative for revealing the relationship between mitochondria, elucidating biological functions of mitochondria and understanding the environmental adaptation mechanisms of plants. Here, the mitochondrial genomes of C. chinensis, C. deltoidea and C. omeiensis were assembled through the Nanopore and Illumina sequencing platform for the first time. The genome organization, gene number, RNA editing sites, repeat sequences, gene migration from chloroplast to mitochondria were compared. The mitogenomes of C. chinensis, C. deltoidea and C. omeiensis have six, two, two circular-mapping molecules with the total length of 1,425,403 bp, 1,520,338 bp and 1,152,812 bp, respectively. The complete mitogenomes harbors 68-86 predicted functional genes including 39-51 PCGs, 26-35 tRNAs and 2-5 rRNAs. C. deltoidea mitogenome host the most abundant repeat sequences, while C. chinensis mitogenome has the largest number of transferred fragments from its chloroplasts. The large repeat sequences and foreign sequences in the mitochondrial genomes of Coptis species were related to substantial rearrangements, changes in relative position of genes and multiple copy genes. Further comparative analysis illustrated that the PCGs under selected pressure in mitochondrial genomes of the three Coptis species mainly belong to the mitochondrial complex I (NADH dehydrogenase). Heat stress adversely affected the mitochondrial complex I and V, antioxidant enzyme system, ROS accumulation and ATP production of the three Coptis species. The activation of antioxidant enzymes, increase of T-AOC and maintenance of low ROS accumulation in C. chinensis under heat stress were suggested as the factors for its thermal acclimation and normal growth at lower altitudes. This study provides comprehensive information on the Coptis mitogenomes and is of great importance to elucidate the mitochondrial functions, understand the different thermal acclimation mechanisms of Coptis plants, and breed heat-tolerant varieties.

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.