Julien Soubrier

Multiple geographic origins of commensalism and complex dispersal history of black rats

The Black Rat (Rattus rattus) spread out of Asia to become one of the worlds worst agricultural and urban pests, and a reservoir or vector of numerous zoonotic diseases, including the devastating plague. Despite the global scale and inestimable cost of their impacts on both human livelihoods and natural ecosystems, little is known of the global genetic diversity of Black Rats, the timing and directions of their historical dispersals, and the risks associated with contemporary movements. We surveyed mitochondrial DNA of Black Rats collected across their global range as a first step towards obtaining an historical genetic perspective on this socioeconomically important group of rodents. We found a strong phylogeographic pattern with well-differentiated lineages of Black Rats native to South Asia, the Himalayan region, southern Indochina, and northern Indochina to East Asia, and a diversification that probably commenced in the early Middle Pleistocene. We also identified two other currently recognised species of Rattus as potential derivatives of a paraphyletic R. rattus. Three of the four phylogenetic lineage units within R. rattus show clear genetic signatures of major population expansion in prehistoric times, and the distribution of particular haplogroups mirrors archaeologically and historically documented patterns of human dispersal and trade. Commensalism clearly arose multiple times in R. rattus and in widely separated geographic regions, and this may account for apparent regionalism in their associated pathogens. Our findings represent an important step towards deeper understanding the complex and influential relationship that has developed between Black Rats and humans, and invite a thorough re-examination of host-pathogen associations among Black Rats.

Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias

Studies of insect assemblages are suited to the simultaneous DNA‐based identification of multiple taxa known as metabarcoding. To obtain accurate estimates of diversity, metabarcoding markers ideally possess appropriate taxonomic coverage to avoid PCR‐amplification bias, as well as sufficient sequence divergence to resolve species. We used in silico PCR to compare the taxonomic coverage and resolution of newly designed insect metabarcodes (targeting 16S) with that of existing markers [16S and cytochrome oxidase c subunit I (COI)] and then compared their efficiency in vitro. Existing metabarcoding primers amplified in silico <75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided >90% coverage. Furthermore, metabarcodes targeting COI appeared to introduce taxonomic PCR‐amplification bias, typically amplifying a greater percentage of Lepidoptera and Diptera species, while failing to amplify certain orders in silico. To test whether bias predicted in silico was observed in vitro, we created an artificial DNA blend containing equal amounts of DNA from 14 species, representing 11 insect orders and one arachnid. We PCR‐amplified the blend using five primer sets, targeting either COI or 16S, with high‐throughput amplicon sequencing yielding more than 6 million reads. In vitro results typically corresponded to in silico PCR predictions, with newly designed 16S primers detecting 11 insect taxa present, thus providing equivalent or better taxonomic coverage than COI metabarcodes. Our results demonstrate that in silico PCR is a useful tool for predicting taxonomic bias in mixed template PCR and that researchers should be wary of potential bias when selecting metabarcoding markers.

Ancient mitochondrial DNA provides high-resolution time scale of the peopling of the Americas.

Native American population history is reexamined using a large data set of pre-Columbian mitochondrial genomes. The exact timing, route, and process of the initial peopling of the Americas remains uncertain despite much research. Archaeological evidence indicates the presence of humans as far as southern Chile by 14.6 thousand years ago (ka), shortly after the Pleistocene ice sheets blocking access from eastern Beringia began to retreat. Genetic estimates of the timing and route of entry have been constrained by the lack of suitable calibration points and low genetic diversity of Native Americans. We sequenced 92 whole mitochondrial genomes from pre-Columbian South American skeletons dating from 8.6 to 0.5 ka, allowing a detailed, temporally calibrated reconstruction of the peopling of the Americas in a Bayesian coalescent analysis. The data suggest that a small population entered the Americas via a coastal route around 16.0 ka, following previous isolation in eastern Beringia for ~2.4 to 9 thousand years after separation from eastern Siberian populations. Following a rapid movement throughout the Americas, limited gene flow in South America resulted in a marked phylogeographic structure of populations, which persisted through time. All of the ancient mitochondrial lineages detected in this study were absent from modern data sets, suggesting a high extinction rate. To investigate this further, we applied a novel principal components multiple logistic regression test to Bayesian serial coalescent simulations. The analysis supported a scenario in which European colonization caused a substantial loss of pre-Columbian lineages.

Neanderthal behaviour, diet, and disease inferred from ancient DNA in dental calculus

Recent genomic data have revealed multiple interactions between Neanderthals and modern humans, but there is currently little genetic evidence regarding Neanderthal behaviour, diet, or disease. Here we describe the shotgun-sequencing of ancient DNA from five specimens of Neanderthal calcified dental plaque (calculus) and the characterization of regional differences in Neanderthal ecology. At Spy cave, Belgium, Neanderthal diet was heavily meat based and included woolly rhinoceros and wild sheep (mouflon), characteristic of a steppe environment. In contrast, no meat was detected in the diet of Neanderthals from El Sidrón cave, Spain, and dietary components of mushrooms, pine nuts, and moss reflected forest gathering. Differences in diet were also linked to an overall shift in the oral bacterial community (microbiota) and suggested that meat consumption contributed to substantial variation within Neanderthal microbiota. Evidence for self-medication was detected in an El Sidrón Neanderthal with a dental abscess and a chronic gastrointestinal pathogen (Enterocytozoon bieneusi). Metagenomic data from this individual also contained a nearly complete genome of the archaeal commensal Methanobrevibacter oralis (10.2× depth of coverage)—the oldest draft microbial genome generated to date, at around 48,000 years old. DNA preserved within dental calculus represents a notable source of information about the behaviour and health of ancient hominin specimens, as well as a unique system that is useful for the study of long-term microbial evolution.

The influence of rate heterogeneity among sites on the time dependence of molecular rates

Molecular evolutionary rate estimates have been shown to depend on the time period over which they are estimated. Factors such as demographic processes, calibration errors, purifying selection, and the heterogeneity of substitution rates among sites (RHAS) are known to affect the accuracy with which rates of evolution are estimated. We use mathematical modeling and Bayesian analyses of simulated sequence alignments to explore how mutational hotspots can lead to time-dependent rate estimates. Mathematical modeling shows that underestimation of molecular rates over increasing time scales is inevitable when RHAS is ignored. Although a gamma distribution is commonly used to model RHAS, we show that when the actual RHAS deviates from a gamma-like distribution, rates can either be under- or overestimated in a time-dependent manner. Simulations performed under different scenarios of RHAS confirm the mathematical modeling and demonstrate the impacts of time-dependent rates on estimates of divergence times. Most notably, erroneous rate estimates can have narrow credibility intervals, leading to false confidence in biased estimates of rates, and node ages. Surprisingly, large errors in estimates of overall molecular rate do not necessarily generate large errors in divergence time estimates. Finally, we illustrate the correlation between time-dependent rate patterns and differential saturation between quickly and slowly evolving sites. Our results suggest that data partitioning or simple nonparametric mixture models of RHAS significantly improve the accuracy with which node ages and substitution rates can be estimated.

Aboriginal mitogenomes reveal 50,000 years of regionalism in Australia

Aboriginal Australians represent one of the longest continuous cultural complexes known. Archaeological evidence indicates that Australia and New Guinea were initially settled approximately 50 thousand years ago (ka); however, little is known about the processes underlying the enormous linguistic and phenotypic diversity within Australia. Here we report 111 mitochondrial genomes (mitogenomes) from historical Aboriginal Australian hair samples, whose origins enable us to reconstruct Australian phylogeographic history before European settlement. Marked geographic patterns and deep splits across the major mitochondrial haplogroups imply that the settlement of Australia comprised a single, rapid migration along the east and west coasts that reached southern Australia by 49–45 ka. After continent-wide colonization, strong regional patterns developed and these have survived despite substantial climatic and cultural change during the late Pleistocene and Holocene epochs. Remarkably, we find evidence for the continuous presence of populations in discrete geographic areas dating back to around 50 ka, in agreement with the notable Aboriginal Australian cultural attachment to their country.

The origins of the enigmatic Falkland Islands wolf

The origins of the extinct Falkland Islands wolf (FIW), Dusicyon australis, have remained a mystery since it was first recorded by Europeans in the seventeenth century. It is the only terrestrial mammal on the Falkland Islands (also known as the Malvinas Islands), which lie ~460 km from Argentina, leading to suggestions of either human-mediated transport or overwater dispersal. Previous studies used ancient DNA from museum specimens to suggest that the FIW diverged from its closest living relative, the South American maned wolf (Chrysocyon brachyurus) around 7 Ma, and colonized the islands ~330 ka by unknown means. Here we retrieve ancient DNA from subfossils of an extinct mainland relative, Dusicyon avus, and reveal the FIW lineage became isolated only 16 ka (8-31 ka), during the last glacial phase. Submarine terraces, formed on the Argentine coastal shelf by low sea-stands during this period, suggest that the FIW colonized via a narrow, shallow marine strait, potentially while it was frozen over.

Early cave art and ancient DNA record the origin of European bison

The two living species of bison (European and American) are among the few terrestrial megafauna to have survived the late Pleistocene extinctions. Despite the extensive bovid fossil record in Eurasia, the evolutionary history of the European bison (or wisent, Bison bonasus) before the Holocene (<11.7 thousand years ago (kya)) remains a mystery. We use complete ancient mitochondrial genomes and genome-wide nuclear DNA surveys to reveal that the wisent is the product of hybridization between the extinct steppe bison (Bison priscus) and ancestors of modern cattle (aurochs, Bos primigenius) before 120 kya, and contains up to 10% aurochs genomic ancestry. Although undetected within the fossil record, ancestors of the wisent have alternated ecological dominance with steppe bison in association with major environmental shifts since at least 55 kya. Early cave artists recorded distinct morphological forms consistent with these replacement events, around the Last Glacial Maximum (LGM, ∼21–18 kya).

Experimental conditions improving in solution target enrichment for ancient DNA

High‐throughput sequencing has dramatically fostered ancient DNA research in recent years. Shotgun sequencing, however, does not necessarily appear as the best‐suited approach due to the extensive contamination of samples with exogenous environmental microbial DNA. DNA capture‐enrichment methods represent cost‐effective alternatives that increase the sequencing focus on the endogenous fraction, whether it is from mitochondrial or nuclear genomes, or parts thereof. Here, we explored experimental parameters that could impact the efficacy of MYbaits in‐solution capture assays of ~5000 nuclear loci or the whole genome. We found that varying quantities of the starting probes had only moderate effect on capture outcomes. Starting DNA, probe tiling, the hybridization temperature and the proportion of endogenous DNA all affected the assay, however. Additionally, probe features such as their GC content, number of CpG dinucleotides, sequence complexity and entropy and self‐annealing properties need to be carefully addressed during the design stage of the capture assay. The experimental conditions and probe molecular features identified in this study will improve the recovery of genetic information extracted from degraded and ancient remains.

Pleistocene Chinese cave hyenas and the recent Eurasian history of the spotted hyena, Crocuta crocuta

The living hyena species (spotted, brown, striped and aardwolf) are remnants of a formerly diverse group of more than 80 fossil species, which peaked in diversity in the Late Miocene (about 7–8 Ma). The fossil history indicates an African origin, and morphological and ancient DNA data have confirmed that living spotted hyenas (Crocuta crocuta) of Africa were closely related to extinct Late Pleistocene cave hyenas from Europe and Asia. The current model used to explain the origins of Eurasian cave hyena populations invokes multiple migrations out of Africa between 3.5–0.35 Ma. We used mitochondrial DNA sequences from radiocarbon‐dated Chinese Pleistocene hyena specimens to examine the origin of Asian populations, and temporally calibrate the evolutionary history of spotted hyenas. Our results support a far more recent evolutionary timescale (430–163 kya) and suggest that extinct and living spotted hyena populations originated from a widespread Eurasian population in the Late Pleistocene, which was only subsequently restricted to Africa. We developed statistical tests of the contrasting population models and their fit to the fossil record. Coalescent simulations and Bayes Factor analysis support the new radiocarbon‐calibrated timescale and Eurasian origins model. The new Eurasian biogeographic scenario proposed for the hyena emphasizes the role of the vast steppe grasslands of Eurasia in contrast to models only involving Africa. The new methodology for combining genetic and geological data to test contrasting models of population history will be useful for a wide range of taxa where ancient and historic genetic data are available.