The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Ancient Clostridium DNA and variants of tetanus neurotoxins associated with human archaeological remains.

The analysis of microbial genomes from human archaeological samples offers a historic snapshot of ancient pathogens and provides insights into the origins of modern infectious diseases. Here, we analyze metagenomic datasets from 38 human archaeological samples and identify bacterial genomic sequences related to modern-day Clostridium tetani, which produces the tetanus neurotoxin (TeNT) and causes the disease tetanus. These genomic assemblies had varying levels of completeness, and a subset of them displayed hallmarks of ancient DNA damage. Phylogenetic analyses revealed known C. tetani clades as well as potentially new Clostridium lineages closely related to C. tetani. The genomic assemblies encode 13 TeNT variants with unique substitution profiles, including a subgroup of TeNT variants found exclusively in ancient samples from South America. We experimentally tested a TeNT variant selected from an ancient Chilean mummy sample and found that it induced tetanus muscle paralysis in mice, with potency comparable to modern TeNT. Thus, our ancient DNA analysis identifies DNA from neurotoxigenic C. tetani in archaeological human samples, and a novel variant of TeNT that can cause disease in mammals.

A 14th century CE Brucella melitensis genome and the recent expansion of the Western Mediterranean clade.

Brucellosis is a disease caused by the bacterium Brucella and typically transmitted through contact with infected ruminants. It is one of the most common chronic zoonotic diseases and of particular interest to public health agencies. Despite its well-known transmission history and characteristic symptoms, we lack a more complete understanding of the evolutionary history of its best-known species-Brucella melitensis. To address this knowledge gap we fortuitously found, sequenced and assembled a high-quality ancient B. melitensis draft genome from the kidney stone of a 14th-century Italian friar. The ancient strain contained fewer core genes than modern B. melitensis isolates, carried a complete complement of virulence genes, and did not contain any indication of significant antimicrobial resistances. The ancient B. melitensis genome fell as a basal sister lineage to a subgroup of B. melitensis strains within the Western Mediterranean phylogenetic group, with a short branch length indicative of its earlier sampling time, along with a similar gene content. By calibrating the molecular clock we suggest that the speciation event between B. melitensis and B. abortus is contemporaneous with the estimated time frame for the domestication of both sheep and goats. These results confirm the existence of the Western Mediterranean clade as a separate group in the 14th CE and suggest that its divergence was due to human and ruminant co-migration.

The DataHarmonizer: a tool for faster data harmonization, validation, aggregation and analysis of pathogen genomics contextual information.

Pathogen genomics is a critical tool for public health surveillance, infection control, outbreak investigations as well as research. In order to make use of pathogen genomics data, they must be interpreted using contextual data (metadata). Contextual data include sample metadata, laboratory methods, patient demographics, clinical outcomes and epidemiological information. However, the variability in how contextual information is captured by different authorities and how it is encoded in different databases poses challenges for data interpretation, integration and their use/re-use. The DataHarmonizer is a template-driven spreadsheet application for harmonizing, validating and transforming genomics contextual data into submission-ready formats for public or private repositories. The tool's web browser-based JavaScript environment enables validation and its offline functionality and local installation increases data security. The DataHarmonizer was developed to address the data sharing needs that arose during the COVID-19 pandemic, and was used by members of the Canadian COVID Genomics Network (CanCOGeN) to harmonize SARS-CoV-2 contextual data for national surveillance and for public repository submission. In order to support coordination of international surveillance efforts, we have partnered with the Public Health Alliance for Genomic Epidemiology to also provide a template conforming to its SARS-CoV-2 contextual data specification for use worldwide. Templates are also being developed for One Health and foodborne pathogens. Overall, the DataHarmonizer tool improves the effectiveness and fidelity of contextual data capture as well as its subsequent usability. Harmonization of contextual information across authorities, platforms and systems globally improves interoperability and reusability of data for concerted public health and research initiatives to fight the current pandemic and future public health emergencies. While initially developed for the COVID-19 pandemic, its expansion to other data management applications and pathogens is already underway.

Monkeypox in Montréal: Epidemiology, Phylogenomics, and Public Health Response to a Large North American Outbreak.

BACKGROUND: Monkeypox, a viral zoonotic disease, is causing a global outbreak outside of endemic areas. OBJECTIVE: To characterize the outbreak of monkeypox in Montréal, the first large outbreak in North America. DESIGN: Epidemiologic and laboratory surveillance data and a phylogenomic analysis were used to describe and place the outbreak in a global context. SETTING: Montréal, Canada. PATIENTS: Probable or confirmed cases of monkeypox. MEASUREMENTS: Epidemiologic, clinical, and demographic data were aggregated. Whole-genome sequencing and phylogenetic analysis were performed for a set of outbreak sequences. The public health response and its evolution are described. RESULTS: Up to 18 October 2022, a total of 402 cases of monkeypox were reported mostly among men who have sex with men (MSM), most of which were suspected to be acquired through sexual contact. All monkeypox genomes nested within the B.1 lineage. Montréal Public Health worked closely with the affected communities to control the outbreak, becoming the first jurisdiction to offer 1 dose of the Modified Vaccinia Ankara-Bavarian Nordic vaccine as preexposure prophylaxis (PrEP) to those at risk in early June 2022. Two peaks of cases were seen in early June and July (43 and 44 cases per week, respectively) followed by a decline toward near resolution of the outbreak in October. Reasons for the biphasic peak are not fully elucidated but may represent the tempo of vaccination and/or several factors related to transmission dynamics and case ascertainment. LIMITATIONS: Clinical data are self-reported. Limited divergence among sequences limited genomic epidemiologic conclusions. CONCLUSION: A large outbreak of monkeypox occurred in Montréal, primarily among MSM. Successful control of the outbreak rested on early and sustained engagement with the affected communities and rapid offer of PrEP vaccination to at-risk persons. PRIMARY FUNDING SOURCE: None.

Examining pathogen DNA recovery across the remains of a 14th century Italian friar (Blessed Sante) infected with Brucella melitensis.

OBJECTIVE: To investigate variation in ancient DNA recovery of Brucella melitensis, the causative agent of brucellosis, from multiple tissues belonging to one individual MATERIALS: 14 samples were analyzed from the mummified remains of the Blessed Sante, a 14 th century Franciscan friar from central Italy, with macroscopic diagnosis of probable brucellosis. METHODS: Shotgun sequencing data from was examined to determine the presence of Brucella DNA. RESULTS: Three of the 14 samples contained authentic ancient DNA, identified as belonging to B. melitensis. A genome (23.81X depth coverage, 0.98 breadth coverage) was recovered from a kidney stone. Nine of the samples contained reads classified as B. melitensis (7-169), but for many the data quality was insufficient to withstand our identification and authentication criteria. CONCLUSIONS: We identified significant variation in the preservation and abundance of B. melitensis DNA present across multiple tissues, with calcified nodules yielding the highest number of authenticated reads. This shows how greatly sample selection can impact pathogen identification. SIGNIFICANCE: Our results demonstrate variation in the preservation and recovery of pathogen DNA across tissues. This study highlights the importance of sample selection in the reconstruction of infectious disease burden and highlights the importance of a holistic approach to identifying disease. LIMITATIONS: Study focuses on pathogen recovery in a single individual. SUGGESTIONS FOR FURTHER RESEARCH: Further analysis of how sampling impacts aDNA recovery will improve pathogen aDNA recovery and advance our understanding of disease in past peoples.

A 16th century Escherichia coli draft genome associated with an opportunistic bile infection.

Escherichia coli - one of the most characterized bacteria and a major public health concern - remains invisible across the temporal landscape. Here, we present the meticulous reconstruction of the first ancient E. coli genome from a 16th century gallstone from an Italian mummy with chronic cholecystitis. We isolated ancient DNA and reconstructed the ancient E. coli genome. It consisted of one chromosome of 4446 genes and two putative plasmids with 52 genes. The E. coli strain belonged to the phylogroup A and an exceptionally rare sequence type 4995. The type VI secretion system component genes appears to be horizontally acquired from Klebsiella aerogenes, however we could not identify any pathovar specific genes nor any acquired antibiotic resistances. A sepsis mouse assay showed that a closely related contemporary E. coli strain was avirulent. Our reconstruction of this ancient E. coli helps paint a more complete picture of the burden of opportunistic infections of the past.

Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments.

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.1 In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.2 Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia3-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,6-12 in paleodietary reconstructions,13-19 and facilitated genomic reconstructions from cave contexts.8,20-22 Murchie et al.6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.24 Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 40,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

SARS-CoV-2 infection and transmission in the North American deer mouse.

Widespread circulation of SARS-CoV-2 in humans raises the theoretical risk of reverse zoonosis events with wildlife, reintroductions of SARS-CoV-2 into permissive nondomesticated animals. Here we report that North American deer mice (Peromyscus maniculatus) are susceptible to SARS-CoV-2 infection following intranasal exposure to a human isolate, resulting in viral replication in the upper and lower respiratory tract with little or no signs of disease. Further, shed infectious virus is detectable in nasal washes, oropharyngeal and rectal swabs, and viral RNA is detectable in feces and occasionally urine. We further show that deer mice are capable of transmitting SARS-CoV-2 to naïve deer mice through direct contact. The extent to which these observations may translate to wild deer mouse populations remains unclear, and the risk of reverse zoonosis and/or the potential for the establishment of Peromyscus rodents as a North American reservoir for SARS-CoV-2 remains unknown.

Response to Brinkmann et al. "Re-assembly of 19th century smallpox vaccine genomes reveals the contemporaneous use of horsepox and horsepox-related viruses in the United States".

We thank Brinkmann and colleagues for their correspondence and their further investigation into these American Civil War Era vaccination strains. Here, we summarize the difficulties and caveats of work with ancient DNA.

American mastodon mitochondrial genomes suggest multiple dispersal events in response to Pleistocene climate oscillations.

Pleistocene glacial-interglacial cycles are correlated with dramatic temperature oscillations. Examining how species responded to these natural fluctuations can provide valuable insights into the impacts of present-day anthropogenic climate change. Here we present a phylogeographic study of the extinct American mastodon (Mammut americanum), based on 35 complete mitochondrial genomes. These data reveal the presence of multiple lineages within this species, including two distinct clades from eastern Beringia. Our molecular date estimates suggest that these clades arose at different times, supporting a pattern of repeated northern expansion and local extirpation in response to glacial cycling. Consistent with this hypothesis, we also note lower levels of genetic diversity among northern mastodons than in endemic clades south of the continental ice sheets. The results of our study highlight the complex relationships between population dispersals and climate change, and can provide testable hypotheses for extant species expected to experience substantial biogeographic impacts from rising temperatures.

The origins and genomic diversity of American Civil War Era smallpox vaccine strains.

Vaccination has transformed public health, most notably including the eradication of smallpox. Despite its profound historical importance, little is known of the origins and diversity of the viruses used in smallpox vaccination. Prior to the twentieth century, the method, source and origin of smallpox vaccinations remained unstandardised and opaque. We reconstruct and analyse viral vaccine genomes associated with smallpox vaccination from historical artefacts. Significantly, we recover viral molecules through non-destructive sampling of historical materials lacking signs of biological residues. We use the authenticated ancient genomes to reveal the evolutionary relationships of smallpox vaccination viruses within the poxviruses as a whole.

Genetic resiliency and the Black Death: No apparent loss of mitogenomic diversity due to the Black Death in medieval London and Denmark.

OBJECTIVES: In the 14th century AD, medieval Europe was severely affected by the Great European Famine as well as repeated bouts of disease, including the Black Death, causing major demographic shifts. This high volatility led to increased mobility and migration due to new labor and economic opportunities, as evidenced by documentary and stable isotope data. This study uses ancient DNA (aDNA) isolated from skeletal remains to examine whether evidence for large-scale population movement can be gleaned from the complete mitochondrial genomes of 264 medieval individuals from England (London) and Denmark. MATERIALS AND METHODS: Using a novel library-conserving approach to targeted capture, we recovered 264 full mitochondrial genomes from the petrous portion of the temporal bones and teeth and compared genetic diversity across the medieval period within and between English (London) and Danish populations and with contemporary populations through population pairwise ΦST analysis. RESULTS: We find no evidence of significant differences in genetic diversity spatially or temporally in our dataset, yet there is a high degree of haplotype diversity in our medieval samples with little exact sequence sharing. DISCUSSION: The mitochondrial genomes of both medieval Londoners and medieval Danes suggest high mitochondrial diversity before, during and after the Black Death. While our mitochondrial genomic data lack geographically correlated signals, these data could be the result of high, continual female migration before and after the Black Death or may simply indicate a large female effective population size unaffected by the upheaval of the medieval period. Either scenario suggests a genetic resiliency in areas of northwestern medieval Europe.

Dorset Pre-Inuit and Beothuk foodways in Newfoundland, ca. AD 500-1829.

Archaeological research on the Canadian island of Newfoundland increasingly demonstrates that the island's subarctic climate and paucity of terrestrial food resources did not restrict past Pre-Inuit (Dorset) and Native American (Beothuk) hunter-gatherer populations to a single subsistence pattern. This study first sought to characterize hunter-gatherer diets over the past 1500 years; and second, to assess the impact of European colonization on Beothuk lifeways by comparing the bone chemistry of Beothuk skeletal remains before and after the intensification of European settlement in the early 18th century. We employed radiocarbon dating and stable carbon and nitrogen isotope ratio analysis of bulk bone collagen from both Dorset (n = 9) and Beothuk (n = 13) cultures, including a naturally mummified 17th century Beothuk individual. Carbon and nitrogen isotope analysis of 108 faunal samples from Dorset and Beothuk archaeological sites around the island were used as a dietary baseline for the humans. We combined our results with previously published isotope data and radiocarbon dates from Dorset (n = 12) and Beothuk (n = 18) individuals and conducted a palaeodietary analysis using Bayesian modelling, cluster analysis and comparative statistical tests. Dorset diets featured more marine protein than those of the Beothuk, and the diets of Beothuk after the 18th century featured less high trophic level marine protein than those of individuals predating the 18th century. Despite inhabiting the same island, Dorset and Beothuk cultures employed markedly different dietary strategies, consistent with interpretations of other archaeological data. Significantly, European colonization had a profound effect on Beothuk lifeways, as in response to the increasing European presence on the coast, the Beothuk relied more extensively on the limited resources of the island's boreal forests and rivers.

Resolving the phylogenetic position of Darwin's extinct ground sloth (Mylodon darwinii) using mitogenomic and nuclear exon data.

Mylodon darwinii is the extinct giant ground sloth named after Charles Darwin, who first collected its remains in South America. We have successfully obtained a high-quality mitochondrial genome at 99-fold coverage using an Illumina shotgun sequencing of a 12 880-year-old bone fragment from Mylodon Cave in Chile. Low level of DNA damage showed that this sample was exceptionally well preserved for an ancient subfossil, probably the result of the dry and cold conditions prevailing within the cave. Accordingly, taxonomic assessment of our shotgun metagenomic data showed a very high percentage of endogenous DNA with 22% of the assembled metagenomic contigs assigned to Xenarthra. Additionally, we enriched over 15 kb of sequence data from seven nuclear exons, using target sequence capture designed against a wide xenarthran dataset. Phylogenetic and dating analyses of the mitogenomic dataset including all extant species of xenarthrans and the assembled nuclear supermatrix unambiguously place Mylodon darwinii as the sister-group of modern two-fingered sloths, from which it diverged around 22 million years ago. These congruent results from both the mitochondrial and nuclear data support the diphyly of the two modern sloth lineages, implying the convergent evolution of their unique suspensory behaviour as an adaption to arboreality. Our results offer promising perspectives for whole-genome sequencing of this emblematic extinct taxon.

Correction: The paradox of HBV evolution as revealed from a 16th century mummy.

[This corrects the article DOI: 10.1371/journal.ppat.1006750.].

The Gateway from Near into Remote Oceania: New Insights from Genome-Wide Data.

A widely accepted two-wave scenario of human settlement of Oceania involves the first out-of-Africa migration circa 50,000 years ago (ya), and the more recent Austronesian expansion, which reached the Bismarck Archipelago by 3,450 ya. Whereas earlier genetic studies provided evidence for extensive sex-biased admixture between the incoming and the indigenous populations, some archaeological, linguistic, and genetic evidence indicates a more complicated picture of settlement. To study regional variation in Oceania in more detail, we have compiled a genome-wide data set of 823 individuals from 72 populations (including 50 populations from Oceania) and over 620,000 autosomal single nucleotide polymorphisms (SNPs). We show that the initial dispersal of people from the Bismarck Archipelago into Remote Oceania occurred in a "leapfrog" fashion, completely by-passing the main chain of the Solomon Islands, and that the colonization of the Solomon Islands proceeded in a bidirectional manner. Our results also support a divergence between western and eastern Solomons, in agreement with the sharp linguistic divide known as the Tryon-Hackman line. We also report substantial post-Austronesian gene flow across the Solomons. In particular, Santa Cruz (in Remote Oceania) exhibits extraordinarily high levels of Papuan ancestry that cannot be explained by a simple bottleneck/founder event scenario. Finally, we use simulations to show that discrepancies between different methods for dating admixture likely reflect different sensitivities of the methods to multiple admixture events from the same (or similar) sources. Overall, this study points to the importance of fine-scale sampling to understand the complexities of human population history.

The paradox of HBV evolution as revealed from a 16th century mummy.

Hepatitis B virus (HBV) is a ubiquitous viral pathogen associated with large-scale morbidity and mortality in humans. However, there is considerable uncertainty over the time-scale of its origin and evolution. Initial shotgun data from a mid-16th century Italian child mummy, that was previously paleopathologically identified as having been infected with Variola virus (VARV, the agent of smallpox), showed no DNA reads for VARV yet did for hepatitis B virus (HBV). Previously, electron microscopy provided evidence for the presence of VARV in this sample, although similar analyses conducted here did not reveal any VARV particles. We attempted to enrich and sequence for both VARV and HBV DNA. Although we did not recover any reads identified as VARV, we were successful in reconstructing an HBV genome at 163.8X coverage. Strikingly, both the HBV sequence and that of the associated host mitochondrial DNA displayed a nearly identical cytosine deamination pattern near the termini of DNA fragments, characteristic of an ancient origin. In contrast, phylogenetic analyses revealed a close relationship between the putative ancient virus and contemporary HBV strains (of genotype D), at first suggesting contamination. In addressing this paradox we demonstrate that HBV evolution is characterized by a marked lack of temporal structure. This confounds attempts to use molecular clock-based methods to date the origin of this virus over the time-frame sampled so far, and means that phylogenetic measures alone cannot yet be used to determine HBV sequence authenticity. If genuine, this phylogenetic pattern indicates that the genotypes of HBV diversified long before the 16th century, and enables comparison of potential pathogenic similarities between modern and ancient HBV. These results have important implications for our understanding of the emergence and evolution of this common viral pathogen.

Successful reconstruction of whole mitochondrial genomes from ancient Central America and Mexico.

The northern and southern peripheries of ancient Mesoamerica are poorly understood. There has been speculation over whether borderland cultures such as Greater Nicoya and Casas Grandes represent Mesoamerican outposts in the Isthmo-Colombian area and the Greater Southwest, respectively. Poor ancient DNA preservation in these regions challenged previous attempts to resolve these questions using conventional genetic techniques. We apply advanced in-solution mitogenome capture and high-throughput sequencing to fourteen dental samples obtained from the Greater Nicoya sites of Jícaro and La Cascabel in northwest Costa Rica (n = 9; A.D. 800-1250) and the Casas Grandes sites of Paquimé and Convento in northwest Mexico (n = 5; A.D. 1200-1450). Full mitogenome reconstruction was successful for three individuals from Jícaro and five individuals from Paquimé and Convento. The three Jícaro individuals belong to haplogroup B2d, a haplogroup found today only among Central American Chibchan-speakers. The five Paquimé and Convento individuals belong to haplogroups C1c1a, C1c5, B2f and B2a which, are found in contemporary populations in North America and Mesoamerica. We report the first successfully reconstructed ancient mitogenomes from Central America, and the first genetic evidence of ancestry affinity of the ancient inhabitants of Greater Nicoya and Casas Grandes with contemporary Isthmo-Columbian and Greater Southwest populations, respectively.

Genetic Discontinuity between the Maritime Archaic and Beothuk Populations in Newfoundland, Canada.

Situated at the furthest northeastern edge of Canada, the island of Newfoundland (approximately 110,000 km2) and Labrador (approximately 295,000 km2) today constitute a province characterized by abundant natural resources but low population density. Both landmasses were covered by the Laurentide ice sheet during the Last Glacial Maximum (18,000 years before present [YBP]); after the glacier retreated, ice patches remained on the island until ca. 9,000 calibrated (cal) YBP [1]. Nevertheless, indigenous peoples, whose ancestors had trekked some 5,000 km from the west coast, arrived approximately 10,000 cal YBP in Labrador and ca. 6,000 cal YBP in Newfoundland [2, 3]. Differential features in material culture indicate at least three settlement episodes by distinct cultural groups, including the Maritime Archaic, Palaeoeskimo, and Beothuk. Newfoundland has remained home to indigenous peoples until present day with only one apparent hiatus (3,400-2,800 YBP). This record suggests abandonment, severe constriction, or local extinction followed by subsequent immigrations from single or multiple source populations, but the specific dynamics and the cultural and biological relationships, if any, among these successive peoples remain enigmatic [4]. By examining the mitochondrial genome diversity and isotopic ratios of 74 ancient remains in conjunction with the archaeological record, we have provided definitive evidence for the genetic discontinuity between the maternal lineages of these populations. This northeastern margin of North America appears to have been populated multiple times by distinct groups that did not share a recent common ancestry, but rather one much deeper in time at the entry point into the continent.