Estimating heritability explained by local ancestry and evaluating stratification bias in admixture mapping from summary statistics.

The heritability explained by local ancestry markers in an admixed population (hγ2) provides crucial insight into the genetic architecture of a complex disease or trait. Estimation of hγ2 can be susceptible to biases due to population structure in ancestral populations. Here, we present heritability estimation from admixture mapping summary statistics (HAMSTA), an approach that uses summary statistics from admixture mapping to infer heritability explained by local ancestry while adjusting for biases due to ancestral stratification. Through extensive simulations, we demonstrate that HAMSTA hγ2 estimates are approximately unbiased and are robust to ancestral stratification compared to existing approaches. In the presence of ancestral stratification, we show a HAMSTA-derived sampling scheme provides a calibrated family-wise error rate (FWER) of ∼5% for admixture mapping, unlike existing FWER estimation approaches. We apply HAMSTA to 20 quantitative phenotypes of up to 15,988 self-reported African American individuals in the Population Architecture using Genomics and Epidemiology (PAGE) study. We observe hˆγ2 in the 20 phenotypes range from 0.0025 to 0.033 (mean hˆγ2 = 0.012 ± 9.2 × 10-4), which translates to hˆ2 ranging from 0.062 to 0.85 (mean hˆ2 = 0.30 ± 0.023). Across these phenotypes we find little evidence of inflation due to ancestral population stratification in current admixture mapping studies (mean inflation factor of 0.99 ± 0.001). Overall, HAMSTA provides a fast and powerful approach to estimate genome-wide heritability and evaluate biases in test statistics of admixture mapping studies.

Ancestral Origins and Admixture History of Kazakhs.

Kazakh people, like many other populations that settled in Central Asia, demonstrate an array of mixed anthropological features of East Eurasian (EEA) and West Eurasian (WEA) populations, indicating a possible scenario of biological admixture between already differentiated EEA and WEA populations. However, their complex biological origin, genomic makeup, and genetic interaction with surrounding populations are not well understood. To decipher their genetic structure and population history, we conducted, to our knowledge, the first whole-genome sequencing study of Kazakhs residing in Xinjiang (KZK). We demonstrated that KZK derived their ancestries from 4 ancestral source populations: East Asian (∼39.7%), West Asian (∼28.6%), Siberian (∼23.6%), and South Asian (∼8.1%). The recognizable interactions of EEA and WEA ancestries in Kazakhs were dated back to the 15th century BCE. Kazakhs were genetically distinctive from the Uyghurs in terms of their overall genomic makeup, although the 2 populations were closely related in genetics, and both showed a substantial admixture of western and eastern peoples. Notably, we identified a considerable sex-biased admixture, with an excess of western males and eastern females contributing to the KZK gene pool. We further identified a set of genes that showed remarkable differentiation in KZK from the surrounding populations, including those associated with skin color (SLC24A5, OCA2), essential hypertension (HLA-DQB1), hypertension (MTHFR, SLC35F3), and neuron development (CNTNAP2). These results advance our understanding of the complex history of contacts between Western and Eastern Eurasians, especially those living or along the old Silk Road.

MultiWaverX: modeling latent sex-biased admixture history.

Sex-biased gene flow has been common in the demographic history of modern humans. However, the lack of sophisticated methods for delineating the detailed sex-biased admixture process prevents insights into complex admixture history and thus our understanding of the evolutionary mechanisms of genetic diversity. Here, we present a novel algorithm, MultiWaverX, for modeling complex admixture history with sex-biased gene flow. Systematic simulations showed that MultiWaverX is a powerful tool for modeling complex admixture history and inferring sex-biased gene flow. Application of MultiWaverX to empirical data of 17 typical admixed populations in America, Central Asia, and the Middle East revealed sex-biased admixture histories that were largely consistent with the historical records. Notably, fine-scale admixture process reconstruction enabled us to recognize latent sex-biased gene flow in certain populations that would likely be overlooked by much of the routine analysis with commonly used methods. An outstanding example in the real world is the Kazakh population that experienced complex admixture with sex-biased gene flow but in which the overall signature has been canceled due to biased gene flow from an opposite direction.

Extensive genetic admixture between Tai-Kadai-speaking people and their neighbours in the northeastern region of the Yungui Plateau inferred from genome-wide variations.

BACKGROUND: Yungui Plateau in Southwest China is characterized by multi-language and multi-ethnic communities and is one of the regions with the wealthiest ethnolinguistic, cultural and genetic diversity in East Asia. There are numerous Tai-Kadai (TK)-speaking populations, but their detailed evolutionary history and biological adaptations are still unclear. RESULTS: Here, we genotyped genome-wide SNP data of 77 unrelated TK-speaking Zhuang and Dong individuals from the Yungui Plateau and explored their detailed admixture history and adaptive features using clustering patterns, allele frequency differentiation and sharing haplotype patterns. TK-speaking Zhuang and Dong people in Guizhou are closely related to geographically close TK and Hmong-Mien (HM)-speaking populations. Besides, we identified that Guizhou TK-speaking people have a close genetic relationship with Austronesian (AN)-speaking Atayal and Paiwan people, which is supported by the common origin of the ancient Baiyue tribe. We additionally found subtle genetic differences among the newly studied TK people and previously reported Dais via the fine-scale genetic substructure analysis based on the shared haplotype chunks. Finally, we identified specific selection candidate signatures associated with several essential human immune systems and neurological disorders, which could provide evolutionary evidence for the allele frequency distribution pattern of genetic risk loci. CONCLUSIONS: Our comprehensive genetic characterization of TK people suggested the strong genetic affinity within TK groups and extensive gene flow with geographically close HM and Han people. We also provided genetic evidence that supported the common origin hypothesis of TK and AN people. The best-fitted admixture models further suggested that ancestral sources from northern millet farmers and southern inland and coastal people contributed to the formation of the gene pool of the Zhuang and Dong people.

Native American ancestry and breast cancer risk in Colombian and Mexican women: ruling out potential confounding through ancestry-informative markers.

BACKGROUND: Latin American and Hispanic women are less likely to develop breast cancer (BC) than women of European descent. Observational studies have found an inverse relationship between the individual proportion of Native American ancestry and BC risk. Here, we use ancestry-informative markers to rule out potential confounding of this relationship, estimating the confounder-free effect of Native American ancestry on BC risk. METHODS AND STUDY POPULATION: We used the informativeness for assignment measure to select robust instrumental variables for the individual proportion of Native American ancestry. We then conducted separate Mendelian randomization (MR) analyses based on 1401 Colombian women, most of them from the central Andean regions of Cundinamarca and Huila, and 1366 Mexican women from Mexico City, Monterrey and Veracruz, supplemented by sensitivity and stratified analyses. RESULTS: The proportion of Colombian Native American ancestry showed a putatively causal protective effect on BC risk (inverse variance-weighted odds ratio [OR] = 0.974 per 1% increase in ancestry proportion, 95% confidence interval [CI] 0.970-0.978, p = 3.1 × 10-40). The corresponding OR for Mexican Native American ancestry was 0.988 (95% CI 0.987-0.990, p = 1.4 × 10-44). Stratified analyses revealed a stronger association between Native American ancestry and familial BC (Colombian women: OR = 0.958, 95% CI 0.952-0.964; Mexican women: OR = 0.973, 95% CI 0.969-0.978), and stronger protective effects on oestrogen receptor (ER)-positive BC than on ER-negative and triple-negative BC. CONCLUSIONS: The present results point to an unconfounded protective effect of Native American ancestry on BC risk in both Colombian and Mexican women which appears to be stronger for familial and ER-positive BC. These findings provide a rationale for personalised prevention programmes that take genetic ancestry into account, as well as for future admixture mapping studies.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island.

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

Type 2 diabetes risks and determinants in second-generation migrants and mixed ethnicity people of South Asian and African Caribbean descent in the UK.

AIMS/HYPOTHESIS: Excess risks of type 2 diabetes in UK South Asians (SA) and African Caribbeans (AC) compared with Europeans remain unexplained. We studied risks and determinants of type 2 diabetes in first- and second-generation (born in the UK) migrants, and in those of mixed ethnicity. METHODS: Data from the UK Biobank, a population-based cohort of ~500,000 participants aged 40-69 at recruitment, were used. Type 2 diabetes was assigned using self-report and HbA1c. Ethnicity was both self-reported and genetically assigned using admixture level scores. European, mixed European/South Asian (MixESA), mixed European/African Caribbean (MixEAC), SA and AC groups were analysed, matched for age and sex to enable comparison. In the frames of this cross-sectional study, we compared type 2 diabetes in second- vs first-generation migrants, and mixed ethnicity vs non-mixed groups. Risks and explanations were analysed using logistic regression and mediation analysis, respectively. RESULTS: Type 2 diabetes prevalence was markedly elevated in SA (599/3317 = 18%) and AC (534/4180 = 13%) compared with Europeans (140/3324 = 4%). Prevalence was lower in second- vs first-generation SA (124/1115 = 11% vs 155/1115 = 14%) and AC (163/2200 = 7% vs 227/2200 = 10%). Favourable adiposity (i.e. lower waist/hip ratio or BMI) contributed to lower risk in second-generation migrants. Type 2 diabetes in mixed populations (MixESA: 52/831 = 6%, MixEAC: 70/1045 = 7%) was lower than in comparator ethnic groups (SA: 18%, AC: 13%) and higher than in Europeans (4%). Greater socioeconomic deprivation accounted for 17% and 42% of the excess type 2 diabetes risk in MixESA and MixEAC compared with Europeans, respectively. Replacing self-reported with genetically assigned ethnicity corroborated the mixed ethnicity analysis. CONCLUSIONS/INTERPRETATION: Type 2 diabetes risks in second-generation SA and AC migrants are a fifth lower than in first-generation migrants. Mixed ethnicity risks were markedly lower than SA and AC groups, though remaining higher than in Europeans. Distribution of environmental risk factors, largely obesity and socioeconomic status, appears to play a key role in accounting for ethnic differences in type 2 diabetes risk.

Genetic Origins and Sex-Biased Admixture of the Huis.

The Hui people are unique among Chinese ethnic minorities in that they speak the same language as Han Chinese (HAN) but practice Islam. However, as the second-largest minority group in China numbering well over 10 million, the Huis are under-represented in both global and regional genomic studies. Here, we present the first whole-genome sequencing effort of 234 Hui individuals (NXH) aged over 60 who have been living in Ningxia, where the Huis are mostly concentrated. NXH are genetically more similar to East Asian than to any other global populations. In particular, the genetic differentiation between NXH and HAN (FST = 0.0015) is only slightly larger than that between northern and southern HAN (FST = 0.0010), largely attributed to the western ancestry in NXH (∼10%). Highly differentiated functional variants between NXH and HAN were identified in genes associated with skin pigmentation (e.g., SLC24A5), facial morphology (e.g., EDAR), and lipid metabolism (e.g., ABCG8). The Huis are also distinct from other Muslim groups such as the Uyghurs (FST = 0.0187), especially, NXH derived much less western ancestry (∼10%) compared with the Uyghurs (∼50%). Modeling admixture history indicated that NXH experienced an episode of two-wave admixture. An ancient admixture occurred ∼1,025 years ago, reflecting the intensive west-east contacts during the late Tang Dynasty, and the Five Dynasties and Ten Kingdoms period. A recent admixture occurred ∼500 years ago, corresponding to the Ming Dynasty. Notably, we identified considerable sex-biased admixture, that is, excess of western males and eastern females contributing to the NXH gene pool. The origins and the genomic diversity of the Hui people imply the complex history of contacts between western and eastern Eurasians.

Genome-wide Significance Thresholds for Admixture Mapping Studies.

Admixture mapping studies have become more common in recent years, due in part to technological advances and growing international efforts to increase the diversity of genetic studies. However, many open questions remain about appropriate implementation of admixture mapping studies, including how best to control for multiple testing, particularly in the presence of population structure. In this study, we develop a theoretical framework to characterize the correlation of local ancestry and admixture mapping test statistics in admixed populations with contributions from any number of ancestral populations and arbitrary population structure. Based on this framework, we develop an analytical approach for obtaining genome-wide significance thresholds for admixture mapping studies. We validate our approach via analysis of simulated traits with real genotype data for 8,064 unrelated African American and 3,425 Hispanic/Latina women from the Women's Health Initiative SNP Health Association Resource (WHI SHARe). In an application to these WHI SHARe data, our approach yields genome-wide significant p value thresholds of 2.1 × 10-5 and 4.5 × 10-6 for admixture mapping studies in the African American and Hispanic/Latina cohorts, respectively. Compared to other commonly used multiple testing correction procedures, our method is fast, easy to implement (using our publicly available R package), and controls the family-wise error rate even in structured populations. Importantly, we note that the appropriate admixture mapping significance threshold depends on the number of ancestral populations, generations since admixture, and population structure of the sample; as a result, significance thresholds are not, in general, transferable across studies.

Influence of Pliocene and Pleistocene climates on hybridization patterns between two closely related oak species in China.

BACKGROUND AND AIMS: Contemporary patterns of genetic admixture reflect imprints of both ancient and recent gene flow, which can provide us with valuable information on hybridization history in response to palaeoclimate change. Here, we examine the relationships between present admixture patterns and past climatic niche suitability of two East Asian Cerris oaks (Quercus acutissima and Q. chenii) to test the hypothesis that the mid-Pliocene warm climate promoted while the Pleistocene cool climate limited hybridization among local closely related taxa. METHODS: We analyse genetic variation at seven nuclear microsatellites (1111 individuals) and three chloroplast intergenic spacers (576 individuals) to determine the present admixture pattern and ancient hybridization history. We apply an information-theoretic model selection approach to explore the associations of genetic admixture degree with past climatic niche suitability at multiple spatial scales. KEY RESULTS: More than 70 % of the hybrids determined by Bayesian clustering analysis and more than 90 % of the individuals with locally shared chloroplast haplotypes are concentrated within a mid-Pliocene contact zone between ~30°N and 35°N. Climatic niche suitabilities for Q. chenii during the mid-Pliocene Warm Period [mPWP, ~3.264-3.025 million years ago (mya)] and during the Last Glacial Maximum (LGM, ~0.022 mya) best explain the admixture patterns across all Q. acutissima populations and across those within the ancient contact zone, respectively. CONCLUSIONS: Our results highlight that palaeoclimate change shapes present admixture patterns by influencing the extent of historical range overlap. Specifically, the mid-Pliocene warm climate promoted ancient contact, allowing widespread hybridization throughout central China. In contrast, the Pleistocene cool climate caused the local extinction of Q. chenii, reducing the probability of interspecific gene flow in most areas except those sites having a high level of ecological stability.

Genomic insight into the population history of Central Han Chinese.

BACKGROUND: In recent decades, considerable attention has been paid to exploring the population genetic characteristics of Han Chinese, mainly documenting a north-south genetic substructure. However, the central Han Chinese have been largely underrepresented in previous studies. AIM: To infer a comprehensive understanding of the homogenisation process and population history of Han Chinese. SUBJECTS AND METHODS: We collected samples from 122 Han Chinese from seven counties of Hubei province in central China and genotyped 534,000 genome-wide SNPs. We compared Hubei Han with both ancient and present-day Eurasian populations using Principal Component Analysis, ADMIXTURE, f statistics, qpWave and qpAdm. RESULTS: We observed Hubei Han Chinese are at a genetically intermediate position on the north-south Han Chinese cline. We have not detected any significant genetic substructure in the studied groups from seven different counties. Hubei Han show significant evidence of genetic admixture deriving about 63% of ancestry from Tai-Kadai or Austronesian-speaking southern indigenous groups and 37% from Tungusic or Mongolic related northern populations. CONCLUSIONS: The formation of Han Chinese has involved extensive admixture with Tai-Kadai or Austronesian-speaking populations in the south and Tungusic or Mongolic speaking populations in the north. The convenient transportation and central location of Hubei make it the key region for the homogenisation of Han Chinese.

Inferring the population structure and admixture history of three Hmong-Mien-speaking Miao tribes from southwest China based on genome-wide SNP genotyping.

BACKGROUND: Hmong-Mien speaking Miao, also called Hmong, is the sixth largest ethnic group in mainland China. However, the fine-scale genetic profiles and population history of Miao populations in southwest China, especially in Guizhou province, remain uncharacterised due to a scarcity of samples of genome-wide data from different tribes. AIM: To further investigate the population substructure and admixture history of the Guizhou Miao minority. SUBJECTS AND METHODS: We collected 29 samples from three Miao tribes of Guizhou province in southwest China and genotyped about 700,000 genome-wide SNPs of each sample. We analysed newly generated data in together with published modern/ancient East Asian populations datasets via a series of population genetic methods, including principal component analysis (PCA), ADMIXTURE, Fst, TreeMix, f-statistics, qpWave, and qpAdm. RESULTS: PCA and ADMIXTURE results showed that the three studied Guizhou Miao groups consistently fell on the Hmong-Mien-related genetic cline and were relatively genetically homogeneous, displaying a genetic affinity with neighbouring Tai-Kadai speaking populations such as Dong. These results were further confirmed by the observed genetic clade in Fst, TreeMix, outgroup-f3-statistics, and f4-statistics. Furthermore, f4-based allele sharing patterns illustrated that compared with Hunan Miao in central China, Guizhou Miao shared more alleles with Hmong-Mien-speaking Vietnam Hmong and Tai-Kadai-speaking CoLao, Dong, while exhibiting less northeast Asian-related ancestry. Admixture-f3 and f4-statistics revealed the North-South admixture pattern for the studied Guizhou Miao. A qpAdm-based two-way admixture model further revealed that the studied Guizhou Miao harboured 44%-55.4% indigenous Austronesian-speaking Atayal-related ancestry and 44.6%-56% Late Neolithic Yellow River farmer-related ancestry. CONCLUSIONS: The population structure within Hmong-Mien-related populations showed a geographic correlation. Hmong-Mien speaking Hunan Miao, Guizhou Miao, and Vietnam Hmong presented close genetic relationships although they dwelt in different regions, suggesting the preservation of the original Hmong-related genetic diversity. The results based on genome-wide SNPs data generally matched the migration history for the Miao population. Our study contributes to a better knowledge of Miao populations and the population structure in southwest China.

Revisiting the genetic background and phylogenetic structure of five Sino-Tibetan-speaking populations: insights from autosomal InDels.

The Tibetan-Yi Corridor, located on the eastern edge of the Tibetan Plateau, is the main route of the people of the plateau. Human settlement and diffusion along the corridor have played a pivotal role in shaping the genetic architecture of Sino-Tibetan-speaking (STs) populations in China. In this study, five STs groups (Chengdu Tibetan, Chengdu Han, Muli Tibetan, Lugu Lake Mosuo and Xichang Yi) settling in the Tibetan-Yi Corridor were genotyped via AGCU InDel 50 kit on the capillary electrophoresis platform to decrypt the genetic landscape and phylogenetic relationship of STs populations and investigate the forensic characteristics. Allele frequency distributions of all autosomal insertion/deletion polymorphisms (InDels) in studied groups comply with Hardy-Weinberg equilibrium. The combined power of discrimination values are 0.9999999999999999998, 0.9999999999999999995, 0.9999999999999999999, 0.999999999999999993 and 0.99999999999999999994, respectively, and all the combined probability of exclusion values exceed 0.9990. Forensically relevant statistics implied that these InDels could be used for individual identification and as a promising alternative to STR profiling in paternity testing. Typical population comparisons showed strikingly high homogeneity among studied STs people, indicating complicated genetic admixture among populations in the Tibetan-Yi Corridor. The STs groups in the Tibetan-Yi Corridor keep close genetic affinity with geographically or linguistically close populations, and the genetic components of investigated populations arose from a mixture of multiple ancestral gene pools (resulting from the admixture from the ancestral Highland Tibetans and ancestral Lowland indigenous populations).

Leveraging eDNA to expand the study of hybrid zones.

Hybrid zones are important windows into ecological and evolutionary processes. Our understanding of the significance and prevalence of hybridization in nature has expanded with the generation and analysis of genome-spanning data sets. That said, most hybridization research still has restricted temporal and spatial resolution, which limits our ability to draw broad conclusions about evolutionary and conservation related outcomes. Here, we argue that rapidly advancing environmental DNA (eDNA) methodology could be adopted for studies of hybrid zones to increase temporal sampling (contemporary and historical), refine and geographically expand sampling density, and collect data for taxa that are difficult to directly sample. Genomic data in the environment offer the potential for near real-time biological tracking of hybrid zones, and eDNA provides broad, but as yet untapped, potential to address eco-evolutionary questions.

Vicariance, dispersal, extinction and hybridization underlie the evolutionary history of Atlantic forest fire-eye antbirds (Aves: Thamnophilidae).

In order to gain insights into the biogeographic processes underlying biotic diversification in the Atlantic Forest (AF), we used a multi-locus approach to examine the evolutionary history of the White-shouldered Fire-eye (Pyriglena leucoptera) and the Fringe-backed Fire-eye (Pyriglena atra), two parapatric sister species endemic to the AF. We sequenced one mitochondrial, three Z chromosome-linked and three anonymous markers of 556 individuals from 66 localities. We recovered four lineages throughout the AF: P. atra and three populations within P. leucoptera. All populations diverged during the late Pleistocene and presented varying levels of admixture. One Z-linked locus showed the highest level of differentiation between the two species. On the other hand, a mitochondrial haplotype was shared extensively between them. Our data supported vicariance driving speciation along with extinction and dispersal as processes underlying intraspecific diversification. Furthermore, signatures of demographic expansion in most populations and areas of genetic admixture were recovered throughout the AF, suggesting that forest fragmentation was also important in differentiation. Genetic admixture areas are located between large rivers suggesting that AF rivers may diminish gene flow. Our results indicated a complex and dynamic biogeographic history of Pyriglena in the AF, with vicariance, extinction, dispersal and secondary contact followed by introgression likely influencing the current patterns of genetic distribution.

Genetic structure and admixture in sheep from terminal breeds in the United States.

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright's FST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright's FST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross-validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.

Genetic diversity of avocado from the southern highlands of Tanzania as revealed by microsatellite markers.

BACKGROUND: Avocado is an important cash crop in Tanzania, however its genetic diversity is not thoroughly investigated. This study was undertaken to explore the genetic diversity of avocado in the southern highlands using microsatellite markers. A total of 226 local avocado trees originating from seeds were sampled in eight districts of the Mbeya, Njombe and Songwe regions. Each district was considered as a population. The diversity at 10 microsatellite loci was investigated. RESULTS: A total of 167 alleles were detected across the 10 loci with an average of 16.7 ± 1.3 alleles per locus. The average expected and observed heterozygosity were 0.84 ± 0.02 and 0.65 ± 0.04, respectively. All but two loci showed a significant deviation from the Hardy-Weinberg principle. Analysis of molecular variance showed that about 6% of the variation was partitioned among the eight geographic populations. Population FST pairwise comparisons revealed lack of genetic differentiation for the seven of 28 population pairs tested. The principal components analysis (PCA) and hierarchical cluster analysis showed a mixing of avocado trees from different districts. The model-based STRUCTURE subdivided the trees samples into four major genetic clusters. CONCLUSION: High diversity detected in the analysed avocado germplasm implies that this germplasm is a potentially valuable source of variable alleles that might be harnessed for genetic improvement of this crop in Tanzania. The mixing of avocado trees from different districts observed in the PCA and dendrogram points to strong gene flow among the avocado populations, which led to population admixture revealed in the STRUCTURE analysis. However, there is still significant differentiation among the tree populations from different districts that can be utilized in the avocado breeding program.

Genetic substructure and admixture of Mongolians and Kazakhs inferred from genome-wide array genotyping.

BACKGROUND: Mongolian populations are widely distributed geographically, showing abundant ethnic diversity with geographic and tribal differences. AIM: To infer the genetic substructure, admixture and ancient genetic sources of Mongolians together with Kazakhs. SUBJECTS AND METHODS: We genotyped more than 690,000 genome-wide SNPs from 33 Mongolian and Chinese Kazakh individuals and compared these with both ancient and present-day Eurasian populations using Principal Component Analysis (PCA), ADMIXTURE, Refine-IBD, f statistics, qpWave and qpAdm. RESULTS: We found genetic substructures within Mongolians corresponding to Ölöd, Chahar, and Inner Mongolian clusters, which was consistent with tribe classifications. Mongolian and Kazakh groups derived about 6-40% of West Eurasian related ancestry, most likely from Bronze Age Steppe populations. The East Asian related ancestry in Mongolian and Kazakh groups was well represented by the Neolithic DevilsCave related nomadic lineage, comprising 42-64% of studied groups. We also detected 10-51% of Han Chinese related ancestry in Mongolian and Kazakh groups, especially in Inner Mongolians. The average admixture times for Inner Mongolian, Mongolian_Chahar, Mongolian_Ölöd and Chinese Kazakh were about 1381, 626, 635 and 632 years ago, respectively, with Han and French as the sources. CONCLUSION: The DevilsCave related ancestry was once widespread westwards covering a wide geographical range from Far East Russia to the Mongolia Plateau. The formation of present-day Mongolic and Turkic-speaking populations has also received genetic influence from agricultural expansion.

Purifying selection does not drive signatures of convergent local adaptation of lodgepole pine and interior spruce.

BACKGROUND: Lodgepole pine (Pinus contorta) and interior spruce (Picea glauca, Picea engelmannii, and their hybrids) are distantly related conifer species. Previous studies identified 47 genes containing variants associated with environmental variables in both species, providing evidence of convergent local adaptation. However, if the intensity of purifying selection varies with the environment, clines in nucleotide diversity could evolve through linked (background) selection that would yield allele frequency-environment signatures resembling local adaptation. If similar geographic patterns in the strength of purifying selection occur in these species, this could result in the convergent signatures of local adaptation, especially if the landscape of recombination is conserved. In the present study, we investigated whether spatially/environmentally varying purifying selection could give rise to the convergent signatures of local adaptation that had previously reported. RESULTS: We analyzed 86 lodgepole pine and 50 interior spruce natural populations spanning heterogeneous environments in western Canada where previous analyses had found signatures of convergent local adaptation. We estimated nucleotide diversity and Tajima's D for each gene within each population and calculated the strength of correlations between nucleotide diversity and environmental variables. Overall, these estimates in the genes with previously identified convergent local adaptation signatures had no similar pattern between pine and spruce. Clines in nucleotide diversity along environmental variables were found for interior spruce, but not for lodgepole pine. In spruce, genes with convergent adaption signatures showed a higher strength of correlations than genes without convergent adaption signatures, but there was no such disparity in pine, which suggests the pattern in spruce may have arisen due to a combination of selection and hybridization. CONCLUSIONS: The results rule out purifying/background selection as a driver of convergent local adaption signatures in lodgepole pine and interior spruce.

Single haplotype admixture models using large scale HLA genotype frequencies to reproduce human admixture.

The human leukocyte antigen (HLA) is the most polymorphic region in humans. Anthropologists use HLA to trace populations' migration and evolution. However, recent admixture between populations can mask the ancestral haplotype frequency distribution. We present a statistical method based on high-resolution HLA haplotype frequencies to resolve population admixture using a non-negative matrix factorization formalism and validated using haplotype frequencies from 56 world populations. The result is a minimal set of source components (SCs) decoding roughly 90% of the total variance in the studied admixtures. These SCs agree with the geographical distribution, phylogenies, and recent admixture events of the studied groups. With the growing population of multi-ethnic individuals, or individuals that do not report race/ethnic information, the HLA matching process for stem-cell and solid organ transplants is becoming more challenging. The presented algorithm provides a framework that facilitates the breakdown of highly admixed populations into SCs, which can be used to better match the rapidly growing population of multi-ethnic individuals worldwide.