Analysis of the Batch Effect Due to Sequencing Center in Population Statistics Quantifying Rare Events in the 1000 Genomes Project

The 1000 Genomes Project (1000G) is one of the most popular whole genome sequencing datasets used in different genomics fields and has boosting our knowledge in medical and population genomics, among other fields. Recent studies have reported the presence of ghost mutation signals in the 1000G. Furthermore, studies have shown that these mutations can influence the outcomes of follow-up studies based on the genetic variation of 1000G, such as single nucleotide variants (SNV) imputation. While the overall effect of these ghost mutations can be considered negligible for common genetic variants in many populations, the potential bias remains unclear when studying low frequency genetic variants in the population. In this study, we analyze the effect of the sequencing center in predicted loss of function (LoF) alleles, the number of singletons, and the patterns of archaic introgression in the 1000G. Our results support previous studies showing that the sequencing center is associated with LoF and singletons independent of the population that is considered. Furthermore, we observed that patterns of archaic introgression were distorted for some populations depending on the sequencing center. When analyzing the frequency of SNPs showing extreme patterns of genotype differentiation among centers for CEU, YRI, CHB, and JPT, we observed that the magnitude of the sequencing batch effect was stronger at MAF < 0.2 and showed different profiles between CHB and the other populations. All these results suggest that data from 1000G must be interpreted with caution when considering statistics using variants at low frequency.

Archaic introgression and variation in pharmacogenes and implications for local adaptation

Modern humans carry Neanderthal and Denisovan (archaic) genome elements which may have been a result of environmental adaptation. These effects may be particularly evident in pharmacogenes - genes responsible for the processing of exogenous substances such as food, pollutants, and medications. However, the health implications and contribution of archaic ancestry in pharmacogenes of modern humans remains understudied. We characterize eleven key cytochrome P450 (CYP450) genes involved in drug metabolizing reactions in three Neanderthal and one Denisovan individuals and examine archaic introgression in modern human populations. We infer the metabolizing efficiency of these eleven genes in archaic individuals and show important genetic differences relative to modern human variants. We identify archaic-specific SNVs in each CYP450 gene, including some that are potentially damaging, which may result in altered metabolism in modern human people carrying these variants. We highlight four genes which display interesting patterns of archaic variation: CYP2B6 - we find a large number of unique variants in the Vindija Neanderthal, some of which are shared with a small subset of African modern humans; CYP2C9 - containing multiple variants that are shared between Europeans and Neanderthals; CYP2A6*12 - a variant defined by a hybridization event that was found in humans and Neanderthals, suggesting the recombination event predates both species; and CYP2J2 - in which we hypothesize a Neanderthal variant was re-introduced in non-African populations by archaic admixture. The genetic variation identified in archaic individuals imply environmental pressures that may have driven CYP450 gene evolution.

Local adaptation and archaic introgression shape global diversity at human structural variant loci

Large genomic insertions, deletions, and inversions are a potent source of functional and fitness-altering variation, but are challenging to resolve with short-read DNA sequencing alone. While recent long-read sequencing technologies have greatly expanded the catalog of structural variants (SVs), their costs have so far precluded their application at population scales. Given these limitations, the role of SVs in human adaptation remains poorly characterized. Here, we used a graph-based approach to genotype 107,866 long-read-discovered SVs in short-read sequencing data from diverse human populations. We then applied an admixture-aware method to scan these SVs for patterns of population-specific frequency differentiation—a signature of local adaptation. We identified 220 SVs exhibiting extreme frequency differentiation, including several SVs that were among the lead variants at their corresponding loci. The top two signatures traced to separate insertion and deletion polymorphisms at the immunoglobulin heavy chain locus, together tagging a 325 Kbp haplotype that swept to high frequency and was subsequently fragmented by recombination. Alleles defining this haplotype are nearly fixed (60-95%) in certain Southeast Asian populations, but are rare or absent from other global populations composing the 1000 Genomes Project. Further investigation revealed that the haplotype closely matches with sequences observed in two of three high-coverage Neanderthal genomes, providing strong evidence of a Neanderthal-introgressed origin. This extraordinary episode of positive selection, which we infer to have occurred between 1700 and 8400 years ago, corroborates the role of immune-related genes as prominent targets of adaptive archaic introgression. Our study demonstrates how combining recent advances in genome sequencing, genotyping algorithms, and population genetic methods can reveal signatures of key evolutionary events that remained hidden within poorly resolved regions of the genome.

Decomposing the admixture statistic, D, suggests a negligible contribution due to archaic introgression into humans

It is widely accepted that non-African humans carry a few percent of Neanderthal DNA due to historical inter-breeding. However, methods used to infer a legacy all assume that mutation rate is constant and that back-mutations can be ignored. Here I decompose the widely used admixture statistic, D, in a way that allows the overall signal to be apportioned to different classes of contributing site. I explore three main characteristics: whether the putative Neanderthal allele is likely derived or ancestral; whether an allele is fixed in one of the two human populations; and the type of mutation that created the polymorphism, defined by the base that mutated and immediately flanking bases. The entire signal used to infer introgression can be attributed to a subset of sites where the putative Neanderthal base is common in Africans and fixed in non-Africans. Moreover, the four triplets containing highly mutable CpG motifs alone contribute 29%. In contrast, sites expected to dominate the signal if introgression has occurred, where the putative Neanderthal allele is absent from Africa and rare outside Africa, contribute negligibly. Together, these observations show that D does not capture a signal due to introgression but instead they support an alternative model in which a higher mutation rate in Africa drives increased divergence from the ancestral state.

A test of the hypothesis that variable mutation rates create signals that have previously been interpreted as evidence of archaic introgression into humans

It is widely accepted that non-African humans carry 1-2% Neanderthal DNA due to historical inter-breeding. However, inferences about introgression rely on a critical assumption that mutation rate is constant and that back-mutations are too rare to be important. Both these assumptions have been challenged, and recent evidence points towards an alternative model where signals interpreted as introgression are driven mainly by higher mutation rates in Africa. In this model, non-Africans appear closer to archaics not because they harbour introgressed fragments but because Africans have diverged more. Here I test this idea by using the density of rare, human-specific variants (RHSVs) as a proxy for recent mutation rate. I find that sites that contribute most to the signal interpreted as introgression tend to occur in tightly defined regions spanning only a few hundred bases in which mutation rate differs greatly between the two human populations being compared. Mutation rate is invariably higher in the population into which introgression is not inferred. I confirmed that RHSV density reflects mutation rate by conducting a parallel analysis looking at the density of RHSVs around sites with three alleles, an independent class of site that also requires recurrent mutations to form. Near-identical peaks in RHSV density are found, suggesting a common cause. Similarly, coalescent simulations confirm that, with constant mutation rate, introgressed fragments do not occur preferentially in regions with a high density of rare, human-specific variants. Together, these observations are difficult to reconcile with a model where excess base-sharing is driven by archaic legacies but instead provide support for a higher mutation rate inside Africa driving increased divergence from the ancestral human state.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with local ancestry and archaic introgression

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity on the region. In a comparison with a European dataset, we identify 166 uniquely Indonesia-specific eQTLs, highlighting the benefits of performing association studies on non-European populations. By combining local ancestry and archaic introgression inference eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits. Our findings illustrate how local ancestry and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations.

Introgression, hominin dispersal and megafaunal survival in Late Pleistocene Island Southeast Asia

The hominin fossil record of Island Southeast Asia (ISEA) indicates that at least two endemic ‘super-archaic’ species – Homo luzonensis and H. floresiensis – were present around the time anatomically modern humans (AMH) arrived in the region >50,000 years ago. Contemporary human populations carry signals consistent with interbreeding events with Denisovans in ISEA – a species that is thought to be more closely related to AMH than the super-archaic endemic ISEA hominins. To query this disparity between fossil and genetic evidence, we performed a comprehensive search for super-archaic introgression in >400 modern human genomes. Our results corroborate widespread Denisovan ancestry in ISEA populations but fail to detect any super-archaic admixture signals. By highlighting local megafaunal survival east of the Wallace Line as a potential signature of deep, pre-H. sapiens hominin-faunal interaction, we propose that this understudied region may hold the key to unlocking significant chapters in Denisovan prehistory.

Signals interpreted as archaic introgression appear to be driven primarily by faster evolution in Africa

Non-African humans appear to carry a few per cent archaic DNA due to ancient inter-breeding. This modest legacy and its likely recent timing imply that most introgressed fragments will be rare and hence will occur mainly in the heterozygous state. I tested this prediction by calculating D statistics, a measure of legacy size, for pairs of humans where one of the pair was conditioned always to be either homozygous or heterozygous. Using coalescent simulations, I confirmed that conditioning the non-African to be heterozygous increased D, while conditioning the non-African to be homozygous reduced D to zero. Repeating with real data reveals the exact opposite pattern. In African–non-African comparisons, D is near-zero if the African individual is held homozygous. Conditioning one of two Africans to be either homozygous or heterozygous invariably generates large values of D, even when both individuals are drawn from the same population. Invariably, the African with more heterozygous sites (conditioned heterozygous > unconditioned > conditioned homozygous) appears less related to the archaic. By contrast, the same analysis applied to pairs of non-Africans always yields near-zero D, showing that conditioning does not create large D without an underlying signal to expose. Large D values in humans are therefore driven almost entirely by heterozygous sites in Africans acting to increase divergence from related taxa such as Neanderthals. In comparison with heterozygous Africans, individuals that lack African heterozygous sites, whether non-African or conditioned homozygous African, always appear more similar to archaic outgroups, a signal previously interpreted as evidence for introgression. I hope these analyses will encourage others to consider increased divergence as well as increased similarity to archaics as mechanisms capable of driving asymmetrical base-sharing.

Analysis of Haplotypic Variation and Deletion Polymorphisms Point to Multiple Archaic Introgression Events, Including from Altai Neanderthal Lineage

The time, extent, and genomic effect of the introgressions from archaic humans into ancestors of extant human populations remain some of the most exciting venues of population genetics research in the past decade. Several studies have shown population-specific signatures of introgression events from Neanderthals, Denisovans, and potentially other unknown hominin populations in different human groups. Moreover, it was shown that these introgression events may have contributed to phenotypic variation in extant humans, with biomedical and evolutionary consequences. In this study, we present a comprehensive analysis of the unusually divergent haplotypes in the Eurasian genomes and show that they can be traced back to multiple introgression events. In parallel, we document hundreds of deletion polymorphisms shared with Neanderthals. A locus-specific analysis of one such shared deletion suggests the existence of a direct introgression event from the Altai Neanderthal lineage into the ancestors of extant East Asian populations. Overall, our study is in agreement with the emergent notion that various Neanderthal populations contributed to extant human genetic variation in a population-specific manner.