Comparison of sequencing data processing pipelines and application to underrepresented African human populations.

BACKGROUND: Population genetic studies of humans make increasing use of high-throughput sequencing in order to capture diversity in an unbiased way. There is an abundance of sequencing technologies, bioinformatic tools and the available genomes are increasing in number. Studies have evaluated and compared some of these technologies and tools, such as the Genome Analysis Toolkit (GATK) and its "Best Practices" bioinformatic pipelines. However, studies often focus on a few genomes of Eurasian origin in order to detect technical issues. We instead surveyed the use of the GATK tools and established a pipeline for processing high coverage full genomes from a diverse set of populations, including Sub-Saharan African groups, in order to reveal challenges from human diversity and stratification. RESULTS: We surveyed 29 studies using high-throughput sequencing data, and compared their strategies for data pre-processing and variant calling. We found that processing of data is very variable across studies and that the GATK "Best Practices" are seldom followed strictly. We then compared three versions of a GATK pipeline, differing in the inclusion of an indel realignment step and with a modification of the base quality score recalibration step. We applied the pipelines on a diverse set of 28 individuals. We compared the pipelines in terms of count of called variants and overlap of the callsets. We found that the pipelines resulted in similar callsets, in particular after callset filtering. We also ran one of the pipelines on a larger dataset of 179 individuals. We noted that including more individuals at the joint genotyping step resulted in different counts of variants. At the individual level, we observed that the average genome coverage was correlated to the number of variants called. CONCLUSIONS: We conclude that applying the GATK "Best Practices" pipeline, including their recommended reference datasets, to underrepresented populations does not lead to a decrease in the number of called variants compared to alternative pipelines. We recommend to aim for coverage of > 30X if identifying most variants is important, and to work with large sample sizes at the variant calling stage, also for underrepresented individuals and populations.

Estimating divergence times from DNA sequences.

The patterns of genetic variation within and among individuals and populations can be used to make inferences about the evolutionary forces that generated those patterns. Numerous population genetic approaches have been developed in order to infer evolutionary history. Here, we present the "Two-Two (TT)" and the "Two-Two-outgroup (TTo)" methods; two closely related approaches for estimating divergence time based in coalescent theory. They rely on sequence data from two haploid genomes (or a single diploid individual) from each of two populations. Under a simple population-divergence model, we derive the probabilities of the possible sample configurations. These probabilities form a set of equations that can be solved to obtain estimates of the model parameters, including population split times, directly from the sequence data. This transparent and computationally efficient approach to infer population divergence time makes it possible to estimate time scaled in generations (assuming a mutation rate), and not as a compound parameter of genetic drift. Using simulations under a range of demographic scenarios, we show that the method is relatively robust to migration and that the TTo method can alleviate biases that can appear from drastic ancestral population size changes. We illustrate the utility of the approaches with some examples, including estimating split times for pairs of human populations as well as providing further evidence for the complex relationship among Neandertals and Denisovans and their ancestors.

Multiple migrations to the Philippines during the last 50,000 years.

Island Southeast Asia has recently produced several surprises regarding human history, but the region's complex demography remains poorly understood. Here, we report ∼2.3 million genotypes from 1,028 individuals representing 115 indigenous Philippine populations and genome-sequence data from two ∼8,000-y-old individuals from Liangdao in the Taiwan Strait. We show that the Philippine islands were populated by at least five waves of human migration: initially by Northern and Southern Negritos (distantly related to Australian and Papuan groups), followed by Manobo, Sama, Papuan, and Cordilleran-related populations. The ancestors of Cordillerans diverged from indigenous peoples of Taiwan at least ∼8,000 y ago, prior to the arrival of paddy field rice agriculture in the Philippines ∼2,500 y ago, where some of their descendants remain to be the least admixed East Asian groups carrying an ancestry shared by all Austronesian-speaking populations. These observations contradict an exclusive "out-of-Taiwan" model of farming-language-people dispersal within the last four millennia for the Philippines and Island Southeast Asia. Sama-related ethnic groups of southwestern Philippines additionally experienced some minimal South Asian gene flow starting ∼1,000 y ago. Lastly, only a few lowlanders, accounting for <1% of all individuals, presented a low level of West Eurasian admixture, indicating a limited genetic legacy of Spanish colonization in the Philippines. Altogether, our findings reveal a multilayered history of the Philippines, which served as a crucial gateway for the movement of people that ultimately changed the genetic landscape of the Asia-Pacific region.

The deep population history in Africa.

Africa is the continent with the greatest genetic diversity among humans and the level of diversity is further enhanced by incorporating non-majority groups, which are often understudied. Many of today's minority populations historically practiced foraging lifestyles, which were the only subsistence strategies prior to the rise of agriculture and pastoralism, but only a few groups practicing these strategies remain today. Genomic investigations of Holocene human remains excavated across the African continent show that the genetic landscape was vastly different compared to today's genetic landscape and that many groups that today are population isolate inhabited larger regions in the past. It is becoming clear that there are periods of isolation among groups and geographic areas, but also genetic contact over large distances throughout human history in Africa. Genomic information from minority populations and from prehistoric remains provide an invaluable source of information on the human past, in particular deep human population history, as Holocene large-scale population movements obscure past patterns of population structure. Here we revisit questions on the nature and time of the radiation of early humans in Africa, the extent of gene-flow among human populations as well as introgression from archaic and extinct lineages on the continent.

Taste perception and lifestyle: insights from phenotype and genome data among Africans and Asians.

Taste is essential for the interaction of animals with their food and has co-evolved with diet. Humans have peopled a large range of environments and present a wide range of diets, but little is known about the diversity and evolution of human taste perception. We measured taste recognition thresholds across populations differing in lifestyles (hunter gatherers and farmers from Central Africa, nomad herders, and farmers from Central Asia). We also generated genome-wide genotype data and performed association studies and selection scans in order to link the phenotypic variation in taste sensitivity with genetic variation. We found that hunter gatherers have lower overall sensitivity as well as lower sensitivity to quinine and fructose than their farming neighbors. In parallel, there is strong population divergence in genes associated with tongue morphogenesis and genes involved in the transduction pathway of taste signals in the African populations. We find signals of recent selection in bitter taste-receptor genes for all four populations. Enrichment analysis on association scans for the various tastes confirmed already documented associations and revealed novel GO terms that are good candidates for being involved in taste perception. Our framework permitted us to gain insight into the genetic basis of taste sensitivity variation across populations and lifestyles.

Khoe-San Genomes Reveal Unique Variation and Confirm the Deepest Population Divergence in Homo sapiens.

The southern African indigenous Khoe-San populations harbor the most divergent lineages of all living peoples. Exploring their genomes is key to understanding deep human history. We sequenced 25 full genomes from five Khoe-San populations, revealing many novel variants, that 25% of variants are unique to the Khoe-San, and that the Khoe-San group harbors the greatest level of diversity across the globe. In line with previous studies, we found several gene regions with extreme values in genome-wide scans for selection, potentially caused by natural selection in the lineage leading to Homo sapiens and more recent in time. These gene regions included immunity-, sperm-, brain-, diet-, and muscle-related genes. When accounting for recent admixture, all Khoe-San groups display genetic diversity approaching the levels in other African groups and a reduction in effective population size starting around 100,000 years ago. Hence, all human groups show a reduction in effective population size commencing around the time of the Out-of-Africa migrations, which coincides with changes in the paleoclimate records, changes that potentially impacted all humans at the time.

Patterns of variation in cis-regulatory regions: examining evidence of purifying selection.

BACKGROUND: With only 2 % of the human genome consisting of protein coding genes, functionality across the rest of the genome has been the subject of much debate. This has gained further impetus in recent years due to a rapidly growing catalogue of genomic elements, based primarily on biochemical signatures (e.g. the ENCODE project). While the assessment of functionality is a complex task, the presence of selection acting on a genomic region is a strong indicator of importance. In this study, we apply population genetic methods to investigate signals overlaying several classes of regulatory elements. RESULTS: We disentangle signals of purifying selection acting directly on regulatory elements from the confounding factors of demography and purifying selection linked to e.g. nearby protein coding regions. We confirm the importance of regulatory regions proximal to coding sequence, while also finding differential levels of selection at distal regions. We note differences in purifying selection among transcription factor families. Signals of constraint at some genomic classes were also strongly dependent on their physical location relative to coding sequence. In addition, levels of selection efficacy across genomic classes differed between African and non-African populations. CONCLUSIONS: In order to assign a valid signal of selection to a particular class of genomic sequence, we show that it is crucial to isolate the signal by accounting for the effects of demography and linked-purifying selection. Our study highlights the intricate interplay of factors affecting signals of selection on functional elements.

Southern African ancient genomes estimate modern human divergence to 350,000 to 260,000 years ago.

Southern Africa is consistently placed as a potential region for the evolution of Homo sapiens We present genome sequences, up to 13x coverage, from seven ancient individuals from KwaZulu-Natal, South Africa. The remains of three Stone Age hunter-gatherers (about 2000 years old) were genetically similar to current-day southern San groups, and those of four Iron Age farmers (300 to 500 years old) were genetically similar to present-day Bantu-language speakers. We estimate that all modern-day Khoe-San groups have been influenced by 9 to 30% genetic admixture from East Africans/Eurasians. Using traditional and new approaches, we estimate the first modern human population divergence time to between 350,000 and 260,000 years ago. This estimate increases the deepest divergence among modern humans, coinciding with anatomical developments of archaic humans into modern humans, as represented in the local fossil record.

Adaptation to infectious disease exposure in indigenous Southern African populations.

Genetic analyses can provide information about human evolutionary history that cannot always be gleaned from other sources. We evaluated evidence of selective pressure due to introduced infectious diseases in the genomes of two indigenous southern African San groups-the ‡Khomani who had abundant contact with other people migrating into the region and the more isolated Ju|'hoansi. We used a dual approach to test for increased selection on immune genes compared with the rest of the genome in these groups. First, we calculated summary values of statistics that measure genomic signatures of adaptation to contrast selection signatures in immune genes and all genes. Second, we located regions of the genome with extreme values of three selection statistics and examined these regions for enrichment of immune genes. We found stronger and more abundant signals of selection in immune genes in the ‡Khomani than in the Ju|'hoansi. We confirm this finding within each population to avoid effects of different demographic histories of the two populations. We identified eight immune genes that have potentially been targets of strong selection in the ‡Khomani, whereas in the Ju|'hoansi, no immune genes were found in the genomic regions with the strongest signals of selection. We suggest that the more abundant signatures of selection at immune genes in the ‡Khomani could be explained by their more frequent contact with immigrant groups, which likely led to increased exposure and adaptation to introduced infectious diseases.

Genome-wide transcriptional and physiological responses to drought stress in leaves and roots of two willow genotypes.

BACKGROUND: Drought is a major environmental stress that can have severe impacts on plant productivity and survival. Understanding molecular mechanisms of drought responses is crucial in order to breed for drought adapted plant cultivars. The aim of the present study was to investigate phenotypic and transcriptional drought responses in two willow genotypes (520 and 592) originating from an experimental cross between S. viminalis × (S. viminalis × S. schwerinii). Willows are woody perennials in the Salicaceae plant family that are grown as bioenergy crops worldwide. METHODS: An experiment was conducted where plants were exposed to drought and different eco-physiological parameters were assessed. RNA-seq data was furthermore generated with the Illumina technology from root tips and leaves from plants grown in drought and well-watered (WW) conditions. The RNA-seq data was assembled de novo with the Trinity assembler to create a reference gene set to which the reads were mapped in order to obtain differentially expressed genes (DEGs) between the drought and WW conditions. To investigate molecular mechanisms involved in the drought response, GO enrichment analyses were conducted. Candidate genes with a putative function in the drought response were also identified. RESULTS: A total of 52,599 gene models were obtained and after filtering on gene expression (FPKM ≥ 1), 35,733 gene models remained, of which 24,421 contained open reading frames. A total of 5,112 unique DEGs were identified between drought and WW conditions, of which the majority were found in the root tips. Phenotypically, genotype 592 displayed less growth reduction in response to drought compared to genotype 520. At the transcriptional level, genotype 520 displayed a greater response in the leaves as more DEGs were found in genotype 520 compared to genotype 592. In contrast, the transcriptional responses in the root tips were rather similar between the two genotypes. A core set of candidate genes encoding proteins with a putative function in drought response was identified, for example MYBs and bZIPs as well as chlorophyll a/b binding proteins. DISCUSSION: We found substantial differences in drought responses between the genotypes, both at the phenotypic and transcriptional levels. In addition to the genotypic variation in several traits, we also found indications for genotypic variation in trait plasticity, which could play a role in drought adaptation. Furthermore, the two genotypes displayed overall similar transcriptional responses in the root tips, but more variation in the leaves. It is thus possible that the observed phenotypic differences could be a result of transcriptional differences mostly at the leaf level. CONCLUSIONS: This study has contributed to a better general understanding of drought responses in woody plants, specifically in willows, and has implications for breeding research towards more drought adapted plants.

Ancient mitochondrial DNA from the northern fringe of the Neolithic farming expansion in Europe sheds light on the dispersion process.

The European Neolithization process started around 12 000 years ago in the Near East. The introduction of agriculture spread north and west throughout Europe and a key question has been if this was brought about by migrating individuals, by an exchange of ideas or a by a mixture of these. The earliest farming evidence in Scandinavia is found within the Funnel Beaker Culture complex (Trichterbecherkultur, TRB) which represents the northernmost extension of Neolithic farmers in Europe. The TRB coexisted for almost a millennium with hunter-gatherers of the Pitted Ware Cultural complex (PWC). If migration was a substantial part of the Neolithization, even the northerly TRB community would display a closer genetic affinity to other farmer populations than to hunter-gatherer populations. We deep-sequenced the mitochondrial hypervariable region 1 from seven farmers (six TRB and one Battle Axe complex, BAC) and 13 hunter-gatherers (PWC) and authenticated the sequences using postmortem DNA damage patterns. A comparison with 124 previously published sequences from prehistoric Europe shows that the TRB individuals share a close affinity to Central European farmer populations, and that they are distinct from hunter-gatherer groups, including the geographically close and partially contemporary PWC that show a close affinity to the European Mesolithic hunter-gatherers.

Investigating population history using temporal genetic differentiation.

The rapid advance of sequencing technology, coupled with improvements in molecular methods for obtaining genetic data from ancient sources, holds the promise of producing a wealth of genomic data from time-separated individuals. However, the population-genetic properties of time-structured samples have not been extensively explored. Here, we consider the implications of temporal sampling for analyses of genetic differentiation and use a temporal coalescent framework to show that complex historical events such as size reductions, population replacements, and transient genetic barriers between populations leave a footprint of genetic differentiation that can be traced through history using temporal samples. Our results emphasize explicit consideration of the temporal structure when making inferences and indicate that genomic data from ancient individuals will greatly increase our ability to reconstruct population history.

Private haplotypes can reveal local adaptation.

BACKGROUND: Genome-wide scans for regions that demonstrate deviating patterns of genetic variation have become common approaches for finding genes targeted by selection. Several genomic patterns have been utilized for this purpose, including deviations in haplotype homozygosity, frequency spectra and genetic differentiation between populations. RESULTS: We describe a novel approach based on the Maximum Frequency of Private Haplotypes--MFPH--to search for signals of recent population-specific selection. The MFPH statistic is straightforward to compute for phased SNP- and sequence-data. Using both simulated and empirical data, we show that MFPH can be a powerful statistic to detect recent population-specific selection, that it performs at the same level as other commonly used summary statistics (e.g. FST, iHS and XP-EHH), and that MFPH in some cases capture signals of selection that are missed by other statistics. For instance, in the Maasai, MFPH reveals a strong signal of selection in a region where other investigated statistics fail to pick up a clear signal that contains the genes DOCK3, MAPKAPK3 and CISH. This region has been suggested to affect height in many populations based on phenotype-genotype association studies. It has specifically been suggested to be targeted by selection in Pygmy groups, which are on the opposite end of the human height spectrum compared to the Maasai. CONCLUSIONS: From the analysis of both simulated and publicly available empirical data, we show that MFPH represents a summary statistic that can provide further insight concerning population-specific adaptation.

Lactase persistence alleles reveal partial East African ancestry of southern African Khoe pastoralists.

The ability to digest milk into adulthood, lactase persistence (LP), as well as specific genetic variants associated with LP, is heterogeneously distributed in global populations. These variants were most likely targets of selection when some populations converted from hunter-gatherer to pastoralist or farming lifestyles. Specific LP polymorphisms are associated with particular geographic regions and populations; however, they have not been extensively studied in southern Africa. We investigate the LP-regulatory region in 267 individuals from 13 southern African populations (including descendants of hunter-gatherers, pastoralists, and agropastoralists), providing the first comprehensive study of the LP-regulatory region in a large group of southern Africans. The "East African" LP single-nucleotide polymorphism (SNP) (14010G>C) was found at high frequency (>20%) in a strict pastoralist Khoe population, the Nama of Namibia, suggesting a connection to East Africa, whereas the "European" LP SNP (13910C>T) was found in populations of mixed ancestry. Using genome-wide data from various African populations, we identify admixture (13%) in the Nama, from an Afro-Asiatic group dating to >1,300 years ago, with the remaining fraction of their genomes being from San hunter-gatherers. We also find evidence of selection around the LCT gene among Khoe-speaking groups, and the substantial frequency of the 14010C variant among the Nama is best explained by adaptation to digesting milk. These genome-local and genome-wide results support a model in which an East African group brought pastoralist practices to southern Africa and admixed with local hunter-gatherers to form the ancestors of Khoe people.

Assessing the maximum contribution from ancient populations.

Ancestral relationships between populations separated by time represent an often neglected dimension in population genetics, a field which historically has focused on analysis of spatially distributed samples from the same point in time. Models are usually straightforward when two time-separated populations are assumed to be completely isolated from all other populations. However, this is usually an unrealistically stringent assumption when there is gene flow with other populations. Here, we investigate continuity in the presence of gene flow from unknown populations. This setup allows a more nuanced treatment of questions regarding population continuity in terms of "level of contribution" from a particular ancient population to a more recent population. We propose a statistical framework which makes use of a biallelic marker sampled at two different points in time to assess population contribution, and present two different interpretations of the concept. We apply the approach to published data from a prehistoric human population in Scandinavia (Malmström H, Gilbert MTP, Thomas MG, Brandström M, Storå J, Molnar P, Andersen PK, Bendixen C, Holmlund G, Götherström A, et al. 2009. Ancient DNA reveals lack of continuity between Neolithic hunter-gatherers and contemporary Scandinavians. Curr Biol. 19:1758-1762) and Pleistocene woolly mammoth (Barnes I, Shapiro B, Lister A, Kuznetsova T, Sher A, Guthrie D, Thomas MG. 2007. Genetic structure and extinction of the woolly mammoth, Mammuthus primigenius. Curr Biol. 17:1072-1075; Debruyne R, Chu G, King CE, Bos K, Kuch M, Schwarz C, Szpak P, Gröcke DR, Matheus P, Zazula G, et al. 2008. Out of America: ancient DNA evidence for a new world origin of late quaternary woolly mammoths. Curr Biol. 18:1320-1326).

Stronger signal of recent selection for lactase persistence in Maasai than in Europeans.

Continued ability to digest lactose after weaning provides a possible selective advantage to individuals who have access to milk as a food source. The lactase persistence (LP) phenotype exists at varying frequencies in different populations and SNPs that modulate the regulation of the LCT gene have been identified in many of these populations. Very strong positive selection for LP has been illustrated for a single SNP (rs4988235) in northwestern European populations, which has become a textbook example of the effect of recent selective sweeps on genetic variation and linkage disequilibrium. In this study, we employed two different methods to detect signatures of positive selection in an East African pastoralist population in the HapMap collection, the Maasai from Kenya, and compared results with other HapMap populations. We found that signatures of recent selection coinciding with the LCT gene are the strongest across the genome in the Maasai population. Furthermore, the genome-wide signal of recent positive selection on haplotypic variation and population differentiation around the LCT gene is greater in the Maasai than in the CEU population (northwestern European descent), possibly due to stronger selection pressure, but it could also be an indication of more recent selection in Maasai compared with the Central European group or more efficient selection in the Maasai due to less genetic drift for their larger effective population size. This signal of recent selection is driven by a putative East African LP haplotype that is different from the haplotype that contributes to the LP phenotype in northwestern Europe.

Anisotropic isolation by distance: the main orientations of human genetic differentiation.

Genetic differentiation among human populations is greatly influenced by geography due to the accumulation of local allele frequency differences. However, little is known about the possibly different increment of genetic differentiation along the different geographical axes (north-south, east-west, etc.). Here, we provide new methods to examine the asymmetrical patterns of genetic differentiation. We analyzed genome-wide polymorphism data from populations in Africa (n = 29), Asia (n = 26), America (n = 9), and Europe (n = 38), and we found that the major orientations of genetic differentiation are north-south in Europe and Africa, and east-west in Asia, but no preferential orientation was found in the Americas. Additionally, we showed that the localization of the individual geographic origins based on single nucleotide polymorphism data was not equally precise along all orientations. Confirming our findings, we obtained that, in each continent, the orientation along which the precision is maximal corresponds to the orientation of maximum differentiation. Our results have implications for interpreting human genetic variation in terms of isolation by distance and spatial range expansion processes. In Europe, for instance, the precise northnorthwest-southsoutheast axis of main European differentiation cannot be explained by a simple Neolithic demic diffusion model without admixture with the local populations because in that case the orientation of greatest differentiation should be perpendicular to the direction of expansion. In addition to humans, anisotropic analyses can guide the description of genetic differentiation for other organisms and provide information on expansions of invasive species or the processes of plant dispersal.

Genomic variation in seven Khoe-San groups reveals adaptation and complex African history.

The history of click-speaking Khoe-San, and African populations in general, remains poorly understood. We genotyped ~2.3 million single-nucleotide polymorphisms in 220 southern Africans and found that the Khoe-San diverged from other populations ≥100,000 years ago, but population structure within the Khoe-San dated back to about 35,000 years ago. Genetic variation in various sub-Saharan populations did not localize the origin of modern humans to a single geographic region within Africa; instead, it indicated a history of admixture and stratification. We found evidence of adaptation targeting muscle function and immune response; potential adaptive introgression of protection from ultraviolet light; and selection predating modern human diversification, involving skeletal and neurological development. These new findings illustrate the importance of African genomic diversity in understanding human evolutionary history.

Resequencing data provide no evidence for a human bottleneck in Africa during the penultimate glacial period.

Based on the accumulation of genetic, climatic, and fossil evidence, a central theory in paleoanthropology stipulates that a demographic bottleneck coincided with the origin of our species Homo Sapiens. This theory proposes that anatomically modern humans--which were only present in Africa at the time--experienced a drastic bottleneck during the penultimate glacial age (130-190 kya) when a cold and dry climate prevailed. Two scenarios have been proposed to describe the bottleneck, which involve either a fragmentation of the range occupied by humans or the survival of one small group of humans. Here, we analyze DNA sequence data from 61 nuclear loci sequenced in three African populations using Approximate Bayesian Computation and numerical simulations. In contrast to the bottleneck theory, we show that a simple model without any bottleneck during the penultimate ice age has the greatest statistical support compared with bottleneck models. Although the proposed bottleneck is ancient, occurring at least 130 kya, we can discard the possibility that it did not leave detectable footprints in the DNA sequence data except if the bottleneck involves a less than a 3-fold reduction in population size. Finally, we confirm that a simple model without a bottleneck is able to reproduce the main features of the observed patterns of genetic variation. We conclude that models of Pleistocene refugium for modern human origins now require substantial revision.