Search

The extinction of the woolly rhinoceros (Coelodonta antiquitatis) at the onset of the Holocene remains an enigma, with conflicting evidence regarding its cause and spatiotemporal dynamics. This partly reflects challenges in determining demographic responses of late Quaternary megafauna to climatic and anthropogenic causal drivers with available genetic and paleontological techniques. Here, we show that elucidating mechanisms of ancient extinctions can benefit from a detailed understanding of fine-scale metapopulation dynamics, operating over many millennia. Using an abundant fossil record, ancient DNA, and high-resolution simulation models, we untangle the ecological mechanisms and causal drivers that are likely to have been integral in the decline and later extinction of the woolly rhinoceros. Our 52,000-y reconstruction of distribution-wide metapopulation dynamics supports a pathway to extinction that began long before the Holocene, when the combination of cooling temperatures and low but sustained hunting by humans trapped woolly rhinoceroses in suboptimal habitats along the southern edge of their range. Modeling indicates that this ecological trap intensified after the end of the last ice age, preventing colonization of newly formed suitable habitats, weakening stabilizing metapopulation processes, triggering the extinction of the woolly rhinoceros in the early Holocene. Our findings suggest that fragmentation and resultant metapopulation dynamics should be explicitly considered in explanations of late Quaternary megafauna extinctions, sending a clarion call to the fragility of the remaining large-bodied grazers restricted to disjunct fragments of poor-quality habitat due to anthropogenic environmental change. [ABSTRACT FROM AUTHOR]

Due to their limited dispersal ability, fossorial species with predominantly belowground activity usually show increased levels of population subdivision across relatively small spatial scales. This may be exacerbated in harsh mountain ecosystems, where landscape geomorphology limits species' dispersal ability and leads to small effective population sizes, making species relatively vulnerable to environmental change. To better understand the environmental drivers of species' population subdivision in remote mountain ecosystems, particularly in understudied high‐elevation systems in Africa, we studied the giant root‐rat (Tachyoryctes macrocephalus), a fossorial rodent confined to the afro‐alpine ecosystem of the Bale Mountains in Ethiopia. Using mitochondrial and low‐coverage nuclear genomes, we investigated 77 giant root‐rat individuals sampled from nine localities across its entire ~1000 km2 range. Our data revealed a distinct division into a northern and southern group, with no signs of gene flow, and higher nuclear genetic diversity in the south. Landscape genetic analyses of the mitochondrial and nuclear genomes indicated that population subdivision was driven by slope and elevation differences of up to 500 m across escarpments separating the north and south, potentially reinforced by glaciation of the south during the Late Pleistocene (~42,000–16,000 years ago). Despite this landscape‐scale subdivision between the north and south, weak geographic structuring of sampling localities within regions indicated gene flow across distances of at least 16 km at the local scale, suggesting high, aboveground mobility for relatively long distances. Our study highlights that despite the potential for local‐scale gene flow in fossorial species, topographic barriers can result in pronounced genetic subdivision. These factors can reduce genetic variability, which should be considered when developing conservation strategies. [ABSTRACT FROM AUTHOR]

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria–Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries. [ABSTRACT FROM AUTHOR]

The impact of post‐divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide‐ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post‐divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding‐window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D‐statistics and f‐branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools. [ABSTRACT FROM AUTHOR]

Stable carbon (δ13C) and nitrogen (δ15N) isotopic compositions of bone and dentine collagen extracted from museum specimens have been widely used to study the paleoecology of past populations. Due to possible systematic differences in stable isotope values between bone and dentine, dentine values need to be transformed into bone-collagen equivalent using a correction factor to allow comparisons between the two collagen sources. Here, we provide correction factors to transform dentine δ13C and δ15N values into bone-collagen equivalent for two toothed whales: narwhal and beluga. We sampled bone and dentine from the skulls of 11 narwhals and 26 belugas. In narwhals, dentine was sampled from tusk and embedded tooth; in belugas, dentine was sampled from tooth. δ13C and δ15N were measured, and intra-individual bone and dentine isotopic compositions were used to calculate correction factors for each species. We detected differences in δ13C and δ15N. In both narwhals and belugas, we found lower average δ13C and δ15N in bone compared with dentine. The correction factors provided by the study enable the combined analysis of stable isotope data from bone and dentine in these species. [ABSTRACT FROM AUTHOR]

The Saiga antelope (Saiga tatarica) is one of the few megafauna species from the mammoth steppe still living today. Currently, saiga are classified as critically endangered, persisting only in small areas of Central Asian steppe and desert ecosystems. The species is divided into two subspecies: Saiga tatarica mongolica and Saiga tatarica tatarica. In this study, we have for the first time characterized the genetic diversity of the Mongolian saiga (S. t. mongolica) using both mitochondrial DNA and microsatellite markers. We also analysed S. t. tatarica specimens in order to genetically compare both subspecies. The mitochondrial control region was sequenced for a total of 89 individuals: 20 skin, 53 umbilical cord, three placentae and a muscle sample from S. t. mongolica , and a total of 12 hair samples from S. t. tatarica. Additionally, 19 microsatellites developed for saiga antelope were also screened. Our results revealed that the Mongolian saiga presents very low genetic diversity at the mitochondrial level, with no shared mitochondrial haplotype between the two subspecies. Low genetic diversity is also present at the autosomal level, with most loci having low heterozygosity (Ho/He) and a low number of alleles per locus. Despite the low genetic diversity, we found no separation between the subpopulations in Mongolia, indicating that conservation corridors are actually promoting contact between different herds. Our results validate current conservation efforts and inform the implementation of new measures to increase the viability of the S. t. mongolica subspecies. [ABSTRACT FROM AUTHOR]

The first complete ancient mitochondrial genome of a vicuña (Vicugna vicugna) (3200–2400 B.P.) recovered from the Tulán-54 site (San Pedro de Atacama, Northern of Chile) was sequenced using mitochondrial enrichment by hybridization capture and subsequent Next Generation Sequencing (NGS). Morphometric analyses place this individual inside of a small size camelid group and Bayesian phylogenetic analysis confirms that this specimen belonged to an ancient vicuña. [ABSTRACT FROM AUTHOR]

The brown bear ( Ursus arctos ) is one of the survivors of the Late Quaternary megafauna extinctions. However, despite being widely distributed across the Holarctic, brown bears have experienced extensive range reductions, and even extirpations in some geographical regions. Previous research efforts using genetic data have provided valuable insights into their evolutionary history. However, most studies have been limited to contemporary individuals or mitochondrial DNA, limiting insights into population processes that preceded the present. Here, we present genomic data from two Late Pleistocene brown bears from Honshu, Japan and eastern Siberia, and combine them with published contemporary and ancient genomes from across the Holarctic range of brown bears to investigate the evolutionary relationships among brown bear populations through time and space. By including genomic data from Late Pleistocene and Holocene individuals sampled outside the current distribution range, we uncover diversity not present in contemporary populations. Notably, although contemporary individuals display geographically structured populations most likely driven by isolation-by-distance, this pattern varies among the ancient samples across different regions. The inclusion of ancient brown bears in our analysis provides novel insights into the evolutionary history of brown bears and contributes to understanding the populations and diversity lost during the Late Quaternary.

Stable carbon ( δ 13 C) and nitrogen ( δ 15 N) isotopic compositions of bone and dentine collagen extracted from museum specimens have been widely used to study the paleoecology of past populations. Due to possible systematic differences in stable isotope values between bone and dentine, dentine values need to be transformed into bone-collagen equivalent using a correction factor to allow comparisons between the two collagen sources. Here, we provide correction factors to transform dentine δ 13 C and δ 15 N values into bone-collagen equivalent for two toothed whales: narwhal and beluga. We sampled bone and dentine from the skulls of 11 narwhals and 26 belugas. In narwhals, dentine was sampled from tusk and embedded tooth; in belugas, dentine was sampled from tooth. δ 13 C and δ 15 N were measured, and intra-individual bone and dentine isotopic compositions were used to calculate correction factors for each species. We detected differences in δ 13 C and δ 15 N. In both narwhals and belugas, we found lower average δ 13 C and δ 15 N in bone compared with dentine. The correction factors provided by the study enable the combined analysis of stable isotope data from bone and dentine in these species.

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer ( Megaloceros giganteus ) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.