Functional dissection and assembly of a small, newly evolved, W chromosome-specific genomic region of the African clawed frog Xenopus laevis.

Genetic triggers for sex determination are frequently co-inherited with other linked genes that may also influence one or more sex-specific phenotypes. To better understand how sex-limited regions evolve and function, we studied a small W chromosome-specific region of the frog Xenopus laevis that contains only three genes (dm-w, scan-w, ccdc69-w) and that drives female differentiation. Using gene editing, we found that the sex-determining function of this region requires dm-w but that scan-w and ccdc69-w are not essential for viability, female development, or fertility. Analysis of mesonephros+gonad transcriptomes during sexual differentiation illustrates masculinization of the dm-w knockout transcriptome, and identifies mostly non-overlapping sets of differentially expressed genes in separate knockout lines for each of these three W-specific gene compared to wildtype sisters. Capture sequencing of almost all Xenopus species and PCR surveys indicate that the female-determining function of dm-w is present in only a subset of species that carry this gene. These findings map out a dynamic evolutionary history of a newly evolved W chromosome-specific genomic region, whose components have distinctive functions that frequently degraded during Xenopus diversification, and evidence the evolutionary consequences of recombination suppression.

Mitogenomics of macaques (Macaca) across Wallace's Line in the context of modern human dispersals.

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.

Manual and automated preparation of single-stranded DNA libraries for the sequencing of DNA from ancient biological remains and other sources of highly degraded DNA.

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.

Mining ancient microbiomes using selective enrichment of damaged DNA molecules.

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.

Developmental Systems Drift and the Drivers of Sex Chromosome Evolution.

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.

Xenopus fraseri: Mr. Fraser, where did your frog come from?

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.

A Method for Single-Stranded Ancient DNA Library Preparation.

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 fourth Denisovan individual.

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.

Single-stranded DNA library preparation from highly degraded DNA using T4 DNA ligase.

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.

Nuclear and mitochondrial DNA sequences from two Denisovan individuals.

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.

Reducing microbial and human contamination in DNA extractions from ancient bones and teeth.

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.

Phylogenomic analyses uncover origin and spread of the Wolbachia pandemic.

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.

Selective enrichment of damaged DNA molecules for ancient genome sequencing.

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.

Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments.

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.

Single-stranded DNA library preparation for the sequencing of ancient or damaged DNA.

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.

A high-coverage genome sequence from an archaic Denisovan individual.

We present a DNA library preparation method that has allowed us to reconstruct a high-coverage (30×) genome sequence of a Denisovan, an extinct relative of Neandertals. The quality of this genome allows a direct estimation of Denisovan heterozygosity indicating that genetic diversity in these archaic hominins was extremely low. It also allows tentative dating of the specimen on the basis of "missing evolution" in its genome, detailed measurements of Denisovan and Neandertal admixture into present-day human populations, and the generation of a near-complete catalog of genetic changes that swept to high frequency in modern humans since their divergence from Denisovans.