In-solution Y-chromosome capture-enrichment on ancient DNA libraries

de Pancorbo M, Calafell F. Y-chromosome target enrichment reveals rapid expansion of haplogroup R1b-DF27 in Iberia during the Bronze Age transition

The Y chromosome can yield a unique perspective into the study of human demographic history. However, due to the repetitive nature of part of its sequence, only a small set of regions are suitable for variant calling and discovery from short-read sequencing data. These regions combined represent 8.9 Mbp or 0.14% of a diploid human genome. Consequently, investing in whole-genome sequencing to resolve Y-chromosome questions is poorly efficient. Here we use, as an alternative, target enrichment technology to greatly increase sequencing effectiveness, validating and applying the technique to 181 males, for 162 of whom we obtained a positive result. Additionally, 75 samples sequenced for the whole genome were also included, for a total sample size of 237. These samples were chosen for their Y chromosome haplogroup: R1b-DF27. In the context of European populations, and particularly in Iberia, this haplogroup stands out for its high frequency and its demographic history. Current evidence indicates that the diffusion of this haplogroup is related to the population movements that mark the cultural Bronze Age transition, making it remarkably interesting for population geneticists. The results of this study show the effects of the rapid radiation of the haplogroup in Spain, as even with the higher discriminating power of whole sequences, most haplotypes still fall within the R1b-DF27* paragroup rather than in the main derived branches. However, we were able to refine the ISOGG 2019-2020 phylogeny, and its two main subbranches, namely L176.2 and Z272, which present geographical differentiation between the Atlantic and Mediterranean coasts of Iberia.

The hazards of genotype imputation in chromosomal regions under selection: A case study using the Lactase gene region

Although imputation of missing SNP results has been widely used in genetic studies, claims about the quality and usefulness of imputation have outnumbered the few studies that have questioned its limitations. But it is becoming clear that these limitations are real-for example, disease association signals can be missed in regions of LD breakdown. Here, as a case study, using the chromosomal region of the well-known lactase gene, LCT, we address the issue of imputation in the context of variants that have become frequent in a limited number of modern population groups only recently, due to selection. We study SNPs in a 500 bp region covering the enhancer of LCT, and compare imputed genotypes with directly genotyped data. We examine the haplotype pairs of all individuals with discrepant and missing genotypes. We highlight the nonrandom nature of the allelic errors and show that most incorrect imputations and missing data result from long haplotypes that are evolutionarily closely related to those carrying the derived alleles, while some relate to rare and recombinant haplotypes. We conclude that bias of incorrectly imputed and missing genotypes can decrease the accuracy of imputed results substantially.

Using Y-chromosome capture enrichment to resolve haplogroup H2 shows new evidence for a two-path Neolithic expansion to Western Europe

Uniparentally-inherited markers on mitochondrial DNA (mtDNA) and the non-recombining regions of the Y chromosome (NRY), have been used for the past 30 years to investigate the history of humans from a maternal and paternal perspective. Researchers have preferred mtDNA due to its abundance in the cells, and comparatively high substitution rate. Conversely, the NRY is less susceptible to back mutations and saturation, and is potentially more informative than mtDNA owing to its longer sequence length. However, due to comparatively poor NRY coverage via shotgun sequencing, and the relatively low and biased representation of Y-chromosome variants on capture assays such as the 1240 k, ancient DNA studies often fail to utilize the unique perspective that the NRY can yield. Here we introduce a new DNA enrichment assay, coined YMCA (Y-mappable capture assay), that targets the "mappable" regions of the NRY. We show that compared to low-coverage shotgun sequencing and 1240 k capture, YMCA significantly improves the mean coverage and number of sites covered on the NRY, increasing the number of Y-haplogroup informative SNPs, and allowing for the identification of previously undiscovered variants. To illustrate the power of YMCA, we show that the analysis of ancient Y-chromosome lineages can help to resolve Y-chromosomal haplogroups. As a case study, we focus on H2, a haplogroup associated with a critical event in European human history: the Neolithic transition. By disentangling the evolutionary history of this haplogroup, we further elucidate the two separate paths by which early farmers expanded from Anatolia and the Near East to western Europe.

Genetic characterization of a collection of Tsantsas from Ecuadorian museums

Tsantsas are shrunken human heads originally made for ceremonial purposes by Amazonian indigenous groups of the Shuar and Achuar family, previously called Jivaroan tribes. A significant demand of these objects during the first half of the 20th century led to the manufacture of counterfeit shrunken heads for commercial purposes. For museums where these collections are held, as well as for the indigenous groups who claim their ownership, it is important to identify the origin and authenticity of these tsantsas. We hypothesized that a collection of 14 tsantsas from 3 different museum collections in Ecuador are human and aimed to characterize their sex and potential origin. We amplified the amelogenin gene and performed a high resolution melting analysis to determine their human origin and characterize their sex. We also analyzed a fragment (16209-16402) from the HVR-1 region to identify the mtDNA haplogroups present in the tsantsa collection. Our exploratory results show that all the tsantsas are human and that the collection is comprised of 13 males and 1 female. A total of seven mtDNA haplogroups were found among the tsantsa collection using the mtDNA EMPOP database. These results show a predominance of the Amerindian mtDNA haplogroups B, C and D. Additional principal component analysis, genetic distance tree and haplotype network analyses suggest a relationship between the tsantsa specimens and Native American groups.

Ancient DNA analysis of rare genetic bone disorders

OBJECTIVE

Review of the current advancements in the field of paleogenetics that provide new opportunities in studying the evolution of rare genetic bone diseases.

MATERIAL AND METHODS

Based on cases from the literature, the genetics of rare bone diseases will be introduced and the main methodological issues will be addressed, focusing on the opportunities presented by the application of aDNA analyses in the field of paleopathology.

RESULTS

Medical literature provides large datasets on the genes responsible for rare bone disorders. These genes, subdivided in functional categories, display important future targets when analyzing rare genetic bone disorders in ancient human remains.

CONCLUSIONS

Knowledge on both phenotype and genotype is required to study rare diseases in ancient human remains.

SIGNIFICANCE

The proposed interdisciplinary research will provide new insight into the occurrence and spread of genetic risk factors in the past and will help in the diagnostics of these rare and often neglected diseases.

LIMITATIONS

The current limitations in ancient DNA research and targeting the disease-causing specific mutations (e.g., somatic or germline).

SUGGESTIONS FOR FURTHER RESEARCH

Methodological advancements and candidate gene lists provide the optimal basis for future interdisciplinary studies of rare genetic bone disorders in ancient human remains.

The evolutionary history of Neanderthal and Denisovan Y chromosomes

Ancient DNA has provided new insights into many aspects of human history. However, we lack comprehensive studies of the Y chromosomes of Denisovans and Neanderthals because the majority of specimens that have been sequenced to sufficient coverage are female. Sequencing Y chromosomes from two Denisovans and three Neanderthals shows that the Y chromosomes of Denisovans split around 700 thousand years ago from a lineage shared by Neanderthals and modern human Y chromosomes, which diverged from each other around 370 thousand years ago. The phylogenetic relationships of archaic and modern human Y chromosomes differ from the population relationships inferred from the autosomal genomes and mirror mitochondrial DNA phylogenies, indicating replacement of both the mitochondrial and Y chromosomal gene pools in late Neanderthals. This replacement is plausible if the low effective population size of Neanderthals resulted in an increased genetic load in Neanderthals relative to modern humans.

Paleogenomic insights into the red complex bacteria Tannerella forsythia in Pre-Hispanic and Colonial individuals from Mexico

The 'red complex' is an aggregate of three oral bacteria (Tannerella forsythia, Porphyromonas gingivalis and Treponema denticola) responsible for severe clinical manifestation of periodontal disease. Here, we report the first direct evidence of ancient T.forsythia DNA in dentin and dental calculus samples from archaeological skeletal remains that span from the Pre-Hispanic to the Colonial period in Mexico. We recovered twelve partial ancient T. forsythia genomes and observed a distinct phylogenetic placement of samples, suggesting that the strains present in Pre-Hispanic individuals likely arrived with the first human migrations to the Americas and that new strains were introduced with the arrival of European and African populations in the sixteenth century. We also identified instances of the differential presence of genes between periods in the T. forsythia ancient genomes, with certain genes present in Pre-Hispanic individuals and absent in Colonial individuals, and vice versa. This study highlights the potential for studying ancient T. forsythia genomes to unveil past social interactions through analysis of disease transmission. Our results illustrate the long-standing relationship between this oral pathogen and its human host, while also unveiling key evidence to understand its evolutionary history in Pre-Hispanic and Colonial Mexico. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Targeted analysis of polymorphic loci from low-coverage shotgun sequence data allows accurate genotyping of HLA genes in historical human populations

The highly polymorphic human leukocyte antigen (HLA) plays a crucial role in adaptive immunity and is associated with various complex diseases. Accurate analysis of HLA genes using ancient DNA (aDNA) data is crucial for understanding their role in human adaptation to pathogens. Here, we describe the TARGT pipeline for targeted analysis of polymorphic loci from low-coverage shotgun sequence data. The pipeline was successfully applied to medieval aDNA samples and validated using both simulated aDNA and modern empirical sequence data from the 1000 Genomes Project. Thus the TARGT pipeline enables accurate analysis of HLA polymorphisms in historical (and modern) human populations.

A roadmap for high-throughput sequencing studies of wild animal populations using noninvasive samples and hybridization capture

Large-scale genomic studies of wild animal populations are often limited by access to high-quality DNA. Although noninvasive samples, such as faeces, can be readily collected, DNA from the sample producers is usually present in low quantities, fragmented, and contaminated by microorganism and dietary DNAs. Hybridization capture can help to overcome these impediments by increasing the proportion of subject DNA prior to high-throughput sequencing. Here we evaluate a key design variable for hybridization capture, the number of rounds of capture, by testing whether one or two rounds are most appropriate, given varying sample quality (as measured by the ratios of subject to total DNA). We used a set of 1,780 quality-assessed wild chimpanzee (Pan troglodytes schweinfurthii) faecal samples and chose 110 samples of varying quality for exome capture and sequencing. We used multiple regression to assess the effects of the ratio of subject to total DNA (sample quality), rounds of capture and sequencing effort on the number of unique exome reads sequenced. We not only show that one round of capture is preferable when the proportion of subject DNA in a sample is above ~2%-3%, but also explore various types of bias introduced by capture, and develop a model that predicts the sequencing effort necessary for a desired data yield from samples of a given quality. Thus, our results provide a useful guide and pave a methodological way forward for researchers wishing to plan similar hybridization capture studies.

The efficacy of whole human genome capture on ancient dental calculus and dentin

OBJECTIVES

Dental calculus is among the richest known sources of ancient DNA in the archaeological record. Although most DNA within calculus is microbial, it has been shown to contain sufficient human DNA for the targeted retrieval of whole mitochondrial genomes. Here, we explore whether calculus is also a viable substrate for whole human genome recovery using targeted enrichment techniques.

MATERIALS AND METHODS

Total DNA extracted from 24 paired archaeological human dentin and calculus samples was subjected to whole human genome enrichment using in-solution hybridization capture and high-throughput sequencing.

RESULTS

Total DNA from calculus exceeded that of dentin in all cases, and although the proportion of human DNA was generally lower in calculus, the absolute human DNA content of calculus and dentin was not significantly different. Whole genome enrichment resulted in up to four-fold enrichment of the human endogenous DNA content for both dentin and dental calculus libraries, albeit with some loss in complexity. Recovering more on-target reads for the same sequencing effort generally improved the quality of downstream analyses, such as sex and ancestry estimation. For nonhuman DNA, comparison of phylum-level microbial community structure revealed few differences between precapture and postcapture libraries, indicating that off-target sequences in human genome-enriched calculus libraries may still be useful for oral microbiome reconstruction.

DISCUSSION

While ancient human dental calculus does contain endogenous human DNA sequences, their relative proportion is low when compared with other skeletal tissues. Whole genome enrichment can help increase the proportion of recovered human reads, but in this instance enrichment efficiency was relatively low when compared with other forms of capture. We conclude that further optimization is necessary before the method can be routinely applied to archaeological samples.

Hybrid Capture-Based Next Generation Sequencing and Its Application to Human Infectious Diseases

This review describes target-enrichment approaches followed by next generation sequencing and their recent application to the research and diagnostic field of modern and past infectious human diseases caused by viruses, bacteria, parasites and fungi.