Search

Sharing genomic variant interpretations across laboratories promotes consistency in variant assertions. A landscape analysis of Australian clinical genetic-testing laboratories in 2017 identified that, despite the national-accreditation-body recommendations encouraging laboratories to submit genotypic data to clinical databases, fewer than 300 variants had been shared to the ClinVar public database. Consultations with Australian laboratories identified resource constraints limiting routine application of manual processes, consent issues, and differences in interpretation systems as barriers to sharing. This information was used to define key needs and solutions required to enable national sharing of variant interpretations. The Shariant platform, using both the GRCh37 and GRCh38 genome builds, was developed to enable ongoing sharing of variant interpretations and associated evidence between Australian clinical genetic-testing laboratories. Where possible, two-way automated sharing was implemented so that disruption to laboratory workflows would be minimized. Terms of use were developed through consultation and currently restrict access to Australian clinical genetic-testing laboratories. Shariant was designed to store and compare structured evidence, to promote and record resolution of inter-laboratory classification discrepancies, and to streamline the submission of variant assertions to ClinVar. As of December 2021, more than 14,000 largely prospectively curated variant records from 11 participating laboratories have been shared. Discrepant classifications have been identified for 11% (28/260) of variants submitted by more than one laboratory. We have demonstrated that co-design with clinical laboratories is vital to developing and implementing a national variant-interpretation sharing effort. This approach has improved inter-laboratory concordance and enabled opportunities to standardize interpretation practices.

Refereed/Peer-reviewed Hotel quarantine for international travelers has been used to prevent coronavirus disease spread into Australia. A quarantine hotel–associated community outbreak was detected in South Australia. Real-time genomic sequencing enabled rapid confirmation tracking the outbreak to a recently returned traveler and linked 2 cases of infection in travelers at the same facility.

The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was the largest carnivorous Australian marsupial to survive into the modern era. Despite last sharing a common ancestor with the eutherian canids ~160 million years ago, their phenotypic resemblance is considered the most striking example of convergent evolution in mammals. The last known thylacine died in captivity in 1936 and many aspects of the evolutionary history of this unique marsupial apex predator remain unknown. Here we have sequenced the genome of a preserved thylacine pouch young specimen to clarify the phylogenetic position of the thylacine within the carnivorous marsupials, reconstruct its historical demography and examine the genetic basis of its convergence with canids. Retroposon insertion patterns placed the thylacine as the basal lineage in Dasyuromorphia and suggest incomplete lineage sorting in early dasyuromorphs. Demographic analysis indicated a long-term decline in genetic diversity starting well before the arrival of humans in Australia. In spite of their extraordinary phenotypic convergence, comparative genomic analyses demonstrated that amino acid homoplasies between the thylacine and canids are largely consistent with neutral evolution. Furthermore, the genes and pathways targeted by positive selection differ markedly between these species. Together, these findings support models of adaptive convergence driven primarily by cis-regulatory evolution.

The Pairwise Sequentially Markov Coalescent (PSMC), and its extension PSMC′, model past population sizes from a single diploid genome. Both models have been widely applied, even to organisms with scaffold-level genome reference assemblies of limited contiguity. However it is unclear how PSMC and PSMC′ perform on short scaffolds. We evaluated psmc and msmc, implementations of the PSMC and PSMC′ models respectively, on simulated genomes with low contiguity, and compared results to those from fully contiguous data. Simulations with scaffolds from 100 Mb to 10 kb revealed that psmc maintains high accuracy down to lengths of 100 kb, while msmc is accurate down to 1 Mb. The discrepancy is not due to differing models, but stems from an implementation detail of msmc—homozygous tracts at the ends of scaffolds are discarded, making msmc unreliable for low contiguity genomes. We recommend excluding data that are aligned to shorter scaffolds when undertaking demographic inference.

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. Analysis of Aboriginal Australian mitochondrial genomes shows geographic patterns and deep splits across the major haplogroups that indicate a single, rapid migration along the coasts around 49–45 ka, followed by longstanding persistence in discrete geographic areas. Aboriginal Australians preserve one of the longest continuous cultural complexes known, with archaeological evidence dating initial settlement of the continent to around 50,000 years ago. Alan Cooper and colleagues have charted the subsequent progress of humanity in and around the continent in the form of 111 mitochondrial genomes from preserved hair samples. The results show that, from landfall in the north of Australia, people spread rapidly around the east and west coasts, meeting in southern Australia as early as 49,000 years ago. Strong regional patterns of mitochondrial DNA variation suggest that when people stopped moving they stayed put, putting down cultural roots that have weathered 50,000 years of significant cultural and climatic change.

Weyrich, Laura S. et al., Recent genomic data have revealed multiple interactions between Neanderthals and modern humans1, 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 gathering2, 3. 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 abscess4 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 Australian Research Council supported this work through the Discovery Project and Fellowship schemes.

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 silico75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided90% 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.

SummaryThe near-simultaneous appearance of most modern animal body plans (phyla) ∼530 million years ago during the Cambrian explosion is strong evidence for a brief interval of rapid phenotypic and genetic innovation, yet the exact speed and nature of this grand adaptive radiation remain debated [1–12]. Crucially, rates of morphological evolution in the past (i.e., in ancestral lineages) can be inferred from phenotypic differences among living organisms—just as molecular evolutionary rates in ancestral lineages can be inferred from genetic divergences [13]. We here employed Bayesian [14] and maximum likelihood [15] phylogenetic clock methods on an extensive anatomical [16] and genomic [17] data set for arthropods, the most diverse phylum in the Cambrian and today. Assuming an Ediacaran origin for arthropods, phenotypic evolution was ∼4 times faster, and molecular evolution ∼5.5 times faster, during the Cambrian explosion compared to all subsequent parts of the Phanerozoic. These rapid evolutionary rates are robust to assumptions about the precise age of arthropods. Surprisingly, these fast early rates do not change substantially even if the radiation of arthropods is compressed entirely into the Cambrian (∼542 mega-annum [Ma]) or telescoped into the Cryogenian (∼650 Ma). The fastest inferred rates are still consistent with evolution by natural selection and with data from living organisms, potentially resolving “Darwin’s dilemma.” However, evolution during the Cambrian explosion was unusual (compared to the subsequent Phanerozoic) in that fast rates were present across many lineages.

High-throughput sequencing (HTS) of PCR amplicons is becoming the method of choice to sequence one or several targeted loci for phylogenetic and DNA barcoding studies. Although the development of HTS has allowed rapid generation of massive amounts of DNA sequence data, preparing amplicons for HTS remains a rate-limiting step. For example, HTS platforms require platform-specific adapter sequences to be present at the 5' and 3' end of the DNA fragment to be sequenced. In addition, short multiplex identifier (MID) tags are typically added to allow multiple samples to be pooled in a single HTS run. Existing methods to incorporate HTS adapters and MID tags into PCR amplicons are either inefficient, requiring multiple enzymatic reactions and clean-up steps, or costly when applied to multiple samples or loci (fusion primers). We describe a method to amplify a target locus and add HTS adapters and MID tags via a linker sequence using a single PCR. We demonstrate our approach by generating reference sequence data for two mitochondrial loci (COI and 16S) for a diverse suite of insect taxa. Our approach provides a flexible, cost-effective and efficient method to prepare amplicons for HTS.

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 (, This research was supported by the Australian Research Council, the European Commission (PIRSES-GA-2009-247652—BIOGEAST), the Polish National Science Centre (N N304 301940 and 2013/11/B/NZ8/00914), the Danish National Research Foundation (DNRF94), the Marie Curie International Outgoing Fellowship (7th European Community Framework Program—MEDITADNA, PIOF-GA-2011-300854, FP7-PEOPLE) and the Russian Foundation for Basic Research (N 15-04-03882).

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.

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.

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.

Author version made available in accordance with Publisher copyright policy., The evolution of the ratite birds has been widely attributed to vicariant speciation, driven by the Cretaceous breakup of the supercontinent Gondwana. The early isolation of Africa and Madagascar implies that the ostrich and extinct Madagascan elephant birds (Aepyornithidae) should be the oldest ratite lineages. We sequenced the mitochondrial genomes of two elephant birds and performed phylogenetic analyses, which revealed that these birds are the closest relatives of the New Zealand kiwi and are distant from the basal ratite lineage of ostriches. This unexpected result strongly contradicts continental vicariance and instead supports flighted dispersal in all major ratite lineages. We suggest that convergence toward gigantism and flightlessness was facilitated by early Tertiary expansion into the diurnal herbivory niche after the extinction of the dinosaurs.

Entre els membres del The Genographic Consortium es troben els autors: Jaume Bertranpetit, David Comas, Marc Haber i Begoña Martínez Cruz The human mitochondrial haplogroup C1 has a broad global distribution but is extremely rare in Europe today. Recent ancient DNA evidence has demonstrated its presence in European Mesolithic individuals. Three individuals from the 7,500 year old Mesolithic site of Yuzhnyy Oleni Ostrov, Western Russia, could be assigned to haplogroup C1 based on mitochondrial hypervariable region I sequences. However, hypervariable region I data alone could not provide enough resolution to establish the phylogenetic relationship of these Mesolithic haplotypes with haplogroup C1 mitochondrial DNA sequences found today in populations of Europe, Asia and the Americas. In order to obtain high-resolution data and shed light on the origin of this European Mesolithic C1 haplotype, we target-enriched and sequenced the complete mitochondrial genome of one Yuzhnyy Oleni Ostrov C1 individual. The updated phylogeny of C1 haplogroups indicated that the Yuzhnyy Oleni Ostrov haplotype represents a new distinct clade, provisionally coined “C1f”. We show that all three C1 carriers of Yuzhnyy Oleni Ostrov belong to this clade. No haplotype closely related to the C1f sequence could be found in the large current database of ancient and present-day mitochondrial genomes. Hence, we have discovered past human mitochondrial diversity that has not been observed in modern-day populations so far. The lack of positive matches in modern populations may be explained by under-sampling of rare modern C1 carriers or by demographic processes, population extinction or replacement, that may have impacted on populations of Northeast Europe since prehistoric times. This research was supported by The Genographic Project, which is supported by funding from the National Geographic Society, IBM, and the Waitt Family Foundation. PB was funded by the Australian Research Council Linkage Project LP0882622 awarded to AC, and the Discovery Projects DP1095782 and DP130102158 awarded to WH and BL. OB was funded by the RAS Programmes ‘‘Molecular and cell biology’’, ‘‘Gene pool dynamics’’ and RFBR grants 13-04-01711, Cover Letter 13-06-0067

BackgroundMitochondrial DNA (mtDNA) typing can be a useful aid for identifying people from compromised samples when nuclear DNA is too damaged, degraded or below detection thresholds for routine short tandem repeat (STR)-based analysis. Standard mtDNA typing, focused on PCR amplicon sequencing of the control region (HVS I and HVS II), is limited by the resolving power of this short sequence, which misses up to 70% of the variation present in the mtDNA genome.MethodsWe used in-solution hybridisation-based DNA capture (using DNA capture probes prepared from modern human mtDNA) to recover mtDNA from post-mortem human remains in which the majority of DNA is both highly fragmented (ResultsWe sequenced whole mitochondrial genomes for human identification from samples where standard nuclear STR typing produced only partial profiles or demonstrably failed and/or where standard mtDNA hypervariable region sequences lacked resolving power. Multiple rounds of enrichment can substantially improve coverage and sequencing depth of mtDNA genomes from highly degraded samples. The application of this method has led to the reliable mitochondrial sequencing of human skeletal remains from unidentified World War Two (WWII) casualties approximately 70 years old and from archaeological remains (up to 2,500 years old).ConclusionsThis approach has potential applications in forensic science, historical human identification cases, archived medical samples, kinship analysis and population studies. In particular the methodology can be applied to any case, involving human or non-human species, where whole mitochondrial genome sequences are required to provide the highest level of maternal lineage discrimination. Multiple rounds of in-solution hybridisation-based DNA capture can retrieve whole mitochondrial genome sequences from even the most challenging samples.

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.

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., Facultad de Ciencias Naturales y Museo

Brotherton, Paul et al.-- The Genographic Consortium, Haplogroup H dominates present-day Western European mitochondrial DNA variability (>40%), yet was less common (∼19%) among Early Neolithic farmers (∼5450 BC) and virtually absent in Mesolithic hunter-gatherers. Here we investigate this major component of the maternal population history of modern Europeans and sequence 39 complete haplogroup H mitochondrial genomes from ancient human remains. We then compare this 'real-time' genetic data with cultural changes taking place between the Early Neolithic (∼5450 BC) and Bronze Age (∼2200 BC) in Central Europe. Our results reveal that the current diversity and distribution of haplogroup H were largely established by the Mid Neolithic (∼4000 BC), but with substantial genetic contributions from subsequent pan-European cultures such as the Bell Beakers expanding out of Iberia in the Late Neolithic (∼2800 BC). Dated haplogroup H genomes allow us to reconstruct the recent evolutionary history of haplogroup H and reveal a mutation rate 45% higher than current estimates for human mitochondria. © 2013 Macmillan Publishers Limited. All rights reserved., We thank the Australian Research Council (grant LP0882622), the Deutsche Forschungsgemeinschaft (Al 287/7-1 and Me 3245/1-1) and National Geographic’s Genographic Project for funding. M.v.O. was supported in part by the Netherlands Forensic Institute (NFI) and a grant from the Netherlands Genomics Initiative (NGI)/Netherlands Organization for Scientific Research (NWO) within the framework of the Forensic Genomics Consortium Netherlands (FGCN).

International audience; 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.

The Black Rat (Rattus rattus) spread out of Asia to become one of the world’s 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. Citation: Aplin KP, Suzuki H, Chinen AA, Chesser RT, ten Have J, et al. (2011) Multiple Geographic Origins of Commensalism and Complex Dispersal History of

The origin of the European bison (Bison bonasus, Linnaeus, 1758) has been widely discussed and investigated in recent years. The species had a wide historic geographic distribution throughout the European continent during the middle and late Holocene, ranging from France in the west to the Caucasus in the east. However, archaeological evidence is needed to resolve the southern extent of the European bison distribution. We discovered one bison skull fragment during archaeological excavations in 2008 in the area of Yenikapi Metro and Marmaray (Turkey). Radiocarbon dating indicated the skull was deposited during the Byzantine period (seventh to eighth century AD). Mitochondrial genome analyses provided clear evidence that the skull was from a European bison. This is the first unambiguous evidence of the presence of this species in southeastern Europe during Byzantine times, which validates the historical written records of a potentially wider range of the European bison in historical times.