Ethnic-difference markers for use in mapping by admixture linkage disequilibrium

Mapping by admixture linkage disequilibrium (MALD) is a potentially powerful technique for the mapping of complex genetic diseases. The practical requirements of this method include (a) a set of markers spanning the genome that have large allele-frequency differences between the parental ethnicities contributing to the admixed population and (b) an understanding of the extent of admixture in the study population. To this end, a DNA-pooling technique was used to screen microsatellite and diallelic insertion/deletion markers for allele-frequency differences between putative representatives of the parental populations of the admixed Mexican American (MA) and African American (AA) populations. Markers with promising pooled differences were then confirmed by individual genotyping in both the parental and admixed populations. For the MA population, screening of >600 markers identified 151 ethnic-difference markers (EDMs) with delta>0.30 (where delta is the absolute value of each allele-frequency difference between two populations, summed over all marker alleles and divided by two) that are likely to be useful for MALD analysis. For the AA population, analysis of >400 markers identified 97 EDMs. In addition, individual genotyping of these markers in Pima Amerindians, Yavapai Amerindians, European American (EA) individuals, Africans from Zimbabwe, MA individuals, and AA individuals, as well as comparison to the CEPH genotyping set, suggests that the differences between subpopulations of an ethnicity are small for many markers with large interethnic differences. Estimates of admixture that are based on individual genotyping of these markers are consistent with a 60% EA:40% Amerindian contribution to MA populations and with a 20% EA:80% African contribution to AA populations. Taken together, these data suggest that EDMs with large interpopulation and small intrapopulation differences can be readily identified for MALD studies in both AA and MA populations.     

Marker selection for the transmission/disequilibrium test, in recently admixed populations.

Recent admixture between genetically differentiated populations can result in high levels of association between alleles at loci that are <=10 cM apart. The transmission/disequilibrium test (TDT) proposed by Spielman et al. (1993) can be a powerful test of linkage between disease and marker loci in the presence of association and therefore could be a useful test of linkage in admixed populations. The degree of association between alleles at two loci depends on the differences in allele frequencies, at the two loci, in the founding populations; therefore, the choice of marker is important. For a multiallelic marker, one strategy that may improve the power of the TDT is to group marker alleles within a locus, on the basis of information about the founding populations and the admixed population, thereby collapsing the marker into one with fewer alleles. We have examined the consequences of collapsing a microsatellite into a two-allele marker, when two founding populations are assumed for the admixed population, and have found that if there is random mating in the admixed population, then typically there is a collapsing for which the power of the TDT is greater than that for the original microsatellite marker. A method is presented for finding the optimal collapsing that has minimal dependence on the disease and that uses estimates either of marker allele frequencies in the two founding populations or of marker allele frequencies in the current, admixed population and in one of the founding populations. Furthermore, this optimal collapsing is not always the collapsing with the largest difference in allele frequencies in the founding populations. To demonstrate this strategy, we considered a recent data set, published previously, that provides frequency estimates for 30 microsatellites in 13 populations.     

Inference of population structure using multilocus genotype data

We describe a model-based clustering method for using multilocus genotype data to infer population structure and assign individuals to populations. We assume a model in which there are K populations (where K may be unknown), each of which is characterized by a set of allele frequencies at each locus. Individuals in the sample are assigned (probabilistically) to populations, or jointly to two or more populations if their genotypes indicate that they are admixed. Our model does not assume a particular mutation process, and it can be applied to most of the commonly used genetic markers, provided that they are not closely linked. Applications of our method include demonstrating the presence of population structure, assigning individuals to populations, studying hybrid zones, and identifying migrants and admixed individuals. We show that the method can produce highly accurate assignments using modest numbers of loci-e.g. , seven microsatellite loci in an example using genotype data from an endangered bird species. The software used for this article is available from http://www.stats.ox.ac.uk/ approximately pritch/home. html.     

Mathematical properties of Fst between admixed populations and their parental source populations

We consider the properties of the F(st) measure of genetic divergence between an admixed population and its parental source populations. Among all possible populations admixed among an arbitrary set of parental populations, we show that the value of F(st) between an admixed population and a specific source population is maximized when the admixed population is simply the most distant of the other source populations. For the case with only two parental populations, as a function of the admixture fraction, we further demonstrate that this F(st) value is monotonic and convex, so that F(st) is informative about the admixture fraction. We illustrate our results using example human population-genetic data, showing how they provide a framework in which to interpret the features of F(st) in admixed populations.     

Polymorphic Admixture Typing in Human Ethnic Populations

A panel of 257 RFLP loci was selected on the basis of high heterozygosity in Caucasian DNA surveys and equivalent spacing throughout the human genome. Probes from each locus were used in a Southern blot survey of allele frequency distribution for four human ethnic groups: Caucasian, African American, Asian (Chinese), and American Indian (Cheyenne). Nearly all RFLP loci were polymorphic in each group, albeit with a broad range of differing allele frequencies (δ). The distribution of frequency differences (δ values) was used for three purposes: (1) to provide estimates for genetic distance (differentiation) among these ethnic groups, (2) to revisit with a large data set the proportion of human genetic variation attributable to differentiation within ethnic groups, and (3) to identify loci with high δ values between recently admixed populations of use in mapping by admixture linkage disequilibrium (MALD). Although most markers display significant allele frequency differences between ethnic groups, the overall genetic distances between ethnic groups were small (.066–.098), and <10% of the measured overall molecular genetic diversity in these human samples can be attributed to “racial” differentiation. The median δ values for pairwise comparisons between groups fell between .15 and .20, permitting identification of highly informative RFLP loci for MALD disease association studies.     

HIV-1 Vaccine-Induced T-Cell Reponses Cluster in Epitope Hotspots that Differ from Those Induced in Natural Infection with HIV-1

Several recent large clinical trials evaluated HIV vaccine candidates that were based on recombinant adenovirus serotype 5 (rAd-5) vectors expressing HIV-derived antigens. These vaccines primarily elicited T-cell responses, which are known to be critical for controlling HIV infection. In the current study, we present a meta-analysis of epitope mapping data from 177 participants in three clinical trials that tested two different HIV vaccines: MRKAd-5 HIV and VRC-HIVAD014-00VP. We characterized the population-level epitope responses in these trials by generating population-based epitope maps, and also designed such maps using a large cohort of 372 naturally infected individuals. We used these maps to address several questions: (1) Are vaccine-induced responses randomly distributed across vaccine inserts, or do they cluster into immunodominant epitope hotspots? (2) Are the immunodominance patterns observed for these two vaccines in three vaccine trials different from one another? (3) Do vaccine-induced hotspots overlap with epitope hotspots induced by chronic natural infection with HIV-1? (4) Do immunodominant hotspots target evolutionarily conserved regions of the HIV genome? (5) Can epitope prediction methods be used to identify these hotspots? We found that vaccine responses clustered into epitope hotspots in all three vaccine trials and some of these hotspots were not observed in chronic natural infection. We also found significant differences between the immunodominance patterns generated in each trial, even comparing two trials that tested the same vaccine in different populations. Some of the vaccine-induced immunodominant hotspots were located in highly variable regions of the HIV genome, and this was more evident for the MRKAd-5 HIV vaccine. Finally, we found that epitope prediction methods can partially predict the location of vaccine-induced epitope hotspots. Our findings have implications for vaccine design and suggest a framework by which different vaccine candidates can be compared in early phases of evaluation.     

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.     

Genetic variation in North Amerindian populations: association with sociocultural complexity

A survey of nine polymorphic loci for 82 North Amerindian populations was undertaken to test the hypothesis that increasing levels of sociocultural complexity are ineluctably accompanied by increased heterozygosity. The data reveal a significant relationship in the predicted direction. Moreover, the significant correlation between average heterozygosity and sociocultural complexity is substantially increased by the removal of 19 highly admixed samples. However, this relationship, at least among North Amerindian populations, may be more apparent than real since both mean heterozygosity and the level of sociocultural organization are significantly negatively correlated with latitude. When this latter variable is controlled for, all correlations between heterozygosity and sociocultural complexity are rendered nonsignificant.     

Hemochromatosis gene variants in three different ethnic populations: effects of admixture for screening programs

Genetic testing for hemochromatosis may have important implications for diagnosis and screening of the disease. However, the relative importance of mutations in the gene for hereditary hemochromatosis, HFE, may vary among populations, when the mutant allele frequencies and their penetrance in a particular genetic and environmental background are taken into account. We present data on the allele and genotype frequencies and population structure of two HFE genetic variants in three different ethnic groups from a highly mixed urban population (S?o Paulo, Brazil). Allele frequencies for both the C282Y and H63D HFE mutations showed significant differences among the studied populations (for the C282Y mutation, Euro-Brazilian 3.7%, admixed 0.7%, Afro-Brazilian 0.5%; and for the H63D mutation, Euro-Brazilian 20.3%, admixed 13.0%, Afro-Brazilian 6.4). The data substantiate a European origin for these mutations. Furthermore, they provide a basis for a more rational strategic planning of population screening programs for the disease.     

Fat patterning in Papua New Guineans: effects of age, sex and acculturation

For the renewed interest in fat patterning to be fruitful, more information on the biology of fat patterning, particularly normative data are required. Nine skinfold thicknesses were measured in 2312 coastal and highland Papua New Guineans of all ages undergoing acculturation. Principal component analysis of the skinfold data revealed components of fatness and fat patterning. First component scores were typical of the changes in fatness with age and the sex differences expected in these populations. Second component scores of limb-trunk fat patterning showed marked changes with age in both sexes but they stabilised in adulthood. This pattern seemed more affected by maturity although it is usually thought to be a sex-specific pattern. The third component, upper-lower body patterning, was more sexually dimorphic in adults. Marked differences were not observed between the two communities but over a 15-year period fatness levels increased slightly and subcutaneous fat tended to be deposited more centrally.     

The Search for Cancer Biomarkers: Assessing the Distribution of INDEL Markers in Different Genetic Ancestries

Cancer is a multifactorial group of diseases, being highly incident and one of the leading causes of death worldwide. In Brazil, there is a great variation in cancer incidence and impact among the different geographic regions, partly due to the genetic heterogeneity of the population in this country, composed mainly by European (EUR), Native American (NAM), African (AFR), and Asian (ASN) ancestries. Among different populations, genetic markers commonly present diverse allelic frequencies, but in admixed populations, such as the Brazilian population, data is still limited, which is an issue that might influence cancer incidence. Therefore, we analyzed the allelic and genotypic distribution of 12 INDEL polymorphisms of interest in populations from the five Brazilian geographic regions and in populations representing EUR, NAM, AFR, and ASN, as well as tissue expression in silico. Genotypes were obtained by multiplex PCR and the statistical analyses were done using R, while data of tissue expression for each marker was extracted from GTEx portal. We highlight that all analyzed markers presented statistical differences in at least one of the population comparisons, and that we found 39 tissues to be differentially expressed depending on the genotype. Here, we point out the differences in genotype distribution and gene expression of potential biomarkers for risk of cancer development and we reinforce the importance of this type of study in populations with different genetic backgrounds.     

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

ABSTRACT: BACKGROUND: Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated ('allopatric') populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations. RESULTS: As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation. CONCLUSIONS: Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of 'genomic islands of speciation'.     

Implications of the Admixture Process in Skin Color Molecular Assessment

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.     

Dermatoglyphic changes during the population admixture between Kam and Han Chinese

Genetic studies and gene localization for human dermatoglyphs are currently ongoing. However, the inheritance modes of various genetic traits are not well understood because of the complexity of dermatoglyph genetics. The study of admixed populations can contribute to the understanding of population genetic traits of dermatoglyphs. Here, we present the dermatoglyphic characteristics of Kam and Liujia Han, and the admixed population consisting of these two parent populations.The characteristics of the admixed population do not always fall in the same ranges as the parent population characters but do seem to be biased to Kam or Liujia parent populations, depending on sex and ethnicity of parents. The total frequencies of different fingerprint types do not differ among these populations, but several of the quantitative traits and the palm true pattern frequencies do significantly differ between admixed and parent populations. The simple arch fingerprint frequency decreases significantly in the admixed population in comparison with parent populations while both simple whorl fingerprint frequency and finger ridge counts increase significantly. True pattern frequency of the span area of interdigital III and IV areas on right hands and the radial-loop frequency of the right index fingers in the admixed populations are consistent with their matrilineal population. These dermatoglyphic changes may result from increased heterozygosity in the admixed population. The genetic modes of these changes may be relatively simple and will be useful for future dermatoglyph genetic studies.      

Estimates of hierarchical variation in flower morphology in natural populations ofScleranthus annuus (Caryophyllaceae), an inbreeding annual

The variation in two sets of morphological characters of the flowers of the highly inbreedingScleranthus annuus (Caryophyllaceae) was assessed using 15–20 plants from each of 20 natural populations from the southernmost region of Sweden. The stamen fertility data set consisted of 10 characters describing the degree of fertility of the ten stamens/staminoids, while the sepal/gynoecium data set comprised 10 sepal characters, together with style and stigma length. Substantial variation was found in both the degree of development and the fertility of the stamens, the degree of variability in stamen fertility being related to stamen position within the flower. Considerable variation was found in the characters of the sepal and the gynoecium. Hierarchical analyses of variance indicated that 29% of the variation in total male reproductive effort was distributed among populations, 28% among plants within populations and 43% represented within-plant variation. The corresponding averages for the characters from the sepal/gynoecium data set are 26, 38, and 35%: a greater proportion of the total variance in female reproductive characters is accounted for by among individual variation than is the case with the male reproductive characters. Significance tests of Mahalanobis distances derived by canonical variate analyses indicated that all populations were significantly separated using the sepal/gynoecium data set, while only 50% of the pairwise comparisons on the basis of the stamen fertility data set were significant. Cluster analysis did not reveal any aggregation of the populations. The incongruence of the two data sets and their ability to discriminate between the populations is 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.     

The Admixture Structure and Genetic Variation of the Archipelago of Cape Verde and Its Implications for Admixture Mapping Studies

Recently admixed populations offer unique opportunities for studying human history and for elucidating the genetic basis of complex traits that differ in prevalence between human populations. Historical records, classical protein markers, and preliminary genetic data indicate that the Cape Verde islands in West Africa are highly admixed and primarily descended from European males and African females. However, little is known about the variation in admixture levels, admixture dynamics and genetic diversity across the islands, or about the potential of Cape Verde for admixture mapping studies. We have performed a detailed analysis of phenotypic and genetic variation in Cape Verde based on objective skin color measurements, socio-economic status (SES) evaluations and data for 50 autosomal, 34 X-chromosome, and 21 non-recombinant Y-chromosome (NRY) markers in 845 individuals from six islands of the archipelago. We find extensive genetic admixture between European and African ancestral populations (mean West African ancestry = 0.57, sd = 0.08), with individual African ancestry proportions varying considerably among the islands. African ancestry proportions calculated with X and Y-chromosome markers confirm that the pattern of admixture has been sex-biased. The high-resolution NRY-STRs reveal additional patterns of variation among the islands that are most consistent with differentiation after admixture. The differences in the autosomal admixture proportions are clearly evident in the skin color distribution across the islands (Pearson r = 0.54, P-value<2e–16). Despite this strong correlation, there are significant interactions between SES and skin color that are independent of the relationship between skin color and genetic ancestry. The observed distributions of admixture, genetic variation and skin color and the relationship of skin color with SES relate to historical and social events taking place during the settlement history of Cape Verde, and have implications for the design of association studies using this population.     

Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.     

Atlas versus range maps: robustness of chorological relationships to distribution data types in European mammals

Aim Chorological relationships describe the patterns of distributional overlap among species. In addition to revealing biogeographical structure, the resulting clusters of species with similar geographical distributions can serve as natural units in conservation planning. Here, we assess the extent to which temporal, methodological and taxonomical differences in the source of species’ distribution data can affect the relationships that are found. Location Western Europe. Methods We used two data sets – the Atlas of European mammals and polygon range maps from the IUCN Global Mammal Assessment – both as presence–absence data for UTM 50 km × 50 km squares. We performed pairwise comparisons among 156 species for each data set to build matrices of the similarity in distribution across species, using both Jaccard’s and Baroni‐Urbani & Buser’s indices. We then compared these similarity matrices (chorological relationships), as well as the species richness and occurrence patterns from the two data sets. Results As expected, range maps increased both the mean prevalence per species and mean species richness per grid cell in comparison to atlas data, reflecting the general view that these data types respectively over‐ and underestimate species occurrence. However, species richness and occurrence patterns in atlas and range map data were positively associated and, most importantly, the chorological relationships underlying the two data sets were highly similar. Main conclusions Despite many methodological, temporal and taxonomical differences between atlas data and range maps, the chorological relationships encountered between species were similar for both data sets. Chorological analyses can thus be robust to the data source used and provide a solid basis for analytical biogeographical studies, even over broad spatial scales.     

A genomewide admixture mapping panel for Hispanic/Latino populations

Admixture mapping (AM) is a promising method for the identification of genetic risk factors for complex traits and diseases showing prevalence differences among populations. Efficient application of this method requires the use of a genomewide panel of ancestry-informative markers (AIMs) to infer the population of origin of chromosomal regions in admixed individuals. Genomewide AM panels with markers showing high frequency differences between West African and European populations are already available for disease-gene discovery in African Americans. However, no such a map is yet available for Hispanic/Latino populations, which are the result of two-way admixture between Native American and European populations or of three-way admixture of Native American, European, and West African populations. Here, we report a genomewide AM panel with 2,120 AIMs showing high frequency differences between Native American and European populations. The average intermarker genetic distance is ~1.7 cM. The panel was identified by genotyping, with the Affymetrix GeneChip Human Mapping 500K array, a population sample with European ancestry, a Mesoamerican sample comprising Maya and Nahua from Mexico, and a South American sample comprising Aymara/Quechua from Bolivia and Quechua from Peru. The main criteria for marker selection were both high information content for Native American/European ancestry (measured as the standardized variance of the allele frequencies, also known as "f value") and small frequency differences between the Mesoamerican and South American samples. This genomewide AM panel will make it possible to apply AM approaches in many admixed populations throughout the Americas.