Search

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an similar to 120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for similar to 500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

Near the end of the Pleistocene epoch, populations of the woolly mammoth (Mammuthus primigenius) were distributed across parts of three continents, from western Europe and northern Asia through Beringia to the Atlantic seaboard of North America. Nonetheless, questions about the connectivity and temporal continuity of mammoth populations and species remain unanswered. We use a combination of targeted enrichment and high-throughput sequencing to assemble and interpret a data set of 143 mammoth mitochondrial genomes, sampled from fossils recovered from across their Holarctic range. Our dataset includes 54 previously unpublished mitochondrial genomes and significantly increases the coverage of the Eurasian range of the species. The resulting global phylogeny confirms that the Late Pleistocene mammoth population comprised three distinct mitochondrial lineages that began to diverge ~1.0-2.0 million years ago (Ma). We also find that mammoth mitochondrial lineages were strongly geographically partitioned throughout the Pleistocene. In combination, our genetic results and the pattern of morphological variation in time and space suggest that male-mediated gene flow, rather than large-scale dispersals, was important in the Pleistocene evolutionary history of mammoths.

Near the end of the Pleistocene epoch, populations of the woolly mammoth (Mammuthus primigenius) were distributed across parts of three continents, from western Europe and northern Asia through Beringia to the Atlantic seaboard of North America. Nonetheless, questions about the connectivity and temporal continuity of mammoth populations and species remain unanswered. We use a combination of targeted enrichment and high-throughput sequencing to assemble and interpret a data set of 143 mammoth mitochondrial genomes, sampled from fossils recovered from across their Holarctic range. Our dataset includes 54 previously unpublished mitochondrial genomes and significantly increases the coverage of the Eurasian range of the species. The resulting global phylogeny confirms that the Late Pleistocene mammoth population comprised three distinct mitochondrial lineages that began to diverge ~1.0-2.0 million years ago (Ma). We also find that mammoth mitochondrial lineages were strongly geographically partitioned throughout the Pleistocene. In combination, our genetic results and the pattern of morphological variation in time and space suggest that male-mediated gene flow, rather than large-scale dispersals, was important in the Pleistocene evolutionary history of mammoths.

BackgroundPigeons and doves (Columbiformes) are one of the oldest and most diverse extant lineages of birds. However, the nature and timing of the group's evolutionary radiation remains poorly resolved, despite recent advances in DNA sequencing and assembly and the growing database of pigeon mitochondrial genomes. One challenge has been to generate comparative data from the large number of extinct pigeon lineages, some of which are morphologically unique and therefore difficult to place in a phylogenetic context.ResultsWe used ancient DNA and next generation sequencing approaches to assemble complete mitochondrial genomes for eleven pigeons, including the extinct Ryukyu wood pigeon (Columba jouyi), the thick-billed ground dove (Alopecoenas salamonis), the spotted green pigeon (Caloenas maculata), the Rodrigues solitaire (Pezophaps solitaria), and the dodo (Raphus cucullatus). We used a Bayesian approach to infer the evolutionary relationships among 24 species of living and extinct pigeons and doves.ConclusionsOur analyses indicate that the earliest radiation of the Columbidae crown group most likely occurred during the Oligocene, with continued divergence of major clades into the Miocene, suggesting that diversification within the Columbidae occurred more recently than has been reported previously.

Smallpox holds a unique position in the history of medicine. It was the first disease for which a vaccine was developed and remains the only human disease eradicated by vaccination. Although there have been claims of smallpox in Egypt, India, and China dating back millennia [1-4], the timescale of emergence of the causative agent, variola virus (VARV), and how it evolved in the context of increasingly widespread immunization, have proven controversial [4-9]. In particular, some molecular-clock-based studies have suggested that key events in VARV evolution only occurred during the last two centuries [4-6] and hence in apparent conflict with anecdotal historical reports, although it is difficult to distinguish smallpox from other pustular rashes by description alone. To address these issues, we captured, sequenced, and reconstructed a draft genome of an ancient strain of VARV, sampled from a Lithuanian child mummy dating between 1643 and 1665 and close to the time of several documented European epidemics [1, 2, 10]. When compared to vaccinia virus, this archival strain contained the same pattern of gene degradation as 20th century VARVs, indicating that such loss of gene function had occurred before ca. 1650. Strikingly, the mummy sequence fell basal to all currently sequenced strains of VARV on phylogenetic trees. Molecular-clock analyses revealed a strong clock-like structure and that the timescale of smallpox evolution is more recent than often supposed, with the diversification of major viral lineages only occurring within the 18th and 19th centuries, concomitant with the development of modern vaccination.

BackgroundPigeons and doves (Columbiformes) are one of the oldest and most diverse extant lineages of birds. However, the nature and timing of the group's evolutionary radiation remains poorly resolved, despite recent advances in DNA sequencing and assembly and the growing database of pigeon mitochondrial genomes. One challenge has been to generate comparative data from the large number of extinct pigeon lineages, some of which are morphologically unique and therefore difficult to place in a phylogenetic context.ResultsWe used ancient DNA and next generation sequencing approaches to assemble complete mitochondrial genomes for eleven pigeons, including the extinct Ryukyu wood pigeon (Columba jouyi), the thick-billed ground dove (Alopecoenas salamonis), the spotted green pigeon (Caloenas maculata), the Rodrigues solitaire (Pezophaps solitaria), and the dodo (Raphus cucullatus). We used a Bayesian approach to infer the evolutionary relationships among 24 species of living and extinct pigeons and doves.ConclusionsOur analyses indicate that the earliest radiation of the Columbidae crown group most likely occurred during the Oligocene, with continued divergence of major clades into the Miocene, suggesting that diversification within the Columbidae occurred more recently than has been reported previously.

The processes leading up to species extinctions are typically characterized by prolonged declines in population size and geographic distribution, followed by a phase in which populations are very small and may be subject to intrinsic threats, including loss of genetic diversity and inbreeding [1]. However, whether such genetic factors have had an impact on species prior to their extinction is unclear [2, 3]; examining this would require a detailed reconstruction of a species' demographic history as well as changes in genome-wide diversity leading up to its extinction. Here, we present high-quality complete genome sequences from two woolly mammoths (Mammuthus primigenius). The first mammoth was sequenced at 17.1-fold coverage and dates to similar to 4,300 years before present, representing one of the last surviving individuals on Wrangel Island. The second mammoth, sequenced at 11.2-fold coverage, was obtained from an similar to 44,800-year-old specimen from the Late Pleistocene population in northeastern Siberia. The demographic trajectories inferred from the two genomes are qualitatively similar and reveal a population bottleneck during the Middle or Early Pleistocene, and a more recent severe decline in the ancestors of the Wrangel mammoth at the end of the last glaciation. A comparison of the two genomes shows that the Wrangel mammoth has a 20% reduction in heterozygosity as well as a 28-fold increase in the fraction of the genome that comprises runs of homozygosity. We conclude that the population on Wrangel Island, which was the last surviving woolly mammoth population, was subject to reduced genetic diversity shortly before it became extinct.

Mammoths (Mammuthus) have been studied extensively at the genetic level. However due to both taphonomic and technological limitations, only one of several late Pleistocene mammoth species, the woolly mammoth (M. primigenius), has been investigated. This limits our impression of mammoth population history to the the northern latitudes, just one of several environments in which mammoths lived and went extinct. It also obscures their evolutionary chronology, which prevents proper climatic and biogeographic contextualization of their history. Fortunately recent technological advances in high-throughput sequencing and targeted enrichment promise to expand Pleistocene faunal population phylogeography to non-permafrost, non-cave burial contexts. However the capacity and behavior of these combined technologies for characterizing ancient DNA is largely unexplored, preventing efficient and routine use for population-level studies. In this thesis I test and apply these technologies to remains of mammoth species from throughout North America. I first demonstrate their potential for poorly-preserved DNA, and then I evaluate their efficient application to large sample sets, as well as for capturing complete nuclear genomes. I then use these technologies to sequence dozens of mitochondrial genomes from Columbian (M. columbi) and other non-woolly mammoths, reconstructing their matrilineal phylogeography south of the ice. The revealed patterns not only imply a deep chronology for mammoth matrilineal diversity, but also that North American mammoth evolution was occurred via separate episodes of interbreeding between resident and invading populations, and between ecotypes. Overall the biological and methodological discoveries afforded by this body of work outline future research avenues on mammoth evolution, behavior, and extinction., Doctor of Philosophy (PhD)

Late Pleistocene North America hosted at least two divergent and ecologically distinct species of mammoth: the periglacial woolly mammoth (Mammuthus primigenius) and the subglacial Columbian mammoth (Mammuthus columbi). To date, mammoth genetic research has been entirely restricted to woolly mammoths, rendering their genetic evolution difficult to contextualize within broader Pleistocene paleoecology and biogeography. Here, we take an interspecific approach to clarifying mammoth phylogeny by targeting Columbian mammoth remains for mitogenomic sequencing. We sequenced the first complete mitochondrial genome of a classic Columbian mammoth, as well as the first complete mitochondrial genome of a North American woolly mammoth. Somewhat contrary to conventional paleontological models, which posit that the two species were highly divergent, the M. columbi mitogenome we obtained falls securely within a subclade of endemic North American M. primigenius. Though limited, our data suggest that the two species interbred at some point in their evolutionary histories. One potential explanation is that woolly mammoth haplotypes entered Columbian mammoth populations via introgression at subglacial ecotones, a scenario with compelling parallels in extant elephants and consistent with certain regional paleontological observations. This highlights the need for multi-genomic data to sufficiently characterize mammoth evolutionary history. Our results demonstrate that the use of next-generation sequencing technologies holds promise in obtaining such data, even from non-cave, non-permafrost Pleistocene depositional contexts.

We have used a systematic protocol for extracting, quantitating, sexing and validating ancient human mitochondrial and nuclear DNA of one male and one female Beothuk, a Native American population from Newfoundland, which became extinct approximately 180 years ago. They carried mtDNA haplotypes, which fall within haplogroups X and C, consistent with Northeastern Native populations today. In addition we have sexed the male using a novel-sexing assay and confirmed the authenticity of his Y chromosome with the presence of the Native American specific Y-QM3 single nucleotide polymorphism (SNP). This is the first ancient nuclear SNP typed from a Native population in the Americas. In addition, using the same teeth we conducted a stable isotopes analysis of collagen and dentine to show that both individuals relied on marine sources (fresh and salt water fish, seals) with no hierarchy seen between them, and that their water sources were pooled or stored water. Both mtDNA sequence data and Y SNP data hint at possible gene flow or a common ancestral population for both the Beothuk and the current day Mikmaq, but more importantly the data do not lend credence to the proposed idea that the Beothuk (specifically, Nonosabasut) were of admixed (European-Native American) descent. We also analyzed patterns of DNA damage in the clones of authentic mtDNA sequences; there is no tendency for DNA damage to occur preferentially at previously defined mutational hotspots, suggesting that such mutational hotspots are not hypervariable because they are more prone to damage. Udgivelsesdato: 2007-Apr

Udgivelsesdato: January 2007, Although ancient DNA (aDNA) miscoding lesions have been studied since the earliest days of the field, their nature remains a source of debate. A variety of conflicting hypotheses exist about which miscoding lesions constitute true aDNA damage as opposed to PCR polymerase amplification error. Furthermore, considerable disagreement and speculation exists on which specific damage events underlie observed miscoding lesions. The root of the problem is that it has previously been difficult to assemble sufficient data to test the hypotheses, and near-impossible to accurately determine the specific strand of origin of observed damage events. With the advent of emulsion-based clonal amplification (emPCR) and the sequencing-by-synthesis technology this has changed. In this paper we demonstrate how data produced on the Roche GS20 genome sequencer can determine miscoding lesion strands of origin, and subsequently be interpreted to enable characterization of the aDNA damage behind the observed phenotypes. Through comparative analyses on 390,965 bp of modern chloroplast and 131,474 bp of ancient woolly mammoth GS20 sequence data we conclusively demonstrate that in this sample at least, a permafrost preserved specimen, Type 2 (cytosine-->thymine/guanine-->adenine) miscoding lesions represent the overwhelming majority of damage-derived miscoding lesions. Additionally, we show that an as yet unidentified guanine-->adenine analogue modification, not the conventionally argued cytosine-->uracil deamination, underpins a significant proportion of Type 2 damage. How widespread these implications are for aDNA will become apparent as future studies analyse data recovered from a wider range of substrates.

Udgivelsesdato: 2007-Jul-6, It is difficult to obtain fossil data from the 10% of Earth's terrestrial surface that is covered by thick glaciers and ice sheets, and hence, knowledge of the paleoenvironments of these regions has remained limited. We show that DNA and amino acids from buried organisms can be recovered from the basal sections of deep ice cores, enabling reconstructions of past flora and fauna. We show that high-altitude southern Greenland, currently lying below more than 2 kilometers of ice, was inhabited by a diverse array of conifer trees and insects within the past million years. The results provide direct evidence in support of a forested southern Greenland and suggest that many deep ice cores may contain genetic records of paleoenvironments in their basal sections.

Udgivelsesdato: 2007-Jul-6, It is difficult to obtain fossil data from the 10% of Earth's terrestrial surface that is covered by thick glaciers and ice sheets, and hence, knowledge of the paleoenvironments of these regions has remained limited. We show that DNA and amino acids from buried organisms can be recovered from the basal sections of deep ice cores, enabling reconstructions of past flora and fauna. We show that high-altitude southern Greenland, currently lying below more than 2 kilometers of ice, was inhabited by a diverse array of conifer trees and insects within the past million years. The results provide direct evidence in support of a forested southern Greenland and suggest that many deep ice cores may contain genetic records of paleoenvironments in their basal sections.

We have used a systematic protocol for extracting, quantitating, sexing and validating ancient human mitochondrial and nuclear DNA of one male and one female Beothuk, a Native American population from Newfoundland, which became extinct approximately 180 years ago. They carried mtDNA haplotypes, which fall within haplogroups X and C, consistent with Northeastern Native populations today. In addition we have sexed the male using a novel-sexing assay and confirmed the authenticity of his Y chromosome with the presence of the Native American specific Y-QM3 single nucleotide polymorphism (SNP). This is the first ancient nuclear SNP typed from a Native population in the Americas. In addition, using the same teeth we conducted a stable isotopes analysis of collagen and dentine to show that both individuals relied on marine sources (fresh and salt water fish, seals) with no hierarchy seen between them, and that their water sources were pooled or stored water. Both mtDNA sequence data and Y SNP data hint at possible gene flow or a common ancestral population for both the Beothuk and the current day Mikmaq, but more importantly the data do not lend credence to the proposed idea that the Beothuk (specifically, Nonosabasut) were of admixed (European-Native American) descent. We also analyzed patterns of DNA damage in the clones of authentic mtDNA sequences; there is no tendency for DNA damage to occur preferentially at previously defined mutational hotspots, suggesting that such mutational hotspots are not hypervariable because they are more prone to damage. Udgivelsesdato: 2007-Apr

Udgivelsesdato: January 2007, Although ancient DNA (aDNA) miscoding lesions have been studied since the earliest days of the field, their nature remains a source of debate. A variety of conflicting hypotheses exist about which miscoding lesions constitute true aDNA damage as opposed to PCR polymerase amplification error. Furthermore, considerable disagreement and speculation exists on which specific damage events underlie observed miscoding lesions. The root of the problem is that it has previously been difficult to assemble sufficient data to test the hypotheses, and near-impossible to accurately determine the specific strand of origin of observed damage events. With the advent of emulsion-based clonal amplification (emPCR) and the sequencing-by-synthesis technology this has changed. In this paper we demonstrate how data produced on the Roche GS20 genome sequencer can determine miscoding lesion strands of origin, and subsequently be interpreted to enable characterization of the aDNA damage behind the observed phenotypes. Through comparative analyses on 390,965 bp of modern chloroplast and 131,474 bp of ancient woolly mammoth GS20 sequence data we conclusively demonstrate that in this sample at least, a permafrost preserved specimen, Type 2 (cytosine-->thymine/guanine-->adenine) miscoding lesions represent the overwhelming majority of damage-derived miscoding lesions. Additionally, we show that an as yet unidentified guanine-->adenine analogue modification, not the conventionally argued cytosine-->uracil deamination, underpins a significant proportion of Type 2 damage. How widespread these implications are for aDNA will become apparent as future studies analyse data recovered from a wider range of substrates.