Genetic variation among world populations: inferences from 100 Alu insertion polymorphisms

A systematic literature review and meta-analysis of circulating nucleic acids as biomarkers in psychiatry

Common mental disorders (CMDs) such as depression, anxiety and post-traumatic stress disorders account for 40% of the global burden of disease. In most psychiatric disorders, both diagnosis and monitoring can be challenging, frequently requiring long-term investigation and follow-up. The discovery of better methods to facilitate accurate and fast diagnosis and monitoring of psychiatric disorders is therefore crucial. Circulating nucleic acids (CNAs) are among these new tools. CNAs (DNA or RNA) can be found circulating in body biofluids, and can be isolated from biological samples such as plasma. They can serve as biomarkers for diagnosis and prognoses. They appear to be promising for disorders (such as psychiatric disorders) that involve organs or structures that are difficult to assess. This review presents an accurate assessment of the current literature about the use of plasma and serum cell-free DNA (cfDNA) as biomarkers for several aspects of psychiatric disorders: diagnosis, prognosis, treatment response, and monitor disease progression. For each psychiatric disorder, we examine the effect sizes to give insights on the efficacy of CNAs as biomarkers. The global effect size for plasma nuclear and mitochondrial cfDNA studies was generally moderate for psychiatric disorders. In addition, we discuss future applications of CNAs and particularly cfDNA as non-invasive biomarkers for these diseases.

Genomic diversity and differentiation of Alu insertion polymorphisms in a native British and four South Asian migrant populations

BACKGROUND

Alu insertions are bi-allelic and primate-specific, this makes them a useful marker for studying genetic variation, migration patterns, forensic analyses, paternity, and evolutionary heritage; however, specific population studies are limited.

AIM

The objective of this study is to document the level and extent of genetic variation at 39 different Alu loci in five populations (British, Indian Punjabi, Indian Gujarati, Pakistani, and Bangladeshi) from the East Midlands region of the UK. Genetic data on migrant populations is currently limited.

SUBJECTS AND METHODS

DNA samples (n = 543) were analysed for 39 Alu insertion polymorphisms using specific primers and standard protocols. Data were analysed for population and forensic genetic parameters.

RESULTS

All studied Alus were polymorphic in the British White population while South Asian migrant populations had a variable number of loci which were monomorphic. Highest heterozygosities and lowest match probabilities were observed in the British sample, while the Bangladeshi sample had the lowest heterozygosity and higher match probability.

CONCLUSION

The analysed Alus insertions (TPA25, Ya5NBC123, Ya5NBC182, Ya5NBC241, and Ya5NBC242) are highly polymorphic and variable among migrant populations. These loci could be useful for population genomic and differentiation studies.

An Alu insertion map of the Indian population: identification and analysis in 1021 genomes of the IndiGen project

Actively retrotransposing primate-specific Alu repeats display insertion-deletion (InDel) polymorphism through their insertion at new loci. In the global datasets, Indian populations remain under-represented and so do their Alu InDels. Here, we report the genomic landscape of Alu InDels from the recently released 1021 Indian Genomes (IndiGen) (available at https://clingen.igib.res.in/indigen). We identified 9239 polymorphic Alu insertions that include private (3831), rare (3974) and common (1434) insertions with an average of 770 insertions per individual. We achieved an 89% PCR validation of the predicted genotypes in 94 samples tested. About 60% of identified InDels are unique to IndiGen when compared to other global datasets; 23% of sites were shared with both SGDP and HGSVC; among these, 58% (1289 sites) were common polymorphisms in IndiGen. The insertions not only show a bias for genic regions, with a preference for introns but also for the associated genes showing enrichment for processes like cell morphogenesis and neurogenesis (P-value < 0.05). Approximately, 60% of InDels mapped to genes present in the OMIM database. Finally, we show that 558 InDels can serve as ancestry informative markers to segregate global populations. This study provides a valuable resource for baseline Alu InDels that would be useful in population genomics.

The Simons Genome Diversity Project: A Global Analysis of Mobile Element Diversity

Ongoing retrotransposition of Alu, LINE-1, and SINE–VNTR–Alu elements generates diversity and variation among human populations. Previous analyses investigating the population genetics of mobile element insertions (MEIs) have been limited by population ascertainment bias or by relatively small numbers of populations and low sequencing coverage. Here, we use 296 individuals representing 142 global populations from the Simons Genome Diversity Project (SGDP) to discover and characterize MEI diversity from deeply sequenced whole-genome data. We report 5,742 MEIs not originally reported by the 1000 Genomes Project and show that high sampling diversity leads to a 4- to 7-fold increase in MEI discovery rates over the original 1000 Genomes Project data. As a result of negative selection, nonreference polymorphic MEIs are underrepresented within genes, and MEIs within genes are often found in the transcriptional orientation opposite that of the gene. Globally, 80% of Alu subfamilies predate the expansion of modern humans from Africa. Polymorphic MEIs show heterozygosity gradients that decrease from Africa to Eurasia to the Americas, and the number of MEIs found uniquely in a single individual are also distributed in this general pattern. The maximum fraction of MEI diversity partitioned among the seven major SGDP population groups (F_ST) is 7.4%, similar to, but slightly lower than, previous estimates and likely attributable to the diverse sampling strategy of the SGDP. Finally, we utilize these MEIs to extrapolate the primary Native American shared ancestry component to back to Asia and provide new evidence from genome-wide identical-by-descent genetic markers that add additional support for a southeastern Siberian origin for most Native Americans.

Mobile element insertions and associated structural variants in longitudinal breast cancer samples

While mobile elements are largely inactive in healthy somatic tissues, increased activity has been found in cancer tissues, with significant variation among different cancer types. In addition to insertion events, mobile elements have also been found to mediate many structural variation events in the genome. Here, to better understand the timing and impact of mobile element insertions and associated structural variants in cancer, we examined their activity in longitudinal samples of four metastatic breast cancer patients. We identified 11 mobile element insertions or associated structural variants and found that the majority of these occurred early in tumor progression. Most of the variants impact intergenic regions; however, we identified a translocation interrupting MAP2K4 involving Alu elements and a deletion in YTHDF2 involving mobile elements that likely inactivate reported tumor suppressor genes. The high variant allele fraction of the translocation, the loss of the other copy of MAP2K4, the recurrent loss-of-function mutations found in this gene in other cancers, and the important function of MAP2K4 indicate that this translocation is potentially a driver mutation. Overall, using a unique longitudinal dataset, we find that most variants are likely passenger mutations in the four patients we examined, but some variants impact tumor progression.

Alu retrotransposons and COVID-19 susceptibility and morbidity

SARS-CoV-2 has spread rapidly across the world and is negatively impacting the global human population. COVID-19 patients display a wide variety of symptoms and clinical outcomes, including those attributed to genetic ancestry. Alu retrotransposons have played an important role in human evolution, and their variants influence host response to viral infection. Intronic Alus regulate gene expression through several mechanisms, including both genetic and epigenetic pathways. With respect to SARS-CoV-2, an intronic Alu within the ACE gene is hypothesized to be associated with COVID-19 susceptibility and morbidity. Here, we review specific Alu polymorphisms that are of particular interest when considering host response to SARS-CoV-2 infection, especially polymorphic Alu insertions in genes associated with immune response and coagulation/fibrinolysis cascade. We posit that additional research focused on Alu-related pathways could yield novel biomarkers capable of predicting clinical outcomes as well as patient-specific treatment strategies for COVID-19 and related infectious diseases.

TypeTE: a tool to genotype mobile element insertions from whole genome resequencing data

Alu retrotransposons account for more than 10% of the human genome, and insertions of these elements create structural variants segregating in human populations. Such polymorphic Alus are powerful markers to understand population structure, and they represent variants that can greatly impact genome function, including gene expression. Accurate genotyping of Alus and other mobile elements has been challenging. Indeed, we found that Alu genotypes previously called for the 1000 Genomes Project are sometimes erroneous, which poses significant problems for phasing these insertions with other variants that comprise the haplotype. To ameliorate this issue, we introduce a new pipeline - TypeTE - which genotypes Alu insertions from whole-genome sequencing data. Starting from a list of polymorphic Alus, TypeTE identifies the hallmarks (poly-A tail and target site duplication) and orientation of Alu insertions using local re-assembly to reconstruct presence and absence alleles. Genotype likelihoods are then computed after re-mapping sequencing reads to the reconstructed alleles. Using a high-quality set of PCR-based genotyping of >200 loci, we show that TypeTE improves genotype accuracy from 83% to 92% in the 1000 Genomes dataset. TypeTE can be readily adapted to other retrotransposon families and brings a valuable toolbox addition for population genomics.

Sensitivity of the polyDetect computational pipeline for phylogenetic analyses

Alu elements are powerful phylogenetic markers. The combination of a recently-developed computational pipeline, polyDetect, with high copy number Alu insertions has previously been utilized to help resolve the Papio baboon phylogeny with high statistical support. Here, the polyDetect method was applied to the highly contentious Cebidae phylogeny within New World monkeys (NWM). The polyDetect method relies on conserved homology/identity of short read sequence data among the species being compared to accurately map predicted shared Alu insertions to each unique flanking sequence. The results of this comprehensive assessment indicate that there were insufficient sequence homology/identity stretches in non-repeated DNA sequences among the four Cebidae genera analyzed in this study to make this strategy phylogenetically viable. The ~20 million years of evolutionary divergence of the Cebidae genera has resulted in random sequence decay within the short read data, obscuring potentially orthologous elements in the species tested. These analyses suggest that the polyDetect pipeline is best suited to resolving phylogenies of more recently diverged lineages when high-quality assembled genomes are not available for the taxa of interest.

A Two-State Model of Tree Evolution and Its Applications to Alu Retrotransposition

Models of tree evolution have mostly focused on capturing the cladogenesis processes behind speciation. Processes that derive the evolution of genomic elements, such as repeats, are not necessarily captured by these existing models. In this article, we design a model of tree evolution that we call the dual-birth model, and we show how it can be useful in studying the evolution of short Alu repeats found in the human genome in abundance. The dual-birth model extends the traditional birth-only model to have two rates of propagation, one for active nodes that propagate often, and another for inactive nodes, that with a lower rate, activate and start propagating. Adjusting the ratio of the rates controls the expected tree balance. We present several theoretical results under the dual-birth model, introduce parameter estimation techniques, and study the properties of the model in simulations. We then use the dual-birth model to estimate the number of active Alu elements and their rates of propagation and activation in the human genome based on a large phylogenetic tree that we build from close to one million Alu sequences.

A computational reconstruction of Papio phylogeny using Alu insertion polymorphisms

Background

Since the completion of the human genome project, the diversity of genome sequencing data produced for non-human primates has increased exponentially. Papio baboons are well-established biological models for studying human biology and evolution. Despite substantial interest in the evolution of Papio, the systematics of these species has been widely debated, and the evolutionary history of Papio diversity is not fully understood. Alu elements are primate-specific transposable elements with a well-documented mutation/insertion mechanism and the capacity for resolving controversial phylogenetic relationships. In this study, we conducted a whole genome analysis of Alu insertion polymorphisms unique to the Papio lineage. To complete these analyses, we created a computational algorithm to identify novel Alu insertions in next-generation sequencing data.

Results

We identified 187,379 Alu insertions present in the Papio lineage, yet absent from M. mulatta [Mmul8.0.1]. These elements were characterized using genomic data sequenced from a panel of twelve Papio baboons: two from each of the six extant Papio species. These data were used to construct a whole genome Alu-based phylogeny of Papio baboons. The resulting cladogram fully-resolved relationships within Papio.

Conclusions

These data represent the most comprehensive Alu-based phylogenetic reconstruction reported to date. In addition, this study produces the first fully resolved Alu-based phylogeny of Papio baboons.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0118-3) contains supplementary material, which is available to authorized users.

Analysis of lineage-specific Alu subfamilies in the genome of the olive baboon, Papio anubis

Background

Alu elements are primate-specific retroposons that mobilize using the enzymatic machinery of L1 s. The recently completed baboon genome project found that the mobilization rate of Alu elements is higher than in the genome of any other primate studied thus far. However, the Alu subfamily structure present in and specific to baboons had not been examined yet.

Results

Here we report 129 Alu subfamilies that are propagating in the genome of the olive baboon, with 127 of these subfamilies being new and specific to the baboon lineage. We analyzed 233 Alu insertions in the genome of the olive baboon using locus specific polymerase chain reaction assays, covering 113 of the 129 subfamilies. The allele frequency data from these insertions show that none of the nine groups of subfamilies are nearing fixation in the lineage.

Conclusions

Many subfamilies of Alu elements are actively mobilizing throughout the baboon lineage, with most being specific to the baboon lineage.

Frequency of Alu insertions within the ACE and PR loci in Northwestern Mexicans

OBJECTIVE

Presently, non-LTR retrotransposons are the most active mobile elements in the human genome. Among these, Alu elements are highly represented in the modern population. Worldwide, distribution of Alu polymorphisms (insertion/deletion; I/D) shows variability between different populations. Two Alu insertion loci, ACE and PR, are significant biomarkers that have served in several genotype-phenotype association studies. In Mexico, studies concerning the frequency of these biomarkers have been conducted mainly in subpopulations from central and southern regions. Here, we screened a population sample of the northwestern region to gain further knowledge regarding the prevalence of Alu polymorphisms within ACE and PR loci.

RESULTS

For ACE locus, the observed genotype frequencies were 26.5, 51.0 and 22.5% for II, ID, and DD, respectively; and allelic frequencies for I and D were 52 and 48%. Whereas respective genotype frequencies for PR locus were 2.7, 26.5 and 70.8%, and the corresponding allele frequencies were 16 and 84%. Furthermore, the insertion frequency within ACE locus was similar between central, western and northwestern subpopulations, and rather higher in southeastern subpopulation (p < 0.05). Although the occurrence of Alu polymorphisms within PR locus has not been widely examined, the insertion frequency was higher in northwestern subpopulation, as compared with western and southeastern subpopulations (p < 0.05). Based on the frequency of Alu insertions found in ACE and PR loci, subpopulations from the northwestern, western and central regions share a common genetic origin, but apparently not with the subpopulation from the southeastern region, in accordance with the notion that assumes the existence of a broad genomic diversity in the Mexican population. In addition, the high prevalence of Alu insertions reveals their potential application as biomarkers with prognostic value for the associated diseases; e.g., as part of the standard protocols for clinical diagnosis.

Genetic dissection of five ethnic groups from Punjab, North-West India-A study based on Autosomal Markers

The present study assessed the applicability of Alu insertion elements and Single Nucleotide Polymorphisms (SNPs) in forensic identification and estimated the extent of genetic variation in five major ethnic groups of Punjab, North-West India. A total of 1012 unrelated samples belonging to Banias, Brahmins, Jat Sikhs, Khatris and Scheduled Castes were genotyped for four Alu elements (ACE, APO, PLAT, D1) and six Single Nucleotide Polymorphisms [ESR (PvuII), LPL (PvuII), HTR2A (MspI), DRD2 Taq1A, Taq1B, Taq1D]. Allele frequencies observed heterozygosity and forensic efficacy parameters were determined. The data on the genetic affinity of the studied populations among themselves and with other populations of India was also analysed using a Neighbor-Joining tree and multidimensional scaling plot respectively. All the 10 loci were polymorphic and their average observed heterozygosity ranged from 0.3872 (Banias) to 0.4311 (Scheduled Castes). Allele frequency variation at the 9 out of 10 loci led to statistically significant pairwise differences among the five study population groups. The result from AMOVA, Structure analysis, and Phylogenetic tree suggests that these populations are homogenous. In the multidimensional scaling plot, the present study populations formed a compact cluster clearly separated from other populations, suggesting a unique genetic identity of the Punjab populations as a whole. All these observations suggest that either a recent common origin of these populations or extensive gene flow across the populations that dissolve the original genetic differences. The data generated in this study will be useful for forensic genetics, molecular anthropological and demographic studies.

Genetic variation of MHC Class I polymorphic Alu insertions (POALINs) in three sub-populations of the East Midlands, UK

BACKGROUND

Alu elements are highly researched due to their useful nature as markers in the study of human population genetics. Recently discovered Major Histocompatibility Complex (MHC) polymorphic Alu insertions (POALINs) have not been examined extensively for genetic variation and their HLA associations.

AIMS

The aim of this study is to assess the genetic variation between three populations using five recently discovered POALINs.

METHODS AND SUBJECTS

The study examined 190 healthy, unrelated subjects from three different populations in the East Midlands (UK) for the presence or absence of five Alu elements (AluHG, AluMICB, AluHJ, AluTF and AluHF) via the polymerase chain reaction followed by gel electrophoresis. Data were analysed for genetic variation and phylogenetic analyses.

RESULTS

All Alus were polymorphic in study populations. Appreciable allele frequency variation was observed at a number of loci. The British population was significantly different from both the Punjabi Jat Sikh and Gujarati Patel populations, although showing a closer genetic relationship to the Punjabi Jat Sikh population than the Gujarati Patel population (Nei's D_A = 0.0031 and 0.0064, respectively).

CONCLUSIONS

MHC POALINs are useful markers in the investigation of genetic variation and the assessment of population relationships, and may have some bearing on disease associations due to their linkage disequilibrium with HLA loci; this warrants further studies.

Population and clinical genetics of human transposable elements in the (post) genomic era

Recent technological developments-in genomics, bioinformatics and high-throughput experimental techniques-are providing opportunities to study ongoing human transposable element (TE) activity at an unprecedented level of detail. It is now possible to characterize genome-wide collections of TE insertion sites for multiple human individuals, within and between populations, and for a variety of tissue types. Comparison of TE insertion site profiles between individuals captures the germline activity of TEs and reveals insertion site variants that segregate as polymorphisms among human populations, whereas comparison among tissue types ascertains somatic TE activity that generates cellular heterogeneity. In this review, we provide an overview of these new technologies and explore their implications for population and clinical genetic studies of human TEs. We cover both recent published results on human TE insertion activity as well as the prospects for future TE studies related to human evolution and health.

Genome-Wide Association Mapping Reveals Multiple QTLs Governing Tolerance Response for Seedling Stage Chilling Stress in Indica Rice

Rice crop is sensitive to cold stress at seedling stage. A panel of population representing 304 shortlisted germplasm lines was studied for seedling stage chilling tolerance in indica rice. Six phenotypic classes were exposed to six low temperature stress regimes under control phenotyping facility to investigate response pattern. A panel of 66 genotypes representing all phenotypic classes was used for ensuring genetic diversity, population structure and association mapping for the trait using 58 simple sequence repeat (SSR) and 2 direct trait linked markers. A moderate level of genetic diversity was detected in the panel population for the trait. Deviation of Hardy-Weinberg's expectation was detected in the studied population using Wright's F statistic. The panel showed 30% variation among population and 70% among individuals. The entire population was categorized into three sub-populations through STRUCTURE analysis. This revealed tolerance for the trait had a common primary ancestor for each sub-population with few admix individuals. The panel population showed the presence of many QTLs for cold stress tolerance in the individuals representing like genome-wide expression of the trait. Nineteen SSR markers were significantly associated at chilling stress of 8°C to 4°C for 7-21 days duration. Thus, the primers linked to the seedling stage cold tolerance QTLs namely qCTS9, qCTS-2, qCTS6.1, qSCT2, qSCT11, qSCT1a, qCTS-3.1, qCTS11.1, qCTS12.1, qCTS-1b, and CTB2 need to be pyramided for development of strongly chilling tolerant variety.

Human Nutrition and Genetic Variation

Human genetic variation is a determinant of nutrient efficacy and of tolerances and intolerances and has the potential to influence nutrient intake values (NIVs). Knowledge derived from the comprehensive identification of human genetic variation offers the potential to predict the physiological and pathological consequences of individual genetic differences and prevent and/or manage adverse outcomes through diet. Nutrients and genomes interact reciprocally; genomes confer differences in nutrient utilization, whereas nutrients effectively modify genome expression, stability, and viability. Understanding the interactions that occur among human genes, including all genetic variants thereof, and environmental exposures is enabling the development of genotype-specific nutritional regimens that prevent disease and promote wellness for individuals and populations throughout the life cycle. Genomic technologies may provide new criteria for establishing NIVs. The impact of a gene variant on NIVs will be dependent on its penetrance and prevalence within a population. Recent experiences indicate that few gene variants are anticipated to be sufficiently penetrant to affect average requirement (AR) values to a greater degree than environmental factors. If highly penetrant gene variants are identified that affect nutrient requirements, the prevalence of the variant in that country or region will determine the feasibility and necessity of deriving more than one AR or upper limit (UL) for affected genetic subgroups.

QTL Detection and Elite Alleles Mining for Stigma Traits in Oryza sativa by Association Mapping

Stigma traits are very important for hybrid seed production in Oryza sativa, which is a self-pollinated crop; however, the genetic mechanism controlling the traits is poorly understood. In this study, we investigated the phenotypic data of 227 accessions across 2 years and assessed their genotypic variation with 249 simple sequence repeat (SSR) markers. By combining phenotypic and genotypic data, a genome-wide association (GWA) map was generated. Large phenotypic variations in stigma length (STL), stigma brush-shaped part length (SBPL) and stigma non-brush-shaped part length (SNBPL) were found. Significant positive correlations were identified among stigma traits. In total, 2072 alleles were detected among 227 accessions, with an average of 8.3 alleles per SSR locus. GWA mapping detected 6 quantitative trait loci (QTLs) for the STL, 2 QTLs for the SBPL and 7 QTLs for the SNBPL. Eleven, 5, and 12 elite alleles were found for the STL, SBPL, and SNBPL, respectively. Optimal cross designs were predicted for improving the target traits. The detected genetic variation in stigma traits and QTLs provides helpful information for cloning candidate STL genes and breeding rice cultivars with longer STLs in the future.

Polymorphic Alu Insertion/Deletion in Different Caste and Tribal Populations from South India

Seven human-specific Alu markers were studied in 574 unrelated individuals from 10 endogamous groups and 2 hill tribes of Tamil Nadu and Kerala states. DNA was isolated, amplified by PCR-SSP, and subjected to agarose gel electrophoresis, and genotypes were assigned for various Alu loci. Average heterozygosity among caste populations was in the range of 0.292-0.468. Among tribes, the average heterozygosity was higher for Paliyan (0.3759) than for Kani (0.2915). Frequency differences were prominent in all loci studied except Alu CD4. For Alu CD4, the frequency was 0.0363 in Yadavas, a traditional pastoral and herd maintaining population, and 0.2439 in Narikuravars, a nomadic gypsy population. The overall genetic difference (Gst) of 12 populations (castes and tribes) studied was 3.6%, which corresponds to the Gst values of 3.6% recorded earlier for Western Asian populations. Thus, our study confirms the genetic similarities between West Asian populations and South Indian castes and tribes and supported the large scale coastal migrations from Africa into India through West Asia. However, the average genetic difference (Gst) of Kani and Paliyan tribes with other South Indian tribes studied earlier was 8.3%. The average Gst of combined South and North Indian Tribes (CSNIT) was 9.5%. Neighbor joining tree constructed showed close proximity of Kani and Paliyan tribal groups to the other two South Indian tribes, Toda and Irula of Nilgiri hills studied earlier. Further, the analysis revealed the affinities among populations and confirmed the presence of North and South India specific lineages. Our findings have documented the highly diverse (micro differentiated) nature of South Indian tribes, predominantly due to isolation, than the endogamous population groups of South India. Thus, our study firmly established the genetic relationship of South Indian castes and tribes and supported the proposed large scale ancestral migrations from Africa, particularly into South India through West Asian corridor.

Mitochondrial Genome and Nuclear Markers Provide New Insight into the Evolutionary History of Macaques

The evolutionary history of macaques, genus Macaca, has been under debate due to the short times of divergence. In this study, maternal, paternal, and biparental genetic systems were applied to infer phylogenetic relationships among macaques and to trace ancient hybridization events in their evolutionary history. Using a PCR display method, 17 newly phylogenetically informative Alu insertions were identified from M. assamensis. We combined presence/absence analysis of 84 Alu elements with mitochondrial genomes as well as nuclear sequences (five autosomal genes, two Y chromosomal genes, and one X chromosomal fragment) to reconstruct a robust macaque phylogeny. Topologies generated from different inherited markers were similar supporting six well defined species groups and a close relationship of M. assamensis and M. thibetana, but differed in the placing of M. arctoides. Both Alu elements and nuclear genes supported that M. arctoides was close to the sinica group, whereas the mitochondrial data clustered it into the fascicularis/mulatta lineage. Our results reveal that a sex-biased hybridization most likely occurred in the evolutionary history of M. arctoides, and suggest an introgressive pattern of male-mediated gene flow from the ancestors of M. arctoides to the M. mulatta population followed by nuclear swamping. According to the estimation of divergence dates, the hybridization occurred around 0.88~1.77 mya (nuclear data) or 1.38~2.56 mya (mitochondrial data). In general, our study indicates that a combination of various molecular markers could help explain complicated evolutionary relationships. Our results have provided new insights into the evolutionary history of macaques and emphasize that hybridization might play an important role in macaque evolution.

Alu SINE analyses of 3,000-year-old human skeletal remains: a pilot study

Background

As Short Interspersed Elements (SINEs), human-specific Alu elements can be used for population genetic studies. Very recent inserts are polymorphic within and between human populations. In a sample of 30 elements originating from three different Alu subfamilies, we investigated whether they are preserved in prehistorical skeletal human remains from the Bronze Age Lichtenstein cave in Lower Saxony, Germany. In the present study, we examined a prehistoric triad of father, mother and daughter.

Results

For 26 of the 30 Alu loci investigated, definite results were obtained. We were able to demonstrate that presence/absence analyses of Alu elements can be conducted on individuals who lived 3,000 years ago. The preservation of the ancient DNA (aDNA) is good enough in two out of three ancient individuals to routinely allow the amplification of 500 bp fragments. The third individual revealed less well-preserved DNA, which results in allelic dropout or complete amplification failures. We here present an alternative molecular approach to deal with these degradation phenomena by using internal Alu subfamily specific primers producing short fragments of approximately 150 bp.

Conclusions

Our data clearly show the possibility of presence/absence analyses of Alu elements in individuals from the Lichtenstein cave. Thus, we demonstrate that our method is reliably applicable for aDNA samples with good or moderate DNA preservation. This method will be very useful for further investigations with more Alu loci and larger datasets. Human population genetic studies and other large-scale investigations would provide insight into Alu SINE-based microevolutionary processes in humans during the last few thousand years and help us comprehend the evolutionary dynamics of our genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0063-y) contains supplementary material, which is available to authorized users.

Human Dispersal Out of Africa: A Lasting Debate

Unraveling the first migrations of anatomically modern humans out of Africa has invoked great interest among researchers from a wide range of disciplines. Available fossil, archeological, and climatic data offer many hypotheses, and as such genetics, with the advent of genome-wide genotyping and sequencing techniques and an increase in the availability of ancient samples, offers another important tool for testing theories relating to our own history. In this review, we report the ongoing debates regarding how and when our ancestors left Africa, how many waves of dispersal there were and what geographical routes were taken. We explore the validity of each, using current genetic literature coupled with some of the key archeological findings.

Transposable element polymorphisms recapitulate human evolution

Background

The human genome contains several active families of transposable elements (TE): Alu, L1 and SVA. Germline transposition of these elements can lead to polymorphic TE (polyTE) loci that differ between individuals with respect to the presence/absence of TE insertions. Limited sets of such polyTE loci have proven to be useful as markers of ancestry in human population genetic studies, but until this time it has not been possible to analyze the full genomic complement of TE polymorphisms in this way.

Results

For the first time here, we have performed a human population genetic analysis based on a genome-wide polyTE data set consisting of 16,192 loci genotyped in 2,504 individuals across 26 human populations. PolyTEs are found at very low frequencies, > 93 % of loci show < 5 % allele frequency, consistent with the deleteriousness of TE insertions. Nevertheless, polyTEs do show substantial geographic differentiation, with numerous group-specific polymorphic insertions. African populations have the highest numbers of polyTEs and show the highest levels of polyTE genetic diversity; Alu is the most numerous and the most diverse polyTE family. PolyTE genotypes were used to compute allele sharing distances between individuals and to relate them within and between human populations. Populations and continental groups show high coherence based on individuals’ polyTE genotypes, and human evolutionary relationships revealed by these genotypes are consistent with those seen for SNP-based genetic distances. The patterns of genetic diversity encoded by TE polymorphisms recapitulate broad patterns of human evolution and migration over the last 60–100,000 years. The utility of polyTEs as ancestry informative markers is further underscored by their ability to accurately predict both ancestry and admixture at the continental level. A genome-wide list of polyTE loci, along with their population group-specific allele frequencies and F_ST values, is provided as a resource for investigators who wish to develop panels of TE-based ancestry markers.

Conclusions

The genetic diversity represented by TE polymorphisms reflects known patterns of human evolution, and ensembles of polyTE loci are suitable for both ancestry and admixture analyses. The patterns of polyTE allelic diversity suggest the possibility that there may be a connection between TE-based genetic divergence and population-specific phenotypic differences.

Graphical Abstract

Unity in diversity: an overview of the genomic anthropology of India

CONTEXT

India is considered a treasure for geneticists and evolutionary biologists due to its vast human diversity, consisting of more than 4500 anthropologically well-defined populations (castes, tribes and religious groups). Each population differs in terms of endogamy, language, culture, physical features, geographic and climatic position and genetic architecture. These factors contributed to India-specific genetic variations which may be responsible for various common diseases in India and its migratory populations. As a result, interpretations of the origins and affinities of Indian populations as well as health and disease conditions require complex and sophisticated genetic analysis. Evidence of ancient human dispersals and settlements is preserved in the genome of Indian inhabitants and this has been extensively analysed in conventional and genomic analyses.

OBJECTIVE AND METHODS

Using genomic analyses of STRs and Alu on a set of populations, this study estimates the level and extent of genetic variation and its implications.

RESULTS

The results show that Indian populations have a higher level of unique genetic diversity which is structured by many social processes and geographical attributes of the country.

CONCLUSION

This overview highlights the need to study the anthropological structure and evolutionary history of Indian populations while designing genomic and epigenomic investigations.

A genomic portrait of haplotype diversity and signatures of selection in indigenous southern African populations

We report a study of genome-wide, dense SNP (∼900K) and copy number polymorphism data of indigenous southern Africans. We demonstrate the genetic contribution to southern and eastern African populations, which involved admixture between indigenous San, Niger-Congo-speaking and populations of Eurasian ancestry. This finding illustrates the need to account for stratification in genome-wide association studies, and that admixture mapping would likely be a successful approach in these populations. We developed a strategy to detect the signature of selection prior to and following putative admixture events. Several genomic regions show an unusual excess of Niger-Kordofanian, and unusual deficiency of both San and Eurasian ancestry, which were considered the footprints of selection after population admixture. Several SNPs with strong allele frequency differences were observed predominantly between the admixed indigenous southern African populations, and their ancestral Eurasian populations. Interestingly, many candidate genes, which were identified within the genomic regions showing signals for selection, were associated with southern African-specific high-risk, mostly communicable diseases, such as malaria, influenza, tuberculosis, and human immunodeficiency virus/AIDs. This observation suggests a potentially important role that these genes might have played in adapting to the environment. Additionally, our analyses of haplotype structure, linkage disequilibrium, recombination, copy number variation and genome-wide admixture highlight, and support the unique position of San relative to both African and non-African populations. This study contributes to a better understanding of population ancestry and selection in south-eastern African populations; and the data and results obtained will support research into the genetic contributions to infectious as well as non-communicable diseases in the region.

Genome-wide analysis of human populations is useful in shedding light on the evolutionary history of the human genome, with a wide range of applications from reconstructing past associations between different population histories to disease mapping. In this manuscript we report on the application of genome-wide data to southern African populations and the identification of genome-wide signatures of selection pre- and post-admixture. Several signals of selection, before and after admixture, were identified, some of which involved loci associated with human diseases, including malaria, influenza, tuberculosis and HIV/AIDS. These results may reflect adaptations of southern African populations to infectious diseases. Consistent with previous studies, this study highlights the significance of the San in the genetics of human populations, as they are distinct from the other populations in many respects i.e. haplotype structure, locations of recombination hotspots, copy number and population structure. Furthermore, our study demonstrates the admixture of the San, Bantu-speaking populations and populations of Eurasian ancestry in some of the southern and eastern African populations. It illustrates the value in correcting for this stratification in future genome-wide association studies, and suggests that a future admixture mapping in these populations would likely be warranted and successful.

Evaluating intra- and inter-individual variation in the human placental transcriptome

BACKGROUND

Gene expression variation is a phenotypic trait of particular interest as it represents the initial link between genotype and other phenotypes. Analyzing how such variation apportions among and within groups allows for the evaluation of how genetic and environmental factors influence such traits. It also provides opportunities to identify genes and pathways that may have been influenced by non-neutral processes. Here we use a population genetics framework and next generation sequencing to evaluate how gene expression variation is apportioned among four human groups in a natural biological tissue, the placenta.

RESULTS

We estimate that on average, 33.2%, 58.9%, and 7.8% of the placental transcriptome is explained by variation within individuals, among individuals, and among human groups, respectively. Additionally, when technical and biological traits are included in models of gene expression they each account for roughly 2% of total gene expression variation. Notably, the variation that is significantly different among groups is enriched in biological pathways associated with immune response, cell signaling, and metabolism. Many biological traits demonstrate correlated changes in expression in numerous pathways of potential interest to clinicians and evolutionary biologists. Finally, we estimate that the majority of the human placental transcriptome exhibits expression profiles consistent with neutrality; the remainder are consistent with stabilizing selection, directional selection, or diversifying selection.

CONCLUSIONS

We apportion placental gene expression variation into individual, population, and biological trait factors and identify how each influence the transcriptome. Additionally, we advance methods to associate expression profiles with different forms of selection.

Identification of species-specific nuclear insertions of mitochondrial DNA (numts) in gorillas and their potential as population genetic markers

The first hyper-variable region (HV1) of the mitochondrial control region (MCR) has been widely used as a molecular tool in population genetics, but inadvertent amplification of nuclear translocated copies of mitochondrial DNA (numts) in gorillas has compromised the use of mitochondrial DNA in population genetic studies. At least three putative classes (I, II, III) of gorilla-specific HV1 MCR numts have been uncovered over the past decade. However, the number, size and location of numt loci in gorillas and other apes are completely unknown. Furthermore, little work to date has assessed the utility of numts as candidate population genetic markers. In the present study, we screened Bacterial Artificial Chromosome (BAC) genomic libraries in the chimpanzee and gorilla to compare patterns of mitochondrial-wide insertion in both taxa. We conducted an intensive BLAST search for numts in the gorilla genome and compared the prevalence of numt loci originating from the MCR with other great ape taxa. Additional gorilla-specific MCR numts were retrieved either through BAC library screens or using an anchored-PCR (A-PCR) amplification using genomic DNA from five unrelated gorillas. Locus-specific primers were designed to identify numt insertional polymorphisms and evaluate their potential as population genetic markers. Mitochondrial-wide surveys of chimpanzee and gorilla BACs showed that the number of numts does not differ between these two taxa. However, MCR numts are more abundant in chimpanzees than in other great apes. We identified and mapped 67 putative gorilla-specific numts, including two that contain the entire HV1 domain, cluster with sequences from two numt classes (I, IIb) and will likely co-amplify with mitochondrial sequences using most published HV1 primers. However, phylogenetic analysis coupled with post-hoc analysis of mitochondrial variation can successfully differentiate nuclear sequences. Insertional polymorphisms were evident in three out of five numts examined, indicating their potential utility as molecular markers. Taken together, these findings demonstrate the potentially powerful insight that numts could make in uncovering population history in gorillas and other mammals.

The population history of endogenous retroviruses in mule deer (Odocoileus hemionus)

Mobile elements are powerful agents of genomic evolution and can be exceptionally informative markers for investigating species and population-level evolutionary history. While several studies have utilized retrotransposon-based insertional polymorphisms to resolve phylogenies, few population studies exist outside of humans. Endogenous retroviruses are LTR-retrotransposons derived from retroviruses that have become stably integrated in the host genome during past infections and transmitted vertically to subsequent generations. They offer valuable insight into host-virus co-evolution and a unique perspective on host evolutionary history because they integrate into the genome at a discrete point in time. We examined the evolutionary history of a cervid endogenous gammaretrovirus (CrERVγ) in mule deer (Odocoileus hemionus). We sequenced 14 CrERV proviruses (CrERV-in1 to -in14), and examined the prevalence and distribution of 13 proviruses in 262 deer among 15 populations from Montana, Wyoming, and Utah. CrERV absence in white-tailed deer (O. virginianus), identical 5' and 3' long terminal repeat (LTR) sequences, insertional polymorphism, and CrERV divergence time estimates indicated that most endogenization events occurred within the last 200000 years. Population structure inferred from CrERVs (F ST = 0.008) and microsatellites (θ = 0.01) was low, but significant, with Utah, northwestern Montana, and a Helena herd being particularly differentiated. Clustering analyses indicated regional structuring, and non-contiguous clustering could often be explained by known translocations. Cluster ensemble results indicated spatial localization of viruses, specifically in deer from northeastern and western Montana. This study demonstrates the utility of endogenous retroviruses to elucidate and provide novel insight into both ERV evolutionary history and the history of contemporary host populations.

Rates and patterns of great ape retrotransposition

We analyzed 83 fully sequenced great ape genomes for mobile element insertions, predicting a total of 49,452 fixed and polymorphic Alu and long interspersed element 1 (L1) insertions not present in the human reference assembly and assigning each retrotransposition event to a different time point during great ape evolution. We used these homoplasy-free markers to construct a mobile element insertions-based phylogeny of humans and great apes and demonstrate their differential power to discern ape subspecies and populations. Within this context, we find a good correlation between L1 diversity and single-nucleotide polymorphism heterozygosity (r(2) = 0.65) in contrast to Alu repeats, which show little correlation (r(2) = 0.07). We estimate that the "rate" of Alu retrotransposition has differed by a factor of 15-fold in these lineages. Humans, chimpanzees, and bonobos show the highest rates of Alu accumulation--the latter two since divergence 1.5 Mya. The L1 insertion rate, in contrast, has remained relatively constant, with rates differing by less than a factor of three. We conclude that Alu retrotransposition has been the most variable form of genetic variation during recent human-great ape evolution, with increases and decreases occurring over very short periods of evolutionary time.

Genomic distribution of SINEs in Entamoeba histolytica strains: implication for genotyping

BACKGROUND

The major clinical manifestations of Entamoeba histolytica infection include amebic colitis and liver abscess. However the majority of infections remain asymptomatic. Earlier reports have shown that some E. histolytica isolates are more virulent than others, suggesting that virulence may be linked to genotype. Here we have looked at the genomic distribution of the retrotransposable short interspersed nuclear elements EhSINE1 and EhSINE2. Due to their mobile nature, some EhSINE copies may occupy different genomic locations among isolates of E. histolytica possibly affecting adjacent gene expression; this variability in location can be exploited to differentiate strains.

RESULTS

We have looked for EhSINE1- and EhSINE2-occupied loci in the genome sequence of Entamoeba histolytica HM-1:IMSS and searched for homologous loci in other strains to determine the insertion status of these elements. A total of 393 EhSINE1 and 119 EhSINE2 loci were analyzed in the available sequenced strains (Rahman, DS4-868, HM1:CA, KU48, KU50, KU27 and MS96-3382. Seventeen loci (13 EhSINE1 and 4 EhSINE2) were identified where a EhSINE1/EhSINE2 sequence was missing from the corresponding locus of other strains. Most of these loci were unoccupied in more than one strain. Some of the loci were analyzed experimentally for SINE occupancy using DNA from strain Rahman. These data helped to correctly assemble the nucleotide sequence at three loci in Rahman. SINE occupancy was also checked at these three loci in 7 other axenically cultivated E. histolytica strains and 16 clinical isolates. Each locus gave a single, specific amplicon with the primer sets used, making this a suitable method for strain typing. Based on presence/absence of SINE and amplification with locus-specific primers, the 23 strains could be divided into eleven genotypes. The results obtained by our method correlated with the data from other typing methods. We also report a bioinformatic analysis of EhSINE2 copies.

CONCLUSIONS

Our results reveal several loci with extensive polymorphism of SINE occupancy among different strains of E. histolytica and prove the principle that the genomic distribution of SINEs is a valid method for typing of E. histolytica strains.

Mobile element scanning (ME-Scan) identifies thousands of novel Alu insertions in diverse human populations

Alu retrotransposons are the most numerous and active mobile elements in humans, causing genetic disease and creating genomic diversity. Mobile element scanning (ME-Scan) enables comprehensive and affordable identification of mobile element insertions (MEI) using targeted high-throughput sequencing of multiplexed MEI junction libraries. In a single experiment, ME-Scan identifies nearly all AluYb8 and AluYb9 elements, with high sensitivity for both rare and common insertions, in 169 individuals of diverse ancestry. ME-Scan detects heterozygous insertions in single individuals with 91% sensitivity. Insertion presence or absence states determined by ME-Scan are 95% concordant with those determined by locus-specific PCR assays. By sampling diverse populations from Africa, South Asia, and Europe, we are able to identify 5799 Alu insertions, including 2524 novel ones, some of which occur in exons. Sub-Saharan populations and a Pygmy group in particular carry numerous intermediate-frequency Alu insertions that are absent in non-African groups. There is a significant dearth of exon-interrupting insertions among common Alu polymorphisms, but the density of singleton Alu insertions is constant across exonic and nonexonic regions. In one case, a validated novel singleton Alu interrupts a protein-coding exon of FAM187B. This implies that exonic Alu insertions are generally deleterious and thus eliminated by natural selection, but not so quickly that they cannot be observed as extremely rare variants.

Human endogenous retrovirus-H insertion screening

Endogenous retroviruses (ERVs) and ERV-like sequences comprise 8% of the human genome. We aimed to analyze genome integration polymorphisms of human endogenous retrovirus (HERV)-H by the inter-retrotransposon amplified polymorphism (IRAP) technique using the sequences of LTR7A (450 bp), LTR7B (445 bp) and LTR7C (471 bp). Blood samples from 20 individuals (10 females and 10 males) of diverse ethnic origins were used for the determination of integration variations at the genomic level. Isolated genomic DNA was screened using 3 pairs of primers corresponding to LTR regions of the HERV-H gene. We observed insertion polymorphism patterns between 0-87% in all subjects. The findings obtained contribute to our understanding of the effects of HERV-H on variations within the human genome.

Mobile element biology: new possibilities with high-throughput sequencing

Mobile elements comprise more than half of the human genome, but until recently their large-scale detection was time consuming and challenging. With the development of new high-throughput sequencing (HTS) technologies, the complete spectrum of mobile element variation in humans can now be identified and analyzed. Thousands of new mobile element insertions (MEIs) have been discovered, yielding new insights into mobile element biology, evolution, and genomic variation. Here, we review several high-throughput methods, with an emphasis on techniques that specifically target MEIs in humans. We highlight recent applications of these methods in evolutionary studies and in the analysis of somatic alterations in human normal and tumor tissues.

Transposable Elements: No More 'Junk DNA'

Since the advent of whole-genome sequencing, transposable elements (TEs), just thought to be 'junk' DNA, have been noticed because of their numerous copies in various eukaryotic genomes. Many studies about TEs have been conducted to discover their functions in their host genomes. Based on the results of those studies, it has been generally accepted that they have a function to cause genomic and genetic variations. However, their infinite functions are not fully elucidated. Through various mechanisms, including de novo TE insertions, TE insertion-mediated deletions, and recombination events, they manipulate their host genomes. In this review, we focus on Alu, L1, human endogenous retrovirus, and short interspersed element/variable number of tandem repeats/Alu (SVA) elements and discuss how they have affected primate genomes, especially the human and chimpanzee genomes, since their divergence.

Introducing the Forensic Research/Reference on Genetics knowledge base, FROG-kb

BACKGROUND

Online tools and databases based on multi-allelic short tandem repeat polymorphisms (STRPs) are actively used in forensic teaching, research, and investigations. The Fst value of each CODIS marker tends to be low across the populations of the world and most populations typically have all the common STRP alleles present diminishing the ability of these systems to discriminate ethnicity. Recently, considerable research is being conducted on single nucleotide polymorphisms (SNPs) to be considered for human identification and description. However, online tools and databases that can be used for forensic research and investigation are limited.

METHODS

The back end DBMS (Database Management System) for FROG-kb is Oracle version 10. The front end is implemented with specific code using technologies such as Java, Java Servlet, JSP, JQuery, and GoogleCharts.

RESULTS

We present an open access web application, FROG-kb (Forensic Research/Reference on Genetics-knowledge base, http://frog.med.yale.edu), that is useful for teaching and research relevant to forensics and can serve as a tool facilitating forensic practice. The underlying data for FROG-kb are provided by the already extensively used and referenced ALlele FREquency Database, ALFRED (http://alfred.med.yale.edu). In addition to displaying data in an organized manner, computational tools that use the underlying allele frequencies with user-provided data are implemented in FROG-kb. These tools are organized by the different published SNP/marker panels available. This web tool currently has implemented general functions possible for two types of SNP panels, individual identification and ancestry inference, and a prediction function specific to a phenotype informative panel for eye color.

CONCLUSION

The current online version of FROG-kb already provides new and useful functionality. We expect FROG-kb to grow and expand in capabilities and welcome input from the forensic community in identifying datasets and functionalities that will be most helpful and useful. Thus, the structure and functionality of FROG-kb will be revised in an ongoing process of improvement. This paper describes the state as of early June 2012.

An Alu-based phylogeny of gibbons (hylobatidae)

Gibbons (Hylobatidae) are small, arboreal apes indigenous to Southeast Asia that diverged from other apes ∼15-18 Ma. Extant lineages radiated rapidly 6-10 Ma and are organized into four genera (Hylobates, Hoolock, Symphalangus, and Nomascus) consisting of 12-19 species. The use of short interspersed elements (SINEs) as phylogenetic markers has seen recent popularity due to several desirable characteristics: the ancestral state of a locus is known to be the absence of an element, rare potentially homoplasious events are relatively easy to resolve, and samples can be quickly and inexpensively genotyped. During radiation of primates, one particular family of SINEs, the Alu family, has proliferated in primate genomes. Nomascus leucogenys (northern white-cheeked gibbon) sequences were analyzed for repetitive content with RepeatMasker using a custom library. The sequences containing Alu elements identified as members of a gibbon-specific subfamily were then compared with orthologous positions in other primate genomes. A primate phylogenetic panel consisting of 18 primate species, including 13 gibbon species representing all four extant genera, was assayed for all loci, and a total of 125 gibbon-specific Alu insertions were identified. The resulting amplification patterns were used to generate a phylogenetic tree. We demonstrate significant support for Symphalangus as the most basal lineage within the family. Our findings also place Nomascus as a derived lineage, sister to Hoolock, with the Nomascus-Hoolock clade sister to Hylobates. Further, our analysis groups N. leucogenys and Nomascus siki as sister taxa to the exclusion of the other Nomascus species assayed. This study represents the first use of SINEs to determine the genus level phylogenetic relationships within the family Hylobatidae. These relationships have been resolved with robust support at most internal nodes, demonstrating the utility of SINE-based phylogenetic analysis. We postulate that hybridization and rapid radiation may have contributed to the complex and contradictory findings of the previous studies. Our findings will aid in the conservation of these threatened primates and inform future studies of the biogeographical history and distribution of modern gibbon species.

Gender-specific association of TSNAX/DISC1 locus for schizophrenia and bipolar affective disorder in South Indian population

Genetic association studies have implicated the TSNAX/DISC1 (disrupted in schizophrenia 1) in schizophrenia (SCZ), bipolar affective disorder (BPAD) and major depression. This study was performed to assess the possible involvement of TSNAX/DISC1 locus in the aetiology of BPAD and SCZ in the Southern Indian population. We genotyped seven single nucleotide polymorphism (SNPs) from TSNAX/DISC1 region in 1252 individuals (419 BPAD patients, 408 SCZ patients and 425 controls). Binary logistic regression revealed a nominal association for rs821616 in DISC1 for BPAD and also combined cases of BPAD or SCZ, but after correcting for multiple testing, these results were non-significant. However, significant association was observed with BPAD, as well as combined cases of BPAD or SCZ, within the female subjects for the rs766288 after applying false discovery rate corrections at the 0.05 level. Two-locus analysis showed C-C (rs766288-rs2812393) as a risk combination in BPAD, and G-T (rs2812393-rs821616) as a protective combination in SCZ and combined cases of BPAD or SCZ. Female-specific associations were observed for rs766288-rs2812393, rs766288-rs821616 and rs8212393-rs821616 in two-locus analysis. Our results provide further evidence for sex-dependent effects of the TSNAX/DISC1 locus in the aetiology of SCZ and BPAD.

Genome-wide patterns of genetic distances reveal candidate Loci contributing to human population-specific traits

Modern humans originated in Africa before migrating across the world with founder effects and adaptations to new environments contributing to their present phenotypic diversity. Determining the genetic basis of differences between populations may provide clues about our evolutionary history and may have clinical implications. Herein, we develop a method to detect genes and biological processes in which populations most differ by calculating the genetic distance between modern populations and a hypothetical ancestral population. We apply our method to large-scale single nucleotide polymorphism (SNP) data from human populations of African, European and Asian origin. As expected, ancestral alleles were more conserved in the African populations and we found evidence of high divergence in genes previously suggested as targets of selection related to skin pigmentation, immune response, senses and dietary adaptations. Our genome-wide scan also reveals novel candidates for contributing to population-specific traits. These include genes related to neuronal development and behavior that may have been influenced by cultural processes. Moreover, in the African populations, we found a high divergence in genes related to UV protection and to the male reproductive system. Taken together, these results confirm and expand previous findings, providing new clues about the evolution and genetics of human phenotypic diversity.

Shared and unique components of human population structure and genome-wide signals of positive selection in South Asia

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.