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

Abstract
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Large genomic insertions and deletions are a potent source of functional variation, but are challenging to resolve with short-read sequencing, limiting knowledge of the role of such structural variants (SVs) in human evolution. Here, we used a graph-based method to genotype long-read-discovered SVs in short-read data from diverse human genomes. We then applied an admixture-aware method to identify 220 SVs exhibiting extreme patterns of frequency differentiation – a signature of local adaptation. The top two variants traced to the immunoglobulin heavy chain locus, tagging a haplotype that swept to near fixation in certain southeast Asian populations, but is rare in other global populations. Further investigation revealed evidence that the haplotype traces to gene flow from Neanderthals, corroborating the role of immune-related genes as prominent targets of adaptive introgression. Our study demonstrates how recent technical advances can help resolve signatures of key evolutionary events that remained obscured within technically challenging regions of the genome.

Introduction

Rapid global dispersal and cultural evolution have exposed modern humans to a striking diversity of environments, to which they have developed numerous genetic adaptations (Fan et al., 2016). The vast majority of genomic research on human adaptation has focused on single-nucleotide polymorphisms (SNPs) and short insertions and deletions (indels) due to their ease of discovery by high-throughput short-read DNA sequencing (e.g., The 1000 Genomes Project Consortium, 2015; The International HapMap Consortium, 2007). This focus overlooks the potential impacts of larger insertions, deletions, and inversions – collectively termed structural variants (SVs) – whose abundance and complexity are now being appreciated with the advent of long-read DNA sequencing (Audano et al., 2019; Chaisson et al., 2015; Sedlazeck et al., 2018a). SVs impact more total sequence per genome than SNPs and may severely disrupt genes and regulatory elements to confer large effects on gene expression (Chiang et al., 2017). The size, modularity, and functional potential of SVs may offer a rich substrate for natural selection, as supported by several specific examples of adaptive insertions and deletions in diet and pigmentation-related genes (Hsieh et al., 2019; Kothapalli et al., 2016; Perry et al., 2007; Saitou and Gokcumen, 2019), as well as genome-wide evidence from primarily short-read data (Almarri et al., 2020; Sudmant et al., 2015).

More comprehensive analysis of SV evolution has been limited by the long and repetitive nature of many SVs. Short-read approaches to SV discovery depend on abnormalities in depth of coverage or other characteristics of read alignments (Kosugi et al., 2019) but achieve sensitivity of only 10–50% (Chaisson et al., 2019; Jakubosky et al., 2020b). In contrast, the recent development of long-read sequencing methods has revolutionized SV discovery, achieving high sensitivity and specificity by spanning SV sequences and their unique flanking regions (Sedlazeck et al., 2018a). This application of long-read sequencing to human samples has dramatically expanded the catalog of known human variation, generating databases of hundreds of thousands of SVs discovered in globally diverse individuals (Audano et al., 2019; Ebert et al., 2021). Yet with rare exceptions (Beyter et al., 2021), these sequencing methods remain impractical for population-scale applications due to their high costs and low throughput.

We addressed this limitation by applying a hybrid approach to integrate SV discovery from long-read sequencing and SV genotyping in short-read sequenced datasets. Using a graph-based method, we genotyped a catalog of long-read discovered SVs in high-coverage short-read sequencing data from the 1000 Genomes Project (1KGP) (Byrska-Bishop et al., 2021). By intersecting these genotype data with existing RNA-seq data from human cell lines, we discovered 1121 significant associations between SVs and the expression of nearby genes, including several components of the human leukocyte antigen (HLA) gene cluster and its interaction partners. We then applied an admixture-aware method to identify 220 SVs that exhibit strong allele frequency differentiation across human populations. Among the top hits, we discovered an extreme signature of local adaptation tagged by separate insertion and deletion polymorphisms at the immunoglobulin heavy chain locus, which encodes key components of the adaptive immune system. Alternative alleles defining this haplotype achieve high frequencies in certain southeast Asian populations, but are rare or absent from other global populations. Searching for signatures of archaic introgression within our set of highly differentiated SVs further revealed evidence that the adaptive haplotype entered the modern human population via ancient hybridization with Neanderthals. Together, our work highlights the role of structurally complex and repetitive regions of the genome as hidden sources of functional diversity and evolutionary innovation in the hominin lineage.

Results

Graph genotyping of structural variation

To assess the role of SVs in human adaptive evolution, we sought to combine the accuracy of long-read sequencing with the scale and population diversity of short-read sequencing data. To this end, we reanalyzed long-read sequencing data of 15 individuals from five continents (Audano et al., 2019) to discover a set of 107,866 SVs (Figure 1—figure supplement 1). The diversity of the long-read sample set (sample ancestries: three African, two American, three East Asian, two European, two South Asian, two hydatidiform moles [likely European]) enables SV discovery across the five continental superpopulations of 1KGP. We found that 89,979 (83.4%) of the long-read-discovered variants, including 30,229 (72.3%) common SVs (allele frequency [AF] ≥ 0.05), are not represented in sets of SVs discovered with short-read sequencing by 1KGP or the Human Genome Diversity Project (HGDP) (Almarri et al., 2020; Sudmant et al., 2015). Despite the much smaller size of the long-read sample set, we were also able to rediscover 66.0% and 17.7% of common SVs found with short reads in the 1KGP and HGDP datasets, respectively, in agreement with results reported in previous studies (Zhao et al., 2021; Figure 1—figure supplement 2 and Figure 1—figure supplement 3; see Materials and methods). We expect that the SVs unique to the short-read datasets reflect both differences in the discovery sample set (i.e., many of the long-read sequenced individuals are in 1KGP, while none are in HGDP) and high rates of false positives in short-read SV detection (Nattestad et al., 2018). We additionally discovered 2.3× and 3.5× more insertions than 1KGP and HGDP, respectively, with an improved insertion-to-deletion ratio (1.46 in the long-read dataset; 0.40 in 1KGP; 0.37 in HGDP). The dearth of insertions in the short-read SV sets reflects the challenge of insertion calling with short reads, which have difficulty mapping sequences not present in the reference genome (Aganezov et al., 2020), while the excess of insertions among long-read SVs reflects a known variant calling bias caused by missing sequences (e.g., collapsed tandem repeats) in the GRCh38 reference (Aganezov et al., 2021).

We then used the graph genotyping software Paragraph (Chen et al., 2019) to genotype this set of long-read discovered SVs in 2504 high-coverage (30×) short-read-sequenced samples from 1KGP (Byrska-Bishop et al., 2021; Figure 1A; see Materials and methods). Results from benchmarking, using ground truth SV data from the Genome in a Bottle Consortium (Zook et al., 2020), show that among both graph-based and non-graph-based SV genotyping tools, Paragraph consistently attains the best balance of precision and recall for both insertions and deletions (Chen et al., 2019; Hickey et al., 2020). Paragraph averages precision and recall rates of 0.72 and 0.70, respectively, across all SV types and genomic regions, and an average precision and recall of 0.86 and 0.79 when repeat regions are excluded (Hickey et al., 2020), underscoring its effectiveness for SV genotyping. Paragraph achieves such accuracy by generating graph representations of SV loci, which include diverging paths for known alternative alleles such as SVs. Short reads are aligned to the graph along the path of best fit, facilitating genotyping even in structurally complex and repetitive regions (Figure 1B). Informed by a large catalog of candidate SV alleles discovered by long-read sequencing, graph genotyping thus permits the study of variants that would be difficult or impossible to discover with short-read data alone (Chen et al., 2019; Hickey et al., 2020; Sibbesen et al., 2018; Sirén et al., 2021).

Figure 1 with 7 supplements see all

Download asset Open asset

Variant graph genotyping of SVs with Paragraph.

(A) Genotyping of SVs was performed using a graph-based approach that represents the reference and alternative alleles of known SVs as edges. The SVs used for graph construction were originally identified from long-read sequencing of 15 individuals (Audano et al., 2019). (B) At candidate SV loci, samples sequenced with short reads are aligned to the graph along the path of best fit, and individuals are genotyped as heterozygous (middle), homozygous for the reference allele (right), or homozygous for the alternative allele (not depicted). We applied this method to the 1KGP dataset to generate population-scale SV genotypes. (C) Allele frequency spectra of SVs genotyped with Paragraph. The left-most bin represents SVs where the alternative allele is absent from the 1KGP sample (AF = 0). Samples are stratified by their 1KGP superpopulation.

Figure 1—source data 1 Superpopulation-level allele frequencies of structural variants used to construct allele frequency spectra in Figure 1C. The third column (‘sum_AC’) provides allele counts, the fourth column (‘sum_NCHROBS’) provides the number of genotyped chromosomes, and the fifth column (‘AF’) provides the alternative allele frequency.: https://cdn.elifesciences.org/articles/67615/elife-67615-fig1-data1-v3.txt
Download elife-67615-fig1-data1-v3.txt

As quality control, we filtered the resulting data based on genotyping rates and adherence to Hardy–Weinberg equilibrium, in accordance with Chen et al., 2019 (see Materials and methods). Specifically, we removed SVs that were not successfully genotyped in ≥50% of samples, as well as SVs that violated one-sided Hardy–Weinberg equilibrium expectations (excess of heterozygotes) in more than half of the 1KGP populations. The latter scenario, whereby nearly all individuals are genotyped as heterozygous, is a common artifact caused by slightly divergent repetitive sequences that are falsely interpreted as alternative alleles at a single locus (Graffelman et al., 2017). These filtering steps removed 15,580 SVs, leaving 92,286 variants for downstream analysis.

Using this remaining set, we examined the allele frequency spectrum both to provide a general description of the data, as well as to quantify the impacts of our unique ascertainment approach (Figure 1C, Figure 1—figure supplement 5). Global alternative allele frequencies were strongly correlated with the number of long-read-sequenced samples in which they were originally discovered (Audano et al., 2019), broadly supporting the accuracy of our graph genotyping results (Figure 1—figure supplement 4). A total of 25,201 (27.3%) alternative SV alleles were absent from all 1KGP samples, likely reflecting ultra-rare variation within the panel of long-read sequenced genomes, as well as the imperfect recall (i.e., true positive rate) of SV genotyping (Chen et al., 2019). This abundance of rare variation is a known feature of human populations, which have experienced rapid demographic expansion over recent history (Keinan and Clark, 2012). In contrast, a total of 1139 (1.2%) alternative SV alleles were observed to be fixed in the 1KGP sample and likely reflect a combination of assembly errors and rare variation in the reference genome. Relaxing the criterion for fixation to 90% alternative allele frequency (thus allowing for a modest rate of false negatives) resulted in a total of 4947 (5.4%) fixed or nearly fixed alternative SV alleles.

The frequency spectrum of segregating variation was meanwhile shaped by an ascertainment bias caused by variant discovery within a small sample and genotyping within a separate and larger sample. This bias, which is similar to that previously described for genotyping microarrays (Lachance and Tishkoff, 2013), is characterized by an apparent depletion of rare variation, given that such variants are not shared with the discovery sample. While important to note, our study should be largely unaffected by this bias, due to our focus on individual variants rather than the shape of the allele frequency spectrum. Moreover, positively selected variants are by definition locally common and thus enriched within a small but globally diverse sample (Audano et al., 2019).

We next sought to test whether different classes of SVs exhibited differences in allele frequencies as a potential consequence of purifying selection. Among SVs that were polymorphic within the 1KGP sample, we observed a negative correlation between SV length and minor allele frequency (Kendall’s τ = –0.140, p-value < 1 × 10^–10), consistent with selection against longer SVs (Figure 1—figure supplement 6). We similarly observed that deletions segregated at lower average minor allele frequencies than insertions (β = –0.555, p-value < 1 ×10^–10).

Quantifying linkage disequilibrium with known SNPs and short indels

One intriguing question regards the extent to which SV genotypes are correlated with those of known SNPs and indels, which has implications for the novelty of any potential discoveries based on the SV data. To address this question, we quantified linkage disequilibrium (LD) between the catalog of long-read discovered SVs and known SNPs and short indels from 1KGP. For each of the 26 1KGP populations, we computed the maximum observed LD between each SV and the nearest 100 variants within a 1 Mb window. Depending on the population subject to measurement, 36–41% of segregating SVs were in strong LD (r² > 0.8) with any known SNP or short indel, and 52–57% were in moderate LD (r² > 0.5) with one of these variants (Figure 1—figure supplement 7). These observations of maximum LD with SNPs are slightly lower than those computed by Jakubosky et al., 2020a, which may reflect our ability to discover SVs in regions that are less accessible to short reads. Levels of LD were lowest for African populations, in accordance with known patterns of haplotype structure (Conrad et al., 2006). As with SNPs, these lower levels of LD suggest that the accuracy of SV genotype imputation based on known SNPs will be lowest in African populations, but that resolution for fine-mapping SVs that are potentially causal in phenotypic associations will also be highest (Wojcik et al., 2019).

We emphasize that our LD observations are strongly affected by the challenge of small variant discovery and genotyping in repetitive regions of the genome that are enriched for SVs. Specifically, 50,917 of all SVs (47.2%) intersect at least partially with one or more genomic intervals deemed inaccessible based on the 1KGP pilot mask, and 5214 SVs (4.8%) occur within regions of the genome identified as problematic by the ENCODE Consortium (Amemiya et al., 2019). Nevertheless, LD with known SNPs and indels is a meaningful metric for our study, because it quantifies the extent to which SVs represent independent and unexplored variation relative to previous evolutionary studies. Low observed levels of LD between a substantial fraction of SVs and known variants presents an opportunity to discover novel functional associations and signatures of adaptation at loci poorly tagged by easily genotyped markers.

Expression quantitative trait locus (eQTL) mapping

Seeking to first test the functional impacts of common structural variation, we intersected the SV genotype data with RNA-seq data from an overlapping set of 447 samples from the Geuvadis Consortium, which was generated from lymphoblastoid cell lines (LCLs) derived from individuals from four European and one African population (Lappalainen et al., 2013). We tested for associations between levels of gene expression and genotypes for SVs within 1 Mb from the transcription start site (TSS). After filtering the data on genotyping call rate, minor allele frequency, and gene expression level (see Materials and methods), we identified a total of 1121 SV-gene pairs with significant gene expression associations at a 10% false discovery rate (FDR; Figure 2A), broadly consistent with expectations from previous eQTL studies when scaled to the number of tested SVs (Chiang et al., 2017). SVs with significant impacts on expression tended to occur near the genes that they regulate, with 62% of significant SV eQTLs occurring within 100 kb of the corresponding TSS (Figure 2B).

Figure 2 with 2 supplements see all

Download asset Open asset

eQTL mapping of SVs.

We used RNA-seq data from the Geuvadis Consortium (Lappalainen et al., 2013), obtained from LCLs derived from individuals from four European and one African population of the 1KGP dataset, to test for associations between SV genotypes and gene expression. SV-eQTL pairs that were significant at a 10% FDR are depicted in purple. (A) Q–Q plot of permutation p-values for all SV-gene pairs tested. (B) Volcano plot of eQTLs and the estimated effect of the alternative allele on expression (β). (C) Distribution of the distance of significant SV eQTLs from the transcription start site (TSS) of their associated genes. (D) Enrichment or depletion of expression associations for SVs that overlap various ChromHMM chromatin state annotations from the Roadmap Epigenomics Project (GM12878 Lymphoblastoid Cells). Chromatin states with ≤0.1% genome-wide representation were omitted for visual clarity.

Figure 2—source data 1 Summarized eQTL data underlying Figure 2. SV IDs, gene IDs and names, distance to the TSS, slopes, p-values, and q-values are provided.: https://cdn.elifesciences.org/articles/67615/elife-67615-fig2-data1-v3.txt
Download elife-67615-fig2-data1-v3.txt

Top gene expression associations include several genes in the HLA complex (HLA-DQB2, HLA-DQA2, HLA-DRB6, and HLA-P), as well as ERAP2, which encodes an endoplasmic reticulum aminopeptidase that processes HLA ligands to lengths suitable for their binding (Figure 2C). While gene expression data from LCLs provides a limited snapshot of the functional implications of SVs, it is notable that these immune loci also exhibit strong gene expression diversity mediated by genetic variation.

We next investigated the mechanisms through which SV eQTLs may impact gene expression. Using epigenetic state annotations generated by the Roadmap Epigenetics Consortium (Roadmap Epigenomics Consortium et al., 2015), we found a strong enrichment of significant gene expression associations for SVs that intersect annotated enhancer elements, as well as a strong depletion of expression associations for SVs that intersect ‘quiescent’ sequences that are devoid of important epigenetic marks (Figure 2D). We also further investigated the 13 significant SV eQTLs that intersected one or more exons of their associated gene. Of the nine such deletions, 8 exhibited associations consistent with direct dosage effects on expression (i.e., a linear association between expression and deletion copy number; Figure 2—figure supplement 1A). Meanwhile, the three exon-intersecting insertions exhibited associations in both directions (one positive and two negative associations; Figure 2—figure supplement 1B). Together, these intuitive results shed light on functional mechanisms through which SVs may mediate impacts on phenotypes and fitness.

Lastly, we performed fine-mapping of SV and SNP eQTLs with CAVIAR (Hormozdiari et al., 2014) with the goal of homing in on candidate causal variants at SV eQTL loci. Of the 1121 loci with a significant SV eQTL, 345 (30.8%) contained the SV within the 90% credible causal set. While the median size of the credible causal sets was 19 variants, we identified 15 loci where an SV occurred among a credible causal set of five or fewer variants (Figure 2—figure supplement 2). These include a 37 Kb deletion (27005_CHM13_del) that is among just two variants in a credible causal set for expression effects on MIF-AS1, a lncRNA partially encompassed by the deletion. The expression of this lncRNA has been previously associated with breast and gastric cancer proliferation (Ding et al., 2019; Li et al., 2018). Similarly, a 3 kb insertion (20290_HX1_ins) is among two variants in a credible causal set for expression effects on CSF2RB, the dysfunction of which has been associated with pulmonary disease (Suzuki et al., 2011). While this approach offers a roadmap for future fine-mapping studies investigating the relative impacts of different variant classes on gene expression, we caution that the present analysis may still underestimate the potential causal role of SVs, due to a modestly higher error rate for genotyping of SVs versus SNPs (Chen et al., 2019).

Admixture-aware scan for signatures of local adaptation

We next sought to use our set of population-scale SV genotypes to search for SVs with signatures of local adaptation. Such locus-specific scans for local adaptation can be broadly classified into frequency differentiation-based and LD-based approaches (Vitti et al., 2013). While powerful, LD-based methods generally require haplotype phasing and/or dense and accurate genotyping – a requirement that is infeasible for many of the complex and repetitive regions of the genome investigated in our study. In contrast, frequency differentiation-based approaches can be applied to individual loci and are based on the logic that positive selection tends to cause a particular allele to increase in frequency only in the population(s) where it became established and was advantageous. However, a limitation of these methods is the requirement that individuals are grouped into pre-defined populations that may not reflect genetic patterns or the substantial admixture exhibited by many human populations.

To overcome these limitations, we used Ohana (Cheng et al., 2017; Ilardo et al., 2018; Cheng et al., 2019), a maximum likelihood-based method that models individuals as possessing ancestry from combinations of k ancestry components, inspired by related methods (Alexander et al., 2009; Falush et al., 2003). The method then tests whether individual variants adhere to this genome-wide null model, or are better explained by an alternative model in which frequencies are allowed to vary in one or more populations by consequence of local adaptation. Following the precedent set by 1KGP (The 1000 Genomes Project Consortium, 2015), we modeled individual genomes as combinations of eight ancestry components, replicating known patterns of population structure at continental scales as well as prominent signatures of admixture within a subset of populations (e.g., African ancestry in Southwest US [ASW] and admixed American populations [AMR]; Figure 3A). The distribution of ancestry components across populations was qualitatively unaffected by the choice of k (Figure 3—figure supplement 1).

Figure 3 with 3 supplements see all

Download asset Open asset

Quantifying allele frequency differentiation among ancestry components.

(A) To conduct an admixture-aware scan for local adaptation, we used Ohana (Cheng et al., 2017; Ilardo et al., 2018) to infer genome-wide patterns of ancestry in the 1KGP samples. This method models each individual as a combination of k ancestry components and then searches for evidence of local adaptation on these component lineages. Admixture proportions (k = 8) for all samples in the 1KGP dataset, grouped by population. Vertical bars represent individual genomes. (B) Using Ohana, we searched for evidence of local adaptation by testing whether the allele frequencies of individual variants were better explained by the genome-wide covariance matrix, or by an alternative covariance matrix where allele frequencies were allowed to vary in one ancestry component. The likelihood ratio statistic (LRS) reflects the relative support for the latter selection hypothesis. For each ancestry component, SVs with LRS > 32 are colored by their maximum linkage disequilibrium (r²) with any known SNP in the corresponding 1KGP superpopulation.

Figure 3—source data 1 Ancestry component results underlying Figure 3A. Ancestry proportions for each individual are provided, along with that individual’s 1KGP population code.: https://cdn.elifesciences.org/articles/67615/elife-67615-fig3-data1-v3.txt
Download elife-67615-fig3-data1-v3.txt

Across all ancestry components, we identified 220 unique SVs exhibiting significant deviation from genome-wide patterns of allele frequency differentiation (99.9th percentile of matched distribution for SNPs and short indels; see Materials and methods; Figure 3B; Supplementary file 1; Figure 3—figure supplement 1). These included 139 SVs with coordinates overlapping with annotated genes, of which 13 intersected with annotated exons. Seven SVs at these frequency-differentiated loci were also significant eQTLs based on our previous analysis of the Geuvadis LCL data (Supplementary file 2). Through comparison to results for SNPs and indels, we found that 88 of the outlier SVs we identified were among the top 10 most frequency-differentiated variants within a 1 Mb window, suggesting that these represent examples where SVs may be causal targets of selection at adaptive loci. Finally, only 119 (54.1%) highly differentiated SVs possessed strong LD with known SNPs or short indels (r² > 0.8; Figure 1—figure supplement 2), indicating that many of these loci constitute novel candidate targets of selection that may have been missed by previous scans.

Known and novel targets of local adaptation

Notable examples of highly differentiated loci included rs333, the Δ32 allele of the chemokine receptor CCR5, which is known to confer resistance to HIV infection and progression (Dean et al., 1996). Among our results, this deletion polymorphism is the 14th most frequency-differentiated SV with respect to ancestry component 3, which is highly represented in Europe. The CCR5-Δ32 allele segregates at moderate frequencies in European populations (MAF = 10.9%) and achieves its highest frequency of 15.6% in the Finnish population, but segregates at low frequencies elsewhere. The case for historical positive selection at this locus has been contentious, with initial studies citing a geographic cline in allele frequencies and strong LD with adjacent microsatellite markers (Stephens et al., 1998), potentially driven by epidemics such as the bubonic plague or smallpox (Galvani and Slatkin, 2003). However, subsequent studies argued that patterns of long-range LD (Sabeti et al., 2005) and temporal allele frequency changes based on ancient DNA (aDNA) samples (Bollback et al., 2008) could not exclude models of neutral evolution.

Our study also identified numerous novel hits, such as a 309 bp intronic insertion in TMEM131 (ancestry component 3), which segregates at a frequency of 63% in Finnish populations, but at a mean frequency of 24% in non-European populations. This gene encodes proteins involved in collagen cargo trafficking from the endoplasmic reticulum to the Golgi (Zhang et al., 2020). Collagen is the most abundant protein in the human body and the major component of human skin. Recurrent positive selection shaping skin pigmentation and other phenotypes is one of the best described examples of local adaptation across human populations (Crawford et al., 2017; Jablonski, 2004).

We also identified a 2.8 kb insertion in an intron of CLEC16A, inferred to be under selection in ancestry component 4, which corresponds to the Peruvian (PEL), Mexican (MXL), and Colombian (CLM) populations of 1KGP. The insertion, which is among the strongest frequency-differentiated variants at its respective locus and segregates in LD with several SNPs and short indels (maximum r² = 0.81; Figure 3—figure supplement 3), occurs at frequencies of 52.4%, 38.3%, and 15.4%, respectively, in these three populations, but is rare in others (AF < 0.04). CLEC16A is thought to impact susceptibility to autoimmune disorders, and SNPs in this gene have been associated with diseases such as type I diabetes, multiple sclerosis, and rheumatoid arthritis (Pandey et al., 2019). Notably, this same SV was also identified in a recent study of structural variation, which similarly reported that it achieves high frequency in Peruvian populations (Ebert et al., 2021).

In rare cases, multiple SVs in LD with one another captured the same underlying signature of frequency differentiation. One such example from South Asian populations (ancestry component 5) involved a linked intronic insertion and deletion in the cellular growth and morphogenesis related gene CSNK1G1. In this case, the reference genome carries the global minor allele, such that the signature of local adaptation presents as a lower frequency of the alternative allele in South Asian populations (60–71% for 22980_HG00514_ins; 60%–72% for 25014_HG02106_del) compared to other global populations (91% and 92%, respectively).

Extreme signatures of adaptation at the immunoglobulin locus in southeast Asian populations

The SV with the strongest evidence of local adaptation across all populations was an insertion polymorphism in an intron of IGHG4, which codes for a constant domain of the immunoglobulin heavy chain. These heavy chains pair with light chains, the latter of which include a domain composed of variable (V), diversity (D), and joining (J) segments. Complementing their substantial germline variation, V(D)J loci experience somatic recombination and hypermutation to generate vast antibody repertoires – the defining feature of the adaptive immune system (Watson et al., 2017). This insertion polymorphism identified by our scan exhibits strong allele frequency differentiation in ancestry component 2, which is highly represented in the Chinese Dai in Xishuangbanna, China (CDX) and Kinh in Ho Chi Minh City, Vietnam (KHV) populations, where it achieves frequencies of 88% and 65%, respectively, while remaining at much lower allele frequencies in other global populations (Figure 4A,B).

Figure 4 with 8 supplements see all

Download asset Open asset

Evidence for Neanderthal introgression of the adaptive *IGH* haplotype.

(A) Local analysis of likelihood ratio statistics (LRS) in the region near the 33 bp insertion (red point) reveals a 325 kb haplotype encompassing 94 SNPs with strong allele frequency differentiation within ancestry component 2. Points where the alternative allele matches an allele observed in the Chagryskaya Neanderthal genome but at a frequency of 1% or less in African populations are highlighted in purple. (B) Individual haplotypes defined by the highly differentiated SNPs (LRS > 450). Four archaic hominin genomes are plotted at the top, while 30 randomly sampled haplotypes from each of 6 populations from 1KGP are plotted below. Archaic hominins samples are colored according to whether they possess more than one aligned read supporting the alternative allele at a given site. ESN refers to the Esan in Nigeria population of 1KGP. Other 1KGP population codes are provided in the main text.

The SV was originally reported as a 34 bp insertion based on long-read sequencing of the genome of a Vietnamese individual (HG02059) (Audano et al., 2019). Based on realignment to a modified version of the reference genome that includes the alternative allele, we revised the sequence of this insertion to 33 bp, but found that it is well supported by patterns of coverage, split reads, and soft clipped alignments at the SV breakpoints (Figure 4—figure supplement 1). The sequence of the insertion itself is repetitive within the human reference genome, with two identical copies of the sequence occurring ~44 and ~117 kb downstream, respectively, also within the IGH gene cluster. This SV is included among sets of indels called by gnomAD (Karczewski et al., 2020), 1KGP, and HGDP, with a slightly altered position and insertion sequence. However, it exhibits a highly skewed allele balance (centered on 0.2, in contrast to the expectation of 0.5) – a symptom of reference bias in read mapping (Chen et al., 2019) – as well as low genotyping rates (31% for CDX in 1KGP, compared with 85% in our study). These biases result in lower estimates of the insertion’s allele frequency in East Asian populations (AF = 0.73 in CDX) and underscore the technical challenges of genotyping at this locus with traditional short-read approaches (Chin et al., 2020; Zhang et al., 2021).

Interestingly, the second strongest signature of adaptation across all populations traced to a nearby 135 bp deletion, which overlaps with two transcription factor binding sites for IGHE, another component of the constant region of immunoglobulin heavy chains. This deletion, which lies only 33 kb upstream of the IGHG4 insertion (Figure 4A), again achieves high frequencies in the CDX and KHV populations (70% and 54%, respectively) but segregates at low frequencies in most other global populations. Despite their genomic proximity and similar patterns of frequency differentiation, these two SVs exhibit only modest levels of LD (r² = 0.17), likely reflecting recombination occurring during or after the episode of selection.

Evidence of Neanderthal-introgressed origin of the high-frequency IGH haplotype

Approximately 2% of the genome of all non-African individuals traces to admixture with Neanderthals between 47 and 65 kya, while Oceanian populations (and Asian populations to a lesser extent) also possess sequences introgressed from Denisovans (Green et al., 2010; Sankararaman et al., 2012; Sankararaman et al., 2016; Vernot et al., 2016). Introgressed alleles, including some SVs (Almarri et al., 2020; Hsieh et al., 2019), are thought to have conferred both beneficial and deleterious effects on modern human populations, especially with respect to immune-related phenotypes (Rotival and Quintana-Murci, 2020). Motivated by these findings, we tested our set of highly differentiated SVs for evidence of archaic hominin introgression. Using results from the method Sprime (Browning et al., 2018), which leverages patterns of divergence and haplotype structure to classify archaic introgressed sequences, we identified SVs that segregate in LD with putative introgressed SNPs (see Materials and methods). Of the 220 highly differentiated SVs, 26 (12%) exhibited moderate LD (r² > 0.5) with putative introgressed haplotypes, while simultaneously exhibiting low allele frequencies (AF < 0.01) within African populations of 1KGP (Figure 4—figure supplement 2; Supplementary file 3). Notably, this set of candidate introgressed SVs included the IGHG4 insertion and nearby deletion, which, despite their low LD with one another, each tag multiple putative introgressed SNPs within the CDX population. Indeed, the original Sprime publication reported that putative introgressed variants of IGHA1, IGHG1, and IGHG3 achieve high frequencies in Eurasian populations (Browning et al., 2018). The variants identified by Browning et al., 2018 are located in a subregion of the broader IGH locus, downstream of the introgressed SVs, and segregate at high allele frequency in East Asian, European, and American populations (Figure 4—figure supplement 3). The southeast Asian-specific haplotype we identify, which includes the IGHG4 insertion and nearby deletion, may have been challenging to discover due to the difficulties of short-read alignment and genotyping in this region of the genome (Zhang et al., 2021). Indeed, 80.7% of the sequence in the broader IGH locus was filtered out by Browning et al., 2018 through strict masking of aDNA genotypes to remove low-coverage, poorly mapping, or repeat-associated reads (Figure 4—figure supplement 4).

Most candidate introgressed SVs fall within complex regions of the genome, such that genotyping directly from fragmented and degraded aDNA data remains infeasible. However, the IGHG4 intronic insertion is short enough to be spanned by individual sequencing reads, thus allowing us to validate our genotyping results and further investigate introgression at this locus by directly searching for an exact match to a diagnostic sequence. Specifically, we identified a 48 bp sequence that is unique to individuals with the insertion, but not observed in the reference genome (see Materials and methods), thereby facilitating rapid searches for the insertion allele in any sequencing dataset. Application of this approach to the 1KGP data broadly validated our graph genotyping results and supported the strong allele frequency differences that we had previously described (Figure 4—figure supplement 5). We then extended this scan to published short-read sequencing data from four high-coverage (~30×) archaic hominin genomes (Mafessoni et al., 2020; Meyer et al., 2012; Prüfer et al., 2017; Prüfer et al., 2014). While the diagnostic sequence was absent from the Denisovan genome, we observed one or more perfect matches in the sequences of each of three Neanderthals, which we found notable given its absence among African populations (Figure 4—figure supplement 5).

We next expanded our analysis to the genomic region around the IGHG4 insertion, investigating the archaic hominin allelic states at each of the 109 highly differentiated variants (2 SVs, 14 short indels, and 93 SNPs) defining the 325 kb LRS peak at the IGH locus. Focusing on the 93 SNPs where archaic alleles are easily compared, we observed that the Denisovan genome exhibited a modest degree of matching (39 shared alleles [42%]), while the Altai, Vindija, and Chagyrskaya Neanderthal genomes exhibited near perfect matching over the entirety of the differentiated region (77 [83%], 88 [95%], and 89 [96%] shared alleles, respectively) (Figure 4B). Additionally, of these 93 highly differentiated SNPs, 64 were called as introgressed in CDX by Sprime (Figure 4—figure supplement 6). Combined with the near absence of these alleles from other global populations and the length of the differentiated haplotype, our observations strongly support the conclusion that this sequence originated via ancient introgression from a Neanderthal population related to the Chagyrskaya and Vindija Neanderthals.

Examination of the haplotype structure at the IGH locus revealed four discrete blocks of LD within the highly differentiated region (Figure 4—figure supplement 7), consistent with substantial recombination after the original haplotype achieved high frequency. The deletion SV (22231_HG02059_del) falls within the largest LD block, while the insertion SV (22237_HG02059_ins) falls within a smaller LD block that exhibits the greatest allele frequency differences. The latter block includes the tag SNP rs150526114, where the global minor allele matches the Neanderthal genomes and segregates at 91% frequency in CDX, 73% frequency in KHV, and 59% frequency in CHS, but is rare or absent in most other populations from 1KGP. Data from HGDP (Bergström et al., 2020) and the Simons Genome Diversity Panel (SGDP) (Mallick et al., 2016) shed additional light on the geographic distribution of the putative Neanderthal-introgressed allele and confirmed its extreme pattern of frequency differentiation specific to southeast Asian populations (Figure 4—figure supplement 8). Notably, the allele is absent from the HGDP populations from the Americas, which are thought to have split from East Asian populations approximately 26 kya (Moreno-Mayar et al., 2018), further supporting the recency and geographically restricted nature of this positive selection event.

Strength and timing of adaptive introgression at the IGH locus

The challenge of dense genotyping and phasing within the broader IGH region (Zhang et al., 2021) hinders haplotype-based approaches for inferring the timing of selection. Nevertheless, the global patterns of allele frequencies we observe (Figure 5A,B) can be interpreted in the context of known population demographic histories, providing a rough estimate of such timing. For example, pairwise divergence times among East Asian populations inferred by Wang et al., 2018 constrain the plausible timing of selection at the IGH locus between approximately 60 generations (1740 years) and 400 generations (11,600 years) ago (assuming a 29-year generation time Tremblay and Vézina, 2000): after the divergence of the Chinese Dai and Vietnamese populations from the Japanese (JPT) population, but before the divergence of CDX and KHV. These divergence time estimates are roughly consistent with those inferred using an alternative approach based on the joint allele frequency spectrum by Jouganous et al., 2017.

Figure 5 with 2 supplements see all

Download asset Open asset

Local adaptation at the *IGH* locus.

(A) Population-specific frequencies of the putative Neanderthal-introgressed insertion allele in each of the 1KGP populations, in the style of the Geography of Genetic Variants browser (Marcus and Novembre, 2017). (B) Tree representation of the best-fit selection hypothesis for Neanderthal-introgressed haplotype tagged by the *IGHG4* insertion polymorphism, as computed by Ohana. (C) Five-population demographic model used for simulation and parameter inference via approximate Bayesian computation (ABC). Population sizes and split times are further described in Materials and methods. (D) Posterior distribution of the selection coefficient (s). (E) Posterior distribution of the timing of the onset of selection (T_adaptive). (F) Posterior distribution of the initial allele frequency at the beginning of the simulation (p₀). (G) Negative relationship between s and T_adaptive for simulations retained by ABC.

Figure 5—source data 1 Parameters and summary statistics for the simulations retained in the ABC analysis, underlying Figure 5D–G. Simulation parameters are provided in the first three columns. Population-specific allele frequencies output by the simulation are provided in the last five columns, labeled by their corresponding 1KGP population codes.: https://cdn.elifesciences.org/articles/67615/elife-67615-fig5-data1-v3.txt
Download elife-67615-fig5-data1-v3.txt

To achieve further insight into this episode of selection, we conducted forward genetic simulations of a simplified five-population demographic model based on parameters obtained from the literature (The 1000 Genomes Project Consortium, 2015; Jouganous et al., 2017; Wang et al., 2018; Figure 5C). For simplicity, we began the simulations 46 kya, after Neanderthal introgression has occurred and the introgressed haplotype is segregating within the modern human population. Our model makes an additional simplification by not including migration, but we note that this omission should make our estimates of the selection coefficient conservative, as stronger selection is required to generate allele frequency differences between populations exchanging migrants. We first performed 100,000 neutral simulations (selection coefficient [s] = 0), finding that the magnitudes of allele frequency differentiation, especially in southeast Asia, far exceed those expected under the null (Figure 5—figure supplement 1). In order to formally estimate the parameters of the selection event, we next applied approximate Bayesian computation (ABC) to an additional 175,000 simulations in which we varied relevant parameters (see Materials and methods). In line with our previous intuition, the results of this analysis indicated a recent onset of extremely strong selection, with the selection coefficient parameter (s) inferred to be 0.06 (95% credible interval [CI] [0.02, 0.12]), the selection onset time parameter (T_adaptive) inferred to be 4400 years ago (95% CI [1700, 8400]), and the initial frequency of the adaptive allele in the ancestral Eurasian population (p₀) to be 0.18 (95% CI [0.04, 0.35]) (Figure 5D–F). We observed a strong correlation between estimates of the selection coefficient and timing of selection. Specifically, older, weaker selection could produce the same frequency differences as stronger and more recent selection (Figure 5G), though even the lower bound on our estimate of s places it among the strongest episodes of positive selection documented in humans.

Discussion

Long-read sequencing is starting to provide more comprehensive views of the landscape of human genetic variation, drawing novel links to phenotypes and diseases. However, long-read sequencing methods remain impractical for most population-scale studies due to their low throughput and high cost. We sought to overcome these limitations by applying variant graph genotyping of a large catalog of long-read-discovered SVs to short-read sequencing data from globally diverse individuals of 1KGP. By mapping eQTLs and scanning for evidence of local adaptation and adaptive introgression, we highlighted the role of SVs as largely unexplored contributors to variation in human genome function and fitness.

Despite the scale and diversity of SVs and samples used in our study, we anticipate that additional SV targets of selection remain undiscovered, either because they were not present in the set of long-read sequenced individuals or because they remain inaccessible to genotyping using graph-based approaches (e.g., tandem repeats) (Chen et al., 2019). Studying the evolutionary impacts of such SVs will require the application of long-read sequencing at population scales (Ebert et al., 2021), which is still infeasible for most genetics studies. Moreover, SV loci exhibiting expression associations or signatures of selection may not themselves be the causal targets, but rather tag nearby causal variation by consequence of LD. Methods such as fine mapping (Schaid et al., 2018) and multiplex reporter assays (Tewhey et al., 2018; van Arensbergen et al., 2019) will be invaluable for disentangling LD to reveal causal relationships and contrast the relative impacts of various forms of genetic variation.

The two strongest signatures of local adaptation in our study traced to the immunoglobulin heavy chain locus. While the precise phenotypic impacts of these variants remain unknown (Figure 5—figure supplement 2), their potential effects on adaptive immunity is intriguing given the established role of immune-related genes as common targets of local adaptation in human populations (Barreiro and Quintana-Murci, 2020). The human IGH locus is highly polymorphic (Mikocziova et al., 2020), with examples of SNPs and copy number variants exhibiting frequency differences between populations (Watson et al., 2013). In developing lymphocytes, this locus undergoes somatic V(D)J recombination and hypermutation to produce antibodies that drive the immune response—processes that may be influenced by nearby germline variation (Watson et al., 2017). The combination of these forms of variation makes the region difficult to probe with traditional sequencing methods (Chin et al., 2020; Zhang et al., 2021), in turn highlighting the power of long-read sequencing and graph genotyping.

Our observation of the IGHG4 insertion within the Neanderthal genomes allowed us to connect and build upon initial evidence of selection and introgression at the IGH locus. Specifically, 1KGP previously reported SNPs in several immunoglobulin genes as allele frequency outliers with respect to the CDX population (The 1000 Genomes Project Consortium, 2015). Browning et al., 2018 later noted that a Neanderthal-introgressed haplotype at the IGH locus achieves high frequencies in Eurasian populations, although they focused on a separate haplotype in the region, split by recombination, that is not specific to southeast Asia. Our findings further expand the evolutionary history of this region of the genome, as well as adding to a growing list of examples of adaptive introgression from archaic hominins (Gittelman et al., 2016; Hsieh et al., 2019; Huerta-Sánchez et al., 2014; Racimo et al., 2017), several of which are thought to have targeted immune-related phenotypes (Abi-Rached et al., 2011; Dannemann et al., 2016; Enard and Petrov, 2018; Gouy and Excoffier, 2020; Mendez et al., 2012a; Mendez et al., 2012b; Sams et al., 2021).

Based on forward genetic simulations, we estimated that selection on the introgressed IGH haplotype initiated between 1700 and 8400 years ago, before the divergence of the CDX and KHV populations and in line with our intuition based on patterns of allele frequencies. This recent onset of selection is intriguing given that introgression from Neanderthals into the ancestors of Eurasian modern human populations dates to 47–65 kya (Sankararaman et al., 2012). Our findings thus suggest that persisting archaic introgressed haplotypes provided a reservoir of functional variation to the ancestors of CDX and KHV that proved adaptive during a period of environmental change, for example in response to local pathogens (Rasmussen et al., 2015). Recent reports that Neanderthal-introgressed sequences mediate individual outcomes of SARS-CoV-2 infections to this day lend plausibility to this hypothesis (Zeberg and Pääbo, 2021; Zeberg and Pääbo, 2020; Zhou et al., 2021), as does polygenic evidence of adaptation in response to ancient viral epidemics, including in East Asia (Souilmi et al., 2021). Moreover, the delayed onset of selection inferred in our study is consistent with recent evidence of time lags between introgression and adaptation favoring other Neanderthal-introgressed alleles (Yair et al., 2021).

Our simulations additionally demonstrated that a selection coefficient between 0.02 and 0.12 best explains the observed frequency differences, comparable to other episodes of strong selection in humans. These include examples such as lactase persistence mutations near LCT (0.01–0.15) (Bersaglieri et al., 2004) and malaria resistance mutations affecting DARC (0.08) (Hamid et al., 2021) and HBB (0.1) (Elguero et al., 2015). We caution against overinterpretation of these parameter estimates given the uncertainty in the underlying demographic model, which relies on estimates of population size and divergence time from two studies (Jouganous et al., 2017; Wang et al., 2018) and does not incorporate admixture between populations. However, our results are broadly consistent with the observation of extreme allele frequency differences among closely related populations in East Asia. Future studies incorporating more complex evolutionary models and fully resolved IGH haplotypes (Rodriguez et al., 2020) will be essential for further refining the complex evolutionary history of selection and recombination at the immunoglobulin locus. Nevertheless, the coexistence of V(D)J recombination, somatic hypermutation, and local adaptation at this locus presents a remarkable example of diversifying selection at multiple scales of biological organization, generating allelic diversity both within individuals and across populations.

Together, our study demonstrates how new sequencing technologies and bioinformatic algorithms are facilitating understanding of complex and repetitive regions of the genome – a new frontier for human population genetics. Our study also demonstrates how the integration of short and long-read datasets allows for investigation of these regions at population scales. Combined with studies of diverse populations, these technologies are providing a more complete picture of human genomes and the evolutionary forces by which they are shaped.

Materials and methods

Structural variant calling and comparison with short-read datasets

Request a detailed protocol

We used published long-read sequencing data from 15 individuals to generate a set of 107,866 SVs for graph genotyping (Audano et al., 2019). Raw reads were downloaded using the accessions provided in the original publication and aligned with NGM-LR (Sedlazeck et al., 2018b) using default PacBio parameters to the main chromosomes of GRCh38. Variants were called with Sniffles (Sedlazeck et al., 2018b), requiring a minimum SV length of 30 bp and a minimum of 10 supporting reads. Resulting VCFs were refined with Iris (Alonge et al., 2020) to polish the reported SV sequences. Variants were then merged with SURVIVOR v1.0.7 (Jeffares et al., 2017), using a merge distance of 50 bp and requiring strand and type to match. For each merged variant, a representative variant was then obtained from the original pre-merged call set to improve accuracy. Such representative variants were selected by first prioritizing homozygous over heterozygous calls, and then by prioritizing variants with greater proportions of reads supporting the non-reference allele. To prepare the variants for input into Paragraph, translocations, mitochondrial DNA variants, inversions and duplications over 5 kb, and variants without a ‘PASS’ filter were removed from the VCF. This resulted in a set of 107,866 SVs.

We compared this set of long-read discovered SVs with SVs discovered from short-read sequencing in 1KGP (Sudmant et al., 2015; http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/phase3/integrated_sv_map/supporting/GRCh38_positions/) and HGDP (Almarri et al., 2020; ftp://ngs.sanger.ac.uk/production/hgdp/hgdp_structural_variation). We merged long-read and short-read SV sets using parameters from Aganezov et al., 2020 (maximum distance of 500 bp between SV start sites, maximum difference of 10% between SV lengths). We chose not to require matching on SV type, strand orientation, or insertion sequence due to inconsistencies in SV representations across the short-read and long-read variant callers (Aganezov et al., 2020). In order to match the minimum SV length used by the short-read datasets, we also removed all long-read discovered SVs shorter than 50 bp. For the purpose of broad allele frequency binning for SVs called in each of these sets, we then calculated dataset-wide allele frequencies with PLINK v1.90b6.4 (Purcell et al., 2007).

Graph genotyping of structural variation

Request a detailed protocol

High-coverage 30× short-read sequencing data for the core 2504 individuals in 1KGP, sequenced by the New York Genome Center (Byrska-Bishop et al., 2021), was obtained from ENA (PRJEB31736). We genotyped SVs in these samples with Paragraph v2.2 (Chen et al., 2019). In accordance with Paragraph’s recommendations, we set the maximum permitted read count for variants to 20 times the mean sample depth in order to limit runtime for repetitive regions. Genotypes from all samples were combined using bcftools v1.9 (Danecek et al., 2021).

To obtain a high-quality set of genotyped SVs, we filtered the resulting data based on dataset-wide genotyping rates and within-population Hardy–Weinberg equilibrium. We determined an SV’s overall genotyping rate with cyvcf2 (Pedersen and Quinlan, 2017) and removed variants that were not genotyped in ≥50% of samples. We additionally calculated one-sided Hardy–Weinberg equilibrium p-values (excess of heterozygotes) for variants within each of the 26 1KGP populations, using the Hardy–Weinberg package from R (Graffelman, 2015). We filtered out SVs that violated equilibrium expectations (Fisher’s exact test, p < 1 × 10^–4) in ≥13 populations. Unfolded, within-population allele frequencies were calculated with PLINK.

Quantifying linkage disequilibrium with known SNPs and short indels

Request a detailed protocol

To calculate linkage disequilibrium (LD) between SVs and SNPs or short indels in 1KGP samples, we used small variant genotypes produced by the 1000 Genomes Consortium. These genotypes were generated by aligning the 1KGP Phase 3 data to GRCh38 and then calling variants against the GRCh38 reference and are restricted to biallelic SNVs and indels (http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/1000_genomes_project/release/20190312_biallelic_SNV_and_INDEL/). Y chromosome genotypes, not included in the former data release, were obtained from Phase 3 variant calls lifted over to GRCh38 (http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/supporting/GRCh38_positions/). We combined the 1KGP SNP and indel genotypes with our SV genotypes and calculated r² within each population, using a 100-variant and 100 Mb window, with PLINK v1.90b6.4.

Genome accessibility masks from 1KGP were obtained from http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/1000_genomes_project/working/20160622_genome_mask_GRCh38/, and ENCODE blacklisted regions were obtained from https://www.encodeproject.org/files/ENCFF356LFX/.

We note peaks of r² around discrete values of 1, 0.5, 0.33, 0.25, etc. (Figure 1—figure supplement 7), which occur due to the properties of the r² statistic for various allele frequency combinations. Specifically, the maximum possible value of r² for a pair of rare alleles with frequencies p_a and p_b is described by VanLiere and Rosenberg, 2008 as:

r_{m a x}^{2} (p_{a}, p_{b}) = \frac{p_{a} (1 - p_{b})}{(1 - p_{a}) p_{b}}

Assuming a sample size of 198 chromosomes (the most common sample size for a given 1KGP population), the maximum possible value of r² for a pair of singletons is therefore:

p_{a} = 1 / 198

p_{b} = 1 / 198

r_{m a x}^{2} = 1

Meanwhile, the maximum possible value of r² for a singleton-doubleton pair is:

p_{a} = 1 / 198

p_{b} = 2 / 198

r_{m a x}^{2} = 0.497

And the maximum possible value of a singleton-tripleton pair is:

p_{a} = 1 / 198

p_{b} = 3 / 198

r_{m a x}^{2} = 0.330

Given that the most common allele frequency of a SNP or SV is 1/N (i.e., singletons), followed by 2/N (doubletons), 3/N (tripletons), etc., the aforementioned allele frequency combinations occur frequently, resulting in some common discrete values of the maximum r².

eQTL mapping

Request a detailed protocol

To conduct eQTL analysis, we used gene expression data generated by the Geuvadis Consortium (Lappalainen et al., 2013), which includes 447 intersecting samples from the following 1KGP populations: CEU, FIN, GBR, TSI, and YRI. Using the recount3 package from R/Bioconductor (Collado-Torres et al., 2017), we extracted gene expression counts for all corresponding samples. Counts were normalized across samples using the trimmed mean of M values (TMM) method from edgeR (Robinson et al., 2010). TPM values were also computed from raw counts. In accordance with the methods employed for eQTL mapping in the GTEx Project (The GTEx Consortium, 2020), we retained all genes with TPM values greater than or equal to 0.1, as well as raw read counts greater than or equal to six in at least 20% of samples. We then applied rank normalization to the TMM values for each remaining gene. We then performed cis-eQTL mapping with a modified version of fastQTL (https://github.com/hall-lab/fastqtl, Colby, 2016, Ongen et al., 2016), which accounts for SV size when determining the appropriate cis window. We conducted nominal and permutation passes with genotype principal components and sex included as covariates. Beta-distribution-approximated permutation p-values from fastQTL were used as input to estimate q-values and control the false discovery rate (FDR) with the qvalue package (http://github.com/jdstorey/qvalue, Storey et al., 2020, Storey and Tibshirani, 2003).

Next, we used the ChromHMM epigenetic state annotations generated by the Roadmap Epigenomics Project (Roadmap Epigenomics Consortium et al., 2015) for the E116 B lymphocyte cell line to determine whether SVs eQTLs were enriched for regulatory elements. For each of the 15 ChromHMM epigenetic states, we identified the SVs that overlapped with its epigenetic state annotations and determined its enrichment for SV eQTLs using Fisher’s exact test.

To identify potential causal variants for the 1121 significant gene-SV pairs, we performed eQTL fine-mapping using CAVIAR v2.2 (Hormozdiari et al., 2014). Nominal pass output of fastQTL was combined for SNPs and SVs, and p-values were back-converted to z-scores while also accounting for the sign of the association (i.e., positive or negative). For each significant eQTL, we extracted all tested variants within a 1 Mb window, excluding any variants with a nominal p-value > 0.05. We then used PLINK to calculate the LD matrix relating this set of variants. Using these two inputs, we then ran CAVIAR for each of the 1121 significant gene-SV pairs, recording the 90% credible causal set for each eQTL.

Admixture-aware scan for signatures of local adaptation

View detailed protocol

Common examples of frequency differentiation-based metrics to search for selection include Wright’s fixation index (F_ST) (Wright, 1949), as well as tree-based extensions of this concept such as the locus-specific branch length (LSBL) (Shriver et al., 2004) and population branch statistic (PBS) (Yi et al., 2010). While useful for polarizing frequency changes on particular lineages, these tree-based tests still require the specification of (typically three) populations for comparison. The number of possible comparisons thus grows in a combinatorial manner with the number of populations in the study. Specifically, for the 26 populations of 1KGP, there are 2600 (26 choose 3) possible comparisons. A second limitation of such tests is the definition of population, which may or may not reflect genetic patterns of population structure that occur at multiple scales. Moreover, many human populations exhibit substantial admixture, which is ignored by, and may confound, some frequency differentiation-based tests.

We overcame these limitations by using the software package Ohana (Cheng et al., 2019) to scan SVs for signatures of positive selection, that is extreme frequency differentiation between populations. Ohana uses allele frequency to model individuals as an admixed combination of k ancestral populations, constructing a genome-wide covariance matrix to describe the relationship between these ancestry components. Variants are then assessed to determine whether their allele frequencies are better explained by the genome-wide ‘neutral’ matrix, or by an alternative matrix that allows allele frequencies to vary in one ancestry component.

In accordance with Ohana’s recommended workflow, we conducted admixture inference on the 1KGP dataset with a set of ~100,000 variants, downsampled from Chromosome 21 of the 1KGP SNV/indel callset used for LD calculations above. Downsampling was performed with PLINK’s variant pruning function. Inference of admixture proportions in the dataset was allowed to continue until increased iterations produced qualitatively similar results (50 iterations). The covariance matrix generated from this downsampled dataset was used as a neutral input for downstream scans for selection on SVs.

In order to generate ‘selection hypothesis’ matrices to search for selection in a specific ancestry component, we modified the neutral covariance matrix by allowing one component at a time to have a greater covariance. A scalar value of 10, representing the furthest possible deviation that a variant could have in a population under selection, was added to elements of the neutral covariance matrix depending on the population of interest. For each variant of interest, Ohana then computes the likelihood of the observed ancestry-component-specific allele frequencies under the selection and neutral models, then compares them by computing a likelihood ratio statistic (LRS), which quantifies relative support for the selection hypothesis. We note that the ‘selection hypothesis’ matrices we generated for our study enable us to perform tests for frequency differentiation in one ancestry component at a time. Testing for selection in multiple ancestry components would require modifying multiple components in the neutral covariance matrix. While some variants have outlier LRS values in multiple ancestry components, this reflects a combined signature of extremely high allele frequency in one ancestry component and extremely low allele frequency in another (e.g., the IGHG4 insertion we highlight is an outlier in both ancestry component 2, where it reaches high allele frequencies, and ancestry component 6, where it is much rarer than expected from the null model [Figure 3B, panel for ancestry component 6]).

We filtered our results to remove extreme outliers in null model log-likelihoods (global log likelihood [LLE] < −1000), which were unremarkable in their patterns of allele frequency and instead indicated a failure of the neutral model to converge for a small subset of rare variants. To calculate p-values, we then compared the LRS to a chi-square distribution with one degree of freedom (Wojcik et al., 2019) and adjusted for multiple hypothesis testing using a Bonferroni correction. To further refine the list of candidate selected loci, we also compared the ancestry component-specific LRS computed for each SV to the observed distribution of LRS computed from SNPs and short indels matched on global minor allele frequency (in 1% frequency bins). Specifically, we identified SVs with LRS exceeding the 99.9 percentile of the empirical LRS distribution for frequency-matched SNPs and short indels as calculated using identical methods. These background SNPs and indels were limited to Chromosome 1 for computational efficiency, but results were qualitatively unaffected by the choice of chromosome.

Assessment of archaic introgression

Request a detailed protocol

To identify candidate archaic introgressed SVs among the set of highly differentiated variants, we tested each SV for LD with putative introgressed SNPs as identified based on published results from the method Sprime, which is based on signatures of LD and divergence from an African outgroup population (Browning et al., 2018). Specifically, we computed pairwise LD between the SV and all SNPs in a 100 kb window, matching the ancestry component to a corresponding population from 1KGP. We additionally computed the allele frequencies of SVs in African populations, excluding the ASW (Americans of African Ancestry in SW USA) and ACB (African Caribbeans in Barbados) populations which exhibit substantial non-African admixture. We reported all highly differentiated SVs possessing r² > 0.5 with any putative introgressed SNP and AF < 0.01 in non-admixed African populations (Supplementary file 2).

We obtained masks for the Altai Neanderthal sequence, which were used in filtering by Browning et al., 2018, from http://cdna.eva.mpg.de/neandertal/Vindija/FilterBed/. These masked regions filter for coverage (stratified by GC content), require a minimum coverage of 10×, require all 35-mers in the region to be unique, and exclude tandem repeats and indels.

To efficiently search for the IGHG4 insertion in additional datasets, we designed a 48 bp sequence (TGGAGAGAGTGGGGGACAGCGTCAGGGACAGGTGGGGACAGCCTGGGG) that spans the insertion breakpoint, extending across the entire 33 bp of the insertion itself as well as 11 bp and 4 bp, into the respective upstream and downstream flanking regions. BLAST searches for this sequence in the hg19 and GRCh38 human reference genomes returned no exact matches, but the diagnostic sequence was sufficiently similar to the reference sequence that reads containing it still mapped to the same locus. Sequence alignments from four high-coverage archaic hominin samples (Altai Neanderthal, Vindija Neanderthal, Chagyrskaya Neanderthal, and Denisovan) were obtained from http://cdna.eva.mpg.de/neandertal and http://cdna.eva.mpg.de/denisova. Forward strand sequences of unique reads were extracted from reads aligned to the hg19 reference genome, and exact matches to the 48 bp query sequence were identified using grep.

Evidence of introgression at the IGH locus was further examined by counting observed alleles at highly differentiated SNPs from sequenced alignments for each high-coverage archaic sample (see above). Sites with two or more reads supporting the alternative allele were used to define matching and color Figure 4. further conditions on alternative allele frequency ≤ 1% within African populations of 1KGP.

Inference of selection parameters with approximate Bayesian computation

Request a detailed protocol

We used a sequential algorithm for approximate Bayesian computation (Lenormand et al., 2013; Pritchard et al., 1999), implemented with the R package EasyABC (Jabot et al., 2013), to infer the strength and timing of selection at the IGHG4 locus, as well as the initial frequency of the adaptive allele. This approach consisted of drawing model parameters from prior distributions as input to the forward evolutionary simulation software package SLiM (Haller and Messer, 2019), computing summary statistics from each simulation, and comparing to those observed from our data. Simulations with summary statistics most closely matching the observed data are then used to construct posterior distributions of the model parameters. The sequential algorithm automatically determines the tolerance level and uses a predetermined stopping criterion (p_accmin = 0.05), thereby reducing the necessary number of simulations and improving estimates of the posterior distribution (Lenormand et al., 2013).

We constructed a simplified five-population demographic model based on parameters obtained from Jouganous et al., 2017 and Wang et al., 2018, as well as population size estimates from 1KGP (The 1000 Genomes Project Consortium, 2015; Figure 5C). Specifically, our model consisted of an initial divergence event between the CEU and East Asian populations at 46 kya, subsequent divergence between the population ancestral to CHB/JPT and the population ancestral to KHV/CHX at 9.8 kya, and final splits between CHB/JPT and KHV/CHX at 9.0 kya and 1.7 kya, respectively. Each new population was drawn from its originating population at the same size of the originating population. The size of the ancestral population was set at 2831 (Jouganous et al., 2017) and allowed to expand exponentially at a rate of 1.25 × 10^–3 per generation, resulting in a final size of 17,883 in each subpopulation, broadly consistent with pairwise sequential Markovian coalescent (PSMC) results presented by 1KGP for these populations (The 1000 Genomes Project Consortium, 2015). We allowed the initial allele frequency of the selected variant (p₀) to vary between 0 and 1, the timing of the onset of selection (T_adaptive) to vary between 1 and 1686 generations ago (i.e., the entire simulated timespan), and the selection coefficient (s) to vary between –0.01 and 0.2 – all drawn from uniform distributions. We used the selected variant’s end-of-simulation allele frequencies as a summary statistic for the model, comparing these frequencies to the observed allele frequencies of the IGHG4 insertion in these five populations.

To obtain a background distribution of the expected frequencies of the introgressed haplotype in the absence of selection, we performed 100,000 of the simulations described above with a selection coefficient of 0. For ease of comparison to the observed data, we used the end-of-simulation population branch statistic (PBS) (Yi et al., 2010), calculated between the CDX, JPT, and CEU populations, as a summary statistic. We compared these PBS values, as well as the inter-population F_ST, to those observed for the IGHG4 insertion in these populations.

Phenotype-wide association analysis

Request a detailed protocol

We examined potential phenotype associations with the putative Neanderthal-introgressed haplotype at the IGH locus by extracting summary statistics from the pan-ancestry analysis of the UK Biobank (https://pan.ukbb.broadinstitute.org/). Specifically, we obtained association p-values for two SNPs, which each tag one of the two major LD blocks (rs115091999 and rs150526114). We restricted analysis to individuals of East Asian ancestry. No variants were significant after Bonferroni correction (Figure 5—figure supplement 2).

Data availability

All code necessary for reproducing our analysis is available on GitHub (https://github.com/mccoy-lab/sv_selection; copy archivedat https://archive.softwareheritage.org/swh:1:rev:2e105bf050fd2acb9f1e0d1833566c0b1ca5b41f). SV genotypes, eQTL results, and selection scan results are available on Zenodo (https://doi.org/10.5281/zenodo.5509980).

The following data sets were generated

1. Yan SM
2. McCoy RC
(2021) Zenodo
Data from: Local adaptation and archaic introgression shape global diversity at human structural variant loci.

https://doi.org/10.5281/zenodo.5509980

The following previously published data sets were used

1. Audano PA
(2019) NCBI BioProject
ID PRJNA300840. Long-read sequencing of a Puerto Rican individual (HG00733).

https://www.ncbi.nlm.nih.gov/bioproject/300840
1. Audano PA
(2019) NCBI BioProject
ID PRJNA288807. Long-read sequencing of a Yoruban individual (NA19240).

https://www.ncbi.nlm.nih.gov/bioproject/288807
1. Audano PA
(2019) NCBI BioProject
ID PRJNA339722. Long-read sequencing of a Gambian individual (HG02818).

https://www.ncbi.nlm.nih.gov/bioproject/339722
1. Audano PA
(2019) NCBI BioProject
ID PRJNA385272. Long-read sequencing of a Luhya individual (NA19434).

https://www.ncbi.nlm.nih.gov/bioproject/385272
1. Audano PA
(2019) NCBI BioProject
ID PRJNA339726. Long-read sequencing of a Vietnamese individual (HG02059).

https://www.ncbi.nlm.nih.gov/bioproject/339726
1. Audano PA
(2019) NCBI BioProject
ID PRJNA323611. Long-read sequencing of a Utah individual with Northern and Western European ancestry (NA12878).

https://www.ncbi.nlm.nih.gov/bioproject/323611
1. Audano PA
(2019) NCBI BioProject
ID PRJNA481794. Long-read sequencing of a Telugu individual (HG04217).

https://www.ncbi.nlm.nih.gov/bioproject/481794
1. Audano PA
(2019) NCBI BioProject
ID PRJNA480858. Long-read sequencing of a Peruvian individual (HG02106).

https://www.ncbi.nlm.nih.gov/bioproject/480858
1. Audano PA
(2019) NCBI BioProject
ID PRJNA480712. Long-read sequencing of a Finnish individual (HG00268).

https://www.ncbi.nlm.nih.gov/bioproject/480712
1. Audano PA
(2019) NCBI BioProject
ID PRJNA300843. Long-read sequencing of a Han Chinese individual (HG00514).

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA300843
1. Byrska-Bishop M
(2019) NCBI BioProject
ID PRJEB31736. 1000 Genomes Project phase 3: 30X coverage whole genome sequencing.

https://www.ncbi.nlm.nih.gov/bioproject/527456
1. Huddleston J
(2017) NCBI BioProject
ID PRJNA246220. Long-read sequencing of a hydatidiform mole (CHM1).

https://www.ncbi.nlm.nih.gov/bioproject/246220
1. Huddleston J
(2017) NCBI BioProject
ID PRJNA269593. Long-read sequencing of a hydatidiform mole (CHM13).

https://www.ncbi.nlm.nih.gov/bioproject/269593
1. Lappalainen T
(2013) ArrayExpress
ID E-GEUV-1. RNA-sequencing of 465 lymphoblastoid cell lines from the 1000 Genomes.

https://www.ebi.ac.uk/arrayexpress/experiments/E-GEUV-1/
1. Seo J-S
(2016) NCBI BioProject
ID PRJNA298944. Long-read sequencing of a Korean individual (AK1).

https://www.ncbi.nlm.nih.gov/bioproject/298944
1. Shi L
(2016) NCBI BioProject
ID PRJNA301527. Long-read sequencing of a Chinese individual (HX1).

https://www.ncbi.nlm.nih.gov/bioproject/301527

References

1. Abi-Rached L
2. Jobin MJ
3. Kulkarni S
4. McWhinnie A
5. Dalva K
6. Gragert L
7. Babrzadeh F
8. Gharizadeh B
9. Luo M
10. Plummer FA
11. Kimani J
12. Carrington M
13. Middleton D
14. Rajalingam R
15. Beksac M
16. Marsh SGE
17. Maiers M
18. Guethlein LA
19. Tavoularis S
20. Little AM
21. Green RE
22. Norman PJ
23. Parham P
(2011) The shaping of modern human immune systems by multiregional admixture with archaic humans
Science 334:89–94.

https://doi.org/10.1126/science.1209202
- PubMed
- Google Scholar
1. Aganezov S
2. Goodwin S
3. Sherman RM
4. Sedlazeck FJ
5. Arun G
6. Bhatia S
7. Lee I
8. Kirsche M
9. Wappel R
10. Kramer M
11. Kostroff K
12. Spector DL
13. Timp W
14. McCombie WR
15. Schatz MC
(2020) Comprehensive analysis of structural variants in breast cancer genomes using single-molecule sequencing
Genome Research 30:1258–1273.

https://doi.org/10.1101/gr.260497.119
- PubMed
- Google Scholar
Preprint
1. Aganezov S
2. Yan SM
3. Soto DC
4. Kirsche M
5. Zarate S
6. Avdeyev P
7. Taylor DJ
8. Shafin K
9. Shumate A
10. Xiao C
11. Wagner J
12. McDaniel J
13. Olson ND
14. Sauria MEG
15. Vollger MR
16. Rhie A
17. Meredith M
18. Martin S
19. Lee J
20. Koren S
21. Rosenfeld JA
22. Paten B
23. Layer R
24. Chin CS
25. Sedlazeck FJ
26. Hansen NF
27. Miller DE
28. Phillippy AM
29. Miga KH
30. McCoy RC
31. Dennis MY
32. Zook JM
33. Schatz MC
(2021) A Complete Reference Genome Improves Analysis of Human Genetic Variation
bioRxiv.

https://doi.org/10.1101/2021.07.12.452063
- Google Scholar
(2009) Fast model-based estimation of ancestry in unrelated individuals
Genome Research 19:1655–1664.

https://doi.org/10.1101/gr.094052.109
- PubMed
- Google Scholar
1. Almarri MA
2. Bergström A
3. Prado-Martinez J
4. Yang F
5. Fu B
6. Dunham AS
7. Chen Y
8. Hurles ME
9. Tyler-Smith C
10. Xue Y
(2020) Population structure, stratification, and introgression of human structural variation
Cell 182:189–199.

https://doi.org/10.1016/j.cell.2020.05.024
- PubMed
- Google Scholar
1. Alonge M
2. Wang X
3. Benoit M
4. Soyk S
5. Pereira L
6. Zhang L
7. Suresh H
8. Ramakrishnan S
9. Maumus F
10. Ciren D
11. Levy Y
12. Harel TH
13. Shalev-Schlosser G
14. Amsellem Z
15. Razifard H
16. Caicedo AL
17. Tieman DM
18. Klee H
19. Kirsche M
20. Aganezov S
21. Ranallo-Benavidez TR
22. Lemmon ZH
23. Kim J
24. Robitaille G
25. Kramer M
26. Goodwin S
27. McCombie WR
28. Hutton S
29. Van Eck J
30. Gillis J
31. Eshed Y
32. Sedlazeck FJ
33. van der Knaap E
34. Schatz MC
35. Lippman ZB
(2020) Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato
Cell 182:145–161.

https://doi.org/10.1016/j.cell.2020.05.021
- PubMed
- Google Scholar
(2019) The ENCODE Blacklist: Identification of Problematic Regions of the Genome
Scientific Reports 9:9354.

https://doi.org/10.1038/s41598-019-45839-z
- PubMed
- Google Scholar
1. Audano PA
2. Sulovari A
3. Graves-Lindsay TA
4. Cantsilieris S
5. Sorensen M
6. Welch AE
7. Dougherty ML
8. Nelson BJ
9. Shah A
10. Dutcher SK
11. Warren WC
12. Magrini V
13. McGrath SD
14. Yi L
15. Wilson RK
16. Eichler EE
(2019) Characterizing the Major Structural Variant Alleles of the Human Genome
Cell 176:663–675.

https://doi.org/10.1016/j.cell.2018.12.019
- PubMed
- Google Scholar
1. Barreiro LB
2. Quintana-Murci L
(2020) Evolutionary and population (epi)genetics of immunity to infection
Human Genetics 139:723–732.

https://doi.org/10.1007/s00439-020-02167-x
- PubMed
- Google Scholar
1. Bergström A
2. McCarthy SA
3. Hui R
4. Almarri MA
5. Ayub Q
6. Danecek P
7. Chen Y
8. Felkel S
9. Hallast P
10. Kamm J
11. Blanché H
12. Deleuze JF
13. Cann H
14. Mallick S
15. Reich D
16. Sandhu MS
17. Skoglund P
18. Scally A
19. Xue Y
20. Durbin R
21. Tyler-Smith C
(2020) Insights into human genetic variation and population history from 929 diverse genomes
Science 367:eaay5012.

https://doi.org/10.1126/science.aay5012
- PubMed
- Google Scholar
(2004) Genetic Signatures of Strong Recent Positive Selection at the Lactase Gene
The American Journal of Human Genetics 74:1111–1120.

https://doi.org/10.1086/421051
- PubMed
- Google Scholar
(2021) Long-read sequencing of 3,622 Icelanders provides insight into the role of structural variants in human diseases and other traits
Nature Genetics 53:779–786.

https://doi.org/10.1038/s41588-021-00865-4
- PubMed
- Google Scholar
(2008) Estimation of 2Nes from temporal allele frequency data
Genetics 179:497–502.

https://doi.org/10.1534/genetics.107.085019
- PubMed
- Google Scholar
1. Browning SR
2. Browning BL
3. Zhou Y
4. Tucci S
5. Akey JM
(2018) Analysis of Human Sequence Data Reveals Two Pulses of Archaic Denisovan Admixture
Cell 173:53–61.

https://doi.org/10.1016/j.cell.2018.02.031
- PubMed
- Google Scholar
Preprint
1. Byrska-Bishop M
2. Evani US
3. Zhao X
4. Basile AO
5. Abel HJ
6. Regier AA
7. Corvelo A
8. Clarke WE
9. Musunuri R
10. Nagulapalli K
11. Fairley S
12. Runnels A
13. Winterkorn L
14. Lowy-Gallego E
15. Flicek P
16. Germer S
17. Brand H
18. Hall IM
19. Talkowski ME
20. Narzisi G
21. Zody MC
(2021) High Coverage Whole Genome Sequencing of the Expanded 1000 Genomes Project Cohort Including 602 Trios
bioRxiv.

https://doi.org/10.1101/2021.02.06.430068
- Google Scholar
(2015) Resolving the complexity of the human genome using single-molecule sequencing
Nature 517:608–611.

https://doi.org/10.1038/nature13907
- PubMed
- Google Scholar
1. Chaisson MJP
2. Sanders AD
3. Zhao X
4. Malhotra A
5. Porubsky D
6. Rausch T
7. Gardner EJ
8. Rodriguez OL
9. Guo L
10. Collins RL
11. Fan X
12. Wen J
13. Handsaker RE
14. Fairley S
15. Kronenberg ZN
16. Kong X
17. Hormozdiari F
18. Lee D
19. Wenger AM
20. Hastie AR
21. Antaki D
22. Anantharaman T
23. Audano PA
24. Brand H
25. Cantsilieris S
26. Cao H
27. Cerveira E
28. Chen C
29. Chen X
30. Chin CS
31. Chong Z
32. Chuang NT
33. Lambert CC
34. Church DM
35. Clarke L
36. Farrell A
37. Flores J
38. Galeev T
39. Gorkin DU
40. Gujral M
41. Guryev V
42. Heaton WH
43. Korlach J
44. Kumar S
45. Kwon JY
46. Lam ET
47. Lee JE
48. Lee J
49. Lee WP
50. Lee SP
51. Li S
52. Marks P
53. Viaud-Martinez K
54. Meiers S
55. Munson KM
56. Navarro FCP
57. Nelson BJ
58. Nodzak C
59. Noor A
60. Kyriazopoulou-Panagiotopoulou S
61. Pang AWC
62. Qiu Y
63. Rosanio G
64. Ryan M
65. Stütz A
66. Spierings DCJ
67. Ward A
68. Welch AE
69. Xiao M
70. Xu W
71. Zhang C
72. Zhu Q
73. Zheng-Bradley X
74. Lowy E
75. Yakneen S
76. McCarroll S
77. Jun G
78. Ding L
79. Koh CL
80. Ren B
81. Flicek P
82. Chen K
83. Gerstein MB
84. Kwok PY
85. Lansdorp PM
86. Marth GT
87. Sebat J
88. Shi X
89. Bashir A
90. Ye K
91. Devine SE
92. Talkowski ME
93. Mills RE
94. Marschall T
95. Korbel JO
96. Eichler EE
97. Lee C
(2019) Multi-platform discovery of haplotype-resolved structural variation in human genomes
Nature Communications 10:1784.

https://doi.org/10.1038/s41467-018-08148-z
- PubMed
- Google Scholar
1. Chen S
2. Krusche P
3. Dolzhenko E
4. Sherman RM
5. Petrovski R
6. Schlesinger F
7. Kirsche M
8. Bentley DR
9. Schatz MC
10. Sedlazeck FJ
11. Eberle MA
(2019) Paragraph: a graph-based structural variant genotyper for short-read sequence data
Genome Biology 20:291.

https://doi.org/10.1186/s13059-019-1909-7
- PubMed
- Google Scholar
(2017) Fast admixture analysis and population tree estimation for SNP and NGS data
Bioinformatics 33:2148–2155.

https://doi.org/10.1093/bioinformatics/btx098
- PubMed
- Google Scholar
Preprint
(2019) Ohana: Detecting Selection in Multiple Populations by Modelling Ancestral Admixture Components
bioRxiv.

https://doi.org/10.1101/546408
- Google Scholar
1. Chiang C
2. Scott AJ
3. Davis JR
4. Tsang EK
5. Li X
6. Kim Y
7. Hadzic T
8. Damani FN
9. Ganel L
10. GTEx Consortium
11. Montgomery SB
12. Battle A
13. Conrad DF
14. Hall IM
(2017) The impact of structural variation on human gene expression
Nature Genetics 49:692–699.

https://doi.org/10.1038/ng.3834
- PubMed
- Google Scholar
1. Chin CS
2. Wagner J
3. Zeng Q
4. Garrison E
5. Garg S
6. Fungtammasan A
7. Rautiainen M
8. Aganezov S
9. Kirsche M
10. Zarate S
11. Schatz MC
12. Xiao C
13. Rowell WJ
14. Markello C
15. Farek J
16. Sedlazeck FJ
17. Bansal V
18. Yoo B
19. Miller N
20. Zhou X
21. Carroll A
22. Barrio AM
23. Salit M
24. Marschall T
25. Dilthey AT
26. Zook JM
(2020) A diploid assembly-based benchmark for variants in the major histocompatibility complex
Nature Communications 11:4794.

https://doi.org/10.1038/s41467-020-18564-9
- PubMed
- Google Scholar
Software
1. Colby C
(2016) FastQTL, version 4d03819
GitHub.

https://github.com/hall-lab/fastqtl
1. Collado-Torres L
2. Nellore A
3. Kammers K
4. Ellis SE
5. Taub MA
6. Hansen KD
7. Jaffe AE
8. Langmead B
9. Leek JT
(2017) Reproducible RNA-seq analysis using recount2
Nature Biotechnology 35:319–321.

https://doi.org/10.1038/nbt.3838
- PubMed
- Google Scholar
1. Conrad DF
2. Jakobsson M
3. Coop G
4. Wen X
5. Wall JD
6. Rosenberg NA
7. Pritchard JK
(2006) A worldwide survey of haplotype variation and linkage disequilibrium in the human genome
Nature Genetics 38:1251–1260.

https://doi.org/10.1038/ng1911
- PubMed
- Google Scholar
1. Crawford NG
2. Kelly DE
3. Hansen MEB
4. Beltrame MH
5. Fan S
6. Bowman SL
7. Jewett E
8. Ranciaro A
9. Thompson S
10. Lo Y
11. Pfeifer SP
12. Jensen JD
13. Campbell MC
14. Beggs W
15. Hormozdiari F
16. Mpoloka SW
17. Mokone GG
18. Nyambo T
19. Meskel DW
20. Belay G
21. Haut J
22. Rothschild H
23. Zon L
24. Zhou Y
25. Kovacs MA
26. Xu M
27. Zhang T
28. Bishop K
29. Sinclair J
30. Rivas C
31. Elliot E
32. Choi J
33. Sa L
34. Hicks B
35. Burgess S
36. Abnet C
37. Watkins-Chow DE
38. Oceana E
39. Song YS
40. Eskin E
41. Brown KM
42. Marks MS
43. Loftus SK
44. Pavan WJ
45. Yeager M
46. Chanock S
47. Tishkoff SA
(2017) Loci associated with skin pigmentation identified in African populations
Science 358:eaan8433.

https://doi.org/10.1126/science.aan8433
- PubMed
- Google Scholar
1. Danecek P
2. Bonfield JK
3. Liddle J
4. Marshall J
5. Ohan V
6. Pollard MO
7. Whitwham A
8. Keane T
9. McCarthy SA
10. Davies RM
11. Li H
(2021) Twelve years of SAMtools and BCFtools
GigaScience 10:giab008.

https://doi.org/10.1093/gigascience/giab008
- PubMed
- Google Scholar
(2016) Introgression of Neandertal- and Denisovan-like Haplotypes Contributes to Adaptive Variation in Human Toll-like Receptors
American Journal of Human Genetics 98:22–23.

https://doi.org/10.1016/j.ajhg.2015.11.015
- PubMed
- Google Scholar
1. Dean M
2. Carrington M
3. Winkler C
4. Huttley GA
5. Smith MW
6. Allikmets R
7. Goedert JJ
8. Buchbinder SP
9. Vittinghoff E
10. Gomperts E
11. Donfield S
12. Vlahov D
13. Kaslow R
14. Saah A
15. Rinaldo C
16. Detels R
17. O’Brien SJ
(1996) Genetic restriction of HIV-1 infection and progression to AIDS by a deletion allele of the CKR5 structural gene Hemophilia Growth and Development Study, Multicenter AIDS Cohort Study, Multicenter Hemophilia Cohort Study, San Francisco City Cohort, ALIVE Study
Science 273:1856–1862.

https://doi.org/10.1126/science.273.5283.1856
- PubMed
- Google Scholar
1. Ding J
2. Wu W
3. Yang J
4. Wu M
(2019) Long non-coding RNA MIF-AS1 promotes breast cancer cell proliferation, migration and EMT process through regulating miR-1249-3p/HOXB8 axis
Pathology, Research and Practice 215:005.

https://doi.org/10.1016/j.prp.2019.03.005
- PubMed
- Google Scholar
1. Ebert P
2. Audano PA
3. Zhu Q
4. Rodriguez-Martin B
5. Porubsky D
6. Bonder MJ
7. Sulovari A
8. Ebler J
9. Zhou W
10. Serra Mari R
11. Yilmaz F
12. Zhao X
13. Hsieh P
14. Lee J
15. Kumar S
16. Lin J
17. Rausch T
18. Chen Y
19. Ren J
20. Santamarina M
21. Höps W
22. Ashraf H
23. Chuang NT
24. Yang X
25. Munson KM
26. Lewis AP
27. Fairley S
28. Tallon LJ
29. Clarke WE
30. Basile AO
31. Byrska-Bishop M
32. Corvelo A
33. Evani US
34. T-Y L
35. Chaisson MJP
36. Chen J
37. Li C
38. Brand H
39. Wenger AM
40. Ghareghani M
41. Harvey WT
42. Raeder B
43. Hasenfeld P
44. Regier AA
45. Abel HJ
46. Hall IM
47. Flicek P
48. Stegle O
49. Gerstein MB
50. Tubio JMC
51. Mu Z
52. Yi L
53. Shi X
54. Hastie AR
55. Ye K
56. Chong Z
57. Sanders AD
58. Zody MC
59. Talkowski ME
60. Mills RE
61. Devine SE
62. Lee C
63. Korbel JO
64. Marschall T
65. Eichler EE
(2021) Haplotype-resolved diverse human genomes and integrated analysis of structural variation
Science 372:eabf7117.

https://doi.org/10.1126/science.abf7117
- PubMed
- Google Scholar
1. Elguero E
2. Délicat-Loembet LM
3. Rougeron V
4. Arnathau C
5. Roche B
6. Becquart P
7. Gonzalez JP
8. Nkoghe D
9. Sica L
10. Leroy EM
11. Durand P
12. Ayala FJ
13. Ollomo B
14. Renaud F
15. Prugnolle F
(2015) Malaria continues to select for sickle cell trait in Central Africa
PNAS 112:7051–7054.

https://doi.org/10.1073/pnas.1505665112
- PubMed
- Google Scholar
1. Enard D
2. Petrov DA
(2018) Evidence that RNA Viruses Drove Adaptive Introgression between Neanderthals and Modern Humans
Cell 175:360–371.

https://doi.org/10.1016/j.cell.2018.08.034
- PubMed
- Google Scholar
(2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies
Genetics 164:1567–1587.

https://doi.org/10.1093/genetics/164.4.1567
- PubMed
- Google Scholar
1. Fan S
2. Hansen MEB
3. Lo Y
4. Tishkoff SA
(2016) Going global by adapting local: A review of recent human adaptation
Science 354:54–59.

https://doi.org/10.1126/science.aaf5098
- PubMed
- Google Scholar
1. Galvani AP
2. Slatkin M
(2003) Evaluating plague and smallpox as historical selective pressures for the CCR5-Delta 32 HIV-resistance allele
PNAS 100:15276–15279.

https://doi.org/10.1073/pnas.2435085100
- PubMed
- Google Scholar
(2016) Archaic Hominin Admixture Facilitated Adaptation to Out-of-Africa Environments
Current Biology 26:3375–3382.

https://doi.org/10.1016/j.cub.2016.10.041
- PubMed
- Google Scholar
1. Gouy A
2. Excoffier L
(2020) Polygenic Patterns of Adaptive Introgression in Modern Humans Are Mainly Shaped by Response to Pathogens
Molecular Biology and Evolution 37:1420–1433.

https://doi.org/10.1093/molbev/msz306
- PubMed
- Google Scholar
1. Graffelman J
(2015) Exploring Diallelic Genetic Markers: The HardyWeinberg Package
Journal of Statistical Software 64:i03.

https://doi.org/10.18637/jss.v064.i03
- Google Scholar
(2017) A genome-wide study of Hardy-Weinberg equilibrium with next generation sequence data
Human Genetics 136:727–741.

https://doi.org/10.1007/s00439-017-1786-7
- PubMed
- Google Scholar
1. Green RE
2. Krause J
3. Briggs AW
4. Maricic T
5. Stenzel U
6. Kircher M
7. Patterson N
8. Li H
9. Zhai W
10. Fritz MH-Y
11. Hansen NF
12. Durand EY
13. Malaspinas A-S
14. Jensen JD
15. Marques-Bonet T
16. Alkan C
17. Prüfer K
18. Meyer M
19. Burbano HA
20. Good JM
21. Schultz R
22. Aximu-Petri A
23. Butthof A
24. Höber B
25. Höffner B
26. Siegemund M
27. Weihmann A
28. Nusbaum C
29. Lander ES
30. Russ C
31. Novod N
32. Affourtit J
33. Egholm M
34. Verna C
35. Rudan P
36. Brajkovic D
37. Kucan Ž
38. Gušic I
39. Doronichev VB
40. Golovanova LV
41. Lalueza-Fox C
42. de la Rasilla M
43. Fortea J
44. Rosas A
45. Schmitz RW
46. Johnson PLF
47. Eichler EE
48. Falush D
49. Birney E
50. Mullikin JC
51. Slatkin M
52. Nielsen R
53. Kelso J
54. Lachmann M
55. Reich D
56. Pääbo S
(2010) A draft sequence of the Neandertal genome
Science 328:710–722.

https://doi.org/10.1126/science.1188021
- PubMed
- Google Scholar
1. Haller BC
2. Messer PW
(2019) SLiM 3: Forward Genetic Simulations Beyond the Wright–Fisher Model
Molecular Biology and Evolution 36:632–637.

https://doi.org/10.1093/molbev/msy228
- PubMed
- Google Scholar
(2021) Rapid adaptation to malaria facilitated by admixture in the human population of Cabo Verde
eLife 10:e63177.

https://doi.org/10.7554/eLife.63177
- PubMed
- Google Scholar
1. Hickey G
2. Heller D
3. Monlong J
4. Sibbesen JA
5. Sirén J
6. Eizenga J
7. Dawson ET
8. Garrison E
9. Novak AM
10. Paten B
(2020) Genotyping structural variants in pangenome graphs using the vg toolkit
Genome Biology 21:35.

https://doi.org/10.1186/s13059-020-1941-7
- PubMed
- Google Scholar
1. Hormozdiari F
2. Kostem E
3. Kang EY
4. Pasaniuc B
5. Eskin E
(2014) Identifying Causal Variants at Loci with Multiple Signals of Association
Genetics 198:497–508.

https://doi.org/10.1534/genetics.114.167908
- PubMed
- Google Scholar
1. Hsieh P
2. Vollger MR
3. Dang V
4. Porubsky D
5. Baker C
6. Cantsilieris S
7. Hoekzema K
8. Lewis AP
9. Munson KM
10. Sorensen M
11. Kronenberg ZN
12. Murali S
13. Nelson BJ
14. Chiatante G
15. Maggiolini FAM
16. Blanché H
17. Underwood JG
18. Antonacci F
19. Deleuze JF
20. Eichler EE
(2019) Adaptive archaic introgression of copy number variants and the discovery of previously unknown human genes
Science 366:eaax2083.

https://doi.org/10.1126/science.aax2083
- PubMed
- Google Scholar
1. Huerta-Sánchez E
2. Jin X
3. Asan B
4. Peter Z
5. Vinckenbosch BM
6. Liang N
7. Yi Y
8. He X
9. Somel M
10. Ni M
11. Wang P
12. Ou B
13. Huasang X
14. Luosang J
15. Zxp C
16. Li K
17. Gao G
18. Yin Y
19. Wang W
20. Zhang X
21. Xu X
22. Yang H
23. Li Y
24. Jian W
25. Jun W
26. Nielsen R
(2014) Altitude adaptation in Tibetans caused by introgression of Denisovan-like DNA
Nature 512:194–197.

https://doi.org/10.1038/nature13408
- PubMed
- Google Scholar
(2018) Physiological and Genetic Adaptations to Diving in Sea Nomads
Cell 173:569–580.

https://doi.org/10.1016/j.cell.2018.03.054
- PubMed
- Google Scholar
1. Jablonski NG
(2004) The Evolution of Human Skin and Skin Color
Annual Review of Anthropology 33:585–623.

https://doi.org/10.1146/annurev.anthro.33.070203.143955
- Google Scholar
(2013) EasyABC: Performing efficient approximate Bayesian computation sampling schemes using R
Methods in Ecology and Evolution 10:12050.

https://doi.org/10.1111/2041-210x.12050
- Google Scholar
(2020a) Properties of structural variants and short tandem repeats associated with gene expression and complex traits
Nature Communications 11:2927.

https://doi.org/10.1038/s41467-020-16482-4
- Google Scholar
(2020b) Discovery and quality analysis of a comprehensive set of structural variants and short tandem repeats
Nature Communications 11:2928.

https://doi.org/10.1038/s41467-020-16481-5
- Google Scholar
1. Jeffares DC
2. Jolly C
3. Hoti M
4. Speed D
5. Shaw L
6. Rallis C
7. Balloux F
8. Dessimoz C
9. Bähler J
10. Sedlazeck FJ
(2017) Transient structural variations have strong effects on quantitative traits and reproductive isolation in fission yeast
Nature Communications 8:14061.

https://doi.org/10.1038/ncomms14061
- PubMed
- Google Scholar
(2017) Inferring the Joint Demographic History of Multiple Populations: Beyond the Diffusion Approximation
Genetics 206:1549–1567.

https://doi.org/10.1534/genetics.117.200493
- PubMed
- Google Scholar
1. Karczewski KJ
2. Francioli LC
3. Tiao G
4. Cummings BB
5. Alföldi J
6. Wang Q
7. Collins RL
8. Laricchia KM
9. Ganna A
10. Birnbaum DP
11. Gauthier LD
12. Brand H
13. Solomonson M
14. Watts NA
15. Rhodes D
16. Singer-Berk M
17. England EM
18. Seaby EG
19. Kosmicki JA
20. Walters RK
21. Tashman K
22. Farjoun Y
23. Banks E
24. Poterba T
25. Wang A
26. Seed C
27. Whiffin N
28. Chong JX
29. Samocha KE
30. Pierce-Hoffman E
31. Zappala Z
32. O’Donnell-Luria AH
33. Minikel EV
34. Weisburd B
35. Lek M
36. Ware JS
37. Vittal C
38. Armean IM
39. Bergelson L
40. Cibulskis K
41. Connolly KM
42. Covarrubias M
43. Donnelly S
44. Ferriera S
45. Gabriel S
46. Gentry J
47. Gupta N
48. Jeandet T
49. Kaplan D
50. Llanwarne C
51. Munshi R
52. Novod S
53. Petrillo N
54. Roazen D
55. Ruano-Rubio V
56. Saltzman A
57. Schleicher M
58. Soto J
59. Tibbetts K
60. Tolonen C
61. Wade G
62. Talkowski ME
63. Neale BM
64. Daly MJ
65. MacArthur DG
(2020) The mutational constraint spectrum quantified from variation in 141,456 humans
Nature 581:434–443.

https://doi.org/10.1038/s41586-020-2308-7
- PubMed
- Google Scholar
1. Keinan A
2. Clark AG
(2012) Recent Explosive Human Population Growth Has Resulted in an Excess of Rare Genetic Variants
Science 336:740–743.

https://doi.org/10.1126/science.1217283
- PubMed
- Google Scholar
1. Kosugi S
2. Momozawa Y
3. Liu X
4. Terao C
5. Kubo M
6. Kamatani Y
(2019) Comprehensive evaluation of structural variation detection algorithms for whole genome sequencing
Genome Biology 20:117.

https://doi.org/10.1186/s13059-019-1720-5
- PubMed
- Google Scholar
1. Kothapalli KSD
2. Ye K
3. Gadgil MS
4. Carlson SE
5. O’Brien KO
6. Zhang JY
7. Park HG
8. Ojukwu K
9. Zou J
10. Hyon SS
11. Joshi KS
12. Gu Z
13. Keinan A
14. Brenna JT
(2016) Positive Selection on a Regulatory Insertion-Deletion Polymorphism in FADS2 Influences Apparent Endogenous Synthesis of Arachidonic Acid
Molecular Biology and Evolution 33:1726–1739.

https://doi.org/10.1093/molbev/msw049
- PubMed
- Google Scholar
1. Lachance J
2. Tishkoff SA
(2013) SNP ascertainment bias in population genetic analyses: Why it is important, and how to correct it
BioEssays 35:780–786.

https://doi.org/10.1002/bies.201300014
- PubMed
- Google Scholar
1. Lappalainen T
2. Sammeth M
3. Friedländer MR
4. ’t Hoen PAC
5. Monlong J
6. Rivas MA
7. Gonzàlez-Porta M
8. Kurbatova N
9. Griebel T
10. Ferreira PG
11. Barann M
12. Wieland T
13. Greger L
14. van Iterson M
15. Almlöf J
16. Ribeca P
17. Pulyakhina I
18. Esser D
19. Giger T
20. Tikhonov A
21. Sultan M
22. Bertier G
23. MacArthur DG
24. Lek M
25. Lizano E
26. Buermans HPJ
27. Padioleau I
28. Schwarzmayr T
29. Karlberg O
30. Ongen H
31. Kilpinen H
32. Beltran S
33. Gut M
34. Kahlem K
35. Amstislavskiy V
36. Stegle O
37. Pirinen M
38. Montgomery SB
39. Donnelly P
40. McCarthy MI
41. Flicek P
42. Strom TM
43. Geuvadis Consortium
44. Lehrach H
45. Schreiber S
46. Sudbrak R
47. Carracedo A
48. Antonarakis SE
49. Häsler R
50. Syvänen A-C
51. van Ommen G-J
52. Brazma A
53. Meitinger T
54. Rosenstiel P
55. Guigó R
56. Gut IG
57. Estivill X
58. Dermitzakis ET
(2013) Transcriptome and genome sequencing uncovers functional variation in humans
Nature 501:506–511.

https://doi.org/10.1038/nature12531
- PubMed
- Google Scholar
(2013) Adaptive approximate Bayesian computation for complex models
Computational Statistics 28:2777–2796.

https://doi.org/10.1007/s00180-013-0428-3
- Google Scholar
1. Li L
2. Li Y
3. Huang Y
4. Ouyang Y
5. Zhu Y
6. Wang Y
7. Guo X
8. Yuan Y
9. Gong K
(2018) Long non‐coding RNA MIF‐AS1 promotes gastric cancer cell proliferation and reduces apoptosis to upregulate NDUFA4
Cancer Science 109:3714–3725.

https://doi.org/10.1111/cas.13801
- PubMed
- Google Scholar
1. Mafessoni F
2. Grote S
3. de Filippo C
4. Slon V
5. Kolobova KA
6. Viola B
7. Markin SV
8. Chintalapati M
9. Peyrégne S
10. Skov L
11. Skoglund P
12. Krivoshapkin AI
13. Derevianko AP
14. Meyer M
15. Kelso J
16. Peter B
17. Prüfer K
18. Pääbo S
(2020) A high-coverage Neandertal genome from Chagyrskaya Cave
PNAS 117:15132–15136.

https://doi.org/10.1073/pnas.2004944117
- PubMed
- Google Scholar
1. Mallick S
2. Li H
3. Lipson M
4. Mathieson I
5. Gymrek M
6. Racimo F
7. Zhao M
8. Chennagiri N
9. Nordenfelt S
10. Tandon A
11. Skoglund P
12. Lazaridis I
13. Sankararaman S
14. Fu Q
15. Rohland N
16. Renaud G
17. Erlich Y
18. Willems T
19. Gallo C
20. Spence JP
21. Song YS
22. Poletti G
23. Balloux F
24. van Driem G
25. de Knijff P
26. Romero IG
27. Jha AR
28. Behar DM
29. Bravi CM
30. Capelli C
31. Hervig T
32. Moreno-Estrada A
33. Posukh OL
34. Balanovska E
35. Balanovsky O
36. Karachanak-Yankova S
37. Sahakyan H
38. Toncheva D
39. Yepiskoposyan L
40. Tyler-Smith C
41. Xue Y
42. Abdullah MS
43. Ruiz-Linares A
44. Beall CM
45. Di Rienzo A
46. Jeong C
47. Starikovskaya EB
48. Metspalu E
49. Parik J
50. Villems R
51. Henn BM
52. Hodoglugil U
53. Mahley R
54. Sajantila A
55. Stamatoyannopoulos G
56. Jts W
57. Khusainova R
58. Khusnutdinova E
59. Litvinov S
60. Ayodo G
61. Comas D
62. Hammer MF
63. Kivisild T
64. Klitz W
65. Winkler CA
66. Labuda D
67. Bamshad M
68. Jorde LB
69. Tishkoff SA
70. Watkins WS
71. Metspalu M
72. Dryomov S
73. Sukernik R
74. Singh L
75. Thangaraj K
76. Pääbo S
77. Kelso J
78. Patterson N
79. Reich D
(2016) The Simons Genome Diversity Project: 300 genomes from 142 diverse populations
Nature 538:201–206.

https://doi.org/10.1038/nature18964
- PubMed
- Google Scholar
1. Marcus JH
2. Novembre J
(2017) Visualizing the geography of genetic variants
Bioinformatics 33:594–595.

https://doi.org/10.1093/bioinformatics/btw643
- PubMed
- Google Scholar
(2012a) A haplotype at STAT2 Introgressed from neanderthals and serves as a candidate of positive selection in Papua New Guinea
American Journal of Human Genetics 91:265–274.

https://doi.org/10.1016/j.ajhg.2012.06.015
- Google Scholar
(2012b) Global genetic variation at OAS1 provides evidence of archaic admixture in Melanesian populations
Molecular Biology and Evolution 29:1513–1520.

https://doi.org/10.1093/molbev/msr301
- Google Scholar
1. Meyer M
2. Kircher M
3. Gansauge MT
4. Li H
5. Racimo F
6. Mallick S
7. Schraiber JG
8. Jay F
9. Prüfer K
10. de Filippo C
11. Sudmant PH
12. Alkan C
13. Fu Q
14. Do R
15. Rohland N
16. Tandon A
17. Siebauer M
18. Green RE
19. Bryc K
20. Briggs AW
21. Stenzel U
22. Dabney J
23. Shendure J
24. Kitzman J
25. Hammer MF
26. Shunkov MV
27. Derevianko AP
28. Patterson N
29. Andrés AM
30. Eichler EE
31. Slatkin M
32. Reich D
33. Kelso J
34. Pääbo S
(2012) A high-coverage genome sequence from an archaic Denisovan individual
Science 338:222–226.

https://doi.org/10.1126/science.1224344
- PubMed
- Google Scholar
1. Mikocziova I
2. Gidoni M
3. Lindeman I
4. Peres A
5. Snir O
6. Yaari G
7. Sollid LM
(2020) Polymorphisms in human immunoglobulin heavy chain variable genes and their upstream regions
Nucleic Acids Research 48:5499–5510.

https://doi.org/10.1093/nar/gkaa310
- PubMed
- Google Scholar
1. Moreno-Mayar JV
2. Potter BA
3. Vinner L
4. Steinrücken M
5. Rasmussen S
6. Terhorst J
7. Kamm JA
8. Albrechtsen A
9. Malaspinas AS
10. Sikora M
11. Reuther JD
12. Irish JD
13. Malhi RS
14. Orlando L
15. Song YS
16. Nielsen R
17. Meltzer DJ
18. Willerslev E
(2018) Terminal Pleistocene Alaskan genome reveals first founding population of Native Americans
Nature 553:203–207.

https://doi.org/10.1038/nature25173
- PubMed
- Google Scholar
1. Nattestad M
2. Goodwin S
3. Ng K
4. Baslan T
5. Sedlazeck FJ
6. Rescheneder P
7. Garvin T
8. Fang H
9. Gurtowski J
10. Hutton E
11. Tseng E
12. Chin CS
13. Beck T
14. Sundaravadanam Y
15. Kramer M
16. Antoniou E
17. McPherson JD
18. Hicks J
19. McCombie WR
20. Schatz MC
(2018) Complex rearrangements and oncogene amplifications revealed by long-read DNA and RNA sequencing of a breast cancer cell line
Genome Research 28:1126–1135.

https://doi.org/10.1101/gr.231100.117
- PubMed
- Google Scholar
1. Ongen H
2. Buil A
3. Brown AA
4. Dermitzakis ET
5. Delaneau O
(2016) Fast and efficient QTL mapper for thousands of molecular phenotypes
Bioinformatics 32:1479–1485.

https://doi.org/10.1093/bioinformatics/btv722
- PubMed
- Google Scholar
1. Pandey R
2. Bakay M
3. Hain HS
4. Strenkowski B
5. Yermakova A
6. Kushner JA
7. Orange JS
8. Hakonarson H
(2019) The Autoimmune Disorder Susceptibility Gene Restrains NK Cell Function in YTS NK Cell Line and Knockout Mice
Frontiers in Immunology 10:68.

https://doi.org/10.3389/fimmu.2019.00068
- PubMed
- Google Scholar
1. Pedersen BS
2. Quinlan AR
(2017) cyvcf2: fast, flexible variant analysis with Python
Bioinformatics 33:1867–1869.

https://doi.org/10.1093/bioinformatics/btx057
- PubMed
- Google Scholar
1. Perry GH
2. Dominy NJ
3. Claw KG
4. Lee AS
5. Fiegler H
6. Redon R
7. Werner J
8. Villanea FA
9. Mountain JL
10. Misra R
11. Carter NP
12. Lee C
13. Stone AC
(2007) Diet and the evolution of human amylase gene copy number variation
Nature Genetics 39:1256–1260.

https://doi.org/10.1038/ng2123
- PubMed
- Google Scholar
(1999) Population growth of human Y chromosomes: a study of Y chromosome microsatellites
Molecular Biology and Evolution 16:1791–1798.

https://doi.org/10.1093/oxfordjournals.molbev.a026091
- PubMed
- Google Scholar
1. Prüfer K
2. Racimo F
3. Patterson N
4. Jay F
5. Sankararaman S
6. Sawyer S
7. Heinze A
8. Renaud G
9. Sudmant PH
10. de Filippo C
11. Li H
12. Mallick S
13. Dannemann M
14. Fu Q
15. Kircher M
16. Kuhlwilm M
17. Lachmann M
18. Meyer M
19. Ongyerth M
20. Siebauer M
21. Theunert C
22. Tandon A
23. Moorjani P
24. Pickrell J
25. Mullikin JC
26. Vohr SH
27. Green RE
28. Hellmann I
29. Johnson PLF
30. Blanche H
31. Cann H
32. Kitzman JO
33. Shendure J
34. Eichler EE
35. Lein ES
36. Bakken TE
37. Golovanova LV
38. Doronichev VB
39. Shunkov MV
40. Derevianko AP
41. Viola B
42. Slatkin M
43. Reich D
44. Kelso J
45. Pääbo S
(2014) The complete genome sequence of a Neanderthal from the Altai Mountains
Nature 505:43–49.

https://doi.org/10.1038/nature12886
- PubMed
- Google Scholar
1. Prüfer K
2. de Filippo C
3. Grote S
4. Mafessoni F
5. Korlević P
6. Hajdinjak M
7. Vernot B
8. Skov L
9. Hsieh P
10. Peyrégne S
11. Reher D
12. Hopfe C
13. Nagel S
14. Maricic T
15. Fu Q
16. Theunert C
17. Rogers R
18. Skoglund P
19. Chintalapati M
20. Dannemann M
21. Nelson BJ
22. Key FM
23. Rudan P
24. Kućan Ž
25. Gušić I
26. Golovanova LV
27. Doronichev VB
28. Patterson N
29. Reich D
30. Eichler EE
31. Slatkin M
32. Schierup MH
33. Andrés AM
34. Kelso J
35. Meyer M
36. Pääbo S
(2017) A high-coverage Neandertal genome from Vindija cave in Croatia
Science 358:655–658.

https://doi.org/10.1126/science.aao1887
- PubMed
- Google Scholar
1. Purcell S
2. Neale B
3. Todd-Brown K
4. Thomas L
5. Ferreira MAR
6. Bender D
7. Maller J
8. Sklar P
9. de Bakker PIW
10. Daly MJ
11. Sham PC
(2007) PLINK: a tool set for whole-genome association and population-based linkage analyses
American Journal of Human Genetics 81:559–575.

https://doi.org/10.1086/519795
- PubMed
- Google Scholar
1. Racimo F
2. Gokhman D
3. Fumagalli M
4. Ko A
5. Hansen T
6. Moltke I
7. Albrechtsen A
8. Carmel L
9. Huerta-Sánchez E
10. Nielsen R
(2017) Archaic adaptive introgression in TBX15/WARS2
Molecular Biology and Evolution 34:509–524.

https://doi.org/10.1093/molbev/msw283
- PubMed
- Google Scholar
1. Rasmussen S
2. Allentoft ME
3. Nielsen K
4. Orlando L
5. Sikora M
6. Sjögren KG
7. Pedersen AG
8. Schubert M
9. Van Dam A
10. Kapel CMO
11. Nielsen HB
12. Brunak S
13. Avetisyan P
14. Epimakhov A
15. Khalyapin MV
16. Gnuni A
17. Kriiska A
18. Lasak I
19. Metspalu M
20. Moiseyev V
21. Gromov A
22. Pokutta D
23. Saag L
24. Varul L
25. Yepiskoposyan L
26. Sicheritz-Pontén T
27. Foley RA
28. Lahr MM
29. Nielsen R
30. Kristiansen K
31. Willerslev E
(2015) Early divergent strains of Yersinia pestis in Eurasia 5,000 years ago
Cell 163:571–582.

https://doi.org/10.1016/j.cell.2015.10.009
- PubMed
- Google Scholar
1. Roadmap Epigenomics Consortium
2. Kundaje A
3. Meuleman W
4. Ernst J
5. Bilenky M
6. Yen A
7. Heravi-Moussavi A
8. Kheradpour P
9. Zhang Z
10. Wang J
11. Ziller MJ
12. Amin V
13. Whitaker JW
14. Schultz MD
15. Ward LD
16. Sarkar A
17. Quon G
18. Sandstrom RS
19. Eaton ML
20. Wu Y-C
21. Pfenning AR
22. Wang X
23. Claussnitzer M
24. Liu Y
25. Coarfa C
26. Harris RA
27. Shoresh N
28. Epstein CB
29. Gjoneska E
30. Leung D
31. Xie W
32. Hawkins RD
33. Lister R
34. Hong C
35. Gascard P
36. Mungall AJ
37. Moore R
38. Chuah E
39. Tam A
40. Canfield TK
41. Hansen RS
42. Kaul R
43. Sabo PJ
44. Bansal MS
45. Carles A
46. Dixon JR
47. Farh K-H
48. Feizi S
49. Karlic R
50. Kim A-R
51. Kulkarni A
52. Li D
53. Lowdon R
54. Elliott G
55. Mercer TR
56. Neph SJ
57. Onuchic V
58. Polak P
59. Rajagopal N
60. Ray P
61. Sallari RC
62. Siebenthall KT
63. Sinnott-Armstrong NA
64. Stevens M
65. Thurman RE
66. Wu J
67. Zhang B
68. Zhou X
69. Beaudet AE
70. Boyer LA
71. De Jager PL
72. Farnham PJ
73. Fisher SJ
74. Haussler D
75. Jones SJM
76. Li W
77. Marra MA
78. McManus MT
79. Sunyaev S
80. Thomson JA
81. Tlsty TD
82. Tsai L-H
83. Wang W
84. Waterland RA
85. Zhang MQ
86. Chadwick LH
87. Bernstein BE
88. Costello JF
89. Ecker JR
90. Hirst M
91. Meissner A
92. Milosavljevic A
93. Ren B
94. Stamatoyannopoulos JA
95. Wang T
96. Kellis M
(2015) Integrative analysis of 111 reference human epigenomes
Nature 518:317–330.

https://doi.org/10.1038/nature14248
- PubMed
- Google Scholar
(2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data
Bioinformatics 26:139–140.

https://doi.org/10.1093/bioinformatics/btp616
- PubMed
- Google Scholar
1. Rodriguez OL
2. Gibson WS
3. Parks T
4. Emery M
5. Powell J
6. Strahl M
7. Deikus G
8. Auckland K
9. Eichler EE
10. Marasco WA
11. Sebra R
12. Sharp AJ
13. Smith ML
14. Bashir A
15. Watson CT
(2020) A Novel Framework for Characterizing Genomic Haplotype Diversity in the Human Immunoglobulin Heavy Chain Locus
Frontiers in Immunology 11:2136.

https://doi.org/10.3389/fimmu.2020.02136
- PubMed
- Google Scholar
1. Rotival M
2. Quintana-Murci L
(2020) Functional consequences of archaic introgression and their impact on fitness
Genome Biology 21:3.

https://doi.org/10.1186/s13059-019-1920-z
- PubMed
- Google Scholar
1. Sabeti PC
2. Walsh E
3. Schaffner SF
4. Varilly P
5. Fry B
6. Hutcheson HB
7. Cullen M
8. Mikkelsen TS
9. Roy J
10. Patterson N
11. Cooper R
12. Reich D
13. Altshuler D
14. O’Brien S
15. Lander ES
(2005) The case for selection at CCR5-Delta32
PLOS Biology 3:e378.

https://doi.org/10.1371/journal.pbio.0030378
- PubMed
- Google Scholar
1. Saitou M
2. Gokcumen O
(2019) Resolving the Insertion Sites of Polymorphic Duplications Reveals a HERC2 Haplotype under Selection
Genome Biology and Evolution 11:1679–1690.

https://doi.org/10.1093/gbe/evz107
- PubMed
- Google Scholar
1. Sams AJ
2. Dumaine A
3. Nédélec Y
4. Yotova V
5. Alfieri C
6. Tanner JE
7. Messer PW
8. Barreiro LB
(2021) Adaptively introgressed Neandertal haplotype at the OAS locus functionally impacts innate immune responses in humans
Genome Biology 17:51466.

https://doi.org/10.1101/051466
- Google Scholar
1. Sankararaman S
2. Patterson N
3. Li H
4. Pääbo S
5. Reich D
(2012) The date of interbreeding between Neandertals and modern humans
PLOS Genetics 8:e1002947.

https://doi.org/10.1371/journal.pgen.1002947
- PubMed
- Google Scholar
(2016) The Combined Landscape of Denisovan and Neanderthal Ancestry in Present-Day Humans
Current Biology 26:1241–1247.

https://doi.org/10.1016/j.cub.2016.03.037
- PubMed
- Google Scholar
(2018) From genome-wide associations to candidate causal variants by statistical fine-mapping
Nature Reviews Genetics 19:491–504.

https://doi.org/10.1038/s41576-018-0016-z
- PubMed
- Google Scholar
1. Sedlazeck FJ
2. Lee H
3. Darby CA
4. Schatz MC
(2018a) Piercing the dark matter: bioinformatics of long-range sequencing and mapping
Nature Reviews Genetics 19:329–346.

https://doi.org/10.1038/s41576-018-0003-4
- Google Scholar
(2018b) Accurate detection of complex structural variations using single-molecule sequencing
Nature Methods 15:461–468.

https://doi.org/10.1038/s41592-018-0001-7
- Google Scholar
1. Shriver MD
2. Kennedy GC
3. Parra EJ
4. Lawson HA
5. Sonpar V
6. Huang J
7. Akey JM
8. Jones KW
(2004) The genomic distribution of population substructure in four populations using 8,525 autosomal SNPs
Human Genomics 1:274–286.

https://doi.org/10.1186/1479-7364-1-4-274
- PubMed
- Google Scholar
(2018) Accurate genotyping across variant classes and lengths using variant graphs
Nature Genetics 50:1054–1059.

https://doi.org/10.1038/s41588-018-0145-5
- PubMed
- Google Scholar
Preprint
1. Sirén J
2. Monlong J
3. Chang X
4. Novak AM
5. Eizenga JM
6. Markello C
7. Sibbesen JA
8. Hickey G
9. Chang PC
10. Carroll A
11. Gupta N
12. Gabriel S
13. Blackwell TW
14. Ratan A
15. Taylor KD
16. Rich SS
17. Rotter JI
18. Haussler D
19. Garrison E
20. Paten B
(2021) Genotyping Common, Large Structural Variations in 5,202 Genomes Using Pangenomes, the Giraffe Mapper, and the vg Toolkit
bioRxiv.

https://doi.org/10.1101/2020.12.04.412486
- Google Scholar
1. Souilmi Y
2. Lauterbur ME
3. Tobler R
4. Huber CD
5. Johar AS
6. Moradi SV
7. Johnston WA
8. Krogan NJ
9. Alexandrov K
10. Enard D
(2021) An ancient viral epidemic involving host coronavirus interacting genes more than 20,000 years ago in East Asia
Current Biology 31:3704.

https://doi.org/10.1016/j.cub.2021.07.052
- PubMed
- Google Scholar
1. Stephens JC
2. Reich DE
3. Goldstein DB
4. Shin HD
5. Smith MW
6. Carrington M
7. Winkler C
8. Huttley GA
9. Allikmets R
10. Schriml L
11. Gerrard B
12. Malasky M
13. Ramos MD
14. Morlot S
15. Tzetis M
16. Oddoux C
17. di Giovine FS
18. Nasioulas G
19. Chandler D
20. Aseev M
21. Hanson M
22. Kalaydjieva L
23. Glavac D
24. Gasparini P
25. Kanavakis E
26. Claustres M
27. Kambouris M
28. Ostrer H
29. Duff G
30. Baranov V
31. Sibul H
32. Metspalu A
33. Goldman D
34. Martin N
35. Duffy D
36. Schmidtke J
37. Estivill X
38. O’Brien SJ
39. Dean M
(1998) Dating the origin of the CCR5-Delta32 AIDS-resistance allele by the coalescence of haplotypes
American Journal of Human Genetics 62:1507–1515.

https://doi.org/10.1086/301867
- PubMed
- Google Scholar
1. Storey JD
2. Tibshirani R
(2003) Statistical significance for genomewide studies
PNAS 100:9440–9445.

https://doi.org/10.1073/pnas.1530509100
- PubMed
- Google Scholar
Software
1. Storey JD
2. Bass AJ
3. Dabney A
4. Robinson D
5. Warnes G
(2020) qvalue: Q-value estimation for false discovery rate control, version 2.21.0
GitHub.

http://github.com/jdstorey/qvalue
1. Sudmant PH
2. Mallick S
3. Nelson BJ
4. Hormozdiari F
5. Krumm N
6. Huddleston J
7. Coe BP
8. Baker C
9. Nordenfelt S
10. Bamshad M
11. Jorde LB
12. Posukh OL
13. Sahakyan H
14. Watkins WS
15. Yepiskoposyan L
16. Abdullah MS
17. Bravi CM
18. Capelli C
19. Hervig T
20. Jts W
21. Tyler-Smith C
22. van Driem G
23. Romero IG
24. Jha AR
25. Karachanak-Yankova S
26. Toncheva D
27. Comas D
28. Henn B
29. Kivisild T
30. Ruiz-Linares A
31. Sajantila A
32. Metspalu E
33. Parik J
34. Villems R
35. Starikovskaya EB
36. Ayodo G
37. Beall CM
38. Di Rienzo A
39. Hammer MF
40. Khusainova R
41. Khusnutdinova E
42. Klitz W
43. Winkler C
44. Labuda D
45. Metspalu M
46. Tishkoff SA
47. Dryomov S
48. Sukernik R
49. Patterson N
50. Reich D
51. Eichler EE
(2015) Global diversity, population stratification, and selection of human copy-number variation
Science 349:aab3761.

https://doi.org/10.1126/science.aab3761
- PubMed
- Google Scholar
1. Suzuki T
2. Maranda B
3. Sakagami T
4. Catellier P
5. Couture CY
6. Carey BC
7. Chalk C
8. Trapnell BC
(2011) Hereditary pulmonary alveolar proteinosis caused by recessive CSF2RB mutations
The European Respiratory Journal 37:201–204.

https://doi.org/10.1183/09031936.00090610
- PubMed
- Google Scholar
1. Tewhey R
2. Kotliar D
3. Park DS
4. Liu B
5. Winnicki S
6. Reilly SK
7. Andersen KG
8. Mikkelsen TS
9. Lander ES
10. Schaffner SF
11. Sabeti PC
(2018) Direct Identification of Hundreds of Expression-Modulating Variants using a Multiplexed Reporter Assay
Cell 172:1132–1134.

https://doi.org/10.1016/j.cell.2018.02.021
- PubMed
- Google Scholar
1. The 1000 Genomes Project Consortium
(2015) A global reference for human genetic variation
Nature 526:68–74.

https://doi.org/10.1038/nature15393
- PubMed
- Google Scholar
1. The GTEx Consortium
(2020) The GTEx Consortium atlas of genetic regulatory effects across human tissues
Science 369:1318–1330.

https://doi.org/10.1126/science.aaz1776
- PubMed
- Google Scholar
1. The International HapMap Consortium
(2007) A second generation human haplotype map of over 3.1 million SNPs
Nature 449:851–861.

https://doi.org/10.1038/nature06258
- PubMed
- Google Scholar
1. Tremblay M
2. Vézina H
(2000) New estimates of intergenerational time intervals for the calculation of age and origins of mutations
American Journal of Human Genetics 66:651–658.

https://doi.org/10.1086/302770
- PubMed
- Google Scholar
(2019) High-throughput identification of human SNPs affecting regulatory element activity
Nature Genetics 51:1160–1169.

https://doi.org/10.1038/s41588-019-0455-2
- PubMed
- Google Scholar
1. VanLiere JM
2. Rosenberg NA
(2008) Mathematical properties of the r2 measure of linkage disequilibrium
Theoretical Population Biology 74:130–137.

https://doi.org/10.1016/j.tpb.2008.05.006
- PubMed
- Google Scholar
1. Vernot B
2. Tucci S
3. Kelso J
4. Schraiber JG
5. Wolf AB
6. Gittelman RM
7. Dannemann M
8. Grote S
9. McCoy RC
10. Norton H
11. Scheinfeldt LB
12. Merriwether DA
13. Koki G
14. Friedlaender JS
15. Wakefield J
16. Pääbo S
17. Akey JM
(2016) Excavating Neandertal and Denisovan DNA from the genomes of Melanesian individuals
Science 352:235–239.

https://doi.org/10.1126/science.aad9416
- PubMed
- Google Scholar
(2013) Detecting natural selection in genomic data
Annual Review of Genetics 47:97–120.

https://doi.org/10.1146/annurev-genet-111212-133526
- PubMed
- Google Scholar
1. Wang Y
2. Lu D
3. Chung YJ
4. Xu S
(2018) Genetic structure, divergence and admixture of Han Chinese, Japanese and Korean populations
Hereditas 155:19.

https://doi.org/10.1186/s41065-018-0057-5
- PubMed
- Google Scholar
1. Watson CT
2. Steinberg KM
3. Huddleston J
4. Warren RL
5. Malig M
6. Schein J
7. Willsey AJ
8. Joy JB
9. Scott JK
10. Graves TA
11. Wilson RK
12. Holt RA
13. Eichler EE
14. Breden F
(2013) Complete haplotype sequence of the human immunoglobulin heavy-chain variable, diversity, and joining genes and characterization of allelic and copy-number variation
American Journal of Human Genetics 92:530–546.

https://doi.org/10.1016/j.ajhg.2013.03.004
- PubMed
- Google Scholar
(2017) The Individual and Population Genetics of Antibody Immunity
Trends in Immunology 38:459–470.

https://doi.org/10.1016/j.it.2017.04.003
- PubMed
- Google Scholar
1. Wojcik GL
2. Graff M
3. Nishimura KK
4. Tao R
5. Haessler J
6. Gignoux CR
7. Highland HM
8. Patel YM
9. Sorokin EP
10. Avery CL
11. Belbin GM
12. Bien SA
13. Cheng I
14. Cullina S
15. Hodonsky CJ
16. Hu Y
17. Huckins LM
18. Jeff J
19. Justice AE
20. Kocarnik JM
21. Lim U
22. Lin BM
23. Lu Y
24. Nelson SC
25. Park SSL
26. Poisner H
27. Preuss MH
28. Richard MA
29. Schurmann C
30. Setiawan VW
31. Sockell A
32. Vahi K
33. Verbanck M
34. Vishnu A
35. Walker RW
36. Young KL
37. Zubair N
38. Acuña-Alonso V
39. Ambite JL
40. Barnes KC
41. Boerwinkle E
42. Bottinger EP
43. Bustamante CD
44. Caberto C
45. Canizales-Quinteros S
46. Conomos MP
47. Deelman E
48. Do R
49. Doheny K
50. Fernández-Rhodes L
51. Fornage M
52. Hailu B
53. Heiss G
54. Henn BM
55. Hindorff LA
56. Jackson RD
57. Laurie CA
58. Laurie CC
59. Li Y
60. Lin DY
61. Moreno-Estrada A
62. Nadkarni G
63. Norman PJ
64. Pooler LC
65. Reiner AP
66. Romm J
67. Sabatti C
68. Sandoval K
69. Sheng X
70. Stahl EA
71. Stram DO
72. Thornton TA
73. Wassel CL
74. Wilkens LR
75. Winkler CA
76. Yoneyama S
77. Buyske S
78. Haiman CA
79. Kooperberg C
80. Le Marchand L
81. Loos RJF
82. Matise TC
83. North KE
84. Peters U
85. Kenny EE
86. Carlson CS
(2019) Genetic analyses of diverse populations improves discovery for complex traits
Nature 570:514–518.

https://doi.org/10.1038/s41586-019-1310-4
- PubMed
- Google Scholar
1. Wright S
(1949) The genetical structure of populations
Annals of Eugenics 15:323–354.

https://doi.org/10.1111/j.1469-1809.1949.tb02451.x
- Google Scholar
1. Yair S
2. Lee KM
3. Coop G
(2021) The timing of human adaptation from Neanderthal introgression
Genetics 218:iyab052.

https://doi.org/10.1093/genetics/iyab052
- PubMed
- Google Scholar
1. Yi X
2. Liang Y
3. Huerta-Sanchez E
4. Jin X
5. Zxp C
6. Pool JE
7. Xu X
8. Jiang H
9. Vinckenbosch N
10. Korneliussen TS
11. Hancheng Z
12. Liu T
13. He W
14. Li K
15. Luo R
16. Nie X
17. Wu H
18. Zhao M
19. Cao H
20. Zou J
21. Shan Y
22. Shuzheng L
23. Yang Q
24. Asan N
25. Tian P
26. Xu G
27. Liu J
28. Jiang X
29. Wu T
30. Zhou R
31. Tang G
32. Qin M
33. Wang J
34. Feng T
35. Li S
36. Huasang G
37. Luosang J
38. Wang W
39. Chen F
40. Wang Y
41. Zheng X
42. Li Z
43. Bianba Z
44. Yang G
45. Wang X
46. Tang S
47. Gao G
48. Chen Y
49. Luo Z
50. Gusang L
51. Cao Z
52. Zhang Q
53. Ouyang W
54. Ren X
55. Liang H
56. Huisong Z
57. Huang Y
58. Li J
59. Bolund L
60. Kristiansen K
61. Li Y
62. Zhang Y
63. Zhang X
64. Li R
65. Songgang L
66. Yang H
67. Nielsen R
68. Jun W
69. Jian W
(2010) Sequencing of 50 human exomes reveals adaptation to high altitude
Science 329:75–78.

https://doi.org/10.1126/science.1190371
- PubMed
- Google Scholar
1. Zeberg H
2. Pääbo S
(2020) The major genetic risk factor for severe COVID-19 is inherited from Neanderthals
Nature 587:610–612.

https://doi.org/10.1038/s41586-020-2818-3
- PubMed
- Google Scholar
1. Zeberg H
2. Pääbo S
(2021) A genomic region associated with protection against severe COVID-19 is inherited from Neandertals
PNAS 118:e2026309118.

https://doi.org/10.1073/pnas.2026309118
- PubMed
- Google Scholar
1. Zhang Z
2. Bai M
3. Barbosa GO
4. Chen A
5. Wei Y
6. Luo S
7. Wang X
8. Wang B
9. Tsukui T
10. Li H
11. Sheppard D
12. Kornberg TB
13. Dk M
(2020) Broadly conserved roles of TMEM131 family proteins in intracellular collagen assembly and secretory cargo trafficking
Science Advances 6:eaay7667.

https://doi.org/10.1126/sciadv.aay7667
- Google Scholar
1. Zhang JY
2. Roberts H
3. Flores DSC
4. Cutler AJ
5. Brown AC
6. Whalley JP
7. Mielczarek O
8. Buck D
9. Lockstone H
10. Xella B
11. Oliver K
12. Corton C
13. Betteridge E
14. Bashford-Rogers R
15. Knight JC
16. Todd JA
17. Band G
(2021) Using de novo assembly to identify structural variation of eight complex immune system gene regions
PLOS Computational Biology 17:e1009254.

https://doi.org/10.1371/journal.pcbi.1009254
- PubMed
- Google Scholar
1. Zhao X
2. Collins RL
3. Lee WP
4. Weber AM
5. Jun Y
6. Zhu Q
7. Weisburd B
8. Huang Y
9. Audano PA
10. Wang H
11. Walker M
12. Lowther C
13. Fu J
14. Gerstein MB
15. Devine SE
16. Marschall T
17. Korbel JO
18. Eichler EE
19. Chaisson MJP
20. Lee C
21. Mills RE
22. Brand H
23. Talkowski ME
(2021) Expectations and blind spots for structural variation detection from long-read assemblies and short-read genome sequencing technologies
American Journal of Human Genetics 108:919–928.

https://doi.org/10.1016/j.ajhg.2021.03.014
- PubMed
- Google Scholar
(2021) A Neanderthal OAS1 isoform protects individuals of European ancestry against COVID-19 susceptibility and severity
Nature Medicine 27:659–667.

https://doi.org/10.1038/s41591-021-01281-1
- PubMed
- Google Scholar
(2020) A robust benchmark for detection of germline large deletions and insertions
Nature Biotechnology 38:1347–1355.

https://doi.org/10.1038/s41587-020-0538-8
- PubMed
- Google Scholar

Article and author information

Author details

Stephanie M Yan

Department of Biology, Johns Hopkins University, Baltimore, Baltimore, United States

Contribution
Conceptualization, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6880-465X
Rachel M Sherman

Department of Computer Science, Johns Hopkins University, Baltimore, United States

Contribution
Data curation, Investigation, Software, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1750-8428
Dylan J Taylor

Department of Biology, Johns Hopkins University, Baltimore, Baltimore, United States

Contribution
Formal analysis, Investigation, Methodology, Software, Visualization, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-5806-4494
Divya R Nair

Department of Biology, Johns Hopkins University, Baltimore, Baltimore, United States

Contribution
Investigation, Software

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3374-8424
Andrew N Bortvin

Department of Biology, Johns Hopkins University, Baltimore, Baltimore, United States

Contribution
Investigation, Software, Visualization

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8784-6786
Michael C Schatz
1. Department of Biology, Johns Hopkins University, Baltimore, Baltimore, United States
2. Department of Computer Science, Johns Hopkins University, Baltimore, United States
Contribution
Conceptualization, Funding acquisition, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4118-4446
Rajiv C McCoy

Department of Biology, Johns Hopkins University, Baltimore, Baltimore, United States

Contribution
Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing

For correspondence
rajiv.mccoy@jhu.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0615-146X

Funding

National Institutes of Health (R35GM133747)

Rajiv C McCoy

National Science Foundation (DBI-1350041)

Michael C Schatz

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Tim O’Connor, Sai Chen, John Kim, and members of the McCoy lab for feedback and helpful discussions. We also thank the staff at the Maryland Advanced Research Computing Center for computing support. This work is supported by the National Institutes of Health (NIH) grant R35GM133747 to RCM and the US National Science Foundation grant DBI-1350041 to MCS.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.