1. College of Life Science, China West Normal University, Nanchong 637002, China
2. Sichuan Agricultural University, Yaan 625014, China
3. Nanchong Agricultural Science Research Institute, Nanchong 637000, China
4. Sichuan Province Traditional Chinese Medicinal Materials Company Dujiangyan Raising Deer Field, Dujiangyan 611800, China
We obtained the complete mitochondrial genome of U.thibetanus mupinensis by DNA sequencing based on the PCR fragments of 18 primers we designed. The results indicate that the mtDNA is 16 868 bp in size, encodes 13 protein genes, 22 tRNA genes, and 2 rRNA genes, with an overall H-strand base composition of 31.2% A, 25.4% C, 15.5% G and 27.9% T. The sequence of the control region (CR) located between tRNA-Pro and tRNA-Phe is 1422 bp in size, consists of 8.43% of the whole genome, GC content is 51.9% and has a 6bp tandem repeat and two 10bp tandem repeats identified by using the Tandem Repeats Finder. U. thibetanus mupinensis mitochondrial genome shares high similarity with those of three other Ursidae: U. americanus (91.46%), U. arctos (89.25%) and U. maritimus (87.66%).
Keywords: Ursus thibetanus mupinensis, mitochondrial genomes, sequencing, sequence analysis
Mitochondria are vital subcellular organelles, responsible for the oxidizing reaction of the tricarboxylic acid cycle, the electron transfer and the energy metabolism in cells and have an independent genetic material called the mitochondria genome (mtDNA). Mammalian mtGenomes exist as closed circular strands and have a set of 13 protein-coding genes (NADH-ubiquinone oxidoreductase chain 1, 2, 3, 4L, 4, 5, 6, cytochrome c oxidase subunit 1, 2 and 3, ATP synthase F0 subunit 8 and ATP synthase F0 subunit 8), two rRNA genes (12s RNA and 16s RNA), and a full set of 22 tRNAgenes. The two strands that make up the genome are commonly known as the heavy strand (H-strand) and the light strand (L-strand) because of molecular weight differences as a result of different base compositions. Of the 13 protein-coding genes, 12 are on the H-strand and only one is on the L-strand. Noncoding regions are mainly limited to areas called the D-loop, thought to have functional roles in replication and transcription, and origin of replication of the L-strand (OL, thought to have a functional role in replication) . The gene order is also highly conserved among most vertebrates . During the last decades, mitochondrial genome sequence and gene arrangement comparisons were employed as powerful new tools for resolving ancient phylogenetic relationships [3, 4-8].
The black bear (Ursus thibetanus mupinensis) is listed in “The Convention on International Trade in Endangered Species of Wild Fauna and Flora “(CITES) the appendix I species, the national 2 levels of key protections wild animals and “China Red Data Book of Endangered Animals” V species . There have been researches for mitochondria genome of Ursidae. Delisle and Strobeck obtained mitochondrial genome sequences of the Americas black bear (Ursus americanus), polar bear (Ursus maritimus) and brown bear (Ursus arctos) in 2002 . Despite some studies of Asian black bear mitochondria individual genes such as the Cytb gene, there are no reports about the mitochondrial genome of the black bear, which live within the boundaries of Sichuan. In the current study, we report the complete mitochondrial genome sequence from a single Asian black bear Sichuan subspecies Ursus thibetanus mupinensis.
Skeletal muscle tissue samples of the Asian black bear Sichuan subspecies (Ursus thibetanus mupinensis) were obtained from black bears collected from Sichuan Province Traditional Chinese Medicinal Materials Company Dujiangyan Raising Deer Field.
Tissue sample was preserved in 70% ethanol until analysis. Total DNA was extracted from approximately 3×3×3 mm of the tissue using the conventional proteinase K/phenol/chloroform method  with some modifications as described by Masuda and Yoshida .
A series of primers were designed (Table 1), based on conserved regions identified from an alignment of published complete mitochondrial genomes from carnivores by Primer Premier 5.0: Ursus americanus (American black bear) (GenBank accession No. AF303109), Ursus maritimus (polar bear) (GenBank accession No. AF303111) and Ursus arctos (brown bear) (GenBank accession No. AF303110) . PCR was using rTaq DNA polymerase (TAKARA) with 1–10μl DNA extracts in a total volume of 50μl.
According to PCR reaction condition exists difference because G, C and A, T basis Content and the base number of primers is different, we lead primers to be divided into 4 sets, the reaction condition of the same of set consistent. And established temperature steps a degree circulation in each set. It set circumstance with the reaction condition for corresponding describe respectively in table 2.
Primers designed with Primer Premier 5.0 for amplifying the complete mitochondrial genome of U.thibetanus mupinensis
|No.||Forward Primer||Reverse Primer||Product length (bp)||Forward Primer Sequence (5'-3')||Reverse Primer Sequence (5'-3')|
Reaction condition of PCR for 18 pairs primerse (Group 1 primers to have: 2, 5, 7, 9, 18; Group 2 primers to have: 1, 3, 4, 6, 8, 10,12,15; Group 3 primers to have: 11,14,16; Group 4 primers to have: 13,17)
|Reaction condition||Group 1||Group 2||Group 3||Group 4|
|Primer extension temperature||72℃||72℃||72℃||72℃|
|Cycleindex for each annealing temperature||2||2||2||2|
|Primer extension temperature||72℃||72℃||72℃||72℃|
The PCR products were fractionated on 1.5% (W/V) agarose gel, and selected bands were purified using a gel extraction kit (Sangon, Shanghai, China). The purified PCR products were ligated into the pMD18-T vector (TakaRa) and transformed into JM109 competent cells. Bacteria were grown in LB-ampicillin agar. Cloned PCR products were sequenced by BGI Life Tech Co., Ltd.
The gene sequences of U. thibetanus mupinensis mitochondrial genome were identified by sequence comparison with published Carnivora gene sequences and similarity analysis using DNASTAR. Start and stop codons were used to help defining the sequences of protein coding genes. Percentage of GC base component and Base skew overall were calculated to analyze the genome characters. Tandem repeat was found by Tandem Repeats Finder .
The complete mtDNA sequence has been determined for U. thibetanus mupinensis mtDNA (GenBank accession no. DQ402478). The mtDNA is 16 868 bp in size, shorter than those of U. maritimus and U. arctos, which are 17017bp and 17020bp in length, yet longer than that of U. americanus which is 16841bp in length. The size differences result from different lengths of the control region among these species. It encodes 13 protein genes, 22 tRNA genes, and 2 rRNA genes (Fig 1). Eight tRNA genes and one protein gene are located on the light strand. The overall base composition of U. thibetanus mupinensis mtDNA for the H-strand is: A, 31.2%; C, 25.4%; G, 15.5%; T, 27.9%. Guanine (G) is the rarest nucleotide and GC content is 40.9%. Nucleotide composition analysis reveals that the U. thibetanus mupinensis genome is biased towards AT rich; such an AT content is lower than that of U. americanus and higher than those of U. maritimus and U. arctos. The gene content and gene order of U. thibetanus mupinensis mtDNA is typical for vertebrates . With respect to the length of intergenic spacers and overlaps, U. thibetanus mupinensis has a rather compact genome similar to the other three Ursidae (Table 3 and Table 4). The U.thibetanus mupinensis mitochondrial genome shares high similarity with those of the other three Ursidae: U. americanus (91.46%), U. arctos (89.25%) and U. maritimus (87.66%).
Vertebrate animal mitochondrial genomes deviate from a random usage of nucleotide. Saccone et al. (1999)  used the formula, base-skew= (A-T/A+T) or (G-C/G+C), to evaluate the degree of the base bias, and found all the values of GC-skew were negative while all the values of AT-skew were positive in amniotes . The base bias overall (GC=-0.242; AT=0.056) of the U.thibetanus mupinensis is most similar to that of U. arctos (GC=-0.238; AT=0.053) (Table 4).
Components of U.thibetanus mupinensis Mitochondrial Genome
|Gene||Nucleotide number||Start codon||Stop codon||Size (bp)||aa||Strand1|
|Origin of L-strand replication||6149-6181||35||-|
1 H=heavy strand;L=light strand
The structure and annotation of the Ursus thibetanas mupinensis mitochondrial genome. (16S: 16S ribosomal RNA gene; 12S: 12S ribosomal RNA gene; ND1: NADH dehydrogenase subunit 1 gene; ND2: NADH dehydrogenase subunit 2 gene; ND3: NADH dehydrogenase subunit 3 gene; ND4L: NADH dehydrogenase subunit 4L gene; ND4: NADH dehydrogenase subunit 4 gene; ND5: NADH dehydrogenase subunit 5 gene; ND6: NADH dehydrogenase subunit 6 gene; COX1: Cytochrome c oxidase subunit 1 gene; COX2: Cytochrome c oxidase subunit 2 gene; COX3: Cytochrome c oxidase subunit 3 gene; ATP8: ATP synthase F0 subunit 8 gene; ATP6: ATP synthase F0 subunit 6 gene; Cyt B: Cytochrome b gene; Phe: Phenylalanine tRNA gene; Val: Valine tRNA gene; Leu: Leucine tRNA gene; Ile:Isoleucine tRNA gene; Gln: Glutamine tRNA gene; Met: Trp: Tryptophan tRNA gene; Ala: Alanine tRNA gene; Tyr: Tyrosine tRNA gene; Ser: Serine tRNA gene; Asn: Asparagine tRNA gene; Cys: Cysteine tRNA gene; Asp: Aspartic acid tRNA gene; Lys: Lysine tRNA gene; Gly: Glycin tRNA gene; Arg: Arginine tRNA gene; His: Histidine tRNA gene; Glu: Glutamic acid tRNA gene; Thr: Threonine tRNA gene; Pro: Proline tRNA gene)
The mtDNA encodes 13 protein genes. Eight tRNA genes and one protein gene are located on the light strand . We found the direction or the encoding-strand selection of the genes of U.thibetanus mupinensis is identical to the typical vertebrates (Table 3). Among 13 protein genes of U.thibetanus mupinensis, ten use ATG as start codon; NADH2, NADH3 and NADH5 use ATA, ATC and ATT as start codon, respectively. Some of these 13 protein genes are terminated with incomplete stop codons: NADH5 and NADH3 are terminated with TA; COXIII (cytochrome c oxidase subunit III) and NADH4 are terminated with T; the rest are terminated with TAA and AGA (Table 1). All these start codons and stop codons are almost the same with those of U. americanus, U. arctos and U. maritimus except that U. americanus, U. arctos and U. maritimus's ND5 use ATC as start codon, U. arctos and U. maritimus's ND4L use GTG as start codon and U. arctos and U. maritimus's ATP8 and COX2 are terminated with TAG. Presumably, these incomplete stop codons are accommodated post-transcriptionally in the mRNA maturation process, i.e. polyadenylation . The 13 protein genes' nucleotide sequence and amino acid sequence of U.t. mupinensis mitochondrial genome shares high similarity with those of the other three Ursidae: U. americanus (88.7% to 93.7% and 80.6% to 100%), U. arctos (84.3% to 92.9% and 83.6% to 98.98%) and U. maritimus (84.8% to 93.4% and 82.1% to 98.98%) (Table 5).
Comparisons of mitochondrial genome features in four Ursidae
|Genome character||Black bear||American black bear||Polar bear||Brown bear|
|GC base component||40.9%||40.6%||41.3%||41.3%|
|Base skew overall||GC=-0.242 AT=0.056||GC=-0.235 AT=0.048||GC=-0.235 AT=0.053||GC=-0.238 AT=0.053|
The 37 gene's similarity comparison at the nt and aa level from between U.thibetanus mupinensis mtDNA with Ursus Americanus, Ursus maritimus and Ursus arctos
|Gene||Ursus Americanus (American black bear)||Ursus maritimus (polar bear)||Ursus arctos (brown bear)|
There are 22 tRNA genes identified in the U. thibetanus mupinensis genome, with length ranges from 59 to 75 bp. The 12S rRNA and 16S rRNA genes are 965 bp and 1 580bp in length, respectively. These are typical for mammalian mitochondrial genomes [16-19]. The 24 RNA genes of U.t. mupinensis mitochondrial genome shares high similarity with those of the other three Ursidae: U. americanus (78.46% to 100%), U. arctos (75.38% to 100%) and U. maritimus (75.38% to 100%), too (Table 5).
The stem-and-1oop structure of origin for the putative light-strand replication of U.thibetanus mupinensis located between tRNA-Asn and tRNA-Cys, is 35bp in size, which is the same size and location detected in U. americanus, U. arctos and U. maritimus. The sequence of the control region (CR) of U. thibetanus mupinensis located between tRNA-Pro and tRNA-Phe is 1422 bp in size and consists of 8.43% of the whole genome. This percentage is slightly lower than that of U. maritimus (9.25%) and U. arctos (9.27%), and higher than that of U. americanus (8.28%). Its GC content is 51.9%, which is lower than U. americanus's (52.6%) and higher than U. maritimus's (48.6%) and U. arctos's (48.2%). By using the Tandem Repeats Finder , we have found a 6bp tandem repeat (ACGTGT) and two 10bp tandem repeats (ACGTGTACGT and TACGTGTACG). But U. maritimus and U. arctos contain two 10bp tandem repeats (U. maritimus and U. arctos have the same two 10bp tandem repeats which are CGTACGCATA and CGCACGTGTA) and no 6bp tandem repeat, while U. americanus contains a 6bp tandem repeat (ACGTGT) and no two 10bp tandem repeats. The sequence of control region of U.thibetanus mupinensis shares 90% identity with that of U. americanus, 71% with that of U. arctos, and 68% with that of U. maritimus.
This work is supported by the Key Chinese National Natural Science Foundation (30570275), Application Technology Project in Sichuan Province (2006J13-057), Key Scientific Research Foundation of Educational Committee of Sichuan Province (2004A101) and Key Discipline Construction Project in Sichuan Province (SZD0420).
The author has declared that no conflict of interest exists.
1. Shadel GS, Clayton DA. Mitochondrial DNA maintenance in vertebrates. Annu. Rev. Biochem. 1997 ;66:409-435
2. Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999 ;27(8):1767-1780
3. Boore JL. Complete mitochondrial genome sequence of Urechis caupo, a representative of the phylum Echiura. BMC Genomics. 2004 ;5(1):67
4. Zardoya R, Meyer A. Complete mitochondrial genome suggests diapsid affinities of turtles. Proc Natl Acad Sci USA. 1998 ;95(24):14226-14231
5. Janke A, Erpenbeck D, Nilsson M. et al. The mitochondrial genomes of the iguana (Iguana iguana) and the caiman (Caiman crocodylus): Implications for amniote phylogeny. Proc. R. Soc. Lond. B Biol. Sci. 2001 ;268(1467):623-631
6. Zhou KY. Molecular phylogenetics of amphibians and reptiles. Zoological Research. 2001 ;22(5):397-405
7. Joshua SR, Jennifer CA, Christopher CA. et al. Molecular systematics of primary reptilian lineages and the tuatara mitochondrial genome. Mol. Phylogenet. Evol. 2003 ;29(2):289-297
8. Wu XB, Wang YQ, Zhou KY. et al. The complete mitochondrial genome sequence of Chinese alligator and the phylogeny of crocodiles [J]. Chinese Science Bulletin. 2003 ;48(18):1954-1958
9. Wang S. China Red Data Book of Endangered Animals: Mammalia. Beijing: Science Press. 1998
10. Delisle I, Strobeck C. Conserved primers for rapid sequencing of the complete mitochondrial genome from carnivores, applied to three species of bears. Mol. Biol. Evol. 2002 ;19(3):357-361
11. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual 2nd edition. New York: Cold Spring Harbor Laboratory Press. 1989
12. Masuda R, Yoshida MC. A molecular phylogeny of the family Mustelidae (Mammalia, Carnivora), based on comparison of mitochondrial cytochrome b nucleotide sequences. Zool. Sci. 1994 ;11:605-612
13. Benson G. Tandem Repeats Finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999 ;27:573-580
14. Saccone S, Giorgi CD, Gissi C. Evolutionary genomics in Metazoa: The mitochondrial DNA as a model system. Gene. 1999 ;238(1):195-209
15. Ojala D. et al. tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981 ;290:470-474
16. Hiendleder S. et al. The complete mitochondrial DNA sequence of the domestic sheep (Ovis aries) and comparison with the other major ovine haplotype. J. Mol. Evol. 1998 ;47:441-448
17. Anderson S. et al. Sequence and organization of the human mitochondrial genome. Nature. 1981 ;290:457-465
18. Lin CS. et al. Complete nucleotide sequence of pig (Sus scrofa) mitochondrial genome and dating evolutionary divergence within Artiodactyla. Gene. 1999 ;236:107-114
19. Parma P. et al. The complete nucleotide sequence of goat (Capra hircus) mitochondrial genome. Goat mitochondrial genome. DNA Seq. 2003 ;14:199-20
Correspondence to: Wan-ru Hou, College of Life Science, China West Normal University, 44# Yuying Road, 637002, Nanchong, China. Phone/Fax: +86-0817-2314311. E-mail: hwr168com.cn.