Detection of Ancestry Informative HLA Alleles Confirms the Admixed Origins of Japanese Population

The polymorphisms in the human leukocyte antigen (HLA) region are powerful tool for studying human evolutionary processes. We investigated genetic structure of Japanese by using five-locus HLA genotypes (HLA-A, -B, -C, -DRB1, and -DPB1) of 2,005 individuals from 10 regions of Japan. We found a significant level of population substructure in Japanese; particularly the differentiation between Okinawa Island and mainland Japanese. By using a plot of the principal component scores, we identified ancestry informative alleles associated with the underlying population substructure. We examined extent of linkage disequilibrium (LD) between pairs of HLA alleles on the haplotypes that were differentiated among regions. The LDs were strong and weak for pairs of HLA alleles characterized by low and high frequencies in Okinawa Island, respectively. The five-locus haplotypes whose alleles exhibit strong LD were unique to Japanese and South Korean, suggesting that these haplotypes had been recently derived from the Korean Peninsula. The alleles characterized by high frequency in Japanese compared to South Korean formed segmented three-locus haplotype that was commonly found in Aleuts, Eskimos, and North- and Meso-Americans but not observed in Korean and Chinese. The serologically equivalent haplotype was found in Orchid Island in Taiwan, Mongol, Siberia, and Arctic regions. It suggests that early Japanese who existed prior to the migration wave from the Korean Peninsula shared ancestry with northern Asian who moved to the New World via the Bering Strait land bridge. These results may support the admixture model for peopling of Japanese Archipelago.     

Haplotypes in SLC24A5 Gene as Ancestry Informative Markers in Different Populations

Ancestry informative markers (AIMs) are human polymorphisms that exhibit substantially allele frequency differences among populations. These markers can be useful to provide information about ancestry of samples which may be useful in predicting a perpetrator’s ethnic origin to aid criminal investigations. Variations in human pigmentation are the most obvious phenotypes to distinguish individuals. It has been recently shown that the variation of a G in an A allele of the coding single-nucleotide polymorphism (SNP) rs1426654 within SLC24A5 gene varies in frequency among several population samples according to skin pigmentation. Because of these observations, the SLC24A5 locus has been evaluated as Ancestry Informative Region (AIR) by typing rs1426654 together with two additional intragenic markers (rs2555364 and rs16960620) in 471 unrelated individuals originating from three different continents (Africa, Asia and Europe). This study further supports the role of human SLC24A5 gene in skin pigmentation suggesting that variations in SLC24A5 haplotypes can correlate with human migration and ancestry. Furthermore, our data do reveal the utility of haplotype and combined unphased genotype analysis of SLC24A5 in predicting ancestry and provide a good example of usefulness of genetic characterization of larger regions, in addition to single polymorphisms, as candidates for population-specific sweeps in the ancestral population.     

Genome-wide analysis of the structure of the South African Coloured Population in the Western Cape

Admixed populations present unique opportunities to discover the genetic factors underlying many multifactorial diseases. The geographical position and complex history of South Africa has led to the establishment of the unique admixed population known as the South African Coloured. Not much is known about the genetic make-up of this population, and the historical record is patchy. We genotyped 959 individuals from the Western Cape area, self-identified as belonging to this population, using the Affymetrix 500k genotyping platform. This resulted in nearly 75,000 autosomal SNPs that could be compared with populations represented in the International HapMap Project and the Human Genome Diversity Project. Analysis by means of both the admixture and linkage models in STRUCTURE revealed that the major ancestral components of this population are predominantly Khoesan (32–43%), Bantu-speaking Africans (20–36%), European (21–28%) and a smaller Asian contribution (9–11%), depending on the model used. This is consistent with historical data. While of great historical and genealogical interest, this information is also essential for future admixture mapping of disease genes in this population.     

A Continuous Correlated Beta Process Model for Genetic Ancestry in Admixed Populations

Admixture and recombination create populations and genomes with genetic ancestry from multiple source populations. Analyses of genetic ancestry in admixed populations are relevant for trait and disease mapping, studies of speciation, and conservation efforts. Consequently, many methods have been developed to infer genome-average ancestry and to deconvolute ancestry into continuous local ancestry blocks or tracts within individuals. Current methods for local ancestry inference perform well when admixture occurred recently or hybridization is ongoing, or when admixture occurred in the distant past such that local ancestry blocks have fixed in the admixed population. However, methods to infer local ancestry frequencies in isolated admixed populations still segregating for ancestry do not exist. In the current paper, I develop and test a continuous correlated beta process model to fill this analytical gap. The method explicitly models autocorrelations in ancestry frequencies at the population-level and uses discriminant analysis of SNP windows to take advantage of ancestry blocks within individuals. Analyses of simulated data sets show that the method is generally accurate such that ancestry frequency estimates exhibited low root-mean-square error and were highly correlated with the true values, particularly when large (±10 or ±20) SNP windows were used. Along these lines, the proposed method outperformed post hoc inference of ancestry frequencies from a traditional hidden Markov model (i.e., the linkage model in structure), particularly when admixture occurred more distantly in the past with little on-going gene flow or was followed by natural selection. The reliability and utility of the method was further assessed by analyzing genetic ancestry in an admixed human population (Uyghur) and three populations from a hybrid zone between Mus domesticus and M. musculus. Considerable variation in ancestry frequencies was detected within and among chromosomes in the Uyghur, with a large region of excess French ancestry harboring a gene with a known disease association. Similar variation was detected in the mouse hybrid zone, with notable constancy in regions of excess ancestry among admixed populations. By filling what has been an analytical gap, the proposed method should be a useful tool for many biologists. A computer program (popanc), written in C++, has been developed based on the proposed method and is available on-line at http://sourceforge.net/projects/popanc/.     

Finding Markers That Make a Difference: DNA Pooling and SNP-Arrays Identify Population Informative Markers for Genetic Stock Identification

Genetic stock identification (GSI) using molecular markers is an important tool for management of migratory species. Here, we tested a cost-effective alternative to individual genotyping, known as allelotyping, for identification of highly informative SNPs for accurate genetic stock identification. We estimated allele frequencies of 2880 SNPs from DNA pools of 23 Atlantic salmon populations using Illumina SNP-chip. We evaluated the performance of four common strategies (global F _ST, pairwise F _ST, Delta and outlier approach) for selection of the most informative set of SNPs and tested their effectiveness for GSI compared to random sets of SNP and microsatellite markers. For the majority of cases, SNPs selected using the outlier approach performed best followed by pairwise F _ST and Delta methods. Overall, the selection procedure reduced the number of SNPs required for accurate GSI by up to 53% compared with randomly chosen SNPs. However, GSI accuracy was more affected by populations in the ascertainment group rather than the ranking method itself. We demonstrated for the first time the compatibility of different large-scale SNP datasets by compiling the largest population genetic dataset for Atlantic salmon to date. Finally, we showed an excellent performance of our top SNPs on an independent set of populations covering the main European distribution range of Atlantic salmon. Taken together, we demonstrate how combination of DNA pooling and SNP arrays can be applied for conservation and management of salmonids as well as other species.     

Using multi-way admixture mapping to elucidate TB susceptibility in the South African Coloured population

Background

The admixed South African Coloured population is ideally suited to the discovery of tuberculosis susceptibility genetic variants and their probable ethnic origins, but previous attempts at finding such variants using genome-wide admixture mapping were hampered by the inaccuracy of local ancestry inference. In this study, we infer local ancestry using the novel algorithm implemented in RFMix, with the emphasis on identifying regions of excess San or Bantu ancestry, which we hypothesize may harbour TB susceptibility genes.

Results

Using simulated data, we demonstrate reasonable accuracy of local ancestry inference by RFMix, with a tendency towards miss-calling San ancestry as Bantu. Regions with either excess San ancestry or excess African (San or Bantu) ancestry are less likely to be affected by this bias, and we therefore proceeded to identify such regions, found in cases but not in controls (642 cases and 91 controls). A number of promising regions were found (overall p-values of 7.19×10^-5 for San ancestry and <2.00×10^-16 for African ancestry), including chromosomes 15q15 and 17q22, which are close to genomic regions previously implicated in TB. Promising immune-related susceptibility genes such as the GADD45A, OSM and B7-H5 genes are also harboured in the identified regions.

Conclusion

Admixture mapping is feasible in the South African Coloured population and a number of novel TB susceptibility genomic regions were uncovered.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1021) contains supplementary material, which is available to authorized users.     

Genomic Ancestry,Self-Rated Health and Its Association with Mortality in an Admixed Population: 10 Year Follow-Up of the Bambui-Epigen (Brazil) Cohort Study of Ageing

Background

Self-rated health (SRH) has strong predictive value for mortality in different contexts and cultures, but there is inconsistent evidence on ethnoracial disparities in SRH in Latin America, possibly due to the complexity surrounding ethnoracial self-classification.

Materials/Methods

We used 370,539 Single Nucleotide Polymorphisms (SNPs) to examine the association between individual genomic proportions of African, European and Native American ancestry, and ethnoracial self-classification, with baseline and 10-year SRH trajectories in 1,311 community dwelling older Brazilians. We also examined whether genomic ancestry and ethnoracial self-classification affect the predictive value of SRH for subsequent mortality.

Results

European ancestry predominated among participants, followed by African and Native American (median = 84.0%, 9.6% and 5.3%, respectively); the prevalence of Non-White (Mixed and Black) was 39.8%. Persons at higher levels of African and Native American genomic ancestry, and those self-identified as Non-White, were more likely to report poor health than other groups, even after controlling for socioeconomic conditions and an array of self-reported and objective physical health measures. Increased risks for mortality associated with worse SRH trajectories were strong and remarkably similar (hazard ratio ~3) across all genomic ancestry and ethno-racial groups.

Conclusions

Our results demonstrated for the first time that higher levels of African and Native American genomic ancestry—and the inverse for European ancestry—were strongly correlated with worse SRH in a Latin American admixed population. Both genomic ancestry and ethnoracial self-classification did not modify the strong association between baseline SRH or SRH trajectory, and subsequent mortality.     

Population structure and morphology of South African deep-water parapagurid hermit crabs

The parapagurid hermit crabs Sympagurus dimorphus (Studer, 1883) and Parapagurus bouvieri (Stebbing, 1910) form dominant components of South African deeper-water benthic invertebrate communities and, as adults, forego the use of “traditional” gastropod shells in favour of symbiotically associated Epizoanthus carcinoecia. A total of 1,122 adult S. dimorphus and 549 adult P. bouvieri derived from demersal research trawl surveys were examined to investigate parameters of population structure, carcinoecia association and morphology. S. dimorphus had a female-biased sex ratio of 1:1.5, and P. bouvieri an equal sex ratio of 1:1. The proportion of males of both species increased with cephalothoracic shield length (CL). The lengths of the major right chelipeds of both S. dimorphus males and females, and P. bouvieri males were positively allometric, and negatively allometric in P. bouvieri females. All P. bouvieri inhabited carcinoecia, while 3.39% (n = 38) of S. dimorphus specimens occupied gastropod shells. There was a positive relationship between S. dimorphus CL and rate of gastropod shell occupation (r² = 0.96). In S. dimorphus, number of zoanthid polyps forming the carcinoecia pseudoshell increased with both total shell volume (R² = 0.73) and colony height (R² = 0.61), but no such relationships were found for P. bouvieri. Euspira napus was the most common original shell around which carcinoecia developed in both S. dimorphus (45.45%) and P. bouvieri (39.00%). We speculate that these hermit species exist along an evolutionary continuum of carcinoecia use. Sympagurus dimorphus retains the ability to occupy both gastropod shells and pseudoshells and, like most other hermit crabs, continue to use shells/carcinoecia as shelters to protect the entire body. In contrast, adult P. bouvieri appear to be totally reliant on carcinoecia, which they use more as a brood pouch covering only the abdomen, rather than as a shelter into which they can fully retract to escape predation.     

Genetic Variability in Markers of HLA-C Expression in Two Diverse South African Populations

An insertion-deletion (indel) polymorphism within the 3′ untranslated region (UTR) of HLA-C has been shown to be involved in the regulation of HLA-C expression. Individuals who carry a deletion at this position exhibit increased HLA-C expression, which associates with lower viral set point in HIV-1 infected individuals. This 263 indel (rs67384697) is reported to be in strong linkage disequilibrium (LD) with a single nucleotide polymorphism (SNP) 35 kilobases upstream of HLA-C (-35T/C; rs9264942) in Caucasian individuals, making this SNP a potential marker for both HLA-C expression and HIV-1 disease progression. We therefore examined genetic variation within the HLA-C 3′ UTR of 265 Black and Caucasian South Africans by direct sequencing and identified haplotypes encompassing the 263 indel and another indel at position 230 in both populations. Concomitant evaluation of variability at the −35 SNP revealed this polymorphism to be an inappropriate marker for the 263 indel in these populations. These findings provide important insights into genetic variability within the regulatory regions of HLA-C that have potential implications for our understanding of the regulation of HLA-C expression and its impact on HIV-1 disease progression.     

A genetic linkage map for the South African Angora goat

Despite their economical importance, relatively few molecular studies have been made on goats compared to other livestock species. The most recent goat map was published in 1998, and lacks complete genome coverage. A large number of discrepancies and especially inter-chromosomal re-assignments were reported between the 1998 goat linkage map and the sheep map. In this study 94 microsatellite markers were amplified in 12 half-sib South African Angora goat families for compilation of a genetic map, aiming to confirm or reject previously reported rearrangements and to improve the alignment between the ovine and caprine maps. The number of informative meiosis per marker ranged from 69 to 836, with an average of 518. The microsatellites were mapped to 23 chromosomes, spanning 1352 cM and resulting in an average marker interval of 23.0 cM. Marker orders were compared to the previously published goat maps, as well as to the ovine map. Six chromosomes (CHI 2, 4, 5, 11, 13 and 19) showed rearrangements in marker order compared to the 1998 Schibler et al. goat map, while nine previously unmapped markers were conclusively assigned to eight chromosomes. Four of the previously reported intra-chromosomal rearrangements between the goat and sheep maps were confirmed to be either population specific or mapping errors. The verification of rearrangements in loci order will lead to improved alignment between the two maps, as well as improved efficiency of genome and fine mapping efforts in goats.     

Unraveling the Complex Hybrid Ancestry and Domestication History of Cultivated Strawberry

Cultivated strawberry (Fragaria × ananassa) is one of our youngest domesticates, originating in early eighteenth-century Europe from spontaneous hybrids between wild allo-octoploid species (Fragaria chiloensis and Fragaria virginiana). The improvement of horticultural traits by 300 years of breeding has enabled the global expansion of strawberry production. Here, we describe the genomic history of strawberry domestication from the earliest hybrids to modern cultivars. We observed a significant increase in heterozygosity among interspecific hybrids and a decrease in heterozygosity among domesticated descendants of those hybrids. Selective sweeps were found across the genome in early and modern phases of domestication—59–76% of the selectively swept genes originated in the three less dominant ancestral subgenomes. Contrary to the tenet that genetic diversity is limited in cultivated strawberry, we found that the octoploid species harbor massive allelic diversity and that F. × ananassa harbors as much allelic diversity as either wild founder. We identified 41.8 M subgenome-specific DNA variants among resequenced wild and domesticated individuals. Strikingly, 98% of common alleles and 73% of total alleles were shared between wild and domesticated populations. Moreover, genome-wide estimates of nucleotide diversity were virtually identical in F. chiloensis,F. virginiana, and F. × ananassa (π = 0.0059–0.0060). We found, however, that nucleotide diversity and heterozygosity were significantly lower in modern F. × ananassa populations that have experienced significant genetic gains and have produced numerous agriculturally important cultivars.     

A Rapid Method to Estimate Population Variables for African Elephants

Abstract: We developed a noninvasive method to estimate reproductive and survival parameters for free-ranging African savannah elephants (Loxodonta africana africana) and used these to estimate finite population growth rates. We used published data from 2 populations with known growth rates and birth and survival histories to validate our technique. Based on body measurements, our method yielded estimates of age at first and last calving, calving interval, and age-specific survival rates that were similar to those determined during long-term studies at both Addo Elephant National Park and Amboseli National Park. Our technique generated population data required to estimate population growth rates. The method may be particularly useful where censuses yield imprecise estimates or where long-term population data are unavailable. (JOURNAL OF WILDLIFE MANAGEMENT 72(3):822–829; 2008)