Quantitating and dating recent gene flow between European and East Asian populations

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Refining models of archaic admixture in Eurasia with ArchaicSeeker 2.0

We developed a method, ArchaicSeeker 2.0, to identify introgressed hominin sequences and model multiple-wave admixture. The new method enabled us to discern two waves of introgression from both Denisovan-like and Neanderthal-like hominins in present-day Eurasian populations and an ancient Siberian individual. We estimated that an early Denisovan-like introgression occurred in Eurasia around 118.8-94.0 thousand years ago (kya). In contrast, we detected only one single episode of Denisovan-like admixture in indigenous peoples eastern to the Wallace-Line. Modeling ancient admixtures suggested an early dispersal of modern humans throughout Asia before the Toba volcanic super-eruption 74 kya, predating the initial peopling of Asia as proposed by the traditional Out-of-Africa model. Survived archaic sequences are involved in various phenotypes including immune and body mass (e.g., ZNF169), cardiovascular and lung function (e.g., HHAT), UV response and carbohydrate metabolism (e.g., HYAL1/HYAL2/HYAL3), while "archaic deserts" are enriched with genes associated with skin development and keratinization.

Population genomics of East Asian ethnic groups

East Asia constitutes one-fifth of the global population and exhibits substantial genetic diversity. However, genetic investigations on populations in this region have been largely under-represented compared with European populations. Nonetheless, the last decade has seen considerable efforts and progress in genome-wide genotyping and whole-genome sequencing of the East-Asian ethnic groups. Here, we review the recent studies in terms of ancestral origin, population relationship, genetic differentiation, and admixture of major East- Asian groups, such as the Chinese, Korean, and Japanese populations. We mainly focus on insights from the whole-genome sequence data and also include the recent progress based on mitochondrial DNA (mtDNA) and Y chromosome data. We further discuss the evolutionary forces driving genetic diversity in East-Asian populations, and provide our perspectives for future directions on population genetics studies, particularly on underrepresented indigenous groups in East Asia.

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.

Ancient Genomes Reveal the Evolutionary History and Origin of Cashmere-Producing Goats in China

Goats are one of the most widespread farmed animals across the world; however, their migration route to East Asia and local evolutionary history remain poorly understood. Here, we sequenced 27 ancient Chinese goat genomes dating from the Late Neolithic period to the Iron Age. We found close genetic affinities between ancient and modern Chinese goats, demonstrating their genetic continuity. We found that Chinese goats originated from the eastern regions around the Fertile Crescent, and we estimated that the ancestors of Chinese goats diverged from this population in the Chalcolithic period. Modern Chinese goats were divided into a northern and a southern group, coinciding with the most prominent climatic division in China, and two genes related to hair follicle development, FGF5 and EDA2R, were highly divergent between these populations. We identified a likely causal de novo deletion near FGF5 in northern Chinese goats that increased to high frequency over time, whereas EDA2R harbored standing variation dating to the Neolithic. Our findings add to our understanding of the genetic composition and local evolutionary process of Chinese goats.

Next generation sequencing of a set of ancestry-informative SNPs: ancestry assignment of three continental populations and estimating ancestry composition for Mongolians

When traditional short tandem repeat profiling fails to provide valuable information to arrest the criminal, forensic ancestry inference of the biological samples left at the crime scene will probably offer investigative leads and facilitate the investigation process of the case. That is why there are consistent efforts in developing panels for ancestry inference in forensic science. Presently, a 30-plex next generation sequencing-based assay was exploited in this study by assembling well-differentiated single nucleotide polymorphisms for ancestry assignment of unknown individuals from three continental populations (African, European and East Asian). And meanwhile, relatively balanced population-specific differentiation values were maintained to avoid the over-estimation or under-estimation of co-ancestry proportions in individuals with admixed ancestry. The principal component analysis and STRUCTURE analysis of reference populations, test populations and the studied Mongolian group indicated that the novel assay was efficient enough to determine the ancestry origin of an unknown individual from the three continental populations. Besides, ancestry membership proportion estimations for the Mongolian group revealed that a large fraction of the ancestry was contributed by East Asian genetic component (approximately 83.9%), followed by European (approximately 12.6%) and African genetic components (approximately 3.5%), respectively. And next generation sequencing technology applied in this study offers possibility to incorporate more single nucleotide polymorphisms for individual identification and phenotype prediction into the same assay to provide as many as possible investigative clues in the future.

Heterogeneity in the distribution of 159 drug-response related SNPs in world populations and their genetic relatedness

Interethnic variability in drug response arises from genetic differences associated with drug metabolism, action and transport. These genetic variations can affect drug efficacy as well as cause adverse drug reactions (ADRs). We retrieved drug-response related single nucleotide polymorphism (SNP) associated data from databases and analyzed to elucidate population specific distribution of 159 drug-response related SNPs in twenty six populations belonging to five super-populations (African, Admixed Americans, East Asian, European and South Asian). Significant interpopulation differences exist in the minor (variant) allele frequencies (MAFs), linkage disequilibrium (LD) and haplotype distributions among these populations. 65 of the drug-response related alleles, which are considered as minor (variant) in global population, are present as the major alleles (frequency ≥0.5) in at least one or more populations. Populations that belong to the same super-population have similar distribution pattern for majority of the variant alleles. These drug response related variant allele frequencies and their pairwise LD measure (r²) can clearly distinguish the populations in a way that correspond to the known evolutionary history of human and current geographic distributions, while D' cannot. The data presented here may aid in identifying drugs that are more appropriate and/or require pharmacogenetic testing in these populations. Our findings emphasize on the importance of distinct, ethnicity-specific clinical guidelines, especially for the African populations, to avoid ADRs and ensure effective drug treatment.

Worldwide genetic variation of the IGHV and TRBV immune receptor gene families in humans

This article presents a comprehensive study of the IGHV and TRBV gene families in a globally diverse sample of humans and shows that the two gene families exhibit starkly different patterns of variation.

The immunoglobulin heavy variable (IGHV) and T cell beta variable (TRBV) loci are among the most complex and variable regions in the human genome. Generated through a process of gene duplication/deletion and diversification, these loci can vary extensively between individuals in copy number and contain genes that are highly similar, making their analysis technically challenging. Here, we present a comprehensive study of the functional gene segments in the IGHV and TRBV loci, quantifying their copy number and single-nucleotide variation in a globally diverse sample of 109 (IGHV) and 286 (TRBV) humans from over a 100 populations. We find that the IGHV and TRBV gene families exhibit starkly different patterns of variation. In addition to providing insight into the different evolutionary paths of the IGHV and TRBV loci, our results are also important to the adaptive immune repertoire sequencing community, where the lack of frequencies of common alleles and copy number variants is hampering existing analytical pipelines.

Whole-genome sequencing of 175 Mongolians uncovers population-specific genetic architecture and gene flow throughout North and East Asia

The genetic variation in Northern Asian populations is currently undersampled. To address this, we generated a new genetic variation reference panel by whole-genome sequencing of 175 ethnic Mongolians, representing six tribes. The cataloged variation in the panel shows strong population stratification among these tribes, which correlates with the diverse demographic histories in the region. Incorporating our results with the 1000 Genomes Project panel identifies derived alleles shared between Finns and Mongolians/Siberians, suggesting that substantial gene flow between northern Eurasian populations has occurred in the past. Furthermore, we highlight that North, East, and Southeast Asian populations are more aligned with each other than these groups are with South Asian and Oceanian populations.

Deep-coverage whole genome sequences and blood lipids among 16,324 individuals

Large-scale deep-coverage whole-genome sequencing (WGS) is now feasible and offers potential advantages for locus discovery. We perform WGS in 16,324 participants from four ancestries at mean depth >29X and analyze genotypes with four quantitative traits-plasma total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol, and triglycerides. Common variant association yields known loci except for few variants previously poorly imputed. Rare coding variant association yields known Mendelian dyslipidemia genes but rare non-coding variant association detects no signals. A high 2M-SNP LDL-C polygenic score (top 5th percentile) confers similar effect size to a monogenic mutation (~30 mg/dl higher for each); however, among those with severe hypercholesterolemia, 23% have a high polygenic score and only 2% carry a monogenic mutation. At these sample sizes and for these phenotypes, the incremental value of WGS for discovery is limited but WGS permits simultaneous assessment of monogenic and polygenic models to severe hypercholesterolemia.

40,000-Year-Old Individual from Asia Provides Insight into Early Population Structure in Eurasia

By at least 45,000 years before present, anatomically modern humans had spread across Eurasia [1-3], but it is not well known how diverse these early populations were and whether they contributed substantially to later people or represent early modern human expansions into Eurasia that left no surviving descendants today. Analyses of genome-wide data from several ancient individuals from Western Eurasia and Siberia have shown that some of these individuals have relationships to present-day Europeans [4, 5] while others did not contribute to present-day Eurasian populations [3, 6]. As contributions from Upper Paleolithic populations in Eastern Eurasia to present-day humans and their relationship to other early Eurasians is not clear, we generated genome-wide data from a 40,000-year-old individual from Tianyuan Cave, China, [1, 7] to study his relationship to ancient and present-day humans. We find that he is more related to present-day and ancient Asians than he is to Europeans, but he shares more alleles with a 35,000-year-old European individual than he shares with other ancient Europeans, indicating that the separation between early Europeans and early Asians was not a single population split. We also find that the Tianyuan individual shares more alleles with some Native American groups in South America than with Native Americans elsewhere, providing further support for population substructure in Asia [8] and suggesting that this persisted from 40,000 years ago until the colonization of the Americas. Our study of the Tianyuan individual highlights the complex migration and subdivision of early human populations in Eurasia.

Modeling Continuous Admixture Using Admixture-Induced Linkage Disequilibrium

Recent migrations and inter-ethnic mating of long isolated populations have resulted in genetically admixed populations. To understand the complex population admixture process, which is critical to both evolutionary and medical studies, here we used admixture induced linkage disequilibrium (LD) to infer continuous admixture events, which is common for most existing admixed populations. Unlike previous studies, we expanded the typical continuous admixture model to a more general scenario with isolation after a certain duration of continuous gene flow. Based on the new models, we developed a method, CAMer, to infer the admixture history considering continuous and complex demographic process of gene flow between populations. We evaluated the performance of CAMer by computer simulation and further applied our method to real data analysis of a few well-known admixed populations.

Inference of multiple-wave population admixture by modeling decay of linkage disequilibrium with polynomial functions

To infer the histories of population admixture, one important challenge with methods based on the admixture linkage disequilibrium (ALD) is to remove the effect of source LD (SLD), which is directly inherited from source populations. In previous methods, only the decay curve of weighted LD between pairs of sites whose genetic distance were larger than a certain starting distance was fitted by single or multiple exponential functions, for the inference of recent single- or multiple-wave admixture. However, the effect of SLD has not been well defined and no tool has been developed to estimate the effect of SLD on weighted LD decay. In this study, we defined the SLD in the formularized weighted LD statistic under the two-way admixture model and proposed a polynomial spectrum (p-spectrum) to study the weighted SLD and weighted LD. We also found that reference populations could be used to reduce the SLD in weighted LD statistics. We further developed a method, iMAAPs, to infer multiple-wave admixture by fitting ALD using a p-spectrum. We evaluated the performance of iMAAPs under various admixture models in simulated data and applied iMAAPs to the analysis of genome-wide single nucleotide polymorphism data from the Human Genome Diversity Project and the HapMap Project. We showed that iMAAPs is a considerable improvement over other current methods and further facilitates the inference of histories of complex population admixtures.

East Eurasian ancestry in the middle of Europe: genetic footprints of Steppe nomads in the genomes of Belarusian Lipka Tatars

Medieval era encounters of nomadic groups of the Eurasian Steppe and largely sedentary East Europeans had a variety of demographic and cultural consequences. Amongst these outcomes was the emergence of the Lipka Tatars—a Slavic-speaking Sunni-Muslim minority residing in modern Belarus, Lithuania and Poland, whose ancestors arrived in these territories via several migration waves, mainly from the Golden Horde. Our results show that Belarusian Lipka Tatars share a substantial part of their gene pool with Europeans as indicated by their Y-chromosomal, mitochondrial and autosomal DNA variation. Nevertheless, Belarusian Lipkas still retain a strong genetic signal of their nomadic ancestry, witnessed by the presence of common Y-chromosomal and mitochondrial DNA variants as well as autosomal segments identical by descent between Lipkas and East Eurasians from temperate and northern regions. Hence, we document Lipka Tatars as a unique example of former Medieval migrants into Central Europe, who became sedentary, changed language to Slavic, yet preserved their faith and retained, both uni- and bi-parentally, a clear genetic echo of a complex population interplay throughout the Eurasian Steppe Belt, extending from Central Europe to northern China.

The Prognostic Value of Decreased LKB1 in Solid Tumors: A Meta-Analysis

Background

Liver kinase B1 (LKB1) is a protein kinase that regulates the growth, integrity and polarity of mammalian cells. Recent studies have reported the prognostic value of decreased LKB1 expression in different tumors. However, the results of these studies remain controversial. Therefore, this meta-analysis was performed to more accurately estimate the role of decreased LKB1 in the prognostication of human solid tumors.

Methods

A systematic literature search in the electronic databases PubMed, Embase, Web of Science and CNKI (updated to October 15, 2015) was performed to identify eligible studies. The overall survival (OS), relapse-free survival (RFS), disease-free survival (DFS) and clinicopathological features data were collected from these studies. The hazard ratios (HRs), odds ratios (ORs) and 95% confidence intervals (CIs) were calculated and pooled with a random-effects models using Stata12.0 software.

Results

A total of 14 studies covering 1915 patients with solid tumors were included in this meta-analysis. Decreased LKB1 was associated with poorer OS in both the univariate (HR: 1.86, 95%CI: 1.42–2.42, P<0.001) and multivariate (HR: 1.55, 95%CI: 1.09–2.21, P = 0.015) analyses. A subgroup analysis revealed that the associations between decreased LKB1 and poor OS were significant within the Asian region (HR 2.18, 95%CI: 1.66–2.86, P<0.001) and obvious for lung cancer (HR: 2.16, 95%CI: 1.47–3.18, P<0.001). However, the articles that involved analyses of both RFS and DFS numbered only 3, and no statistically significant correlations of decreased LKB1 with RFS or DFS were observed in this study. Additionally, the pooled odds ratios (ORs) indicated that decreased LKB1 was associated with larger tumor size (OR: 1.60, 95%CI: 1.09–2.36, P = 0.017), lymph node metastasis (OR: 2.41, 95%CI: 1.53–3.78, P<0.001) and a higher TNM stage (OR: 3.35, 95%CI: 2.20–5.09, P<0.001).

Conclusion

These results suggest that decreased LKB1 expression in patients with solid tumors might be related to poor prognosis and serve as a potential predictive marker of poor clinicopathological prognostic factors. Additional studies are required to verify the clinical utility of decreased LKB1 in solid tumors.

Premature termination codons in modern human genomes

The considerable range of genetic variation in human populations may partly reflect distinctive processes of adaptation to variable environmental conditions. However, the adaptive genomic signatures remain to be completely elucidated. This research explores candidate loci under selection at the population level by characterizing recently arisen premature termination codons (PTCs), some of which indicate a human knockout. From a total of 7595 participants from two population exome projects, 246 PTCs were found where natural selection has resulted in new alleles with a high frequency (from 1% to 96%) of derived alleles and various levels of population differentiation (F_ST = 0.00139–0.626). The PTC genes formed protein and regulatory networks limited to 15 biological processes or gene families, of which seven categories were previously unreported. PTC mutations have a strong tendency to be introduced into members of the same gene family, even during modern human evolution, although the exact nature of the selection is not fully known. The findings here suggest the ongoing evolutionary plasticity of modern humans at the genetic level and also partly provide insights into common human knockouts.

DXS998-DXS548-FRAXAC1 represents a novel informative haplotype at the FMR1 locus in the Iranian population

Fragile X syndrome, which is caused by mutation in the FMR1 gene region, is one of the most prevalent forms of mental retardation. Direct diagnosis of the disease is based on PCR and southern blot analysis, but because of technical problems, use of polymorphic DNA markers can be helpful for carrier detection and prenatal diagnosis in families with an affected individual. The polymorphic markers usually show a population-based haplotype frequency and heterozygosity. In the present study, genotyping and analysis of haplotype frequency of three microsatellite markers including DXS998, DXS548 and FRAXAC1 at the FMR1 gene region were carried out in 140 unrelated healthy women and 26 families from the Iranian population. The data indicated the presence of a novel allele for DXS998 in the Iranian population. Estimation of haplotype frequency using Arlequin program showed 50 different DXS998-DXS548-FRAXAC1 haplotypes for the input data of 5, 7 and 4 alleles, respectively. Among these haplotypes five of them showed relatively high frequencies (≥0.05). Analysis of linkage disequilibrium (LD) for the unrelated individuals using the PowerMarker computer program, showed that this haplotype combination can be an informative haplotype for linkage analysis in carrier detection and possible molecular diagnosis of fragile X in the Iranian population.