Mapping of disease-associated variants in admixed populations

Recent developments in high-throughput genotyping and whole-genome sequencing will enhance the identification of disease loci in admixed populations. We discuss how a more refined estimation of ancestry benefits both admixture mapping and association mapping, making disease loci identification in admixed populations more powerful.     

Case-control studies of association in structured or admixed populations

Case-control tests for association are an important tool for mapping complex-trait genes. But population structure can invalidate this approach, leading to apparent associations at markers that are unlinked to disease loci. Family-based tests of association can avoid this problem, but such studies are often more expensive and in some cases--particularly for late-onset diseases--are impractical. In this review article we describe a series of approaches published over the past 2 years which use multilocus genotype data to enable valid case-control tests of association, even in the presence of population structure. These tests can be classified into two categories. "Genomic control" methods use the independent marker loci to adjust the distribution of a standard test statistic, while "structured association" methods infer the details of population structure en route to testing for association. We discuss the statistical issues involved in the different approaches and present results from simulations comparing the relative performance of the methods under a range of models.     

Mathematical properties and bounds on haplotyping populations by pure parsimony

Although the haplotype data can be used to analyze the function of DNA, due to the significant efforts required in collecting the haplotype data, usually the genotype data is collected and then the population haplotype inference (PHI) problem is solved to infer haplotype data from genotype data for a population. This paper investigates the PHI problem based on the pure parsimony criterion (HIPP), which seeks the minimum number of distinct haplotypes to infer a given genotype data. We analyze the mathematical structure and properties for the HIPP problem, propose techniques to reduce the given genotype data into an equivalent one of much smaller size, and analyze the relations of genotype data using a compatible graph. Based on the mathematical properties in the compatible graph, we propose a maximal clique heuristic to obtain an upper bound, and a new polynomial-sized integer linear programming formulation to obtain a lower bound for the HIPP problem.     

Linkage analysis of a complex disease through use of admixed populations

Linkage disequilibrium arising from the recent admixture of genetically distinct populations can be potentially useful in mapping genes for complex diseases. McKeigue has proposed a method that conditions on parental admixture to detect linkage. We show that this method tests for linkage only under specific assumptions, such as equal admixture in the parental generation and admixture that occurs in a single generation. In practice, these assumptions are unlikely to hold for natural populations, resulting in an inflation of the type I error rate when testing for linkage by this method. In this article, we generalize McKeigue's approach of testing for linkage to allow two different admixture models: (1) intermixture admixture and (2) continuous gene flow. We calculate the sample size required for a genomewide search by this method under different disease models: multiplicative, additive, recessive, and dominant. Our results show that the sample size required to obtain 90% power to detect a putative mutant allele at a genomewide significance level of 5% can usually be achieved in practice if informative markers are available at a density of 2 cM.     

Genetic relationships between Indians and their neighboring populations

Using gene frequency data for 18 protein and blood group loci, we studied the genetic relationships of four Indian subcontinent populations (peoples from Punjab, Gujarati, Andhra Pradesh, and Bangladesh) with their neighboring populations (Iranians, Afghans, Sinhalese in Sri Lanka, Nepalese, Bhutanese, Malays, Bataks in northern Sumatra, and Chinese). The results obtained indicate that the four Indian subcontinent populations and the Sinhalese are genetically closer to Iranians and Afghans (Caucasoid) than to the other neighboring Mongoloid populations. Genetic distance analysis shows a clear-cut dichotomy between the Caucasoid and Mongoloid populations.     

Adjustment for local ancestry in genetic association analysis of admixed populations

MOTIVATION: Admixed populations offer a unique opportunity for mapping diseases that have large disease allele frequency differences between ancestral populations. However, association analysis in such populations is challenging because population stratification may lead to association with loci unlinked to the disease locus. Methods and results: We show that local ancestry at a test single nucleotide polymorphism (SNP) may confound with the association signal and ignoring it can lead to spurious association. We demonstrate theoretically that adjustment for local ancestry at the test SNP is sufficient to remove the spurious association regardless of the mechanism of population stratification, whether due to local or global ancestry differences among study subjects; however, global ancestry adjustment procedures may not be effective. We further develop two novel association tests that adjust for local ancestry. Our first test is based on a conditional likelihood framework which models the distribution of the test SNP given disease status and flanking marker genotypes. A key advantage of this test lies in its ability to incorporate different directions of association in the ancestral populations. Our second test, which is computationally simpler, is based on logistic regression, with adjustment for local ancestry proportion. We conducted extensive simulations and found that the Type I error rates of our tests are under control; however, the global adjustment procedures yielded inflated Type I error rates when stratification is due to local ancestry difference.     

Mapping genes that underlie ethnic differences in disease risk: methods for detecting linkage in admixed populations, by conditioning on parental admixture.

Genes that underlie ethnic differences in disease risk can be mapped in affected individuals of mixed descent if the ancestry of the alleles at each marker locus can be assigned to one of the two founding populations. Linkage can be detected by testing for association of the disease with the ancestry of alleles at the marker locus, by conditioning on the admixture (defined as the proportion of genes that have ancestry from the high-risk population) of both parents. With regard to exploiting the effects of admixture, this test is more flexible and powerful than the transmission-disequilibrium test. Under the assumption of a multiplicative model, the statistical power for a given sample size depends only on parental admixture and the risk ratio r between populations that is generated by the locus. The most informative families are those in which mean parental admixture is .2-.7 and in which admixture is similar in both parents. The number of markers required for a genome search depends on the number of generations since admixture and on the information content for ancestry (f) of the markers, defined as a function of allele frequencies in the two founding populations. Simulations using a hidden Markov model suggest that, when admixture has occurred 2-10 generations earlier, a multipoint analysis using 2,000 biallelic markers, with f values of 30%, can extract 70%-90% of the ancestry information for each locus. Sets of such markers could be selected from libraries of single-nucleotide polymorphisms, when these become available.     

Ancestral traits, parental populations, and hybrids

The distance traveled over the centuries by any human population is a rocky road; the fusion of the paths of two meandering populations often leads to an even rougher one. Tracing such a wayward child of history back to its long-lost parents presents complex problems for the physical anthropologist. For those parents may share unequally in their legacy to their offspring at the start of the journey. Many genes may disappear into the gutters of genetic drift; some byways may stop at the barricades of selection; fresh migrants may join by side roads to swell the stream; and inbreeding may narrow the roadbed. Searching for the genes that the hybrid population holds in common with alleged parents provides one clue to the relative role of the ancestors. Comparison of measurements of the living hybrid with the nearest known relatives of the presumed parents yields another clue. Resuscitation of ancestors from their bones and teeth supplies yet another link. Distinctive features of such skeletal remains are most valuable where their inheritance is understood. Bones may yield another clue from substances left within them of known genetic control. And finally all the evidence of history, every glimpse caught of the population along the route, gives knowledge of the contribution of the forebears to the living people. The ideal estimate will utilize a wide variety of information about the known groups.     

Distribution of a 9.1-kb insertion-deletion polymorphism among native and admixed populations from Argentina

We studied five population samples from Argentina, four drawn from Native American groups of the northeast region (Wichí, Pilagá, Toba, and Mbyá-Guaraní) and one from two small villages of the Córdoba province. In this study we report genotypes and allele frequencies of the 9.1-kb insertion-deletion polymorphism located on chromosome 22. The frequency of the deletion allele ranges from 0.276 in the Mbyá-Guaraní to 0.470 in the Pilagá. The coefficient of population differentiation is fairly low (F(ST) = 0.013), does not reflect any geographic or linguistic pattern, and seems to be more related to stochastic processes than to directional forces.     

Analysis of the region V mitochondrial marker in two Black communities of Ecuador,and in their parental populations

Data on the frequency of the mtDNA region V deletion were used to estimate the relative maternal contribution from the parental populations to the gene pools of the two Black communities of Rio Cayapas and Viche in northern Ecuador. Ethnohistorical records and nuclear DNA data indicate that these populations are hybrids of West African and Amerindian populations. The unique distribution of the DNA marker in these parental groups provided good admixture estimates. The fraction of mtDNA of Amerindian origin in the population of Rio Cayapas is quite small (8%±5%), whereas in the community of Viche the native Americans contributed the major portion of the gene pool (51%±15). The mtDNA estimate for Rio Cayapas is similar to that of some protein polymorphisms, which confirms the cultural and genetic isolation of this community from the neighboring native population. On the other hand, the admixture value obtained from nuclear genes in Viche is statistically different from the estimate obtained from mtDNA data. This supports the traditional belief, gathered from historical records and cultural data, that the contribution from Indian females was higher than that of Indian males, at least in the primary settlements of the African-American population of Esmeraldas.     

Robust estimation of local genetic ancestry in admixed populations using a nonparametric bayesian approach

We present a new haplotype-based approach for inferring local genetic ancestry of individuals in an admixed population. Most existing approaches for local ancestry estimation ignore the latent genetic relatedness between ancestral populations and treat them as independent. In this article, we exploit such information by building an inheritance model that describes both the ancestral populations and the admixed population jointly in a unified framework. Based on an assumption that the common hypothetical founder haplotypes give rise to both the ancestral and the admixed population haplotypes, we employ an infinite hidden Markov model to characterize each ancestral population and further extend it to generate the admixed population. Through an effective utilization of the population structural information under a principled nonparametric Bayesian framework, the resulting model is significantly less sensitive to the choice and the amount of training data for ancestral populations than state-of-the-art algorithms. We also improve the robustness under deviation from common modeling assumptions by incorporating population-specific scale parameters that allow variable recombination rates in different populations. Our method is applicable to an admixed population from an arbitrary number of ancestral populations and also performs competitively in terms of spurious ancestry proportions under a general multiway admixture assumption. We validate the proposed method by simulation under various admixing scenarios and present empirical analysis results from a worldwide-distributed dataset from the Human Genome Diversity Project.     

Starting the invasion pathway: the interaction between source populations and human transport vectors

Human transport hubs, such as shipping ports, airports and mail centers are important foci for the spread of non-indigenous species. High relative abundance in a transport hub has been proposed as a correlate of invasion success, since abundant species are thought more likely to colonize vectors and to be transported more frequently. We here evaluate the relative importance of vector characteristics and local source assemblages in determining the pool of species that is transported by hull fouling on recreational boats. We compared the resident fouling communities of three recreational boat harbors in Australia with the assemblages on the hulls of boats that travel between them. We used data on the recent travel and maintenance history of the boats to evaluate correlates of transport probability and the potential for intra-coastal spread of fouling organisms. Invertebrate assemblages on heavily fouled vessels reflected the composition of biotic assemblages within the marina in which they were moored, but by itself, relative abundance in the source port was not a reliable predictor of transport probability. More important was the age of the antifouling paint on the vessels’ hulls, which acted selectively on some groups of organisms. Movements of vessels were characterized by “fidelity” (vessels remaining close to homeport) interspersed with “promiscuity” (vessels traveling to a diverse pool of destinations). In an infested harbor, measures taken to increase the resistance of vectors to colonization by the invader should be effective in slowing the rate of spread to other locations, by decreasing the overall frequency of transport.     

Gene flow in admixed populations and implications for the conservation of the Western honeybee, Apis mellifera

Anthropogenic activity, especially modern apiculture, has considerable impact on the natural distribution of the Western honeybee, Apis mellifera, leading to the spread, replacement and fragmentation of many subspecies. This creates demand for the conservation of some subspecies, in particular, Apis mellifera mellifera, which once was widely distributed in Western Europe and nowadays is endangered through habitat loss and fragmentation. Moreover, A. m. mellifera may be further endangered by hybridisation in populations that now occur in artificial sympatry with other subspecies. Here, we quantify and compare individual hybridisation between sympatric and allopatric honeybee populations of A. m. mellifera and A. m. carnica using microsatellite markers and a Bayesian model-based approach. We had a special focus on pure breeding populations, which are a major tool in honeybee conservation. Our results demonstrate that subspecies are still highly differentiated, but gene flow is not prevented by the current management strategies, creating urgent demand for an improved conservation management of A. m. mellifera. However, the occurrence of a high number of pure individuals might suggest that some sort of hybrid barrier acts against the complete admixture of the two subspecies.     

Population structure in admixed populations: effect of admixture dynamics on the pattern of linkage disequilibrium

Gene flow between genetically distinct populations creates linkage disequilibrium (admixture linkage disequilibrium [ALD]) among all loci (linked and unlinked) that have different allele frequencies in the founding populations. We have explored the distribution of ALD by using computer simulation of two extreme models of admixture: the hybrid-isolation (HI) model, in which admixture occurs in a single generation, and the continuous-gene-flow (CGF) model, in which admixture occurs at a steady rate in every generation. Linkage disequilibrium patterns in African American population samples from Jackson, MS, and from coastal South Carolina resemble patterns observed in the simulated CGF populations, in two respects. First, significant association between two loci (FY and AT3) separated by 22 cM was detected in both samples. The retention of ALD over relatively large (>10 cM) chromosomal segments is characteristic of a CGF pattern of admixture but not of an HI pattern. Second, significant associations were also detected between many pairs of unlinked loci, as observed in the CGF simulation results but not in the simulated HI populations. Such a high rate of association between unlinked markers in these populations could result in false-positive linkage signals in an admixture-mapping study. However, we demonstrate that by conditioning on parental admixture, we can distinguish between true linkage and association resulting from shared ancestry. Therefore, populations with a CGF history of admixture not only are appropriate for admixture mapping but also have greater power for detection of linkage disequilibrium over large chromosomal regions than do populations that have experienced a pattern of admixture more similar to the HI model, if methods are employed that detect and adjust for disequilibrium caused by continuous admixture.     

Genetic signatures of parental contribution in black and white populations in Brazil

Two hundred and three individuals classified as white were tested for 11 single nucleotide polymorphisms plus two insertion/deletions in their Y-chromosomes. A subset of these individuals (n = 172) was also screened for sequences in the first hypervariable segment of their mitochondrial DNA (mtDNA). In addition, complementary studies were done for 11 of the 13 markers indicated above in 54 of 107 black subjects previously investigated in this southern Brazilian population. The prevalence of Y-chromosome haplogroups among whites was similar to that found in the Azores (Portugal) or Spain, but not to that of other European countries. About half of the European or African mtDNA haplogroups of these individuals were related to their places of origin, but not their Amerindian counterparts. Persons classified in these two categories of skin color and related morphological traits showed distinct genomic ancestries through the country. These findings emphasize the need to consider in Brazil, despite some general trends, a notable heterogeneity in the pattern of admixture dynamics within and between populations/groups.