Capturing the Resistome: a Targeted Capture Method To Reveal Antibiotic Resistance Determinants in Metagenomes.

Identification of the nucleotide sequences encoding antibiotic resistance elements and determination of their association with antibiotic resistance are critical to improve surveillance and monitor trends in antibiotic resistance. Current methods to study antibiotic resistance in various environments rely on extensive deep sequencing or laborious culturing of fastidious organisms, both of which are heavily time-consuming operations. An accurate and sensitive method to identify both rare and common resistance elements in complex metagenomic samples is needed. Referencing the sequences in the Comprehensive Antibiotic Resistance Database, we designed a set of 37,826 probes to specifically target over 2,000 nucleotide sequences associated with antibiotic resistance in clinically relevant bacteria. Testing of this probe set on DNA libraries generated from multidrug-resistant bacteria to selectively capture resistance genes reproducibly produced higher numbers of reads on target at a greater length of coverage than shotgun sequencing. We also identified additional resistance gene sequences from human gut microbiome samples that sequencing alone was not able to detect. Our method to capture the resistome enables a sensitive means of gene detection in diverse environments where genes encoding antibiotic resistance represent less than 0.1% of the metagenome.

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.

Shotgun Mitogenomics Provides a Reference Phylogenetic Framework and Timescale for Living Xenarthrans.

Xenarthra (armadillos, sloths, and anteaters) constitutes one of the four major clades of placental mammals. Despite their phylogenetic distinctiveness in mammals, a reference phylogeny is still lacking for the 31 described species. Here we used Illumina shotgun sequencing to assemble 33 new complete mitochondrial genomes, establishing Xenarthra as the first major placental clade to be fully sequenced at the species level for mitogenomes. The resulting data set allowed the reconstruction of a robust phylogenetic framework and timescale that are consistent with previous studies conducted at the genus level using nuclear genes. Incorporating the full species diversity of extant xenarthrans points to a number of inconsistencies in xenarthran systematics and species definition. We propose to split armadillos into two distinct families Dasypodidae (dasypodines) and Chlamyphoridae (euphractines, chlamyphorines, and tolypeutines) to better reflect their ancient divergence, estimated around 42 Ma. Species delimitation within long-nosed armadillos (genus Dasypus) appeared more complex than anticipated, with the discovery of a divergent lineage in French Guiana. Diversification analyses showed Xenarthra to be an ancient clade with a constant diversification rate through time with a species turnover driven by high but constant extinction. We also detected a significant negative correlation between speciation rate and past temperature fluctuations with an increase in speciation rate corresponding to the general cooling observed during the last 15 My. Biogeographic reconstructions identified the tropical rainforest biome of Amazonia and the Guiana Shield as the cradle of xenarthran evolutionary history with subsequent dispersions into more open and dry habitats.

Second-pandemic strain of Vibrio cholerae from the Philadelphia cholera outbreak of 1849.

In the 19th century, there were several major cholera pandemics in the Indian subcontinent, Europe, and North America. The causes of these outbreaks and the genomic strain identities remain a mystery. We used targeted high-throughput sequencing to reconstruct the Vibrio cholerae genome from the preserved intestine of a victim of the 1849 cholera outbreak in Philadelphia, part of the second cholera pandemic. This O1 biotype strain has 95 to 97% similarity with the classical O395 genome, differing by 203 single-nucleotide polymorphisms (SNPs), lacking three genomic islands, and probably having one or more tandem cholera toxin prophage (CTX) arrays, which potentially affected its virulence. This result highlights archived medical remains as a potential resource for investigations into the genomic origins of past pandemics.

Ancient whole genome enrichment using baits built from modern DNA.

We report metrics from complete genome capture of nuclear DNA from extinct mammoths using biotinylated RNAs transcribed from an Asian elephant DNA extract. Enrichment of the nuclear genome ranged from 1.06- to 18.65-fold, to an apparent maximum threshold of ∼80% on-target. This projects an order of magnitude less costly complete genome sequencing from long-dead organisms, even when a reference genome is unavailable for bait design.

Correction to 'Ancient DNA and the tropics: a rodent's tale'.

The present erratum is in regards to our article entitled 'Ancient DNA and the tropics: a rodent's tale'. We were made aware of problems with some of the ancient sequences submitted to GenBank and conducted a systematic review of all the files used in our study. We discovered that, unfortunately, an incorrect file was sent to GenBank and was also used in some of our downstream analyses. We immediately contacted GenBank, explained the situation and corrected the file. We have redone some analyses with the correct file and describe these changes below.

Ancient DNA and the tropics: a rodent's tale.

Most genetic studies of Holocene fauna have been performed with ancient samples from dry and cold regions, in which preservation of fossils is facilitated and molecular damage is reduced. Ancient DNA work from tropical regions has been precluded owing to factors that limit DNA preservation (e.g. temperature, hydrolytic damage). We analysed ancient DNA from rodent jawbones identified as Ototylomys phyllotis, found in Holocene and Late Pleistocene stratigraphic layers from Loltún, a humid tropical cave located in the Yucatan peninsula. We extracted DNA and amplified six short overlapping fragments of the cytochrome b gene, totalling 666 bp, which represents an unprecedented success considering tropical ancient DNA samples. We performed genetic, phylogenetic and divergence time analyses, combining sequences from ancient and modern O. phyllotis, in order to assess the ancestry of the Loltún samples. Results show that all ancient samples fall into a unique clade that diverged prior to the divergence of the modern O. phyllotis, supporting it as a distinct Pleistocene form of the Ototylomys genus. Hence, this rodent's tale suggests that the sister group to modern O. phyllotis arose during the Miocene-Pliocene, diversified during the Pleistocene and went extinct in the Holocene.

Erratum: Probe design for simultaneous, targeted capture of diverse metagenomic targets.

[This corrects the article DOI: 10.1016/j.crmeth.2021.100069.].

Ancient DNA from lake sediments: bridging the gap between paleoecology and genetics.

BACKGROUND: Quaternary plant ecology in much of the world has historically relied on morphological identification of macro- and microfossils from sediments of small freshwater lakes. Here, we report new protocols that reliably yield DNA sequence data from Holocene plant macrofossils and bulk lake sediment used to infer ecological change. This will allow changes in census populations, estimated from fossils and associated sediment, to be directly associated with population genetic changes. RESULTS: We successfully sequenced DNA from 64 samples (out of 126) comprised of bulk sediment and seeds, leaf fragments, budscales, and samaras extracted from Holocene lake sediments in the western Great Lakes region of North America. Overall, DNA yields were low. However, we were able to reliably amplify samples with as few as 10 copies of a short cpDNA fragment with little detectable PCR inhibition. Our success rate was highest for sediments < 2000 years old, but we were able to successfully amplify DNA from samples up to 4600 years old. DNA sequences matched the taxonomic identity of the macrofossil from which they were extracted 79% of the time. Exceptions suggest that DNA molecules from surrounding nearby sediments may permeate or adhere to macrofossils in sediments. CONCLUSIONS: An ability to extract ancient DNA from Holocene sediments potentially allows exciting new insights into the genetic consequences of long-term environmental change. The low DNA copy numbers we found in fossil material and the discovery of multiple sequence variants from single macrofossil extractions highlight the need for careful experimental and laboratory protocols. Further application of these protocols should lead to better understanding of the ecological and evolutionary consequences of environmental change.

Out of America: ancient DNA evidence for a new world origin of late quaternary woolly mammoths.

Although the iconic mammoth of the Late Pleistocene, the woolly mammoth (Mammuthus primigenius), has traditionally been regarded as the end point of a single anagenetically evolving lineage, recent paleontological and molecular studies have shown that successive allopatric speciation events must have occurred within Pleistocene Mammuthus in Asia, with subsequent expansion and hybridization between nominal taxa [1, 2]. However, the role of North American mammoth populations in these events has not been adequately explored from an ancient-DNA standpoint. To undertake this task, we analyzed mtDNA from a large data set consisting of mammoth samples from across Holarctica (n = 160) and representing most of radiocarbon time. Our evidence shows that, during the terminal Pleistocene, haplotypes originating in and characteristic of New World populations replaced or succeeded those endemic to Asia and western Beringia. Also, during the Last Glacial Maximum, mammoth populations do not appear to have suffered an overall decline in diversity, despite differing responses on either side of the Bering land bridge. In summary, the "Out-of-America" hypothesis holds that the dispersal of North American woolly mammoths into other parts of Holarctica created major phylogeographic structuring within Mammuthus primigenius populations, shaping the last phase of their evolutionary history before their demise.

New insights from old bones: DNA preservation and degradation in permafrost preserved mammoth remains.

Despite being plagued by heavily degraded DNA in palaeontological remains, most studies addressing the state of DNA degradation have been limited to types of damage which do not pose a hindrance to Taq polymerase during PCR. Application of serial qPCR to the two fractions obtained during extraction (demineralization and protein digest) from six permafrost mammoth bones and one partially degraded modern elephant bone has enabled further insight into the changes which endogenous DNA is subjected to during diagenesis. We show here that both fractions exhibit individual qualities in terms of the prevailing type of DNA (i.e. mitochondrial versus nuclear DNA) as well as the extent of damage, and in addition observed a highly variable ratio of mitochondrial to nuclear DNA among the six mammoth samples. While there is evidence suggesting that mitochondrial DNA is better preserved than nuclear DNA in ancient permafrost samples, we find the initial DNA concentration in the bone tissue to be as relevant for the total accessible mitochondrial DNA as the extent of DNA degradation post-mortem. We also evaluate the general applicability of indirect measures of preservation such as amino-acid racemization, bone crystallinity index and thermal age to these exceptionally well-preserved samples.

Probe design for simultaneous, targeted capture of diverse metagenomic targets.

The compounding challenges of low signal, high background, and uncertain targets plague many metagenomic sequencing efforts. One solution has been DNA capture, wherein probes are designed to hybridize with target sequences, enriching them in relation to their background. However, balancing probe depth with breadth of capture is challenging for diverse targets. To find this balance, we have developed the HUBDesign pipeline, which makes use of sequence homology to design probes at multiple taxonomic levels. This creates an efficient probe set capable of simultaneously and specifically capturing known and related sequences. We validated HUBDesign by generating probe sets targeting the breadth of coronavirus diversity, as well as a suite of bacterial pathogens often underlying sepsis. In separate experiments demonstrating significant, simultaneous enrichment, we captured SARS-CoV-2 and HCoV-NL63 in a human RNA background and seven bacterial strains in human blood. HUBDesign (https://github.com/zacherydickson/HUBDesign) has broad applicability wherever there are multiple organisms of interest.

Comparison of methods in the recovery of nucleic acids from archival formalin-fixed paraffin-embedded autopsy tissues.

Archival formalin-fixed paraffin-embedded (FFPE) human tissue collections are typically in poor states of storage across the developing world. With advances in biomolecular techniques, these extraordinary and virtually untapped resources have become an essential part of retrospective epidemiological studies. To successfully use such tissues in genomic studies, scientists require high nucleic acid yields and purity. In spite of the increasing number of FFPE tissue kits available, few studies have analyzed their applicability in recovering high-quality nucleic acids from archived human autopsy samples. Here we provide a study involving 10 major extraction methods used to isolate total nucleic acid from FFPE tissues ranging in age from 3 to 13years. Although all 10 methods recovered quantifiable amounts of DNA, only 6 recovered quantifiable RNA, varying considerably and generally yielding lower DNA concentrations. Overall, we show quantitatively that TrimGen's WaxFree method and our in-house phenol-chloroform extraction method recovered the highest yields of amplifiable DNA, with considerable polymerase chain reaction (PCR) inhibition, whereas Ambion's RecoverAll method recovered the most amplifiable RNA.

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.

Ancient Mitogenomes Reveal the Evolutionary History and Biogeography of Sloths.

Living sloths represent two distinct lineages of small-sized mammals that independently evolved arboreality from terrestrial ancestors. The six extant species are the survivors of an evolutionary radiation marked by the extinction of large terrestrial forms at the end of the Quaternary. Until now, sloth evolutionary history has mainly been reconstructed from phylogenetic analyses of morphological characters. Here, we used ancient DNA methods to successfully sequence 10 extinct sloth mitogenomes encompassing all major lineages. This includes the iconic continental ground sloths Megatherium, Megalonyx, Mylodon, and Nothrotheriops and the smaller endemic Caribbean sloths Parocnus and Acratocnus. Phylogenetic analyses identify eight distinct lineages grouped in three well-supported clades, whose interrelationships are markedly incongruent with the currently accepted morphological topology. We show that recently extinct Caribbean sloths have a single origin but comprise two highly divergent lineages that are not directly related to living two-fingered sloths, which instead group with Mylodon. Moreover, living three-fingered sloths do not represent the sister group to all other sloths but are nested within a clade of extinct ground sloths including Megatherium, Megalonyx, and Nothrotheriops. Molecular dating also reveals that the eight newly recognized sloth families all originated between 36 and 28 million years ago (mya). The early divergence of recently extinct Caribbean sloths around 35 mya is consistent with the debated GAARlandia hypothesis postulating the existence at that time of a biogeographic connection between northern South America and the Greater Antilles. This new molecular phylogeny has major implications for reinterpreting sloth morphological evolution, biogeography, and diversification history.

Nuclear gene sequences from a late pleistocene sloth coprolite.

The determination of nuclear DNA sequences from ancient remains would open many novel opportunities such as the resolution of phylogenies, the sexing of hominid and animal remains, and the characterization of genes involved in phenotypic traits. However, to date, single-copy nuclear DNA sequences from fossils have been determined only from bones and teeth of woolly mammoths preserved in the permafrost. Since the best preserved ancient nucleic acids tend to stem from cold environments, this has led to the assumption that nuclear DNA would be retrievable only from frozen remains. We have previously shown that Pleistocene coprolites stemming from the extinct Shasta sloth (Nothrotheriops shastensis, Megatheriidae) contain mitochondrial (mt) DNA from the animal that produced them as well as chloroplast (cp) DNA from the ingested plants. Recent attempts to resolve the phylogeny of two families of extinct sloths by using strictly mitochondrial DNA has been inconclusive. We have prepared DNA extracts from a ground sloth coprolite from Gypsum Cave, Nevada, and quantitated the number of mtDNA copies for three different fragment lengths by using real-time PCR. We amplified one multicopy and three single-copy nuclear gene fragments and used the concatenated sequence to resolve the phylogeny. These results show that ancient single-copy nuclear DNA can be recovered from warm, arid climates. Thus, nuclear DNA preservation is not restricted to cold climates.

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.

A molecular portrait of maternal sepsis from Byzantine Troy.

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman's remains. Scanning electron microscopy of the tissue revealed 'ghost cells', resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections.

Plasmodium falciparum malaria in 1st-2nd century CE southern Italy.

The historical record attests to the devastation malaria exacted on ancient civilizations, particularly the Roman Empire [1]. However, evidence for the presence of malaria during the Imperial period in Italy (1st-5th century CE) is based on indirect sources, such as historical, epigraphic, or skeletal evidence. Although these sources are crucial for revealing the context of this disease, they cannot establish the causative species of Plasmodium. Importantly, definitive evidence for the presence of malaria is now possible through the implementation of ancient DNA technology. As malaria is presumed to have been at its zenith during the Imperial period [1], we selected first or second molars from 58 adults from three cemeteries from this time: Isola Sacra (associated with Portus Romae, 1st-3rd century CE), Velia (1st-2nd century CE), and Vagnari (1st-4th century CE). We performed hybridization capture using baits designed from the mitochondrial (mtDNA) genomes of Plasmodium spp. on a prioritized subset of 11 adults (informed by metagenomic sequencing). The mtDNA sequences generated provided compelling phylogenetic evidence for the presence of P. falciparum in two individuals. This is the first genomic data directly implicating P. falciparum in Imperial period southern Italy in adults.

The phylogenetic affinities of the extinct glyptodonts.

Among the fossils of hitherto unknown mammals that Darwin collected in South America between 1832 and 1833 during the Beagle expedition were examples of the large, heavily armored herbivores later known as glyptodonts. Ever since, glyptodonts have fascinated evolutionary biologists because of their remarkable skeletal adaptations and seemingly isolated phylogenetic position even within their natural group, the cingulate xenarthrans (armadillos and their allies). In possessing a carapace comprised of fused osteoderms, the glyptodonts were clearly related to other cingulates, but their precise phylogenetic position as suggested by morphology remains unresolved. To provide a molecular perspective on this issue, we designed sequence-capture baits using in silico reconstructed ancestral sequences and successfully assembled the complete mitochondrial genome of Doedicurus sp., one of the largest glyptodonts. Our phylogenetic reconstructions establish that glyptodonts are in fact deeply nested within the armadillo crown-group, representing a distinct subfamily (Glyptodontinae) within family Chlamyphoridae. Molecular dating suggests that glyptodonts diverged no earlier than around 35 million years ago, in good agreement with their fossil record. Our results highlight the derived nature of the glyptodont morphotype, one aspect of which is a spectacular increase in body size until their extinction at the end of the last ice age.