International Journal of Biological Sciences

Impact factor
3.873

ISSN 1449-2288

News feeds of IJBS published articles
My Manuscript
My Account

Journal of Biomedicinenew

Theranostics

International Journal of Medical Sciences

Journal of Cancer

Oncomedicine

Journal of Genomics

Journal of Bone and Joint Infection (JBJI)

Nanotheranostics

PubMed Central Indexed in Journal Impact Factor

Int J Biol Sci 2014; 10(5):479-489. doi:10.7150/ijbs.8739

Research Paper

Complete Mitochondrial Genomes of Carcinoscorpius rotundicauda and Tachypleus tridentatus (Xiphosura, Arthropoda) and Implications for Chelicerate Phylogenetic Studies

Su Youn Baek1,#, Eun Hwa Choi 2,#, Kuem Hee Jang2, Shi Hyun Ryu2, Sang Myeon Park3, Ho Young Suk4, Cheon Young Chang5, Ui Wook Hwang2, Corresponding address

1. School of Life Sciences, Graduate School, Kyungpook National University, Daegu 702-701, South Korea,
2. Department of Biology, Teachers College & Institute for Phylogenomics and Evolution, Kyungpook National University, Daegu 702-701, South Korea,
3. Department of Science Education, Graduate School, Kyungpook National University, Daegu 702-701, South Korea,
4. Department of Life Sciences, Yeungnam University, Gyeongsan, Gyeongsangbuk-do 705-717, South Korea,
5. Department of Biology, College of Natural Science, Daegu University, Gyeongsan, Gyeongsangbuk-do 712-714, South Korea
# Equally contributed.

This is an open access article distributed under the terms of the Creative Commons Attribution (CC BY-NC) License. See http://ivyspring.com/terms for full terms and conditions.
How to cite this article:
Baek SY, Choi EH, Jang KH, Ryu SH, Park SM, Suk HY, Chang CY, Hwang UW. Complete Mitochondrial Genomes of Carcinoscorpius rotundicauda and Tachypleus tridentatus (Xiphosura, Arthropoda) and Implications for Chelicerate Phylogenetic Studies. Int J Biol Sci 2014; 10(5):479-489. doi:10.7150/ijbs.8739. Available from http://www.ijbs.com/v10p0479.htm

Abstract

Horseshoe crabs (order Xiphosura) are often referred to as an ancient order of marine chelicerates and have been considered as keystone taxa for the understanding of chelicerate evolution. However, the mitochondrial genome of this order is only available from a single species, Limulus polyphemus. In the present study, we analyzed the complete mitochondrial genomes from two Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus to offer novel data for the evolutionary relationship within Xiphosura and their position in the chelicerate phylogeny. The mitochondrial genomes of C. rotundicauda (15,033 bp) and T. tridentatus (15,006 bp) encode 13 protein-coding genes, two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Overall sequences and genome structure of two Asian species were highly similar to that of Limulus polyphemus, though clear differences among three were found in the stem-loop structure of the putative control region. In the phylogenetic analysis with complete mitochondrial genomes of 43 chelicerate species, C. rotundicauda and T. tridentatus were recovered as a monophyly, while L. polyphemus solely formed an independent clade. Xiphosuran species were placed at the basal root of the tree, and major other chelicerate taxa were clustered in a single monophyly, clearly confirming that horseshoe crabs composed an ancestral taxon among chelicerates. By contrast, the phylogenetic tree without the information of Asian horseshoe crabs did not support monophyletic clustering of other chelicerates. In conclusion, our analyses may provide more robust and reliable perspective on the study of evolutionary history for chelicerates than earlier analyses with a single Atlantic species.

Keywords: Carcinoscorpius rotundicauda, Tachypleus tridentatus, Horseshoe crabs, Xiphosura, Mitochondrial genome, Phylogenetics

Introduction

Horseshoe crabs (Xiphosura, Arthorpoda) are marine chelicerates dwelling shallow waters on sandy or muddy bottoms. The Limulidae, only living family of Xiphosura, contains four horseshoe crab species, Tachypleus tridentatus Leach, Tachypleus gigas Müller, Carcinoscorpius rotundicauda Latreille and Limulus polyphemus Linnaeus [1]. The Atlantic horseshoe crab, L. polyphemus, occurs along the eastern coast of North America, while the remaining three species are widely distributed across Asia from Southeast to Northeast [1]. The horseshoe crabs have been dramatically decreased in population sizes and many populations have extirpated throughout much of their ranges in the world, probably due to the global degradation of coastal areas.

Despite its name, horseshoe crabs are rather more closely related to arachnids (subphylum Chelicerata) than to crabs (subphylum Crustacea) [1-3]. They are also referred to as living fossil, since little change has been made in their external morphologies over the last 450 million years [1-3]. Phylogenetic relationship among four living horseshoe crabs remains scientifically unresolved yet. Two Tachypleus species have been known to compose a monophyletic assemblage on the basis of morphological evaluation [4]. Amino acid sequencing [5] and interspecific crossing experiments [6] showed rather different results suggesting that T. tridentatus should be more closely related to C. rotundicauta than to T. gigas. However, studies based on different mtDNA regions failed to provide a congruent pattern; for example, T. gigas and C. rotundicauda were clustered as a monophyly in the combined analyses using 16S ribosomal RNA (16S rRNA) and cytochrome oxidase subunit I (COI) [7], while T. gigas formed the closest relationship with T. tridentatus in the combined analysis using 18S rRNA, 28S rRNA and COI [8].

Xiphosura has been considered as an ancestral taxon in chelicerates as well as arthropods [9, 10], and the complete mitochondrial genome of L. polyphemus has frequently been used to resolve the phylogenetic relationships among chelicerates [3, 11]. However, additional mitochondrial genomic information of horseshoe crab species is critically required to reconstruct the robust evolutionary relationship among chelicerates as well as arthropods, because the phylogenetic relationship between L. polyphemus and the remaining Asian species is still vague. Here, the complete mitochondrial genomes of two Asian horseshoe crabs, C. rotundicauda and T. tridentatus were analyzed to provide novel data for the relationship among species within the order and the placement of horseshoe crabs in the phylogenetic tree of chelicerates.

Materials and methods

Sampling

Specimens of Carcinoscorpius rotundicauda was purchased from a pet shop (http://www.hanqua.co.kr, Republic of Korea). The tissue sample of Tachypleus tridentatus for analyses was given from Department of Biology, Daegu University (Daegu, Republic of Korea). The remaining tissue specimens were deposited in the Department of Biology, Teacher's College, Kyungpook National University (voucher numbers, HC-13-01: C. rotundicauda; HC-13-02: T. tridentatus). Total cellular DNA was extracted using the DNeasy Tissue Kit (QIAGEN Co., Germany) according to the manufacturer's protocol.

Amplification and sequencing

The mitochondrial genomes were amplified using long-range PCR [12] with primer sets shown in Table 1 and the Expand Long Template PCR Kit (Roche Co., Germany). The PCR setting was as followed: [92℃ for 2 min], [92℃ for 30 sec, 52℃ for 30 sec, 68℃ for 9 min] × 14, [92℃ for 10 sec, 52℃ for 30 sec, 68℃ for 9 min (+ 20 sec per cycle)] × 24, [68℃ for 5 min]. PCR reactants were loaded on a 1.0 % agarose gel and stained with ethidium bromide to visualize the bands on ultraviolet transilluminator. The PCR products were purified using the PCR Purification Kit (QIAGEN Co., Germany) and sequenced using the primer-walking strategy with the ABI PRISM BigDye terminator system and the ABI3700 model automatic sequencer (Genotech Co., Korea).

 Table 1 

Primers used in the long-range PCR for the mitochondrial genome analyses of Carcinoscorpius rotundicauda and Tachypleus tridentatus.

SegmentPrimersSequenceCitation
a. Carcinoscorpius rotundicauda
cox1 - rrnLHCO24985´-TAA ACT TCA GGG TGA CCA AAA AAT CA-3´34
16SB5´-CCG GTY TGA ACT CAR ATC A-3´35
rrnL - cox1LCO14905´-GGT CAA CAA ATC ATA AAG ATA TTG G-3'34
16SA5´-CGC CTG TTT AHC AAA AAC AT-3´35
b. Tachypleus tridentatus
cox1 - rrnLHCO24985´-TAA ACT TCA GGG TGA CCA AAA AAT CA-3´34
16SB5´-CCG GTY TGA ACT CAR ATC A-3´35
rrnL - cob16SA5´-CGC CTG TTT AHC AAA AAC AT-3´35
COB-F15'-CGA GTA ATT CAT GCA AAC GGA GC-3'13
cob - cox 1Tachy-CytbL5´-GCA GGA ACA GGA TGA ACA GT-3´This study
Tachy-CO1H5´-GCA GGA ACA GGA TGA ACA GT-3´This study

Sequence analysis

Thirteen protein-coding genes and two ribosomal RNA genes were identified based on sequence similarity under BLAST searches in the NCBI database. The boundary of each gene was determined by alignments with the sequences of the American horseshoe crab, L. polyphemus [13] in the Clustal X 2.1 [14]. The position and secondary structure of tRNA genes were determined by using tRNAscan-SE Search Server [15]. Ribosomal RNA (rRNA) genes were identified based on sequence similarity in BLAST searches. Alignments of protein coding genes were translated into amino acid sequences to confirm whether the amplified domains are functional with no frame-shifting or no premature stop codons. The complete mitochondrial genome sequences were deposited in the NCBI GenBank under the accession number JQ178358 (C. rotundicauda) and JQ739210 (T. tridentatus). The CG-skew values in both coding and non-coding regions (Table 2) were calculated based on CG - skew = (C - G)/(C+G) [11, 16].

Phylogenetic analysis

The mitochondrial genomes of C. rotundicauta and T. tridentatus were analyzed in the phylogenetic tree using a total of 41 other chelicerate mitochondrial genomes retrieved from NCBI GenBank (Supplementary Table 1). The nucleic acid sequences of 12 protein-coding genes were aligned using Bio-Edit sequence alignment editor (Ibis Biosciences, USA). The atp8 gene was excluded from the analyses, because the substitution rate was too high to recover the true phylogeny and too short to provide enough numbers of phylogenetically informative characters. Using the annotated gene boundary information, each protein coding gene was excised from the genomic sequence and put into an individual file. The program EMBOSS Transeq [17] was used to translate the nucleotide sequences into amino acid sequences based on the invertebrate mitochondrial genetic code. CLUSTAL X was finally used to align amino acid sequences of 12 genes. Only well-aligned and conserved alignment sites were extracted from each alignment subset using Gblock ver. 0.91b program [18] with the default setting. The extracted conserved blocks of amino acids were subsequently concatenated into a unified, single large alignment set.

The refined alignment was subjected to two different tree-making algorithms: the maximum likelihood (ML) and Bayesian inference (BI) methods. The best fitting model of sequence evolution was tested by ProtTest ver. 1.3 [19] under Akaike information criterion for the amino acid data. The mtREV [20] + I + Γ was consequently selected as the best fitting model and employed for the phylogenetic reconstruction. The ML analysis was carried out using PHYML v2.4.4 [21]. The bootstrap proportions (BPML; 1000 replicates) of the ML tree were obtained by the fast-ML method using PHYML. The BI analyses implemented in MrBayes v.3.1 [22] with 1,000,000 generations, four MCMC chains (one hot and three cold) and burn-in step of the first 1,000. Node confidence values of the tree were presented with Bayesian posterior probabilities (BPP).

Results

Genome composition

The mitochondrial genomes of C. rotundicauda (15,033 bp) and T. tridentatus (15,006 bp) included 13 protein coding genes (nad1 - 6, nad4L, cox1 - 3, cob, atp6 and atp8), two ribosomal RNA genes (rrnL and rrnS), 22 tRNA genes and one large non-coding region (putative control region, CR) (Fig. 1; Table 2). In both genomes, nine protein coding genes (nad2, cox1, cox2, atp8, atp6, cox3, nad3, nad6, and cob) were encoded on the heavy strand along with 13 tRNAs (trnI, trnM, trnW, trnK, trnD, trnG, trnA, trnR, trnN, trnS2, trnE, trnT, and trnS2), while the remaining four protein coding genes (nad5, nad4, nad4L, nad1) with nine tRNAs (trnQ, trnC, trnY, trnF, trnH, trnP, trnL1, trnL2, and trnV) and two rRNAs (rrnL and rrnS) were encoded on the light strand (Fig. 1; Table 2). Overall gene arrangement of both mitochondrial genomes was completely identical to that of L. polyphemus [13].

The gene components were rather loosely juxtaposed with 51/23 (C. rotundicauda) and 40/24 (T. tridentatus) of gap/overlapping nucleotides, considering that of L. polyphemus (23/19; Table 2) [13]. Although the overall A + T contents of 73.8 % in C. rotundicauda and 74.0 % in T. tridentatus were relatively higher than that of L. polyphemus (67.57 %), those values are within the range (60.2 - 80.4 %) of chelicerates (Supplementary Table 1). The overall pattern of nucleotide skew was highly similar among three mitochondrial genomes including that of L. polyphemus, with only an exception found on the putative control region (Table 3).

Protein coding genes

The inferred start/stop codons for protein coding genes of both C. rotundicauda and T. tridentatus are listed in Table 2. All of the protein coding genes in both mitochondrial genomes were initiated by ATN, with the exceptions in cox1, nad1 and nad5 (Table 2). The open reading frame of cox1 in both Asian horseshoe crabs started with TTA, while that of nad1 and nad5 started with TTG. The canonical stop codon (TAA or TAG) occurs in seven protein coding genes (nad1, nad2, cox1, atp6, nad3, nad4 and nad6; Table 2), while the remaining six (nad5, nad4L, cob, cox2, cox3 and atp8) had incomplete T- or TA- stop codons (Table 2). The codon usage pattern of the 13 protein coding genes is provided in Table 4. A + T rich codons, such as Leu (UUA), Ile (AUU) and Phe (UUU), are frequent in both C. rotundicauda and T. tridentatus (Table 4).

Transfer and ribosomal RNA genes

A total of 22 tRNA genes ranging from 62 to 73 bp in length were identified in mitochondrial genomes of both C. rotundicauda and T. tridentatus. The putative secondary structure for each tRNA gene could be predicted based on the sequences (Fig. 2). The sequences, anticodon nucleotides, and secondary structures of tRNA genes from both C. rotundicauda and T. tridentatus were very similar to those observed in L. polyphemus [13]. All of the tRNA genes were capable of forming typical cloverleaf secondary structures (dihydrouridine, DHU), with an exception of trnS1; trnS1 lacked DHU arm in T. tridentatus (Fig. 2B) and had a shortened stem (2 bp) in C. rotundicauda (Fig. 2A). Two rRNA genes (rrnS and rrnL) were encoded on the light strand and were separated by a trnV in both mitochondrial genomes analyzed. The sizes of rrnS and rrnL were estimated to be 816 bp (rrnS) and 1,301 bp (rrnL) for C. rotundicauda and 800 bp (rrnS) and 1,294 bp for T. tridentatus.

 Fig 1 

Gene map of the mitochondrial genome of Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus. Thirteen protein-coding genes and 2 ribosomal RNA genes are abbreviated as follows: nad1-6 and nad4L, NADH dehydrogenase subunits 1-6 and 4L; cox 1-3, cytochrome C oxidase subunits I-III; cob, cytochrome b apoenzyme; atp 6 and 8, ATPase subunits 6 and 8. The 22 transfer RNA genes are identified by the IUPAC amino acid single-letter codes.

Int J Biol Sci Image (Click on the image to enlarge.)
 Table 2 

The mitochondrial genome profile of two Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus

GeneStrandCarcinoscorpius rotundicauda (15,033 bp)Tachypleus tridentatus I (15,006 bp)
positionSizeCodonIntergenic
bp*
positionSizeCodonIntergenic
bp*
FromtobpStartStopFromtobpStartStop
trnIH16767-216767-2
trnQL6613267-26613267-2
trnMH131199690131199690
nad 2H20012191020ATTTAA-220012191020ATTTAA-2
trnWH1218128669-112181286691
trnCL1286135065012881351640
trnYL1351141565-51352141564-5
cox 1H141129461536TTATAA3141129461536TTATAA3
cox 2H29503633684ATGT_ _029503633684ATGTAA1
trnKH363537057103635370571-1
trnDH3706376964037053767630
atp 8H37703923154ATTT_ _-537683923156ATTTAA-7
atp 6H39194593675ATGTAA-139174591675ATGTAA-1
cox 3H45935376784ATGT_ _045915374784ATGT_ _0
trnGH5377544266053755438640
nad 3H54435787345ATATAA654395783345ATTTAA6
trnAH5794586370057905857680
trnRH5864592562-15858591962-1
trnNH5925599470-15919598870-1
trnS2H5994605966-15988605467-1
trnEH60596121636605461176411
trnFL6128619265061296192640
nad 5L619379061714TTGT_ _0619379061714TTGT_ _0
trnHL7907797064179077970641
nad 4L797293091338ATGTAG4797293151344ATTTAA1
nad 4LL93039602300ATGT_ _293179602286ATGT_ _2
trnTH9605966864196059669651
trnPL9670973566396719735653
nad 6H973910200462ATTTAA-1973910200462ATCTAA-1
cobH10200113281129ATGT_ _010200113281129ATGT_ _0
trnS1H113291139668241132911397699
nad1L1142112356936TTGTAA01140712339933TTGTAA0
trnL1L12354124196611234012405661
trnL2L12421124876701240712473670
rrnLL124881378713000124741376712940
trnVL13788138577001376813837700
rrnSL1385814673816013838146378000
CR1467415033360014638150063690

*Intergenic bp indicates gap nucleotides (positive value) or overlapped nucleotides (negative value) between two adjacent genes

 Table 3 

AT/CG skews in the mitochondrial protein coding genes (PCG), 2 rRNA genes, CR and the entire mitochondrial genome from three horseshoe crabs, Carcinoscorpius rotundicauda, Tachypleus tridentatus, and Limulus polyphemus

GeneAT-skewCG-skew
C. rotundicaudaT. tridentatusL. polyphemusC. rotundicaudaT. tridentatusL. polyphemus
cox1-0.065-0.0890.0440.1530.1460.210
cox2-0.002-0.0350.0360.3050.2810.331
cox3-0.130-0.141-0.0730.2300.1830.333
atp8-0.108-0.0880.0970.7500.5480.721
atp6-0.107-0.088-0.0290.4560.3980.526
nad1-0.243-0.200-0.388-0.438-0.414-0.511
nad2-0.140-0.143-0.0520.4960.4410.570
nad3-0.183-0.254-0.0420.4620.3410.524
nad4-0.257-0.251-0.374-0.432-0.393-0.485
nad4L-0.162-0.139-0.217-0.697-0.657-0.636
nad5-0.269-0.250-0.319-0.401-0.348-0.396
nad6-0.119-0.070-0.0060.6090.6120.629
cob-0.113-0.159-0.0810.3170.3060.427
rrnL0.0440.0300.1480.4560.4290.447
rrnS0.0880.1040.1450.3520.3330.331
CR-0.026-0.023-0.046-0.111-0.1150.077
13 PCG-0.162-0.163-0.1600.0460.0420.088
overall0.0360.0240.1110.3460.3150.399
 Table 4 

Codon usage pattern of the 13 mitochondrial protein-coding genes from three horseshoe crab species, Carcinoscorpius rotundicauda, Tachypleus tridentatus, and Limulus polyphemus

Amino acidsCodonNo.Amino acidsCodonNo.Amino acidsCodonNo.
CTLCTLCTL
NonpolarPheUUC466685SerAGA758366
AlaGCA484555UUU333310246AGC548
GCC271639TrpUGA10110294AGG2712
GCG496UGG10617AGU242122
GCU738476UCA11611498
IleAUC4342107PolarUCC493165
AUU346335241AsnAAU126139103UCG236
LeuCUA414294AAC352044UCU106122110
CUC171855CysUGC779Acidic
CUG214UGU434740AspGAC141124
CUU8580106GlnCAA596352GAU485136
UUA348363226CAG10315GluGAA686964
UUG443884GlyGGA110119109GAG141423
MetAUA215229171GGC81121Basic
AUG402343GGG322866ArgCGA323534
ProCCA606452GGU675941CGC337
CCC14827ThrACA787784CGG427
CCG224ACC293140CGU202014
CCU606869ACG102HisCAC191341
ValGUA766282ACU596553CAU566437
GUC6514TyrUAC262034LysAAA858466
GUG111936UAU10510887AAG16819
GUU738690

Bold numbers indicate strong differences (+/-25%) to Limulus polyphemus.

 Fig 2 

Putative tRNA secondary structures predicted from the 22 tRNA gene sequences found in the (A) Carcinoscorpius rotundicauda and (B) Tachypleus tridentatus mitochondrial genome.

Int J Biol Sci Image (Click on the image to enlarge.)
Int J Biol Sci Image (Click on the image to enlarge.)

Non-coding regions

The major non-coding regions (putative control region, CR) of 360 bp (C. rotundicauda) and 369 bp (T. tridentatus) were found between rrnS and trnQ (Table 2). The A + T contents in this region were higher, with 85% and 83.47% in C. rotundicauda and T. tridentatus, respectively, compared to other regions of the mitochondrial genomes (supplementary Table 1). The sequence of this region is highly conservative among xiphosuran species (Fig. 3A). The stem-loop structure found in C. rotundicauda (Fig. 3B) was composed of the 13 bp stem and 11 bp loop. However, the structure of T. tridentatus was rather complicated with the 10 bp stem and 13 bp loop followed by the alternative 8 bp stem and 12 bp loop (Fig. 3C). The mitochondrial genome of L. polyphemus had such an alternative structure [13], although it does not seem to be stable to build the complete stem-loop structure (Fig. 3D). A small (24 bp) non-coding fragment was found between the trnS (UGA) and nad1 gene in C. rotundicauda (Fig. 3E and Table 2. T. tridentatus and L. polyphemus had just a short (9 bp) and highly conserved non-coding fragment (TTTCTAAA) in this region.

 Fig 3 

Comparison in the primary and secondary structures of non-coding regions found in the mitochondrial genomes of Carcinoscorpius rotundicauda, Tachypleus tridentatus and Limulus polyphemus. The CR pair-wise alignment of the three xiphosuran species (A): the dots indicate nts identical to those of the first line in the alignment. Hairpin-like structure (60 bp and 57 bp) from the control region observed in C. rotundicauda (B), T. tridentatus (C), and L. polyphemus (D). An alignment of the gap region between trnS1 (UCN) and nad1 in three xiphosuran species (E): the box indicates a conserved region across three horseshoe crabs. L; L. polyphemus, C; C. rotundicauda, T; T. tridentatus.

Int J Biol Sci Image (Click on the image to enlarge.)

Phylogenetic analyses

The phylogenetic trees reconstructed from ML and BI algorithms were completely identical in clustering pattern and overall topology (Fig. 4). C. rotundicauda and T. tridentatus were recovered as a monophyly (BPML = 100, BPP = 1.00), whereas L. polyphemus solely formed a separate clade from two Asian species (BPML = 100, BPP = 1.00). Horseshoe crabs were placed as a sister relationship to the Arachnids by strong node-supporting values (Fig. 4A). The tree also supported the monophyletic clusterings of 'Opiliones + Scorpions' (BPML = 63, BPP = 0.83), 'Uropigy + Amblypygi' (BPML = 62, BPP = 1.0), Acari (Parasotoformes; BPML = 62, BPP = 1.00), 'Richinulei + Araneae' (BPML = 36, BPP = 0.84), Pantopoda (BPML = 60, BPP = 0.97) and 'Pseudoscorpiones + Acarifomes' (BPML = 94, BPP = 1.00; Fig. 4A). The monophyletic grouping of 'Opiliones + Scorpions' and 'Uropygi + Amblypygi' are compatible with the result from several previous analyses conducted based on morphological characters [23-25]. When the Asian horseshoe crabs were not considered in the tree, however, the tree was not resolved to support the monophyly of Arachnida, with Solifugae being independently placed at the basal root (Fig. 4B). In addition, such tree did not support the clustering of 'Opiliones + Scorpions' and 'Uropygi + Amblypygi' (Fig. 4B).

Discussion

The present study was designed to provide critical mitochondrial genome information of C. rotundicauda and T. tridentatus necessary for the concrete inference of phylogenetic relationships among species within Xiphosura and their placement within Chelicerata as well as Arthropoda. The overall architecture of both mitochondrial genomes is highly consistent with those previously reported in the Atlantic horseshoe crab, L. polyphemus. Those three mitochondrial genomes share the identical genome content (13 protein coding genes, 2 ribosomal RNA genes, 22 tRNA genes and a single large control region) and overall gene arrangement, as typically shown in many metazoan mitochondrial genomes [13, 26-28]. The mitochondrial genomes of C. rotundicauda and T. tridentatus are slightly larger than that of L. polyphemus, due to more frequent occurrence of gap nucleotides between genes.

In two Asian horseshoe crabs, ten of the protein coding genes (atp6, atp8, nad2, nad3, nad4. nad6, nad4L, cob, cox2 and cox3) have either Met or Ile as the start codon (ATA, ATG or ATT), while the remaining three (cox1, nad1, and nad5) start with Leu (TTA or TTG). Such start codon pattern was frequently found in a few other chelicerates [25, 26, 29, 30] as well as in several arthropods [31, 32]. Six (nad5, nad4L, cob, cox2, cox3 and atp8) protein coding genes have incomplete stop codon. Such truncated stop codons are commonly observed among many arthropod mitochondrial genomes [30-32] and are expected to be converted into a fully functional TAA via post-translational polyadenylation [29].

The sequences, anticodon nucleotides and secondary structures of tRNA genes are highly consistent with those previously reported in the Atlantic horseshoe crab. Most of the tRNA genes form dihydrouridine (DHU) arm. However, such structure is not found in trnS1 of T. tridentatus, and incomplete conformation with shortened stem of 2 bp is observed in C. rotundicauda. Missing of complete DHU arm in trnS1 has been reported among a number of metazoan mitochondrial genomes [29-32].

The A + T rich non-coding regions (putative control regions) were found in both two Asian horseshoe crab mitochondrial genomes. The sequence of this region is highly conservative among xiphosuran species as well as among a variety of arthropods [33]. However, the stem-loop structure is highly variable among three xiphosuran species. An interesting feature of stem-loop structure in T. tridentatus is its ability to form an alternative secondary structure that appears to be complete and stable. Although the mitochondrial genome of L. polyphemus has an additional 8 nucleotides complementary to the main stem, it does not appear to be stable. It is also surprising that there is a small non-coding fragment (24 bp) between trnS and nad1 in C. rotundicauda.

In the phylogenetic analysis among 43 chelicerates, C. rotundicauda and T. tridentatus were recovered as a monophyly, whereas L. polyphemus solely forms a separate clade from two Asian species. Such result supports previous findings showing clear genetic differentiation between L. polyphemus and Asian species [7, 8]. Major other chelicerate taxa were clustered in a single monophyletic assemblage, three xiphosuran species being placed at the basal root, suggesting that horseshoe crabs should be considered as an ancestral taxon in arthropods as well as in chelicerates. The phylogenetic tree without the information of Asian horseshoe crabs, by contrast, was not resolved to support monophyletic clustering of chelicerates, since solifuges formed a separate cluster. In conclusion, additional information of Xiphosuran mitochondrial genomes provide more robust and reliable perspective on the evolutionary history of chelicerate as well as arthropods.

 Fig 4 

Maximum likelihood phylogenetic trees of chelicerates based on the amino acid sequences of concatenated 12 mitochondrial protein coding genes. Tree was reconstructed either (A) with three horseshoe species (C. rotundicauda, T. tridentatus and L. polyphemus) or (B) with only L. polyphemus. Numbers at the branch indicate the percentages from ML bootstrapping (left) and Bayesian posterior probabilities (right).

Int J Biol Sci Image (Click on the image to enlarge.)
Int J Biol Sci Image (Click on the image to enlarge.)

Supplementary Material

Attachment Supplementary Table 1 

Acknowledgements

The present work was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science (NRF-2012R1A2A2A01046531 to U. W. H., NRF-2011-0025314 to S. H. R., and NRF-2011-0014704 to H. Y. S.).

Competing Interests

The authors have declared that no competing interest exists.

References

1. Sekiguchi K, Shuster CN. Limits on the global distribution of horseshoe crabs (Limulacea): lessons learned from two lifetimes of observations: Asia and America. In: (ed.) Tanacredi JT, Botton ML Smith DR. Biology and Conservation of Horseshoe Crabs. Dordrecht: Springer Publishing. 2009:5-24

2. Wheeler WC, Hayashi CY. The phylogeny of the extant chelicerate orders. Cladistics. 1998;14:173- 192

3. Masta SE, Longhorn SJ, Boore JL. Arachnid relationships based on mitochondrial genomes: asymmetric nucleotide and amino acid bias affects phylogenetic analyses. Mol Phylogenet Evol. 2009;50:117- 128

4. Shuster CN, Anderson LI. A history of skeletal structure: Clues to relationships among species. In: (ed.) Shuster CN, Barlow RB, Brockmann J. The American Horseshoe Crab. Cambridge, MA: Harvard University Press. 2003:154-188

5. Shishikura F, Nakamura S, Takahashi K, Sekiguchi K. Horseshoe crab phylogeny based on amino acid sequences of the fibrinopeptide-like peptide. J Exp Zool. 1982;223:89- 91

6. Sekiguchi K, Sugita H. Systematics and hybridization in the 4 living species of horseshoe crabs. Evolution. 1980;34:712- 718

7. Xia XH. Phylogenetic relationship among horseshoe crab species: effect of substitution models on phylogenetic analyses. Syst Biol. 2000;49:87- 100

8. Obst M, Faurby S, Bussarawit S, Funch P. Molecular phylogeny of extant horseshoe crabs (Xiphosura, Limulidae) indicates Paleogene diversification of Asian species. Mol Phylogenet Evol. 2012;62:21- 26

9. Shultz JW. A phylogenetic analysis of the arachnid orders based on morphological characters. Zool J Linn Soc. 2007;150:221- 265

10. Staton JL, Daehler LL, Brown WM. Mitochondrial gene arrangement of the horseshoe crab Limulus polyphemus L.: conservation of major features among arthropod classes. Mol Biol Evol. 1997;14:867- 874

11. Ovchinnikov S, Masta SE. Pseudoscorpion mitochondria show rearranged genes and genome-wide reductions of RNA gene sizes and inferred structures, yet typical nucleotide composition bias. BMC Evol Biol. 2012;12:31

12. Cheng S, Chang SY, Gravitt P. Long PCR. Nature. 1994;369:684- 685

13. Lavrov DV, Boore JL, Brown WN. The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus. Mol Biol Evol. 2000;17:813- 824

14. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL - X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876- 4882

15. Schattner P, Brooks AN, Lowe TM. The tRNAscan-SE: sonscan and snoGPS web servers for the detection of tRNAs and sonRNAs. Nucleic Acids Res. 2005;33(2):W686- 689

16. Arabi J, Judson ML, Deharvenq L, Lourenço WR, Cruaud C, Hassanin A. Nucleotide composition of CO1 sequences in Chelicerata (Arthropoda): detecting new mitogenomic rearrangements. J Mol Evol. 2012;74(1-2):81- 95

17. McWillian H, Li W, Squizzato S, Park YM, Buso N, Cowley AP, Lopex R. Analysis tool web services from the EMBL-EBI. Nucleic Acids Res. 2013;41:W597- 600

18. Castresana J. Selection of conserved block from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540- 552

19. Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21(9):2104- 2105

20. Abascal F, Posada D, Zardoya R. MtArt: a new model of amino acid replacement for Arthropoda. Mol Biol Evol. 2007;24(1):1- 5

21. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696- 704

22. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572- 1574

23. Giribet G, Edgecombe GD, Wheeler WC, Babbitt C. Phylogeny and systematic position of Opiliones: a combined analysis of chelicerate relationships using morphological and molecular data. Cladistics. 2002;18:5- 70

24. Shultz JW. A phylogenetic analysis of the arachnid orders based on morphological characters. Zool J Linn Soc. 2007;150:221- 265

25. Masta SE, McCall A, Longhorn SJ. Rare genomic changes and mitochondrial sequences provide independent support for congruent relationships among the sea spiders (Arthropoda, Pycnogonida). Mol Phylogenet Evol. 2010;57:59- 70

26. Masta SE. Mitochondrial rRNA secondary structures and genome arrangements distinguish chelicerates: comparisons with a harvestman (Arachnida: Opiliones: Phalangium opilio). Gene. 2010;449:1- 9

27. Fahrein K, Talarico G, Braband A, Podsiadlowski L. The complete mitochondrial genome of Pseudocellus pearsei (Chelicerata: Ricinulei) and a comparison of mitochondrial gene rearrangements in Arachnida. BMC Genomics. 2007;8:386

28. Castellana S, Vicario S, Saccone S. Evolutionary pattern of the mitochondrial genome in metazoa: Exploring the role of mutation and selection in mitochondrial protein-coding genes. Genome Biol. Evol. 2011;3:1067-1079

29. Ernsting BR, Edwards DD, Aldred KJ, Fites JS, Neff CR. Mitochondrial genome sequence of Unionicola foili (Acari: Unionicolidae): a unique gene order with implications for phylogenetic inference. Exp Appl Acarol. 2009;49:305- 16

30. Liu GH, Chen F, Chen YZ, Song HQ, Lin RQ, Zhou DH, Zhu XQ. Complete mitochondrial genome sequence data provides genetic evidence that the brown dog tick Rhipicephalus sanguineus (Acar: Ixodidae) represents a species complex. Int J Biol Sci. 2013;9(4):361- 369

31. Li H, Gao JY, Liu HY, Liu H, Liang AP, Zhou XG, Cai W. The architecture and complete sequence of mitochondrial genome of an assassin bug Agiosphodrus dohrni (Hemiptera: Reduviidae). Int J Biol Sci. 2011;7(6):792- 804

32. Wang YY, Liu XY, Yang D. The first mitochondrial genome for caddisfly (Insecta: Trichoptera) with phylogenetic implications. Int J Biol Sci. 2014;10(1):53- 63

33. Montagna M, Sassera D, Griggio F, Epis S, Bandi C, Gissi C. Tick-box for 3'-end formation of mitochondrial transcripts in Ixodida, basal chelicerates and Drosophila. PLoS ONE. 2012;7(10):e47538

34. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294- 299

35. Kambhampati S, Smith PT. PCR primers for the amplification of four insect mitochondrial gene fragments. Insect Mol Biol. 1995;4:233- 236

Author contact

Corresponding address Corresponding author: Ui Wook Hwang; Department of Biology, Teachers College, Kyungpook National University, Daegu 702-701, South Korea. E-mail: uwhwangac.kr; Tel: +82-53-950-5911; Fax: +82-53-950-6809; Cell: +82-10-4064-0375


Received 2014-2-1
Accepted 2014-4-2
Published 2014-4-16