Search

Background The identification of bona fide microbial taxa in microbiomes derived from ancient and historical samples is complicated by the unavoidable mixture between DNA from ante- and post-mortem microbial colonizers. One possibility to distinguish between these sources of microbial DNA is querying for the presence of age-associated degradation patterns typical of ancient DNA (aDNA). The presence of uracils, resulting from cytosine deamination, has been detected ubiquitously in aDNA retrieved from diverse sources, and used as an authentication criterion. Here, we employ a library preparation method that separates molecules that carry uracils from those that do not for a set of samples that includes Neandertal remains, herbarium specimens and archaeological plant remains. Results We show that sequencing DNA libraries enriched in molecules carrying uracils effectively amplifies age associated degradation patterns in microbial mixtures of ancient and historical origin. This facilitates the discovery of authentic ancient microbial taxa in cases where degradation patterns are difficult to detect due to large sequence divergence in microbial mixtures. Additionally, the relative enrichment of taxa in the uracil enriched fraction can help to identify bona fide ancient microbial taxa that could be missed using a more targeted approach. Conclusions Our experiments show, that in addition to its use in enriching authentic endogenous DNA of organisms of interest, the selective enrichment of damaged DNA molecules can be a valuable tool in the discovery of ancient microbial taxa.

SARS-CoV-2 causes substantial morbidity and mortality in elderly and immunocompromised individuals, particularly in retirement homes, where transmission from asymptomatic staff and visitors may introduce the infection. Here we present a cheap and fast screening method based on direct RT-qPCR to detect SARS-CoV-2 in single or pooled gargle lavages (“mouthwashes”). This method detects individuals with large viral loads (Ct≤29) and we use it to test all staff at a nursing home daily over a period of three weeks in order to reduce the risk that the infection penetrates the facility. This or similar approaches can be implemented to protect hospitals, nursing homes and other institutions in this and future viral epidemics.

Wallace's Line demarcates a biogeographical boundary between the Indomalaya and Australasian ecoregions. Most placental mammalian genera, for example, occur to the west of this line, whereas most marsupial genera occur to the east. However, macaque monkeys are unusual because they naturally occur on both western and eastern sides. To further explore this anomalous distribution, we analyzed 222 mitochondrial genomes from ∼20 macaque species, including new genomes from 60 specimens. These comprise a population sampling of most Sulawesi macaques, Macaca fascicularis (long-tailed macaques) specimens that were collected by Alfred R. Wallace and specimens that were recovered during archaeological excavations at Liang Bua, a cave on the Indonesian island of Flores. In M. fascicularis, three mitochondrial lineages span the southernmost portion of Wallace's Line between Bali and Lombok, and divergences within these lineages are contemporaneous with, and possibly mediated by, past dispersals of modern human populations. Near the central portion of Wallace's Line between Borneo and Sulawesi, a more ancient dispersal of macaques from mainland Asia to Sulawesi preceded modern human colonization, which was followed by rapid dispersal of matrilines and was subsequently influenced by recent interspecies hybridization. In contrast to previous studies, we find no strong signal of recombination in most macaque mitochondrial genomes. These findings further characterize macaque evolution before and after modern human dispersal throughout Southeast Asia and point to possible effects on biodiversity of ancient human cultural diasporas.

Additional file 3: Figure S3. Taxonomic tree of reads from a Solanum tuberosum library assigned to different taxonomic levels. The size of the circle represents the amount of reads assigned to the node displayed in the tree or to any taxonomic level below it. Reads assigned to some species Phytophthora infestans, Pseudomonas syringae and Pseudomonas rhizosphaerae, as well as the genera Pseudomonas and Solanum are named in the taxonomic tree.

Additional file 2: Figure S2. Length and GC content of plant historical specimens. A. Distributions of mean length for non-selected and U-selected libraries (U-enriched and U-depleted fractions). Median values are denoted as black lines and points show the original value for each individual sample. B. Length distribution of Arabidopsis thaliana sample NY1365375 for a non-selected and U-selected library (U-enriched and U-depleted fractions). C. Distributions of mean GC content for non-selected and U-selected libraries (U-enriched and U-depleted fractions).

Additional file 1: Figure S1. Mapped and taxonomically assigned reads of plant historical specimens. A. Distributions of percentage of mapped reads for non-selected and U-selected libraries (U-enriched and U-depleted fractions). B. Distributions of percentage of taxonomically assigned reads for non-selected and U-selected libraries (U-enriched and U-depleted fractions). C. Correlation of the percentage of mapped reads between the U-enriched and the non-selected library D. Correlation of the percentage of taxonomically assigned reads between the U-enriched and the non-selected library E. Relation between percentages of mapped and taxonomically assigned reads from U-selected libraries (U-enriched fraction).

Additional file 4: Figure S4. Analysis of de novo assembled Pseudomonas contigs from a Solanum tuberosum sample. Contig coverage of Pseudomonas syringae (A.) and Pseudomonas rhizosphaerae (B.). The circles represent each circular bacterial reference genome, and the black lines depict genomic regions covered by alignments of de novo assembled contigs. C. Percentage of reference genome of P. syringae and P. rhizosphaerae covered by de novo assembled contigs from A. and B., respectively.

Additional file 5: Figure S5. Analysis of uniquely mapped Pseudomonas contigs from a Solanum tuberosum sample. Uniquely mapped contig coverage of Pseudomonas syringae (A.) and Pseudomonas rhizosphaerae (B.). The circles represent each circular bacterial reference genome, and the black lines depict genomic regions covered by alignments of uniquely mapped contigs. Histograms of contig k-mer coverage from de novo assembled contigs uniquely mapping to P. syringae (C.) and P. rhizosphaerae (D.).

It has been shown that highly fragmented DNA is most efficiently converted into DNA libraries for sequencing if both strands of the DNA fragments are processed independently. We present an updated protocol for library preparation from single-stranded DNA, which is based on the splinted ligation of an adapter oligonucleotide to the 3′ ends of single DNA strands, the synthesis of a complementary strand using a DNA polymerase and the addition of a 5′ adapter via blunt-end ligation. The efficiency of library preparation is determined individually for each sample using a spike-in oligonucleotide. The whole workflow, including library preparation, quantification and amplification, requires two work days for up to 16 libraries. Alternatively, we provide documentation and electronic protocols enabling automated library preparation of 96 samples in parallel on a Bravo NGS Workstation (Agilent Technologies). After library preparation, molecules with uninformative short inserts (shorter than ~30−35 base pairs) can be removed by polyacrylamide gel electrophoresis if desired. Here the authors describe an updated protocol for single-stranded sequencing library preparation suitable for highly degraded DNA from ancient remains or other sources. The procedure can be performed manually or in an automated fashion.

Phenotypic invariance—the outcome of purifying selection—is a hallmark of biological importance. However, invariant phenotypes might be controlled by diverged genetic systems in different species. Here, we explore how an important and invariant phenotype—the development of sexually differentiated individuals—is controlled in over two dozen species in the frog family Pipidae. We uncovered evidence in different species for 1) an ancestral W chromosome that is not found in many females and is found in some males, 2) independent losses and 3) autosomal segregation of this W chromosome, 4) changes in male versus female heterogamy, and 5) substantial variation among species in recombination suppression on sex chromosomes. We further provide evidence of, and evolutionary context for, the origins of at least seven distinct systems for regulating sex determination among three closely related genera. These systems are distinct in their genomic locations, evolutionary origins, and/or male versus female heterogamy. Our findings demonstrate that the developmental control of sexual differentiation changed via loss, sidelining, and empowerment of a mechanistically influential gene, and offer insights into novel factors that impinge on the diverse evolutionary fates of sex chromosomes.

Genomic library preparation from highly degraded DNA is more efficient when library molecules are prepared separately from the complementary strands of DNA fragments. We describe a protocol in which libraries are constructed from single DNA strands in a three-step procedure: single-stranded ligation of the first adapter with T4 DNA ligase in the presence of a splinter oligonucleotide, copying of the DNA strand with a proofreading polymerase, and blunt-end ligation of the second double-stranded adapter with T4 DNA ligase.

Genomic library preparation from highly degraded DNA is more efficient when library molecules are prepared separately from the complementary strands of DNA fragments. We describe a protocol in which libraries are constructed from single DNA strands in a three-step procedure: single-stranded ligation of the first adapter with T4 DNA ligase in the presence of a splinter oligonucleotide, copying of the DNA strand with a proofreading polymerase, and blunt-end ligation of the second double-stranded adapter with T4 DNA ligase.

A comprehensive, accurate, and revisable alpha taxonomy is crucial for biodiversity studies, but is challenging when data from reference specimens are difficult to collect or observe. However, recent technological advances can overcome some of these challenges. To illustrate this, we used modern approaches to tackle a centuries-old taxonomic enigma presented by Fraser's Clawed Frog, Xenopus fraseri, including whether X. fraseri is different from other species, and if so, where it is situated geographically and phylogenetically. To facilitate these inferences, we used high-resolution techniques to examine morphological variation, and we generated and analyzed complete mitochondrial genome sequences from all Xenopus species, including >150-year-old type specimens. Our results demonstrate that X. fraseri is indeed distinct from other species, firmly place this species within a phylogenetic context, and identify its minimal geographic distribution in northern Ghana and northern Cameroon. These data also permit novel phylogenetic resolution into this intensively studied and biomedically important group. Xenopus fraseri was formerly thought to be a rainforest endemic placed alongside species in the amieti species group; in fact this species occurs in arid habitat on the borderlands of the Sahel, and is the smallest member of the muelleri species group. This study illustrates that the taxonomic enigma of Fraser's frog was a combined consequence of sparse collection records, interspecies conservation and intraspecific polymorphism in external anatomy, and type specimens with unusual morphology.

Denisovans, a sister group of Neandertals, have been described on the basis of a nuclear genome sequence from a finger phalanx (Denisova 3) found in Denisova Cave in the Altai Mountains. The only other Denisovan specimen described to date is a molar (Denisova 4) found at the same site. This tooth carries a mtDNA sequence similar to that of Denisova 3. Here we present nuclear DNA sequences from Denisova 4 and a morphological description, as well as mitochondrial and nuclear DNA sequence data, from another molar (Denisova 8) found in Denisova Cave in 2010. This new molar is similar to Denisova 4 in being very large and lacking traits typical of Neandertals and modern humans. Nuclear DNA sequences from the two molars form a clade with Denisova 3. The mtDNA of Denisova 8 is more diverged and has accumulated fewer substitutions than the mtDNAs of the other two specimens, suggesting Denisovans were present in the region over an extended period. The nuclear DNA sequence diversity among the three Denisovans is comparable to that among six Neandertals, but lower than that among present-day humans.

Although great progress has been made in improving methods for generating DNA sequences from ancient biological samples, many, if not most, samples are still not amenable for analyses due to overwhelming contamination with microbial or modern human DNA. Here we explore different DNA decontamination procedures for ancient bones and teeth for use prior to DNA library preparation and high-throughput sequencing. Two procedures showed promising results: (i) the release of surface-bound DNA by phosphate buffer and (ii) the removal of DNA contamination by sodium hypochlorite treatment. Exposure to phosphate removes on average 64% of the microbial DNA from bone powder but only 37% of the endogenous DNA (from the organism under study), increasing the percentage of informative sequences by a factor of two on average. An average 4.6-fold increase, in one case reaching 24-fold, is achieved by sodium hypochlorite treatment, albeit at the expense of destroying 63% of the endogenous DNA preserved in the bone. While both pretreatment methods described here greatly reduce the cost of genome sequencing from ancient material due to efficient depletion of microbial DNA, we find that the removal of human DNA contamination remains a challenging problem.

Contamination by present-day human and microbial DNA is one of the major hindrances for large-scale genomic studies using ancient biological material. We describe a new molecular method, U selection, which exploits one of the most distinctive features of ancient DNA—the presence of deoxyuracils—for selective enrichment of endogenous DNA against a complex background of contamination during DNA library preparation. By applying the method to Neanderthal DNA extracts that are heavily contaminated with present-day human DNA, we show that the fraction of useful sequence information increases ~10-fold and that the resulting sequences are more efficiently depleted of human contamination than when using purely computational approaches. Furthermore, we show that U selection can lead to a four- to fivefold increase in the proportion of endogenous DNA sequences relative to those of microbial contaminants in some samples. U selection may thus help to lower the costs for ancient genome sequencing of nonhuman samples also. [Supplemental material is available for this article.]

Although an inverse relationship is expected in ancient DNA samples between the number of surviving DNA fragments and their length, ancient DNA sequencing libraries are strikingly deficient in molecules shorter than 40 bp. We find that a loss of short molecules can occur during DNA extraction and present an improved silica-based extraction protocol that enables their efficient retrieval. In combination with single-stranded DNA library preparation, this method enabled us to reconstruct the mitochondrial genome sequence from a Middle Pleistocene cave bear (Ursus deningeri) bone excavated at Sima de los Huesos in the Sierra de Atapuerca, Spain. Phylogenetic reconstructions indicate that the U. deningeri sequence forms an early diverging sister lineage to all Western European Late Pleistocene cave bears. Our results prove that authentic ancient DNA can be preserved for hundreds of thousand years outside of permafrost. Moreover, the techniques presented enable the retrieval of phylogenetically informative sequences from samples in which virtually all DNA is diminished to fragments shorter than 50 bp.

DNA library preparation for high-throughput sequencing of genomic DNA usually involves ligation of adapters to double-stranded DNA fragments. However, for highly degraded DNA, especially ancient DNA, library preparation has been found to be more efficient if each of the two DNA strands are converted into library molecules separately. We present a new method for single-stranded library preparation, ssDNA2.0, which is based on single-stranded DNA ligation with T4 DNA ligase utilizing a splinter oligonucleotide with a stretch of random bases hybridized to a 3΄ biotinylated donor oligonucleotide. A thorough evaluation of this ligation scheme shows that single-stranded DNA can be ligated to adapter oligonucleotides in higher concentration than with CircLigase (an RNA ligase that was previously chosen for end-to-end ligation in single-stranded library preparation) and that biases in ligation can be minimized when choosing splinters with 7 or 8 random nucleotides. We show that ssDNA2.0 tolerates higher quantities of input DNA than CircLigase-based library preparation, is less costly and better compatible with automation. We also provide an in-depth comparison of library preparation methods on degraded DNA from various sources. Most strikingly, we find that single-stranded library preparation increases library yields from tissues stored in formalin for many years by several orders of magnitude.

Of all obligate intracellular bacteria, Wolbachia is probably the most common. In general, Wolbachia are either widespread, opportunistic reproductive parasites of arthropods or essential mutualists in a single group of filarial nematodes, including many species of medical significance. To date, a robust phylogenetic backbone of Wolbachia is lacking and consequently, many Wolbachia-related phenomena cannot be discussed in a broader evolutionary context. Here we present the first comprehensive phylogenomic analysis of Wolbachia supergroup relationships based on new whole-genome-shotgun data. Our results suggest that Wolbachia has switched between its two major host groups at least twice. The ability of some arthropod-infecting Wolbachia to universally infect and to adapt to a broad range of hosts quickly is restricted to a single monophyletic lineage (containing supergroups A and B). Thus, the currently observable pandemic has likely a single evolutionary origin and is unique within the radiation of Wolbachia strains.

This protocol describes a method for converting short single-stranded and double-stranded DNA into libraries compatible with high-throughput sequencing using Illumina technology. This method has primarily been developed to improve sequence retrieval from ancient DNA, but it is also applicable to the sequencing of short or degraded DNA from other sources, and it can also be used for sequencing oligonucleotides. Single-stranded library preparation is performed by ligating a biotinylated adapter oligonucleotide to the 3′ ends of heat-denatured DNA. The resulting strands are then immobilized on streptavidin-coated beads and copied with a polymerase. A second adapter is attached by blunt-end ligation, and library preparation is completed by PCR amplification. We estimate that intact DNA strands are recovered in the library with ∼50% efficiency. Libraries can be generated from up to 12 DNA or oligonucleotide samples in parallel within 2 d.

Ancient Genomics The Denisovans were archaic humans closely related to Neandertals, whose populations overlapped with the ancestors of modern-day humans. Using a single-stranded library preparation method, Meyer et al. (p. 222 , published online 30 August) provide a detailed analysis of a high-quality Denisovan genome. The genomic sequence provides evidence for very low rates of heterozygosity in the Denisova, probably not because of recent inbreeding, but instead because of a small population size. The genome sequence also illuminates the relationships between humans and archaics, including Neandertals, and establishes a catalog of genetic changes within the human lineage.

DNA retrieved from a tooth discovered deep in Denisova Cave allows us to assign it to the Denisovans, a group of archaic hominins., The presence of Neandertals in Europe and Western Eurasia before the arrival of anatomically modern humans is well supported by archaeological and paleontological data. In contrast, fossil evidence for Denisovans, a sister group of Neandertals recently identified on the basis of DNA sequences, is limited to three specimens, all of which originate from Denisova Cave in the Altai Mountains (Siberia, Russia). We report the retrieval of DNA from a deciduous lower second molar (Denisova 2), discovered in a deep stratigraphic layer in Denisova Cave, and show that this tooth comes from a female Denisovan individual. On the basis of the number of “missing substitutions” in the mitochondrial DNA determined from the specimen, we find that Denisova 2 is substantially older than two of the other Denisovans, reinforcing the view that Denisovans were likely to have been present in the vicinity of Denisova Cave over an extended time period. We show that the level of nuclear DNA sequence diversity found among Denisovans is within the lower range of that of present-day human populations.