Mitochondrial genomes reveal an explosive radiation of extinct and extant bears near the Miocene-Pliocene boundary

Genomic evidence for island population conversion resolves conflicting theories of polar bear evolution

Despite extensive genetic analysis, the evolutionary relationship between polar bears (Ursus maritimus) and brown bears (U. arctos) remains unclear. The two most recent comprehensive reports indicate a recent divergence with little subsequent admixture or a much more ancient divergence followed by extensive admixture. At the center of this controversy are the Alaskan ABC Islands brown bears that show evidence of shared ancestry with polar bears. We present an analysis of genome-wide sequence data for seven polar bears, one ABC Islands brown bear, one mainland Alaskan brown bear, and a black bear (U. americanus), plus recently published datasets from other bears. Surprisingly, we find clear evidence for gene flow from polar bears into ABC Islands brown bears but no evidence of gene flow from brown bears into polar bears. Importantly, while polar bears contributed <1% of the autosomal genome of the ABC Islands brown bear, they contributed 6.5% of the X chromosome. The magnitude of sex-biased polar bear ancestry and the clear direction of gene flow suggest a model wherein the enigmatic ABC Island brown bears are the descendants of a polar bear population that was gradually converted into brown bears via male-dominated brown bear admixture. We present a model that reconciles heretofore conflicting genetic observations. We posit that the enigmatic ABC Islands brown bears derive from a population of polar bears likely stranded by the receding ice at the end of the last glacial period. Since then, male brown bear migration onto the island has gradually converted these bears into an admixed population whose phenotype and genotype are principally brown bear, except at mtDNA and X-linked loci. This process of genome erosion and conversion may be a common outcome when climate change or other forces cause a population to become isolated and then overrun by species with which it can hybridize.

MHC class II DQB diversity in the Japanese black bear, Ursus thibetanus japonicus

Background

The major histocompatibility complex (MHC) genes are one of the most important genetic systems in the vertebrate immune response. The diversity of MHC genes may directly influence the survival of individuals against infectious disease. However, there has been no investigation of MHC diversity in the Asiatic black bear (Ursus thibetanus). Here, we analyzed 270-bp nucleotide sequences of the entire exon 2 region of the MHC DQB gene by using 188 samples from the Japanese black bear (Ursus thibetanus japonicus) from 12 local populations.

Results

Among 185 of 188 samples, we identified 44 MHC variants that encoded 31 different amino acid sequences (allotypes) and one putative pseudogene. The phylogenetic analysis suggests that MHC variants detected from the Japanese black bear are derived from the DQB locus. One of the 31 DQB allotypes, Urth-DQB*01, was found to be common to all local populations. Moreover, this allotype was shared between the black bear on the Asian continent and the Japanese black bear, suggesting that Urth-DQB*01 might have been maintained in the ancestral black bear population for at least 300,000 years. Our findings, from calculating the ratio of non-synonymous to synonymous substitutions, indicate that balancing selection has maintained genetic variation of peptide-binding residues at the DQB locus of the Japanese black bear. From examination of genotype frequencies among local populations, we observed a considerably lower level of observed heterozygosity than expected.

Conclusions

The low level of observed heterozygosity suggests that genetic drift reduced DQB diversity in the Japanese black bear due to a bottleneck event at the population or species level. The decline of DQB diversity might have been accelerated by the loss of rare variants that have been maintained by negative frequency-dependent selection. Nevertheless, DQB diversity of the black bear appears to be relatively high compared with some other endangered mammalian species. This result suggests that the Japanese black bears may also retain more potential resistance against pathogens than other endangered mammalian species. To prevent further decline of potential resistance against pathogens, a conservation policy for the Japanese black bear should be designed to maintain MHC rare variants in each local population.

Kretzoiarctos gen. nov., the oldest member of the giant panda clade

The phylogenetic position of the giant panda, Ailuropoda melanoleuca (Carnivora: Ursidae: Ailuropodinae), has been one of the most hotly debated topics by mammalian biologists and paleontologists during the last century. Based on molecular data, it is currently recognized as a true ursid, sister-taxon of the remaining extant bears, from which it would have diverged by the Early Miocene. However, from a paleobiogeographic and chronological perspective, the origin of the giant panda lineage has remained elusive due to the scarcity of the available Miocene fossil record. Until recently, the genus Ailurarctos from the Late Miocene of China (ca. 8–7 mya) was recognized as the oldest undoubted member of the Ailuropodinae, suggesting that the panda lineage might have originated from an Ursavus ancestor. The role of the purported ailuropodine Agriarctos, from the Miocene of Europe, in the origins of this clade has been generally dismissed due to the paucity of the available material. Here, we describe a new ailuropodine genus, Kretzoiarctos gen. nov., based on remains from two Middle Miocene (ca. 12–11 Ma) Spanish localities. A cladistic analysis of fossil and extant members of the Ursoidea confirms the inclusion of the new genus into the Ailuropodinae. Moreover, Kretzoiarctos precedes in time the previously-known, Late Miocene members of the giant panda clade from Eurasia (Agriarctos and Ailurarctos). The former can be therefore considered the oldest recorded member of the giant panda lineage, which has significant implications for understanding the origins of this clade from a paleobiogeographic viewpoint.

Paleopopulation genetics

Paleopopulation genetics is a new field that focuses on the population genetics of extinct groups and ancestral populations (i.e., populations ancestral to extant groups). With recent advances in DNA sequencing technologies, we now have unprecedented ability to directly assay genetic variation from fossils. This allows us to address issues, such as past population structure, changes in population size, and evolutionary relationships between taxa, at a much greater resolution than can traditional population genetics studies. In this review, we discuss recent developments in this emerging field as well as prospects for the future.

Extra-Mediterranean refugia: The rule and not the exception?

Some decades ago, biogeographers distinguished three major faunal types of high importance for Europe: (i) Mediterranean elements with exclusive glacial survival in the Mediterranean refugia, (ii) Siberian elements with glacial refugia in the eastern Palearctic and only postglacial expansion to Europe and (iii) arctic and/or alpine elements with large zonal distributions in the periglacial areas and postglacial retreat to the North and/or into the high mountain systems. Genetic analyses have unravelled numerous additional refugia both of continental and Mediterranean species, thus strongly modifying the biogeographical view of Europe. This modified notion is particularly true for the so-called Siberian species, which in many cases have not immigrated into Europe during the postglacial period, but most likely have survived the last, or even several glacial phases, in extra-Mediterranean refugia in some climatically favourable but geographically limited areas of southern Central and Eastern Europe. Recently, genetic analyses revealed that typical Mediterranean species have also survived the Last Glacial Maximum in cryptic northern refugia (e.g. in the Carpathians or even north of the Alps) in addition to their Mediterranean refuge areas.

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

Background

The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs).

Results

To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are monophyletic in the mtgenome phylogeny, respectively, but a H. pileatus individual falls into the H. lar Y chromosome clade.

Conclusions

Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.

Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaska's Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5-10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4-5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.

Functional annotation from the genome sequence of the giant panda

The giant panda is one of the most critically endangered species due to the fragmentation and loss of its habitat. Studying the functions of proteins in this animal, especially specific trait-related proteins, is therefore necessary to protect the species. In this work, the functions of these proteins were investigated using the genome sequence of the giant panda. Data on 21,001 proteins and their functions were stored in the Giant Panda Protein Database, in which the proteins were divided into two groups: 20,179 proteins whose functions can be predicted by GeneScan formed the known-function group, whereas 822 proteins whose functions cannot be predicted by GeneScan comprised the unknown-function group. For the known-function group, we further classified the proteins by molecular function, biological process, cellular component, and tissue specificity. For the unknown-function group, we developed a strategy in which the proteins were filtered by cross-Blast to identify panda-specific proteins under the assumption that proteins related to the panda-specific traits in the unknown-function group exist. After this filtering procedure, we identified 32 proteins (2 of which are membrane proteins) specific to the giant panda genome as compared against the dog and horse genomes. Based on their amino acid sequences, these 32 proteins were further analyzed by functional classification using SVM-Prot, motif prediction using MyHits, and interacting protein prediction using the Database of Interacting Proteins. Nineteen proteins were predicted to be zinc-binding proteins, thus affecting the activities of nucleic acids. The 32 panda-specific proteins will be further investigated by structural and functional analysis.

Ancient DNA studies: new perspectives on old samples

In spite of past controversies, the field of ancient DNA is now a reliable research area due to recent methodological improvements. A series of recent large-scale studies have revealed the true potential of ancient DNA samples to study the processes of evolution and to test models and assumptions commonly used to reconstruct patterns of evolution and to analyze population genetics and palaeoecological changes. Recent advances in DNA technologies, such as next-generation sequencing make it possible to recover DNA information from archaeological and paleontological remains allowing us to go back in time and study the genetic relationships between extinct organisms and their contemporary relatives. With the next-generation sequencing methodologies, DNA sequences can be retrieved even from samples (for example human remains) for which the technical pitfalls of classical methodologies required stringent criteria to guaranty the reliability of the results. In this paper, we review the methodologies applied to ancient DNA analysis and the perspectives that next-generation sequencing applications provide in this field.

Resurrecting ancient animal genomes: the extinct moa and more

Recently two developments have had a major impact on the field of ancient DNA (aDNA). First, new advances in DNA sequencing, in combination with improved capture/enrichment methods, have resulted in the recovery of orders of magnitude more DNA sequence data from ancient animals. Second, there has been an increase in the range of tissue types employed in aDNA. Hair in particular has proven to be very successful as a source of DNA because of its low levels of contamination and high level of ancient endogenous DNA. These developments have resulted in significant advances in our understanding of recently extinct animals: namely their evolutionary relationships, physiology, and even behaviour. Hair has been used to recover the first complete ancient nuclear genome, that of the extinct woolly mammoth, which then facilitated the expression and functional analysis of haemoglobins. Finally, we speculate on the consequences of these developments for the possibility of recreating extinct animals.

Multiplex PCR amplification of ancient DNA

Multiplex PCR allows the simultaneous amplification of up to dozens of target fragments in a single PCR. It is therefore a powerful tool to obtain many kilobases of continuous sequence from minute amounts of ancient DNA (aDNA), which usually must be amplified in multiple short and overlapping fragments. Because significantly less template is required compared to amplifying each fragment separately, multiplex PCR is particularly beneficial when the fossil material itself, or access to the fossil material, is limited. The recently refined two-step multiplex PCR protocol consists of a first-step reaction (the actual multiplex PCR) that then acts as the template for the second-step PCR. During the second step, nested primers are used in individual amplification reactions. Although the same set of primers can be used in both steps, using a nested set in the second step adds an additional level of selectivity and specificity, minimizing PCR artifacts. This is particularly important when complex mixtures of template DNA, such as aDNA extracts, are amplified.

Setting up an ancient DNA laboratory

Entering into the world of ancient DNA research is nontrivial. Because the DNA in most ancient specimens is degraded to some extent, the potential for contamination of ancient samples and DNA extracts with modern DNA is considerable. To minimize the risk associated with working with ancient DNA, experimental protocols specific to handling ancient specimens have been introduced. Here, I outline the challenges associated with working with ancient DNA and describe guidelines for setting up a new ancient DNA laboratory. I also discuss steps that can be taken at the sample collection and preparation stage to minimize the potential for contamination with exogenous sources of DNA.

Targeted enrichment of ancient pathogens yielding the pPCP1 plasmid of Yersinia pestis from victims of the Black Death

Although investigations of medieval plague victims have identified Yersinia pestis as the putative etiologic agent of the pandemic, methodological limitations have prevented large-scale genomic investigations to evaluate changes in the pathogen's virulence over time. We screened over 100 skeletal remains from Black Death victims of the East Smithfield mass burial site (1348-1350, London, England). Recent methods of DNA enrichment coupled with high-throughput DNA sequencing subsequently permitted reconstruction of ten full human mitochondrial genomes (16 kb each) and the full pPCP1 (9.6 kb) virulence-associated plasmid at high coverage. Comparisons of molecular damage profiles between endogenous human and Y. pestis DNA confirmed its authenticity as an ancient pathogen, thus representing the longest contiguous genomic sequence for an ancient pathogen to date. Comparison of our reconstructed plasmid against modern Y. pestis shows identity with several isolates matching the Medievalis biovar; however, our chromosomal sequences indicate the victims were infected with a Y. pestis variant that has not been previously reported. Our data reveal that the Black Death in medieval Europe was caused by a variant of Y. pestis that may no longer exist, and genetic data carried on its pPCP1 plasmid were not responsible for the purported epidemiological differences between ancient and modern forms of Y. pestis infections.

The complete mitochondrial genome of an 11,450-year-old aurochsen (Bos primigenius) from Central Italy

Background

Bos primigenius, the aurochs, is the wild ancestor of modern cattle breeds and was formerly widespread across Eurasia and northern Africa. After a progressive decline, the species became extinct in 1627. The origin of modern taurine breeds in Europe is debated. Archaeological and early genetic evidence point to a single Near Eastern origin and a subsequent spread during the diffusion of herding and farming. More recent genetic data are instead compatible with local domestication events or at least some level of local introgression from the aurochs. Here we present the analysis of the complete mitochondrial genome of a pre-Neolithic Italian aurochs.

Results

In this study, we applied a combined strategy employing both multiplex PCR amplifications and 454 pyrosequencing technology to sequence the complete mitochondrial genome of an 11,450-year-old aurochs specimen from Central Italy. Phylogenetic analysis of the aurochs mtDNA genome supports the conclusions from previous studies of short mtDNA fragments - namely that Italian aurochsen were genetically very similar to modern cattle breeds, but highly divergent from the North-Central European aurochsen.

Conclusions

Complete mitochondrial genome sequences are now available for several modern cattle and two pre-Neolithic mtDNA genomes from very different geographic areas. These data suggest that previously identified sub-groups within the widespread modern cattle mitochondrial T clade are polyphyletic, and they support the hypothesis that modern European breeds have multiple geographic origins.

Potential efficacy of mitochondrial genes for animal DNA barcoding: a case study using eutherian mammals

BACKGROUND

A well-informed choice of genetic locus is central to the efficacy of DNA barcoding. Current DNA barcoding in animals involves the use of the 5' half of the mitochondrial cytochrome oxidase 1 gene (CO1) to diagnose and delimit species. However, there is no compelling a priori reason for the exclusive focus on this region, and it has been shown that it performs poorly for certain animal groups. To explore alternative mitochondrial barcoding regions, we compared the efficacy of the universal CO1 barcoding region with the other mitochondrial protein-coding genes in eutherian mammals. Four criteria were used for this comparison: the number of recovered species, sequence variability within and between species, resolution to taxonomic levels above that of species, and the degree of mutational saturation.

RESULTS

Based on 1,179 mitochondrial genomes of eutherians, we found that the universal CO1 barcoding region is a good representative of mitochondrial genes as a whole because the high species-recovery rate (> 90%) was similar to that of other mitochondrial genes, and there were no significant differences in intra- or interspecific variability among genes. However, an overlap between intra- and interspecific variability was still problematic for all mitochondrial genes. Our results also demonstrated that any choice of mitochondrial gene for DNA barcoding failed to offer significant resolution at higher taxonomic levels.

CONCLUSIONS

We suggest that the CO1 barcoding region, the universal DNA barcode, is preferred among the mitochondrial protein-coding genes as a molecular diagnostic at least for eutherian species identification. Nevertheless, DNA barcoding with this marker may still be problematic for certain eutherian taxa and our approach can be used to test potential barcoding loci for such groups.

The European Nucleotide Archive

The European Nucleotide Archive (ENA; http://www.ebi.ac.uk/ena) is Europe's primary nucleotide-sequence repository. The ENA consists of three main databases: the Sequence Read Archive (SRA), the Trace Archive and EMBL-Bank. The objective of ENA is to support and promote the use of nucleotide sequencing as an experimental research platform by providing data submission, archive, search and download services. In this article, we outline these services and describe major changes and improvements introduced during 2010. These include extended EMBL-Bank and SRA-data submission services, extended ENA Browser functionality, support for submitting data to the European Genome-phenome Archive (EGA) through SRA, and the launch of a new sequence similarity search service.

Mitochondrial genome sequences effectively reveal the phylogeny of Hylobates gibbons

Background

Uniquely among hominoids, gibbons exist as multiple geographically contiguous taxa exhibiting distinctive behavioral, morphological, and karyotypic characteristics. However, our understanding of the evolutionary relationships of the various gibbons, especially among Hylobates species, is still limited because previous studies used limited taxon sampling or short mitochondrial DNA (mtDNA) sequences. Here we use mtDNA genome sequences to reconstruct gibbon phylogenetic relationships and reveal the pattern and timing of divergence events in gibbon evolutionary history.

Methodology/Principal Findings

We sequenced the mitochondrial genomes of 51 individuals representing 11 species belonging to three genera (Hylobates, Nomascus and Symphalangus) using the high-throughput 454 sequencing system with the parallel tagged sequencing approach. Three phylogenetic analyses (maximum likelihood, Bayesian analysis and neighbor-joining) depicted the gibbon phylogenetic relationships congruently and with strong support values. Most notably, we recover a well-supported phylogeny of the Hylobates gibbons. The estimation of divergence times using Bayesian analysis with relaxed clock model suggests a much more rapid speciation process in Hylobates than in Nomascus.

Conclusions/Significance

Use of more than 15 kb sequences of the mitochondrial genome provided more informative and robust data than previous studies of short mitochondrial segments (e.g., control region or cytochrome b) as shown by the reliable reconstruction of divergence patterns among Hylobates gibbons. Moreover, molecular dating of the mitogenomic divergence times implied that biogeographic change during the last five million years may be a factor promoting the speciation of Sundaland animals, including Hylobates species.

Mitochondrial genome phylogeny among Asiatic black bear Ursus thibetanus subspecies and comprehensive analysis of their control regions

The complete mitochondrial genome (16,824 bp) of an Asiatic black bear Ursus thibetanus ussuricus (Mammalia, Carnivora, Ursidae) was newly sequenced and characterized in detail. It is the second mitochondrial genome from this subspecies which has been completely sequenced. The two U. t. ussuricus individuals were compared with each other and then with individuals from the other four U. thibetanus subspecies and the other nine ursid species, focusing especially on the control regions in the 14 mitochondrial genomes. Within these control regions, tandem repeats of basically 10 bp (5'-ACGCACGTGT-3' or its derivatives) were found in Domain II. Plausible secondary structures of the repeat region were compared between the North and South Korean individuals of U. t. ussuricus. According to the maximum likelihood and Bayesian inference trees inferred from the nucleotide sequences of 13 protein-coding and two rRNA genes, the ursine members within the monophyletic ursid clade can be divided into at least three groups: A, B, and C. According to this analysis, U. thibetanus subspecies were found with Ursus americanus and Ursus malayanus within Group A, showing the following relationships with nodal bootstrap values above 91% and Bayesian posterior probabilities of 1.00: ([(U. t. thibetanus, U. t. formosanus), U. t. spp.], U. t. ussuricus), U. t. mupinensis. In addition, we present a hypothetical scenario of the evolution of the major repeat motifs in the control region.

Biomechanical consequences of rapid evolution in the polar bear lineage

The polar bear is the only living ursid with a fully carnivorous diet. Despite a number of well-documented craniodental adaptations for a diet of seal flesh and blubber, molecular and paleontological data indicate that this morphologically distinct species evolved less than a million years ago from the omnivorous brown bear. To better understand the evolution of this dietary specialization, we used phylogenetic tests to estimate the rate of morphological specialization in polar bears. We then used finite element analysis (FEA) to compare the limits of feeding performance in the polar bear skull to that of the phylogenetically and geographically close brown bear. Results indicate that extremely rapid evolution of semi-aquatic adaptations and dietary specialization in the polar bear lineage produced a cranial morphology that is weaker than that of brown bears and less suited to processing tough omnivorous or herbivorous diets. Our results suggest that continuation of current climate trends could affect polar bears by not only eliminating their primary food source, but also through competition with northward advancing, generalized brown populations for resources that they are ill-equipped to utilize.

Next Generation Sequencing of Ancient DNA: Requirements, Strategies and Perspectives

The invention of next-generation-sequencing has revolutionized almost all fields of genetics, but few have profited from it as much as the field of ancient DNA research. From its beginnings as an interesting but rather marginal discipline, ancient DNA research is now on its way into the centre of evolutionary biology. In less than a year from its invention next-generation-sequencing had increased the amount of DNA sequence data available from extinct organisms by several orders of magnitude. Ancient DNA research is now not only adding a temporal aspect to evolutionary studies and allowing for the observation of evolution in real time, it also provides important data to help understand the origins of our own species. Here we review progress that has been made in next-generation-sequencing of ancient DNA over the past five years and evaluate sequencing strategies and future directions.

Pseudogenization of the umami taste receptor gene Tas1r1 in the giant panda coincided with its dietary switch to bamboo

Although it belongs to the order Carnivora, the giant panda is a vegetarian with 99% of its diet being bamboo. The draft genome sequence of the giant panda shows that its umami taste receptor gene Tas1r1 is a pseudogene, prompting the proposal that the loss of the umami perception explains why the giant panda is herbivorous. To test this hypothesis, we sequenced all six exons of Tas1r1 in another individual of the giant panda and five other carnivores. We found that the open reading frame (ORF) of Tas1r1 is intact in all these carnivores except the giant panda. The rate ratio (ω) of nonsynonymous to synonymous substitutions in Tas1r1 is significantly higher for the giant panda lineage than for other carnivore lineages. Based on the ω change and the observed number of ORF-disrupting substitutions, we estimated that the functional constraint on the giant panda Tas1r1 was relaxed ∼ 4.2 Ma, with its 95% confidence interval between 1.3 and 10 Ma. Our estimate matches the approximate date of the giant panda's dietary switch inferred from fossil records. It is probable that the giant panda's decreased reliance on meat resulted in the dispensability of the umami taste, leading to Tas1r1 pseudogenization, which in turn reinforced its herbivorous life style because of the diminished attraction of returning to meat eating in the absence of Tas1r1. Nonetheless, additional factors are likely involved because herbivores such as cow and horse still retain an intact Tas1r1.

A multi-locus phylogeny of Nectogalini shrews and influences of the paleoclimate on speciation and evolution

Nectogaline shrews are a major component of the small mammalian fauna of Europe and Asia, and are notable for their diverse ecology, including utilization of aquatic habitats. So far, molecular phylogenetic analyses including nectogaline species have been unable to infer a well-resolved, well-supported phylogeny, thus limiting the power of comparative evolutionary and ecological analyses of the group. Here, we employ Bayesian phylogenetic analyses of eight mitochondrial and three nuclear genes to infer the phylogenetic relationships of nectogaline shrews. We subsequently use this phylogeny to assess the genetic diversity within the genus Episoriculus, and determine whether adaptation to aquatic habitats evolved independently multiple times. Moreover, we both analyze the fossil record and employ Bayesian relaxed clock divergence dating analyses of DNA to assess the impact of historical global climate change on the biogeography of Nectogalini. We infer strong support for the polyphyly of the genus Episoriculus. We also find strong evidence that the ability to heavily utilize aquatic habitats evolved independently in both Neomys and Chimarrogale+Nectogale lineages. Our Bayesian molecular divergence analysis suggests that the early history of Nectogalini is characterized by a rapid radiation at the Miocene/Pliocene boundary, thus potentially explaining the lack of resolution at the base of the tree. Finally, we find evidence that nectogalines once inhabited northern latitudes, but the global cooling and desiccating events at the Miocene/Pliocene and Pliocene/Pleistocene boundaries and Pleistocene glaciation resulted in the migration of most Nectogalini lineages to their present day southern distribution.

Genomic analysis of expressed sequence tags in American black bear Ursus americanus

Background

Species of the bear family (Ursidae) are important organisms for research in molecular evolution, comparative physiology and conservation biology, but relatively little genetic sequence information is available for this group. Here we report the development and analyses of the first large scale Expressed Sequence Tag (EST) resource for the American black bear (Ursus americanus).

Results

Comprehensive analyses of molecular functions, alternative splicing, and tissue-specific expression of 38,757 black bear EST sequences were conducted using the dog genome as a reference. We identified 18 genes, involved in functions such as lipid catabolism, cell cycle, and vesicle-mediated transport, that are showing rapid evolution in the bear lineage Three genes, Phospholamban (PLN), cysteine glycine-rich protein 3 (CSRP3) and Troponin I type 3 (TNNI3), are related to heart contraction, and defects in these genes in humans lead to heart disease. Two genes, biphenyl hydrolase-like (BPHL) and CSRP3, contain positively selected sites in bear. Global analysis of evolution rates of hibernation-related genes in bear showed that they are largely conserved and slowly evolving genes, rather than novel and fast-evolving genes.

Conclusion

We provide a genomic resource for an important mammalian organism and our study sheds new light on the possible functions and evolution of bear genes.

Complete mitochondrial genome of a Pleistocene jawbone unveils the origin of polar bear

The polar bear has become the flagship species in the climate-change discussion. However, little is known about how past climate impacted its evolution and persistence, given an extremely poor fossil record. Although it is undisputed from analyses of mitochondrial (mt) DNA that polar bears constitute a lineage within the genetic diversity of brown bears, timing estimates of their divergence have differed considerably. Using next-generation sequencing technology, we have generated a complete, high-quality mt genome from a stratigraphically validated 130,000- to 110,000-year-old polar bear jawbone. In addition, six mt genomes were generated of extant polar bears from Alaska and brown bears from the Admiralty and Baranof islands of the Alexander Archipelago of southeastern Alaska and Kodiak Island. We show that the phylogenetic position of the ancient polar bear lies almost directly at the branching point between polar bears and brown bears, elucidating a unique morphologically and molecularly documented fossil link between living mammal species. Molecular dating and stable isotope analyses also show that by very early in their evolutionary history, polar bears were already inhabitants of the Artic sea ice and had adapted very rapidly to their current and unique ecology at the top of the Arctic marine food chain. As such, polar bears provide an excellent example of evolutionary opportunism within a widespread mammalian lineage.

Sequencing, annotation and comparative analysis of nine BACs of giant panda (Ailuropoda melanoleuca)

A 10-fold BAC library for giant panda was constructed and nine BACs were selected to generate finish sequences. These BACs could be used as a validation resource for the de novo assembly accuracy of the whole genome shotgun sequencing reads of giant panda newly generated by the Illumina GA sequencing technology. Complete sanger sequencing, assembly, annotation and comparative analysis were carried out on the selected BACs of a joint length 878 kb. Homologue search and de novo prediction methods were used to annotate genes and repeats. Twelve protein coding genes were predicted, seven of which could be functionally annotated. The seven genes have an average gene size of about 41 kb, an average coding size of about 1.2 kb and an average exon number of 6 per gene. Besides, seven tRNA genes were found. About 27 percent of the BAC sequence is composed of repeats. A phylogenetic tree was constructed using neighbor-join algorithm across five species, including giant panda, human, dog, cat and mouse, which reconfirms dog as the most related species to giant panda. Our results provide detailed sequence and structure information for new genes and repeats of giant panda, which will be helpful for further studies on the giant panda.

The sequence and de novo assembly of the giant panda genome

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

Revising the recent evolutionary history of equids using ancient DNA

The rich fossil record of the family Equidae (Mammalia: Perissodactyla) over the past 55 MY has made it an icon for the patterns and processes of macroevolution. Despite this, many aspects of equid phylogenetic relationships and taxonomy remain unresolved. Recent genetic analyses of extinct equids have revealed unexpected evolutionary patterns and a need for major revisions at the generic, subgeneric, and species levels. To investigate this issue we examine 35 ancient equid specimens from four geographic regions (South America, Europe, Southwest Asia, and South Africa), of which 22 delivered 87-688 bp of reproducible aDNA mitochondrial sequence. Phylogenetic analyses support a major revision of the recent evolutionary history of equids and reveal two new species, a South American hippidion and a descendant of a basal lineage potentially related to Middle Pleistocene equids. Sequences from specimens assigned to the giant extinct Cape zebra, Equus capensis, formed a separate clade within the modern plain zebra species, a phenotypicically plastic group that also included the extinct quagga. In addition, we revise the currently recognized extinction times for two hemione-related equid groups. However, it is apparent that the current dataset cannot solve all of the taxonomic and phylogenetic questions relevant to the evolution of Equus. In light of these findings, we propose a rapid DNA barcoding approach to evaluate the taxonomic status of the many Late Pleistocene fossil Equidae species that have been described from purely morphological analyses.

Dogs, cats, and kin: a molecular species-level phylogeny of Carnivora

Phylogenies underpin comparative biology as high-utility tools to test evolutionary and biogeographic hypotheses, inform on conservation strategies, and reveal the age and evolutionary histories of traits and lineages. As tools, most powerful are those phylogenies that contain all, or nearly all, of the taxa of a given group. Despite their obvious utility, such phylogenies, other than summary 'supertrees', are currently lacking for most mammalian orders, including the order Carnivora. Carnivora consists of about 270 extant species including most of the world's large terrestrial predators (e.g., the big cats, wolves, bears), as well as many of man's favorite wild (panda, cheetah, tiger) and domesticated animals (dog, cat). Distributed globally, carnivores are highly diverse ecologically, having occupied all major habitat types on the planet and being diverse in traits such as sociality, communication, body/brain size, and foraging ecology. Thus, numerous studies continue to address comparative questions within the order, highlighting the need for a detailed species-level phylogeny. Here we present a phylogeny of Carnivora that increases taxon sampling density from 28% in the most detailed primary-data study to date, to 82% containing 243 taxa (222 extant species, 17 subspecies). In addition to extant species, we sampled four extinct species: American cheetah, saber-toothed cat, cave bear and the giant short-faced bear. Bayesian analysis of cytochrome b sequences data-mined from GenBank results in a phylogenetic hypothesis that is largely congruent with prior studies based on fewer taxa but more characters. We find support for the monophyly of Carnivora, its major division into Caniformia and Feliformia, and for all but one family within the order. The only exception is the placement of the kinkajou outside Procyonidae, however, prior studies have already cast doubt on its family placement. In contrast, at the subfamily and genus level, our results indicate numerous problems with current classification. Our results also propose new, controversial hypotheses, such as the possible placement of the red panda (Ailuridae) sister to canids (Canidae). Our results confirm previous findings suggesting that the dog was domesticated from the Eurasian wolf (Canis lupus lupus) and are congruent with the Near East domestication of the cat. In sum, this study presents the most detailed species-level phylogeny of Carnivora to date and a much needed tool for comparative studies of carnivoran species. To demonstrate one such use, we perform a phylogenetic analysis of evolutionary distinctiveness (EDGE), which can be used to help establish conservation priorities. According with those criteria, and under one of the many possible sets of parameters, the highest priority Carnivora species for conservation of evolutionary diversity include: monk seals, giant and red panda, giant otter, otter civet, Owston's palm civet, sea otter, Liberian mongoose, spectacled bear, walrus, binturong, and the fossa.

Ancient mitogenomics

The mitochondrial genome has been the traditional focus of most research into ancient DNA, owing to its high copy number and population-level variability. Despite this long-standing interest in mitochondrial DNA, it was only in 2001 that the first complete ancient mitogenomic sequences were obtained. As a result of various methodological developments, including the introduction of high-throughput sequencing techniques, the total number of ancient mitogenome sequences has increased rapidly over the past few years. In this review, we present a brief history of ancient mitogenomics and describe the technical challenges that face researchers in the field. We catalogue the diverse sequencing methods and source materials used to obtain ancient mitogenomic sequences, summarise the associated genetic and phylogenetic studies that have been conducted, and evaluate the future prospects of the field.

Predormancy omnivory in European cave bears evidenced by a dental microwear analysis of Ursus spelaeus from Goyet, Belgium

Previous morphological and isotopic studies indicate that Late Pleistocene cave bear (Ursus spelaeus) diet ranged from mostly vegetarian to omnivory or even carnivory. However, such analyses do not provide information on seasonal diets, and only provide an average record of diet. A dental microwear analysis of 43 young and adult individuals demonstrate that, during the predormancy period, cave bears from Goyet (Late Pleistocene, Belgium) were not strictly herbivorous, but had a mixed diet composed of hard items (e.g., possibly bone), invertebrates (e.g., insects), meat (ungulates, small vertebrates), and/or plant matter (hard mast, seeds, herbaceous vegetations, and fruits). Therefore, our results indicate that cave bears at Goyet were generalist omnivores during the predormancy period, which is consistent with current data on the dietary ecology of extant bears during this season. These data also raise questions about the ecological role and causes of the extinction of cave bears.

Direct multiplex sequencing (DMPS)--a novel method for targeted high-throughput sequencing of ancient and highly degraded DNA

Although the emergence of high-throughput sequencing technologies has enabled whole-genome sequencing from extinct organisms, little progress has been made in accelerating targeted sequencing from highly degraded DNA. Here, we present a novel and highly sensitive method for targeted sequencing of ancient and degraded DNA, which couples multiplex PCR directly with sample barcoding and high-throughput sequencing. Using this approach, we obtained a 96% complete mitochondrial genome data set from 31 cave bear (Ursus spelaeus) samples using only two 454 Life Sciences (Roche) GS FLX runs. In contrast to previous studies relying only on short sequence fragments, the overlapping portion of our data comprises almost 10 kb of replicated mitochondrial genome sequence, allowing for the unambiguous differentiation of three major cave bear clades. Our method opens up the opportunity to simultaneously generate many kilobases of overlapping sequence data from large sets of difficult samples, such as museum specimens, medical collections, or forensic samples. Embedded in our approach, we present a new protocol for the construction of barcoded sequencing libraries, which is compatible with all current high-throughput technologies and can be performed entirely in plate setup.

Deciphering the complete mitochondrial genome and phylogeny of the extinct cave bear in the Paleolithic painted cave of Chauvet

Retrieving a large amount of genetic information from extinct species was demonstrated feasible, but complete mitochondrial genome sequences have only been deciphered for the moa, a bird that became extinct a few hundred years ago, and for Pleistocene species, such as the woolly mammoth and the mastodon, both of which could be studied from animals embedded in permafrost. To enlarge the diversity of mitochondrial genomes available for Pleistocene species, we turned to the cave bear (Ursus spelaeus), whose only remains consist of skeletal elements. We collected bone samples from the Paleolithic painted cave of Chauvet-Pont d'Arc (France), which displays the earliest known human drawings, and contains thousands of bear remains. We selected a cave bear sternebra, radiocarbon dated to 32,000 years before present, from which we generated overlapping DNA fragments assembling into a 16,810-base pair mitochondrial genome. Together with the first mitochondrial genome for the brown bear western lineage, this study provides a statistically secured molecular phylogeny assessing the cave bear as a sister taxon to the brown bear and polar bear clade, with a divergence inferred to 1.6 million years ago. With the first mitochondrial genome for a Pleistocene carnivore to be delivered, our study establishes the Chauvet-Pont d'Arc Cave as a new reservoir for Paleogenetic studies. These molecular data enable establishing the chronology of bear speciation, and provide a helpful resource to rescue for genetic analysis archeological samples initially diagnosed as devoid of amplifiable DNA.