African and Asian leopards are highly differentiated at the genomic level.

Leopards are the only big cats still widely distributed across the continents of Africa and Asia. They occur in a wide range of habitats and are often found in close proximity to humans. But despite their ubiquity, leopard phylogeography and population history have not yet been studied with genomic tools. Here, we present population-genomic data from 26 modern and historical samples encompassing the vast geographical distribution of this species. We find that Asian leopards are broadly monophyletic with respect to African leopards across almost their entire nuclear genomes. This profound genetic pattern persists despite the animals' high potential mobility, and despite evidence of transfer of African alleles into Middle Eastern and Central Asian leopard populations within the last 100,000 years. Our results further suggest that Asian leopards originated from a single out-of-Africa dispersal event 500-600 thousand years ago and are characterized by higher population structuring, stronger isolation by distance, and lower heterozygosity than African leopards. Taxonomic categories do not take into account the variability in depth of divergence among subspecies. The deep divergence between the African subspecies and Asian populations contrasts with the much shallower divergence among putative Asian subspecies. Reconciling genomic variation and taxonomy is likely to be a growing challenge in the genomics era.

Ancestral mitogenome capture of the Southeast Asian banded linsang.

Utilising a reconstructed ancestral mitochondrial genome of a clade to design hybridisation capture baits can provide the opportunity for recovering mitochondrial sequences from all its descendent and even sister lineages. This approach is useful for taxa with no extant close relatives, as is often the case for rare or extinct species, and is a viable approach for the analysis of historical museum specimens. Asiatic linsangs (genus Prionodon) exemplify this situation, being rare Southeast Asian carnivores for which little molecular data is available. Using ancestral capture we recover partial mitochondrial genome sequences for seven banded linsangs (P. linsang) from historical specimens, representing the first intraspecific genetic dataset for this species. We additionally assemble a high quality mitogenome for the banded linsang using shotgun sequencing for time-calibrated phylogenetic analysis. This reveals a deep divergence between the two Asiatic linsang species (P. linsang, P. pardicolor), with an estimated divergence of ~12 million years (Ma). Although our sample size precludes any robust interpretation of the population structure of the banded linsang, we recover two distinct matrilines with an estimated tMRCA of ~1 Ma. Our results can be used as a basis for further investigation of the Asiatic linsangs, and further demonstrate the utility of ancestral capture for studying divergent taxa without close relatives.

Molecular insights into an ancient form of Paget's disease of bone.

Paget's disease of bone (PDB) is a chronic skeletal disorder that can affect one or several bones in individuals older than 55 y of age. PDB-like changes have been reported in archaeological remains as old as Roman, although accurate diagnosis and natural history of the disease is lacking. Six skeletons from a collection of 130 excavated at Norton Priory in the North West of England, which dates to medieval times, show atypical and extensive pathological changes resembling contemporary PDB affecting as many as 75% of individual skeletons. Disease prevalence in the remaining collection is high, at least 16% of adults, with age at death estimations as low as 35 y. Despite these atypical features, paleoproteomic analysis identified sequestosome 1 (SQSTM1) or p62, a protein central to the pathological milieu of PDB, as one of the few noncollagenous human sequences preserved in skeletal samples. Targeted proteomic analysis detected >60% of the ancient p62 primary sequence, with Western blotting indicating p62 abnormalities, including in dentition. Direct sequencing of ancient DNA excluded contemporary PDB-associated SQSTM1 mutations. Our observations indicate that the ancient p62 protein is likely modified within its C-terminal ubiquitin-associated domain. Ancient miRNAs were remarkably preserved in an osteosarcoma from a skeleton with extensive disease, with miR-16 expression consistent with that reported in contemporary PDB-associated bone tumors. Our work displays the use of proteomics to inform diagnosis of ancient diseases such as atypical PDB, which has unusual features presumably potentiated by yet-unidentified environmental or genetic factors.

Double-Stranded Library Preparation for Ancient and Other Degraded Samples.

High-throughput sequencing (HTS) allows fast and cost-efficient sequencing of ancient DNA (aDNA) without prior information about what sequences should be targeted. One necessary step for HTS is the preparation of a sequencing library. Commercial kits are available for this purpose, but many of these are not suitable for aDNA or other types of damaged DNA. Here, we outline a protocol for HTS library preparation that is optimized for ancient DNA. We report the library conversion rate for a range of input template and adapter concentrations. Our results show that the protocol performs at a high efficiency.

Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations.

BACKGROUND: Resolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts? RESULTS: In this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these. CONCLUSIONS: The phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.

Optimized DNA sampling of ancient bones using Computed Tomography scans.

The prevalence of contaminant microbial DNA in ancient bone samples represents the principal limiting factor for palaeogenomic studies, as it may comprise more than 99% of DNA molecules obtained. Efforts to exclude or reduce this contaminant fraction have been numerous but also variable in their success. Here, we present a simple but highly effective method to increase the relative proportion of endogenous molecules obtained from ancient bones. Using computed tomography (CT) scanning, we identify the densest region of a bone as optimal for sampling. This approach accurately identifies the densest internal regions of petrous bones, which are known to be a source of high-purity ancient DNA. For ancient long bones, CT scans reveal a high-density outermost layer, which has been routinely removed and discarded prior to DNA extraction. For almost all long bones investigated, we find that targeted sampling of this outermost layer provides an increase in endogenous DNA content over that obtained from softer, trabecular bone. This targeted sampling can produce as much as 50-fold increase in the proportion of endogenous DNA, providing a directly proportional reduction in sequencing costs for shotgun sequencing experiments. The observed increases in endogenous DNA proportion are not associated with any reduction in absolute endogenous molecule recovery. Although sampling the outermost layer can result in higher levels of human contamination, some bones were found to have more contamination associated with the internal bone structures. Our method is highly consistent, reproducible and applicable across a wide range of bone types, ages and species. We predict that this discovery will greatly extend the potential to study ancient populations and species in the genomics era.

eQTL mapping of transposon silencing reveals a position-dependent stable escape from epigenetic silencing and transposition of AtMu1 in the Arabidopsis lineage.

Transposons are massively abundant in all eukaryotic genomes and are suppressed by epigenetic silencing. Transposon activity contributes to the evolution of species; however, it is unclear how much transposition-induced variation exists at a smaller scale and how transposons are targeted for silencing. Here, we exploited differential silencing of the AtMu1c transposon in the Arabidopsis thaliana accessions Columbia (Col) and Landsberg erecta (Ler). The difference persisted in hybrids and recombinant inbred lines and was mapped to a single expression quantitative trait locus within a 20-kb interval. In Ler only, this interval contained a previously unidentified copy of AtMu1c, which was inserted at the 3' end of a protein-coding gene and showed features of expressed genes. By contrast, AtMu1c(Col) was intergenic and associated with heterochromatic features. Furthermore, we identified widespread natural AtMu1c transposition from the analysis of over 200 accessions, which was not evident from alignments to the reference genome. AtMu1c expression was highest for insertions within 3' untranslated regions, suggesting that this location provides protection from silencing. Taken together, our results provide a species-wide view of the activity of one transposable element at unprecedented resolution, showing that AtMu1c transposed in the Arabidopsis lineage and that transposons can escape epigenetic silencing by inserting into specific genomic locations, such as the 3' end of genes.