Human Papilloma Virus-Independent/p53abnormal Keratinizing Squamous Cell Carcinoma of the Uterine Cervix Associated With Uterine Prolapse.

Knowledge about the morphologic and molecular characteristics of cervical squamous cell carcinomas (CSCCs) associated with uterine prolapse is very limited. Detailed histopathological and immunohistochemical (p16, p53, and cytokeratin 17), as well as molecular evaluation for human papillomavirus (HPV)-DNA and p53-mutational analyses in 4 consecutive CSCCs associated with uterine prolapse with definition of a hitherto not well-described HPV-independent/p53abnormal precursor lesion (HPV-independent cervical intraepithelial neoplasia [CIN; differentiated CIN]) and molecular tumorigenetic pathway. Cases diagnosed within 7 years with a mean age of 75 (range: 69-83) years and a mean tumor size of 7.3 cm (range: 5.2-9.4 cm). All patients presented with locally advanced disease, and 1 woman died of the disease within 4, and another within 14 months of follow-up. All CSCCs and their adjacent precursor lesions were negative for p16, with aberrant p53-expression and diffuse and strong staining for cytokeratin 17. Both the CSCCs and their precursors were negative for HPV-DNA but harbored a TP53 mutation. The precursor lesions were characterized by epithelial thickening with superficial keratinization, and the presence of basal and parabasal keratinocytes with mitotic figures beyond the basal layer, thus showing features similar to those seen in differentiated types of vulvar intraepithelial lesions (vulvar intraepithelial neoplasia [VIN] syn. HPV-independent/p53abn VIN), suggesting the terminology of differentiated CIN or HPV-independent/p53abn CIN. An HPV-independent pathogenetic pathway with a p53-alteration was identified for these cases. CSCC associated with uterine prolapse represents HPV-independent tumors harboring a TP53 mutation. For the first time, a precursor lesion of HPV-independent CSCC of the uterine cervix is described with a differentiated VIN-like morphology, and a separate tumorigenic pathway defined.

The genetic landscape of intellectual disability and epilepsy in adults and the elderly: a systematic genetic work-up of 150 individuals.

PURPOSE: Genetic diagnostics of neurodevelopmental disorders with epilepsy (NDDE) are predominantly applied in children, thus limited information is available regarding adults or elderly. METHODS: We investigated 150 adult/elderly individuals with NDDE by conventional karyotyping, FMR1 testing, chromosomal microarray, panel sequencing, and for unresolved cases, also by exome sequencing (nsingle = 71, ntrios = 24). RESULTS: We identified (likely) pathogenic variants in 71 cases (47.3%) comprising fragile X syndrome (n = 1), disease-causing copy number (n = 23), and single-nucleotide variants (n = 49). Seven individuals displayed multiple independent genetic diagnoses. The diagnostic yield correlated with the severity of intellectual disability. Individuals with anecdotal evidence of exogenic early-life events (e.g., nuchal cord, complications at delivery) with alleged/unproven association to the disorder had a particularly high yield of 58.3%. Screening for disease-specific comorbidities was indicated in 45.1% and direct treatment consequences arose in 11.8% of diagnosed individuals. CONCLUSION: Panel/exome sequencing displayed the highest yield and should be considered as first-tier diagnostics in NDDE. This high yield and the numerous indications for additional screening or treatment modifications arising from genetic diagnoses indicate a current medical undersupply of genetically undiagnosed adult/elderly individuals with NDDE. Moreover, knowledge of the course of elderly individuals will ultimately help in counseling newly diagnosed individuals with NDDE.

Ancient horse genomes reveal the timing and extent of dispersals across the Bering Land Bridge.

The Bering Land Bridge (BLB) last connected Eurasia and North America during the Late Pleistocene. Although the BLB would have enabled transfers of terrestrial biota in both directions, it also acted as an ecological filter whose permeability varied considerably over time. Here we explore the possible impacts of this ecological corridor on genetic diversity within, and connectivity among, populations of a once wide-ranging group, the caballine horses (Equus spp.). Using a panel of 187 mitochondrial and eight nuclear genomes recovered from present-day and extinct caballine horses sampled across the Holarctic, we found that Eurasian horse populations initially diverged from those in North America, their ancestral continent, around 1.0-0.8 million years ago. Subsequent to this split our mitochondrial DNA analysis identified two bidirectional long-range dispersals across the BLB ~875-625 and ~200-50 thousand years ago, during the Middle and Late Pleistocene. Whole genome analysis indicated low levels of gene flow between North American and Eurasian horse populations, which probably occurred as a result of these inferred dispersals. Nonetheless, mitochondrial and nuclear diversity of caballine horse populations retained strong phylogeographical structuring. Our results suggest that barriers to gene flow, currently unidentified but possibly related to habitat distribution across Beringia or ongoing evolutionary divergence, played an important role in shaping the early genetic history of caballine horses, including the ancestors of living horses within Equus ferus.

Fossil and genomic evidence constrains the timing of bison arrival in North America.

The arrival of bison in North America marks one of the most successful large-mammal dispersals from Asia within the last million years, yet the timing and nature of this event remain poorly determined. Here, we used a combined paleontological and paleogenomic approach to provide a robust timeline for the entry and subsequent evolution of bison within North America. We characterized two fossil-rich localities in Canada's Yukon and identified the oldest well-constrained bison fossil in North America, a 130,000-y-old steppe bison, Bison cf. priscus We extracted and sequenced mitochondrial genomes from both this bison and from the remains of a recently discovered, ∼120,000-y-old giant long-horned bison, Bison latifrons, from Snowmass, Colorado. We analyzed these and 44 other bison mitogenomes with ages that span the Late Pleistocene, and identified two waves of bison dispersal into North America from Asia, the earliest of which occurred ∼195-135 thousand y ago and preceded the morphological diversification of North American bison, and the second of which occurred during the Late Pleistocene, ∼45-21 thousand y ago. This chronological arc establishes that bison first entered North America during the sea level lowstand accompanying marine isotope stage 6, rejecting earlier records of bison in North America. After their invasion, bison rapidly colonized North America during the last interglaciation, spreading from Alaska through continental North America; they have been continuously resident since then.

Parental mosaicism in epilepsies due to alleged de novo variants.

Severe early onset epilepsies are often caused by de novo pathogenic variants. Few studies have reported the frequency of somatic mosaicism in parents of children with severe epileptic encephalopathies. Here we aim to investigate the frequency of mosaicism in the parents of children with epilepsy caused by alleged de novo variants. We tested parental genomic DNA derived from different tissues for 75 cases using targeted next-generation sequencing. Five parents (6.6%) showed mosaicism at minor allele frequencies of 0.8%-29% for the pathogenic variant detected in their offspring. Parental mosaicism was observed in the following genes: SCN1A, SCN2A, SCN8A, and STXBP1. One of the identified parents had epilepsy himself. Our results show that de novo events can occur already in parental tissue and in some cases can be detected in peripheral blood. Consequently, parents affected by low-grade mosaicism are faced with an increased recurrence risk for transmitting the pathogenic variant, compared to the overall recurrence risk for a second affected child estimated at approximately 1%. However, testing for parental somatic mosaicism will help identifying those parents who truly are at higher risk and will significantly improve genetic counseling in the respective families.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse.

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560-780 thousand years before present (kyr BP). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr BP), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski's horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0-4.5 million years before present (Myr BP), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski's and domestic horse populations diverged 38-72 kyr BP, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski's horse investigated. This supports the contention that Przewalski's horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski's and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski's horse. Such regions could correspond to loci selected early during domestication.

Bison phylogeography constrains dispersal and viability of the Ice Free Corridor in western Canada.

The Ice Free Corridor has been invoked as a route for Pleistocene human and animal dispersals between eastern Beringia and more southerly areas of North America. Despite the significance of the corridor, there are limited data for when and how this corridor was used. Hypothetical uses of the corridor include: the first expansion of humans from Beringia into the Americas, northward postglacial expansions of fluted point technologies into Beringia, and continued use of the corridor as a contact route between the north and south. Here, we use radiocarbon dates and ancient mitochondrial DNA from late Pleistocene bison fossils to determine the chronology for when the corridor was open and viable for biotic dispersals. The corridor was closed after ∼23,000 until 13,400 calendar years ago (cal y BP), after which we find the first evidence, to our knowledge, that bison used this route to disperse from the south, and by 13,000 y from the north. Our chronology supports a habitable and traversable corridor by at least 13,000 cal y BP, just before the first appearance of Clovis technology in interior North America, and indicates that the corridor would not have been available for significantly earlier southward human dispersal. Following the opening of the corridor, multiple dispersals of human groups between Beringia and interior North America may have continued throughout the latest Pleistocene and early Holocene. Our results highlight the utility of phylogeographic analyses to test hypotheses about paleoecological history and the viability of dispersal routes over time.

Spatiotemporal Dynamics of Genetic Variation in the Iberian Lynx along Its Path to Extinction Reconstructed with Ancient DNA.

There is the tendency to assume that endangered species have been both genetically and demographically healthier in the past, so that any genetic erosion observed today was caused by their recent decline. The Iberian lynx (Lynx pardinus) suffered a dramatic and continuous decline during the 20th century, and now shows extremely low genome- and species-wide genetic diversity among other signs of genomic erosion. We analyze ancient (N = 10), historical (N = 245), and contemporary (N = 172) samples with microsatellite and mitogenome data to reconstruct the species' demography and investigate patterns of genetic variation across space and time. Iberian lynx populations transitioned from low but significantly higher genetic diversity than today and shallow geographical differentiation millennia ago, through a structured metapopulation with varying levels of diversity during the last centuries, to two extremely genetically depauperate and differentiated remnant populations by 2002. The historical subpopulations show varying extents of genetic drift in relation to their recent size and time in isolation, but these do not predict whether the populations persisted or went finally extinct. In conclusion, current genetic patterns were mainly shaped by genetic drift, supporting the current admixture of the two genetic pools and calling for a comprehensive genetic management of the ongoing conservation program. This study illustrates how a retrospective analysis of demographic and genetic patterns of endangered species can shed light onto their evolutionary history and this, in turn, can inform conservation actions.

Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch.

Viruses preserved in ancient materials provide snapshots of past viral diversity and a means to trace viral evolution through time. Here, we use a metagenomics approach to identify filterable and nuclease-resistant nucleic acids preserved in 700-y-old caribou feces frozen in a permanent ice patch. We were able to recover and characterize two viruses in replicated experiments performed in two different laboratories: a small circular DNA viral genome (ancient caribou feces associated virus, or aCFV) and a partial RNA viral genome (Ancient Northwest Territories cripavirus, or aNCV). Phylogenetic analysis identifies aCFV as distantly related to the plant-infecting geminiviruses and the fungi-infecting Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1 and aNCV as within the insect-infecting Cripavirus genus. We hypothesize that these viruses originate from plant material ingested by caribou or from flying insects and that their preservation can be attributed to protection within viral capsids maintained at cold temperatures. To investigate the tropism of aCFV, we used the geminiviral reverse genetic system and introduced a multimeric clone into the laboratory model plant Nicotiana benthamiana. Evidence for infectivity came from the detection of viral DNA in newly emerged leaves and the precise excision of the viral genome from the multimeric clones in inoculated leaves. Our findings indicate that viral genomes may in some circumstances be protected from degradation for centuries.

Targeted next generation sequencing reveals unique mutation profile of primary melanocytic tumors of the central nervous system.

Melanocytic tumors originating in the central nervous system (MT-CNS) are rare tumors that generally have a favorable prognosis, however malignant tumors do occur. Pathogenetically MT-CNS are not well characterized. Similar to uveal melanoma and blue nevi, they frequently harbor activating GNAQ or GNA11 mutations. Rare NRAS mutations have also been reported. Other mutations have not yet been described. We analyzed 19 MT-CNS, 7 uveal melanomas and 19 cutaneous melanomas using a targeted next generation sequencing approach analyzing 29 genes known to be frequently mutated in other melanocytic tumors (in particular uveal and cutaneous melanomas). In concordance with previous studies, cutaneous melanoma samples showed frequent NRAS or BRAF mutations, as well as mutations in other genes (e.g. NF1, RAC1, PIK3CA, ARID1A). Metastasized uveal melanomas exhibited mutations in GNAQ, GNA11 and BAP1. In contrast, MT-CNS almost exclusively demonstrated mutations in GNAQ (71 %) or GNA11 (12 %). Interestingly both GNA11 mutations identified were detected in MT-CNS diagnosed as intermediate grade melanocytomas which also recurred. One of these recurrent cases also harbored an inactivating BAP1 mutation and was found to have lost one copy of chromosome 3. Our findings show that while MT-CNS do have GNAQ or GNA11 mutations, they rarely harbor other recurrent mutations found in uveal or cutaneous melanomas. Considering chromosome 3 and BAP1 loss are robust markers of poor prognosis in uveal melanoma, it will prove interesting to determine whether these genomic alterations are also of prognostic significance in MT-CNS.

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.

Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears.

Polar bears are an arctic, marine adapted species that is closely related to brown bears. Genome analyses have shown that polar bears are distinct and genetically homogeneous in comparison to brown bears. However, these analyses have also revealed a remarkable episode of polar bear gene flow into the population of brown bears that colonized the Admiralty, Baranof and Chichagof islands (ABC islands) of Alaska. Here, we present an analysis of data from a large panel of polar bear and brown bear genomes that includes brown bears from the ABC islands, the Alaskan mainland and Europe. Our results provide clear evidence that gene flow between the two species had a geographically wide impact, with polar bear DNA found within the genomes of brown bears living both on the ABC islands and in the Alaskan mainland. Intriguingly, while brown bear genomes contain up to 8.8% polar bear ancestry, polar bear genomes appear to be devoid of brown bear ancestry, suggesting the presence of a barrier to gene flow in that direction.

Ancient mitochondrial DNA and the genetic history of Eurasian beaver (Castor fiber) in Europe.

After centuries of human hunting, the Eurasian beaver Castor fiber had disappeared from most of its original range by the end of the 19th century. The surviving relict populations are characterized by both low genetic diversity and strong phylogeographical structure. However, it remains unclear whether these attributes are the result of a human-induced, late Holocene bottleneck or already existed prior to this reduction in range. To investigate genetic diversity in Eurasian beaver populations during the Holocene, we obtained mitochondrial control region DNA sequences from 48 ancient beaver samples and added 152 modern sequences from GenBank. Phylogeographical analyses of the data indicate a differentiation of European beaver populations into three mitochondrial clades. The two main clades occur in western and eastern Europe, respectively, with an early Holocene contact zone in eastern Europe near a present-day contact zone. A divergent and previously unknown clade of beavers from the Danube Basin survived until at least 6000 years ago, but went extinct during the transition to modern times. Finally, we identify a recent decline in effective population size of Eurasian beavers, with a stronger bottleneck signal in the western than in the eastern clade. Our results suggest that the low genetic diversity and the strong phylogeographical structure in recent beavers are artefacts of human hunting-associated population reductions. While beaver populations have been growing rapidly since the late 19th century, genetic diversity within modern beaver populations remains considerably reduced compared to what was present prior to the period of human hunting and habitat reduction.

Interactive webtool for analyzing drug sensitivity and resistance associated with genetic signatures of cancer cell lines.

PURPOSE: A wide therapeutic repertoire has become available to oncologists including radio- and chemotherapy, small molecules and monoclonal antibodies. However, drug efficacy can be limited by genetic heterogeneity. Here, we designed a webtool that facilitates the data analysis of the in vitro drug sensitivity data on 265 approved compounds from the GDSC database in association with a plethora of genetic changes documented for 1001 cell lines in the CCLE data. METHODS: The webtool computes odds ratios of drug resistance for a queried set of genetic alterations. It provides results on the efficacy of single compounds or groups of compounds assigned to cellular signaling pathways. Webtool availability: https://tools.hornlab.org/GDSC/ . RESULTS: We first replicated established associations of genetic driver mutations in BRAF, RAS genes and EGFR with drug response. We then tested the 'BRCAness' hypothesis and did not find increased sensitivity to the assayed PARP inhibitors. Analyzing specific PIK3CA mutations related to cancer and mendelian overgrowth, we found support for the described sensitivity of H1047 mutants to GSK690693 targeting the AKT pathway. Testing a co-mutated gene pair, GATA3 activation abolished PTEN-related sensitivity to PI3K/mTOR inhibition. Finally, the pharmacogenomic modifier ABCB1 was associated with olaparib resistance. CONCLUSIONS: This tool could identify potential drug candidates in the presence of custom sets of genetic changes and moreover, improve the understanding of signaling pathways. The underlying computer code can be adapted to larger drug response datasets to help structure and accommodate the increasingly large biomedical knowledge base.

TERT expression is associated with metastasis from thin primaries, exhausted CD4+ T cells in melanoma and with DNA repair across cancer entities.

Telomerase reverse transcriptase (TERT) promoter mutations occur frequently in cancer, have been associated with increased TERT expression and cell proliferation, and could potentially influence therapeutic regimens for melanoma. As the role of TERT expression in malignant melanoma and the non-canonical functions of TERT remain understudied, we aimed to extend the current knowledge on the impact of TERT promoter mutations and expression alterations in tumor progression by analyzing several highly annotated melanoma cohorts. Using multivariate models, we found no consistent association for TERT promoter mutations or TERT expression with the survival rate in melanoma cohorts under immune checkpoint inhibition. However, the presence of CD4+ T cells increased with TERT expression and correlated with the expression of exhaustion markers. While the frequency of promoter mutations did not change with Breslow thickness, TERT expression was increased in metastases arising from thinner primaries. As single-cell RNA-sequencing (RNA-seq) showed that TERT expression was associated with genes involved in cell migration and dynamics of the extracellular matrix, this suggests a role of TERT during invasion and metastasis. Co-regulated genes found in several bulk tumors and single-cell RNA-seq cohorts also indicated non-canonical functions of TERT related to mitochondrial DNA stability and nuclear DNA repair. This pattern was also evident in glioblastoma and across other entities. Hence, our study adds to the role of TERT expression in cancer metastasis and potentially also immune resistance.

Database for queries of melanoma cohorts with molecular and clinical data.

Large genome-scale studies are deposited in various public sequence repositories. However, their access and analysis can be non-trivial to infrequent users. Here, we present a new database connecting whole transcriptomes with clinical data for straight-forward access and analysis of patient-specific samples. Users can perform association tests of survival and gene expression across different cohorts, identify cell-type expressions, or correlate the presence of immune cells. In summary, we present a new data hub for bench scientists to perform replication and discovery studies.

Comparative bioinformatic analysis of KRAS, STK11 and KEAP1 (co-)mutations in non-small cell lung cancer with a special focus on KRAS G12C.

OBJECTIVES: Mutations in STK11 (STK11MUT) and KEAP1 (KEAP1MUT) occur frequently in non-small cell lung cancer (NSCLC) and are often co-mutated with KRAS. Several studies linked the co-occurrence of KRASMUT + STK11MUT, as well as KRASMUT + KEAP1MUT to reduced response to immune checkpoint inhibitors (ICI) and even a negative impact on survival. Data focusing STK11 + KEAP1 co-mutations or the triple mutation (KRAS + STK11 + KEAP1) are scarce. The recent availability of KRAS-G12C inhibitors increases the clinical relevance of those co-mutations in KRAS-mutated NSCLC. MATERIALS AND METHODS: We present a comprehensive bioinformatic analysis encompassing six datasets retrieved from cBioPortal. RESULTS: Independent of the treatment, triple mutations and STK11MUT + KEAP1MUT were significantly associated with a reduced overall survival (OS). Across treatments, OS of patients with a KRAS G12C triple mutation was significantly reduced compared to patients with KRAS G12C-only. Under ICI-therapy, there was no significant difference in OS between patients harboring the KRAS G12C-only and patients with the KRAS G12C triple mutation, but a significant difference between patients harboring KRAS non-G12C and KRAS non-G12C triple mutations. Triple mutated primary tumors showed a significantly increased frequency of distant metastases to bone and adrenal glands compared to KRAS-only mutated tumors. Additionally, our drug response analysis in cancer cell lines harboring the triple mutations revealed the WNT pathway inhibitor XAV-939 as a potential future drug candidate for this mutational situation. CONCLUSION: The triple mutation status may serve as a negative prognostic and predictive factor across treatments compared to KRASMUT-only. KRAS G12C generally seems to be a negative predictive marker for ICI-therapy.

Improving the performance of true single molecule sequencing for ancient DNA.

BACKGROUND: Second-generation sequencing technologies have revolutionized our ability to recover genetic information from the past, allowing the characterization of the first complete genomes from past individuals and extinct species. Recently, third generation Helicos sequencing platforms, which perform true Single-Molecule DNA Sequencing (tSMS), have shown great potential for sequencing DNA molecules from Pleistocene fossils. Here, we aim at improving even further the performance of tSMS for ancient DNA by testing two novel tSMS template preparation methods for Pleistocene bone fossils, namely oligonucleotide spiking and treatment with DNA phosphatase. RESULTS: We found that a significantly larger fraction of the horse genome could be covered following oligonucleotide spiking however not reproducibly and at the cost of extra post-sequencing filtering procedures and skewed %GC content. In contrast, we showed that treating ancient DNA extracts with DNA phosphatase improved the amount of endogenous sequence information recovered per sequencing channel by up to 3.3-fold, while still providing molecular signatures of endogenous ancient DNA damage, including cytosine deamination and fragmentation by depurination. Additionally, we confirmed the existence of molecular preservation niches in large bone crystals from which DNA could be preferentially extracted. CONCLUSIONS: We propose DNA phosphatase treatment as a mechanism to increase sequence coverage of ancient genomes when using Helicos tSMS as a sequencing platform. Together with mild denaturation temperatures that favor access to endogenous ancient templates over modern DNA contaminants, this simple preparation procedure can improve overall Helicos tSMS performance when damaged DNA templates are targeted.

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.

Withering away--25,000 years of genetic decline preceded cave bear extinction.

The causes of the late Pleistocene megafaunal extinctions are still enigmatic. Although the fossil record can provide approximations for when a species went extinct, the timing of its disappearance alone cannot resolve the causes and mode of the decline preceding its extinction. However, ancient DNA analyses can reveal population size changes over time and narrow down potential causes of extinction. Here, we present an ancient DNA study comparing late Pleistocene population dynamics of two closely related species, cave and brown bears. We found that the decline of cave bears started approximately 25,000 years before their extinction, whereas brown bear population size remained constant. We conclude that neither the effects of climate change nor human hunting alone can be responsible for the decline of the cave bear and suggest that a complex of factors including human competition for cave sites lead to the cave bear's extinction.