Complete Nucleotide Sequence and Gene Rearrangement of the Mitochondrial Genome of the Bell-Ring Frog, Buergeria Buergeri (Family Rhacophoridae)
Total Page:16
File Type:pdf, Size:1020Kb
Genes Genet. Syst. (2004) 79, p. 151–163 Complete nucleotide sequence and gene rearrangement of the mitochondrial genome of the bell-ring frog, Buergeria buergeri (family Rhacophoridae) Naomi Sano1, Atsushi Kurabayashi1, Tamotsu Fujii2, Hiromichi Yonekawa3, and Masayuki Sumida1* 1Institute for Amphibian Biology, Hiroshima University, Higashihiroshima, Hiroshima 739-8526, Japan 2Department of Health Science, Hiroshima Prefectural Women’s University, Hiroshima 734-8558, Japan 3Department of Laboratory Animal Science, The Tokyo Metropolitan Institute of Medical Science, Tokyo 113-8613, Japan (Received 22 March 2004, accepted 26 May 2004) In this study we determined the complete nucleotide sequence (19,959 bp) of the mitochondrial DNA of the rhacophorid frog Buergeria buergeri. The gene content, nucleotide composition, and codon usage of B. buergeri conformed to those of typ- ical vertebrate patterns. However, due to an accumulation of lengthy repetitive sequences in the D-loop region, this species possesses the largest mitochondrial genome among all the vertebrates examined so far. Comparison of the gene organ- izations among amphibian species (Rana, Xenopus, salamanders and caecilians) revealed that the positioning of four tRNA genes and the ND5 gene in the mtDNA of B. buergeri diverged from the common vertebrate gene arrangement shared by Xenopus, salamanders and caecilians. The unique positions of the tRNA genes in B. buergeri are shared by ranid frogs, indicating that the rearrangements of the tRNA genes occurred in a common ancestral lineage of ranids and rhacophorids. On the other hand, the novel position of the ND5 gene seems to have arisen in a lineage leading to rhacophorids (and other closely related taxa) after ranid divergence. Phylogenetic analysis based on nucleotide sequence data of all mito- chondrial genes also supported the gene rearrangement pathway. Key words: Buergeria buergeri, complete mitochondrial genome, gene rearrange- ment, Japanese bell-ring frog, ND5 gene 1981; Roe et al., 1985; Tzeng et al., 1992; Boore, 1999). INTRODUCTION Although minor rearrangements have been reported for Vertebrate mitochondrial DNA (mtDNA) is a closed cir- marsupials (Pääbo et al., 1991), birds (Desjardins and cular molecule. The genome organization is highly com- Morais, 1990, 1991; Quinn and Wilson, 1993; Mindell et pact and the genome is approximately 16 kbp in length al., 1998), reptiles (Kumazawa and Nishida, 1995; Quinn (Wolstenholme, 1992). This genome typically contains and Mindell, 1996; Macey et al., 1997), lampreys (Lee and 37 genes for 2 ribosomal (r)RNAs, 22 transfer (t)RNAs Kocher, 1995), and teleost fishes (Miya and Nishida, and 13 proteins, and a long noncoding region called the 1999; Inoue et al., 2001; Miya et al., 2001), most of these D-loop that includes signals required for regulating and cases involve only a few rearrangements of tRNA genes initiating mtDNA replication and transcription (Wolsten- and/or the D-loop region. holme, 1992; Janke et al., 1994). Mitochondrial gene arrangements have attracted the Mitochondrial gene arrangements are generally con- attention of evolutionary biologists as novel phylogenetic served in vertebrates. All 37 genes are arranged in the markers (Smith et al., 1993; Kumazawa and Nishida, same relative order in almost all vertebrate species from 1995; Quinn and Mindell, 1996; Macey et al., 1997; Boore teleost fishes to eutherian mammals (Anderson et al., and Brown, 1998; Boore, 1999; Kurabayashi and Ueshima, 2000). The complete mtDNA sequences have Edited by Toshihiko Shiroishi been published for only six amphibian species, including * Corresponding author. E-mail: [email protected] a caecilian (Zardoya and Meyer, 2000), three sala- 152 N. SANO et al. manders, (Zardoya et al., 2003; Zhang et al., 2003a, b), ity to the family Ranidae (true frogs), and the genus Buer- and two anurans, the clawed frog Xenopus laevis (Roe et geria is regarded as the most basal group in the former al., 1985) and the Japanese pond frog Rana nigromacu- (Channing, 1989; Jiang et al., 1987; Liem, 1970; Richards lata (Sumida et al., 2001). The mitochondrial gene and Moore, 1998; Wilkinson and Drewes, 2000; Wilkinson arrangements of five of these six amphibian species are et al., 2002). identical with those of typical vertebrates. However, in In this report, we present the first data on the complete the R. nigromaculata mtDNA, the positions of four tRNA mtDNA sequence of a rhacophorid frog and describe fea- genes (tRNA-Leu(CUN), tRNA-Thr, tRNA-Pro, and tures of the genome. The evolutionary implications of tRNA-Phe) differ from those of typical vertebrates (see our findings are also discussed. Fig. 3). The same gene arrangement is also found in other ranid frogs so far investigated (Rana catesbeiana, MATERIALS AND METHODS Yoneyama, 1987; Rana limnocharis, Macey et al., 1997; Rana porosa, Sumida et al., 2000). A lack of available DNA sources. Bell-ring frogs (Buergeria buergeri) were information on the mitochondrial genomic organization of collected from Ota River, Hiroshima prefecture, western other anuran family members makes it difficult to deter- Japan. The total genomic DNAs were extracted from a mine whether this unique arrangement is shared by other clipped toe of a living frog using the DNeasy Tissue Kit anuran groups besides the ranids. (QUIAGEN) according to the manufacturer’s protocol. In order to elucidate various aspects of mitochondrial gene rearrangement, we determined the complete PCR and sequencing. The total length of B. buergeri mtDNA sequence of the bell-ring frog Buergeria buergeri, mtDNA was amplified by polymerase chain reaction a representative of the family Rhacophoridae (tree frogs). (PCR), beginning with partial mitochondrial segments Rhacophoridae is generally considered to have close affin- and finishing with the remaining mtDNA region. Two Fig. 1. Sequencing and cloning strategy for Buergeria buergeri mtDNA. Localizations and directions of primers used in the LA-PCR amplification and DNA sequencing are denoted by arrowheads. The sequences of these primers are available from the www site, http://home.hiroshimau.ac.jp/~amphibia/syukeisei/usedprimers.html. LA-PCR products are shown as bold lines below the gene map. Cloned restriction fragments and their lengths are indicated above the gene map. Table 1. Primers used for PCR amplification in the present study. Primer name Sequence (5’-3’) F20N7a TGAATCGGGGGCCAACCAG R16b ATAGTGGGGTATCTAATCCCAGTTTGTTTT FR16c AAAACAAACTGGGATTAGATACCCCACTAT R51d GGTCTGAACTCAGATCACGTA F70e GGGTATCCCAGTGGTGCAGCCGCTACTAAT R71e CGAATGTCTTGTTCGTCATTGAGGTTATGA ND5Fow1e TTYATHGGHTGRGARGGVGTNGG R16M1c GGGTATCTAATCCCAGTTTG Positions with mixed bases are labeled with their IUB codes: H = A / T / G, N = A / T / G / C, R = A / G, V = A / G / C, Y = C / T. aDesigned based on Sumida et al. (2001), Roe et al. (1985), Zardoya and Meyer (2000) and Zardoya and Meyer (1997). bSumida et al. (1998). cModified from Sumida et al. (1998). dSumida et al. (2002). ePresent study. Mitochondrial genome of rhacophorid frog 153 partial mitochondrial segments were amplified from B. employed to sequence almost all portions of the mtDNA buergeri total DNA by long-and-accurate PCR (LA-PCR) (Fig. 1). The sequencing was performed using 373A and using two primer sets (F20N7 and R16; FR16 and R51) 3100-Avant automated DNA sequencers (ABI) with DYE- (Fig. 1 and Table 1). PCR mixtures were prepared using namic ET Terminator cycle sequencing reagent a TaKaRa LA Taq™ Kit as recommended by the manu- (Amersham). Sequencing primers for internal portions facturer (TaKaRa). LA-PCR reactions consisted of an of the long PCR fragment were designed by cloning five initial denaturation at 94°C for 1 min, 30 cycles of dena- restriction fragments (see Fig. 1) into pUC 118 E. coli vec- turation at 94°C for 20 sec plus annealing and extension tor and sequencing them. The D-loop region was diffi- at 60°C for 3 min, and a final extension at 72°C for 10 cult to sequence by primer walking due to an abundance min. The resultant PCR fragments were electrophore- of lengthy tandem repeat units. To determine the pre- sed on a 0.7% agarose gel, and the DNAs were purified cise sequence of this region, a series of deleted subclones from excised pieces of gel using Wizard SV Gel and the was made from the clones of an Sma I/EcoR I-digested PCR Clean-UP System (Promega) and used for sequenc- fragment (see Fig. 1) using the Exonuclease III deletion ing (see below). The remaining mtDNA fragment was method (Henikoff, 1987). The nucleotide sequence of the amplified with the primer set F70 and R71 (Fig. 1 and B. buergeri mtDNA was deposited in the DDBJ database Table 1) designed based on the ND5 gene and 16S rRNA accession number AB127977. gene sequences determined above. LA-PCR reactions consisted of 1 cycle of 1 min at 94°C, 14 cycles of 20 sec Phylogenetic analysis. For the phylogenetic analysis, at 98°C followed by 20 min at 68°C, 17 cycles of 20 sec at we created an alignment dataset using CLUSTAL W 98°C followed by 20 min 20 sec at 68°C, and 1 final cycle (Thompson et al., 1994). The data set consisted of all 37 of 10 min at 72°C. The amplified mtDNA fragment of mitochondrial gene sequences from 9 vertebrates, includ- approximately 16.5 kbp was then purified using ing 7 amphibians, 1 coelacanth (Latimeria chalumnae), MicroSpin™ S-300 HR Columns (Pharmacia Biotech) and and 1 lungfish (Protopterus dolloi). The latter two were used for the sequencing reaction. used as outgroups. The alignment was checked by eye, The entire mtDNA genome of B. buergeri was sequenced and all positions with gaps and ambiguous sites were from both strands. The primer walking method was excluded. Based on the alignment data (12,865 nucle- Fig. 2. Organization of the B. buergeri mitochondrial genome. All protein-coding genes are encoded by the H strand, with the exception of ND6, which is encoded by the L strand. Transfer RNA genes are represented by the standard one-letter amino acid code and those encoded by the H and L strands are shown outside and inside the circle, respectively.