Phylogenetic study of complete cytochrome b genes in musk deer (genus Moschus) using museum samples

Predicting the current habitat refugia of Himalayan Musk deer (Moschus chrysogaster) across Nepal

Himalayan Musk deer, Moschus chrysogaster is widely distributed but one of the least studied species in Nepal. In this study, we compiled a total of 429 current presence points of direct observation of the species, pellets droppings, and hoofmarks based on field-based surveys during 2018-2021 and periodic data held by the Department of National Park and Wildlife Conservation. We developed the species distribution model using an ensemble modeling approach. We used a combination of bioclimatic, anthropogenic, topographic, and vegetation-related variables to predict the current suitable habitat for Himalayan Musk deer in Nepal. A total of 16 predictor variables were used for habitat suitability modeling after the multicollinearity test. The study shows that the 6973.76 km2 (5%) area of Nepal is highly suitable and 8387.11 km2 (6%) is moderately suitable for HMD. The distribution of HMD shows mainly by precipitation seasonality, precipitation of the warmest quarter, temperature ranges, distance to water bodies, anthropogenic variables, and land use and land cover change (LULC). The probability of occurrence is less in habitats with low forest cover. The response curves indicate that the probability of occurrence of HMD decreases with an increase in precipitation seasonality and remains constant with an increase in precipitation of the warmest quarter. Thus, the fortune of the species distribution will be limited by anthropogenic factors like poaching, hunting, habitat fragmentation and habitat degradation, and long-term forces of climate change.

Population genomics reveals moderate genetic differentiation between populations of endangered Forest Musk Deer located in Shaanxi and Sichuan

BACKGROUND

Many endangered species exist in small, genetically depauperate, or inbred populations, hence promoting genetic differentiation and reducing long-term population viability. Forest Musk Deer (Moschus berezovskii) has been subject to illegal hunting for hundreds of years due to the medical and commercial values of musk, resulting in a significant decline in population size. However, it is still unclear to what extent the genetic exchange and inbreeding levels are between geographically isolated populations. By using whole-genome data, we reconstructed the demographic history, evaluated genetic diversity, and characterized the population genetic structure of Forest Musk Deer from one wild population in Sichuan Province and two captive populations from two ex-situ centers in Shaanxi Province.

RESULTS

SNP calling by GATK resulted in a total of 44,008,662 SNPs. Principal component analysis (PCA), phylogenetic tree (NJ tree), ancestral component analysis (ADMIXTURE) and the ABBA-BABA test separated Sichuan and Shaanxi Forest Musk Deer as two genetic clusters, but no obvious genetic differentiation was observed between the two captive populations. The average pairwise FST value between the populations in Sichuan and Shaanxi ranged from 0.05-0.07, suggesting a low to moderate genetic differentiation. The mean heterozygous SNPs rate was 0.14% (0.11%-0.15%) for Forest Musk Deer at the genomic scale, and varied significantly among three populations (Chi-square = 1.22, p < 0.05, Kruskal-Wallis Test), with the Sichuan population having the lowest (0.11%). The nucleotide diversity of three populations varied significantly (p < 0.05, Kruskal-Wallis Test), with the Sichuan population having the lowest genetic θπ (1.69 × 10-3).

CONCLUSIONS

Genetic diversity of Forest Musk Deer was moderate at the genomic scale compared with other endangered species. Genetic differentiation between populations in Sichuan and Shaanxi may not only result from historical biogeographical factors but also be associated with contemporary human disturbances. Our findings provide scientific aid for the conservation and management of Forest Musk Deer. They can extend the proposed measures at the genomic level to apply to other musk deer species worldwide.

Mitochondrial capture and subsequent genetic divergence generates a novel haplogroup: evidence from ancient and ongoing hybridization in mule and white-tailed deer

Odocoileus virginianus (white-tailed deer) and O. hemionus (mule deer) are sympatric across much of North America. Molecular evidence suggests that up to 24% of individuals in some populations are a product of hybrid ancestry. Several studies have alluded to ancient and recent introgression between Odocoileus spp.; however, no divergence dates were proposed. Herein, phylogenetic analyses of DNA sequences obtained from the mitochondrial Cytochrome b gene in 690 individuals identified three clades corresponding to black-tailed deer, white-tailed deer, or a unique combination of both white-tailed deer and mule deer. White-tailed deer and mule deer diverged from a common ancestor of approximately 3.13 mya followed by an ancient hybridization event of approximately 1.32 mya, in which the white-tailed deer mitochondrial genome was “captured” by mule deer. This hybridization event produced a novel haplogroup for white-tailed deer and mule deer located west of the Appalachian Mountains and east of the Cascade Range, south to Veracruz, Mexico, and north to the Yukon Territory, Canada. The ancestral mule deer-like mitochondrial genome appears to be restricted to black-tailed deer distributed along the western portion of the Cascade and Sierra Nevada Ranges of the United States and Canada, whereas the ancestral white-tailed deer-like mitochondrial genome is restricted to the eastern United States and portions of Latin America and Caribbean regions. The “captured mitochondrial genome” has continued on an independent evolutionary trajectory and represents a unique and broadly distributed haplogroup that is 7.25% and 2.84% different from the ancestral mule deer and ancestral white-tailed deer haplogroups, respectively.

Mitochondrial genome of captive Alpine musk deer, Moschus chrysogaster (Moschidae), and phylogenetic analyses with its coordinal species

Alpine musk deer, Moschus chrysogaster, a solitary, primitive ungulate inhabiting high elevation areas (3000–4500 m) is an endangered species facing threat of extinction globally due to excessive hunting for its musk. In this study, we determined the complete mitochondrial genome of M. chrysogaster, which was 16,354 bp in length, and revealed the same gene order and genomic organization as typical Moschidae mitochondrial DNA. Start codons in 13 protein-coding genes (PCGs) were all typical ATGs except ATA for ND2 and ND3 and ATT for ND5. Stop codons were all typical types except an incomplete stop codon T for COX3, ND2, ND3, and ND4. Secondary structures in 22 transfer RNA genes all showed typical cloverleaf except tRNA-Ser (AGY), in which the dihydrouridine arm formed a simple loop. No repeat units were found in the control region. The topology structure indicated that M. cupreus was primitive and located at the root of the Moschidae clade. Phylogenetic reconstruction placed M. chrysogaster as a distinct lineage, closely related to the branch of M. leucogaster, M. berezovskii (wild) and predicted a sister relationship with M. moschiferus, M. anhuiensis, and M. berezovskii (captive). However, we suggested that the genetic resources of M. chrysogaster_JQ608470 should be further investigated.

Whole-genome sequencing of wild Siberian musk deer (Moschus moschiferus) provides insights into its genetic features

Background

Siberian musk deer, one of the seven species, is distributed in coniferous forests of Asia. Worldwide, the population size of Siberian musk deer is threatened by severe illegal poaching for commercially valuable musk and meat, habitat losses, and forest fire. At present, this species is categorized as Vulnerable on the IUCN Red List. However, the genetic information of Siberian musk deer is largely unexplored.

Results

Here, we produced 3.10 Gb draft assembly of wild Siberian musk deer with a contig N50 of 29,145 bp and a scaffold N50 of 7,955,248 bp. We annotated 19,363 protein-coding genes and estimated 44.44% of the genome to be repetitive. Our phylogenetic analysis reveals that wild Siberian musk deer is closer to Bovidae than to Cervidae. Comparative analyses showed that the genetic features of Siberian musk deer adapted in cold and high-altitude environments. We sequenced two additional genomes of Siberian musk deer constructed demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Finally, we identified several candidate genes that may play a role in the musk secretion based on transcriptome analysis.

Conclusions

Here, we present a high-quality draft genome of wild Siberian musk deer, which will provide a valuable genetic resource for further investigations of this economically important musk deer.

European minnows through time: museum collections aid genetic assessment of species introductions in freshwater fishes (Cyprinidae: Phoxinus species complex)

Massive fish introductions have taken place throughout much of the world, mostly over the last 70 years, and present a major threat to the genetic diversity of native fishes. Introductions have been reported for European Phoxinus, a ubiquitous small cyprinid that populates a wide variety of habitats. Species delineation in European Phoxinus has proven difficult with one reason being ranges of distribution that often traverse drainage boundaries. The present study combines recent samples with museum samples to better understand the current distribution of Phoxinus species and their distributions prior to the massive introductions of fishes in Europe, and to evaluate the use of museum specimens for species distribution studies. For these purposes, genetic lineages from sites collected prior to 1900 (n = 14), and between 1900 and 1950 (n = 8), were analysed using two mitochondrial and nuclear markers. Although possible fish introductions were detected, our results show that the distribution of genetic lineages of museum samples is comparable to that of the extant lineages of European Phoxinus present in those areas. These observations suggest that in the studied ranges the distribution of Phoxinus lineages has been driven by natural processes.

mtDNA analysis confirms the endangered Kashmir musk deer extends its range to Nepal

Musk deer Moschus spp. are endemic to the high mountain forests of central Asia. The taxonomic status of musk deer in the central and western Himalayas is poorly understood. We investigated the phylogenetic relationship of musk deer from the central and western Himalayas based on mitochondrial genomic data of Cytochrome b (380 bps) and D-loop (1000 bps). Our results distinguished two divergent lineages using higher bootstrap support (bs) values from the Maximum likelihood and Bayesian posterior probabilities (bpp). Both the Manang and Kaski lineages from central Nepal are confirmed as Himalayan musk deer Moschus leucogaster and represent a species complex widespread throughout the central and eastern Himalayan region. The musk deer Mustang lineage was confirmed as Kashmir musk deer Moschus cupreus and has wide distribution in the western Himalayas (from central Nepal to Afghanistan). Our analysis validates that Kashmir musk deer is a genetically distinct species and it clarifies that Himalayan musk deer and Kashmir musk deer are confirmed instead of Alpine musk deer Moschus chrysogaster which has been previously described from the southern parts of Himalayas of Nepal, India and Pakistan.

Influence of pairwise genetic distance computation and reference sample size on the reliability of species identification using Cyt b and COI gene fragments in a group of native passerines

Species identification is fundamental to wildlife forensic practice. The desirability of molecular genetic methods is increasing rapidly. The sequence of a marker, rather than its particular diagnostic nucleotides, provides greater safety through comparisons between intra- and inter-specific pairwise genetic distances. However, it has not been well described how reliability of species assignment is influenced by distance computing methods and reference sample sizes. In this study, the influences were tested using 12 species from 4 genera of passerine birds and the sequences of partial Cytochrome b (Cyt b) and Cytochrome Oxidase subunit I (COI) genes. Results showed that different substitution types have different outcomes of pairwise genetic distance estimation and this influences the risk of false inclusion and exclusion. Transition (Ts) is the most effective substitution type to reveal optimal species resolution for both Cyt b and COI gene fragments no matter whether K2P and p-distance are used. Sample size required to accurately estimate pairwise distance is essentially determined by the genetic diversity of a species in reference to a given strictness of predefined acceptable accuracy. These findings suggest that for future forensic work on birds by use of Cyt b and COI gene fragments, transition should be used exclusively for marker validation and identification practice when targeting closely related species. Meanwhile, the reference database should sufficiently represent overall genetic diversity of the species. The minimum sample size should be estimated based on existing knowledge of genetic diversity. Special caution should be used for species assignment when only several reference data are available for animals that are considered likely to have high genetic diversity.

Species dilemma of musk deer (Moschus spp) in India: molecular data on cytochrome c oxidase I suggests distinct genetic lineage in Uttarakhand compared to other Moschus species

Musk deer are of high conservation priority owing to poaching pressure because of its musk pod. Representation of musk deer status using genetics is poorly documented in India, and it is not confirmed as to how many species of musk deer are present. We characterize for the first time, the genetic diversity of musk deer from Uttarakhand using Cytochrome Oxidase sub-unit (COI) gene (486 bp) and compared with the data available for other species. Results revealed the presence of six haplotypes in the Uttarakhand population amongst 17 sequences. Of these, 12 sequences shared the single haplotype. The intra-species sequences divergence was 0.003-0.017, whereas divergence with other species of musk deer was 0.071-0.081. Bayesian phylogenetic tree revealed that samples from Uttarakhand formed a separate clade with respect to other species of musk deer, whereas three species distributed in China clustered in the same clade and showed low sequences divergence, i.e., 0.002-0.061. Because of different ecomorph reported, we suggest using the barcoding based approach for inter and intra-species distinction and delineating species boundaries across the range for effective conservation. Besides, systematic classification, DNA barcoding would also help in dealing wildlife offence cases for disposal of the legal report in court.

Next generation sequencing yields the complete mitogenome of captive forest musk deer, Moschus berezovskii (Ruminantia: Moschidae)

Moschus berezovskii is an endangered species, but its captive populations are valuable on musk secretions in traditional Chinese medicine and perfume manufacture. The mitogenome of M. berezovskii was 16,353 bp in size. Stop codons in 13 PCGs were all typical types except incomplete stop codon T for COX3, ND2 and ND4, and TA for ND3. No tandem repeat was found in control region. Phylogenetic analysis indicated that Moschidae has the closest relationship with Bovidae. We supported that M. berezovskii should be categorized into two subspecies, and suggested that the status of M. chrysogaster JQ608470 should be further investigated.

Micromeryx? eiselei-A new moschid species from Steinheim am Albuch, Germany, and the first comprehensive description of moschid cranial material from the Miocene of Central Europe

Moschids are enigmatic pecoran ruminants whose phylogeny is still not fully understood. So far we know only little of the family’s early evolutionary history and the origin of the modern genus, Moschus. Here we present a comprehensive description of cranial material, including the ear region and the dentition, of fossil moschid material from the Middle Miocene locality Steinheim am Albuch (13.5 Ma; Germany). This study provides the first exhaustive dataset for the cranial osteology of Micromeryx flourensianus, the most likely oldest true moschid. It furthermore reveals the presence of a second, so far undescribed moschid species, we here name Micromeryx? eiselei, in the abundant material from the locality. The two taxa can be clearly distinguished by characters of the skull, the ear region, the dentition, as well as by size. This evidences the sympatric occurrence of two moschid species in the locality Steinheim am Albuch.

Species Delimitation in the Genus Moschus (Ruminantia: Moschidae) and Its High-Plateau Origin

The authenticity of controversial species is a significant challenge for systematic biologists. Moschidae is a small family of musk deer in the Artiodactyla, composing only one genus, Moschus. Historically, the number of species in the Moschidae family has been debated. Presently, most musk deer species were restricted in the Tibetan Plateau and surrounding/adjacent areas, which implied that the evolution of Moschus might have been punctuated by the uplift of the Tibetan Plateau. In this study, we aimed to determine the evolutionary history and delimit the species in Moschus by exploring the complete mitochondrial genome (mtDNA) and other mitochondrial gene. Our study demonstrated that six species, M. leucogaster, M. fuscus, M. moschiferus, M. berezovskii, M. chrysogaster and M. anhuiensis, were authentic species in the genus Moschus. Phylogenetic analysis and molecular dating showed that the ancestor of the present Moschidae originates from Tibetan Plateau which suggested that the evolution of Moschus was prompted by the most intense orogenic movement of the Tibetan Plateau during the Pliocene age, and alternating glacial-interglacial geological eras.

Comparative molecular cytogenetics in Cetartiodactyla

Cetartiodactyla comprises Artiodactyla (even-toed ungulates) and Cetacea (whales, dolphins and porpoises). Artiodactyla is a large taxon represented by about 200 living species ranked in 10 families. Cetacea are classified into 13 families with almost 80 species. Many publications concerning karyotypic relationships in Cetartiodactyla have been published in previous decades. Formerly, the karyotypes of closely related species were compared by chromosome banding. Introduction of molecular cytogenetic methods facilitated comparative mapping between species with highly rearranged karyotypes and distantly related species. Such information is a prerequisite for the understanding of karyotypic phylogeny and the reconstruction of the karyotypes of common ancestors. This study summarizes the data on chromosome evolution in Cetartiodactyla, mainly derived from molecular cytogenetic studies. Traditionally, phylogenetic relationships of most groups have been estimated using morphological data. However, the results of some molecular studies of mammalian phylogeny are discordant with traditional conceptions of phylogeny. Cetartiodactyls provide several examples of incongruence between traditional morphological and molecular data. Such cases of conflict include the relationships of the major clades of artiodactyls, the relationships among the extant families of the suborder Ruminantia or the phylogeny of the family Bovidae. The most unexpected aspect of the molecular phylogeny was the recognition that Cetacea is a deeply nested member of Artiodactyla. The largest living order of terrestrial hoofed mammals is the even-toed hoofed mammals, or Artiodactyla. The artiodactyls are composed of over 190 living species including pigs, peccaries, hippos, camels, llamas, deer, pronghorns, giraffes, sheep, goats, cattle and antelopes. Cetacea is an order of wholly aquatic mammals, which include whales, dolphins and porpoises. Cetartiodactyla has become the generally accepted name for the clade containing both of these orders.

Phylogeographic analysis reveals significant spatial genetic structure of Incarvillea sinensis as a product of mountain building

BACKGROUND

Incarvillea sinensis is widely distributed from Southwest China to Northeast China and in the Russian Far East. The distribution of this species was thought to be influenced by the uplift of the Qinghai-Tibet Plateau and Quaternary glaciation. To reveal the imprints of geological events on the spatial genetic structure of Incarvillea sinensis, we examined two cpDNA segments ( trnH- psbA and trnS- trnfM) in 705 individuals from 47 localities.

RESULTS

A total of 16 haplotypes was identified, and significant genetic differentiation was revealed (GST =0.843, NST = 0.975, P < 0.05). The survey detected two highly divergent cpDNA lineages connected by a deep gap with allopatric distributions: the southern lineage with higher genetic diversity and differentiation in the eastern Qinghai-Tibet Plateau, and the northern lineage in the region outside the Qinghai-Tibet Plateau. The divergence between these two lineages was estimated at 4.4 MYA. A correlation between the genetic and the geographic distances indicates that genetic drift was more influential than gene flow in the northern clade with lower diversity and divergence. However, a scenario of regional equilibrium between gene flow and drift was shown for the southern clade. The feature of spatial distribution of the genetic diversity of the southern lineage possibly indicated that allopatric fragmentation was dominant in the collections from the eastern Qinghai-Tibet Plateau.

CONCLUSIONS

The results revealed that the uplift of the Qinghai-Tibet Plateau likely resulted in the significant divergence between the lineage in the eastern Qinghai-Tibet Plateau and the other one outside this area. The diverse niches in the eastern Qinghai-Tibet Plateau created a wide spectrum of habitats to accumulate and accommodate new mutations. The features of genetic diversity of populations outside the eastern Qinghai-Tibet Plateau seemed to reveal the imprints of extinction during the Glacial and the interglacial and postglacial recolonization. Our study is a typical case of the significance of the uplift of the Qinghai-Tibet Plateau and the Quaternary Glacial in spatial genetic structure of eastern Asian plants, and sheds new light on the evolution of biodiversity in the Qinghai-Tibet Plateau at the intraspecies level.

Evaluation of the 16S and 12S rRNA genes as universal markers for the identification of commercial fish species in South Africa

The development of DNA-based methods for the identification of fish species is important for fisheries research and control, as well as for the detection of unintentional or fraudulent species substitutions in the marketplace. The aim of this study was to generate a comprehensive reference database of DNA sequences from the mitochondrial 16S and 12S ribosomal RNA (rRNA) genes for 53 commercial fish species in South Africa and to evaluate the applicability of these genetic markers for the identification of fish at the species level. The DNA extracted from all target species was readily amplified using universal primers targeting both rRNA gene regions. Sequences from the 16S and 12S rRNA genes were submitted to GenBank for the first time for 34% and 53% of the fish species, respectively. Cumulative analysis of the 16S rRNA gene sequences revealed mean conspecific, congeneric and confamilial Kimura two parameter (K2P) distances of 0.03%, 0.70% and 5.10% and the corresponding values at the 12S level were 0.03%, 1.00% and 5.57%. K2P neighbour-joining trees based on both sequence datasets generally clustered species in accordance with their taxonomic classifications. The nucleotide variation in both the 16S and 12S sequences was suitable for identifying the large majority of the examined fish specimens to at least the level of genus, but was found to be less useful for the explicit differentiation of certain congeneric fish species. It is recommended that one or more faster-evolving DNA regions be analysed to confirm the identities of closely-related fish species in South Africa.

Plant molecular phylogeography in China and adjacent regions: Tracing the genetic imprints of Quaternary climate and environmental change in the world's most diverse temperate flora

The Sino-Japanese Floristic Region (SJFR) of East Asia harbors the most diverse of the world's temperate flora, and was the most important glacial refuge for its Tertiary representatives ('relics') throughout Quaternary ice-age cycles. A steadily increasing number of phylogeographic studies in the SJFR of mainland China and adjacent areas, including the Qinghai-Tibetan-Plateau (QTP) and Sino-Himalayan region, have documented the population histories of temperate plant species in these regions. Here we review this current literature that challenges the oft-stated view of the SJFR as a glacial sanctuary for temperate plants, instead revealing profound effects of Quaternary changes in climate, topography, and/or sea level on the current genetic structure of such organisms. There are three recurrent phylogeographic scenarios identified by different case studies that broadly agree with longstanding biogeographic or palaeo-ecological hypotheses: (i) postglacial re-colonization of the QTP from (south-)eastern glacial refugia; (ii) population isolation and endemic species formation in Southwest China due to tectonic shifts and river course dynamics; and (iii) long-term isolation and species survival in multiple localized refugia of (warm-)temperate deciduous forest habitats in subtropical (Central/East/South) China. However, in four additional instances, phylogeographic findings seem to conflict with a priori predictions raised by palaeo-data, suggesting instead: (iv) glacial in situ survival of some hardy alpine herbs and forest trees on the QTP platform itself; (v) long-term refugial isolation of (warm-)temperate evergreen taxa in subtropical China; (vi) 'cryptic' glacial survival of (cool-)temperate deciduous forest trees in North China; and (vii) unexpectedly deep (Late Tertiary/early-to-mid Pleistocene) allopatric-vicariant differentiation of disjunct lineages in the East China-Japan-Korea region due to past sea transgressions. We discuss these and other consequences of the main phylogeographic findings in light of palaeo-environmental evidence, emphasize notable gaps in our knowledge, and outline future research prospects for disentangling the evolution and biogeographic history of the region's extremely diverse temperate flora.

An overview to the investigative approach to species testing in wildlife forensic science

The extent of wildlife crime is unknown but it is on the increase and has observable effects with the dramatic decline in many species of flora and fauna. The growing awareness of this area of criminal activity is reflected in the increase in research papers on animal DNA testing, either for the identification of species or for the genetic linkage of a sample to a particular organism. This review focuses on the use of species testing in wildlife crime investigations. Species identification relies primarily on genetic loci within the mitochondrial genome; focusing on the cytochrome b and cytochrome oxidase 1 genes. The use of cytochrome b gained early prominence in species identification through its use in taxonomic and phylogenetic studies, while the gene sequence for cytochrome oxidase was adopted by the Barcode for Life research group. This review compares how these two loci are used in species identification with respect to wildlife crime investigations. As more forensic science laboratories undertake work in the wildlife area, it is important that the quality of work is of the highest standard and that the conclusions reached are based on scientific principles. A key issue in reporting on the identification of a particular species is a knowledge of both the intraspecies variation and the possible overlap of sequence variation from one species to that of a closely related species. Recent data showing this degree of genetic separation in mammalian species will allow greater confidence when preparing a report on an alleged event where the identification of the species is of prime importance. The aim of this review is to illustrate aspects of species testing in wildlife forensic science and to explain how a knowledge of genetic variation at the genus and species level can aid in the reporting of results.

Reconstructing mammalian phylogenies: a detailed comparison of the cytochrome B and cytochrome oxidase subunit I mitochondrial genes

The phylogeny and taxonomy of mammalian species were originally based upon shared or derived morphological characteristics. However, genetic analyses have more recently played an increasingly important role in confirming existing or establishing often radically different mammalian groupings and phylogenies. The two most commonly used genetic loci in species identification are the cytochrome oxidase I gene (COI) and the cytochrome b gene (cyt b). For the first time this study provides a detailed comparison of the effectiveness of these two loci in reconstructing the phylogeny of mammals at different levels of the taxonomic hierarchy in order to provide a basis for standardizing methodologies in the future. Interspecific and intraspecific variation is assessed and for the first time, to our knowledge, statistical confidence is applied to sequence comparisons. Comparison of the DNA sequences of 217 mammalian species reveals that cyt b more accurately reconstructs their phylogeny and known relationships between species based on other molecular and morphological analyses at Super Order, Order, Family and generic levels. Cyt b correctly assigned 95.85% of mammal species to Super Order, 94.31% to Order and 98.16% to Family compared to 78.34%, 93.36% and 96.93% respectively for COI. Cyt b also gives better resolution when separating species based on sequence data. Using a Kimura 2-parameter p-distance (x100) threshold of 1.5-2.5, cyt b gives a better resolution for separating species with a lower false positive rate and higher positive predictive value than those of COI.

DNA typing in wildlife crime: recent developments in species identification

Species identification has become a tool in the investigation of acts of alleged wildlife crimes. This review details the steps required in DNA testing in wildlife crime investigations and highlights recent developments where not only can individual species be identified within a mixture of species but multiple species can be identified simultaneously. 'What species is this?' is a question asked frequently in wildlife crime investigations. Depending on the material being examined, DNA analysis may offer the best opportunity to answer this question. Species testing requires the comparison of the DNA type from the unknown sample to DNA types on a database. The areas of DNA tested are on the mitochondria and include predominantly the cytochrome b gene and the cytochrome oxidase I gene. Standard analysis requires the sequencing of part of one of these genes and comparing the sequence to that held on a repository of DNA sequences such as the GenBank database. Much of the DNA sequence of either of these two genes is conserved with only parts being variable. A recent development is to target areas of those sequences that are specific to a species; this can increase the sensitivity of the test with no loss of specificity. The benefit of targeting species specific sequences is that within a mixture of two of more species, the individual species within the mixture can be identified. This identification would not be possible using standard sequencing. These new developments can lead to a greater number of samples being tested in alleged wildlife crimes.

Molecular identification of Indian crocodile species: PCR-RFLP method for forensic authentication*

South East Asian countries are known for illegal poaching and trade of crocodiles clandestinely, to be used in skin, medicinal, and cosmetic industries. Besides crocodiles being listed in the Convention on International Trade in Endangered Species of Wild Fauna and Flora, India has its Wildlife Protection Act, 1972 for conservation of crocodile species. Hitherto, lack of any rapid and reliable technique for examinations of crocodile-based crime exhibits such as skin, bones, etc. has been a major problem for an effective promulgation of law on illegal trade. DNA-based identification of species using PCR-RFLP technique for an apt identification of all the three Indian crocodile species namely, Crocodylus porosus, Crocodylus palustris and Gavialis gangeticus is presented here. A 628 bp segment of cytochrome b gene was amplified using novel primers followed by restriction digestion with three enzymes i.e., HaeIII, MboI, and MwoI, separately and in combination. The technique has produced a species-specific pattern for identifying the three crocodile species individually, which fulfills the requirement for its forensic application. It is expected that the technique will prove handy in identification of all the three Indian crocodile species and strengthen conservation efforts.

Cross-species chromosome painting in Cetartiodactyla: reconstructing the karyotype evolution in key phylogenetic lineages

Recent molecular and morphological studies place Artiodactyla and Cetacea into the order Cetartiodactyla. Within the Cetartiodactyla such families as Bovidae, Cervidae, and Suidae are well studied by comparative chromosome painting, but many taxa that are crucial for understanding cetartiodactyl phylogeny remain poorly studied. Here we present the genome-wide comparative maps of five cetartiodactyl species obtained by chromosome painting with human and dromedary paint probes from four taxa: Cetacea, Hippopotamidae, Giraffidae, and Moschidae. This is the first molecular cytogenetic report on pilot whale, hippopotamus, okapi, and Siberian musk deer. Our results, when integrated with previously published comparative chromosome maps allow us to reconstruct the evolutionary pathway and rates of chromosomal rearrangements in Cetartiodactyla. We hypothesize that the putative cetartiodactyl ancestral karyotype (CAK) contained 25-26 pairs of autosomes, 2n = 52-54, and that the association of human chromosomes 8/9 could be a cytogenetic signature that unites non-camelid cetartiodactyls. There are no unambiguous cytogenetic landmarks that unite Hippopotamidae and Cetacea. If we superimpose chromosome rearrangements on the supertree generated by Price and colleagues, several homoplasy events are needed to explain cetartiodactyl karyotype evolution. Our results apparently favour a model of non-random breakpoints in chromosome evolution. Cetariodactyl karyotype evolution is characterized by alternating periods of low and fast rates in various lineages. The highest rates are found in Suina (Suidae+Tayasuidae) lineage (1.76 rearrangements per million years (R/My)) and the lowest in Cetaceans (0.07 R/My). Our study demonstrates that the combined use of human and camel paints is highly informative for revealing evolutionary karyotypic rearrangements among cetartiodactyl species.

Genetic diversity of captive forest musk deer (Moschus berezovskii) inferred from the mitochondrial DNA control region

Forest musk deer (Moschus berezovskii) were once distributed widely in China. However, wild populations have declined dramatically because of poaching and habitat loss. Captive breeding populations have been established for several decades, but the genetic backgrounds of most captive populations were unclear and the population sizes increased very slowly. To provide useful information for conservation and management of this species, we investigated the genetic diversity and population structure of forest musk deer by analysing a 582-bp fragment of the mitochondrial DNA (mtDNA) control region (CR) in three captive breeding populations in Sichuan Province, China. Ninety-four variable sites and 27 haplotypes were observed in 109 individuals, and the nucleotide and haplotype diversities were relatively high compared with those of other endangered mammals. Of the three investigated populations, the Maerkang population had the highest nucleotide diversity (pi=0.0568), haplotype diversity (h=0.836) and average intra-population genetic distance (0.062). The analysis of molecular variance demonstrated that most variation occurred within samples and that there was significant differentiation of the three populations. Estimates of gene flow indicated that there were few genetic exchanges among the three populations. Building pedigree records and increasing gene flow between populations will be helpful for conserving these populations and this species.

Species identification using the cytochrome b gene of commercial turtle shells

Turtle shells and their gelled products are familiar in some countries as foods, tonics and medicines. These shells may come from endangered and protected species, requiring the identification of the species present to enforce national and international legislation. We report on the design of five combinations of primer pairs for the identification of turtle shells and shell fragments used as ornaments, food products and medicines. The types of samples used are those encountered frequently and will typically contain highly degraded DNA. The success rate for species identification using the test described is dependent upon the choice of primer sets used and the length of the expected amplification product. Gelled products were simulated by the process of decoction for up to 12 h, after which all the turtle species could be identified from the liquid samples. This study establishes a method for the identification of commercial turtle shells and illustrates a simulated case using gelled products.

Molecular evolution study in China: progress and future promise

China has a large land area with highly diverse topography, climate and vegetation, and animal resources and is ranked eighth in the world and first in the Northern Hemisphere on richness of biodiversity. Even though little work on molecular evolution had been reported a decade ago, studies on both the evolution of macromolecules and the molecular phylogeny have become active in China in recent years. This review highlights some of the interesting and important developments in molecular evolution study in China. Chinese scientists have made significant contribution on the methods inferring phylogeny and biogeography of animals and plants in East Asia using molecular data. Studies on population and conservation genetics of animals and plants, such as Golden monkey and Chinese sturgeon, provided useful information for conserving the endangered species. East and South Asia has been demonstrated to be one of the centres of domestication. Origin and evolution of genes and gene families have been explored, which shed new insight on the genetic mechanism of adaptation. In the genomic era, Chinese researchers also made a transition from single-gene to a genomic investigation approach. Considering the fact that amazing progress has been made in the past few years, and more and more talented young scientists are entering field, the future of molecular evolution study in China holds much promise.

Molecular phylogeny of musk deer: a genomic view with mitochondrial 16S rRNA and cytochrome b gene

The phylogenetic status of the infra order Pecora is controversial, even though it is supported by paleontological, morphological, and molecular evidence. We analyzed two mitochondrial genes (i.e., 16S rRNA and cytochrome b) to resolve the phylogenetic position of pecoran species, i.e., the Bovidae, Cervidae, and Moschidae endemic to the Indian subcontinent. We used phylogenetic analysis based on different algorithms, including neighbor joining, maximum parsimony, Bayesian inference, maximum likelihood, minimum evolution, median joining network, along with multidimensional scaling, and DNA word analysis. Our results established the basal position of Tragulidae and the monophyly of the infra order Pecora within the Suborder Ruminantia. Our results also demonstrated that Bovidae, Cervidae, and Moschidae are allied with the placement of musk deer as more closely related to bovids than to cervids. Molecular dating based on sequence analysis shows that the radiation of Pecora occurred during the early Oligocene and that the majority of the pecoran families radiated and dispersed rapidly during the Oligocene/Miocene transition.

Critical factors for assembling a high volume of DNA barcodes

Large-scale DNA barcoding projects are now moving toward activation while the creation of a comprehensive barcode library for eukaryotes will ultimately require the acquisition of some 100 million barcodes. To satisfy this need, analytical facilities must adopt protocols that can support the rapid, cost-effective assembly of barcodes. In this paper we discuss the prospects for establishing high volume DNA barcoding facilities by evaluating key steps in the analytical chain from specimens to barcodes. Alliances with members of the taxonomic community represent the most effective strategy for provisioning the analytical chain with specimens. The optimal protocols for DNA extraction and subsequent PCR amplification of the barcode region depend strongly on their condition, but production targets of 100K barcode records per year are now feasible for facilities working with compliant specimens. The analysis of museum collections is currently challenging, but PCR cocktails that combine polymerases with repair enzyme(s) promise future success. Barcode analysis is already a cost-effective option for species identification in some situations and this will increasingly be the case as reference libraries are assembled and analytical protocols are simplified.

A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla)

Despite the biological and economic importance of the Cetartiodactyla, the phylogeny of this clade remains controversial. Using the supertree approach of matrix representation with parsimony, we present the first phylogeny to include all 290 extant species of the Cetacea (whales and dolphins) and Artiodactyla (even-toed hoofed mammals). At the family-level, the supertree is fully resolved. For example, the relationships among the Ruminantia appear as (((Cervidae, Moschidae) Bovidae) (Giraffidae, Antilocapridae) Tragulidae). However, due to either lack of phylogenetic study or contradictory information, polytomies occur within the clades Sus, Muntiacus, Cervus, Delphinidae, Ziphiidae and Bovidae. Complete species-level phylogenies are necessary for both illustrating and analysing biological, geographical and ecological patterns in an evolutionary framework. The present species-level tree of the Cetartiodactyla provides the first opportunity to examine comparative hypotheses across entirely aquatic and terrestrial species within a single mammalian order.

A complete estimate of the phylogenetic relationships in Ruminantia: a dated species-level supertree of the extant ruminants

This paper presents the first complete estimate of the phylogenetic relationships among all 197 species of extant and recently extinct ruminants combining morphological, ethological and molecular information. The composite tree is derived by applying matrix representation using parsimony analysis to 164 previous partial estimates, and is remarkably well resolved, containing 159 nodes (> 80 % of the potential nodes in the completely resolved phylogeny). Bremer decay index has been used to indicate the degree of certainty associated with each clade. The ages of over 80% of the clades in the tree have been estimated from information in the literature. The supertree for Ruminantia illustrates which areas of ruminant phylogeny are still only roughly known because of taxa with controversial relationships (e.g. Odocoileini, Antilopinae) or not studied in great detail (e.g. Muntiacus). It supports the monophyly of the ruminant families and Pecora. According to this analysis Antilocapridae and Giraffidae constitute the superfamily Giraffoidea, which is the sister group of a clade clustering Bovoidea and Cervoidea. The position of several taxa whose systematic positions have remained controversial in the past (Saiga, Pelea, Aepycerus, Pantholops, Ammotragus, Pseudois) is unambiguously established. Nevertheless, the position of Neotragus and Oreotragus within the original radiation of the non-bovine bovids remains unresolved in the present analysis. It also shows that six successive rapid cladogenesis events occurred within the infraorder Pecora during the Oligocene to middle Pliocene, which coincided with periods of global climatic change. Finally, the presented supertree will be a useful framework for comparative and evolutionary biologists interested in studies involving the ruminants.

New insights into the karyotypic relationships of Chinese muntjac (Muntiacus reevesi), forest musk deer (Moschus berezovskii) and gayal (Bos frontalis)

To investigate the karyotypic relationships between Chinese muntjac (Muntiacus reevesi), forest musk deer (Moschus berezovskii) and gayal (Bos frontalis), a complete set of Chinese muntjac chromosome-specific painting probes has been assigned to G-banded chromosomes of these three species. Sixteen autosomal probes (i.e. 6-10, 12-22) of the Chinese muntjac each delineated one pair of conserved segments in the forest musk deer and gayal, respectively. The remaining six autosomal probes (1-5, and 11) each delineated two to five pairs of conserved segments. In total, the 22 autosomal painting probes of Chinese muntjac delineated 33 and 34 conserved chromosomal segments in the genomes of forest musk deer and gayal, respectively. The combined analysis of comparative chromosome painting and G-band comparison reveals that most interspecific homologous segments show a high degree of conservation in G-banding patterns. Eleven chromosome fissions and five chromosome fusions differentiate the karyotypes of Chinese muntjac and forest musk deer; twelve chromosome fissions and six fusions are required to convert the Chinese muntjac karyotype to that of gayal; one chromosome fission and one fusion separate the forest musk deer and gayal. The musk deer has retained a highly conserved karyotype that closely resembles the proposed ancestral pecoran karyotype but shares none of the rearrangements characteristic for the Cervidae and Bovidae. Our results substantiate that chromosomes 1-5 and 11 of Chinese muntjac originated through exclusive centromere-to-telomere fusions of ancestral acrocentric chromosomes.

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.

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.

Rapid and parallel chromosomal number reductions in muntjac deer inferred from mitochondrial DNA phylogeny

Muntjac deer (Muntiacinae, Cervidae) are of great interest in evolutionary studies because of their dramatic chromosome variations and recent discoveries of several new species. In this paper, we analyze the evolution of karyotypes of muntjac deer in the context of a phylogeny which is based on 1,844-bp mitochondrial DNA sequences of seven generally recognized species in the muntjac subfamily. The phylogenetic results support the hypothesis that karyotypic evolution in muntjac deer has proceeded via reduction in diploid number. However, the reduction in number is not always linear, i.e., not strictly following the order: 46-->14/13-->8/9-->6/7. For example, Muntiacus muntjak (2n = 6/7) shares a common ancestor with Muntiacus feae (2n = 13/14), which indicates that its karyotype was derived in parallel with M. feae's from an ancestral karyotype of 2n >/= 13/14. The newly discovered giant muntjac (Muntiacus vuquangensis) may represent another parallel reduction lineage from the ancestral 2n = 46 karyotype. Our phylogenetic results indicate that the giant muntjac is relatively closer to Muntiacus reevesi than to other muntjacs and may be placed in the genus Muntiacus Analyses of sequence divergence reveal that the rate of change in chromosome number in muntjac deer is one of the fastest in vertebrates. Within the muntjac subfamily, the fastest evolutionary rate is found in the Fea's lineage, in which two species with different karyotypes diverged in around 0.5 Myr.