A mitochondrial control region and cytochrome b phylogeny of sika deer (Cervus nippon) and report of tandem repeats in the control region

First description of the mitochondrial genomes of the Central American brocket deer Mazama temama (Kerr, 1792) and the Yucatán Peninsula brocket deer Odocoileus pandora Merriam, 1901

BACKGROUND

The Central American (Mazama temama) and the Yucatán Peninsula brocket deer (Odocoileus pandora) are deer species with cryptic habits, and little is known about their biology. Odocoileus pandora is listed as Vulnerable on the 2015 IUCN Red List of Threatened Species, while M. temama is considered Data Deficient; however, it currently faces a decreasing population trend.

METHODS AND RESULTS

We assembled the complete mitochondrial genome for two M. temama specimens and one complete and one partial for O. pandora from Illumina 150 bp paired-end reads. The mitogenomes of M. temama and O. pandora have a length of 16,479-16,480 and 16,419 bp, respectively, AT-biased; they consist of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and one non-coding control region, most of them follow a transcription direction in the heavy strand of the molecule. The mitochondrial genome of O. pandora shows some particularities compared to other deer species, like a shorter control region of 987-990 bp and a cytochrome b gene with a length of 1,143 bp. Our phylogenetic analyses confirm the close affinity of M. temama to South American M. americana and the nested position of the genus Odocoileus, including O. pandora, into the genus Mazama.

CONCLUSIONS

Here, we described for the first time the complete mitochondrial genome for these two species. While our study provides additional information about the taxonomic status of the northern neotropical brocket deer, further research is needed to solve the complicated taxonomy of neotropical deer.

Molecular Evidence Reveals the Sympatric Distribution of Cervus nippon yakushimae and Cervus nippon taiouanus on Jeju Island, South Korea

Non-native species threaten native ecosystems and species, particularly on islands where rates of endemism and vulnerability to threats are high. Understanding species invasion will aid in providing insights into ecological and evolutionary processes. To identify the non-native sika deer (Cervus nippon) population in Jeju, South Korea, and their phylogenetic affinities, we collected tissue samples from roadkill and the World Natural Heritage Headquarters in Jeju. Mitochondrial DNA cytochrome B (CytB) gene sequences were analyzed to determine two distinct CytB haplotypes. Phylogenetic analysis using maximum likelihood tree revealed two haplotypes of CytB clustered into two different groups representing two subspecies: C. n. yakushimae, native to Japan, and C. n. taiouanus, native to Taiwan. The tentative divergence time between the two subspecies was estimated at 1.81 million years. Our study confirmed that the two subspecies of sika deer are sympatric in the natural ecosystem of Jeju Island. This study provides valuable information to help government and conservation agencies understand alien species and determine control policies for conserving native biodiversity in South Korea.

Characterization of the complete mitochondrial genome of Myrmuslateralis (Heteroptera, Rhopalidae) and its implication for phylogenetic analyses

Mitochondrial genomes (mitogenomes) are widely used in research studies on phylogenetic relationships and evolutionary history. Here, we sequenced and analyzed the mitogenome of the scentless plant bug Myrmuslateralis Hsiao, 1964 (Heteroptera, Rhopalidae). The complete 17,309 bp genome encoded 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region. The mitogenome revealed a high A+T content (75.8%), a positive AT-skew (0.092), and a negative GC-skew (-0.165). All 13 PCGs were found to start with ATN codons, except for cox1, in which TTG was the start codon. The Ka/Ks ratios of 13 PCGs were all lower than 1, indicating that purifying selection evolved in these genes. All tRNAs could be folded into the typical cloverleaf secondary structure, except for trnS1 and trnV, which lack dihydrouridine arms. Phylogenetic trees were constructed and analyzed based on the PCG+rRNA from 38 mitogenomes, using maximum likelihood and Bayesian inference methods, showed that M.lateralis and Chorosomamacilentum Stål, 1858 grouped together in the tribe Chorosomatini. In addition, Coreoidea and Pyrrhocoroidea were sister groups among the superfamilies of Trichophora, and Rhopalidae was a sister group to Alydidae + Coreidae.

Development and validation of a 1 K sika deer (Cervus nippon) SNP Chip

Background

China is the birthplace of the deer family and the country with the most abundant deer resources. However, at present, China’s deer industry faces the problem that pure sika deer and hybrid deer cannot be easily distinguished. Therefore, the development of a SNP identification chip is urgently required.

Results

In this study, 250 sika deer, 206 red deer, 23 first-generation hybrid deer (F1), 20 s-generation hybrid deer (F2), and 20 third-generation hybrid deer (F3) were resequenced. Using the chromosome-level sika deer genome as the reference sequence, mutation detection was performed on all individuals, and a total of 130,306,923 SNP loci were generated. After quality control filtering was performed, the remaining 31,140,900 loci were confirmed. From molecular-level and morphological analyses, the sika deer reference population and the red deer reference population were established. The Fst values of all SNPs in the two reference populations were calculated. According to customized algorithms and strict screening principles, 1000 red deer-specific SNP sites were finally selected for chip design, and 63 hybrid individuals were determined to contain red deer-specific SNP loci. The results showed that the gene content of red deer gradually decreased in subsequent hybrid generations, and this decrease roughly conformed to the law of statistical genetics. Reaction probes were designed according to the screening sites. All candidate sites met the requirements of the Illumina chip scoring system. The average score was 0.99, and the MAF was in the range of 0.3277 to 0.3621. Furthermore, 266 deer (125 sika deer, 39 red deer, 56 F1, 29 F2,17 F3) were randomly selected for 1 K SNP chip verification. The results showed that among the 1000 SNP sites, 995 probes were synthesized, 4 of which could not be typed, while 973 loci were polymorphic. PCA, random forest and ADMIXTURE results showed that the 1 K sika deer SNP chip was able to clearly distinguish sika deer, red deer, and hybrid deer and that this 1 K SNP chip technology may provide technical support for the protection and utilization of pure sika deer species resources.

Conclusion

We successfully developed a low-density identification chip that can quickly and accurately distinguish sika deer from their hybrid offspring, thereby providing technical support for the protection and utilization of pure sika deer germplasm resources.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00994-z.

Variation of craniometrical characters in an artificial population of Cervus nippon from Askania-Nova, Ukraine

Craniometrical characters of adult sika deer from an artificial population, which was created by import of animals from a natural population from the Primorsky Krai of Russia, are analysed. It was revealed that the mean value of variation (by 11 characters) of skulls of sika deer is 3.56 ± 0.35 in males and 3.89 ± 0.33 in females. Males are significantly larger than females (p < 0.001) by all characters. The study of craniometrical characters of the sika deer using principal component analysis showed that 93.3 % of total variance is described by the first principal component (PC1). The highest, though quite equal, factor loadings on PC1 have characters of length (full length, condylobasal length, basic length). Study results suggest that sex-related differences in C. nippon by craniometrical characters are expressed by general linear dimensions and by proportions of the skull. The correlation matrix of craniometrical characters of adult sika deer from the Askanian subpopulation shows a relatively weak dependence between different parameters. It was established that the coefficient of correlation in most pairs of characters varies from 0.02 to 0.96 in males and from 0.01 to 0.97 in females. The highest coefficients of correlation in both males and females (r = 0.96 and r = 0.97, respectively) are revealed between the condylobasal and basal lengths of the skull. Correlations between cranial indices are significantly lower compared to craniometrical characters. The highest coefficient of correlation in both males and females has the variation of the index of mastoid width in relation to the relative largest skull width (r = 0.72). Skulls of sika deer from the Askanian subpopulation (of both males and females) have significantly larger dimensions compared to those in animals from the natural population from the Primorsky Krai.

Increased genetic marker density reveals high levels of admixture between red deer and introduced Japanese sika in Kintyre, Scotland

Hybridization is a natural process at species range boundaries, but increasing numbers of species are hybridizing due to direct or indirect human activities. In such cases of anthropogenic hybridization, subsequent introgression can threaten the survival of native species. To date, many such systems have been studied with too few genetic markers to assess the level of threat resulting from advanced backcrossing. Here, we use 44,999 single nucleotide polymorphisms (SNPs) and the ADMIXTURE program to study two areas of Scotland where a panel of 22 diagnostic microsatellites previously identified introgression between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon). In Kintyre, we reclassify 26% of deer from the pure species categories to the hybrid category whereas in the NW Highlands we only reclassify 2%. As expected, the reclassified individuals are mostly advanced backcrosses. We also investigate the ability of marker panels selected on different posterior allele frequency criteria to find hybrids assigned by the full marker set and show that in our data, ancestry informative markers (i.e. those that are highly differentiated between the species, but not fixed) are better than diagnostic markers (those markers that are fixed between the species) because they are more evenly distributed in the genome. Diagnostic loci are concentrated on the X chromosome to the detriment of autosomal coverage.

Population level mitogenomics of long-lived bats reveals dynamic heteroplasmy and challenges the Free Radical Theory of Ageing

Bats are the only mammals capable of true, powered flight, which drives an extremely high metabolic rate. The "Free Radical Theory of Ageing" (FTRA) posits that a high metabolic rate causes mitochondrial heteroplasmy and the progressive ageing phenotype. Contrary to this, bats are the longest-lived order of mammals given their small size and high metabolic rate. To investigate if bats exhibit increased mitochondrial heteroplasmy with age, we performed targeted, deep sequencing of mitogenomes and measured point heteroplasmy in wild, long lived Myotis myotis. Blood was sampled from 195 individuals, aged between <1 and at 6+ years old, and whole mitochondria deep-sequenced, with a subset sampled over multiple years. The majority of heteroplasmies were at a low frequency and were transitions. Oxidative mutations were present in only a small number of individuals, suggesting local oxidative stress events. Cohort data showed no significant increase in heteroplasmy with age, while longitudinal data from recaptured individuals showed heteroplasmy is dynamic, and does not increase uniformly over time. We show that bats do not suffer from the predicted, inevitable increase in heteroplasmy as posited by the FRTA, instead heteroplasmy was found to be dynamic, questioning its presumed role as a primary driver of ageing.

Introgression of exotic Cervus (nippon and canadensis) into red deer (Cervus elaphus) populations in Scotland and the English Lake District

Since the mid-19th century, multiple introductions of Japanese sika deer (Cervus nippon nippon) and North American wapiti (C. canadensis) have taken place in the British Isles. While wapiti have generally been unsuccessful, sika have been very successful, especially in Scotland where they now overlap at least 40% of the range of native red deer (C. elaphus). Hybridization between these two species and red deer has been demonstrated in captivity and in the wild. Using a panel of 22 microsatellite loci that are highly diagnostic between red deer and sika, and moderately diagnostic between red deer and wapiti, we investigated the extent of introgression between these species in 2,943 deer sampled from around Scotland and from the English Lake District using the Bayesian clustering software STRUCTURE. We also used a diagnostic mitochondrial marker for red deer and sika. Our survey extends previous studies indicating little introgression of wapiti nuclear alleles into red deer, in particular in Northern Scotland, Kintyre, and the Lake District. We found a new area of extensive sika introgression in South Kintyre. In the North Highlands, we show for the first time geographically scattered evidence of past hybridization followed by extensive backcrossing, including one red-like individual with sika introgression, two sika-like individuals with red deer introgression, and six individuals that were apparently pure sika at the nuclear markers assessed but which carried red deer mitochondria. However, there has not been a collapse of assortative mating in this region. Similarly, in the English Lake District red deer, we found only traces of past sika introgression. No sika alleles were detected in the Central Highlands or the Hebridean red deer refugia. We make suggestions for management to prevent further spread of sika alleles into red deer and vice versa.

Phylogenetic relationships between Dicrocoelium chinensis populations in Japan and China based on mitochondrial nad1 gene sequences

We carried out phylogenetic analyses of the relationships between Dicrocoelium chinensis populations in Japan and China using molecular markers. One hundred nine lancet flukes collected from Japan and China were identified as D. chinensis based on their testis orientation and the nucleotide sequences of their ribosomal ITS2. These flukes were analyzed phylogenetically using mitochondrial nad1 gene sequences. An analysis of molecular variance found that the percentage of variation between the countries was extremely high, indicating that the D. chinensis populations in Japan and China are differentiated genetically. D. chinensis mainly parasitizes wild sika deer, which is thought to originate in northeast Asia and to have colonized into Japan from the Eurasia continent in the Pleistocene glaciations. In addition, phylogenic analyses indicated that Japanese sika deer is genetically differentiated from Chinese population; therefore, we hypothesize that D. chinensis might have been introduced into Japan along with the migration of infected wild ruminants in the Pleistocene, and then the population became differentiated from the Chinese population. This study provides the nucleotide sequences of the nad1 gene of D. chinensis in Japan for the first time and shows that these sequences are useful for elucidating the phylogenetic relationships of the Dicrocoelium species prevalent in Asia.

PCR-fingerprint profiles of mitochondrial and genomic DNA extracted from Fetus cervi using different extraction methods

The use of Fetus cervi, which is derived from the embryo and placenta of Cervus Nippon Temminck or Cervs elaphus Linnaeus, has been documented for a long time in China. There are abundant species of deer worldwide. Those recorded by China Pharmacopeia (2010 edition) from all the species were either authentic or adulterants/counterfeits. Identification of their origins or authenticity became a key in the preparation of the authentic products. The traditional SDS alkaline lysis and salt-outing methods were modified to extract mt DNA and genomic DNA from fresh and dry Fetus cervi in addition to Fetus from false animals, respectively. A set of primers were designed by bioinformatics to target the intra-and inter-variation. The mt DNA and genomic DNA extracted from Fetus cervi using the two methods meet the requirement for authenticity. Extraction of mt DNA by SDS alkaline lysis is more practical and accurate than extraction of genomic DNA by salt-outing method. There were differences in length and number of segments amplified by PCR between mt DNA from authentic Fetus cervi and false animals Fetus. The distinctive PCR-fingerprint patterns can distinguish the Fetus cervi from adulterants and counterfeit animal Fetus.

Revisiting Coleoptera a + T-rich region: structural conservation, phylogenetic and phylogeographic approaches in mitochondrial control region of bioluminescent Elateridae species (Coleoptera)

The control region (CR) or A + T-rich region in Coleoptera mt genome is poorly characterized, including the Elateroidea bioluminescent species. Here, we provided the first attempt to characterize and compare the structure and organization of the CR of different species within Elateridae. We also revisited some sequenced Coleoptera CR and observed consensus T-stretches, non-conserved sequences near the stem-loop and unusual inner tRNAs-like sequences. All these features are probably involved in the replication start of the mt genome. The phylogenetic relationships in Elateridae bioluminescent groups using partial sequence of CR showed the monophyly of Pyrearinus pumilus group and Pyrearinus as a polyphyletic genus, corroborating our previous results. The wider genetic variation obtained by CR analysis could separate two different lineages that occur within P. termitilluminans populations. In Elateridae, the CR exhibited high polymorphism within and between populations, which was also observed in other Coleoptera species, suggesting that the CR could be described as a suitable molecular marker to be applied in phylogenetic and phylogeographic studies.

The complete mitochondrial genome sequence of the tubeworm Lamellibrachia satsuma and structural conservation in the mitochondrial genome control regions of Order Sabellida

The control region of the mitochondrial genomes shows high variation in conserved sequence organizations, which follow distinct evolutionary patterns in different species or taxa. In this study, we sequenced the complete mitochondrial genome of Lamellibrachia satsuma from the cold-seep region of Kagoshima Bay, as a part of whole genome study and extensively studied the structural features and patterns of the control region sequences. We obtained 15,037 bp of mitochondrial genome using Illumina sequencing and identified the non-coding AT-rich region or control region (354 bp, AT=83.9%) located between trnH and trnR. We found 7 conserved sequence blocks (CSB), scattered throughout the control region of L. satsuma and other taxa of Annelida. The poly-TA stretches, which commonly form the stem of multiple stem-loop structures, are most conserved in the CSB-I and CSB-II regions. The mitochondrial genome of L. satsuma encodes a unique repetitive sequence in the control region, which forms a unique secondary structure in comparison to Lamellibrachia luymesi. Phylogenetic analyses of all protein-coding genes indicate that L. satsuma forms a monophyletic clade with L. luymesi along with other tubeworms found in cold-seep regions (genera: Lamellibrachia, Escarpia, and Seepiophila). In general, the control region sequences of Annelida could be aligned with certainty within each genus, and to some extent within the family, but with a higher rate of variation in conserved regions.

Detection of 40 bp insertion-deletion (INDEL) in mitochondrial control region among sambar (Rusa unicolor) populations in India

Background

The mitochondrial DNA (mtDNA) control region is extensively used in the phylogeography of species. We examined sequence variations in the mtDNA control region of sambar (Rusa unicolor) populations from the South, Central and North India.

Results

Most of the samples collected from the south India exhibited a 40 bp insertion in the mtDNA control region. This insertion was not observed in the North and Central Indian populations.

Conclusion

This study provided a potential marker for molecular screening and identification of sambar populations in the form of a distinct 40 bp insertion. Some populations in South India did not exhibit this insertion. It indicates that there could be an ecological barrier that might be preventing the expansion of insertion-positive sambar population.

Development of multiplex PCR assay for authentication of Cornu Cervi Pantotrichum in traditional Chinese medicine based on cytochrome b and C oxidase subunit 1 genes

This study describes a method for discriminating the true Cervus antlers from its counterfeits using multiplex PCR. Bioinformatics were carried out to design the specific alleles primers for mitochondrial (mt) cytochrome b (Cyt b) and cytochrome C oxidase subunit 1 (Cox 1) genes. The mt DNA and genomic DNA were extracted from Cervi Cornu Pantotrichum through the modified alkaline and the salt-extracting method in addition to its counterfeits, respectively. Sufficient DNA templates were extracted from all samples used in two methods, and joint fragments of 354 bp and 543 bp that were specifically amplified from both of true Cervus antlers served as a standard control. The data revealed that the multiplex PCR-based assays using two primer sets can be used for forensic and quantitative identification of original Cervus deer products from counterfeit antlers in a single step.

An examination of the origin and evolution of additional tandem repeats in the mitochondrial DNA control region of Japanese sika deer (Cervus Nippon)

Tandem repeat units are only detected in the left domain of the mitochondrial DNA control region in sika deer. Previous studies showed that Japanese sika deer have more tandem repeat units than its cousins from the Asian continent and Taiwan, which often have only three repeat units. To determine the origin and evolution of these additional repeat units in Japanese sika deer, we obtained the sequence of repeat units from an expanded dataset of the control region from all sika deer lineages. The functional constraint is inferred to act on the first repeat unit because this repeat has the least sequence divergence in comparison to the other units. Based on slipped-strand mispairing mechanisms, the illegitimate elongation model could account for the addition or deletion of these additional repeat units in the Japanese sika deer population. We also report that these additional repeat units could be occurring in the internal positions of tandem repeat regions, possibly via coupling with a homogenization mechanism within and among these lineages. Moreover, the increased number of repeat units in the Japanese sika deer population could reflect a balance between mutation and selection, as well as genetic drift.

Faunal record identifies Bering isthmus conditions as constraint to end-Pleistocene migration to the New World

Human colonization of the New World is generally believed to have entailed migrations from Siberia across the Bering isthmus. However, the limited archaeological record of these migrations means that details of the timing, cause and rate remain cryptic. Here, we have used a combination of ancient DNA, 14C dating, hydrogen and oxygen isotopes, and collagen sequencing to explore the colonization history of one of the few other large mammals to have successfully migrated into the Americas at this time: the North American elk (Cervus elaphus canadensis), also known as wapiti. We identify a long-term occupation of northeast Siberia, far beyond the species's current Old World distribution. Migration into North America occurred at the end of the last glaciation, while the northeast Siberian source population became extinct only within the last 500 years. This finding is congruent with a similar proposed delay in human colonization, inferred from modern human mitochondrial DNA, and suggestions that the Bering isthmus was not traversable during parts of the Late Pleistocene. Our data imply a fundamental constraint in crossing Beringia, placing limits on the age and mode of human settlement in the Americas, and further establish the utility of ancient DNA in palaeontological investigations of species histories.

Dybowski's sika deer (Cervus nippon hortulorum): genetic divergence between natural primorian and introduced Czech populations

Dybowski's sika deer (Cervus nippon hortulorum) originally inhabited the majority of the Primorsky Krai in Far Eastern Russia, north-eastern China, and Korean Peninsula. At present, only the Russian population seems to be stable, even though this taxon is still classified as endangered by the Russian Federation. Almost 100 years ago, this subspecies, among others, was imported to several European countries including the Czech Republic. We used both mitochondrial (mtDNA; the cytochrome b gene and the control region) and nuclear DNA markers to examine the actual taxonomic status of modern Czech Dybowski's sika population and to compare the genetic diversity between the introduced and the native populations. Altogether, 124 Czech samples and 109 Primorian samples were used in the analyses. Within the samples obtained from individuals that were all morphologically classified as Dybowski's sika, we detected mtDNA haplotypes of Dybowski's sika (84 samples), as well as those belonging to other sika subspecies: northern Japanese sika (25 samples), southern Japanese sika (6 samples), and south-eastern Chinese sika (8 samples). Microsatellite analysis revealed a certain level of heterozygote deficiency and a high level of inbreeding in both populations. The high number of private alleles, factorial correspondence analysis, and Bayesian clustering analysis indicate a high level of divergence between both populations. The large degree of differentiation and the high number of population-specific alleles could be a result of a founder effect, could be a result of a previously suggested bottleneck within the Primorian population, and could also be affected by the crossbreeding of captive individuals with other sika subspecies.

The complete mitochondrial genome of the domestic red deer (Cervus elaphus) of New Zealand and its phylogenic position within the family Cervidae

We determined the complete nucleotide sequence of the mitochondrial genome of the semidomestic red deer (Cervus elaphus) of New Zealand. The genome was 16,357 bp long and contained 13 protein-coding genes, 12SrRNA, 16SrRNA, 22 tRNAs and a D-loop as found in other mammals. Database homology searches showed that the mitochondrial DNA (mtDNA) sequence from the New Zealand semidomestic deer was similar to partial mtDNA sequences from the European, Norwegian (C. e. atlanticus) and Spanish red deer (C. e. hispanicus). Phylogenetic analysis of the mitochondrial protein-coding regions revealed two well-defined monophyletic clades in subfamilies Cervinae and Muntiacinae. However, red deer and Sika deer were not found to be close relatives. The analysis did identify the red deer as a sister taxon of a Samber/Sika deer clade, although it was more closely related to the Samber than the Sika group.

Mitochondrial cytochrome b sequence variation and phylogenetics of the highly specialized Schizothoracine fishes (Teleostei: Cyprinidae) in the Qinghai-Tibet plateau

The complete 1140 bp mitochondrial cytochrome b sequences were obtained from 39 individuals representing five species of all four genera of highly specialized schizothoracine fishes distributed in the Qinghai-Tibet plateau. Sequence variation of the cytochrome b gene was surveyed among the 39 individuals as well as three primitive schizothoracines and one outgroup. Phylogenetic analysis suggested that the group assignment based on 1140 bp of the cytochrome b sequence is obviously different from previous assignments, and the highly specialized schizothoracine fishes (Schizopygopsis pylzovi, Gymnocypris przewalskii, G. eckloni, Chuanchia labiosa, and Platypharodon extremus) form a monophyletic group that is sister to the clade formed by the primitive schizothoracine fishes (Schizothorax prenanti, S. pseudaksaiensis, and S. argentatus). The haplotypes of Schizopygopsis pylzovi and G. przewalskii were paraphyletic based on cytochrome b data, which most likely reflected incomplete sorting of mitochondrial DNA lineages. The diploid chromosome numbers of Schizothoracinae were considered in phylogenetic analysis and provided a clear pattern of relationships. Molecular dating estimated for highly specialized schizothoracine fishes suggested that the highly specialized schizothoracine fishes diverged in the late Miocene Pliocene to Pleistocene (4.5 x 10(4)-4.05 x 10(6 ) years BP). The relationship between the cladogenesis of highly specialized schizothoracine fishes and geographical events of the Qinghai-Tibet plateau is discussed.

Evolution and phylogeny of old world deer

The phylogenetic pattern and timing of the radiation of Old World deer was determined based on the complete mitochondrial cytochrome b gene from 33 Cervinae taxa. Using rooted and unrooted phylogenies derived from distinct theoretical approaches, strong support was achieved for monophyly of the Old World deer with muntjacs as sister group as well as for the divergence of at least three distinct genera: Rucervus, Dama, and Cervus. The latter clade comprises what have previously been regarded as the genera or subgenera Panolia, Rusa, Cervus, Sika, and probably Przewalskium. Our data also consistently confirmed paraphyly of nominate C. elaphus and did not support the monophyly of Axis. We used these molecular phylogenies to assess the homoplastic evolution of morphological, geographical, ecological, and selected behavioural character state differences within the Cervinae. Reliable fossil calibrations, large molecular data sets, and improved dating methods are shaping a molecular time scale for the evolutionary radiation of Old World deer that occurred at the Miocene/Pliocene transition and is largely compatible with existing palaeontological evidence. Using node ages estimated from sequence data, we estimated an average per-lineage diversification rate of 0.51+/-0.1 species per million years (my) over roughly the last 6 mya.

Mitochondrial DNA phylogeography of red deer (Cervus elaphus)

In order to understand the origin, phylogeny, and phylogeography of the species Cervus elaphus, we examined the DNA sequence variation of the mitochondrial cytochrome b gene of 51 populations of deer from the entire distribution area of Cervinae with an emphasis on Europe and Asia. Several methods, including maximum parsimony, maximum likelihood, and nested clade analysis, revealed that red deer originated from the area between Kyrgyzstan and Northern India. We found two distinct groups of red deer: a western group consisting of four subgroups and an eastern group consisting of three subgroups. Our mtDNA data do not support the traditional classification of red deer as only one species nor its division into numerous subspecies. The discrepancies between the geographical pattern of differentiation based on mtDNA cytochrome b and the existing specific and subspecific taxonomy based on morphology are discussed.

Development of PCR assay for differentiation of some important wild animal meat of Sri Lanka

A polymerase chain reaction (PCR) assay was developed to differentiate meat of Ceylon spotted deer (Axis axis ceylonensis), Ceylon hog deer (A. porcius oryzus), Ceylon sambhur (Cervus unicolor unicolor) and barking deer (Muntiacus muntijak malabaricus) from meat of cattle, goat, buffalo, pig, dog and sheep. A set of primers was designed according to the sequence of the mitochondrial cytochrome b gene of C. elaphus canadensis and by PCR amplification about 450 bp band was observed for all four animal species and these primers were not cross reacted with DNA of other animal species tested in the study under the tested reaction conditions. A band of 649 bp size was observed for all animal species when DNA was amplified with the universal primers and that indicated the presence of mitochondrial DNA in the samples. Further, the results indicated that this technique was sensitive enough to differentiate rotten meat, at least 5 days after the killing of an animal. Under these PCR conditions, the DNA of bacteria, which is involved in decomposition of meat, was not amplified with both universal and specific primers. However, the method was not sensitive enough in differentiating cooked meat of these species. Slaughtering of these four wild animal species is banned, but the animals are being killed illegally. Lack of meat identification methods has been identified as one of the major constraints to implement legal procedures and conserve biodiversity in the country.

A phylogenetic comparison of red deer and wapiti using mitochondrial DNA

A phylogeny was constructed for red deer/wapiti (Cervus elaphus) subspecies using sequence data from the control region of mitochondrial DNA (mtDNA). The tree was rooted using Cervus nippon (sika deer), Cervus albirostris (Thorold's white-lipped deer), and several Odocoileinae species. A division between the mtDNA haplotypes of red deer (European) and wapiti (Asian/North American) corresponds to subspecies found on opposite sides of the Himalayan Mountains and Gobi, which suggests wapiti should be reconsidered for the status of C. canadensis. Using parsimony and distance analysis, red deer and wapiti are derived from a single recent common ancestor, which is consistent with current taxonomy that recognizes the subspecies of Cervus elaphus as monophyletic group. However, maximum-likelihood analysis using weighted transitional substitutions caused red deer to form a sister group to sika deer (Cervus nippon) and wapiti. A phenetic comparison revealed wapiti also share more nucleotide similarities with sika deer, although approximately 5% sequence divergence separates wapiti, sika, and red deer. Phylogenetic evidence from the cytochrome b sequences corroborated observations from the control region. Observations from this study suggest that the species status of wapiti should be reinstated.

Cytochrome b gene for species identification of the conservation animals

A partial DNA sequence of cytochrome b gene was used to identify the remains of endangered animals and species endemic to Taiwan. The conservation of animals species included in this study were: the formosan gem-faced civets, leopard cats, tigers, clouded leopards, lion, formosan muntjacs, formosan sika deers, formosan sambars, formosan serows, water buffalo, formosan pangolins and formosan macaques. The control species used included domestic cats, domestic dogs, domestic sheeps, domestic cattles, domestic pigs and humans. Heteroplasmy was detected in the formosan macaque, domestic pig and domestic cats. The frequencies of heteroplasmy in these animals were about 0.25% (1 in 402bp). Sequences were aligned by Pileup program of GCG computer package, and the phylogenetic tree was constructed by the neighbor-joining method. The results of sequence comparison showed that the percentage range of sequence diversity in the same species was from 0.25 to 2.74%, and that between the different species was from 5.97 to 34.83%. The results of phylogenetic analysis showed that the genetic distance between the different species was from 6.33 to 40.59. Animals of the same species, both the endangered animal species and domestic animals, were clustered together in the neighbor-joining tree. Three unknown samples of animal remains were identified by this system. The partial sequence of cytochrome b gene adopted in this study proved to be usable for animal identification.

Bottlenecks, drift and differentiation: the population structure and demographic history of sika deer (Cervus nippon) in the Japanese archipelago

We assessed genetic differentiation and diversity in 14 populations of sika deer (Cervus nippon) from Japan and four populations of sika deer introduced to the UK, using nine microsatellite loci. We observed extreme levels of differentiation and significant differences in diversity between populations. Our results do not support morphological subspecies designations, but are consistent with previous mitochondrial DNA analyses which suggest the existence of two genetically distinct lineages of sika deer in Japan. The source of sika introduced to the UK was identified as Kyushu. The underlying structure of Japanese populations probably derives from drift in separate glacial refugia and male dispersal limited by distance. This structure has been perturbed by bottlenecks and habitat fragmentation, resulting from human activity from the mid-nineteenth century. Most current genetic differentiation and differences in diversity among populations probably result from recent drift. Coalescent model analysis suggests sika on each of the main Japanese islands have experienced different recent population histories. Hokkaido, which has large areas of continuous habitat, has maintained high levels of gene flow. In Honshu the population is highly fragmented and is likely to have been evolving by drift alone. In Kyushu there has been a balance between gene flow and drift but all the populations have experienced high levels of drift. Habitat fragment size was not significantly associated with genetic diversity in populations but there was a significant correlation between habitat fragment size and effective population size.