Genetic Consequences of the Transatlantic Slave Trade in the Americas.

According to historical records of transatlantic slavery, traders forcibly deported an estimated 12.5 million people from ports along the Atlantic coastline of Africa between the 16th and 19th centuries, with global impacts reaching to the present day, more than a century and a half after slavery's abolition. Such records have fueled a broad understanding of the forced migration from Africa to the Americas yet remain underexplored in concert with genetic data. Here, we analyzed genotype array data from 50,281 research participants, which-combined with historical shipping documents-illustrate that the current genetic landscape of the Americas is largely concordant with expectations derived from documentation of slave voyages. For instance, genetic connections between people in slave trading regions of Africa and disembarkation regions of the Americas generally mirror the proportion of individuals forcibly moved between those regions. While some discordances can be explained by additional records of deportations within the Americas, other discordances yield insights into variable survival rates and timing of arrival of enslaved people from specific regions of Africa. Furthermore, the greater contribution of African women to the gene pool compared to African men varies across the Americas, consistent with literature documenting regional differences in slavery practices. This investigation of the transatlantic slave trade, which is broad in scope in terms of both datasets and analyses, establishes genetic links between individuals in the Americas and populations across Atlantic Africa, yielding a more comprehensive understanding of the African roots of peoples of the Americas.

Characterization of Prevalence and Health Consequences of Uniparental Disomy in Four Million Individuals from the General Population.

Meiotic nondisjunction and resulting aneuploidy can lead to severe health consequences in humans. Aneuploidy rescue can restore euploidy but may result in uniparental disomy (UPD), the inheritance of both homologs of a chromosome from one parent with no representative copy from the other. Current understanding of UPD is limited to ∼3,300 case subjects for which UPD was associated with clinical presentation due to imprinting disorders or recessive diseases. Thus, the prevalence of UPD and its phenotypic consequences in the general population are unknown. We searched for instances of UPD across 4,400,363 consented research participants from the personal genetics company 23andMe, Inc., and 431,094 UK Biobank participants. Using computationally detected DNA segments identical-by-descent (IBD) and runs of homozygosity (ROH), we identified 675 instances of UPD across both databases. We estimate that UPD is twice as common as previously thought, and we present a machine-learning framework to detect UPD using ROH. While we find a nominally significant association between UPD of chromosome 22 and autism risk, we do not find significant associations between UPD and deleterious traits in the 23andMe database.

Diversity of KIR, HLA Class I, and Their Interactions in Seven Populations of Sub-Saharan Africans.

HLA class I and KIR sequences were determined for Dogon, Fulani, and Baka populations of western Africa, Mbuti of central Africa, and Datooga, Iraqw, and Hadza of eastern Africa. Study of 162 individuals identified 134 HLA class I alleles (41 HLA-A, 60 HLA-B, and 33 HLA-C). Common to all populations are three HLA-C alleles (C1+C*07:01, C1+C*07:02, and C2+C*06:02) but no HLA-A or -B Unexpectedly, no novel HLA class I was identified in these previously unstudied and anthropologically distinctive populations. In contrast, of 227 KIR detected, 22 are present in all seven populations and 28 are novel. A high diversity of HLA A-C-B haplotypes was observed. In six populations, most haplotypes are represented just once. But in the Hadza, a majority of haplotypes occur more than once, with 2 having high frequencies and 10 having intermediate frequencies. The centromeric (cen) part of the KIR locus exhibits an even balance between cenA and cenB in all seven populations. The telomeric (tel) part has an even balance of telA to telB in East Africa, but this changes across the continent to where telB is vestigial in West Africa. All four KIR ligands (A3/11, Bw4, C1, and C2) are present in six of the populations. HLA haplotypes of the Iraqw and Hadza encode two KIR ligands, whereas the other populations have an even balance between haplotypes encoding one and two KIR ligands. Individuals in these African populations have a mean of 6.8-8.4 different interactions between KIR and HLA class I, compared with 2.9-6.5 for non-Africans.

The genetic ancestry of African Americans, Latinos, and European Americans across the United States.

Over the past 500 years, North America has been the site of ongoing mixing of Native Americans, European settlers, and Africans (brought largely by the trans-Atlantic slave trade), shaping the early history of what became the United States. We studied the genetic ancestry of 5,269 self-described African Americans, 8,663 Latinos, and 148,789 European Americans who are 23andMe customers and show that the legacy of these historical interactions is visible in the genetic ancestry of present-day Americans. We document pervasive mixed ancestry and asymmetrical male and female ancestry contributions in all groups studied. We show that regional ancestry differences reflect historical events, such as early Spanish colonization, waves of immigration from many regions of Europe, and forced relocation of Native Americans within the US. This study sheds light on the fine-scale differences in ancestry within and across the United States and informs our understanding of the relationship between racial and ethnic identities and genetic ancestry.

Germ line variants predispose to both JAK2 V617F clonal hematopoiesis and myeloproliferative neoplasms.

We conducted a genome-wide association study (GWAS) to identify novel predisposition alleles associated with Philadelphia chromosome-negative myeloproliferative neoplasms (MPNs) and JAK2 V617F clonal hematopoiesis in the general population. We recruited a web-based cohort of 726 individuals with polycythemia vera, essential thrombocythemia, and myelofibrosis and 252 637 population controls unselected for hematologic phenotypes. Using a single-nucleotide polymorphism (SNP) array platform with custom probes for the JAK2 V617F mutation (V617F), we identified 497 individuals (0.2%) among the population controls who were V617F carriers. We performed a combined GWAS of the MPN cases plus V617F carriers in the control population (n = 1223) vs the remaining controls who were noncarriers for V617F (n = 252 140). For these MPN cases plus V617F carriers, we replicated the germ line JAK2 46/1 haplotype (rs59384377: odds ratio [OR] = 2.4, P = 6.6 × 10(-89)), previously associated with V617F-positive MPN. We also identified genome-wide significant associations in the TERT gene (rs7705526: OR = 1.8, P = 1.1 × 10(-32)), in SH2B3 (rs7310615: OR = 1.4, P = 3.1 × 10(-14)), and upstream of TET2 (rs1548483: OR = 2.0, P = 2.0 × 10(-9)). These associations were confirmed in a separate replication cohort of 446 V617F carriers vs 169 021 noncarriers. In a joint analysis of the combined GWAS and replication results, we identified additional genome-wide significant predisposition alleles associated with CHEK2, ATM, PINT, and GFI1B All SNP ORs were similar for MPN patients and controls who were V617F carriers. These data indicate that the same germ line variants endow individuals with a predisposition not only to MPN, but also to JAK2 V617F clonal hematopoiesis, a more common phenomenon that may foreshadow the development of an overt neoplasm.

Genome-wide analysis points to roles for extracellular matrix remodeling, the visual cycle, and neuronal development in myopia.

Myopia, or nearsightedness, is the most common eye disorder, resulting primarily from excess elongation of the eye. The etiology of myopia, although known to be complex, is poorly understood. Here we report the largest ever genome-wide association study (45,771 participants) on myopia in Europeans. We performed a survival analysis on age of myopia onset and identified 22 significant associations ([Formula: see text]), two of which are replications of earlier associations with refractive error. Ten of the 20 novel associations identified replicate in a separate cohort of 8,323 participants who reported if they had developed myopia before age 10. These 22 associations in total explain 2.9% of the variance in myopia age of onset and point toward a number of different mechanisms behind the development of myopia. One association is in the gene PRSS56, which has previously been linked to abnormally small eyes; one is in a gene that forms part of the extracellular matrix (LAMA2); two are in or near genes involved in the regeneration of 11-cis-retinal (RGR and RDH5); two are near genes known to be involved in the growth and guidance of retinal ganglion cells (ZIC2, SFRP1); and five are in or near genes involved in neuronal signaling or development. These novel findings point toward multiple genetic factors involved in the development of myopia and suggest that complex interactions between extracellular matrix remodeling, neuronal development, and visual signals from the retina may underlie the development of myopia in humans.

Replicability and robustness of genome-wide-association studies for behavioral traits.

A recent genome-wide-association study of educational attainment identified three single-nucleotide polymorphisms (SNPs) whose associations, despite their small effect sizes (each R (2) ≈ 0.02%), reached genome-wide significance (p < 5 × 10(-8)) in a large discovery sample and were replicated in an independent sample (p < .05). The study also reported associations between educational attainment and indices of SNPs called "polygenic scores." In three studies, we evaluated the robustness of these findings. Study 1 showed that the associations with all three SNPs were replicated in another large (N = 34,428) independent sample. We also found that the scores remained predictive (R (2) ≈ 2%) in regressions with stringent controls for stratification (Study 2) and in new within-family analyses (Study 3). Our results show that large and therefore well-powered genome-wide-association studies can identify replicable genetic associations with behavioral traits. The small effect sizes of individual SNPs are likely to be a major contributing factor explaining the striking contrast between our results and the disappointing replication record of most candidate-gene studies.

Rapid, global demographic expansions after the origins of agriculture.

The invention of agriculture is widely assumed to have driven recent human population growth. However, direct genetic evidence for population growth after independent agricultural origins has been elusive. We estimated population sizes through time from a set of globally distributed whole mitochondrial genomes, after separating lineages associated with agricultural populations from those associated with hunter-gatherers. The coalescent-based analysis revealed strong evidence for distinct demographic expansions in Europe, southeastern Asia, and sub-Saharan Africa within the past 10,000 y. Estimates of the timing of population growth based on genetic data correspond neatly to dates for the initial origins of agriculture derived from archaeological evidence. Comparisons of rates of population growth through time reveal that the invention of agriculture facilitated a fivefold increase in population growth relative to more ancient expansions of hunter-gatherers.

Web-based genome-wide association study identifies two novel loci and a substantial genetic component for Parkinson's disease.

Although the causes of Parkinson's disease (PD) are thought to be primarily environmental, recent studies suggest that a number of genes influence susceptibility. Using targeted case recruitment and online survey instruments, we conducted the largest case-control genome-wide association study (GWAS) of PD based on a single collection of individuals to date (3,426 cases and 29,624 controls). We discovered two novel, genome-wide significant associations with PD-rs6812193 near SCARB2 (p = 7.6 × 10(-10), OR = 0.84) and rs11868035 near SREBF1/RAI1 (p = 5.6 × 10(-8), OR = 0.85)-both replicated in an independent cohort. We also replicated 20 previously discovered genetic associations (including LRRK2, GBA, SNCA, MAPT, GAK, and the HLA region), providing support for our novel study design. Relying on a recently proposed method based on genome-wide sharing estimates between distantly related individuals, we estimated the heritability of PD to be at least 0.27. Finally, using sparse regression techniques, we constructed predictive models that account for 6%-7% of the total variance in liability and that suggest the presence of true associations just beyond genome-wide significance, as confirmed through both internal and external cross-validation. These results indicate a substantial, but by no means total, contribution of genetics underlying susceptibility to both early-onset and late-onset PD, suggesting that, despite the novel associations discovered here and elsewhere, the majority of the genetic component for Parkinson's disease remains to be discovered.

Hunter-gatherer genomic diversity suggests a southern African origin for modern humans.

Africa is inferred to be the continent of origin for all modern human populations, but the details of human prehistory and evolution in Africa remain largely obscure owing to the complex histories of hundreds of distinct populations. We present data for more than 580,000 SNPs for several hunter-gatherer populations: the Hadza and Sandawe of Tanzania, and the ≠Khomani Bushmen of South Africa, including speakers of the nearly extinct N|u language. We find that African hunter-gatherer populations today remain highly differentiated, encompassing major components of variation that are not found in other African populations. Hunter-gatherer populations also tend to have the lowest levels of genome-wide linkage disequilibrium among 27 African populations. We analyzed geographic patterns of linkage disequilibrium and population differentiation, as measured by F(ST), in Africa. The observed patterns are consistent with an origin of modern humans in southern Africa rather than eastern Africa, as is generally assumed. Additionally, genetic variation in African hunter-gatherer populations has been significantly affected by interaction with farmers and herders over the past 5,000 y, through both severe population bottlenecks and sex-biased migration. However, African hunter-gatherer populations continue to maintain the highest levels of genetic diversity in the world.

Assessing genetic contributions to phenotypic differences among 'racial' and 'ethnic' groups.

Descriptions of human genetic variation given thirty years ago have held up well, considering the substantial accrual of DNA sequence data in the interim. Most importantly, estimates of between-group genetic variation have remained relatively low. Despite the low average level of between-group variation, clusters recently inferred from multilocus genetic data coincide closely with groups defined by self-identified race or continental ancestry. This correspondence implies that genetic factors might contribute to unexplained between-group phenotypic variation. Current understanding of the contribution of genes to variation in most complex traits is limited, however. Under these circumstances, assumptions about genetic contributions to group differences are unfounded. In the absence of detailed understanding, 'racial' and 'ethnic' categories will remain useful in biomedical research. Further, we suggest approaches and guidelines for assessing the contribution of genetic factors to between-group phenotypic differences, including studies of candidate genes and analyses of recently admixed populations.

The genetic legacy of African Americans from Catoctin Furnace.

Few African Americans have been able to trace family lineages back to ancestors who died before the 1870 United States Census, the first in which all Black people were listed by name. We analyzed 27 individuals from Maryland's Catoctin Furnace African American Cemetery (1774-1850), identifying 41,799 genetic relatives among consenting research participants in 23andMe, Inc.'s genetic database. One of the highest concentrations of close relatives is in Maryland, suggesting that descendants of the Catoctin individuals remain in the area. We find that many of the Catoctin individuals derived African ancestry from the Wolof or Kongo groups and European ancestry from Great Britain and Ireland. This study demonstrates the power of joint analysis of historical DNA and large datasets generated through direct-to-consumer ancestry testing.

Navigating a research partnership between academia and industry to assess the impact of personalized genetic testing.

PURPOSE: To describe the process of structuring a partnership between academic researchers and two personalized genetic testing companies that would manage conflicts of interest while allowing researchers to study the impact of this nascent industry. METHODS: We developed a transparent process of ongoing communication about the interests of all research partners to address challenges in establishing study goals, survey development, data collection, analysis, and manuscript preparation. Using the existing literature on conflicts of interest and our experience, we created a checklist for academic and industry researchers seeking to structure research partnerships. RESULTS: Our checklist includes questions about the risk to research participants, sponsorship of the study, control of data analysis, freedom to publish results, the impact of the research on industry customers, openness to input from all partners, sharing results before publication, and publication of industry-specific data. Transparency is critical to building trust between partners. Involving all partners in the research development enhanced the quality of our research and provided an opportunity to manage conflicts early in the research process. CONCLUSION: Navigating relationships between academia and industry is complex and requires strategies that are transparent and responsive to the concerns of all. Employing a checklist of questions prior to beginning a research partnership may help to manage conflicts of interest.

Genetic variants associated with breast size also influence breast cancer risk.

BACKGROUND: While some factors of breast morphology, such as density, are directly implicated in breast cancer, the relationship between breast size and cancer is less clear. Breast size is moderately heritable, yet the genetic variants leading to differences in breast size have not been identified. METHODS: To investigate the genetic factors underlying breast size, we conducted a genome-wide association study (GWAS) of self-reported bra cup size, controlling for age, genetic ancestry, breast surgeries, pregnancy history and bra band size, in a cohort of 16,175 women of European ancestry. RESULTS: We identified seven single-nucleotide polymorphisms (SNPs) significantly associated with breast size (p<5.10(-8)): rs7816345 near ZNF703, rs4849887 and (independently) rs17625845 flanking INHBB, rs12173570 near ESR1, rs7089814 in ZNF365, rs12371778 near PTHLH, and rs62314947 near AREG. Two of these seven SNPs are in linkage disequilibrium (LD) with SNPs associated with breast cancer (those near ESR1 and PTHLH), and a third (ZNF365) is near, but not in LD with, a breast cancer SNP. The other three loci (ZNF703, INHBB, and AREG) have strong links to breast cancer, estrogen regulation, and breast development. CONCLUSIONS: These results provide insight into the genetic factors underlying normal breast development and show that some of these factors are shared with breast cancer. While these results do not directly support any possible epidemiological relationships between breast size and cancer, this study may contribute to a better understanding of the subtle interactions between breast morphology and breast cancer risk.

Y-chromosomal evidence of a pastoralist migration through Tanzania to southern Africa.

Although geneticists have extensively debated the mode by which agriculture diffused from the Near East to Europe, they have not directly examined similar agropastoral diffusions in Africa. It is unclear, for example, whether early instances of sheep, cows, pottery, and other traits of the pastoralist package were transmitted to southern Africa by demic or cultural diffusion. Here, we report a newly discovered Y-chromosome-specific polymorphism that defines haplogroup E3b1f-M293. This polymorphism reveals the monophyletic relationship of the majority of haplotypes of a previously paraphyletic clade, E3b1-M35*, that is widespread in Africa and southern Europe. To elucidate the history of the E3b1f haplogroup, we analyzed this haplogroup in 13 populations from southern and eastern Africa. The geographic distribution of the E3b1f haplogroup, in association with the microsatellite diversity estimates for populations, is consistent with an expansion through Tanzania to southern-central Africa. The data suggest this dispersal was independent of the migration of Bantu-speaking peoples along a similar route. Instead, the phylogeography and microsatellite diversity of the E3b1f lineage correlate with the arrival of the pastoralist economy in southern Africa. Our Y-chromosomal evidence supports a demic diffusion model of pastoralism from eastern to southern Africa approximately 2,000 years ago.

Characterizing the time dependency of human mitochondrial DNA mutation rate estimates.

Previous research has established a discrepancy of nearly an order of magnitude between pedigree-based and phylogeny-based (human vs. chimpanzee) estimates of the mitochondrial DNA (mtDNA) control region mutation rate. We characterize the time dependency of the human mitochondrial hypervariable region one mutation rate by generating 14 new phylogeny-based mutation rate estimates using within-human comparisons and archaeological dates. Rate estimates based on population events between 15,000 and 50,000 years ago are at least 2-fold lower than pedigree-based estimates. These within-human estimates are also higher than estimates generated from phylogeny-based human-chimpanzee comparisons. Our new estimates establish a rapid decay in evolutionary mutation rate between approximately 2,500 and 50,000 years ago and a slow decay from 50,000 to 6 Ma. We then extend this analysis to the mtDNA-coding region. Our within-human coding region mutation rate estimates display a similar, though less rapid, time-dependent decay. We explore the possibility that multiple hits explain the discrepancy between pedigree-based and phylogeny-based mutation rates. We conclude that whereas nucleotide substitution models incorporating multiple hits do provide a possible explanation for the discrepancy between pedigree-based and human-chimpanzee mutation rate estimates, they do not explain the rapid decline of within-human rate estimates. We propose that demographic processes such as serial bottlenecks prior to the Holocene could explain the difference between rates estimated before and after 15,000 years ago. Our findings suggest that human mtDNA estimates of dates of population and phylogenetic events should be adjusted in light of this time dependency of the mutation rate estimates.

Direct-to-consumer genetic testing for factor V Leiden and prothrombin 20210G>A: the consumer experience.

BACKGROUND: Clinical genetic testing for inherited predisposition to venous thromboembolism (VTE) is common among patients and their families. However, there is incomplete consensus about which individuals should receive testing, and the relative risks and benefits. METHODS: We assessed outcomes of receiving direct-to-consumer (DTC) results for the two most common genetic risk factors for VTE, factor V Leiden in the F5 gene (FVL) and prothrombin 20210G>A in the F2 gene (PT). Two thousand three hundred fifty-four customers (1244 variant-positive and 1110 variant-negative individuals) of the personal genetics company 23andMe, Inc., who had received results online for F5 and F2 variants, participated in an online survey-based study. Participants responded to questions about perception of VTE risk, discussion of results with healthcare providers (HCPs) and recommendations received, actions taken to control risk, emotional responses to receiving risk results, and perceived value of the information. RESULTS: Most participants (90% of variant-positive individuals, 99% of variant-negative individuals) had not previously been tested for F5 and/or F2 variants. The majority of variant-positive individuals correctly perceived that they were at higher than average risk for developing VTE. These individuals reported moderate rates of discussing results with HCPs (41%); receiving prevention advice from HCPs (31%), and making behavioral changes to control risk (e.g., exercising more, 30%). A minority (36%) of variant-positive individuals worried more after receiving VTE results. Nevertheless, most participants reported that knowing their risk had been an advantage (78% variant-positive and 58% variant-negative) and were satisfied knowing their genetic probability for VTE (81% variant-positive and 67% variant-negative). CONCLUSION: Consumers reported moderate rates of behavioral change and perceived personal benefit from receiving DTC genetic results for VTE risk.

REJECTOR: software for population history inference from genetic data via a rejection algorithm.

UNLABELLED: We introduce REJECTOR, software for parameter estimation and comparison of alternate models of population history from genetic data via a rejection algorithm. Through coalescent simulation, REJECTOR generates numerous gene genealogies, and hence simulated data, under a model of population history specified by the user. Summary statistics derived from such simulated data are compared with observed statistics, leading to acceptance or rejection of a given set of parameter values. We performed tests of the software using known parameter values in order to assess the inferential power provided by each summary statistic. The tests demonstrate the precision and accuracy of estimation made possible using this approach. AVAILABILITY: http://www.rejector.org