Widespread male sex bias in mammal fossil and museum collections

A label-free quantification method for assessing sex from modern and ancient bovine tooth enamel

Identification of the sex of modern, fossil and archaeological animal remains offers many insights into their demography, mortality profiles and domestication pathways. However, due to many-factors, sex determination of osteological remains is often problematic. To overcome this, we have developed an innovative protocol to determine an animal's sex from tooth enamel, by applying label-free quantification (LFQ) of two unique AmelY peptides 'LRYPYP' (AmelY;[M+2]2 + 404.7212 m/z) and 'LRYPYPSY' (AmelY;[M+2]2 + 529.7689 m/z) that are only present in the enamel of males. We applied this method to eight modern cattle (Bos taurus) of known sex, and correctly assigned them to sex. We then applied the same protocol to twelve archaeological Bos teeth from the Neolithic site of Beisamoun, Israel (8-th-7-th millennium BC) and determined the sex of the archaeological samples. Since teeth are usually better preserved than bones, this innovative protocol has potential to facilitate sex determination in ancient and modern bovine remains that currently cannot be sexed.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

Shrouded in history: Unveiling the ways of life of an early Muslim population in Santarém, Portugal (8th- 10th century AD)

In around 716 AD, the city of Santarém, Portugal, was conquered by the Berber and Arab armies that swept the Iberian Peninsula and went on to rule the region until the 12th century. Archaeological excavations in 2007/08 discovered an Islamic necropolis (Avenida 5 de Outubro #2-8) that appears to contain the remains of an early Muslim population in Santarém (8th- 10th century). In this study, skeletal material from 58 adult individuals was analysed for stable carbon (δ13Ccol; δ13Cap), nitrogen (δ15N) and sulphur (δ34S) isotope ratios in bones, and stable oxygen (δ18O), carbon (δ13Cen) and radiogenic strontium (87Sr/86Sr) isotopes in tooth enamel. The results of this study revealed a dietary pattern of predominantly C3-plant and domestic C3-fed herbivore consumption during adulthood (δ13Ccol and δ15N, respectively) but a higher proportion of C4-plant input during childhood (δ13Cen) for some individuals-interpreted as possible childhood consumption of millet porridge, a common practice in North Africa-in those with unorthodox burial types (Groups 1 and 2) that was not practiced in the individuals with canonical burials (Group 3). In this first mobility study of a medieval Muslim population in Portugal, δ18ODW values revealed greater heterogeneity in Groups 1 and 2, consistent with diverse origins, some in more humid regions than Santarém when compared to regional precipitation δ18O data, contrasting the more homogenous Group 3, consistent with the local precipitation δ18O range. Ancient DNA analysis conducted on three individuals revealed maternal (mtDNA) and paternal (Y-chromosome) lineages compatible with a North African origin for (at least) some of the individuals. Additionally, mobility of females in this population was higher than males, potentially resulting from a patrilocal social system, practiced in Berber and Arab communities. These results serve to offer a more detailed insight into the ancestry and cultural practices of early Muslim populations in Iberia.

Sex-biased sampling may influence Homo naledi tooth size variation

A frequent source of debate in paleoanthropology concerns the taxonomic unity of fossil assemblages, with many hominin samples exhibiting elevated levels of variation that can be interpreted as indicating the presence of multiple species. By contrast, the large assemblage of hominin fossils from the Rising Star cave system, assigned to Homo naledi, exhibits a remarkably low degree of variation for most skeletal elements. Many factors can contribute to low sample variation, including genetic drift, strong natural selection, biased sex ratios, and sampling of closely related individuals. In this study, we tested for potential sex-biased sampling in the Rising Star dental sample. We compared coefficients of variation for the H. naledi teeth to those for eight extant hominoid samples. We used a resampling procedure that generated samples from the extant taxa that matched the sample size of the fossil sample for each possible Rising Star dental sex ratio. We found that variation at four H. naledi tooth positions-I2, M1, P4, M1-is so low that the possibility that one sex is represented by few or no individuals in the sample cannot be excluded. Additional evidence is needed to corroborate this inference, such as ancient DNA or enamel proteome data, and our study design does not address other potential factors that would account for low sample variation. Nevertheless, our results highlight the importance of considering the taphonomic history of a hominin assemblage and suggest that sex-biased sampling is a plausible explanation for the low level of phenotypic variation found in some aspects of the current H. naledi assemblage.

Evidence for early domestic yak, taurine cattle, and their hybrids on the Tibetan Plateau

Domestic yak, cattle, and their hybrids are fundamental to herder survival at high altitudes on the Tibetan Plateau. However, little is known about their history. Bos remains are uncommon in this region, and ancient domestic yak have not been securely identified. To identify Bos taxa and investigate their initial management, we conducted zooarchaeological analyses of 193 Bos specimens and sequenced five nuclear genomes from recently excavated assemblages at Bangga. Morphological data indicated that more cattle than yak were present. Ancient mitochondrial DNA and nuclear genome sequences identified taurine cattle and provided evidence for domestic yak and yak-cattle hybridization ~2500 years ago. Reliance on diverse Bos species and their hybrid has increased cattle adaptation and herder resilience to plateau conditions. Ancient cattle and yak at Bangga were closely related to contemporary livestock, indicating early herder legacies and the continuity of cattle and yak husbandry on the Tibetan Plateau.

The macroevolutionary dynamics of mammalian sexual size dimorphism

Sexual size dimorphism (SSD) is a common phenomenon across the animal kingdom. Mammals are unusual in primarily displaying male-biased SSD, where males of a species are typically larger than females. The driving factors behind the evolution of this SSD have been much debated, with popular hypotheses invoking the influence of mating system and social organization via sexual selection, dietary niche divergence and broad-scale correlations with body size (Rensch's rule). Here, we investigate the macroevolutionary origins and maintenance of SSD among mammals, using phylogenetic general mixed linear models and a comprehensive global dataset to evaluate correlations of diet, body mass, seasonality, social organization and mating system with SSD type. We find that SSD as a whole is lost at a greater rate than it is gained, with female-biased SSD being particularly unstable. Non-monogamous mating systems, vertebrate prey consumption and temperature seasonality correlate with male-biased SSD, while polyandry correlates with female-biased SSD, and both types of SSD are positively correlated with body mass. This is in partial contrast to the predictions of Rensch's rule, which predicts that female-biased SSD would correlate negatively with body size. Taken together, our results highlight the importance of considering multiple ecological and social drivers when evaluating the macroevolutionary trajectory of sex differences in body size.

How low can you go? Introducing SeXY: sex identification from low-quantity sequencing data despite lacking assembled sex chromosomes

Accurate sex identification is crucial for elucidating the biology of a species. In the absence of directly observable sexual characteristics, sex identification of wild fauna can be challenging, if not impossible. Molecular sexing offers a powerful alternative to morphological sexing approaches. Here, we present SeXY, a novel sex‐identification pipeline, for very low‐coverage shotgun sequencing data from a single individual. SeXY was designed to utilize low‐effort screening data for sex identification and does not require a conspecific sex‐chromosome assembly as reference. We assess the accuracy of our pipeline to data quantity by downsampling sequencing data from 100,000 to 1000 mapped reads and to reference genome selection by mapping to a variety of reference genomes of various qualities and phylogenetic distance. We show that our method is 100% accurate when mapping to a high‐quality (highly contiguous N50 > 30 Mb) conspecific genome, even down to 1000 mapped reads. For lower‐quality reference assemblies (N50 < 30 Mb), our method is 100% accurate with 50,000 mapped reads, regardless of reference assembly quality or phylogenetic distance. The SeXY pipeline provides several advantages over previously implemented methods; SeXY (i) requires sequencing data from only a single individual, (ii) does not require assembled conspecific sex chromosomes, or even a conspecific reference assembly, (iii) takes into account variation in coverage across the genome, and (iv) is accurate with only 1000 mapped reads in many cases.

Grey wolf genomic history reveals a dual ancestry of dogs

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived^1–8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

Bias assessments to expand research harnessing biological collections

Biological collections are arguably the most important resources for investigations into the impacts of human activities on biodiversity. However, the apparent opportunities presented by museum-derived datasets have not resulted in consistent or widespread use of specimens in ecology outside phenological research and species distribution modeling. We attribute this gap between opportunity and application to biases introduced by collectors, curators, and preservation practices and an imperfect understanding of these biases and how to mitigate them. To facilitate broader use of specimen-based data, we characterize collection biases across key axes and explore interactions among them. We then present a framework for determining the bias assessments needed when extracting data from biological collections. We show that bias assessments required by particular ecological studies will depend on the response variables being measured and the predictor axes of interest. We argue that quantification of biases in specimen-derived datasets is needed to facilitate the widespread application of these data.

The ultimate database to (re)set the evolutionary history of primate genital bones

Scientific literature concerning genital bones in primates consists of both ancient works (dating back to the nineteenth century) and more recent revisions/meta-analyses, which, however, are not always so detailed or exhaustive. Based on a thorough analysis, several conflicting data, inaccurate references, and questionable claims have emerged. We generated a binary matrix of genital bone occurrence data, considering only data at the species level, based on (1) a rigorous literature search protocol, (2) raw data (collected exclusively from primary literature), (3) an updated taxonomy (often tracing back to the species taxonomic history) and (4) new occurrence data from scanned genitals of fresh and museum specimens (using micro-computed tomography-micro-CT). Thanks to this methodological approach, we almost doubled available occurrence data so far, avoiding any arbitrary extension of generic data to conspecific species. This practice, in fact, has been recently responsible for an overestimation of the occurrence data, definitively flattening the interspecific variability. We performed the ancestral state reconstruction analysis of genital bone occurrence and results were mapped onto the most updated phylogeny of primates. As for baculum, we definitively demonstrated its simplesiomorphy for the entire order. As for baubellum, we interpreted all scattered absences as losses, actually proposing (for the first time) a simplesiomorphic state for the clitoral bone as well. The occurrence data obtained, while indirectly confirming the baculum/baubellum homology (i.e., for each baubellum a baculum was invariably present), could also directly demonstrate an intra-specific variability affecting ossa genitalia occurrence. With our results, we established a radically improved and updated database about the occurrence of genital bones in primates, available for further comparative analyses.

Model Systems in Ecology, Evolution, and Behavior: A Call for Diversity in Our Model Systems and Discipline

AbstractEcologists and evolutionary biologists are fascinated by life's variation but also seek to understand phenomena and mechanisms that apply broadly across taxa. Model systems can help us extract generalities from amid all the wondrous diversity, but only if we choose and develop them carefully, use them wisely, and have a range of model systems from which to choose. In this introduction to the Special Feature on Model Systems in Ecology, Evolution, and Behavior (EEB), we begin by grappling with the question, What is a model system? We then explore where our model systems come from, in terms of the skills and other attributes required to develop them and the historical biases that influence traditional model systems in EEB. We emphasize the importance of communities of scientists in the success of model systems-narrow scientific communities can restrict the model organisms themselves. We also consider how our discipline was built around one type of "model scientist"-a history still reflected in the field. This lack of diversity in EEB is unjust and also narrows the field's perspective, including by restricting the questions asked and talents used to answer them. Increasing diversity, equity, and inclusion will require acting at many levels, including structural changes. Diversity in EEB, in both model systems and the scientists who use them, strengthens our discipline.

Relative demographic susceptibility does not explain the extinction chronology of Sahul's megafauna

The causes of Sahul’s megafauna extinctions remain uncertain, although several interacting factors were likely responsible. To examine the relative support for hypotheses regarding plausible ecological mechanisms underlying these extinctions, we constructed the first stochastic, age-structured models for 13 extinct megafauna species from five functional/taxonomic groups, as well as 8 extant species within these groups for comparison. Perturbing specific demographic rates individually, we tested which species were more demographically susceptible to extinction, and then compared these relative sensitivities to the fossil-derived extinction chronology. Our models show that the macropodiformes were the least demographically susceptible to extinction, followed by carnivores, monotremes, vombatiform herbivores, and large birds. Five of the eight extant species were as or more susceptible than the extinct species. There was no clear relationship between extinction susceptibility and the extinction chronology for any perturbation scenario, while body mass and generation length explained much of the variation in relative risk. Our results reveal that the actual mechanisms leading to the observed extinction chronology were unlikely related to variation in demographic susceptibility per se, but were possibly driven instead by finer-scale variation in climate change and/or human prey choice and relative hunting success.

A Multidisciplinary Review of the Inka Imperial Resettlement Policy and Implications for Future Investigations

The rulers of the Inka empire conquered approximately 2 million km2 of the South American Andes in just under 100 years from 1438-1533 CE. Inside the empire, the elite conducted a systematic resettlement of the many Indigenous peoples in the Andes that had been rapidly colonised. The nature of this resettlement phenomenon is recorded within the Spanish colonial ethnohistorical record. Here we have broadly characterised the resettlement policy, despite the often incomplete and conflicting details in the descriptions. We then review research from multiple disciplines that investigate the empirical reality of the Inka resettlement policy, including stable isotope analysis, intentional cranial deformation morphology, ceramic artefact chemical analyses and genetics. Further, we discuss the benefits and limitations of each discipline for investigating the resettlement policy and emphasise their collective value in an interdisciplinary characterisation of the resettlement policy.

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.

Origins and genetic legacy of prehistoric dogs

Dogs were the first domestic animal, but little is known about their population history and to what extent it was linked to humans. We sequenced 27 ancient dog genomes and found that all dogs share a common ancestry distinct from present-day wolves, with limited gene flow from wolves since domestication but substantial dog-to-wolf gene flow. By 11,000 years ago, at least five major ancestry lineages had diversified, demonstrating a deep genetic history of dogs during the Paleolithic. Coanalysis with human genomes reveals aspects of dog population history that mirror humans, including Levant-related ancestry in Africa and early agricultural Europe. Other aspects differ, including the impacts of steppe pastoralist expansions in West and East Eurasia and a near-complete turnover of Neolithic European dog ancestry.

Biomolecular analyses reveal the age, sex and species identity of a near-intact Pleistocene bird carcass

Ancient remains found in permafrost represent a rare opportunity to study past ecosystems. Here, we present an exceptionally well-preserved ancient bird carcass found in the Siberian permafrost, along with a radiocarbon date and a reconstruction of its complete mitochondrial genome. The carcass was radiocarbon dated to approximately 44–49 ka BP, and was genetically identified as a female horned lark. This is a species that usually inhabits open habitat, such as the steppe environment that existed in Siberia at the time. This near-intact carcass highlights the potential of permafrost remains for evolutionary studies that combine both morphology and ancient nucleic acids.

Nicolas Dussex et al. identify a 44,000–49,000 year old bird found in Siberian permafrost as a female horned lark using ancient DNA. This exceptionally well-preserved specimen illustrates the potential contribution to science of permafrost deposits, such as the study of ecology and evolution of ancient ecosystems, calibration of molecular clocks, and furthering our understanding of processes such as biological regulation and gene expression in relation to climate change.

Decadal heat and drought drive body size of North American bison (Bison bison) along the Great Plains

Large grazers are visible and valuable indicators of the effects of projected changes in temperature and drought on grasslands. The grasslands of the Great Plains have supported the greatest number of bison (Bison bison; Linnaeus, 1758) since prehistoric times. We tested the hypothesis that body mass (BM, kg) and asymptotic body mass (ABM, kg) of Bison decline with rising temperature and increasing drought over both temporal and spatial scales along the Great Plains. Temporally, we modeled the relationship of annual measures of BM and height (H, m) of 5,781 Bison at Wind Cave National Park (WICA) from 1966 to 2015. We used Gompertz equations of BM against age to estimate ABM in decadal cohorts; both females and males decreased from the 1960s to the 2010s. Male ABM was variable but consistently larger (699 vs. 441 kg) than female ABM. We used local mean decadal temperature (MDT) and local mean decadal Palmer Drought Severity Index (dPDSI) to model the effects of climate on ABM. Drought decreased ABM temporally (-16 kg/local dPDSI) at WICA. Spatially, we used photogrammetry to measure body height (HE ) of 773 Bison to estimate BME in 19 herds from Saskatchewan to Texas, including WICA. Drought also decreased ABM spatially (-16 kg/local dPDSI) along the Great Plains. Temperature decreased ABM both temporally at WICA (-115 kg/°C local MDT) and spatially (-1 kg/°C local MDT) along the Great Plains. Our data indicate that temperature and drought drive Bison ABM presumably by affecting seasonal mass gain. Bison body size is likely to decline over the next five decades throughout the Great Plains due to projected increases in temperatures and both the frequency and intensity of drought.