Open-access Complete sequence and organization of the mitochondrial genome of Cyclemys atripons (Testudines, Geoemydidae)

Abstract

The Black Bridged Leaf Turtle, Cyclemys atripons (Testudines; Cryptodira; Geoemydidae), is a poorly known species within the genus Cyclemys. We determined the complete nucleotide sequence of the Cyclemys atripons mitochondrial genome (mtDNA) and found it to be 16,500 base pairs (bp) in length, with the genome organization, gene order and base composition being identical to that of the typical vertebrate. However, unlike for most turtle mtDNA so far reported, an extra base was not found in the NADH3 gene. The C. atripons control region of mtDNA was 981 bp long. Comparisons with three other geoemydids showed that the C. atripons control region contained a highly variable region at the 3' end composed of AT enriched tandem repeats containing a fifteen-unit 5'-A (AT)3-3' variable number of tandem repeats (VNTRs).

control region; Cyclemys atripons; mitochondrial genome; tandem repeats


EVOLUTIONARY GENETICS

SHORT COMMUNICATION

Complete sequence and organization of the mitochondrial genome of Cyclemys atripons (Testudines, Geoemydidae)

Li Zhang; Liu-wang Nie; Ying Zhang; Jin-long Rui; Yan-yun Zhang

Life Science College, Anhui Normal University, The Provincial Key Lab of the Conservation and Exploitation Research of Biological Resources in Anhui, Wuhu, China

Send correspondence to Send correspondence to: Prof. Liu-wang Nie Life Science College, Anhui Normal University 1 East Beijing Road, Wuhu 241000, Anhui, China E-mail: lwnie@mail.ahnu.edu.cn

ABSTRACT

The Black Bridged Leaf Turtle, Cyclemys atripons (Testudines; Cryptodira; Geoemydidae), is a poorly known species within the genus Cyclemys. We determined the complete nucleotide sequence of the Cyclemys atripons mitochondrial genome (mtDNA) and found it to be 16,500 base pairs (bp) in length, with the genome organization, gene order and base composition being identical to that of the typical vertebrate. However, unlike for most turtle mtDNA so far reported, an extra base was not found in the NADH3 gene. The C. atripons control region of mtDNA was 981 bp long. Comparisons with three other geoemydids showed that the C. atripons control region contained a highly variable region at the 3' end composed of AT enriched tandem repeats containing a fifteen-unit 5'-A (AT)3-3' variable number of tandem repeats (VNTRs).

Key words: control region, Cyclemys atripons, mitochondrial genome, tandem repeats.

Vertebrate mitochondrial (mt) DNA forms a double-stranded circular molecule of about 15-20 kb which generally contain 37 genes encoding 13 proteins, 22 tRNAs, 2 rRNAs and a major noncoding region bearing signals for mitochondrial replication and transcription (Wolstenholme, 1992). Due to its maternal inheritance and relative lack of recombination, the mitochondrial genome has been widely employed as a marker in vertebrate phylogenetic analyses, and have been often used in turtle science. Turtles are easily recognizable by the public, with approximately 270 species in the world (Iverson, 1992) and have been widely studied, with many earlier studies of mtDNA having concentrated on polymorphism analysis using restriction fragment length polymorphism (RFLP) and determining partial sequences (Lamb et al., 1989; 1994). However, the trend in more and more studies is to move to direct sequencing of the complete mtDNA (Peng et al., 2005, 2006; Parham et al., 2006). As of March 2007, complete mitochondrial genomes have been released from GenBank for only 17 turtle species, including 16 cryptodiran turtles and one side-necked turtle, which is far from being sufficient for reliable turtle studies.

The Black Bridged Leaf Turtle, Cyclemys atripons (Testudines; Cryptodira; Geoemydidae) is a poorly understood cryptodiran turtle species within the genus Cyclemys (Guicking et al., 2002). Previously, only the cytochrome b (Cyt b) gene of C. atripons mtDNA has been published (Spinks et al., 2004), clearly, further studies on this species are necessary. In our study, described in this paper, we sequenced and characterized the complete mitochondrial genome of C. atripons, which has laid the foundation for the further comparative analyses between C. atripons and other turtles.

In 2005 a C. atripons specimen was obtained from the suburbia of Longzhou city in the Chinese region of Guangxi, after natural death frozen at -80 ºC for preservation. Total DNA was extracted from the liver and muscle tissue using the proteinase K method (Sambrook and Russell, 2001) and kept at -20 ºC until needed for polymerase chain reaction (PCR) amplification.

Based on partial sequences reported by Spinks et al. (2004) and the similarity of mtDNA sequences of the painted turtle (Chrysemys picta; GenBank NC_002073) and Reeve's Turtle (Chinemys reevesii; GenBank AY676201) we designed 16 pairs of primers for PCR amplification (Table 1). The PCR was carried out in a total volume of 25 µL containing 100 ng of sample genomic DNA, 2.5 µL of 10×Buffer (TaKaRa, Japan), 2 µL of 2.5 mol L-1 of MgCl2 , 1.5 µL of each dNTP, 0.25 µL of each primer (25 µmol L-1) and 1 unit of Taq DNA polymerase (TaKaRa). The thermal cycles were 95 ºC pre-denaturing for 2 min, followed by 35 cycles of 94 ºC for 40 s, 51 ºC to 58 ºC for 45 s and 72 ºC for 1 min, plus a final extension at 72 ºC for 10 min. The resultant PCR fragments were first resolved on 1% (w/v) agarose gel (Promega, USA). After electrophoreses, gels were stained with ethidium bromide and bands were visualized under ultraviolet. Bands of intended size were excised and recovered with Gel Extract Purification Kit (TaKaRa, Japan). The cleaned PCR products were sequenced in both directions on an ABI3730 automated sequencer (Invitrogen Biotechnology). The sequences obtained from each sequencing reaction averaged 1000 bp in length and each segment overlapped the next contig by roughly 150 bp. The whole mtDNA genome sequence was read at least twice.

Sequence data were analyzed with the EditSeq (DNASTAR) and ClustalX1.8 (Thompson et al., 1997) programs. The locations of protein-coding, rRNA and tRNA genes were identified by the tRNA Scan-SE1.21 and SQUEIN v. 5.35 programs, which were also used for the comparisons with the corresponding sequences from the other turtles cited above. The analysis of the control region sequence was carried out with the DNAsis and BioEdit programs. The resultant complete mitochondrial genome of C. atripons (16,500 bp) was deposited in GenBank under accession number EF067858.

The structural organization (Table 2) and gene order (Figure 1) of the complete C. atripons mitochondrial genome was identical that of other typical vertebrates, with the genome containing the following: 13 protein-coding genes, all of which except NADH6 being encoded on the H-strand; 22 tRNA genes, 14 on the H-strand and 8 on the L-strand; 2 rRNA genes, 12S and 16S, both on the H-strand; and one control region. There were few, or small, noncoding intergenic spacer nucleotides, with intervening sequences of 8 bp between cytochrome c oxidase mitochondrial subunit II gene (COII) and tRNALys plus 13 bp between NADH4 and tRNAHis (Table 2). We found that the base composition of the major coding strand of C. atripons mtDNA was A = 34.42%, G = 13.01%, C = 25.36% and T = 27.20%, demonstrating the low G and high A+T bias seen in most other turtles (Pu et al., 2005; Peng et al., 2005, 2006). We also found that in C. atripons mtDNA three protein genes (COIII, NADH6 and Cyt b) have an incomplete stop codon such as T, while the cytochrome c oxidase mitochondrial subunit I gene (COI) has GTG instead of ATG as a start codon. As in other vertebrate mitochondrial genomes, we found three instances of reading frame overlap, 10 nucleotides for ATP8 and ATP6, 7 for NADH4L and NADH4, and 5 for NADH5 and NADH6 (Table 2).


However, in C. atripons our analysis did not find the extra base usually found at a specific position in NADH3 of most other turtles (Mindell et al., 1998, 1999; Pu et al., 2005; Parham et al., 2006). Such an insertion in the NADH3 gene has been reported in most turtles with the exception of Pelodiscus sinensis and Kinosternon flavescens (Peng et al., 2005; Pu et al., 2005; Parham et al., 2006). It is generally thought the base could be related to a TAA stop codon frameshift prematurely terminating protein translation if not corrected by RNA editing or other mechanisms (Mindell et al., 1998, 1999). However, since the additional base is apparently absent from C. atripons more studies are needed to ascertain whether or not the extra base in NADH3 is a common characteristic of turtles or is specific to certain species and genera.

We found that the C. atripons mitochondrial genome contained 22 tRNA genes, ranging in size from 66 nucleotides to 76 nucleotides (Table 2), interspersed between the rRNA and protein-coding genes, which is typical of the mtDNAs of other vertebrates. Most of the tRNA genes could be folded into the canonical cloverleaf secondary structure, the exception being the tRNASer (AGY) gene which lacks the dihydrouridine arm (D arm). The length of the C. atripons 12S rRNA gene was 969 nucleotides and the 16S rRNA gene was 1,601 nucleotides long, these genes, as in other vertebrates, being separated by the tRNAVal gene and positioned between the of tRNAPhe and tRNALeu (UUR) genes (Table 2).

In the C. atripons mitochondrial DNA, we also found that the light-strand replication origin (31 nucleotides) was located between the tRNAAsn and tRNACys genes inside the WANCY tRNA gene cluster (Figure 1). This region has also been discovered in mtDNAs of all other cryptodiran turtles investigated (Pu et al., 2005; Peng et al., 2006), contrasting with its apparent disappearance from Pelomedusa subrufa (Zardoya and Meyer, 1998b). This C. atripons sequence may potentially fold into a stable stem-loop secondary structure with a stem comprised of 10 bp and a loop of 10 bp. The secondary structures of the origin of light strand replication (OL) for 19 cryptodiran turtles (Figure 2) shows that the OL sequence nucleotides are rather conserved and the secondary structures of these sequences are also similar because 9 bp are identical in the stems, possibly a common characteristic of cryptodiran turtles.


We found that the C. atripons D-loop control region was 981 bp long, 69.52% A+T rich and flanked by tRNAPro and tRNAPhe genes (Table 2). A comparison of the complete control region sequences of four geoemydid turtles is given in the online edition of this paper (Figure S1), Similar to three other geoemydid turtles, three conserved sequence blocks (CSBs) 1-3 (Walberg and Clayton, 1981) were identified in the C. atripons control region. The whole lengths of four control region sequences ranged from 981 bp in C. atripons to 1,379 bp in Cuora aurocapitata, mainly resulting from sequences positioned at the 3' end. Interestingly, a large number of AT enriched tandem repeats containing variable number tandem repeats (VNTRs) were revealed at the 3' end (right domain) of the control regions. Furthermore, the composition and number of these tandem units were different for the different species, with the C. atripons VNTR being composed of fifteen 5'-A (AT) 3-3' units (Figure S1).

The control region is usually considered to be the most variable parts of mtDNAs in terms of nucleotide substitutions, short insertions/deletions and VNTRs dynamics. However, these variations are not distributed randomly across the whole region but occur in particular hypervariable sites and domains at the 5' and 3' ends (Su, 2005). Previous studies of turtles utilizing control region sequences were primarily focused on the 5' end adjacent to tRNAPro and several regulatory motifs (Lamb et al., 1994; Walker et al., 1997; Walker and Avise, 1998). However, at present most work focuses on the 3' end close to tRNAPhe, especially tandemly repeated sequences, including VNTRs (Serb et al., 2001). The length difference between mitochondrial genomes among species is caused mainly by the divergent tandem repeats, which are thought to be generated by strand slippage and mispairing during replication (Fumagalli et al., 1996).

Tandemly repeated control region DNA has been reported from an ever-growing number of taxa. Zardoya and Meyer (1998a) characterized six tandem repeats (containing VNTRs) in the 3' domain of the P. subrufa control region and suggested that this sequence might be a potentially informative molecular marker for population studies by its unique localization in the maternally inherited mitochondrial molecule. What is remarkable is that the tandem repeats are present in the four geoemydid turtles discussed in our present paper. The repeat consists of two different repeat cores, "ATTATATC" followed by "AT" in Pyxidea mouhtii (DQ659152) than the one "A (AT) 3" in C. atripons; in C. reevesii (AY676201) it is over ten 5'-ATATATC-3' units succeeded by AT-rich repeat; whereas in C. aurocapitata (AY874540) an approximately 490 bp AT-rich repeat is in the 3' of the CR , with only a few "G" nucleotides and no "C" nucleotide.

In the mitochondrial genomes of C. atripons and other turtles, AT enriched tandem repeats (containing VNTRs) reflect heteroplasmy, suggesting interspecies turtle genetic diversity. The occurrence of tandemly repeated mtDNA in the control region could be regarded as a special molecular marker in turtle species researches. However, to confirm this issue and the presence or absence of variation between specimens of the same species (intraspecies polymorphism) more informative characters need to be obtained from more turtle species and specimens. Taken as a whole, our study highlights the need for further work focusing on the 3' domain of the mitochondrial control region.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (grant number 30770296), the Natural and Science Key Project of Anhui Education Department (grant number KJ2007A022) and the Key Laboratory of Biotic Environment and Ecological Safety of Anhui province. We thank three anonymous reviewers and the technical editor for critically reviewing the manuscript.

Internet Resources

Received: March 19, 2007; Accepted: November 13, 2007.

Associate Editor: João S. Morgante

Supplementary Material

The following online material is available for this article:

- Figure S1. The complete control regions of four geoemydids. Dots and dashes represent the same nucleotides as Pyxidea and gaps, respectively. The potential regulatory elements (CSBs) 1-3 are boxed. Shaded region designates tandem repeat sequences at the 3' end of CRs.


This material is made available as part of the online article from http://www.scielo.br.gmb.

References

  • Fumagalli L, Taberlet P, Favre L and Hausser J (1996) Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews. Mol Biol Evol 13:31-46.
  • Guicking D, Fritz U, Wink M and Lehr E (2002) New data on the diversity of the Southeast Asian leaf turtle genus Cyclemys Bell, 1834. Molecular results (Reptilia, Testudines, Geoemydidae*). Faunist Abh Staatl Mus Tierkunde Dresden 23:75-86.
  • Iverson JB (1992) A Revised Checklist with Distribution Maps of the Turtles of the World. Privately printed, Richmond, 363 pp.
  • Lamb T, Lydeard C and Gibbons JW (1994) Molecular systematics of map turtles (Graptemys), a comparison of mitochondrial restriction site vs. sequence data. Syst Biol 43:543-559.
  • Lamb T, Avise JC and Gibbons JW (1989) Phylo-geographical patterns in mitochondrial DNA of the desert tortoise (Xerobates agassizi), and evolution relationships among the North American gother tortoise. Evolution 43:76-87.
  • Mindell DP, Sorenson MD and Dimcheff DE (1998) An extra nucleotide is not translated in mitochondrial ND3 of some birds and turtles. Mol Biol Evol 15:1568-1571.
  • Mindell DP, Sorenson MD, Dimcheff DE, Hasegawa M, Ast JC and Yuri T (1999) Interordinal relationships of birds and other reptiles based on whole mitochondrial genomes. Syst Biol 48:138-152.
  • Parham JF, Feldman CR and Boore J (2006) The complete mitochondrial genome of the enigmatic bigheaded turtle (Platysternon), description of unusual genomic features and the reconciliation of phylogenetic hypotheses based on mitochondrial and nuclear DNA. BMC Evol Biol 6:11.
  • Parham JF, Macey JR, Papenfuss TJ, Feldman CR, Türkzan O, Polymeni R and Boore J (2006) The phylogeny of Mediterranean tortoises and their close relatives based on complete mitochondrial genome sequences from museum specimens. Mol Phylogenet Evol 38:50-64.
  • Peng QL, Nie LW and Pu YG (2006) Complete mitochondrial genome of Chinese big-headed turtle, Platysternon megacephalum, with a novel gene organization in vertebrate mtDNA. Gene 380:14-20.
  • Peng QL, Pu YG, Wang ZF and Nie LW (2005) Complete mitochondrial DNA sequence of Chinese softshell turtle (Pelodiscus sinensis). Chin J Biochem Mol Biol 21:591-596.
  • Pu YG, Peng QL, Wang ZF and Nie LW (2005) Sequence and organization of the complete mitochondrial genome of the Chinese three-keeled pond turtle (Chinemys reevesii). Acta Zool Sin 51:691-696.
  • Sambrook J and Russell DW (2001) Molecular Cloning, A Laboratory Manual. 3rd ed. Cold Spring Harbor Laboratory Press, New York.
  • Serb JM, Phillips CA and Iverson JB (2001) Molecular phylogeny and biogeography of Kinosternon flavescens based on complete mitochondrial control region sequences. Mol Phylogenet Evol 18:149-162.
  • Spinks PQ, Shaffer HB, Iverson JB and McCord WP (2004) Phylogenetic hypotheses for the turtle family Geoemydidae. Mol Phylogenet Evol 32:164-182.
  • Su Y (2005) Conserved and heteroplasmy on mitochondrial DNA control region in animal. Sichuan J Zool 24:669-672.
  • Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F and Higgins DG (1997) The clustalx windows interface: Exible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 24:4876-4882.
  • Walker D, Nelson WS, Buhlmann KA and Avise JC (1997) Mitochondrial DNA phylogeography and subspecies issues in the monotypic freshwater turtle (Sternotherus odoratus). Copeia 1997:16-21.
  • Walker D and Avise JC (1998) Principles of phylogeography as illustrated by freshwater and terrestrial turtles in the southeastern United States. Annu Rev Ecol Syst 29:23-58.
  • Walberg MW and Clayton DA (1981) Sequence and properties of the human KB cell and mouse L cell D-loop regions of mitochondrial DNA. Nucleic Acids Res 9:5411-5421.
  • Wolstenholme DR (1992) Animal mitochondrial DNA: Structure cursor for mitochondrial genomic rearrangement. Int Rev Cytol 141:173-216.
  • Zardoya R and Meyer A (1998a) Cloning and characterization of a microsatellite in the mitochondrial control region of the African side-necked turtle, Pelomedusa subrufa. Gene 216:149-153.
  • Zardoya R and Meyer A (1998b) Complete mitochondrial genome suggests diapsid affinities of turtles. Proc Natl Acad Sci USA 95:14226-14231.
  • Scan-SE1.21 program is available at http://lowelab.ucsc.edu/ tRNA
  • EditSeq, DNAsis and BioEdit programs are available at http:// www.bioon.com/Soft/
  • Send correspondence to:
    Prof. Liu-wang Nie
    Life Science College, Anhui Normal University
    1 East Beijing Road, Wuhu
    241000, Anhui, China
    E-mail:
  • Publication Dates

    • Publication in this collection
      19 Aug 2008
    • Date of issue
      2008

    History

    • Received
      19 Mar 2007
    • Accepted
      13 Nov 2007
    location_on
    Sociedade Brasileira de Genética Rua Cap. Adelmio Norberto da Silva, 736, 14025-670 Ribeirão Preto SP Brazil, Tel.: (55 16) 3911-4130 / Fax.: (55 16) 3621-3552 - Ribeirão Preto - SP - Brazil
    E-mail: editor@gmb.org.br
    rss_feed Stay informed of issues for this journal through your RSS reader
    Acessibilidade / Reportar erro