William Klitz

Reconstructing Native American population history.

The peopling of the Americas has been the subject of extensive genetic, archaeological and linguistic research; however, central questions remain unresolved. One contentious issue is whether the settlement occurred by means of a single migration or multiple streams of migration from Siberia. The pattern of dispersals within the Americas is also poorly understood. To address these questions at a higher resolution than was previously possible, we assembled data from 52 Native American and 17 Siberian groups genotyped at 364,470 single nucleotide polymorphisms. Here we show that Native Americans descend from at least three streams of Asian gene flow. Most descend entirely from a single ancestral population that we call ‘First American’. However, speakers of Eskimo–Aleut languages from the Arctic inherit almost half their ancestry from a second stream of Asian gene flow, and the Na-Dene-speaking Chipewyan from Canada inherit roughly one-tenth of their ancestry from a third stream. We show that the initial peopling followed a southward expansion facilitated by the coast, with sequential population splits and little gene flow after divergence, especially in South America. A major exception is in Chibchan speakers on both sides of the Panama isthmus, who have ancestry from both North and South America.

Geographic patterns of genome admixture in Latin American Mestizos.

The large and diverse population of Latin America is potentially a powerful resource for elucidating the genetic basis of complex traits through admixture mapping. However, no genome-wide characterization of admixture across Latin America has yet been attempted. Here, we report an analysis of admixture in thirteen Mestizo populations (i.e. in regions of mainly European and Native settlement) from seven countries in Latin America based on data for 678 autosomal and 29 X-chromosome microsatellites. We found extensive variation in Native American and European ancestry (and generally low levels of African ancestry) among populations and individuals, and evidence that admixture across Latin America has often involved predominantly European men and both Native and African women. An admixture analysis allowing for Native American population subdivision revealed a differentiation of the Native American ancestry amongst Mestizos. This observation is consistent with the genetic structure of pre-Columbian populations and with admixture having involved Natives from the area where the Mestizo examined are located. Our findings agree with available information on the demographic history of Latin America and have a number of implications for the design of association studies in population from the region.

Ancestral Asian Source(s) of New World Y-Chromosome Founder Haplotypes

Haplotypes constructed from Y-chromosome markers were used to trace the origins of Native Americans. Our sample consisted of 2,198 males from 60 global populations, including 19 Native American and 15 indigenous North Asian groups. A set of 12 biallelic polymorphisms gave rise to 14 unique Y-chromosome haplotypes that were unevenly distributed among the populations. Combining multiallelic variation at two Y-linked microsatellites (DYS19 and DXYS156Y) with the unique haplotypes results in a total of 95 combination haplotypes. Contra previous findings based on Y- chromosome data, our new results suggest the possibility of more than one Native American paternal founder haplotype. We postulate that, of the nine unique haplotypes found in Native Americans, haplotypes 1C and 1F are the best candidates for major New World founder haplotypes, whereas haplotypes 1B, 1I, and 1U may either be founder haplotypes and/or have arrived in the New World via recent admixture. Two of the other four haplotypes (YAP+ haplotypes 4 and 5) are probably present because of post-Columbian admixture, whereas haplotype 1G may have originated in the New World, and the Old World source of the final New World haplotype (1D) remains unresolved. The contrasting distribution patterns of the two major candidate founder haplotypes in Asia and the New World, as well as the results of a nested cladistic analysis, suggest the possibility of more than one paternal migration from the general region of Lake Baikal to the Americas.

A Genomewide Admixture Map for Latino Populations

Admixture mapping is an economical and powerful approach for localizing disease genes in populations of recently mixed ancestry and has proven successful in African Americans. The method holds equal promise for Latinos, who typically inherit a mix of European, Native American, and African ancestry. However, admixture mapping in Latinos has not been practical because of the lack of a map of ancestry-informative markers validated in Native American and other populations. To address this, we screened multiple databases, containing millions of markers, to identify 4,186 markers that were putatively informative for determining the ancestry of chromosomal segments in Latino populations. We experimentally validated each of these markers in at least 232 new Latino, European, Native American, and African samples, and we selected a subset of 1,649 markers to form an admixture map. An advantage of our strategy is that we focused our map on markers distinguishing Native American from other ancestries and restricted it to markers with very similar frequencies in Europeans and Africans, which decreased the number of markers needed and minimized the possibility of false disease associations. We evaluated the effectiveness of our map for localizing disease genes in four Latino populations from both North and South America.

Association mapping of disease loci, by use of a pooled DNA genomic screen.

Genomic screening to map disease loci by association requires automation, pooling of DNA samples, and 3,000-6,000 highly polymorphic, evenly spaced microsatellite markers. Case-control samples can be used in an initial screen, followed by family-based data to confirm marker associations. Association mapping is relevant to genetic studies of complex diseases in which linkage analysis may be less effective and to cases in which multigenerational data are difficult to obtain, including rare or late-onset conditions and infectious diseases. The method can also be used effectively to follow up and confirm regions identified in linkage studies or to investigate candidate disease loci. Study designs can incorporate disease heterogeneity and interaction effects by appropriate subdivision of samples before screening. Here we report use of pooled DNA amplifications-the accurate determination of marker-disease associations for both case-control and nuclear family-based data-including application of correction methods for stutter artifact and preferential amplification. These issues, combined with a discussion of both statistical power and experimental design to define the necessary requirements for detecting of disease loci while virtually eliminating false positives, suggest the feasibility and efficiency of association mapping using pooled DNA screening.

High-Resolution Patterns of Meiotic Recombination across the Human Major Histocompatibility Complex

Definitive characteristics of meiotic recombination events over large (i.e., >1 Mb) segments of the human genome remain obscure, yet they are essential for establishing the haplotypic structure of the genome and for efficient mapping of complex traits. We present a high-resolution map of recombination at the kilobase level across a 3.3-Mb interval encompassing the major histocompatibility complex (MHC). Genotyping of 20,031 single sperm from 12 individuals resulted in the identification and fine mapping of 325 recombinant chromosomes within genomic intervals as small as 7 kb. Several principal characteristics of recombination in this region were observed: (1) rates of recombination can differ significantly between individuals; (2) intense hot spots of recombination occur at least every 0.8 Mb but are not necessarily evenly spaced; (3) distribution in the location of recombination events can differ significantly among individuals; (4) between hot spots, low levels of recombination occur fairly evenly across 100-kb segments, suggesting the presence of warm spots of recombination; and (5) specific sequence motifs associate significantly with recombination distribution. These data provide a plausible model for recombination patterns of the human genome overall.

The Simons Genome Diversity Project: 300 genomes from 142 diverse populations

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans.

Antibiotic-refractory Lyme arthritis is associated with HLA-DR molecules that bind a Borrelia burgdorferi peptide

An association has previously been shown between antibiotic-refractory Lyme arthritis, the human histocompatibility leukocyte antigen (HLA)–DR4 molecule, and T cell recognition of an epitope of Borrelia burgdorferi outer-surface protein A (OspA163–175). We studied the frequencies of HLA-DRB1-DQA1-DQB1 haplotypes in 121 patients with antibiotic-refractory or antibiotic-responsive Lyme arthritis and correlated these frequencies with in vitro binding of the OspA163–175 peptide to 14 DRB molecules. Among the 121 patients, the frequencies of HLA-DRB1-DQA1-DQB1 haplotypes were similar to those in control subjects. However, when stratified by antibiotic response, the frequencies of DRB1 alleles in the 71 patients with antibiotic-refractory arthritis differed significantly from those in the 50 antibiotic-responsive patients (log likelihood test, P = 0.006; exact test, P = 0.008; effect size, Wn = 0.38). 7 of the 14 DRB molecules (DRB1*0401, 0101, 0404, 0405, DRB5*0101, DRB1*0402, and 0102) showed strong to weak binding of OspA163–175, whereas the other seven showed negligible or no binding of the peptide. Altogether, 79% of the antibiotic-refractory patients had at least one of the seven known OspA peptide–binding DR molecules compared with 46% of the antibiotic-responsive patients (odds ratio = 4.4; P < 0.001). We conclude that binding of a single spirochetal peptide to certain DRB molecules is a marker for antibiotic-refractory Lyme arthritis and might play a role in the pathogenesis of the disease.

Global diversity, population stratification, and selection of human copy-number variation.

Duplications and deletions in the human genome Duplications and deletions can lead to variation in copy number for genes and genomic loci among humans. Such variants can reveal evolutionary patterns and have implications for human health. Sudmant et al. examined copy-number variation across 236 individual genomes from 125 human populations. Deletions were under more selection, whereas duplications showed more population-specific structure. Interestingly, Oceanic populations retain large duplications postulated to have originated in an ancient Denisovan lineage. Science, this issue 10.1126/science.aab3761 Copy-number variation reveals how selection affects the human genome across the globe. INTRODUCTION Most studies of human genetic variation have focused on single-nucleotide variants (SNVs). However, copy-number variants (CNVs) affect more base pairs of DNA among humans, and yet our understanding of CNV diversity among human populations is limited. RATIONALE We aimed to understand the pattern, selection, and diversity of copy-number variation by analyzing deeply sequenced genomes representing the diversity of all humans. We compared the selective constraints of deletions versus duplications to understand population stratification in the context of the ancestral human genome and to assess differences in CNV load between African and non-African populations. RESULTS We sequenced 236 individual genomes from 125 distinct human populations and identified 14,467 autosomal CNVs and 545 X-linked CNVs with a sequence read-depth approach. Deletions exhibit stronger selective pressure and are better phylogenetic markers of population relationships than duplication polymorphisms. We identified 1036 population-stratified copy-number–variable regions, 295 of which intersect coding regions and 199 of which exhibit extreme signatures of differentiation. Duplicated loci were 1.8-fold more likely to be stratified than deletions but were poorly correlated with flanking genetic diversity. Among these, we highlight a duplication polymorphism restricted to modern Oceanic populations yet also present in the genome of the archaic Denisova hominin. This 225–kilo–base pair (kbp) duplication includes two microRNA genes and is almost fixed among human Papuan-Bougainville genomes. The data allowed us to reconstruct the ancestral human genome and create a more accurate evolutionary framework for the gain and loss of sequences during human evolution. We identified 571 loci that segregate in the human population and another 2026 loci of fixed-copy 2 in all human genomes but absent from the reference genome. The total deletion and duplication load between African and non-African population groups showed no difference after we account for ancestral sequences missing from the human reference. However, we did observe that the relative number of base pairs affected by CNVs compared to single-nucleotide polymorphisms is higher among non-Africans than Africans. CONCLUSION Deletions, duplications, and CNVs have shaped, to different extents, the genetic diversity of human populations by the combined forces of mutation, selection, and demography. Figure Global human CNV diversity and archaic introgression of a chromosome 16 duplication. (Left) The geographic coordinates of populations sampled are indicated on a world map (colored dots). The pie charts show the continental population allele frequency of a single ~225-kbp duplication polymorphism found exclusively among Oceanic populations and an archaic Denisova. (Right) The ancestral structure of this duplication locus (1) and the Denisova duplication structure (2) are shown in relation to their position on chromosome 16. We estimate that the duplication emerged ~440 thousand years ago (ka) in the Denisova and then introgressed into ancestral Papuan populations ~40 ka. In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide–variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Confirmation of an Association Between the TNF(−308) Promoter Polymorphism and Stroke Risk in Children With Sickle Cell Anemia

Background and Purpose— The etiology of stroke in children with sickle cell anemia (SCA) is complex and poorly understood. Growing evidence suggests that genetic factors beyond the sickle cell mutation influence stroke risk in SCA. We previously reported risk associations with polymorphisms in several proinflammatory genes in SCA children with ischemic stroke. The aim of this replication study was to confirm our previous findings of associations between the TNF(−308) G/A, IL4R 503 S/P, and ADRB2 27 Q/E polymorphisms and large vessel stroke risk. Methods— Using previously collected MRA data, we assessed an independent population of SCA children from the multicenter Stroke Prevention Trial in Sickle Cell Anemia (STOP) for the presence or absence of large vessel stenosis. Samples were genotyped for 104 polymorphisms among 65 candidate vascular disease genes. Genotypic associations with risk of large vessel stroke were screened using univariable analysis and compared with results from our original study. Joint analysis of the 2 study populations combined was performed using multivariable logistic regression. Results— A total of 96 children (49 MRA-positive, 47 MRA-negative) were included in this study. Of the SNP associations previously identified in the original study, the TNF(−308) G/A association with large vessel stroke remained significant and the IL4R 503 S/P variant approached significance in the joint analysis of the combined study populations. Consistent with our original findings, the TNF(−308) GG genotype was associated with a >3-fold increased risk of large vessel disease (OR=3.27; 95% CI=1.6, 6.9; P=0.006). Unadjusted analyses also revealed a previously unidentified association between the LTC4S(−444) A/C variant and large vessel stroke risk. Conclusions— Similar findings in 2 independent study populations strongly suggest that the TNF(−308) G/A promoter polymorphism is a clinically important risk factor for large vessel stroke in children with SCA. The previously observed association with the IL4R 503 S/P variant and the novel association with the LTC4S(−444) A/C variant suggest that these loci may also contribute to large vessel stroke risk in children with SCA.