The Complete Maternal and Paternal Mitochondrial Genomes of the Mediterranean Mussel Mytilus galloprovincialis: Implications for the Doubly Uniparental Inheritance Mode of mtDNA

Mitogenome sequence comparison in the endangered congeneric Pinna nobilis and Pinna rudis bivalves

BACKGROUND

The pen shells Pinna nobilis and Pinna rudis are large wedge-shaped bivalve molluscs. Both species are threatened by different anthropogenic pressures. In the last few years, P. nobilis populations have significantly reduced due to massive mortality events. The complete mitochondrial DNA sequences of these congeneric species have been determined and compared for the first time.

RESULTS

The mitogenome sequences of P. nobilis and P. rudis were 18,919 bp and 18,264 bp in length, respectively. Each mitogenome is composed of 12 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA (tRNAs) genes and non-coding regions. A putative Adenosine Triphosphate synthase subunit 8 gene could only be proposed for P. nobilis. Both newly sequenced mitogenomes present a conserved gene order between them, comparable to the closely related Atrina pectinata, but global arrangement greatly differs from other available bivalve mitochondrial sequences. Multiple copies of tRNA-Cys were identified, located in different positions probably due to mechanisms of mitochondrial genome rearrangements, and detected 2 and 3 times in P. rudis and in P. nobilis, respectively.

CONCLUSION

A close relationship was shown between Pinna species and Atrina pectinata and a consistent clustering showing a monophyletic origin of Pinnidae family sequences was evidenced. The mitochondrial genomes will provide a valuable genetic resource for further studies on population genetics and species identification.

A mitochondrial genome phylogeny of Mytilidae (Bivalvia: Mytilida)

The family Mytilidae is a family of bivalve mussels that are distributed worldwide in diverse marine habitats. Within the family, classification systems and phylogenetic relationships among subfamilies remain not yet fully resolved. In this study, we newly determined 9 mitochondrial genome sequences from 7 subfamilies: Bathymodiolus thermophilus (Bathymodiolinae), Modiolus nipponicus (Modiolinae), Lithophaga curta (the first representative of Lithophaginae), Brachidontes mutabilis (Brachidontinae), Mytilisepta virgata (Brachidontinae), Mytilisepta keenae (Brachidontinae), Crenomytilus grayanus (Mytilinae), Gregariella coralliophaga (Crenellinae), and Septifer bilocularis (the first representative of Septiferinae). Phylogenetic trees using maximum likelihood and Bayesian inference methods for 28 mitochondrial genomes (including 19 previously published sequences) showed two major clades with high support values: Clade 1 ((Bathymodiolinae + Modiolinae) + (Lithophaginae + Limnoperninae)) and Clade 2 (((Mytilinae + Crenellinae) + Septiferinae) + Brachidontinae). The position of the genus Lithophaga (representing Lithophaginae) differed from a previously published molecular phylogeny. Divergence time analysis with a molecular clock indicated that lineage splitting among the major subfamilies of Mytilidae (including the habitat transition from marine to freshwater environments by ancestral Limnoperninae) occurred in the Mesozoic period, coinciding with high diversification rates of marine fauna during that time. This is the first mitochondrial genome-based phylogenetic study of the Mytilidae that covers nearly all subfamily members, excluding the subfamily Dacrydiinae.

Transoceanic dispersal of the mussel Mytilus galloprovincialis on Japanese tsunami marine debris: An approach for evaluating rafting of a coastal species at sea

Biofouled debris from the 2011 Great East Japan earthquake and tsunami has landed in the Northeast Pacific and along the Hawaiian Islands since 2012. As of 2017, >630 biofouled debris items with >320 living species of algae, invertebrates, and fish have been examined. The invasive mussel Mytilus galloprovincialis was present on >50% of those items. Size, reproduction, and growth of this filter-feeding species were examined to better understand long-distance rafting of a coastal species. The majority of mussels (79%) had developing or mature gametes, and growth rates averaged 0.075±0.018 SE mm/day. Structural and elemental (barium/calcium) analysis of mussel shells generated estimates of growth in coastal waters (mean=1.3 to 25mm total length), which provides an indication of residence times in waters along North America and the Hawaiian Islands prior to landing. Detailed studies of individual species contribute to our understanding of debris as a transport vector and aid efforts to evaluate potential risks associated with marine debris.

Actively transcribed and expressed atp8 gene in Mytilus edulis mussels

Background

Animal mitochondrial genomes typically encode 37 genes: 13 proteins, 22 tRNAs and two rRNAs. However, many species represent exceptions to that rule. Bivalvia along with Nematoda and Platyhelminthes are often suspected to fully or partially lack the ATP synthase subunit 8 (atp8) gene. This raises the question as to whether they are really lacking this gene or is this maybe an annotation problem? Among bivalves, Mytilus edulis has been inferred to lack an ATP8 gene since the characterization of its mitochondrial genome in 1992. Even though recent bioinformatic analyses suggested that atp8 is present in Mytilus spp., due to high divergence in predicted amino acid sequences, the existence of a functional atp8 gene in this group remains controversial.

Results

Here we demonstrate that M. edulis mitochondrial open reading frames suggested to be atp8 (in male and female mtDNAs) are actively translated proteins. We also provide evidence that both proteins are an integral part of the ATP synthase complex based on in-gel detection of ATP synthase activity and two-dimensional Blue-Native and SDS polyacrylamide electrophoresis.

Conclusion

Many organisms (e.g., Bivalvia along with Nematoda and Platyhelminthes) are considered to be lacking certain mitochondrial genes often only based on poor similarity between protein coding gene sequences in genetically closed species. In some situations, this may lead to the inference that the ATP8 gene is absent, when it is in fact present, but highly divergent. This shows how important complementary role protein-based approaches, such as those in the present study, can provide to bioinformatic, genomic studies (i.e., ability to confirm the presence of a gene).

Distribution and frequency of mitochondrial DNA polymorphisms in blue mussel (Mytilus edulis) populations of southwestern Nova Scotia (Canada)

The Atlantic blue mussel (Mytilus edulis Linnaeus, 1758) exhibits doubly uniparental inheritance of mitochondrial (mt) DNA. Females are usually homoplasmic for a female-transmitted mt genome (the F type) and males are heteroplasmic for an F type and a male-transmitted mt genome (the M type). F types can undergo “role-reversal” events, resulting in new male-transmitted mtDNA genomes known as recently masculinized (RM) types that co-occur in populations with evolutionarily older standard-male (SM) types. Phylogenetic analyses have shown that RM types periodically replace SM types. It has also been shown that sperm with RM mtDNA have greater swimming velocity and more efficient components of the electron transport chain compared to sperm with SM mtDNA, thus leading to the hypothesis that RM sperm may have a selective advantage over SM sperm. The present study examines the distribution of RM and SM mitotypes in male M. edulis (n = 225) from 13 localities in southwestern Nova Scotia (Canada). The SM type was more common in all populations, with the proportion of RM types ranging from 0% to 24.1%. The highest proportion of RM types was observed in an aquaculture operation. Analyses of additional populations are required to evaluate the selective pressures affecting the geographic distribution of RM and SM mitotypes in M. edulis.

Gender-Associated Mitochondrial DNA Heteroplasmy in Somatic Tissues of the Endangered Freshwater Mussel Unio crassus (Bivalvia: Unionidae): Implications for Sex Identification and Phylogeographical Studies

Some bivalve species possess two independent mitochondrial DNA lineages: maternally (F-type) and paternally (M-type) inherited. This phenomenon is called doubly uniparental inheritance. It is generally agreed that F-type mtDNA is typically present in female somatic and gonadal tissues as well as in male somatic tissues, whereas the M-type mtDNA occurs only in male germ line and gonadal tissue. In the present study, the mtDNA heteroplasmy (for both F and M genomes) in male somatic tissues of Unio crassus (Philipsson, 1788), species threatened with extinction, has been confirmed. Taking advantage from the presence of Mcox1 marker only in male somatic tissues, we developed a new method of sex identification in this endangered species, using nondestructive tissue sampling. Probability of correct sex identification was estimated at 97.5%. The present study is the first report on gender-associated mitochondrial DNA heteroplasmy in male somatic tissues of thick-shelled river mussel and first approach to U. crassus sex identification at molecular level. Our study also confirmed the utility of paternally inherited Mcox1 gene fragment as a complementary molecular tool for resolving phylogeographical relationships among populations of thick-shelled river mussel.

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

Mitogenomic trees for Bivalvia have proved problematic in the past, but several highly divergent lineages were missing from these analyses and increased representation of these groups may yet improve resolution. Here, we add seven new sequences from the Anomalodesmata and one unidentified semelid species (Bryopa lata, Euciroa cf. queenslandica, Laternula elliptica, Laternula truncata, Lyonsia norwegica, Myadora brevis, Tropidomya abbreviata, "Abra" sp.). We show that relationships in a mitogenomic tree for the Class are improved by the addition of seven anomalodesmatans from this highly divergent clade, but are still not completely consistent with relationships recovered in studies of nuclear genes. We suggest that some anomalous relationships (for instance the non-monophyly of Bivalvia) may be partially explained by compositional heterogeneity in the mitogenome and suggest that the addition of more taxa may help resolve both this effect and possible instances of long branch attraction. We also identify several curious features about anomalodesmatan mitogenomes. For example, many protein-coding gene boundaries are poorly defined in marine bivalves, but particularly so in anomalodesmatans, primarily due to non-conserved boundary sequences. The use of transcriptomic and genomic data together enabled better definition of gene boundaries, the identification of possible pseudogenes and suggests that most genes are translated monocistronically, which contrasts with many other studies. We also identified a possible case of gene duplication of ND5 in Myadora brevis (Myochamidae). Mitogenome size in the Anomalodesmata ranges from very small compact molecules, with the smallest for Laternula elliptica (Laternulidae) only 14,622bp, to Bryopa lata (Clavagellidae) which is at least 31,969bp long and may be >40,000bp. Finally, sampled species show a high degree of sequence divergence and variable gene order, although intraspecific variation in Laternula elliptica is very low.

Complete female mitochondrial genome of Mytilus chilensis

The controversy surrounding the origin of antitropical distribution of Mytilus mussels and the taxonomic status of southern hemisphere populations remain unsolved, despite the efforts. One of the limiting factors remains the lack of the complete sequences of the representative mitochondrial genomes which would allow their proper comparison with the relatively well-represented northern hemisphere congeneric mussels. To fill this gap we sequenced the representative maternal (F) genome of a native Chilean mussel. The genome is 16,748bp long and structurally identical to the northern hemisphere M. edulis and M. galloprovincialis F genomes. However, the genetic distance from them (≈5%) is twice as high as the maximum distance between them (<2.5%). Thus, the notion that the name M. chilensis should be used for native Chilean Mytilus mussels, with the same rank as M. galloprovincialis and M. edulis is supported.

Sequence motifs associated with paternal transmission of mitochondrial DNA in the horse mussel, Modiolus modiolus (Bivalvia: Mytilidae)

In the majority of metazoans paternal mitochondria represent evolutionary dead-ends. In many bivalves, however, this paradigm does not hold true; both maternal and paternal mitochondria are inherited. Herein, we characterize maternal and paternal mitochondrial control regions of the horse mussel, Modiolus modiolus (Bivalvia: Mytilidae). The maternal control region is 808bp long, while the paternal control region is longer at 2.3kb. We hypothesize that the size difference is due to a combination of repeated duplications within the control region of the paternal mtDNA genome, as well as an evolutionarily ancient recombination event between two sex-associated mtDNA genomes that led to the insertion of a second control region sequence in the genome that is now transmitted via males. In a comparison to other mytilid male control regions, we identified two evolutionarily Conserved Motifs, CMA and CMB, associated with paternal transmission of mitochondrial DNA. CMA is characterized by a conserved purine/pyrimidine pattern, while CMB exhibits a specific 13bp nucleotide string within a stem and loop structure. The identification of motifs CMA and CMB in M. modiolus extends our understanding of Sperm Transmission Elements (STEs) that have recently been identified as being associated with the paternal transmission of mitochondria in marine bivalves.

Comparative population structure of two dominant species, Shinkaia crosnieri (Munidopsidae: Shinkaia) and Bathymodiolus platifrons (Mytilidae: Bathymodiolus), inhabiting both deep-sea vent and cold seep inferred from mitochondrial multi-genes

Deep-sea hydrothermal vents and cold seeps, limited environments without sunlight, are two types of extreme habitat for marine organisms. The differences between vents and cold seeps may facilitate genetic isolation and produce population heterogeneity. However, information on such chemosynthetic fauna taxa is rare, especially regarding the population diversity of species inhabiting both vents and cold seeps. In this study, three mitochondrial DNA fragments (the cytochrome c oxidase submit I (COI), cytochrome b gene (Cytb), and 16S) were concatenated as a mitochondrial concatenated dataset (MCD) to examine the genetic diversity, population structure, and demographic history of Shinkaia crosnieri and Bathymodiolus platifrons. The genetic diversity differences between vent and seep populations were statistically significant for S. crosnieri but not for B. platifrons. S. crosnieri showed less gene flow and higher levels of genetic differentiation between the vent and seep populations than B. platifrons. In addition, the results suggest that all the B. platifrons populations, but only the S. crosnieri vent populations, passed through a recent expansion or bottleneck. Therefore, different population distribution patterns for the two dominant species were detected; a pattern of population differentiation for S. crosnieri and a homogeneity pattern for B. platifrons. These different population distribution patterns were related to both extrinsic restrictive factors and intrinsic factors. Based on the fact that the two species were collected in almost identical or adjacent sampling sites, we speculated that the primary factors underlying the differences in the population distribution patterns were intrinsic. The historical demographics, dispersal ability, and the tolerance level of environmental heterogeneity are most likely responsible for the different distribution patterns.

Mitogenomics of southern hemisphere blue mussels (Bivalvia: Pteriomorphia): Insights into the evolutionary characteristics of the Mytilus edulis complex

Marine blue mussels (Mytilus spp.) are widespread species that exhibit an antitropical distribution with five species occurring in the Northern Hemisphere (M. trossulus, M. edulis, M. galloprovincialis, M. californianus and M. coruscus) and three in the Southern Hemisphere (M. galloprovincialis, M. chilensis and M. platensis). Species limits in this group remain controversial, in particular for those forms that live in South America. Here we investigated structural characteristics of marine mussels mitogenomes, based on published F mtDNA sequences of Northern Hemisphere species and two newly sequenced South American genomes, one from the Atlantic M. platensis and another from the Pacific M. chilensis. These mitogenomes exhibited similar architecture to those of other genomes of Mytilus, including the presence of the Atp8 gene, which is missing in most of the other bivalves. Our evolutionary analysis of mitochondrial genes indicates that purifying selection is the predominant force shaping the evolution of the coding genes. Results of our phylogenetic analyses supported the monophyly of Pteriomorphia and fully resolved the phylogenetic relationships among its five orders. Finally, the low genetic divergence of specimens assigned to M. chilensis and M. platensis suggests that these South American marine mussels represent conspecific variants rather than distinct species.

A First Insight into the Genome of the Filter-Feeder Mussel Mytilus galloprovincialis

Mussels belong to the phylum Mollusca, one of the largest and most diverse taxa in the animal kingdom. Despite their importance in aquaculture and in biology in general, genomic resources from mussels are still scarce. To broaden and increase the genomic knowledge in this family, we carried out a whole-genome sequencing study of the cosmopolitan Mediterranean mussel (Mytilus galloprovincialis). We sequenced its genome (32X depth of coverage) on the Illumina platform using three pair-end libraries with different insert sizes. The large number of contigs obtained pointed out a highly complex genome of 1.6 Gb where repeated elements seem to be widespread (~30% of the genome), a feature that is also shared with other marine molluscs. Notwithstanding the limitations of our genome sequencing, we were able to reconstruct two mitochondrial genomes and predict 10,891 putative genes. A comparative analysis with other molluscs revealed a gene enrichment of gene ontology categories related to multixenobiotic resistance, glutamate biosynthetic process, and the maintenance of ciliary structures.

The complete mitochondrial genome of the golden mussel Limnoperna fortunei and comparative mitogenomics of Mytilidae

Here we describe the mitochondrial genome of the golden mussel Limnoperna fortunei, an Asian bivalve which has become one of the most aggressive invasive species in Japan and South America. The mitochondrial genome of L. fortunei does not present conserved gene arrangement when compared to the other Mytilidae species suggesting a high degree of gene recombination in the mitochondria of this clade. In addition, the golden mussel mitogenome encodes two copies of tRNA(Lys) and presents a putative pseudogene for the atp8 gene sequence that encodes a 27 amino acid peptide containing an in-frame stop codon. The presence of this pseudogene raises the question as to whether atp8 is encoded in some bivalve mitochondrial genomes or not. The phylogenetic analysis of all complete mitochondrial genomes available from Mytilidae mussels confirmed the close evolutionary relationships among bivalves from the genus Mytilys and placed L. fortunei coming from a more ancestral branch on the family. The supermatrix phylogeny described used the concatenation of all 12 genes from the mitochondria and disputed the monophyly of the genus Perna, as Perna perna was shown to be more closely related to Brachidontes exustus than to Perna viridis. The comparative analysis of mitogenome synteny also confirmed the polyphyly of the genus Perna. The complete and annotated mitogenome has been published in GenBank under the accession number KP756905.

A Distinct Mitochondrial Genome with DUI-Like Inheritance in the Ocean Quahog Arctica islandica

Mitochondrial DNA (mtDNA) is strictly maternally inherited in metazoans. The major exception to this rule has been found in many bivalve species which allow the presence of different sex-linked mtDNA molecules. This mechanism, named doubly uniparental inheritance (DUI), is characterized by the presence of two mtDNAs: The female mtDNA is found in somatic tissue and female gonads, whereas the male mtDNA is usually found in male gonads and sperm. In this study we highlight the existence of two divergent mitochondrial haplotypes with a low genetic difference around 6–8% in Arctica islandica, a long-lived clam belonging to the Arcticidae, a sister group to the Veneridae in which DUI has been found. Phylogenetic analysis on cytochrome b and 16S sequences from somatic and gonadic tissues of clams belonging to different populations reveals the presence of the “divergent” type in male gonads only and the “normal” type in somatic tissues and female gonads. This peculiar segregation of divergent mtDNA types speaks for the occurrence of the DUI mechanism in A. islandica. This example also highlights the difficulties to assess the presence of such particular mitochondrial inheritance system and underlines the possible misinterpretations in phylogeographic and phylogenetic studies of bivalve species linked to the presence of two poorly differentiated mitochondrial genomes.

Presence of two mitochondrial genomes in the mytilid Perumytilus purpuratus: Phylogenetic evidence for doubly uniparental inheritance

This study presents evidence, using sequences of ribosomal 16S and COI mtDNA, for the presence of two mitochondrial genomes in Perumytilus purpuratus. This may be considered evidence of doubly uniparental mtDNA inheritance. The presence of the two types of mitochondrial genomes differentiates females from males. The F genome was found in the somatic and gonadal tissues of females and in the somatic tissues of males; the M genome was found in the gonads and mantle of males only. For the mitochondrial 16S region, ten haplotypes were found for the F genome (nucleotide diversity 0.004), and 7 haplotypes for the M genome (nucleotide diversity 0.001), with a distance Dxy of 0.125 and divergence Kxy of 60.33%. For the COI gene 17 haplotypes were found for the F genome (nucleotide diversity 0.009), and 10 haplotypes for the M genome (nucleotide diversity 0.010), with a genetic distance Dxy of 0.184 and divergence Kxy of 99.97%. Our results report the presence of two well-differentiated, sex-specific types of mitochondrial genome (one present in the male gonad, the other in the female gonad), implying the presence of DUI in P. purpuratus. These results indicate that care must be taken in phylogenetic comparisons using mtDNA sequences of P. purpuratus without considering the sex of the individuals.

A protein binding site in the M mitochondrial genome of Mytilus galloprovincialis may be responsible for its paternal transmission

Sea mussels (genus Mytilus) have two mitochondrial genomes in obligatory co-existence, one that is transmitted through the egg and the other through the sperm. The phenomenon, known as Doubly Uniparental Inheritance (DUI) of mitochondrial DNA (mtDNA), is presently known to occur in more than 40 molluscan bivalve species. Females and the somatic tissues of males contain mainly the maternal (F) genome. In contrast, the sperm contains only the paternal (M) genome. Through electrophoretic mobility shift assay (EMSA) experiments we have identified a sequence element in the control region (CR) of the M genome that acts as a binding site for the formation of a complex with a protein factor that occurs in the male gonad. An adenine tract upstream to the element is also essential for the formation of the complex. The reaction is highly specific. It does not occur with protein extracts from the female gonad or from a male or female somatic tissue. Further experiments showed that the interaction takes place in mitochondria surrounding the nucleus of the cells of male gonads, suggesting a distinct role of perinuclear mitochondria. We propose that at a certain point during spermatogenesis mitochondria are subject to degradation and that perinuclear mitochondria with the M mtDNA-protein complex are protected from this degradation with the result that mature spermatozoa contain only the paternal mitochondrial genome.

Mitogenomics of recombinant mitochondrial genomes of Baltic Sea Mytilus mussels

Recombination in the control region (CR) of Mytilus mitochondrial DNA (mtDNA) was originally reported based on the relatively short, sequenced fragments of mitochondrial genomes. Recombination outside the CR has been reported recently with the suggestion that such processes are common in Mytilus. We have fully sequenced a set of 11 different mitochondrial haplotypes representing the high diversity of paternally inherited mitochondrial genomes of Baltic Sea Mytilus mussels, including the haplotype close to the native Mytilus trossulus mitochondrial genome, which was thought to have been entirely eliminated from this population. Phylogenetic and comparative analysis showed that the recombination is limited to the vicinity of the CR in all sequenced genomes. Coding sequence comparison indicated that all paternally inherited genomes showed increased accumulation of nonsynonymous substitutions, including the genomes which switched their transmission route very recently. The acquisition of certain CR sequences through recombination with highly divergent paternally inherited genomes seems to precede and favor the switch, but it is not a prerequisite for this process. Interspecies hybridization in the Baltic Sea during the recent 10,000 years created conditions for both structural and evolutionary mitochondrial instability which resulted in the observed variation and dynamics of mtDNA in Baltic Sea Mytilus mussels. In conclusion, the data shows that the effects of mitochondrial recombination are limited to the CR of few phylogenetic lineages.

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

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

Complete mitochondrial genome of the giant African snail, Achatina fulica (Mollusca: Achatinidae): a novel location of putative control regions (CR) in the mitogenome within Pulmonate species

The whole sequence (15,057 bp) of the mitochondrial DNA (mtDNA) of the terrestrial snail Achatina fulica (order Stylommatophora) was determined. The mitogenome, as the typical metazoan mtDNA, contains 13 protein-coding genes (PCG), 2 ribosomal RNA genes (rRNA) and 22 transfer RNA genes (tRNA). The tRNA genes include two trnS without standard secondary structure. Interestingly, among the known mitogenomes of Pulmonata species, we firstly characterized an unassigned lengthy sequence (551 bp) between the cox1 and the trnV which may be the CR for the sake of its AT bases usage bias (65.70%) and potential hairpin structure.

Mitochondrial control region variability in Mytilus galloprovincialis populations from the central-Eastern Mediterranean Sea

The variable domain 1 (VD1) domain of the control region and a small segment of the rrnaL gene of the F mtDNA type were sequenced and analyzed in 174 specimens of Mytilus galloprovincialis. Samples were collected from eight locations in four Central-Eastern (CE) Mediterranean countries (Italy, Croatia, Greece and Turkey). A new primer, specific for the F mtDNA type, was designed for the sequencing procedure. In total 40 different haplotypes were recorded, 24 of which were unique. Aside from the two populations situated in Thermaikos gulf (Northern Aegean, Greece), relatively high levels of haplotype and nucleotide diversity were estimated for both Central and Eastern Mediterranean populations. Eight out of the 40 haplotypes were shared by at least three populations while two of them were found in all populations. ΦST and cluster analysis revealed lack of structuring among CE Mediterranean populations with the exception of those located at the Sea of Marmara and Croatian coast which were highly differentiated. Apart from the species' inherit dispersal ability, anthropogenic activities, such as the repeated translocations of mussel spat, seem to have played an important role in shaping the current genetic population structure of CE M. galloprovincialis mussels.

The F type mitochondrial genome of hard-shelled mussel: Mytilus coruscus (Mytiloida, Mytilidae)

We determined the complete mitochondrial genome of the female hard-shelled mussel Mytilus coruscus (Gould, 1869) (F type). The F type genome is composed of 16,642 bp in length including 12 protein coding genes (PCGs), 2 rRNA and 23 tRNA with the same gene content and order as the other Mytilus species which is characterized by the absence of the ATPase8 gene and addition of tRNA(Met) (AUA). The nucleotide composition of the genome shows that the percentage of A+T (63.2%) is higher than those of M. edulis complex species (M. edulis, M. galloprovincialis, M. trossulus) but lower than that of M. californianus. The F type mitochondrial genome of M. coruscus will provide useful information on the evolutionary history of the Mytilus species in the family Mitilidae.

Co-expressed mitochondrial genomes: recently masculinized, recombinant mitochondrial genome is co-expressed with the female-transmitted mtDNA genome in a male Mytilus trossulus mussel from the Baltic Sea

BACKGROUND

Few exceptions have been described from strict maternal inheritance of mitochondrial DNA in animals, including sea mussels (Mytilidae), clams (Donacidae, Veneridae and Solenidae) and freshwater mussels (Unionoidae) order. In these bivalves mitochondria and their DNA are transferred through two separate routes. The females inherit only the maternal mitochondrial DNA whereas the males inherit maternal as well as paternal mitochondrial DNA, which is usually present only in gonads and sperm. The mechanism controlling this phenomenon is unclear but leads to the existence of two separate mitochondrial DNA lineages in a single species. The lineages are usually well differentiated: up to 20-50% divergence in nucleotide sequence. Occasionally, a maternal mitochondrial DNA can invade the paternal transmission route, eventually replacing the diverged M-type and lowering the divergence. Such role reversal (masculinization) event has happened recently in the Mytilus population of the Baltic Sea which consists of M. edulis × M. trossulus hybrids, but the functional status of the resulting mitochondrial genome was unknown.

RESULTS

In this paper we sequenced transcripts from one specimen that was identified as male carrying both the female mitochondrial genome and a recently masculinized mitochondrial genome. Additionally, the analysis of the control region has showed that the recently masculinized, recombinant genome, not only has an M-type control region and all coding regions derived from the F-type, but also is transcriptionally active along side the maternally inherited F-type genome. In the comparative analysis, the two genomes exhibit different substitution patterns, typical for the M vs. F genome comparisons. The genetic distances and ratios of non-synonymous substitutions also suggest that one of the genomes is transitioning from the maternal to the paternal inheritance mode, consistent with its recent masculinization.

CONCLUSION

We have shown, for the first time, that the recently masculinized mitochondrial genome is active and that it accumulates excess of non-synonymous substitutions across its coding sequence. This suggests, that, under certain cytonuclear incompatibility conditions, masculinization may serve to restore the endangered functionality of the paternally inherited genome. This is also another example of a mitochondrial genome in which the recombination in the control region predated its transition from paternal to maternal transmission route.

The rRNA and tRNA transcripts of maternally and paternally inherited mitochondrial DNAs of Mytilus galloprovincialis suggest presence of a "degradosome" in mussel mitochondria and necessitate the re-annotation of the l-rRNA/CR boundary

Species of the genus Mytilus carry two mitochondrial genomes in obligatory coexistence; one transmitted though the eggs (the F type) and one through the sperm (the M type). We have studied the 3' and 5' ends of rRNA and tRNA transcripts using RT-PCR and RNA circularization techniques in both the F and M genomes of Mytilus galloprovincialis. We have found polyadenylated and non-adenylated transcripts for both ribosomal and transfer RNAs. In all these genes the 5' ends of the transcripts coincided with the first nucleotide of the annotated genes, but the 3' ends were heterogeneous. The l-rRNA 3' end is 47 or 48 nucleotides upstream from the one assigned by a previous annotation, which makes the adjacent first domain (variable domain one, VD1) of the main control region (CR) correspondingly longer. We have observed s-rRNA and l-rRNA transcripts with truncated 3' end and polyadenylated tRNA transcripts carrying the CCA trinucleotide. We have also detected polyadenylated RNA remnants carrying the sequences of the control region, which strongly suggests RNA degradation activity and thus presence of degradosomes in Mytilus mitochondria.

The complete maternally and paternally inherited mitochondrial genomes of the endangered freshwater mussel Solenaia carinatus (Bivalvia: Unionidae) and implications for Unionidae taxonomy

Doubly uniparental inheritance (DUI) is an exception to the typical maternal inheritance of mitochondrial (mt) DNA in Metazoa, and found only in some bivalves. In species with DUI, there are two highly divergent gender-associated mt genomes: maternal (F) and paternal (M), which transmit independently and show different tissue localization. Solenaia carinatus is an endangered freshwater mussel species exclusive to Poyang Lake basin, China. Anthropogenic events in the watershed greatly threaten the survival of this species. Nevertheless, the taxonomy of S. carinatus based on shell morphology is confusing, and the subfamilial placement of the genus Solenaia remains unclear. In order to clarify the taxonomic status and discuss the phylogenetic implications of family Unionidae, the entire F and M mt genomes of S. carinatus were sequenced and compared with the mt genomes of diverse freshwater mussel species. The complete F and M mt genomes of S. carinatus are 16716 bp and 17102 bp in size, respectively. The F and M mt genomes of S. carinatus diverge by about 40% in nucleotide sequence and 48% in amino acid sequence. Compared to F counterparts, the M genome shows a more compact structure. Different gene arrangements are found in these two gender-associated mt genomes. Among these, the F genome cox2-rrnS gene order is considered to be a genome-level synapomorphy for female lineage of the subfamily Gonideinae. From maternal and paternal mtDNA perspectives, the phylogenetic analyses of Unionoida indicate that S. carinatus belongs to Gonideinae. The F and M clades in freshwater mussels are reciprocal monophyly. The phylogenetic trees advocate the classification of sampled Unionidae species into four subfamilies: Gonideinae, Ambleminae, Anodontinae, and Unioninae, which is supported by the morphological characteristics of glochidia.

The largest unassigned regions of the male- and female-transmitted mitochondrial DNAs in Musculista senhousia (Bivalvia Mytilidae)

Musculista senhousia is a marine mussel with doubly uniparental inheritance (DUI) of mitochondria. In this study we analyzed the largest unassigned region (LUR) of its female- and male-transmitted mitochondrial genomes, described their fine characteristics and searched for shared features. Our results suggest that both LURs contain the control region of their respective mitochondrial genomes. The female-transmitted control region is duplicated in tandem, with the two copies evolving in concert. This makes the F-mtDNA of M. senhousia the first Bivalve mitochondrial genome with this feature. We also compared M. senhousia control regions to that of other Mytilidae, and demonstrated that signals for basic mtDNA functions are retained over evolutionary times even among the fast-evolving mitochondrial genomes of DUI species. Finally, we discussed how similarities between female and male LURs may be explained in the context of DUI evolution and if the duplicated female control region might have influenced the DUI system in this species.

Complete female mitochondrial genome of Anodonta anatina (Mollusca: Unionidae): confirmation of a novel protein-coding gene (F ORF)

Freshwater mussels are among animals having two different, gender-specific mitochondrial genomes. We sequenced complete female mitochondrial genomes from five individuals of Anodonta anatina, a bivalve species common in palearctic ecozone. The length of the genome was variable: 15,637-15,653 bp. This variation was almost entirely confined to the non-coding parts, which constituted approximately 5% of the genome. Nucleotide diversity was moderate, at 0.3%. Nucleotide composition was typically biased towards AT (66.0%). All genes normally seen in animal mtDNA were identified, as well as the ORF characteristic for unionid mitochondrial genomes, bringing the total number of genes present to 38. If this additional ORF does encode a protein, it must evolve under a very relaxed selection since all substitutions within this gene were non-synonymous. The gene order and structure of the genome were identical to those of all female mitochondrial genomes described in unionid bivalves except the Gonideini.

Diversity of mitochondrial genome organization

In this review, we discuss types of mitochondrial genome structural organization (architecture), which includes the following characteristic features: size and the shape of DNA molecule, number of encoded genes, presence of cryptogenes, and editing of primary transcripts.

Structure, transcription, and variability of metazoan mitochondrial genome: perspectives from an unusual mitochondrial inheritance system

Despite its functional conservation, the mitochondrial genome (mtDNA) presents strikingly different features among eukaryotes, such as size, rearrangements, and amount of intergenic regions. Nonadaptive processes such as random genetic drift and mutation rate play a fundamental role in shaping mtDNA: the mitochondrial bottleneck and the number of germ line replications are critical factors, and different patterns of germ line differentiation could be responsible for the mtDNA diversity observed in eukaryotes. Among metazoan, bivalve mollusc mtDNAs show unusual features, like hypervariable gene arrangements, high mutation rates, large amount of intergenic regions, and, in some species, an unique inheritance system, the doubly uniparental inheritance (DUI). The DUI system offers the possibility to study the evolutionary dynamics of mtDNAs that, despite being in the same organism, experience different genetic drift and selective pressures. We used the DUI species Ruditapes philippinarum to study intergenic mtDNA functions, mitochondrial transcription, and polymorphism in gonads. We observed: 1) the presence of conserved functional elements and novel open reading frames (ORFs) that could explain the evolutionary persistence of intergenic regions and may be involved in DUI-specific features; 2) that mtDNA transcription is lineage-specific and independent from the nuclear background; and 3) that male-transmitted and female-transmitted mtDNAs have a similar amount of polymorphism but of different kinds, due to different population size and selection efficiency. Our results are consistent with the hypotheses that mtDNA evolution is strongly dependent on the dynamics of germ line formation, and that the establishment of a male-transmitted mtDNA lineage can increase male fitness through selection on sperm function.

[Phylogenetic application and analysis of horizontal transfer based on the prokaryote eno gene]

The phenomenon of conflicting gene trees has become a remarkable and difficult problem. Application of multiple genes has been a widespread practice to reconstruct phylogenies in phylogenetic studies. Enolase is a key glycolytic enzyme, The enzymes from a large variety of organisms, including archaebacteria, eubacteria and eukaryotes, were studied. We downloaded eno sequences from the genomes of bacteria and archaea that have been completely sequenced. The comprehensive homology search and phylogenetic analysis of the eno were used, and nineteen horizontally transferred genes were identified. The results of analysis showed lots of differences between the features of horizontal transferred genes and the ones of whole genomic genes, such as nucleotide composition, gene combination, codon usage bias, and selection pressure. These results reconfirmed that the horizontally transferred genes were exogenous. The result revealed that prokaryote eno genes were highly conserved, medium-sized, is a good material in the phylogenetic. This paper can provide a reference in study of life habit and evolutionary history of donor and receptor, and enolase structure and function.

Mitogenomics does not resolve deep molluscan relationships (yet?)

The origin of molluscs among lophotrochozoan metazoans is unresolved and interclass relationships are contradictory between morphology-based, multi-locus, and recent phylogenomic analyses. Within the "Deep Metazoan Phylogeny" framework, all available molluscan mitochondrial genomes were compiled, covering 6 of 8 classes. Genomes were reannotated, and 13 protein coding genes (PCGs) were analyzed in various taxon settings, under multiple masking and coding regimes. Maximum Likelihood based methods were used for phylogenetic reconstructions. In all cases, molluscs result mixed up with lophotrochozoan outgroups, and most molluscan classes with more than single representatives available are non-monophyletic. We discuss systematic errors such as long branch attraction to cause aberrant, basal positions of fast evolving ingroups such as scaphopods, patellogastropods and, in particular, the gastropod subgroup Heterobranchia. Mitochondrial sequences analyzed either as amino acids or nucleotides may perform well in some (Cephalopoda) but not in other palaeozoic molluscan groups; they are not suitable to reconstruct deep (Cambrian) molluscan evolution. Supposedly "rare" mitochondrial genome level features have long been promoted as phylogenetically informative. In our newly annotated data set, features such as genome size, transcription on one or both strands, and certain coupled pairs of PCGs show a homoplastic, but obviously non-random distribution. Apparently congruent (but not unambiguous) signal for non-trivial subclades, e.g. for a clade composed of pteriomorph and heterodont bivalves, needs confirmation from a more comprehensive bivalve sampling. We found that larger clusters not only of PCGs but also of rRNAs and even tRNAs can bear local phylogenetic signal; adding trnG-trnE to the end of the ancestral cluster trnM-trnC-trnY-trnW-trnQ might be synapomorphic for Mollusca. Mitochondrial gene arrangement and other genome level features explored and reviewed herein thus failed as golden bullets, but are promising as additional characters or evidence supporting deep molluscan clades revealed by other data sets. A representative and dense sampling of molluscan subgroups may contribute to resolve contentious interclass relationships in the future, and is vital for exploring the evolution of especially diverse mitochondrial genomes in molluscs.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.