Complete mitochondrial DNA sequence of Crassostrea nippona: comparative and phylogenomic studies on seven commercial Crassostrea species

Phylogenetic relationships and adaptation in deep-sea carideans revealed by mitogenomes

The deep-sea environment is characterized by extreme and inhospitable conditions, including oxygen depletion, low temperatures, high pressure, absence of light, and limited food availability. Mitochondria and mitogenomes play a crudial role in aerobic respiration to generate energy for eukaryotes. Here, using the Illumina Hiseq 4000 platform, we performed mitogenome sequencing for five deep-sea caridean species: Lebbeus shinkaiae, Lebbeus Formosus, Glyphocrangon regalis, Heterocarpus dorsalis, and Heterocarpus laevigatus, and five deep-sea caridean mitogenomes were assembled and identified. Each of the five mitogenomes contained 13 protein-coding genes, 2 rRNAs and 22 tRNAs. Specific elements, such as tandem repeats and AT-rich sequences, were observed in the control regions of Lebbeus formosus and Lebbeus shinkaiae, potentially take a role in regulating mitochondrial genome replication and transcription. The gene order of all obtained mitogenomes follows caridean ancestral type organization. Phylogenetic analysis shows a robustly supported phylogenetic tree for the infraorder Caridea. The monophyly of the families included in this study was strongly supported. This study supports the monophyly of Oplophoroidea, but rejects the monophyletic status of Nematocarcinoidea, Crangonoidea, and Alpheoidea. At the genus level, Plesionika is polyphyletic and Rimicaris is paraphyletic in our analysis. Furthermore, Paralebbeus may be considered invalid and synonymous with Lebbeus. Positive selection analysis reveals evidence for adaptive changes in the mitogenome of different deep-sea caridean lineages. Nine residues located in cox1, cox3, atp6, nad1, nad2, nad4, nad5, nad6 and cytb were determined to have undergone positive selection. Mitogenome of different deep-sea lineages experienced different positive selection, and the lineage represented by Alvinocarididae living in deep-sea hydrothermal vents experienced the strongest positive selection. This study provides valuable insights into the adaptive evolution of deep-sea shrimps at the mitochondrial, highlighting the mitogenomic strategy that contribute to their unique adaptations in the deep-sea environment.

The First Complete Mitochondrial Genome of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) and Phylogenetic Relationships within Infraorder Brachyura

Characterizing the complete mitochondrial genome (mitogenome) of an organism is useful for genomic studies in taxonomy and evolution. The mitogenomic characteristics of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) have never been studied. The present study decodes the first mitogenome of E. crenata by high-throughput sequencing (HTS). The length of the mitogenome is 15,597 bp, and it contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL), and 22 transfer RNA genes. There are 14 and 23 genes observed on the heavy and light strands, respectively. E. crenata possesses a trnH-cac translocation, with the trnH-cac shifted between trnE-gaa and trnF-ttc instead of the usual location between nad5 and nad4 in decapods. Phylogenetic analyses based on the current dataset of 33 Brachyuran mitogenomes indicate that E. crenata. is closely related to Ashtoret lunaris of Matutidae. The similar codon usage and rearrangements in the two species provide evidence for their close phylogenetic relationship. Positive selection analysis showed that one residue located in cox1 was identified as a positively selected site with high BEB value (>95%), indicating that this gene was under positive selection pressure. This study is the first complete mitogenome record for the family Goneplacidae, and the results obtained may improve the understanding of the phylogeny of Goneplacidae in Brachyura.

Complete mitochondrial genome of Cultellus attenuatus and its phylogenetic implications

BACKGROUND

The mitochondrial genomes of three species in Solenoidea of Heterodonta have been reported, but the mitochondrial genes and phylogenetic relationships of Cultellus attenuatus, which also belongs to this superfamily and has high economic value, are unknown.

METHODS AND RESULTS

The complete mitochondrial genome of C. attenuatus was sequenced and compared with mitogenomes of seven species of Heterodonta bivalve mollusks in GenBank. The mitochondrial genome of C. attenuatus has a length of 16,888 bp and contains 36 genes, including 12 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. In comparison with C. attenuates, the mitochondrial genes of Sinonovacula constricta from the same family were not rearranged, but those of six other species from different families were rearranged to different degrees. The location, size, and composition of the largest noncoding regions in eight species suggested a closer relationship between C. attenuatus and S. constricta. The phylogenetic analysis showed that C. attenuatus and S. constricta belonging to Cultellidae cluster into one branch and that two species of Solenidae (Solen grandis and Solen strictus) clustered as their sister taxa.

CONCLUSIONS

Overall, we used mitochondrial genome data to demonstrate that C. attenuatus and S. constricta exhibit the closest relationship in Heterodonta. These data and analyses provide new insights into the phylogenetic relationships in Heterodonta.

Adaptive divergence and underlying mechanisms in response to salinity gradients between two Crassostrea oysters revealed by phenotypic and transcriptomic analyses

Comparing the responses of closely related species to environmental changes is an efficient method to explore adaptive divergence, for a better understanding of the adaptive evolution of marine species under rapidly changing climates. Oysters are keystone species thrive in intertidal and estuarine areas where frequent environmental disturbance occurs including fluctuant salinity. The evolutionary divergence of two sister species of sympatric estuarine oysters, Crassostrea hongkongensis and Crassostrea ariakensis, in response to euryhaline habitats on phenotypes and gene expression, and the relative contribution of species effect, environment effect, and their interaction to the divergence were explored. After a 2-month outplanting at high- and low-salinity locations in the same estuary, the high growth rate, percent survival, and high tolerance indicated by physiological parameters suggested that the fitness of C. ariakensis was higher under high-salinity conditions and that of C. hongkongensis was higher under low-salinity conditions. Moreover, a transcriptomic analysis showed the two species exhibited differentiated transcriptional expression in high- and low-salinity habitats, largely caused by the species effect. Several of the important pathways enriched in divergent genes between species were also salinity-responsive pathways. Specifically, the pyruvate and taurine metabolism pathway and several solute carriers may contribute to the hyperosmotic adaptation of C. ariakensis, and some solute carriers may contribute to the hypoosmotic adaptation of C. hongkongensis. Our findings provide insights into the phenotypic and molecular mechanisms underlying salinity adaptation in marine mollusks, which will facilitate the assessment of the adaptive capacity of marine species in the context of climate change and will also provide practical information for marine resource conservation and aquaculture.

Immunomodulatory activity of aqueous extract from Crassostrea sikamea in the splenocytes of Sprague-Dawley rats

Crassostrea sikamea (C. sikamea) is used as an important edible and medicinal seafood in China. In the present study, an aqueous extract of C. sikamea (AECs) was prepared, and its immunomodulatory effects on rat splenocytes were studied. 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium (MTS) assay revealed that AECs was able to promote splenocyte proliferation. Moreover, flow cytometry revealed that AECs treatment markedly altered the populations of splenic lymphocyte subtypes. Data from real-time quantitative PCR (RT-qPCR) and enzyme-linked immunosorbent assay (ELISA) showed that AECs promoted the mRNA expression and secretion of TNF-α, IL-2, IL-6, IL-12, and IFN-γ. Mechanistically, p38 MAPK phosphorylation in splenocytes was significantly upregulated under AECs treatment and p38 MAPK inhibitor reversed the promoting effect of AECs on the expression of inflammatory cytokines. Collectively, our novel evidence suggests that AECs exhibits immunomodulatory activity in vitro, supporting the further application of C. sikamea as a potential functional food.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Immunomodulatory effects of flazin from Crassostrea sikamea on splenic lymphocytes of Sprague-Dawley rats

Crassostrea sikamea (C.sikamea) is an important edible and medicinal seafood in China. In the present study, a compound named flazin was separated and identified from the ethyl acetate extract of C.sikamea (EAECs) for the first time. In addition, the 3-(4, 5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetra zolium (MTS) assay revealed that EAECs and flazin inhibited the transformation of splenic lymphocytes in vitro. Moreover, flazin (20 μg·mL^-1) altered the populations of splenic lymphocyte subtypes. Real-time quantitative PCR (RT-qPCR) analysis and enzyme-linked immunosorbent assay (ELISA) showed that flazin suppressed the mRNA expression and secretion of TNF-α and IL-2, and reversed Concanavalin A (ConA)-induced mRNA up-regulation and protein secretion of TNF-α and IL-2. Western blot results showed that flazin reversed ConA-induced increases in p-ERK1/2 and p-p38 in splenocytes. In conclusion, flazin exhibits effective immunomodulatory function and may be useful for treating immune-related disorders, which indicates the application potential of C.sikamea as a functional food or immunomodulator.

The complete mitogenome of Phymorhynchus sp. (Neogastropoda, Conoidea, Raphitomidae) provides insights into the deep-sea adaptive evolution of Conoidea

The deep-sea environment is characterized by darkness, hypoxia, and high hydrostatic pressure. Mitochondria play a vital role in energy metabolism; thus, they may endure the selection process during the adaptive evolution of deep-sea organisms. In the present study, the mitogenome of Phymorhynchus sp. from the Haima methane seep was completely assembled and characterized. This mitogenome is 16,681 bp in length and contains 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of most sequenced conoidean gastropods. Some special elements, such as tandem repeat sequences and AT-rich sequences, which are involved in the regulation of the replication and transcription of the mitogenome, were observed in the control region. Phylogenetic analysis revealed that Conoidea is divided into two separate clades with high nodal support. Positive selection analysis revealed evidence of adaptive changes in the mitogenomes of deep-sea conoidean gastropods. Eight residues located in atp6, cox1, cytb, nad1, nad4, and nad5 were determined to have undergone positive selection. This study explores the adaptive evolution of deep-sea conoidean gastropods and provides valuable clues at the mitochondrial level regarding the exceptional adaptive ability of organisms in deep-sea environments.

Revision shock in Pacific oysters taxonomy: the genus Magallana (formerly Crassostrea in part) is well-founded and necessary

The description of the genus Magallana provoked taxonomic ‘revision shock’ among Bayne and colleagues, who criticised: (1) the supporting evidence for this taxonomic change, (2) the procedure used and (3) its impact on the nomenclatural stability of oysters. Here, we demonstrate that the description of the genus Magallana fulfils the norms of taxonomic and indexing revisions, and is well-founded on a scientific basis. The clade named Magallana is supported by a robust phylogeny based on comprehensive taxon sampling, independent datasets and varied analytical methods. A new maximum likelihood phylogenetic analysis of oysters, based on mitochondrial gene order data (representing molecular phenotypes above the sequence level), has provided additional support for this classification. These results are fully consistent with all previously published phylogenetic studies, thus providing an unambiguous indication of the stability of the clade that meets all the currently accepted criteria for naming clades as taxa. We show that, while the criterion of ‘morphological diagnosability’ is not applicable at any level of oyster classification, several striking molecular phenotypes are diagnostic of Magallana, both at the nuclear and the mitochondrial genome level. The classification with three genera Magallana, Crassostrea and Talonostrea reflects the evolutionary diversity of Crassostreinae and is consistent with taxonomic ranking criteria adopted for other oyster subfamilies.

First de novo transcriptome assembly of Iwagaki oyster, Crassostrea nippona, and comparative evolutionary analysis of salinity-stress response genes in Crassostrea oysters

Crassostrea nippona is a commercially important oyster species in East Asia for it is edible during the summer when the other oyster species are unavailable. Salinity is one of the important limiting factors to the survival and distribution of this stenohaline species. In this study, 535 million reads (74G data) from C. nippona were produced and assembled into 66,742 transcripts. The number of 19,253 differentially expressed genes (DEGs) under salinity stress were identified as salinity stress-response genes. Through comparative evolutionary analysis in five Crassostrea species from East Asia, salinity stress-response genes were noticed to have higher adaptive evolution rates than other genes. This study presents the first de novo transcriptome of C. nippona. Furthermore, comparative evolutionary analysis implies that salinity plays an important role in speciation of Crassostrea species.

The complete mitochondrial genome of Calyptogena marissinica (Heterodonta: Veneroida: Vesicomyidae): Insight into the deep-sea adaptive evolution of vesicomyids

The deep-sea chemosynthetic environment is one of the most extreme environments on the Earth, with low oxygen, high hydrostatic pressure and high levels of toxic substances. Species of the family Vesicomyidae are among the dominant chemosymbiotic bivalves found in this harsh habitat. Mitochondria play a vital role in oxygen usage and energy metabolism; thus, they may be under selection during the adaptive evolution of deep-sea vesicomyids. In this study, the mitochondrial genome (mitogenome) of the vesicomyid bivalve Calyptogena marissinica was sequenced with Illumina sequencing. The mitogenome of C. marissinica is 17,374 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes. All of these genes are encoded on the heavy strand. Some special elements, such as tandem repeat sequences, “G(A)_nT” motifs and AT-rich sequences, were observed in the control region of the C. marissinica mitogenome, which is involved in the regulation of replication and transcription of the mitogenome and may be helpful in adjusting the mitochondrial energy metabolism of organisms to adapt to the deep-sea chemosynthetic environment. The gene arrangement of protein-coding genes was identical to that of other sequenced vesicomyids. Phylogenetic analyses clustered C. marissinica with previously reported vesicomyid bivalves with high support values. Positive selection analysis revealed evidence of adaptive change in the mitogenome of Vesicomyidae. Ten potentially important adaptive residues were identified, which were located in cox1, cox3, cob, nad2, nad4 and nad5. Overall, this study sheds light on the mitogenomic adaptation of vesicomyid bivalves that inhabit the deep-sea chemosynthetic environment.

Complete mitochondrial genome of the first deep-sea spongicolid shrimp Spongiocaris panglao (Decapoda: Stenopodidea): Novel gene arrangement and the phylogenetic position and origin of Stenopodidea

Stenopodidea Claus, 1872 (Crustacea: Decapoda) is one of the major groups of decapods crustaceans. Hitherto, only one complete mitochondrial genome (mitogenome) from the family Stenopodidae is available for the infraorder Stenopodidea. Here, we determined the complete mitogenome of Spongiocaris panglao de Grave and Saito, 2016 using Illumina sequencing, representing the first species from the family Spongicolidae. The 15,909 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. Although the overall genome organization is typical for metazoans, the mitogenome of S. panglao shows some derived characters. A + T content of 77.42% in S. pamglao mitogenome is second-highest among the dacapods described to date. The trnR gene exhibit modified secondary structure with the TψC loop completely missing, which might be a putative autapomorphy of S. pamglao mitogenome. Compared with the shallow-water stenopodidean species S. hispidus, the control region of S. pamglao exhibits three characteristics: larger size, higher A + T content, and more tandem repeat sequences. The gene order exhibited difference from the ancestral mitogenome pattern of the Pancrustacea, with 5 tRNA genes rearrangement. The result from BI was agreed with most morphological characters and molecular evidences, revealing that Stenopodidea and Reptantia had the closest relationship, as the sister group of Caridea. Still, the alternative hypothesis supported from ML topology cannot be completely rejected based on the current data. Estimated times revealed that the two stenopodideans families Stenopodidae and Spongicolidae diverged from each other around 122 Mya. The divergence time of spongicolid shrimp is in good agreement with the origin of their hexactinellid hosts (78-144 Mya).

The complete mitochondrial genome of the alvinocaridid shrimp Shinkaicaris leurokolos (Decapoda, Caridea): Insight into the mitochondrial genetic basis of deep-sea hydrothermal vent adaptation in the shrimp

Deep-sea hydrothermal vent is one of the most extreme environments on Earth with low oxygen and high levels of toxins. Decapod species from the family Alvinocarididae have colonized and successfully adapted to this extremely harsh environment. Mitochondria plays a vital role in oxygen usage and energy metabolism, thus it may be under selection in the adaptive evolution of the hydrothermal vent shrimps. In this study, the mitochondrial genome (mitogenome) of alvinocaridid shrimp Shinkaicaris leurokolos (Kikuchi & Hashimoto, 2000) was determined through Illumina sequencing. The mitogenome of S. leurokolos was 15,903bp in length, containing 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of sequenced alvinocaridids. It has the longest concatenated sequences of protein-coding genes, tRNAs and shortest pooled rRNAs among the alvinocaridids. The control regions (CRs) of alvinocaridid were significantly longer (p<0.01) than those of the other caridaen. The alignment of the alvinocaridid CRs revealed two conserved sequence blocks (CSBs), and each of the CSBs included a noncanonical open reading frame (ORF), which may be involved in adjusting mitochondrial energy metabolism to adapt to the hydrothermal environment. Phylogenetic analysis supported that the deep-sea hydrothermal vent shrimps may have originated from those living in shallow area. Positive selection analysis reveals the evidence of adaptive change in the mitogenome of Alvinocarididae. Thirty potentially important adaptive residues were identified, which were located in atp6, cox1, cox3, cytb and nad1-5. This study explores the mitochondrial genetic basis of hydrothermal vent adaptation in alvinocaridid for the first time, and provides valuable clues regarding the adaptation.

Studying the genetic basis of speciation in high gene flow marine invertebrates

A growing number of genes responsible for reproductive incompatibilities between species (barrier loci) exhibit the signals of positive selection. However, the possibility that genes experiencing positive selection diverge early in speciation and commonly cause reproductive incompatibilities has not been systematically investigated on a genome-wide scale. Here, I outline a research program for studying the genetic basis of speciation in broadcast spawning marine invertebrates that uses a priori genome-wide information on a large, unbiased sample of genes tested for positive selection. A targeted sequence capture approach is proposed that scores single-nucleotide polymorphisms (SNPs) in widely separated species populations at an early stage of allopatric divergence. The targeted capture of both coding and non-coding sequences enables SNPs to be characterized at known locations across the genome and at genes with known selective or neutral histories. The neutral coding and non-coding SNPs provide robust background distributions for identifying FST-outliers within genes that can, in principle, identify specific mutations experiencing diversifying selection. If natural hybridization occurs between species, the neutral coding and non-coding SNPs can provide a neutral admixture model for genomic clines analyses aimed at finding genes exhibiting strong blocks to introgression. Strongylocentrotid sea urchins are used as a model system to outline the approach but it can be used for any group that has a complete reference genome available.

The complete mitochondrial genome of Crassostrea gasar (Bivalvia: Ostreidae)

The complete mitochondrial genome of Crassostrea gasar was sequenced using the Ion Proton technology in combination with 454 Roche GS-FLX plataform data. We assembled a 17,686 bp complete circular mitochondrial genome, containing 13 protein-coding genes, a major non-coding region (MNR), two ribosomal RNA genes and 24 transfer RNA genes. Phylogenetic analysis of concatenated amino acid sequences from mitochondria showed monophyletic clades formed with high bootstrap values. This is the first complete mitochondrial sequence of an oyster from South America. Mitogenome sequence was deposited in GenBank under the accession number KR856227.

Intraspecific Variation in Mitogenomes of Five Crassostrea Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China's coast, which are shaped by unknown mechanisms in a north-south divide.

The complete mitochondrial genome of Neptunea arthritica cumingii Crosse, (Gastropoda: Buccinidae)

The complete mitochondrial (mt) genome of the Neptune whelk, Neptunea arthritica cumingii, was determined using genome walking techniques in this study. The total length of the mt genome sequence of N. arthritica cumingii was 15 256 bp, including 13 protein-coding genes, 21 transfer RNA genes and 2 ribosomal RNA genes. The overall composition of the mitogenome was estimated to be 30.85% for A, 38.59% for T, 15.15% for C and 15.40% for G, indicating that an A + T (69.44%)-rich feature occurs in the N. arthritica cumingii mitogenome. The phylogenetic relationships of 11 mollusc species were constructed based on the complete mtDNA sequences by the neighbour-joining method using MEGA 5.0 software (MEGA Inc., Englewood, NJ).

Complete mitochondrial DNA sequences of Saccostrea mordax and Saccostrea cucullata: genome organization and phylogeny analysis

Classified in the phylum mollusks, oysters are bivalves which are found in estuaries and coastal zones. Because of their plastic shell, mitochondrial DNA analysis of this species becomes an interesting field, necessary to investigate their phylogenetic and evolution of relations. In our study, two oyster species: Saccostrea mordax and Saccostrea cucullata from Indian Ocean (Madagascar) were investigated. The complete sequence of Saccostrea mordax (16 512 bp) and Saccostrea cucullata (16 396 bp) were described and determined, with their mitogenomes deposited in the GenBank with accession number KP769562 and KP967577 respectively. Both mitochondrial genome sequences contained 12 protein-coding genes, 23 tRNAs, and two rRNAs, all encoded in the same heavy strand. High levels of similarity in the gene arrangement of the two Saccostrea species were evident. The phylogenetic analysis shows a closer relationship between the two Saccostrea species and confirms the strong relationship within Saccostrea, Crassostrea and Ostrea genus in taxonomy of Ostreidae family.

SNP mining in Crassostrea gigas EST data: transferability to four other Crassostrea species, phylogenetic inferences and outlier SNPs under selection

Oysters, with high levels of phenotypic plasticity and wide geographic distribution, are a challenging group for taxonomists and phylogenetics. Our study is intended to generate new EST-SNP markers and to evaluate their potential for cross-species utilization in phylogenetic study of the genus Crassostrea. In the study, 57 novel SNPs were developed from an EST database of C. gigas by the HRM (high-resolution melting) method. Transferability of 377 SNPs developed for C. gigas was examined on four other Crassostrea species: C. sikamea, C. angulata, C. hongkongensis and C. ariakensis. Among the 377 primer pairs tested, 311 (82.5%) primers showed amplification in C. sikamea, 353 (93.6%) in C. angulata, 254 (67.4%) in C. hongkongensis and 253 (67.1%) in C. ariakensis. A total of 214 SNPs were found to be transferable to all four species. Phylogenetic analyses showed that C. hongkongensis was a sister species of C. ariakensis and that this clade was sister to the clade containing C. sikamea, C. angulata and C. gigas. Within this clade, C. gigas and C. angulata had the closest relationship, with C. sikamea being the sister group. In addition, we detected eight SNPs as potentially being under selection by two outlier tests (fdist and hierarchical methods). The SNPs studied here should be useful for genetic diversity, comparative mapping and phylogenetic studies across species in Crassostrea and the candidate outlier SNPs are worth exploring in more detail regarding association genetics and functional studies.

The complete mitochondrial genome of Rapana venosa (Gastropoda, Muricidae)

The complete mitochondrial (mt) genome of the veined rapa whelk, Rapana venosa, was determined using genome walking techniques in this study. The total length of the mt genome sequence of R. venosa was 15,271 bp, which is comparable to the reported Muricidae mitogenomes to date. It contained 13 protein-coding genes, 21 transfer RNA genes, and two ribosomal RNA genes. A bias towards a higher representation of nucleotides A and T (69%) was detected in the mt genome of R. venosa. A small number of non-coding nucleotides (302 bp) was detected, and the largest non-coding region was 74 bp in length.

Complete mitochondrial genome of Ostrea denselamellosa (Bivalvia, Ostreidae)

The complete mitochondrial (mt) genome of the flat oyster, Ostrea denselamellosa, was determined using Long-PCR and genome walking techniques in this study. The total length of the mt genome sequence of O. denselamellosa was 16,227 bp, which is the smallest reported Ostreidae mt genome to date. It contained 12 protein-coding genes (lacking of ATP8), 23 transfer RNA genes, and two ribosomal RNA genes. A bias towards a higher representation of nucleotides A and T (60.7%) was detected in the mt genome of O. denselamellosa. The rrnL was split into two fragments (3' half, 711 bp; 5' half, 509 bp), which seems to be the unique characteristics of Ostreidae mt genomes.

A unique tRNA gene family and a novel, highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS(-AGY) and trnQ genes, a trnM isomerism, but it lacks trnS(-CUN). We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP→trnS(-AGY)) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, "Oxidored_q1_N", which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.

New features of Asian Crassostrea oyster mitochondrial genomes: a novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.