Skin pigmentation,biogeographical ancestry and admixture mapping

Ancestry informative markers (AIMs) are genetic loci showing alleles with large frequency differences between populations. AIMs can be used to estimate biogeographical ancestry at the level of the population, subgroup (e.g. cases and controls) and individual. Ancestry estimates at both the subgroup and individual level can be directly instructive regarding the genetics of the phenotypes that differ qualitatively or in frequency between populations. These estimates can provide a compelling foundation for the use of admixture mapping (AM) methods to identify the genes underlying these traits. We present details of a panel of 34 AIMs and demonstrate how such studies can proceed, by using skin pigmentation as a model phenotype. We have genotyped these markers in two population samples with primarily African ancestry, viz. African Americans from Washington D.C. and an African Caribbean sample from Britain, and in a sample of European Americans from Pennsylvania. In the two African population samples, we observed significant correlations between estimates of individual ancestry and skin pigmentation as measured by reflectometry (R(2)=0.21, P<0.0001 for the African-American sample and R(2)=0.16, P<0.0001 for the British African-Caribbean sample). These correlations confirm the validity of the ancestry estimates and also indicate the high level of population structure related to admixture, a level that characterizes these populations and that is detectable by using other tests to identify genetic structure. We have also applied two methods of admixture mapping to test for the effects of three candidate genes (TYR, OCA2, MC1R) on pigmentation. We show that TYR and OCA2 have measurable effects on skin pigmentation differences between the west African and west European parental populations. This work indicates that it is possible to estimate the individual ancestry of a person based on DNA analysis with a reasonable number of well-defined genetic markers. The implications and applications of ancestry estimates in biomedical research are discussed.     

Maximum-likelihood estimation of admixture proportions from genetic data

For an admixed population, an important question is how much genetic contribution comes from each parental population. Several methods have been developed to estimate such admixture proportions, using data on genetic markers sampled from parental and admixed populations. In this study, I propose a likelihood method to estimate jointly the admixture proportions, the genetic drift that occurred to the admixed population and each parental population during the period between the hybridization and sampling events, and the genetic drift in each ancestral population within the interval between their split and hybridization. The results from extensive simulations using various combinations of relevant parameter values show that in general much more accurate and precise estimates of admixture proportions are obtained from the likelihood method than from previous methods. The likelihood method also yields reasonable estimates of genetic drift that occurred to each population, which translate into relative effective sizes (N(e)) or absolute average N(e)'s if the times when the relevant events (such as population split, admixture, and sampling) occurred are known. The proposed likelihood method also has features such as relatively low computational requirement compared with previous ones, flexibility for admixture models, and marker types. In particular, it allows for missing data from a contributing parental population. The method is applied to a human data set and a wolflike canids data set, and the results obtained are discussed in comparison with those from other estimators and from previous studies.     

Recent admixture in an Indian population of African ancestry

Identification and study of genetic variation in recently admixed populations not only provides insight into historical population events but also is a powerful approach for mapping disease loci. We studied a population (OG-W-IP) that is of African-Indian origin and has resided in the western part of India for 500 years; members of this population are believed to be descendants of the Bantu-speaking population of Africa. We have carried out this study by using a set of 18,534 autosomal markers common between Indian, CEPH-HGDP, and HapMap populations. Principal-components analysis clearly revealed that the African-Indian population derives its ancestry from Bantu-speaking west-African as well as Indo-European-speaking north and northwest Indian population(s). STRUCTURE and ADMIXTURE analyses show that, overall, the OG-W-IPs derive 58.7% of their genomic ancestry from their African past and have very little inter-individual ancestry variation (8.4%). The extent of linkage disequilibrium also reveals that the admixture event has been recent. Functional annotation of genes encompassing the ancestry-informative markers that are closer in allele frequency to the Indian ancestral population revealed significant enrichment of biological processes, such as ion-channel activity, and cadherins. We briefly examine the implications of determining the genetic diversity of this population, which could provide opportunities for studies involving admixture mapping.     

Informativeness of genetic markers for inference of ancestry

Inference of individual ancestry is useful in various applications, such as admixture mapping and structured-association mapping. Using information-theoretic principles, we introduce a general measure, the informativeness for assignment (I(n)), applicable to any number of potential source populations, for determining the amount of information that multiallelic markers provide about individual ancestry. In a worldwide human microsatellite data set, we identify markers of highest informativeness for inference of regional ancestry and for inference of population ancestry within regions; these markers, which are listed in online-only tables in our article, can be useful both in testing for and in controlling the influence of ancestry on case-control genetic association studies. Markers that are informative in one collection of source populations are generally informative in others. Informativeness of random dinucleotides, the most informative class of microsatellites, is five to eight times that of random single-nucleotide polymorphisms (SNPs), but 2%-12% of SNPs have higher informativeness than the median for dinucleotides. Our results can aid in decisions about the type, quantity, and specific choice of markers for use in studies of ancestry.     

Phylogenetic inference under varying proportions of indel-induced alignment gaps

Background  

The effect of alignment gaps on phylogenetic accuracy has been the subject of numerous studies. In this study, we investigated the relationship between the total number of gapped sites and phylogenetic accuracy, when the gaps were introduced (by means of computer simulation) to reflect indel (insertion/deletion) events during the evolution of DNA sequences. The resulting (true) alignments were subjected to commonly used gap treatment and phylogenetic inference methods.     

The population genetic effects of ancestry and admixture in a subdivided cattle breed

The genetic structure of the Dexter, a minority cattle breed with complex demographic history, was investigated using microsatellite markers and a range of statistical approaches designed to detect both admixture and genetic drift. Modern representatives of two putative ancestral populations, the Devon and Kerry, together with the different populations of the Dexter, which have experienced different demographic histories, were analysed. Breed units showed comparatively high levels of genetic variability ( H _E = 0.63–0.68); however, distinct genetic subgroups were detected within the Dexter, which could be attributed to known demographic events. Much lower diversity was identified in three small, isolated Dexter populations ( H _E = 0.52–0.55) and higher differentiation ( F _ST > 0.13) was found. For one of these populations, where strong selection has taken place, we also found evidence of a demographic bottleneck. Three methods for quantifying breed admixture were applied and substantial method-based variation in estimates for the genetic contribution of the two proposed ancestral populations for each subdivision of the Dexter was found. Results were consistent only in the case of a group consisting of selected Traditional Dexter animals, where the ancestor of the modern Kerry breed was also determined as the greater parental contributor to the Dexter. The inconsistency of estimation of admixture proportions between the methods highlights the potentially confounding role of genetic drift in shaping small population structure, and the consequences of accurately describing population histories from contemporary genetic data.     

Estimating African American admixture proportions by use of population-specific alleles.

We analyzed the European genetic contribution to 10 populations of African descent in the United States (Maywood, Illinois; Detroit; New York; Philadelphia; Pittsburgh; Baltimore; Charleston, South Carolina; New Orleans; and Houston) and in Jamaica, using nine autosomal DNA markers. These markers either are population-specific or show frequency differences >45% between the parental populations and are thus especially informative for admixture. European genetic ancestry ranged from 6.8% (Jamaica) to 22.5% (New Orleans). The unique utility of these markers is reflected in the low variance associated with these admixture estimates (SEM 1.3%-2.7%). We also estimated the male and female European contribution to African Americans, on the basis of informative mtDNA (haplogroups H and L) and Y Alu polymorphic markers. Results indicate a sex-biased gene flow from Europeans, the male contribution being substantially greater than the female contribution. mtDNA haplogroups analysis shows no evidence of a significant maternal Amerindian contribution to any of the 10 populations. We detected significant nonrandom association between two markers located 22 cM apart (FY-null and AT3), most likely due to admixture linkage disequilibrium created in the interbreeding of the two parental populations. The strength of this association and the substantial genetic distance between FY and AT3 emphasize the importance of admixed populations as a useful resource for mapping traits with different prevalence in two parental populations.     

Genetic variation in Arizona Mexican Americans: estimation and interpretation of admixture proportions

Mexican Americans are a numerous and fast growing ethnic population in the United States. Yet little is known about their genetic structure. Since they are a hybrid, it is of interest to identify their parental populations and to estimate the relative contributions of these groups. This information is relevant to historical, biomedical, and evolutionary concerns. New genetic typings on 730 Arizona Mexican Americans for the HLA-A, HLA-B, ABO, Rh, MNSs, Duffy, Kidd, and Kell loci are presented here and they are used to estimate ancestral contributions. We considered both a dihybrid model with Amerindians and Spaniards as proposed ancestors, and a trihybrid model with Amerindians, Spaniards, and Africans as proposed ancestors. A modified weighted least squares method that allows for linkage disequilibrium was used to estimate ancestral contributions for each model. The following admixture estimates were obtained: Amerindian, 0.29 +/- 0.04; Spaniard, 0.68 +/- 0.05; and African, 0.03 +/- 0.02. The interpretation of these results with respect to Amerindian and Spanish ancestry is straightforward. African ancestry is strongly supported by the presence of a marker of African descent, Fy, despite the fact that the standard error of the estimate is as large as the estimated admixture proportion. An evaluation of the sensitivity of these results to a number of variables is presented: 1) our choices of ancestral allele frequencies, 2) the possibility of selection at HLA and the blood groups, and 3) genetic drift in Mexican Americans.     

Estimating admixture proportions with microsatellites: comparison of methods based on simulated data

Several methods have been developed to estimate the parental contributions in the genetic pool of an admixed population. Some pair-comparisons have been performed on real data but, to date, no systematic comparison of a large number of methods has been attempted. In this study, we performed a simulated data-based comparison of six of the most cited methods in the literature of the last 20 years. Five of these methods use allele frequencies and differ in the statistical treatment of the data. The last one also considers the degree of molecular divergence by estimating the coalescence times. Comparisons are based on the frequency at which the method can be applied, the bias and the mean square error of the estimation, and the frequency at which the true value is within the confidence interval. Eventually, each method was applied to a real data set of variously introgressed honeybee populations. In optimal conditions (highly differentiated parental populations, recent hybridization event), all methods perform equally well. When conditions are not optimal, the methods perform differently, but no method is always better or worse than all others. Some guidelines are given for the choice of the method.     

mtDna and the islands of the North Atlantic: estimating the proportions of Norse and Gaelic ancestry

A total of 1,664 new mtDNA control-region sequences were analyzed in order to estimate Gaelic and Scandinavian matrilineal ancestry in the populations of Iceland, Orkney, the Western Isles, and the Isle of Skye and to investigate other aspects of their genetic history. A relative excess of private lineages in the Icelanders is indicative of isolation, whereas the scarcity of private lineages in Scottish island populations may be explained by recent gene flow and population decline. Differences in the frequencies of lineage clusters are observed between the Scandinavian and the Gaelic source mtDNA pools, and, on a continent-wide basis, such differences between populations seem to be associated with geography. A multidimensional scaling analysis of genetic distances, based on mtDNA lineage-cluster frequencies, groups the North Atlantic islanders with the Gaelic and the Scandinavian populations, whereas populations from the central, southern, and Baltic regions of Europe are arranged in clusters in broad agreement with their geographic locations. This pattern is highly significant, according to a Mantel correlation between genetic and geographic distances (r=.716). Admixture analyses indicate that the ancestral contributions of mtDNA lineages from Scandinavia to the populations of Iceland, Orkney, the Western Isles, and the Isle of Skye are 37.5%, 35.5%, 11.5%, and 12.5%, respectively.     

Patterns of insertion and deletion in contrasting chromatin domains

Transposable elements (TEs) play a fundamental role in the evolution of genomes. In Drosophila they are disproportionately represented in regions of low recombination, such as in heterochromatin. This pattern has been attributed to selection against repeated elements in regions of normal recombination, owing to either (1) the slightly deleterious position effects of TE insertions near or into genes, or (2) strong selection against chromosomal abnormalities arising from ectopic exchange between TE repeats. We have used defective non-long-terminal repeat (LTR) TEs that are "dead-on-arrival" (DOA) and unable to transpose in order to estimate spontaneous deletion rates in different constituents of chromatin. These elements have previously provided evidence for an extremely high rate of spontaneous deletion in Drosophila as compared with mammals, potentially explaining at least part of the differences in the genome sizes in these organisms. However, rates of deletion could be overestimated due to positive selection for a smaller likelihood of ectopic exchange. In this article, we show that rates of spontaneous deletion in DOA repeats are as high in heterochromatin and regions of euchromatin with low recombination as they are in regions of euchromatin with normal recombination. We have also examined the age distribution of five non-LTR families throughout the genome. We show that there is substantial variation in the historical pattern of transposition of these TEs. The overrepresentation of TEs in the heterochromatin is primarily due to their longer retention time in heterochromatin, as evidenced by the average time since insertion. Fragments inserted recently are much more evenly distributed in the genome. This contrast demonstrates that the accumulation of TEs in heterochromatin and in euchromatic regions of low recombination is not due to biased transposition but by greater probabilities of fixation in these regions relative to regions of normal recombination.     

Development of a panel of genome-wide ancestry informative markers to study admixture throughout the Americas

Most individuals throughout the Americas are admixed descendants of Native American, European, and African ancestors. Complex historical factors have resulted in varying proportions of ancestral contributions between individuals within and among ethnic groups. We developed a panel of 446 ancestry informative markers (AIMs) optimized to estimate ancestral proportions in individuals and populations throughout Latin America. We used genome-wide data from 953 individuals from diverse African, European, and Native American populations to select AIMs optimized for each of the three main continental populations that form the basis of modern Latin American populations. We selected markers on the basis of locus-specific branch length to be informative, well distributed throughout the genome, capable of being genotyped on widely available commercial platforms, and applicable throughout the Americas by minimizing within-continent heterogeneity. We then validated the panel in samples from four admixed populations by comparing ancestry estimates based on the AIMs panel to estimates based on genome-wide association study (GWAS) data. The panel provided balanced discriminatory power among the three ancestral populations and accurate estimates of individual ancestry proportions (R2 > 0.9 for ancestral components with significant between-subject variance). Finally, we genotyped samples from 18 populations from Latin America using the AIMs panel and estimated variability in ancestry within and between these populations. This panel and its reference genotype information will be useful resources to explore population history of admixture in Latin America and to correct for the potential effects of population stratification in admixed samples in the region.     

Estimation of admixture proportions: a likelihood-based approach using Markov chain Monte Carlo.

When populations are separated for long periods and then brought into contact for a brief episode in part of their range, this can result in genetic admixture. To analyze this type of event we considered a simple model under which two parental populations (P1 and P2) mix and create a hybrid population (H). After that event, the three populations evolve under pure drift without exchange during T generations. We developed a new method, which allows the simultaneous estimation of the time since the admixture event (scaled by the population size t(i) = T/N(i), where N(i) is the effective population size of population i) and the contribution of one of two parental populations (which we call p1). This method takes into account drift since the admixture event, variation caused by sampling, and uncertainty in the estimation of the ancestral allele frequencies. The method is tested on simulated data sets and then applied to a human data set. We find that (i) for single-locus data, point estimates are poor indicators of the real admixture proportions even when there are many alleles; (ii) biallelic loci provide little information about the admixture proportion and the time since admixture, even for very small amounts of drift, but can be powerful when many loci are used; (iii) the precision of the parameters' estimates increases with sample size n = 50 vs. n = 200 but this effect is larger for the t(i)'s than for p1; and (iv) the increase in precision provided by multiple loci is quite large, even when there is substantial drift (we found, for instance, that it is preferable to use five loci than one locus, even when drift is 100 times larger for the five loci). Our analysis of a previously studied human data set illustrates that the joint estimation of drift and p1 can provide additional insights into the data.     

Statistical alignment based on fragment insertion and deletion models

Motivation: The topic of this paper is the estimation of alignments and mutation rates based on stochastic sequence-evolution models that allow insertions and deletions of subsequences ('fragments') and not just single bases. The model we propose is a variant of a model introduced by Thorne et al., (J. Mol. Evol., 34, 3-16, 1992). The computational tractability of the model depends on certain restrictions in the insertion/deletion process; possible effects we discuss. Results: The process of fragment insertion and deletion in the sequence-evolution model induces a hidden Markov structure at the level of alignments and thus makes possible efficient statistical alignment algorithms. As an example we apply a sampling procedure to assess the variability in alignment and mutation parameter estimates for HVR1 sequences of human and orangutan, improving results of previous work. Simulation studies give evidence that estimation methods based on the proposed model also give satisfactory results when applied to data for which the restrictions in the insertion/deletion process do not hold. Availability: The source code of the software for sampling alignments and mutation rates for a pair of DNA sequences according to the fragment insertion and deletion model is freely available from http://www.math.uni-frankfurt.de/~stoch/software/mcmcsalut under the terms of the GNU public license (GPL, 2000).     

Inferring admixture proportions from molecular data: extension to any number of parental populations

The relative contribution of two parental populations to a hybrid group (the admixture proportions) can be estimated using not only the frequencies of different alleles, but also the degree of molecular divergence between them. In this paper, we extend this possibility to the case of any number of parental populations. The newly derived multiparental estimator is tested by Monte Carlo simulations and by generating artificial hybrid groups by pooling mtDNA samples from human populations. The general properties (including the variance) of the two-parental estimator seem to be retained by the multiparental estimator. When mixed human populations are considered and hypervariable single-locus data are analyzed (mtDNA control region), errors in the estimated contributions appear reasonably low only when highly differentiated parental populations are involved. Finally, the method applied to the hybrid Canary Island population points to a much lower female contribution from Spain than has previously been estimated.     

基因组水平的插入和缺失变异研究进展

插入和缺失(insertion and deletion)是DNA和蛋白质在进化过程中发生的序列长度上的改变,由于缺乏祖先序列的信息,不能肯定其到底是插入事件还是缺失事件,故统称之为增减(indel)。indel是分子水平进化变异的主要来源之一,近年来对这种进化事件的研究已经涵盖了其发生频率、大小、分布模式、序列进化模型及应用等各个方面。该文总结了基因组水平上插入和缺失的研究进展和发生机制;介绍了已经提出的插入和缺失进化模型,包括TKF91、TKF92、Long Indel模型和序列环境模型;讨论了插入和缺失作为分子标记在分子进化、基因分型和药物设计等方面的应用。     

Ancestry informative DIP loci for dissecting genetic structure and ancestry proportions of Qinghai Tibetan and Tibet Tibetan groups

Molecular Biology Reports - Tibetans living in the Qing-Tibet plateau show unique genetic features since they are exposed to the high altitude environment. Accordingly, it is necessary for us to...     

Male ancestry structure and interethnic admixture in African-descent communities from the Amazon as revealed by Y-chromosome Strs

Some genetic markers on both the Y chromosome and mtDNA are highly polymorphic and population‐specific in humans, representing useful tools for reconstructing the past history of populations with poor historical records. Such lack of information is usually true in the case of recent African‐descent populations of the New World founded by fugitive slaves throughout the slavery period in the Americas, particularly in Brazil, where those communities are known as quilombos. Aiming to recover male‐derived ethnic structure of nine quilombos from the Brazilian Amazon, a total of 300 individuals, belonging to Mazagão Velho (N = 24), Curiaú (N = 48), Mazagão (N = 36), Trombetas (N = 20), Itacoã (N = 22), Saracura (N = 46), Marajó (N = 58), Pitimandeua (N = 26), and Pontal (N = 20), were investigated for nine Y‐STRs (DYS393, DYS19, DYS390, DYS389 I, DYS389 II, DYS392, DYS391, DYS385 I/II). From the 169 distinct haplotypes obtained, 120 were singletons. The results suggest the West African coast as the main origin of slaves brought to Brazil (54% of male contribution); the European contribution was high (41%), while the Amerindian's was low (5%). Those results contrast with previous mtDNA data that showed high Amerindian female contribution (46.6%) in African‐descent populations. AMOVA suggests that the genetic differentiation among the quilombos is mainly influenced by admixture with European. However, when restricting AMOVA to African‐specific haplotypes, low differentiation was detected, suggesting great genetic homogeneity of the African founding populations and/or a later homogenization by intense slave trade inside Brazil. Am J Phys Anthropol, 2011. © 2010 Wiley‐Liss, Inc.