Complete sequences of maternally inherited mitochondrial genomes in mussels Unio pictorum (Bivalvia, Unionidae)

Investigating the Impact of a Curse: Diseases, Population Isolation, Evolution and the Mother's Curse

The mitochondrion was characterized for years as the energy factory of the cell, but now its role in many more cellular processes is recognized. The mitochondrion and mitochondrial DNA (mtDNA) also possess a set of distinct properties, including maternal inheritance, that creates the Mother's Curse phenomenon. As mtDNA is inherited from females to all offspring, mutations that are harmful to males tend to accumulate more easily. The Mother's Curse is associated with various diseases, and has a significant effect on males, in many cases even affecting their reproductive ability. Sometimes, it even leads to reproductive isolation, as in crosses between different populations, the mitochondrial genome cannot cooperate effectively with the nuclear one resulting in a mito-nuclear incompatibility and reduce the fitness of the hybrids. This phenomenon is observed both in the laboratory and in natural populations, and have the potential to influence their evolution and speciation. Therefore, it turns out that the study of mitochondria is an exciting field that finds many applications, including pest control, and it can shed light on the molecular mechanism of several diseases, improving successful diagnosis and therapeutics. Finally, mito-nuclear co-adaptation, paternal leakage, and kin selection are some mechanisms that can mitigate the impact of the Mother's Curse.

No evidence of DUI in the Mediterranean alien species Brachidontes pharaonis (P. Fisher, 1870) despite mitochondrial heteroplasmy

Two genetically different mitochondrial haplogroups of Brachidontes pharaonis (p-distance 6.8%) have been identified in the Mediterranean Sea. This hinted at a possible presence of doubly uniparental inheritance in this species. To ascertain this possibility, we sequenced two complete mitogenomes of Brachidontes pharaonis mussels and performed a qPCR analysis to measure the relative mitogenome copy numbers of both mtDNAs. Despite the presence of two very similar regions composed entirely of repetitive sequences in the two haplogroups, no recombination between mitogenomes was detected. In heteroplasmic individuals, both mitogenomes were present in the generative tissues of both sexes, which argues against the presence of doubly uniparental inheritance in this species.

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher d_N /d_S ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

Complete paternally inherited mitogenomes of two freshwater mussels Unio pictorum and Sinanodonta woodiana (Bivalvia: Unionidae)

Freshwater bivalves from the family Unionidae usually have two very divergent mitogenomes, inherited according to the doubly uniparental model. The early divergence of these two mitogenomic lineages gives a unique opportunity to use two mitogenomic data sets in a single phylogenetic context. However, the number of complete sequences of the maternally inherited mitogenomes of these animals available in GenBank greatly exceeds that of the paternally inherited mitogenomes. This is a problem for phylogenetic reconstruction because it limits the use of both mitogenomic data sets. Moreover, since long branch attraction phenomenon can bias reconstructions if only a few but highly divergent taxa are considered, the shortage of the faster evolving paternally inherited mitogenome sequences is a real problem. Here we provide, for the first time, complete sequences of the M mitogenomes sampled from Polish populations of two species: native Unio pictorum and invasive Sinanodonta woodiana. It increases the available set of mitogenomic pairs to 18 species per family, and allows unambiguous reconstruction of phylogenetic relationships among them. The reconstructions based on M and F mitogenomes which were separated for many millions of years, and subject to differing evolutionary dynamics, are fully congruent.

Testing the utility of DNA barcodes and a preliminary phylogenetic framework for Chinese freshwater mussels (Bivalvia: Unionidae) from the middle and lower Yangtze River

The middle and lower portions of the Yangtze River basin is the most species-rich region for freshwater mussels in Asia. The management and conservation of the taxa in this region has been greatly hampered by the lack of a well-developed phylogeny and species-level taxonomic framework. In this study, we tested the utility of two mitochondrial genes commonly used as DNA barcodes: the first subunit of the cytochrome oxidase c gene (COI) and the first subunit of the NADH dehydrogenase gene (ND1) for 34 putative species representing 15 genera, and also generated phylogenetic hypotheses for Chinese unionids based on the combined dataset of the two mitochondrial genes. The results showed that both loci performed well as barcodes for species identification, but the ND1 sequences provided better resolution when compared to COI. Based on the two-locus dataset, Bayesian Inference (BI) and Maximum Likelihood (ML) phylogenetic analyses indicated 3 of the 15 genera of Chinese freshwater mussels examined were polyphyletic. Additionally, the analyses placed the 15 genera into 3 subfamilies: Unioninae (Aculamprotula, Cuneopsis, Nodularia and Schistodesmus), Gonideninae (Lamprotula, Solenaia and Ptychorhychus) and Anodontinae (Cristaria, Arconaia, Acuticosta, Lanceolaria, Anemina and Sinoanodonta). Our results contradict previous taxonomic classification that placed the genera Arconaia, Acuticosta and Lanceolaria in the Unioninae. This study represents one of the first attempts to develop a molecular phylogenetic framework for the Chinese members of the Unionidae and will provide a basis for future research on the evolution, ecology, and conservation of Chinese freshwater mussels.

Hermaphroditic freshwater mussel Anodonta cygnea does not have supranumerary open reading frames in the mitogenome

The complete mitogenome of Anodonta cygnea is 15,613 bp long. This compact, circular molecule contains the set of 37 genes, typical for invertebrate mitogenomes, in the same order and orientation as in maternally inherited genomes of other bivalves from the same subfamily. There are only two unassigned regions longer than 200 bp (266 bp and 274 bp) and no indication of any supranumerary open reading frames.

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.

Evolution of sex-dependent mtDNA transmission in freshwater mussels (Bivalvia: Unionida)

Doubly uniparental inheritance (DUI) describes a mode of mtDNA transmission widespread in gonochoric freshwater mussels (Bivalvia: Palaeoheterodonta: Unionida). In this system, both female- and male-transmitted mtDNAs, named F and M respectively, coexist in the same species. In unionids, DUI is strictly correlated to gonochorism and to the presence of the atypical open reading frames (ORFans) F-orf and M-orf, respectively inside F and M mtDNAs, which are hypothesized to participate in sex determination. However, DUI is not found in all three Unionida superfamilies (confirmed in Hyrioidea and Unionoidea but not in Etherioidea), raising the question of its origin in these bivalves. To reconstruct the co-evolution of DUI and of ORFans, we sequenced the mtDNAs of four unionids (two gonochoric with DUI, one gonochoric and one hermaphroditic without DUI) and of the related gonochoric species Neotrigonia margaritacea (Palaeoheterodonta: Trigoniida). Our analyses suggest that rearranged mtDNAs appeared early during unionid radiation, and that a duplicated and diverged atp8 gene evolved into the M-orf associated with the paternal transmission route in Hyrioidea and Unionoidea, but not in Etherioidea. We propose that novel mtDNA-encoded genes can deeply influence bivalve sex determining systems and the evolution of the mitogenomes in which they occur.

The Complete Maternally and Paternally Inherited Mitochondrial Genomes of a Freshwater Mussel Potamilus alatus (Bivalvia: Unionidae)

Doubly uniparental inheritance (DUI) of mitochondrial DNA, found only in some bivalve families and characterized by the existence of gender-associated mtDNA lineages that are inherited through males (M-type) or females (F-type), is one of the very few exceptions to the general rule of strict maternal mtDNA inheritance in animals. M-type sequences are often undetected and hence still underrepresented in the GenBank, which hinders the progress of the understanding of the DUI phenomenon. We have sequenced and analyzed the complete M and F mitogenomes of a freshwater mussel, Potamilus alatus. The M-type was 493 bp longer (M = 16 560, F = 16 067 bp). Gene contents, order and the distribution of genes between L and H strands were typical for unionid mussels. Candidates for the two ORFan genes (forf and morf) were found in respective mitogenomes. Both mitogenomes had a very similar A+T bias: F = 61% and M = 62.2%. The M mitogenome-specific cox2 extension (144 bp) is much shorter than in other sequenced unionid mitogenomes (531-576 bp), which might be characteristic for the Potamilus genus. The overall topology of the phylogenetic tree is in very good agreement with the currently accepted phylogenetic relationships within the Unionidae: both studied sequences were placed within the Ambleminae subfamily clusters in the corresponding M and F clades.

The complete maternal and paternal mitochondrial genomes of Unio crassus: Mitochondrial molecular clock and the overconfidence of molecular dating

The availability of a rapidly growing number of complete mitochondrial genome sequences provokes high confidence dating approaches. However, even if the congruence between mitochondrial and nuclear markers is reasonable, the resulting topologies are frequently questionable. The unique opportunity to study the evolutionary history of two independent mitochondrial genomes in one phylogenetic context exists in the freshwater mussels family Unionidae. The two lineages function under doubly uniparental inheritance since well before the emergence of the family. Despite the relatively high number of available complete sequences of maternally inherited genomes, comparative analyses are limited by the small number of sequences of counterpart paternally inherited genomes. We have sequenced for the first time the representative set of five sequences (two maternal and three paternal) from the species Unio crassus. Comparative analysis of the phylogenies reconstructed using relevant mitogenomic data available in GenBank (13 species in total) reveal that single - genome inferences are congruent only if the relaxed clock is assumed.

The F-type complete mitochondrial genome of Chinese freshwater mussels Cuneopsis pisciculus

The complete F-type mitochondrial genome of Cuneopsis pisciculus was first determined in this research. This circular genome (15,907 bp in size) contained 13 protein-coding genes, 22 tRNA genes and 2rRNA genes. The protein-coding genes initiated with orthodox ATG or ATA start codon. Twenty-five non-coding regions were found throughout the mitochondrial genome of C. pisciculus, ranging in size from 1 to 326 bp. The maximum was between tRNA(Gln) and gene ND5 (326 bp).

Complete F-type mitochondrial genome of freshwater mussels Unio douglasiae

The complete F-type mitochondrial genome of Unio douglasiae is reported for the first time in this research. It is 15,767 bp in size and contains 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. Except for COI, ND2, CYTB, ND1, ND6, ND4L, ATP8 with ATA start codon, and ND4 with ATT start codon, the remaining protein-coding genes initiate with the orthodox ATG start codon. There are 25 noncoding regions found throughout the mitogenome of Unio douglasiae, ranging in the size from 1 to 312 bp, the largest of which is between tRNA^Glu and ND2 (312 bp).

The complete mitochondrial genome of Chinese pond mussel Sinanodonta woodiana (Unionoida: Unionidae)

The complete mitochondrial genome of Sinanodonta woodiana is a circular molecule of 16256 bp in length, containing 14 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and 2 control regions. The A + T content of the overall base composition of H-strand is 65.9% (T: 28.0%; C: 22.3%; A: 37.9%; G: 11.8%). F ORF (Female-specific open reading frame) begins with ATA, Cyt b begins with ATC, ATP6, ATP8, COII, COIII, ND1, ND2, ND3 and ND5 begin with ATG, ND4L begins with GTG, COI begins with TTG, and other two protein-coding genes begin with ATT as start codon. COII, COIII, F ORF, ND1, ND3, ND5 and ND6 genes are terminated with TAA as stop codon, ATP6, ATP8, COI, Cyt b, ND4 and ND4L end with TAG, and ND2 gene ends with T.

Complete F-type mitochondrial genome of freshwater mussel Lanceolaria glayana

Lanceolaria glayana is widely distributed in China as a common species of freshwater pearl mussel. The complete F-type mitochondrial genome of L. glayana was firstly reported by this research. This circle genome (15,736 bp in size) contains 13 protein-coding genes, 22 tRNA genes, 2rRNA genes. Most of the protein-coding genes initiated with the orthodox ATG start codon while COI and ND4 were corresponding to ATA start codon. Further, 28 non-coding regions were found throughout the mitogenome of L. glayana, ranging from 1 to 310 bp. The maximum is between ND5 and tRNA Gln (310 bp), which is longer than the control region sequences (about 280 bp) of freshwater mussels from Poland and South Korea.

Complete F-type mitochondrial genome of Chinese freshwater mussels Lamprotula gottschei

Lamprotula gottschei (Von Martens, 1894) is one of the freshwater pearl mussels which have great meaning on economy in China. The complete F-type mitochondrial genome of L. gottschei was firstly determined. This circle genome (15,915 bp in size) comprises 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, except for COI, CYTB, ND4, ND6 and ND4L, with ATA start codon, the remaining protein-coding genes initiated with the orthodox ATG start codon. There were 25 non-coding regions found throughout the mitogenome of L. gottschei, ranging in size from 1 to 305 bp, the largest of which is between ND2 and tRNA(Gln) (305 bp), and is longer than the control region sequences (about 280 bp) of freshwater mussels from Poland and South Korea. Strictly maternal inheritance (SMI) of mitochondrial DNA (mtDNA) is the standard rule in nearly all anisogamic organisms (SMI), (Roze et al., 2005; White et al., 2008). However, bivalves of the marine family Mytilidae and Veneroida as well as the freshwater family Unionidae possess two independently inheritance mode of mitochondrial DNA (mtDNA) (Boyle & Etter, 2013), which involves maternal (F) and paternal (M) transmission routes concomitant with highly divergent gender-associated mtDNA genomes. This system of mtDNA transmission was called doubly uniparental inheritance (DUI) (Zouros, 2000).

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.

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.

Complete F-type mitochondrial genome of Chinese freshwater mussel Lamprotula tortuosa

Lamprotula tortuosa is the endemic species of freshwater pearl mussel in China. The complete F-type mitochondrial genome of L. tortuosa was first determined. This circle genome (15,722 bp in size) comprises 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, except for COI, CYTB, ND6, and ND4L, with TTG, ATT, ATT, and ATA start codon, respectively; the remaining protein-coding genes initiated with the orthodox ATG start codon. There were 25 noncoding regions found throughout the mitogenome of L. tortuosa, ranging in size from 1 to 241 bp, the largest of which is between ND5 and tRNA(Gln) (241 bp), but is shorter than the control region sequences(about 280 bp) of freshwater mussels from Poland and South Korea.

Complete mitochondrial genome of Coelomactra antiquata (Mollusca: Bivalvia): The first representative from the family Mactridae with novel gene order and unusual tandem repeats

The complete mitochondrial genome plays an important role in the accurate inference of phylogenetic relationships among metazoans. Mactridae, also known as trough shells or duck clams, is an important family of marine bivalve clams in the order Veneroida. Here we present the complete mitochondrial genome sequence of the Xishishe Coelomactra antiquata (Mollusca: Bivalvia), which is the first representative from the family Mactridae. The mitochondrial genome of C. antiquata is of 17,384bp in length, and encodes 35 genes, including 12 protein-coding, 21 transfer RNA, and 2 ribosomal RNA genes. Compared with the typical gene content of animal mitochondrial genomes, atp8 and tRNAS(2) are missing. Gene order of the mitochondrial genome of C. antiquata is unique compared with others from Veneroida. In the mitochondrial genome of the C. antiquata, a total of 2189bp of non-coding nucleotides are scattered among 26 non-coding regions. The largest non-coding region contains one section of tandem repeats (99 bp×11), which is the second largest tandem repeats found in the mitochondrial genomes from Veneroida. The phylogenetic trees based on mitochondrial genomes support the monophyly of Veneridae and Lucinidae, and the relationship at the family level: ((Veneridae+Mactridae)+(Cardiidae+Solecurtidae))+Lucinidae. The phylogenetic result is consistent with the morphological classification. Meanwhile, bootstrap values are very high (BP=94-100), suggesting that the evolutionary relationship based on mitochondrial genomes is very reliable.

Sequencing technologies and genome sequencing

The high-throughput - next generation sequencing (HT-NGS) technologies are currently the hottest topic in the field of human and animals genomics researches, which can produce over 100 times more data compared to the most sophisticated capillary sequencers based on the Sanger method. With the ongoing developments of high throughput sequencing machines and advancement of modern bioinformatics tools at unprecedented pace, the target goal of sequencing individual genomes of living organism at a cost of $1,000 each is seemed to be realistically feasible in the near future. In the relatively short time frame since 2005, the HT-NGS technologies are revolutionizing the human and animal genome researches by analysis of chromatin immunoprecipitation coupled to DNA microarray (ChIP-chip) or sequencing (ChIP-seq), RNA sequencing (RNA-seq), whole genome genotyping, genome wide structural variation, de novo assembling and re-assembling of genome, mutation detection and carrier screening, detection of inherited disorders and complex human diseases, DNA library preparation, paired ends and genomic captures, sequencing of mitochondrial genome and personal genomics. In this review, we addressed the important features of HT-NGS like, first generation DNA sequencers, birth of HT-NGS, second generation HT-NGS platforms, third generation HT-NGS platforms: including single molecule Heliscope™, SMRT™ and RNAP sequencers, Nanopore, Archon Genomics X PRIZE foundation, comparison of second and third HT-NGS platforms, applications, advances and future perspectives of sequencing technologies on human and animal genome research.