Thirty years of ancient DNA and the faunal biogeography of Aotearoa New Zealand: lessons and future directions.

Thirty years ago, DNA sequences were obtained from an extinct Aotearoa New Zealand animal for the first time. Since then, ancient DNA research has provided many - often unexpected - insights into the origins of New Zealand's terrestrial and marine vertebrate fauna. Because recent human activities in New Zealand have caused the decline or extinction of many endemic plant, bird, reptile, and marine mammal species, ancient DNA has been instrumental in reconstructing their identities and origins. However, most ancient DNA studies focusing on New Zealand species have been restricted to vertebrates, with small sample sizes, and/or relatively few genetic markers. This has limited their power to infer fine-scale biogeographic patterns, including (pre)historic distributions and range-shifts driven by past climate and environmental change. Recently, 'next-generation' methodological and technological advances have broadened the range of hypotheses that can feasibly be tested with ancient DNA. These advances represent an exciting opportunity for further exploring New Zealand biogeography using ancient DNA, but their promise has not yet been fully realised. In this review, we summarise the last 30 years of ancient DNA research into New Zealand faunal biogeography and highlight key objectives, challenges, and possibilities for the next 30 years and beyond.

Ancient genomes reveal over two thousand years of dingo population structure.

Dingoes are culturally and ecologically important free-living canids whose ancestors arrived in Australia over 3,000 B.P., likely transported by seafaring people. However, the early history of dingoes in Australia-including the number of founding populations and their routes of introduction-remains uncertain. This uncertainty arises partly from the complex and poorly understood relationship between modern dingoes and New Guinea singing dogs, and suspicions that post-Colonial hybridization has introduced recent domestic dog ancestry into the genomes of many wild dingo populations. In this study, we analyzed genome-wide data from nine ancient dingo specimens ranging in age from 400 to 2,746 y old, predating the introduction of domestic dogs to Australia by European colonists. We uncovered evidence that the continent-wide population structure observed in modern dingo populations had already emerged several thousand years ago. We also detected excess allele sharing between New Guinea singing dogs and ancient dingoes from coastal New South Wales (NSW) compared to ancient dingoes from southern Australia, irrespective of any post-Colonial hybrid ancestry in the genomes of modern individuals. Our results are consistent with several demographic scenarios, including a scenario where the ancestry of dingoes from the east coast of Australia results from at least two waves of migration from source populations with varying affinities to New Guinea singing dogs. We also contribute to the growing body of evidence that modern dingoes derive little genomic ancestry from post-Colonial hybridization with other domestic dog lineages, instead descending primarily from ancient canids introduced to Sahul thousands of years ago.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago.

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

Genomic ancestry and social dynamics of the last hunter-gatherers of Atlantic France.

Since the early Holocene, western and central Europe was inhabited by a genetically distinct group of Western Hunter-Gatherers (WHGs). This group was eventually replaced and assimilated by the incoming Neolithic farmers. The western Atlantic façade was home to some of the last Mesolithic sites of mainland Europe, represented by the iconic open-air sites at Hoedic and Téviec in southern Brittany, France. These sites are known for the unusually well-preserved and rich burials. Genomic studies of Mesolithic European hunter-gatherers have been limited to single or a few individuals per site and our understanding of the social dynamics of the last Mesolithic hunter-gatherers of Europe and their interactions with incoming farmers is limited. We sequenced and analyzed the complete genomes of 10 individuals from the Late Mesolithic sites of Hoedic, Téviec, and Champigny, in France, four of which sequenced to between 23- and 8-times genome coverage. The analysis of genomic, chronological and dietary data revealed that the Late Mesolithic populations in Brittany maintained distinct social units within a network of exchanging mates. This resulted in low intra-group biological relatedness that prevented consanguineous mating, despite the small population size of the Late Mesolithic groups. We found no genetic ancestry from Neolithic farmers in the analyzed hunter-gatherers, even though some of them may have coexisted with the first farming groups in neighboring regions. Hence, contrary to previous conclusions based on stable isotope data from the same sites, the Late Mesolithic forager community was limited in mate-exchange to neighboring hunter-gatherer groups, to the exclusion of Neolithic farmers.

Beringia and the peopling of the Western Hemisphere.

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

Ancient mitochondrial genomes recovered from small vertebrate bones through minimally destructive DNA extraction: Phylogeography of the New Zealand gecko genus Hoplodactylus.

Methodological and technological improvements are continually revolutionizing the field of ancient DNA. Most ancient DNA extraction methods require the partial (or complete) destruction of finite museum specimens, which disproportionately impacts small or fragmentary subfossil remains, and future analyses. We present a minimally destructive ancient DNA extraction method optimized for small vertebrate remains. We applied this method to detect lost mainland genetic diversity in the large New Zealand diplodactylid gecko genus Hoplodactylus, which is presently restricted to predator-free island and mainland sanctuaries. We present the first mitochondrial genomes for New Zealand diplodactylid geckos, recovered from 19 modern, six historical/archival (1898-2011) and 16 Holocene Hoplodactylus duvaucelii sensu latu specimens, and one modern Woodworthia sp. specimen. No obvious damage was observed in post-extraction micro-computed tomography reconstructions. All "large gecko" specimens examined from extinct populations were found to be conspecific with extant Hoplodactylus species, suggesting their large relative size evolved only once in the New Zealand diplodactylid radiation. Phylogenetic analyses of Hoplodactylus samples recovered two genetically (and morphologically) distinct North and South Island clades, probably corresponding to distinct species. Finer phylogeographical structuring within Hoplodactylus spp. highlighted the impacts of Late Cenozoic biogeographical barriers, including the opening and closure of Pliocene marine straits, fluctuations in the size and suitability of glacial refugia, and eustatic sea-level change. Recent mainland extinction obscured these signals from the modern tissue-derived data. These results highlight the utility of minimally destructive DNA extraction in genomic analyses of less well studied small vertebrate taxa, and the conservation of natural history collections.

Facilitating accessible, rapid, and appropriate processing of ancient metagenomic data with AMDirT.

Background: Access to sample-level metadata is important when selecting public metagenomic sequencing datasets for reuse in new biological analyses. The Standards, Precautions, and Advances in Ancient Metagenomics community (SPAAM, https://spaam-community.org) has previously published AncientMetagenomeDir, a collection of curated and standardised sample metadata tables for metagenomic and microbial genome datasets generated from ancient samples. However, while sample-level information is useful for identifying relevant samples for inclusion in new projects, Next Generation Sequencing (NGS) library construction and sequencing metadata are also essential for appropriately reprocessing ancient metagenomic data. Currently, recovering information for downloading and preparing such data is difficult when laboratory and bioinformatic metadata is heterogeneously recorded in prose-based publications. Methods: Through a series of community-based hackathon events, AncientMetagenomeDir was updated to provide standardised library-level metadata of existing and new ancient metagenomic samples. In tandem, the companion tool 'AMDirT' was developed to facilitate rapid data filtering and downloading of ancient metagenomic data, as well as improving automated metadata curation and validation for AncientMetagenomeDir. Results: AncientMetagenomeDir was extended to include standardised metadata of over 6000 ancient metagenomic libraries. The companion tool 'AMDirT' provides both graphical- and command-line interface based access to such metadata for users from a wide range of computational backgrounds. We also report on errors with metadata reporting that appear to commonly occur during data upload and provide suggestions on how to improve the quality of data sharing by the community. Conclusions: Together, both standardised metadata reporting and tooling will help towards easier incorporation and reuse of public ancient metagenomic datasets into future analyses.

Population dynamics and demographic history of Eurasian collared lemmings.

BACKGROUND: Ancient DNA studies suggest that Late Pleistocene climatic changes had a significant effect on population dynamics in Arctic species. The Eurasian collared lemming (Dicrostonyx torquatus) is a keystone species in the Arctic ecosystem. Earlier studies have indicated that past climatic fluctuations were important drivers of past population dynamics in this species. RESULTS: Here, we analysed 59 ancient and 54 modern mitogenomes from across Eurasia, along with one modern nuclear genome. Our results suggest population growth and genetic diversification during the early Late Pleistocene, implying that collared lemmings may have experienced a genetic bottleneck during the warm Eemian interglacial. Furthermore, we find multiple temporally structured mitogenome clades during the Late Pleistocene, consistent with earlier results suggesting a dynamic late glacial population history. Finally, we identify a population in northeastern Siberia that maintained genetic diversity and a constant population size at the end of the Pleistocene, suggesting suitable conditions for collared lemmings in this region during the increasing temperatures associated with the onset of the Holocene. CONCLUSIONS: This study highlights an influence of past warming, in particular the Eemian interglacial, on the evolutionary history of the collared lemming, along with spatiotemporal population structuring throughout the Late Pleistocene.

An extinct and deeply divergent tiger lineage from northeastern China recognized through palaeogenomics.

Tigers (Panthera tigris) are flagship big cats and attract extensive public attention due to their charismatic features and endangered status. Despite this, little is known about their prehistoric lineages and detailed evolutionary histories. Through palaeogenomic analyses, we identified a Pleistocene tiger from northeastern China, dated to beyond the limits of radiocarbon dating (greater than 43 500 years ago). We used a simulated dataset and different reads processing pipelines to test the validity of our results and confirmed that, in both mitochondrial and nuclear phylogenies, this ancient individual belongs to a previously unknown lineage that diverged prior to modern tiger diversification. Based on the mitochondrial genome, the divergence time of this ancient lineage was estimated to be approximately 268 ka (95% CI: 187-353 ka), doubling the known age of tigers' maternal ancestor to around 125 ka (95% CI: 88-168 ka). Furthermore, by combining our findings with putative mechanisms underlying the discordant mito-nuclear phylogenetic placement for the South China tigers, we proposed a more complex scenario of tiger evolution that would otherwise be missed using data from modern tigers only. Our study provides the first glimpses of the genetic antiquity of tigers and demonstrates the utility of aDNA-based investigation for further understanding tiger evolution.

Paleogenomics reveals independent and hybrid origins of two morphologically distinct wolf lineages endemic to Japan.

Little is known about the spatiotemporal dynamics of gray wolves in the Pleistocene across low-latitude regions of Eurasia. In Japan, a small-bodied endemic subspecies of Japanese wolves existed and went extinct in the early 1900s. The fossil record indicates that a giant wolf, which reached 70 cm in body height, inhabited Japan during the Pleistocene, but its evolutionary relationship, if any, with the Japanese wolf remains uncertain. Here, to reveal the genetic origin of the Japanese wolf, we analyzed ancient DNA from remains (recovered in Japan) of one Pleistocene wolf that lived 35,000 years ago and one Holocene wolf from 5,000 years ago. The analysis of the mitochondrial DNA revealed that the Pleistocene wolf was not part of the Japanese wolf clade but rather an earlier-diverging lineage. The analysis of the nuclear DNA of the Holocene Japanese wolf revealed that it was an admixture of the Japanese Pleistocene wolf and continental wolf lineages. These findings suggest that the Japanese wolf originated via waves of colonization of multiple Pleistocene wolf populations at 57-35 and 37-14 ka, respectively, followed by interpopulation hybridization.

Modern Siberian dog ancestry was shaped by several thousand years of Eurasian-wide trade and human dispersal.

Dogs have been essential to life in the Siberian Arctic for over 9,500 y, and this tight link between people and dogs continues in Siberian communities. Although Arctic Siberian groups such as the Nenets received limited gene flow from neighboring groups, archaeological evidence suggests that metallurgy and new subsistence strategies emerged in Northwest Siberia around 2,000 y ago. It is unclear if the Siberian Arctic dog population was as continuous as the people of the region or if instead admixture occurred, possibly in relation to the influx of material culture from other parts of Eurasia. To address this question, we sequenced and analyzed the genomes of 20 ancient and historical Siberian and Eurasian Steppe dogs. Our analyses indicate that while Siberian dogs were genetically homogenous between 9,500 to 7,000 y ago, later introduction of dogs from the Eurasian Steppe and Europe led to substantial admixture. This is clearly the case in the Iamal-Nenets region (Northwestern Siberia) where dogs from the Iron Age period (∼2,000 y ago) possess substantially less ancestry related to European and Steppe dogs than dogs from the medieval period (∼1,000 y ago). Combined with findings of nonlocal materials recovered from these archaeological sites, including glass beads and metal items, these results indicate that Northwest Siberian communities were connected to a larger trade network through which they acquired genetically distinctive dogs from other regions. These exchanges were part of a series of major societal changes, including the rise of large-scale reindeer pastoralism ∼800 y ago.

How a Paleogenomic Approach Can Provide Details on Bioarchaeological Reconstruction: A Case Study from the Globular Amphorae Culture.

Ancient human remains have the potential to explain a great deal about the prehistory of humankind. Due to recent technological and bioinformatics advances, their study, at the palaeogenomic level, can provide important information about population dynamics, culture changes, and the lifestyles of our ancestors. In this study, mitochondrial and nuclear genome data obtained from human bone remains associated with the Neolithic Globular Amphorae culture, which were recovered in the Megalithic barrow of Kierzkowo (Poland), were reanalysed to gain insight into the social organisation and use of the archaeological site and to provide information at the individual level. We were able to successfully estimate the minimum number of individuals, sex, kin relationships, and phenotypic traits of the buried individuals, despite the low level of preservation of the bone samples and the intricate taphonomic conditions. In addition, the evaluation of damage patterns allowed us to highlight the presence of "intruders"-that is, of more recent skeletal remains that did not belong to the original burial. Due to its characteristics, the study of the Kierzkowo barrow represented a challenge for the reconstruction of the biological profile of the human community who exploited it and an excellent example of the contribution that ancient genomic analysis can provide to archaeological reconstruction.

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.

Genome of Pestera Muierii skull shows high diversity and low mutational load in pre-glacial Europe.

Few complete human genomes from the European Early Upper Palaeolithic (EUP) have been sequenced. Using novel sampling and DNA extraction approaches, we sequenced the genome of a woman from "Pestera Muierii," Romania who lived ∼34,000 years ago to 13.5× coverage. The genome shows similarities to modern-day Europeans, but she is not a direct ancestor. Although her cranium exhibits both modern human and Neanderthal features, the genome shows similar levels of Neanderthal admixture (∼3.1%) to most EUP humans but only half compared to the ∼40,000-year-old Pestera Oase 1. All EUP European hunter-gatherers display high genetic diversity, demonstrating that the severe loss of diversity occurred during and after the Last Glacial Maximum (LGM) rather than just during the out-of-Africa migration. The prevalence of genetic diseases is expected to increase with low diversity; however, pathogenic variant load was relatively constant from EUP to modern times, despite post-LGM hunter-gatherers having the lowest diversity ever observed among Europeans.

Reproducible, portable, and efficient ancient genome reconstruction with nf-core/eager.

The broadening utilisation of ancient DNA to address archaeological, palaeontological, and biological questions is resulting in a rising diversity in the size of laboratories and scale of analyses being performed. In the context of this heterogeneous landscape, we present an advanced, and entirely redesigned and extended version of the EAGER pipeline for the analysis of ancient genomic data. This Nextflow pipeline aims to address three main themes: accessibility and adaptability to different computing configurations, reproducibility to ensure robust analytical standards, and updating the pipeline to the latest routine ancient genomic practices. The new version of EAGER has been developed within the nf-core initiative to ensure high-quality software development and maintenance support; contributing to a long-term life-cycle for the pipeline. nf-core/eager will assist in ensuring that a wider range of ancient DNA analyses can be applied by a diverse range of research groups and fields.

Genomes of Pleistocene Siberian Wolves Uncover Multiple Extinct Wolf Lineages.

Extant Canis lupus genetic diversity can be grouped into three phylogenetically distinct clades: Eurasian and American wolves and domestic dogs.1 Genetic studies have suggested these groups trace their origins to a wolf population that expanded during the last glacial maximum (LGM)1-3 and replaced local wolf populations.4 Moreover, ancient genomes from the Yana basin and the Taimyr peninsula provided evidence of at least one extinct wolf lineage that dwelled in Siberia during the Pleistocene.35 Previous studies have suggested that Pleistocene Siberian canids can be classified into two groups based on cranial morphology. Wolves in the first group are most similar to present-day populations, although those in the second group possess intermediate features between dogs and wolves.67 However, whether this morphological classification represents distinct genetic groups remains unknown. To investigate this question and the relationships between Pleistocene canids, present-day wolves, and dogs, we resequenced the genomes of four Pleistocene canids from Northeast Siberia dated between >50 and 14 ka old, including samples from the two morphological categories. We found these specimens cluster with the two previously sequenced Pleistocene wolves, which are genetically more similar to Eurasian wolves. Our results show that, though the four specimens represent extinct wolf lineages, they do not form a monophyletic group. Instead, each Pleistocene Siberian canid branched off the lineage that gave rise to present-day wolves and dogs. Finally, our results suggest the two previously described morphological groups could represent independent lineages similarly related to present-day wolves and dogs.

Competitive mapping allows for the identification and exclusion of human DNA contamination in ancient faunal genomic datasets.

BACKGROUND: After over a decade of developments in field collection, laboratory methods and advances in high-throughput sequencing, contamination remains a key issue in ancient DNA research. Currently, human and microbial contaminant DNA still impose challenges on cost-effective sequencing and accurate interpretation of ancient DNA data. RESULTS: Here we investigate whether human contaminating DNA can be found in ancient faunal sequencing datasets. We identify variable levels of human contamination, which persists even after the sequence reads have been mapped to the faunal reference genomes. This contamination has the potential to affect a range of downstream analyses. CONCLUSIONS: We propose a fast and simple method, based on competitive mapping, which allows identifying and removing human contamination from ancient faunal DNA datasets with limited losses of true ancient data. This method could represent an important tool for the ancient DNA field.

Rapid size change associated with intra-island evolutionary radiation in extinct Caribbean "island-shrews".

BACKGROUND: The Caribbean offers a unique opportunity to study evolutionary dynamics in insular mammals. However, the recent extinction of most Caribbean non-volant mammals has obstructed evolutionary studies, and poor DNA preservation associated with tropical environments means that very few ancient DNA sequences are available for extinct vertebrates known from the region's Holocene subfossil record. The endemic Caribbean eulipotyphlan family Nesophontidae ("island-shrews") became extinct ~ 500 years ago, and the taxonomic validity of many Nesophontes species and their wider evolutionary dynamics remain unclear. Here we use both morphometric and palaeogenomic methods to clarify the status and evolutionary history of Nesophontes species from Hispaniola, the second-largest Caribbean island. RESULTS: Principal component analysis of 65 Nesophontes mandibles from late Quaternary fossil sites across Hispaniola identified three non-overlapping morphometric clusters, providing statistical support for the existence of three size-differentiated Hispaniolan Nesophontes species. We were also able to extract and sequence ancient DNA from a ~ 750-year-old specimen of Nesophontes zamicrus, the smallest non-volant Caribbean mammal, including a whole-mitochondrial genome and partial nuclear genes. Nesophontes paramicrus (39-47 g) and N. zamicrus (~ 10 g) diverged recently during the Middle Pleistocene (mean estimated divergence = 0.699 Ma), comparable to the youngest species splits in Eulipotyphla and other mammal groups. Pairwise genetic distance values for N. paramicrus and N. zamicrus based on mitochondrial and nuclear genes are low, but fall within the range of comparative pairwise data for extant eulipotyphlan species-pairs. CONCLUSIONS: Our combined morphometric and palaeogenomic analyses provide evidence for multiple co-occurring species and rapid body size evolution in Hispaniolan Nesophontes, in contrast to patterns of genetic and morphometric differentiation seen in Hispaniola's extant non-volant land mammals. Different components of Hispaniola's mammal fauna have therefore exhibited drastically different rates of morphological evolution. Morphological evolution in Nesophontes is also rapid compared to patterns across the Eulipotyphla, and our study provides an important new example of rapid body size change in a small-bodied insular vertebrate lineage. The Caribbean was a hotspot for evolutionary diversification as well as preserving ancient biodiversity, and studying the surviving representatives of its mammal fauna is insufficient to reveal the evolutionary patterns and processes that generated regional diversity.

Three-dimensional characterization of osteocyte volumes at multiple scales, and its relationship with bone biology and genome evolution in ray-finned fishes.

Osteocytes, cells embedded within the bone mineral matrix, inform on key aspects of vertebrate biology. In particular, a relationship between volumes of the osteocytes and bone growth and/or genome size has been proposed for several tetrapod lineages. However, the variation in osteocyte volume across different scales is poorly characterized and mostly relies on incomplete, two-dimensional information. In this study, we characterize the variation of osteocyte volumes in ray-finned fishes (Actinopterygii), a clade including more than half of modern vertebrate species in which osteocyte biology is poorly known. We use X-ray synchrotron micro-computed tomography (SRµCT) to achieve a three-dimensional visualization of osteocyte lacunae and direct measurement of their size (volumes). Our specimen sample is designed to characterize variation in osteocyte lacuna morphology at three scales: within a bone, among the bones of one individual and among species. At the intra-bone scale, we find that osteocyte lacunae vary noticeably in size between zones of organized and woven bone (being up to six times larger in woven bone), and across cyclical bone deposition. This is probably explained by differences in bone deposition rate, with larger osteocyte lacunae contained in bone that deposits faster. Osteocyte lacuna volumes vary 3.5-fold among the bones of an individual, and this cannot readily be explained by variation in bone growth rate or other currently observable factors. Finally, we find that genome size provides the best explanation of variation in osteocyte lacuna volume among species: actinopterygian taxa with larger genomes (polyploid taxa in particular) have larger osteocyte lacunae (with a ninefold variation in median osteocyte volume being measured). Our findings corroborate previous two-dimensional studies in tetrapods that also observed similar patterns of intra-individual variation and found a correlation with genome size. This opens new perspectives for further studies on bone evolution, physiology and palaeogenomics in actinopterygians, and vertebrates as a whole.

Consensify: A Method for Generating Pseudohaploid Genome Sequences from Palaeogenomic Datasets with Reduced Error Rates.

A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes.