# Expanding the Search for Sperm Transmission Elements in the Mitochondrial Genomes of Bivalve Mollusks

Donald Stewart, Brent Robicheau, Noor Youssef, Manuel Garrido-Ramos,

Emily Chase, Sophie Breton

► To cite this version:

Donald Stewart, Brent Robicheau, Noor Youssef, Manuel Garrido-Ramos, Emily Chase, et al. Expanding the Search for Sperm Transmission Elements in the Mitochondrial Genomes of Bivalve Mollusks. Genes, 2021, 12 (8), pp.1211. 10.3390/genes12081211. hal-03663192

HAL Id: hal-03663192

https://amu.hal.science/hal-03663192

Submitted on 10 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Expanding the Search for Sperm Transmission Elements in the Mitochondrial Genomes of Bivalve Mollusks

Donald T. Stewart ¹,*,† D, Brent M. Robicheau ²,†, Noor Youssef ², Manuel A. Garrido-Ramos ³ D, Emily E. Chase ⁴,⁵
and Sophie Breton ⁶

¹ Department of Biology, Acadia University, Wolfville, NS B4P 2R6, Canada

² Department of Biology, Dalhousie University, Halifax, NS B3H 4R2, Canada;
brentrobicheau@acadiau.ca (B.M.R.); N.Youssef@dal.ca (N.Y.)

³ Departamento de Genética, Facultad de Ciencias, Universidad de Granada, 18071 Granada, Spain;
mgarrido@ugr.es

⁴ Institut Méditerranéen d'Océanologie, Aix-Marseille University, 13288 Marseille, France;
emilychase@acadiau.ca

5 Institut Hospitalo-Universitaire (IHU) Méditerranée Infection, 19-21 Boulevard Jean Moulin,
13005 Marseille, France

6 Département de Sciences Biologiques, Université de Montréal, Montréal, QC H3C 3J7, Canada;
s.breton@umontreal.ca

* Correspondence: don.stewart@acadiau.ca

† Authors contributed equally.

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Stewart, D.T.; Robicheau, B.M.; Youssef, N.; Garrido-Ramos, M.A.; Chase, E.E.; Breton, S. Expanding the Search for Sperm Transmission Elements in the Mitochondrial Genomes of Bivalve Mollusks. *Genes* **2021**, *12*, 1211.
https://doi.org/10.3390/genes12081211

Academic Editor: Antonio Figueras

Received: 25 June 2021

Accepted: 3 August 2021

Published: 5 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA) in bivalve mollusks is one of the most notable departures from the paradigm of strict maternal inheritance of mtDNA among metazoans. Recently, work on the Mediterranean mussel *Mytilus galloprovincialis* suggested that a nucleotide motif in the control region of this species, known as the sperm transmission element (STE), helps protect male-transmitted mitochondria from destruction during spermatogenesis. Subsequent studies found similar, yet divergent, STE motifs in other marine mussels. Here, we extend the in silico search for mtDNA signatures resembling known STEs. This search is carried out for the large unassigned regions of 157 complete mitochondrial genomes from within the Mytiloida, Veneroida, Unionoida, and Ostreoida bivalve orders. Based on a sliding window approach, we present evidence that there are additional putative STE signatures in the large unassigned regions of several marine clams and freshwater mussels with DUI. We discuss the implications of this finding for interpreting the origin of doubly uniparental inheritance in ancestral bivalve mollusks, as well as potential future in vitro and in silico studies that could further refine our understanding of the early evolution of this unusual system of mtDNA inheritance.

**Keywords:** mitochondrial DNA inheritance; sperm transmission elements; bivalves; DUI

## 1. Introduction

Over 100 species of bivalved mollusks are known to exhibit the unusual system of doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA) [1]. DUI is an inheritance pattern characterized by two sex-associated lineages by which mtDNA is transmitted, a male lineage and a female lineage [2]. Typically, the male-transmitted mtDNA type (or M-type) is the exclusive type present in sperm whereas the female-transmitted type (or F-type) is the exclusive (or at least vastly predominant) type present in eggs and in somatic tissues of females [3]. The F-type is also the more abundant type in male somatic tissues, whereas the M-type is typically present in lower amounts [4]. DUI was first described in marine mussels (genus *Mytilus*) approximately 30 years ago [5,6]. In addition to marine mussels, DUI has been found in marine clams of the order Veneroida, freshwater mussels of the order Unionoida, and nutshell of the order Nuculanoida; however, the phylogenetic position of Nuculanoida (e.g., *Yoldia*) within or outside of the Protobranchia
is controversial [1,7]. In a recent review of current research in the area of DUI, Stewart and colleagues [8] discussed the following topics: whether DUI evolved once or multiple origins, the challenges of annotating mtDNA genomes of DUI species, whether there is a link between the male- and female-transmitted mtDNA genomes and sex determination, what possible role(s) non-oxidative phosphorylation mitochondrial genes (i.e., open reading frames or ORFs) may play in freshwater mussels, and whether or not there are obvious evolutionary advantages of DUI, among other questions. One of the other points of discussion briefly considered was the possible function of conserved sequence motifs and sperm transmission elements (STEs) in the mitochondrial genomes of marine mussels. In this paper, we focus on this point by extending the analysis of Robicheau et al. [9] to search for conserved DNA sequence motifs that could represent possible sperm transmission elements as first described in Kyriakou et al. [10].

Within mitochondrial genomes, the search for molecular signals that may play a role in their sex-associated mode of transmission is facilitated by a particularly curious phenomenon best studied in marine mussels. As directly documented in several species within the blue mussel genus *Mytilus*, and inferred in several other members of the Mytilidae through genome and phylogenetic analyses, the male- and female-transmitted genomes in DUI positive taxa occasionally recombine (reviewed in [11]). While all manner of combinations of recombinant genomes are likely possibly, several known examples consist of recombinants composed of a mostly F-type genome (including the F-type protein coding genes, 12S and 16S ribosomal RNA genes, and most or all F-type transfer RNA genes) and an inserted control region from an M-type genome (e.g., [12]). Although there is variation in terms of exactly which additional adjacent portions of the male genome may also be recombined (e.g., [13]), this type of recombination event can result in a mostly female-type mitochondrial genome that is transmitted in a father-son pathway via sperm [14]. Such recombinant genomes have been referred to as recently masculinized (RM) genomes and the event is described as a “role-reversal,” which refers to the switch in sex-associated transmission for the majority of genes in these recombinant mitochondrial genomes [11]. An obvious conclusion to be drawn from the study of these recombinant genomes is that there could be a molecular signature (for e.g., a DNA motif) present in the M-type control region that facilitates its transmission from male parent to male offspring.

As summarized by Kyriakou et al. [10], there are three possible mechanisms that could result in the M-type dominating transmission in sperm while the F-type remains the predominant type found in somatic tissues of males: (i) only mitochondria from the male parent enter the primordial germ cells (PCGs) that will go on to populate the spermatogonia; (ii) mitochondria from both parents enter the PCGs but the M-type mitochondria experience more rapid replication and, in effect, outcompete the F-type mitochondria in this arena; and (iii) the F-type mitochondria are actively eliminated from the PCGs leaving only the M-type, which must also have entered the relevant blastomere in early development. Using electrophoretic mobility shift assays of gonad extracts from Mediterranean sea mussels (*Mytilus galloprovincialis*), Kyriakou et al. [10] provided evidence that there is a DNA sequence motif in the control region of the M-type that forms a DNA-protein complex. These authors also demonstrated that this DNA-protein interaction occurs in the perinuclear space within male gonad cells, and they proposed that these perinuclear mitochondria are protected from general degradation of cytosolic mitochondria that results in the destruction of the F-type mitochondria and their mtDNA genomes. Kyriakou and colleagues thus argued that the third scenario above, i.e., the destruction of the F-type and the preferential survival of the M-type, is the most likely of the three possible mechanisms underpinning the DUI system, at least in marine mussels [10].

Relatively little experimental or in silico work has been performed on mytilids over the past several years to follow up on the aforementioned work by Kyriakou et al. [10]. Robicheau et al. [9] used the STE signature of Kyriakou et al. [10], a 22-23 bp conserved DNA sequence string, to search for putative STE motifs in the highly divergent male and female-transmitted mitochondrial genomes of the horse mussel, *Modiolus modiolus* [15].
Overall, this work suggested that divergent, yet still recognizable, STEs may be found when less stringent search parameters are implemented during alignment procedures (e.g., by converting sequences to simply purines and pyrimidines). Indeed, a putative STE in the control region of the male mitochondrial genome of *M. modiolus* was identified using this more relaxed search approach [9]. Hence, the identification of homologous DNA sequence motifs that could be involved in the protection and (or) sequestration of the M-type mitochondrial genomes remains a feasible area for discovery despite the STEs' rather short nucleotide signature.

In addition to the obvious relevance for understanding the evolution and maintenance of DUI in bivalve mollusks, studies of the molecular signatures involved in precisely controlling the nature and extent of mitochondrial heteroplasmy in animal cells have broad applicability in cell and molecular biology as noted by [16]. Mitochondrial heteroplasmy is implicated in human diseases, some quite severe, affecting a number of organ systems [17,18]. Teasing apart the mechanism(s) of the mitochondrial dynamics of the male- and female-transmitted mitochondrial genomes in bivalves could shed light on pathways of mitochondrial biogenesis (e.g., of the role of the nuclear-encoded gene PCG-1α [19,20] and mitophagy (e.g., the PINK1/Parkin axis that ultimately activates mitophagy of select mitochondria in animals [21]).

In this paper, we extend the search for putative STE signatures throughout the Bivalvia using a sliding window approach, reminiscent of the methodologies used in Robicheau et al. [9], to scan the large unassigned regions of available whole mitochondrial genomes for various bivalves. Collectively, our analyses provide evidence of putative STE motifs in the large unassigned regions of all three major lineages of bivalves that exhibit DUI, the Mytiloida, the Veneroida, and the Unionoida. Based on this finding, we then discuss strategies for further characterizing these molecular signatures/structures moving forward, as well as the broader implications of finding STE-like signatures throughout these three major lineages with respect to the evolution of doubly uniparental inheritance.

## 2. Materials and Methods

The definitive/known "STE motif" corresponds to the 22 bp nucleotide string reported by Kyriakou et al. [10] for an *M. galloprovincialis* RM-type that was shown to be functional via EMSAs: its sequence is 1bp-CCATAAATGTTTGAAAATAAGG-22bp. A diagram reiterating the homology between the definitive STE motif to other known STE signatures within closely related mytilids, as was originally proposed by Kyriakou et al. [10], is given in Figure 1a. As discussed elsewhere [8], the extremely short nucleotide length of the STE motif causes a significant challenge for identifying homologous STE signatures across Bivalvia. As an extreme example, the similarity between a previously described *M. trossulus* male-type STE versus the functional/known *M. galloprovincialis* RM-type STE was reported as ~63% (Figure 1a; also see [10] and the figures given therein). To overcome the issue of short sequence similarities, in this study, we made use of a sliding window approach implemented in Geneious Prime 2020.1.2 [22] to conduct a widespread search for putative STEs. The sliding window procedure in Geneious is based on EMBOSS dot [23]. From within Geneious, the sliding window parameters selected included a "high sensitivity/slow" algorithm and an exact scoring matrix that assumed no ambiguous matches (e.g., C vs. Y; [22]). The tile size for visualizations was set to 5000, and the sliding window size was set to the entire STE motif length of 22bp. The threshold selected for returning a positive match to the known STE motif was manually determined through calibration using known mytilid STE sequences along with their complementary large unassigned regions (LURs) from the sequences given in Kyriakou et al. [10]. We refer to this "threshold selected" value as Ts. The relationship between the sliding window Ts value and the number of observed matches (or alignment hits) during the sliding window procedure is shown in Figure 1b. The relationship is such that lower T values will lead to more spurious matches and potentially false positives, while higher T values will lead to fewer spurious matches but potential false negatives. Using the definitive
STE motif (STE-Mgallo) and with prior knowledge that the quantity of STE alignment matches in *M. galloprovincialis* RM (DQ399833) should be *n* = 3, and in *M. galloprovincialis* M (AY363687) and *M. trossulus* M (GU936626) the STE quantity should be *n* = 1 each [10], a final threshold of *Ts* = 56 was selected. This threshold would allow for the observation of putative STEs in silico using the 22 bp STE motif, while also attempting to minimize false positive matches. We also found that the homologous 24 bp STE signature in the *M. trossulus* M mitotype (STE-Mtross) was unable to capture the definitive STE motif and its known duplications in *M. gallprovincialis* RM mitotype [10] when the sliding window was equivalent to the STE-Mtross motif length and above *T* = 57. Furthermore, although comparison to *M. gallprovincialis* M and *M. trossulus* M would suggest that a lower *T* value could have been selected for the known 22 bp STE motif (STE-Mgallo), we decided to select a more conservative *T* value under the assumption that it would be more informative to identify fewer potentially "true" STE matches, as opposed to capturing more numerous dissimilar STEs that may also have a higher likelihood of being false positives. Hence, we assume that our approach would still be unable to capture extremely ancient STE motifs that may fall below the threshold selected.

Given that STEs have already been discussed elsewhere for Mytiloida, we aimed to screen as many Veneroida and Unionoida bivalves as possible for putative STE motifs. However, to provide additional context we also screened twenty mitochondrial genomes for Ostreoida, as well as several mitochondrial genomes for Mytiloida. Complete mitochondrial genomes were retrieved from NCBI [24] using the Genome resource database therein (https://www.ncbi.nlm.nih.gov/genome/; accessed on 24 May 2021). Additional mitochondrial genomes (especially for M-types) were also retrieved from the nucleotide database of NCBI [24,25]. As best as possible, mitotype sex was inferred from NCBI sequence descriptors, for example, sometimes a DUI mitotype is specifically listed for a given sequence, while in other instances, tissue-type can suggest that a particular mt genome is an F-type based on the general principle that female-transmitted mtDNA is typically present in somatic tissues of both sexes [3,4]. Any mt genomes with an unresolved or uncertain transmission route based on the sequence descriptors are denoted by question marks. Overall, 157 mitochondrial genomes were retrieved for analysis (Supplementary Table S1). To facilitate the manual search of STE matches using the sliding window approach described above, STE motif searches were first limited to the large unassigned regions in each genome—these were manually extracted from complete mtDNA from within Geneious [22]. LURs are the positions most likely to house mitochondrial control regions [26] and are known to be implemented in the formation of recently masculinized mitotypes via mitochondrial recombination events [11]. The STE identified by Kyriakou et al. [10] was also found within a large unassigned control region. The positions of putative LURs were taken as the longest stretches of unannotated DNA between genes identified as either protein-coding (oxidative phosphorylation and ORF proteins) or ribosomal RNA. If the exact position of a DUI-related H-/F-/M-ORF was uncertain due to annotation inconsistencies (e.g., in *Unio pictorum*, the F-ORF was annotated as a "miscellaneous feature" rather than a gene), then it was not factored into the LUR boundaries. Transfer RNAs were not masked during the downstream sliding window analysis, nor were they excluded in the assessment of LURs. In total, 215 LURs were screened against the known 22 bp STE motif (STE-Mgallo). Any genome found to contain a putative STE motif within a LUR was then re-screened; for this second iteration, STE motifs were searched for within the entire circular mt genome. Final gene maps showing putative STE motif positions relative to other mitochondrial genes were generated in Geneious Prime 2020.1.2 [22], and then manually assembled into a composite display. In addition, the 30 bp upstream from putative STEs, as well as the putative STE motifs themselves (for those associated with LURs) were further retrieved (via Geneious) and inspected for sequence conservation. The 30 bp upstream was left unaligned and adenines highlighted; meanwhile, the putative STE nucleotides were aligned with the known STE motif (STE-Mgallo) using the MUSCLE algorithm [27] and then conserved nucleotide positions highlighted. The sequence similarity between
each putative STE and the known 22 bp STE motif was further plotted from this alignment using the *ggplot2* [28] and *ggrepel* [29] packages in R version 4.0.3 [30] via RStudio version 1.3.1093 [31] (note *ggplot2* was also used for Figure 1b).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Known sperm transmission elements (STEs) from Kyriakou et al. [10]. (b) Calibration of the sliding window procedure used to identify novel putative STE motifs. (c) Alignments of STE motifs from Kyriakou et al. [10] (accessions shown in shaded grey boxes) versus additional putative STE motifs identified in this study; the characterized STE motif from Kyriakou et al. [10] is marked with a ★. Nucleotides of the STE motif region were aligned using a MUSCLE algorithm [27], while the 30 bp upstream remains unaligned and simply shows adenine residues. H = Hermaphroditic mitotype, F = Female mitotype, M = Male mitotype, RM = Recently masculinized mitotype.</figcaption></figure>
The distribution of putative STEs that could be attributed to LURs was further displayed on a bivalve phylogeny reconstructed in MEGAX version 10.1.8 [32,33] for the mitochondrial genomes screened. Since our objective was to visually display the putative STE motifs across the taxonomic orders investigated, as opposed to conducting a comprehensive Bivalvia evolutionary tree, we chose to use the simpler Neighbour-joining phylogenetic reconstruction method, as well as p-distances and pair-wise deletion [32,33]. The tree was also bootstrapped using 1000 replicates and made use of complete *col* sequences that were retrieved from each of the complete mitochondrial genomes indicated earlier (accomplished through Geneious; [22]). Furthermore, said *col* sequences were aligned via MUSCLE [27] and then trimmed to a conserved start/end nucleotide prior to constructing the phylogeny (the alignment and sequence trimming were both completed within MEGAX [32,33]). Onto this phylogeny, the major bivalve clades, as well as the lineages with putative STEs present, were shown through additional/manual annotation.

Given the small size of the known STE motif (22bp), the final stage of our analysis was to assess whether the putative STEs identified herein are likely to represent biologically significant matches between the known STE motif and LURs. This was accomplished by comparing the known STE motif to pseudo-sequences that were randomly generated. To this end, we generated 1000 random nucleotide sequences with lengths equivalent to the length of the LUR in the original analysis. The random sequences were generated using the empirical frequencies estimated from the original LUR (for either the forward or reverse DNA strand depending on putative STE orientation). For comparison, the analysis was repeated assuming equal nucleotide frequencies (i.e., 0.25 each). For each randomly generated sequence, we used a sliding window of 22 positions to identify potential STE motifs. A match was identified if the similarity between the window and the known STE was greater than or equal to the percent similarity observed for a newly identified putative STE motif (i.e., calculating potential false positive matches). Multiple matches could potentially be identified within a single sequence. We therefore report both the total number of false STE motifs found and the number of sequences with ≥1 false matches. This procedure allows us to estimate a pseudo *p*-value of the number of motif matches we would expect to find by random chance and in the absence of any duplication events; the estimate is calculated as the proportion of randomly generated sequences (out of 1000 replicates) with at least one motif match. This is similar to the ideas presented in [34]. Overall, the probabilities reported allow for a heuristic sense of whether the putative STE motifs identified during our analysis were likely to have originated from a mitochondrial region simply prone to displaying random matches to the known STE motif, or whether matches are more likely to have arisen from DUI relevant homologies. The statistical analyses above were carried out in Python [35] using the *NumPy* package [36] (see Supplemental Methods 1 for code).

**3. Results and Discussion**

The first extensive analysis of putative sperm transmission elements (STEs) in bivalve mollusks with the DUI system for mitochondrial DNA comes from the study by Kyriakou et al. [10], which focused on *Mytilus* spp. Accordingly, we used the 22 bp motif identified by those authors as a starting point to search for STEs in other bivalve mitochondrial genomes. As a critical first step, we needed to specify the search parameters to find a similarity threshold range that would minimize both false positive matches (i.e., searches that falsely identify non-homologous regions as matches due to low stringency) as well as minimizing false negative searches (i.e., searches with a stringency set so high that no or virtually no homologous regions were identified). As outlined earlier, a threshold of T = 56 for the known motif signature along with a sliding window length of 22 bp met these criteria (Figure 1a,b).

Using the calibrated sliding window method above, we were able to identify a total of 23 STE-like signatures associated with LURs from the 157 mitochondrial genomes screened. Figure 1c shows the nucleotides identified as having homology to the known STE motif
(pink). The 30 bp upstream from each signature (with adenines highlighted yellow) is also shown in Figure 1c. Out of the total LUR-associated putative STEs, four are previously known from Kyriakou et al. [10] (Map Codes 3, 4, 5, and 8) and one is the known STE motif (★), these are highlighted grey (Figure 1c). The STE signature from *M. trossulus* M (Map Code 8) was manually added to the analysis, as its relatively low sequence similarity makes it interesting for comparative purposes, yet it could not be identified using our search strategy (hence, this motif represents a true positive that falls below the threshold of detection). The degree of sequence similarity between the putative STE motifs of Figure 1c and the known STE is shown in Figure 2. Putative motifs range in sequence similarity from 62.5% (15/24bp) to 100% or 22/22bp.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Degree of sequence similarity between STE-like signatures and the known 22 bp STE motif from Kyriakou et al. [10]. Nucleotides were aligned using MUSCLE [27].</figcaption></figure>

We further investigated whether the mitochondrial genomes housing putative STE motifs in their LURs also displayed evidence of additional STE homology to other parts of their mitochondrial genomes (i.e., outside the LURs). The full genome maps for protein-coding genes, tRNAs, and putative STE motif signatures are given (Figure 3). As illustrated, for many of the genomes [12/20 mtDNAs] these motifs were only localized to the unassigned/control region (Figure 3). Interestingly, for the eight genomes with additional motifs outside the LURs, these were in many cases on the complementary DNA strand to the protein-coding or rRNA gene they were identified in and were usually associated with either 16S (often), as well as *nd5*, *atp6*, or *co1* (less often; Figure 3). The localization of many putative STE motifs to LURs is interesting given that nearly two decades ago, Burzyński et al. [37], Zbawicka et al. [38] and Cao et al. [26] suggested that genetic signatures associated with DUI signaling were probably to be found in the non-coding portion of mtDNA as these unassigned regions or control regions, as they are variously called, are the mtDNA regions primarily associated with recombination events leading to the formation of masculinized mitochondrial genomes in the family Mytilidae.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Mitochondrial genome maps for bivalve species exhibiting putative STE motifs localized to large unassigned regions in mtDNA. STE-like signatures are identified by map codes within boxes (these correspond to codes in Figure 1c) and can either be found on the Reverse (Rev) or Forward (For) mtDNA orientation. *M. gallo* refers to *M. galloprovincialis*.</figcaption></figure>
The phylogenetic distribution of bivalve mitotypes exhibiting putative STE signatures within their LURs is shown in the tree diagram presented in Figure 4. It is important to note that the Neighbour-joining tree in Figure 4 is based on the alignment of complete *col* sequences. Although the tree is computationally quite simple, it nevertheless does an adequate job of representing the known broader evolutionary relationships previously indicated in the literature [1,39]. Overall, the phylogeny indicates that several distantly related clades from across the Unionida, Veneroida and Mytiloida have signatures similar to the known 22 bp STE motif (Figure 4). For those in Unionoida and Veneroida, signatures were found in thirteen species: three in DUI paternal genomes, two in DUI maternal genomes, one in an SMI hermaphroditic genome, and seven in genomes of unknown inheritance pattern for which we suppose are maternally-transmitted mtDNA (Figure 1). Despite our widespread finding of this signature, we did not observe their presence among taxa as a monophyletic group, nor were the signatures limited to male-transmitted genomes (Figure 4). For both trends, we must consider more closely that recombination between DUI mitotypes may obfuscate evolutionary patterns. As indicated earlier, recombination and masculinization of mtDNA genomes were demonstrated in the Mytilidae [13,40,41], and recombination and role-reversal have been events implicated in the other DUI positive bivalve orders [7,11]. Some of the male-transmitted mitochondrial genomes that were used during our analysis were known to be chimeric and, indeed, all M-types could be chimeric or derived from a chimeric ancestor. This means that the putative STE motifs are being mapped onto a genome tree rather than a gene tree. If we could confidently determine that recombination had occurred repeatedly and could reconstruct the recombination events throughout the hundreds of millions of years of bivalve evolutionary history, then we could determine whether the STE motif is monophyletic. Monophyly of the STE motifs across all bivalves would be strong evidence that the phenomenon evolved once. However, given the extremely long time periods involved, and the extremely small length of the putative STEs, reconstructing the evolutionary history of such events with greater confidence than we present here is highly unlikely. On a similar note, it is also possible that the putative STE motifs localized to F-type genomes could be the product of failed masculinization events (as suggested by the mt genomes analyzed in [13,42]), or alternatively, they may instead represent non-functional STEs that have yet to acquire enough nucleotide differences via mutation to appear dissimilar to the functional STE (as implied by sequences in Kyriakou et al. [10] for the *M. galloprovincialis* F-type genome). Accordingly, the presence of STE-like signatures in F-types does not allow us to definitively argue that DUI arose only once, in an ancestor to the Mytiloida, Unionoida, and Veneroida orders. If the DNA sequence motif of the STE was only found in sperm transmitted genomes and could be shown to be present in the Mytiloida, Unionoida, and Veneroida, then the most reasonable explanation would be that the phenomenon had evolved only once. While this is still the most parsimonious explanation, it is possible that the M-type sequestration mechanism evolved independently three or more times in the history of the Bivalvia. Analysis of additional complete sequences of M and F mitochondrial genomes in the future should help test hypotheses of putative cases of recombination and masculinization within particular bivalve orders, which should, in turn, shed more light on the roles of STEs and uncharacterized mitochondrial open reading frames (ORFs) in the origin and maintenance of DUI across the Bivalvia [8].

Despite these aforementioned challenges in understanding the evolutionary history of STEs, our search for putative STE motifs amongst a larger set of bivalve mitochondrial genomes has allowed for two major steps of progress in our understanding: (i) our analyses suggest that the STE signature may represent a unifying feature that is present within LURs from all three of the major lineages of bivalves exhibiting DUI (clams, freshwater mussels, and marine mussels) despite its origins still remaining unclear; and (ii) using the putative STEs identified herein, we are now positioned to design targeted electrophoretic mobility shift assays to detect DNA-protein complexes in Unionoida and Veneroida as was done for Mytiloida [10]. The most promising lineage currently appears to be *Limecola balthica* given
its 81% similarity to the 22 bp STE motif (Figure 2). Moving forward, if the same DNA-protein complexes, and localization of perinuclear mitochondria can be demonstrated in the Unionoida as well as the Veneroida, this would be further strong evidence that the DUI phenomenon evolved once, but that there were subsequent role-reversal/recombination events in the ancestors to these two orders of bivalves as well.

<figure><img src="image_7.png" /></figure>

**Figure 4.** Taxonomic distribution of STE-like signatures discussed in this study. The phylogeny is a Neighbor-joining tree with 1000 bootstrap replicates reconstructed in MEGAX [32,33]. Full-length *co1* genes were aligned using MUSCLE [27] and trimmed to a conserved start/end nucleotide. Bootstrap values are shown in the expanded version of this phylogeny in Figure S1. ⭐ = Known STE motif from Kyriakou et al. [10]. The *Hyriopsis schlegelii* male sequence is likely an F-type genome derived from a male individual [43].
Until relatively recent, sex-linked bivalve heteroplasmy was a central tenant of the DUI framework [44]. In DUI-positive species, males are typically heteroplasmic with the M-type dominating in sperm and male gonad, while the F-type dominates in all tissues of females as well as the somatic tissues of males. A recent analysis of the marine mussel *Semimytilus algosus*, which is hermaphroditic, has demonstrated that the sperm-producing gonad located on one half of the body contains one type of mtDNA while the egg-producing gonad located on the opposite side of the body contains another type of mtDNA [45]. This exception to the rule is still consistent with the hypothesis that some signal is present within the sperm-transmitted mitochondrial genome that is associated with a genetic mechanism that influences its fate (see [46], and particularly [47], for more advanced discussions on genetic signatures/models of DUI).

As a final component of our study, we further assessed whether STE matches are expected to represent true homology as opposed to false homology. False homology may arise as a product of random chance when the nucleotides in the short 22 bp motif are found stochastically, or if the frequency of nucleotides in the mtDNA being searched has a high overall similarity to the nucleotide frequency of the STE. The probabilities associated with identifying an STE match of comparable sequence similarity by random chance in the LURs examined are explored in Table 1; as the trends suggest, higher spurious matches are to be expected when the LUR size is larger and if the nucleotide frequencies in random sequence are equivalent to the LUR frequencies. There is also the possibility that a single random sequence may have more than one STE match; however, this is much rarer than individual sequences having at least one match (Table 1). Based on these analyses of 1000 replicates each, most of the novel LUR-associated STE matches except those for *F. mutica* and *M. lyrata* were significant (or near significant for Paphia textile, Arctica islandica, and Sinonovacula constricta). This is because out of 1000 searches there was a ≤5% probability (at *α* = 0.05) that at least one comparable (or better) match would be found by random chance. For less significant STE-like signatures the probability of observing a similar (or better) match at least once in 1000 random sequences was ~5–14% (Table 1). In contrast, we found that several putative STE motifs from Veneroida, and nearly all Unionoida matches, were very unlikely to have a similar (or better) STE signature by random chance; these probabilities were closer to ~0–1.5% (Table 1). Consequently, the random sequence tests suggest that the widespread occurrence of putative STE signatures across the three major DUI-exhibiting lineages is likely a true DUI-relevant feature.

In terms of future studies to further characterize the putative STEs presented herein, in addition to repeating the EMSA and DNA-protein complex analyses pioneered by Kyriakou et al. [10], further in silico analyses of the conserved DNA motifs of the putative STEs from M-type genomes from Veneroida and Unionoida should be focused on characterizing potential DNA bending properties that were hypothesized by Kyriakou and colleagues as crucial to the interaction of this STE motif with a dimerized protein complex. Additionally, in a follow-up paper [48], Kyriakou and colleagues showed that the homologous control region from the F-type mitochondrial genome of *M. galloprovincialis* was not capable of forming a stable protein structure comparable to the M-derived STE likely due in part to differences in the folding parameters. The ideas surrounding secondary structural folding may be even more relevant given that we did not observe any characteristic strings of adenine residues upstream from the putative STE motifs as was observed for various mytilids ([10]; Figure 1c).
<table><caption>**Table 1.** The probability of observing a comparable or better alignment match as those identified in the sliding window analysis. Probabilities are assessed based on generating 1000 random sequences of equal lengths to each Large Unassigned Region (LUR) under investigation. Reported are the total number of matches (Hits) found in the 1000 random sequences, the number of replicates (Seqs) where at least one or more matches were found, and the final probability of each LUR having a similar (or better) match as those found in our earlier analyses (/1000 replicates). Sequences were generated considering both the nucleotide frequency for each LUR ($\pi^*$) or assuming equal frequencies ($\pi_i = 0.25$, for *i* equal to A, C, G, or T). Frequencies are those of the forward (F) or reverse strand (R) depending on the orientation of the putative STE motif. Significant values (bold) are those where the number of replicates with one or more hits is $\le 0.05$.</caption><thead><tr><th rowspan="2">Map Code (Strand)</th><th rowspan="2">ORDER Taxon</th><th rowspan="2">%</th><th rowspan="2">LUR bp</th><th rowspan="2">LUR π<sup>*</sup> = (π<sub>T</sub>, π<sub>C</sub>, π<sub>G</sub>, π<sub>A</sub>)</th><th colspan="3">If π<sup>*</sup></th><th colspan="3">If π<sub>i</sub> = 0.25</th></tr><tr><th>Total Hits</th><th>Seqs. with ≥1 Hit</th><th>/1000</th><th>Total Hits</th><th>Seqs. with ≥1 Hit</th><th>/1000</th></tr></thead><tbody><tr><td colspan="10" style="text-align: left; font-weight: bold;">MYTILOIDA</td></tr><tr><td>1 (F)</td><td>Mytilus galloprovincialis RM</td><td>72.7</td><td>3590</td><td>(28.2, 15.8, 20.9, 35.0)</td><td>85</td><td>81</td><td>0.081</td><td>14</td><td>14</td><td><strong>0.014</strong></td></tr><tr><td>6 (F)</td><td>Mytilus galloprovincialis M</td><td>77.3 <sup>a</sup></td><td>1529</td><td>(30.1, 15.2, 18.9, 35.8)</td><td>7</td><td>7</td><td><strong>0.007</strong></td><td>0</td><td>0</td><td><strong>0.000</strong></td></tr><tr><td>C (F)</td><td>Mytilus galloprovincialis F</td><td>72.7</td><td>1225</td><td>(28.0, 14.6, 25.5, 31.9)</td><td>11</td><td>11</td><td><strong>0.011</strong></td><td>3</td><td>3</td><td><strong>0.003</strong></td></tr><tr><td>9 (F)</td><td>Mytilus trossulus RM</td><td>100</td><td>1070</td><td>(27.9, 15.3, 19.5, 37.3)</td><td>0</td><td>0</td><td><strong>0.000</strong></td><td>0</td><td>0</td><td><strong>0.000</strong></td></tr><tr><td>12 (F)</td><td>Mytilus edulis M</td><td>95.5</td><td>993</td><td>(28.9, 15.7, 19.4, 35.9)</td><td>0</td><td>0</td><td><strong>0.000</strong></td><td>0</td><td>0</td><td><strong>0.000</strong></td></tr><tr><td>D (F)</td><td>Mytilus edulis F</td><td>72.7</td><td>1226</td><td>(28.1, 14.3, 25.1, 32.5)</td><td>17</td><td>17</td><td><strong>0.017</strong></td><td>4</td><td>4</td><td><strong>0.004</strong></td></tr><tr><td colspan="10" style="text-align: left; font-weight: bold;">VENERIDA</td></tr><tr><td>14 (F)</td><td>Sinonovacula constricta</td><td>72.7</td><td>1602</td><td>(33.1, 8.40, 24.7, 33.8)</td><td>36</td><td>35</td><td><strong>0.035</strong></td><td>6</td><td>6</td><td><strong>0.006</strong></td></tr><tr><td>16 (F)</td><td>Paphia textile</td><td>72.7</td><td>1986</td><td>(26.8, 11.9, 25.4, 35.9)</td><td>53</td><td>51</td><td>0.051</td><td>9</td><td>9</td><td><strong>0.009</strong></td></tr><tr><td>17 (R)</td><td>Coelomactra antiquata</td><td>72.7</td><td>1285</td><td>(31.4, 27.6, 13.3, 27.6)</td><td>7</td><td>7</td><td><strong>0.007</strong></td><td>4</td><td>4</td><td><strong>0.004</strong></td></tr><tr><td>19 (F)</td><td>Fulvia mutica</td><td>72.7</td><td>4368</td><td>(30.7, 12.3, 22.1, 34.9)</td><td>108</td><td>106</td><td>0.106</td><td>9</td><td>9</td><td><strong>0.009</strong></td></tr><tr><td>21 (F)</td><td>Meretrix lyrata</td><td>72.7</td><td>4620</td><td>(41.5, 5.80, 20.8, 31.9)</td><td>142</td><td>137</td><td>0.137</td><td>12</td><td>12</td><td><strong>0.012</strong></td></tr><tr><td>23 (F)</td><td>Limecola balthica</td><td>81.8</td><td>3968</td><td>(35.8, 10.6, 20.9, 32.7)</td><td>0</td><td>0</td><td><strong>0.000</strong></td><td>0</td><td>0</td><td><strong>0.000</strong></td></tr><tr><td>A (F)</td><td>Hippopus hippopus</td><td>77.3</td><td>3027</td><td>(33.4, 12.1, 25.7, 28.8)</td><td>4</td><td>4</td><td><strong>0.004</strong></td><td>1</td><td>1</td><td><strong>0.001</strong></td></tr><tr><td>B (F)</td><td>Arctica islandica</td><td>72.7</td><td>1497</td><td>(37.4, 10.2, 17.4, 35.0)</td><td>52</td><td>52</td><td>0.052</td><td>6</td><td>6</td><td><strong>0.006</strong></td></tr><tr><td colspan="10" style="text-align: left; font-weight: bold;">UNIONIDA</td></tr><tr><td>24 (R)</td><td>Leptodea leptodon</td><td>72.7</td><td>830</td><td>(37.5, 8.40, 26.5, 27.6)</td><td>16</td><td>15</td><td><strong>0.015</strong></td><td>3</td><td>3</td><td><strong>0.003</strong></td></tr><tr><td>25 (R)</td><td>Lampilsilis siliquoidea</td><td>72.7</td><td>575</td><td>(34.8, 11.1, 23.8, 30.3)</td><td>10</td><td>10</td><td><strong>0.010</strong></td><td>4</td><td>4</td><td><strong>0.004</strong></td></tr><tr><td>26 (R)</td><td>Toxolasma parvus</td><td>72.7</td><td>647</td><td>(32.0, 10.8, 25.8, 31.4)</td><td>13</td><td>13</td><td><strong>0.013</strong></td><td>0</td><td>0</td><td><strong>0.000</strong></td></tr><tr><td>27 (R)</td><td>Pyganodon grandis</td><td>72.7</td><td>575</td><td>(33.7, 11.7, 19.8, 34.8)</td><td>9</td><td>9</td><td><strong>0.009</strong></td><td>2</td><td>2</td><td><strong>0.002</strong></td></tr><tr><td>29 (R)</td><td>Unio crassus</td><td>72.7</td><td>557</td><td>(31.2, 13.5, 20.3, 35.0)</td><td>12</td><td>12</td><td><strong>0.012</strong></td><td>0</td><td>0</td><td><strong>0.000</strong></td></tr></tbody></table>

## 4. Conclusions

Overall, our analyses have suggested that STE-like motifs appear to be widespread among several species of marine and freshwater mussels as well as clams. In agreement with previous findings in Kyriakou et al. [10], STEs within the mytilids remain an important feature of this group of bivalves. Based on the entire set of putative STEs identified herein, at present *Limecola balthica* appears a particularly promising study species for future functional experiments regarding STE motifs within the Veneroida. Although the analyses we provided in this study were rather extensive (screening many previously published mitochondrial bivalve genomes for STE-like signatures), our survey still understandably remains somewhat focused on maternally inherited mitochondrial genomes given that F genomes are often the ones reported by researchers. Consequently, as has been the case in past years, there is an ongoing need to sequence novel paternally inherited mitochondrial genomes. Finally, there of course remains the exciting prospect of characterizing the precise protein(s) that interact/bind with STE motif(s), and how/if STE signatures may fit into the framework of other DUI relevant bivalve mitochondrial features (e.g., smithRNAs discussed in [49]). Moving forward, future studies aimed at addressing these topics will be important for progressing our overall understanding of DUI in bivalves.

<footnote><sup>a</sup> STE motif is likely 23bp, however, the 22 bp identified by the sliding window was a 17/22 bp match.</footnote>
**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/genes12081211/s1, Figure S1: Expanded view of Figure 4, Table S1: List of all mitochondrial genomes used in this study, Supplemental Methods 1: Code for random sequence tests.

**Author Contributions:** Conceptualization, D.T.S., B.M.R., E.E.C., M.A.G.-R. and S.B.; conducted initial data analyses, B.M.R., D.T.S. and N.Y.; writing—original draft preparation, B.M.R. and D.T.S. All authors further contributed to critically revising the work for important intellectual content. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Discovery grants from the Natural Sciences and Engineering Research Council of Canada (NSERC) to D.T.S. (grant number 217175) and S.B. (grant number 435656). E.E.C. was funded by the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No713750, carried out with the financial support of the Regional Council of Provence-Alpes-Côte d'Azur and with the financial support of the A*MIDEX (n° ANR- 11-IDEX-0001-02), funded by the Investissements d'Avenir project funded by the French Government, managed by the French National Research Agency [ANR]. B.M.R. was supported by an NSERC CGS-D award, Killam Predoctoral Scholarship, and a NS Graduate Scholarship. M.A.G.-R. was supported by a Harrison McCain Visiting Professorship Award at Acadia University from the Harrison McCain Foundation.

**Institutional Review Board Statement:** This study was conducted with existing datasets. Accordingly, no animal care committee approval was required.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The datasets generated in this study represent a form of meta-analysis; original data sources are referenced within the manuscript itself and (or) listed as sequence accession codes provided in Table S1.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Gusman, A.; Lecomte, S.; Stewart, D.T.; Passamonti, M.; Breton, S. Pursuing the quest for better understanding the taxonomic distribution of the system of doubly uniparental inheritance of mtDNA. *PeerJ* **2016**, *4*, e2760. [CrossRef]

2. Zouros, E. Biparental inheritance through uniparental transmission: The Doubly Uniparental Inheritance (DUI) of mitochondrial DNA. *Evol. Biol.* **2013**, *40*, 1–31. [CrossRef]

3. Breton, S.; Beaupré, H.D.; Stewart, D.T.; Hoeh, W.R.; Blier, P.U. The unusual system of doubly uniparental inheritance of mtDNA: Isn't one enough? *Trends Genet.* **2007**, *23*, 465–474. [CrossRef]

4. Garrido-Ramos, M.A.; Stewart, D.T.; Sutherland, B.W.; Zouros, E. The distribution of male-transmitted and female-transmitted mitochondrial DNA types in somatic tissues of blue mussels: Implications for the operation of Doubly Uniparental Inheritance of mitochondrial DNA. *Genome* **1998**, *41*, 818–824. [CrossRef]

5. Fisher, C.; Skibinski, D.O.F. Sex-biased mitochondrial DNA heteroplasmy in the marine mussel *Mytilus*. *Proc. R. Soc. Lond. B Biol. Sci.* **1990**, *242*, 149–156. [CrossRef]

6. Zouros, E.; Freeman, K.R.; Ball, A.O.; Pogson, G.H. Direct evidence for extensive paternal mitochondrial DNA inheritance in the marine mussel *Mytilus*. *Nature* **1992**, *359*, 412–414. [CrossRef] [PubMed]

7. González, V.L.; Andrade, S.C.S.; Bieler, R.; Collins, T.M.; Dunn, C.W.; Mikkelsen, P.M.; Taylor, J.D.; Giribet, G. A phylogenetic backbone for Bivalvia: An RNA-seq approach. *Proc. R. Soc. B Biol. Sci.* **2015**, *282*, 20142332. [CrossRef] [PubMed]

8. Stewart, D.T.; Breton, S.; Chase, E.E.; Robicheau, B.M.; Bettinazzi, S.; Pante, E.; Youssef, N.; Garrido-Ramos, M.A. An unusual evolutionary strategy: The origins, genetic repertoire, and implications of Doubly Uniparental Inheritance of mitochondrial DNA in bivalves. In *Evolutionary Biology—A Transdisciplinary Approach*; Pontarotti, P., Ed.; Springer International Publishing: Cham, Germany, 2020; pp. 301–323. ISBN 978-3-030-57246-4.

9. Robicheau, B.M.; Breton, S.; Stewart, D.T. Sequence motifs associated with paternal transmission of mitochondrial DNA in the horse mussel, *Modiolus Modiolus* (Bivalvia: Mytilidae). *Gene* **2017**, *605*, 32–42. [CrossRef]

10. Kyriakou, E.; Kravariti, L.; Vasilopoulos, T.; Zouros, E.; Rodakis, G.C. A protein binding site in the M mitochondrial genome of *Mytilus galloprovincialis* may be responsible for its paternal transmission. *Gene* **2015**, *562*, 83–94. [CrossRef]

11. Stewart, D.T.; Breton, S.; Blier, P.U.; Hoeh, W.R. Masculinization events and Doubly Uniparental Inheritance of mitochondrial DNA: A model for understanding the evolutionary dynamics of gender-associated mtDNA in mussels. In *Evolutionary Biology*; Pontarotti, P., Ed.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 163–173. ISBN 978-3-642-00951-8.

12. Śmietanka, B.; Burzyński, A. Disruption of Doubly Uniparental Inheritance of mitochondrial DNA associated with hybridization area of European *Mytilus edulis* and *Mytilus trossulus* in Norway. *Mar. Biol.* **2017**, *164*, 209. [CrossRef]
13. Zbawicka, M.; Wenne, R.; Burzyński, A. Mitogenomics of recombinant mitochondrial genomes of Baltic Sea *Mytilus* mussels. *Mol. Genet. Genom.* **2014**, *289*, 1275–1287. [CrossRef]

14. Venetis, C.; Theologidis, I.; Zouros, E.; Rodakis, G.C. A mitochondrial genome with a reversed transmission route in the Mediterranean Mussel *Mytilus galloprovincialis*. *Gene* **2007**, *406*, 79–90. [CrossRef] [PubMed]

15. Robicheau, B.M.; Powell, A.E.; Del Bel, L.; Breton, S.; Stewart, D.T. Evidence for extreme sequence divergence between the male- and female-transmitted mitochondrial genomes in the bivalve mollusc, *Modiolus modiolus* (Mytilidae). *J. Zool. Syst. Evol. Res.* **2017**, *55*, 89–97. [CrossRef]

16. Zouros, E.; Rodakis, G.C. Doubly Uniparental Inheritance of mtDNA: An unappreciated defiance of a general rule. In *Cellular and Molecular Basis of Mitochondrial Inheritance: Mitochondrial Disease and Fitness*; Sutovsky, P., Ed.; Advances in Anatomy, Embryology and Cell Biology; Springer International Publishing: Cham, Germany, 2019; pp. 25–49. ISBN 978-3-030-04570-8.

17. Nandakumar, P.; Tian, C.; O’Connell, J.; 23andMe Research Team; Hinds, D.; Paterson, A.D.; Sondheimer, N. Nuclear genome-wide associations with mitochondrial heteroplasmy. *Sci. Adv.* **2021**, *7*, eabe7520. [CrossRef]

18. Stewart, J.B.; Chinnery, P.F. The dynamics of mitochondrial DNA heteroplasmy: Implications for human health and disease. *Nat. Rev. Genet.* **2015**, *16*, 530–542. [CrossRef]

19. Garcia, I.; Jones, E.; Ramos, M.; Innis-Whitehouse, W.; Gilkerson, R. The little big genome: The organization of mitochondrial DNA. *Front. Biosci. Landmark Ed.* **2017**, *22*, 710–721.

20. Puigserver, P.; Wu, Z.; Park, C.W.; Graves, R.; Wright, M.; Spiegelman, B.M. A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. *Cell* **1998**, *92*, 829–839. [CrossRef]

21. Carelli, V.; Maresca, A.; Caporali, L.; Trifunov, S.; Zanna, C.; Rugolo, M. Mitochondria: Biogenesis and mitophagy balance in segregation and clonal expansion of mitochondrial DNA mutations. *Int. J. Biochem. Cell Biol.* **2015**, *63*, 21–24. [CrossRef]

22. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. *Bioinformatics* **2012**, *28*, 1647–1649. [CrossRef]

23. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. *Nucleic Acids Res.* **2019**, *47*, W636–W641. [CrossRef]

24. NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. *Nucleic Acids Res.* **2018**, *46*, D8–D13. [CrossRef]

25. Benson, D.A.; Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. *Nucleic Acids Res.* **2015**, *43*, D30–D35. [CrossRef]

26. Cao, L.; Kenchington, E.; Zouros, E.; Rodakis, G.C. Evidence that the large noncoding sequence is the main control region of maternally and paternally transmitted mitochondrial genomes of the marine mussel (*Mytilus* spp.). *Genetics* **2004**, *167*, 835–850. [CrossRef]

27. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* **2004**, *32*, 1792–1797. [CrossRef] [PubMed]

28. Wickham, H. *ggplot2: Elegant Graphics for Data Analysis*; Springer International Publishing: Cham, Switzerland, 2016.

29. Slowikowski, K.; Schep, A.; Hughes, S.; Lukauskas, S.; Irisson, J.-O.; Kamvar, Z.N.; Ryan, T.; Christophe, D.; Hiroaki, Y.; Gramme, P. *Ggrepel: Automatically Position Non-Overlapping Text Labels with 'ggplot2'*. R Package. 2018. Available online: https://rdrr.io/github/slowkow/ggrepel/ (accessed on 25 June 2021).

30. R Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2020. Available online: https://www.R-project.org/ (accessed on 31 May 2021).

31. RStudio Team. *RStudio: Integrated Development Environment for R*. RStudio; PBC: Boston, MA, USA, 2020. Available online: http://www.rstudio.com/ (accessed on 31 May 2021).

32. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549. [CrossRef] [PubMed]

33. Stecher, G.; Tamura, K.; Kumar, S. Molecular Evolutionary Genetics Analysis (MEGA) for MacOS. *Mol. Biol. Evol.* **2020**, *37*, 1237–1239. [CrossRef]

34. Noble, W.S. How does multiple testing correction work? *Nat. Biotechnol.* **2009**, *27*, 1135–1137. [CrossRef] [PubMed]

35. Welcome to Python.Org. Available online: https://www.python.org/ (accessed on 1 June 2021).

36. Harris, C.R.; Millman, K.J.; van der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array programming with NumPy. *Nature* **2020**, *585*, 357–362. [CrossRef]

37. Burzyński, A.; Zbawicka, M.; Skibinski, D.O.F.; Wenne, R. Evidence for recombination of mtDNA in the Marine mussel *Mytilus trossulus* from the Baltic. *Mol. Biol. Evol.* **2003**, *20*, 388–392. [CrossRef]

38. Zbawicka, M.; Skibinski, D.; Wenne, R. Doubly Uniparental Transmission of mitochondrial DNA length variants in the mussel *Mytilus trossulus*. *Mar. Biol.* **2003**, *142*, 455–460. [CrossRef]

39. Lemer, S.; Bieler, R.; Giribet, G. Resolving the relationships of clams and cockles: Dense transcriptome sampling drastically improves the bivalve tree of life. *Proc. R. Soc. B Biol. Sci.* **2019**, *286*, 20182684. [CrossRef]

40. Rawson, P.D. Nonhomologous recombination between the large unassigned region of the male and female mitochondrial genomes in the Mussel, *Mytilus trossulus*. *J. Mol. Evol.* **2005**, *61*, 717–732. [CrossRef]
41. Cao, L.; Ort, B.S.; Mizi, A.; Pogson, G.; Kenchington, E.; Zouros, E.; Rodakis, G.C. The control region of maternally and paternally inherited mitochondrial genomes of three species of the sea mussel genus *Mytilus*. *Genet. 2009*, *181*, 1045–1056. [CrossRef] [PubMed]

42. Ladoukakis, E.D.; Theologidis, I.; Rodakis, G.C.; Zouros, E. Homologous recombination between highly diverged mitochondrial sequences: Examples from maternally and paternally transmitted genomes. *Mol. Biol. Evol.* **2011**, *28*, 1847–1859. [CrossRef]

43. Fukata, Y.; Iigo, M. F-type complete mitochondrial genomes of two *Hyriposis* species, *H. schlegelii* and *H. cumingii* (Bivalvia; Unionoida; Unionidae) from Lake Kasumigaura, Japan, and its phylogenetic analysis. *Mitochondrial DNA Part B* **2020**, *5*, 1218–1219. [CrossRef] [PubMed]

44. Ghiselli, F.; Maurizii, M.G.; Reunov, A.; Ariño-Bassols, H.; Cifaldi, C.; Pecci, A.; Alexandrova, Y.; Bettini, S.; Passamonti, M.; Franceschini, V.; et al. Natural heteroplasmy and mitochondrial inheritance in bivalve molluscs. *Integr. Comp. Biol.* **2019**, *59*, 1016–1032. [CrossRef] [PubMed]

45. Lubośny, M.; Przyłucka, A.; Śmietanka, B.; Burzyński, A. *Semimytilus algosus*: First known hermaphroditic mussel with Doubly Uniparental Inheritance of mitochondrial DNA. *Sci. Rep.* **2020**, *10*, 11256. [CrossRef]

46. Plazzi, F.; Passamonti, M. Footprints of unconventional mitochondrial inheritance in bivalve phylogeny: Signatures of positive selection on clades with doubly uniparental inheritance. *J. Zool. Syst. Evol. Res.* **2019**, *57*, 258–271. [CrossRef]

47. Zouros, E. Doubly uniparental inheritance of mitochondrial DNA: Might it be simpler than we thought? *J. Zool. Syst. Evol. Res.* **2020**, *58*, 624–631. [CrossRef]

48. Kyriakou, E.; Kravariti, L.; Zouros, E.; Rodakis, G.C. No sex-specific protein-binding site in the VD1 of the F mitochondrial genome of the mussel *Mytilus galloprovincialis*. *Gene Rep.* **2016**, *5*, 148–150. [CrossRef]

49. Pozzi, A.; Plazzi, F.; Milani, L.; Ghiselli, F.; Passamonti, M. SmithRNAs: Could mitochondria "bend" nuclear regulation? *Mol. Biol. Evol.* **2017**, *34*, 1960–1973. [CrossRef] [PubMed]