Application of palaeogenetic techniques to historic mollusc shells reveals phylogeographic structure in a New Zealand abalone.

Natural history collections worldwide contain a plethora of mollusc shells. Recent studies have detailed the sequencing of DNA extracted from shells up to thousands of years old and from various taphonomic and preservational contexts. However, previous approaches have largely addressed methodological rather than evolutionary research questions. Here, we report the generation of DNA sequence data from mollusc shells using such techniques, applied to Haliotis virginea Gmelin, 1791, a New Zealand abalone, in which morphological variation has led to the recognition of several forms and subspecies. We successfully recovered near-complete mitogenomes from 22 specimens including 12 dry-preserved shells up to 60 years old. We used a combination of palaeogenetic techniques that have not previously been applied to shell, including DNA extraction optimized for ultra-short fragments and hybridization-capture of single-stranded DNA libraries. Phylogenetic analyses revealed three major, well-supported clades comprising samples from: (1) The Three Kings Islands; (2) the Auckland, Chatham and Antipodes Islands; and (3) mainland New Zealand and Campbell Island. This phylogeographic structure does not correspond to the currently recognized forms. Critically, our nonreliance on freshly collected or ethanol-preserved samples enabled inclusion of topotypes of all recognized subspecies as well as additional difficult-to-sample populations. Broader application of these comparatively cost-effective and reliable methods to modern, historical, archaeological and palaeontological shell samples has the potential to revolutionize invertebrate genetic research.

Genetic evidence for post-glacial expansion from a southern refugium in the eastern moa (Emeus crassus).

Cycles of glacial expansion and contraction throughout the Pleistocene drove increases and decreases, respectively, in the geographical range and population size of many animal species. Genetic data have revealed that during glacial maxima the distribution of many Eurasian animals was restricted to small refugial areas, from which species expanded to reoccupy parts of their former range as the climate warmed. It has been suggested that the extinct eastern moa (Emeus crassus)-a large, flightless bird from New Zealand-behaved analogously during glacial maxima, possibly surviving only in a restricted area of lowland habitat in the southern South Island of New Zealand during the Last Glacial Maximum (LGM). However, previous studies have lacked the power and geographical sampling to explicitly test this hypothesis using genetic data. Here we analyse 46 ancient mitochondrial genomes from Late Pleistocene and Holocene bones of the eastern moa from across their post-LGM distribution. Our results are consistent with a post-LGM increase in the population size and genetic diversity of eastern moa. We also demonstrate that genetic diversity was higher in eastern moa from the southern extent of their range, supporting the hypothesis that they expanded from a single glacial refugium following the LGM.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge.

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Mitogenomes Reveal Two Major Influxes of Papuan Ancestry across Wallacea Following the Last Glacial Maximum and Austronesian Contact.

The tropical archipelago of Wallacea contains thousands of individual islands interspersed between mainland Asia and Near Oceania, and marks the location of a series of ancient oceanic voyages leading to the peopling of Sahul-i.e., the former continent that joined Australia and New Guinea at a time of lowered sea level-by 50,000 years ago. Despite the apparent deep antiquity of human presence in Wallacea, prior population history research in this region has been hampered by patchy archaeological and genetic records and is largely concentrated upon more recent history that follows the arrival of Austronesian seafarers ~3000-4000 years ago (3-4 ka). To shed light on the deeper history of Wallacea and its connections with New Guinea and Australia, we performed phylogeographic analyses on 656 whole mitogenomes from these three regions, including 186 new samples from eight Wallacean islands and three West Papuan populations. Our results point to a surprisingly dynamic population history in Wallacea, marked by two periods of extensive demographic change concentrated around the Last Glacial Maximum ~15 ka and post-Austronesian contact ~3 ka. These changes appear to have greatly diminished genetic signals informative about the original peopling of Sahul, and have important implications for our current understanding of the population history of the region.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote.

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

Revision of the Litoria watjulumensis (Anura: Pelodryadidae) group from the Australian monsoonal tropics, including the resurrection of L. spaldingi.

We show that the Wotjulum frog, Litoria watjulumensis (Copland, 1957), comprises two deeply divergent mitochondrial DNA lineages that are also reciprocally monophyletic for a nuclear gene locus and have discrete distributions. The taxa are differentiated in multivariate analysis of shape but show no appreciable differences in colour and pattern. The two taxa differ substantially in the degree of female biased sexual size dimorphism, with the western taxon showing considerably more pronounced dimorphism. We subsequently resurrect Litoria (Hyla) spaldingi (Hosmer, 1964) for populations from east of the Daly River system in the Northern Territory through to western Queensland and restrict L. watjulumensis to populations from the Kimberley region of north-western Australia and the Victoria River system of the western Northern Territory. The complex advertisement call of L. spaldingi is described for the first time.

Examining Natural History through the Lens of Palaeogenomics.

The many high-resolution tools that are uniquely applicable to specimens from the Quaternary period (the past ~2.5 Ma) provide an opportunity to cross-validate data and test hypotheses based on the morphology and distribution of fossils. Among these tools is palaeogenomics - the genome-scale sequencing of genetic material from ancient specimens - that can provide direct insight into ecology and evolution, potentially improving the accuracy of inferences about past ecological communities over longer timescales. Palaeogenomics has revealed instances of over- and underestimation of extinct diversity, detected cryptic faunal migration and turnover, allowed quantification of widespread sex biases and sexual dimorphism in the fossil record, revealed past hybridisation events and hybrid individuals, and has highlighted previously unrecognised routes of zoonotic disease transfer.

Dire wolves were the last of an ancient New World canid lineage.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently.

Oligocene divergence of frogmouth birds (Podargidae) across Wallace's Line.

Wallace's Line demarcates the transition between the differentiated regional faunas of Asia and Australia. However, while patterns of biotic differentiation across these two continental landmasses and the intervening island groups (Wallacea) have been extensively studied, patterns of long-term dispersal and diversification across this region are less well understood. Frogmouths (Aves: Podargidae) are a relictual family of large nocturnal birds represented by three extant genera occurring, respectively, in Asia, 'Sahul' (Australia and New Guinea) and the Solomon Islands, thus spanning Wallace's Line. We used new mitochondrial genomes from each of the extant frogmouth genera to estimate the timeline of frogmouth evolution and dispersal across Wallace's Line. Our results suggest that the three genera diverged and dispersed during the mid-Cenozoic between approximately 30 and 40 Mya. These divergences are among the oldest inferred for any trans-Wallacean vertebrate lineage. In addition, our results reveal that the monotypic Solomons frogmouth (Rigidipenna inexpectata) is one of the most phylogenetically divergent endemic bird lineages in the southwest Pacific. We suggest that the contemporary distribution of exceptionally deep divergences among extant frogmouth lineages may be explained by colonization of, and subsequent long-term persistence on, island arcs in the southwest Pacific during the Oligocene. These island arcs may have provided a pathway for biotic dispersal out of both Asia and Australia that preceded the formation of extensive emergent landmasses in Wallacea by at least 10 million years.

The Dogma of Dingoes-Taxonomic status of the dingo: A reply to Smith et al.

Adopting the name Canis dingo for the Dingo to explicitly denote a species-level taxon separate from other canids was suggested by Crowther et al.  (2014) as a means to eliminate taxonomic instability and contention. However, Jackson et al.  (2017), using standard taxonomic and nomenclatural approaches and principles, called instead for continued use of the nomen C. familiaris for all domestic dogs and their derivatives, including the Dingo. (This name, C. familiaris, is applied to all dogs that derive from the domesticated version of the Gray Wolf, Canis lupus, based on nomenclatural convention.) The primary reasons for this call by Jackson et al.  (2017) were: (1) a lack of evidence to show that recognizing multiple species amongst the dog, including the Dingo and New Guinea Singing Dog, was necessary taxonomically, and (2) the principle of nomenclatural priority (the name familiaris Linnaeus, 1758, antedates dingo Meyer, 1793). Overwhelming current evidence from archaeology and genomics indicates that the Dingo is of recent origin in Australia and shares immediate ancestry with other domestic dogs as evidenced by patterns of genetic and morphological variation. Accordingly, for Smith et al.  (2019) to recognise Canis dingo as a distinct species, the onus was on them to overturn current interpretations of available archaeological, genomic, and morphological datasets and instead show that Dingoes have a deeply divergent evolutionary history that distinguishes them from other named forms of Canis (including C. lupus and its domesticated version, C. familiaris). A recent paper by Koepfli et al.  (2015) demonstrates exactly how this can be done in a compelling way within the genus Canis-by demonstrating deep evolutionary divergence between taxa, on the order of hundreds of thousands of years, using data from multiple genetic systems. Smith et al.  (2019) have not done this; instead they have misrepresented the content and conclusions of Jackson et al.  (2017), and contributed extraneous arguments that are not relevant to taxonomic decisions. Here we dissect Smith et al.  (2019), identifying misrepresentations, to show that ecological, behavioural and morphological evidence is insufficient to recognise Dingoes as a separate species from other domestic dogs. We reiterate: the correct binomial name for the taxon derived from Gray Wolves (C. lupus) by passive and active domestication, including Dingoes and other domestic dogs, is Canis familiaris. We are strongly sympathetic to arguments about the historical, ecological, cultural, or other significance of the Dingo, but these are issues that will have to be considered outside of the more narrow scope of taxonomy and nomenclature.

Widespread male sex bias in mammal fossil and museum collections.

A recent study of mammoth subfossil remains has demonstrated the potential of using relatively low-coverage high-throughput DNA sequencing to genetically sex specimens, revealing a strong male-biased sex ratio [P. Pecnerová et al., Curr. Biol. 27, 3505-3510.e3 (2017)]. Similar patterns were predicted for steppe bison, based on their analogous female herd-based structure. We genetically sexed subfossil remains of 186 Holarctic bison (Bison spp.), and also 91 brown bears (Ursus arctos), which are not female herd-based, and found that ∼75% of both groups were male, very close to the ratio observed in mammoths (72%). This large deviation from a 1:1 ratio was unexpected, but we found no evidence for sex differences with respect to DNA preservation, sample age, material type, or overall spatial distribution. We further examined ratios of male and female specimens from 4 large museum mammal collections and found a strong male bias, observable in almost all mammalian orders. We suggest that, in mammals at least, 1) wider male geographic ranges can lead to considerably increased chances of detection in fossil studies, and 2) sexual dimorphic behavior or appearance can facilitate a considerable sex bias in fossil and modern collections, on a previously unacknowledged scale. This finding has major implications for a wide range of studies of fossil and museum material.

Low-cost cross-taxon enrichment of mitochondrial DNA using in-house synthesised RNA probes.

Hybridization capture with in-solution oligonucleotide probes has quickly become the preferred method for enriching specific DNA loci from degraded or ancient samples prior to high-throughput sequencing (HTS). Several companies synthesize sets of probes for in-solution hybridization capture, but these commercial reagents are usually expensive. Methods for economical in-house probe synthesis have been described, but they do not directly address one of the major advantages of commercially synthesised probes: that probe sequences matching many species can be synthesised in parallel and pooled. The ability to make "phylogenetically diverse" probes increases the cost-effectiveness of commercial probe sets, as they can be used across multiple projects (or for projects involving multiple species). However, it is labour-intensive to replicate this with in-house methods, as template molecules must first be generated for each species of interest. While it has been observed that probes can be used to enrich for phylogenetically distant targets, the ability of this effect to compensate for the lack of phylogenetically diverse probes in in-house synthesised probe sets has not been tested. In this study, we present a refined protocol for in-house RNA probe synthesis and evaluated the ability of probes generated using this method from a single species to successfully enrich for the target locus in phylogenetically distant species. We demonstrated that probes synthesized using long-range PCR products from a placental mammal mitochondrion (Bison spp.) could be used to enrich for mitochondrial DNA in birds and marsupials (but not plants). Importantly, our results were obtained for approximately a third of the cost of similar commercially available reagents.

Mitogenomes Uncover Extinct Penguin Taxa and Reveal Island Formation as a Key Driver of Speciation.

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Correction: Low-cost cross-taxon enrichment of mitochondrial DNA using in-house synthesised RNA probes.

[This corrects the article DOI: 10.1371/journal.pone.0209499.].

Reconstructing the Evolution of Giant Extinct Kangaroos: Comparing the Utility of DNA, Morphology, and Total Evidence.

Combined "total evidence" analysis of molecular and morphological data offers the opportunity to objectively merge fossils into the tree of life, and challenges the primacy of solely DNA based phylogenetic and dating inference, even among modern taxa. To investigate the relative utility of DNA, morphology, and total evidence for evolutionary inference, we sequenced the first near-complete mitochondrial genomes from extinct Australian megafauna: a 40-50 thousand year old giant short-faced kangaroo (Simosthenurus occidentalis) and giant wallaby (Protemnodon anak). We analyzed the ancient DNA and fossil data alongside comparable data from extant species to infer phylogeny, divergence times, and ancestral body mass among macropods (kangaroos and wallabies). Our results confirm a close relationship between Protemnodon and the iconic kangaroo genus complex "Macropus", and unite the giant Simothenurus with the hare-sized Lagostrophus fasciatus (banded hare-wallaby), suggesting that the latter is the closest living link to the once diverse sthenurine kangaroo radiation. We find that large body size evolved multiple times among kangaroos, coincident with expansion of open woodland habitats beginning in the Late Miocene. In addition, our results suggest that morphological data mislead macropod phylogeny reconstruction and in turn can distort total evidence estimation of divergence dates. However, a novel result with potentially broad application is that the accuracy and precision of reconstructing ancestral body mass was improved by tracing body mass on morphological branch lengths. This is likely due to positive allometric correlation between morphological and body size variation-a relationship that may be masked or even misleadingly inverted with the temporal or molecular branch lengths that typically underpin ancestral body size reconstruction. Our study supports complementary roles for DNA and morphology in evolutionary inference, and opens a new window into the evolution of Australia's unique marsupial fauna.

Evolution and extinction of the giant rhinoceros Elasmotherium sibiricum sheds light on late Quaternary megafaunal extinctions.

Understanding extinction events requires an unbiased record of the chronology and ecology of victims and survivors. The rhinoceros Elasmotherium sibiricum, known as the 'Siberian unicorn', was believed to have gone extinct around 200,000 years ago-well before the late Quaternary megafaunal extinction event. However, no absolute dating, genetic analysis or quantitative ecological assessment of this species has been undertaken. Here, we show, by accelerator mass spectrometry radiocarbon dating of 23 individuals, including cross-validation by compound-specific analysis, that E. sibiricum survived in Eastern Europe and Central Asia until at least 39,000 years ago, corroborating a wave of megafaunal turnover before the Last Glacial Maximum in Eurasia, in addition to the better-known late-glacial event. Stable isotope data indicate a dry steppe niche for E. sibiricum and, together with morphology, a highly specialized diet that probably contributed to its extinction. We further demonstrate, with DNA sequencing data, a very deep phylogenetic split between the subfamilies Elasmotheriinae and Rhinocerotinae that includes all the living rhinoceroses, settling a debate based on fossil evidence and confirming that the two lineages had diverged by the Eocene. As the last surviving member of the Elasmotheriinae, the demise of the 'Siberian unicorn' marked the extinction of this subfamily.

Molecular resolution to a morphological controversy: The case of North American fossil muskoxen Bootherium and Symbos.

The musk ox (Ovibos moschatus) is the only surviving member of a group of Pleistocene North American musk ox genera (Praeovibos, Ovibos, Bootherium, Euceratherium, and Soergelia) whose taxonomy is uncertain. The helmeted musk ox (Bootherium bombifrons) and the woodland musk ox (Symbos cavifrons) have been synonymised as male and female forms of a single Nearctic species found from Alaska, in the north, to Texas, in the south. However, this reclassification has not been tested using molecular data, despite the potential to use ancient DNA to examine these late Pleistocene taxa. In the present study, we sequenced mitochondrial genomes from seven subfossil musk ox specimens (originally identified as Bootherium and/or Symbos), allowing us to evaluate the identity of these muskoxen, explore their phylogeography, and estimate the timeline for their evolution. We also used nuclear genomic data to determine the sex of six of our seven samples. Ultimately, our molecular data support the synonymisation of the North American muskoxen Bootherium and Symbos.

Genetic diversity and drivers of dwarfism in extinct island emu populations.

Australia's iconic emu (Dromaius novaehollandiae novaehollandiae) is the only living representative of its genus, but fossil evidence and reports from early European explorers suggest that three island forms (at least two of which were dwarfs) became extinct during the nineteenth century. While one of these-the King Island emu-has been found to be conspecific with Australian mainland emus, little is known about how the other two forms-Kangaroo Island and Tasmanian emus-relate to the others, or even the size of Tasmanian emus. We present a comprehensive genetic and morphological analysis of Dromaius diversity, including data from one of the few definitively genuine Tasmanian emu specimens known. Our genetic analyses suggest that all the island populations represent sub-populations of mainland Dnovaehollandiae Further, the size of island emus and those on the mainland appears to scale linearly with island size but not time since isolation, suggesting that island size-and presumably concomitant limitations on resource availability-may be a more important driver of dwarfism in island emus, though its precise contribution to emu dwarfism remains to be confirmed.