Positive Selection in the Chromosome 16 VKORC1 Genomic Region Has Contributed to the Variability of Anticoagulant Response in Humans

VKORC1 (vitamin K epoxide reductase complex subunit 1, 16p11.2) is the main genetic determinant of human response to oral anticoagulants of antivitamin K type (AVK). This gene was recently suggested to be a putative target of positive selection in East Asian populations. In this study, we genotyped the HGDP-CEPH Panel for six VKORC1 SNPs and downloaded chromosome 16 genotypes from the HGDP-CEPH database in order to characterize the geographic distribution of footprints of positive selection within and around this locus. A unique VKORC1 haplotype carrying the promoter mutation associated with AVK sensitivity showed especially high frequencies in all the 17 HGDP-CEPH East Asian population samples. VKORC1 and 24 neighboring genes were found to lie in a 505 kb region of strong linkage disequilibrium in these populations. Patterns of allele frequency differentiation and haplotype structure suggest that this genomic region has been submitted to a near complete selective sweep in all East Asian populations and only in this geographic area. The most extreme scores of the different selection tests are found within a smaller 45 kb region that contains VKORC1 and three other genes (BCKDK, MYST1 (KAT8), and PRSS8) with different functions. Because of the strong linkage disequilibrium, it is not possible to determine if VKORC1 or one of the three other genes is the target of this strong positive selection that could explain present-day differences among human populations in AVK dose requirement. Our results show that the extended region surrounding a presumable single target of positive selection should be analyzed for genetic variation in a wide range of genetically diverse populations in order to account for other neighboring and confounding selective events and the hitchhiking effect.     

Dissecting the genetic structure of Korean population using genome-wide SNP arrays

Genome-wide SNP arrays have generated unprecedented quantities of data allow the detection of human evolutionary history and dense genome-wide data also enable the identification of distance ancestry among individuals or ethnic groups. To explain wider aspects of the genetic structure of Koreans and the East Asian population, we analyzed 79 individuals from the Korean HapMap project at 555,352 common single-nucleotide polymorphism loci, and compared this data with the worldwide population groups with the 53 ethnic groups from Human Genome Diversity Panel (HGDP-CEPH). Population differentiation (F_ST), Principal Component Analyses, STRUCTURE and ADMIXTURE are examined. In general, all the individual samples studies here were classified into subset of ethnic groups according to their geographical origins. Korean HapMap individuals were grouped together with East Asian populations from HGDP panel. Recently, a sub-population structure within Korean population has been reported. Our result, however, revealed the genetic homogeneity of Korean population. The ADMIXTURE analysis showed that, overall the Korean populations derive 79 % of their genomic ancestry from southern Asia and have relatively little northern Asian ancestry (21 %). The present work, therefore, provide the evidence that the male-biased southern-to-northern migration influenced not only for the genetic make up of the Y chromosome in the Korean population but also, its autosomal composition.     

Distinctive KIR and HLA diversity in a panel of north Indian Hindus

HLA and KIR are diverse and rapidly evolving gene complexes that work together in human immunity mediated by cytolytic lymphocytes. Understanding their complex immunogenetic interaction requires study of both HLA and KIR diversity in the same human population. Here a panel of 72 unrelated north Indian Hindus was analyzed. HLA- A, B, C, DRB1, DQA1, and DQB1 alleles and their frequencies were determined by sequencing or high-resolution typing of genomic DNA; KIR genotypes were determined by gene-specific typing and by allele-level DNA typing for KIR2DL1, 2DL3, 2DL5, 3DL1, and 3DL2. From HLA analysis, the north Indian population is seen to have several characteristics shared either with Caucasian or East Asian populations, consistent with the demographic history of north India, as well as specific features, including several alleles at high frequency that are rare or absent in other populations. A majority of the north Indian KIR gene profiles have not been seen in Caucasian and Asian populations. Most striking is a higher frequency of the B group of KIR haplotypes, resulting in equal frequencies for A and B group haplotypes in north Indians. All 72 members of the north Indian panel have different HLA genotype and different KIR genotype.     

Diversity of killer cell immunoglobulin-like receptor genes in four ethnic groups in China

Killer cell immunoglobulin-like receptors (KIRs) show extensive variation in terms of gene content and allelic polymorphisms among different populations. The aim of this study was to analyze the distribution of KIR genes in the Bulang, Nu, Yugu, and Zhuang ethnic groups, which belong to four different language families in China, and thus to provide basic KIR gene and genotype data for these Chinese ethnic groups. Genotyping of 16 KIR genes was performed in 425 unrelated individuals using the polymerase chain reaction-sequence-specific oligonucleotide probe method with the Luminex MultiAnalyte Profiling System. The four framework KIR genes were detected in all four ethnic groups. The activating KIR genes as well as the inhibitory KIR genes showed extreme diversity among these four populations. A total of 35 distinct KIR genotypes were identified, one of which was previously unknown. The four most common genotypes were identified in all four populations and comprised 66.1~91.1% of all the genotypes. The group A haplotype occurred more frequently than the group B haplotype in the Nu, Yugu, and Zhuang populations, as in other East Asian populations. In contrast, the group A and group B haplotypes occurred equally in the Bulang population. The results of the present study suggested that the KIR genes and genotypes are diverse in these four ethnic groups, and each ethnic group has its own characteristic KIR distribution. The findings with respect to KIR gene diversity in these four populations should provide relevant genomic diversity data for the future study of viral infections, autoimmune diseases, and reproductive fitness.     

Gene admixture in the silk road region of China: evidence from mtDNA and melanocortin 1 receptor polymorphism

Mitochondrial DNA control region segment I sequences and melanocortin 1 receptor (MC1R) gene polymorphism were examined in ethnic populations in the silk road region of China. Both the frequencies of the MC1R variants and the results of mtDNA data in this region presented intermediate values between those of Europe and East and Southeast Asia, which suggested extensive gene admixture in this area and was in general agreement with previous studies. Phylogenetic analysis of the ethnic populations in the Silk Road region that based on mtDNA data didn't show expected cluster pattern according to their ethnogenesis. We suspect that a high migration rate in female among these closely related populations and other three demographic events might account for it.     

Features of evolution and expansion of modern humans,inferred from genomewide microsatellite markers

We study data on variation in 52 worldwide populations at 377 autosomal short tandem repeat loci, to infer a demographic history of human populations. Variation at di-, tri-, and tetranucleotide repeat loci is distributed differently, although each class of markers exhibits a decrease of within-population genetic variation in the following order: sub-Saharan Africa, Eurasia, East Asia, Oceania, and America. There is a similar decrease in the frequency of private alleles. With multidimensional scaling, populations belonging to the same major geographic region cluster together, and some regions permit a finer resolution of populations. When a stepwise mutation model is used, a population tree based on TD estimates of divergence time suggests that the branches leading to the present sub-Saharan African populations of hunter-gatherers were the first to diverge from a common ancestral population (approximately 71-142 thousand years ago). The branches corresponding to sub-Saharan African farming populations and those that left Africa diverge next, with subsequent splits of branches for Eurasia, Oceania, East Asia, and America. African hunter-gatherer populations and populations of Oceania and America exhibit no statistically significant signature of growth. The features of population subdivision and growth are discussed in the context of the ancient expansion of modern humans.     

Phylogeography and genetic diversity of the red squirrel (Sciurus vulgaris) in China: Implications for the species' postglacial expansion history

The geographical distribution of Sciurus vulgaris spans much of the Palearctic from western Europe and the UK eastward to the pacific coast of East Asia. S. vulgaris occurs in China only in the far northwest and the northeast. Understanding of the species’ postglacial expansion history in East Asia has been limited by the paucity of molecular data. In this study, we used partial D-loop and cytochrome b gene sequences to assess mitochondrial DNA variation in S. vulgaris in China. Our objectives were to (1) determine phylogeographical patterns of S. vulgaris in China; (2) understand the species’ postglacial expansion history in this region; and (3) quantify genetic diversity levels within S. vulgaris populations in China. We identified a supported phylogenetic group in S. vulgaris from China, and found no tendency for haplotypes to cluster by geographic region. Our analysis of S. vulgaris from China and other regions supports the hypothesis that the Calabria region of southern Italy is a glacial refugium for the species. We tentatively propose a postglacial expansion pattern for the squirrels: migrating from Calabria via Central and Eastern Europe to Russia and from there to China, and firstly to the northwest and then to northeast in China. We found high levels of genetic diversity in S. vulgaris populations across China as a whole, and discussed its influential factors.     

Straightforward inference of ancestry and admixture proportions through ancestry-informative insertion deletion multiplexing

Ancestry-informative markers (AIMs) show high allele frequency divergence between different ancestral or geographically distant populations. These genetic markers are especially useful in inferring the likely ancestral origin of an individual or estimating the apportionment of ancestry components in admixed individuals or populations. The study of AIMs is of great interest in clinical genetics research, particularly to detect and correct for population substructure effects in case-control association studies, but also in population and forensic genetics studies. This work presents a set of 46 ancestry-informative insertion deletion polymorphisms selected to efficiently measure population admixture proportions of four different origins (African, European, East Asian and Native American). All markers are analyzed in short fragments (under 230 basepairs) through a single PCR followed by capillary electrophoresis (CE) allowing a very simple one tube PCR-to-CE approach. HGDP-CEPH diversity panel samples from the four groups, together with Oceanians, were genotyped to evaluate the efficiency of the assay in clustering populations from different continental origins and to establish reference databases. In addition, other populations from diverse geographic origins were tested using the HGDP-CEPH samples as reference data. The results revealed that the AIM-INDEL set developed is highly efficient at inferring the ancestry of individuals and provides good estimates of ancestry proportions at the population level. In conclusion, we have optimized the multiplexed genotyping of 46 AIM-INDELs in a simple and informative assay, enabling a more straightforward alternative to the commonly available AIM-SNP typing methods dependent on complex, multi-step protocols or implementation of large-scale genotyping technologies.     

Molecular identification, phylogenetic status, and geographic distribution of Culicoides oxystoma (Diptera: Ceratopogonidae) in Israel

Culicoides oxystoma (Diptera: Ceratopogonidae) is an important vector species, reported mainly from Asia, with high potential to transmit viral diseases affecting livestock. In Japan, many arboviruses have been isolated from C. oxystoma, suggesting it as a key player in the epidemiology of several Culicoides-borne diseases. Over the years, C. oxystoma has also been reported in the Middle East region, including Israel. In this region, however, C. oxystoma cannot be easily distinguished morphologically from its sibling species included in the Culicoides schultzei complex. We therefore used genomic data for species identification and phylogeny resolution. Phylogenetic analyses based on internal transcribed spacer 1 (ITS-1) of ribosomal DNA and the mitochondrial gene encoding cytochrome oxidase subunit I (COI) showed that C. oxystoma from Israel is closely related to C. oxystoma from Japan. Using differential probing PCR, we showed that C. oxystoma is distributed all over the country, especially in Mediterranean climate regions. Culicoides oxystoma is less common or even absent in arid regions, while the other genetic cluster of C. schultzei complex was found only in the east of the country (mostly arid and semiarid regions). The molecular finding of C. oxystoma in wide geographical regions, together with its high proportion in the general Culicoides population and its vectoring potential, imply that it may be an important vector species in the Middle East.     

Phylogeography and Demographic History of Chinese Black-Spotted Frog Populations (<Emphasis Type="Italic">Pelophylax nigromaculata</Emphasis>): Evidence for Independent Refugia Expansion and Secondary Contact

Background  

Pleistocene glaciations had considerable impact on phylogeographic patterns within and among closely related species of many vertebrates. Compared to Europe and North America, research on the phylogeography of vertebrates in East Asia, particularly in China, remains limited. The black-spotted frog (Pelophylax nigromaculata) is a widespread species in East Asia. The wide distribution of this species in China makes it an ideal model for the study of palaeoclimatic effects on vertebrates in East Asia. Our previous studies of P. nigromaculata revealed significant subdivisions between the northeast China populations and populations in other regions of the mainland. In the present study, we aim to see whether the deepest splits among lineages and perhaps subsequent genealogical divisions are temporally consistent with a Pleistocene origin and whether clade geographic distributions, with insight into expansion patterns, are similarly spatially consistent with this model.     

Diversity of killer cell immunoglobulin-like receptor genes in Southern Turkey

Killer cell immunoglobulin-like receptors (KIRs) are a family of inhibitory and activating receptors expressed by natural killer (NK) cells and regulate NK cells’ activity. KIR genes are highly polymorphic markers, characterized by a wide diversity, and can therefore be considered as good population genetic markers. The aim of this study was to determine KIR gene frequencies, ratios of haplotypes and genotypes in Southern Turkey and also to compare the data with other worldwide populations studied previously. The study group consisted of 200 non-related individuals from Southern Turkey. The percentage of each KIR gene in the population group was determined by direct counting. Differences between populations in the distribution of each KIR gene and genotype profile were estimated by two-tailed Fisher Exact test. The most frequent non-framework KIR genes detected in Southern Turkey population were: KIR 2DL1 (97%), KIR 3DL1 (91%), KIR 2DS4 (92%) and the pseudogene 2DP1 (96%). Fourty different genotypes were found in 200 subjects and AA1 genotype was the most frequent (27%). Among 40 different genotypes, ten of these were described for the first time in this study and were added to the database () numerized as genotype ID from 400 to 409. Gene frequencies and found genotypes demonstrated similarity of Southern Turkey’s KIR repertoire with the KIR repertoires of Middle East and European population. High variability seen in KIR genome in this region is thought to be formed as a result of migration and settlement of different civilizations in this region and heterogenity formed in time.     

Genetic structure of Hmong-Mien speaking populations in East Asia as revealed by mtDNA lineages

Hmong-Mien (H-M) is a major language family in East Asia, and its speakers distribute primarily in southern China and Southeast Asia. To date, genetic studies on H-M speaking populations are virtually absent in the literature. In this report, we present the results of an analysis of genetic variations in the mitochondrial DNA (mtDNA) hypervariable segment 1 (HVS1) region and diagnostic variants in the coding regions in 537 individuals sampled from 17 H-M populations across East Asia. The analysis showed that the haplogroups that are predominant in southern East Asia, including B, R9, N9a, and M7, account for 63% (ranging from 45% to 90%) of mtDNAs in H-M populations. Furthermore, analysis of molecular variance (AMOVA), phylogenetic tree analysis, and principal component (PC) analysis demonstrate closer relatedness between H-M and other southern East Asians, suggesting a general southern origin of maternal lineages in the H-M populations. The estimated ages of the mtDNA lineages that are specific to H-M coincide with those based on archeological cultures that have been associated with H-M. Analysis of genetic distance and phylogenetic tree indicated some extent of difference between the Hmong and the Mien populations. Together with the higher frequency of north-dominating lineages observed in the Hmong people, our results indicate that the Hmong populations had experienced more contact with the northern East Asians, a finding consistent with historical evidence. Moreover, our data defined some new (sub-)haplogroups (A6, B4e, B4f, C5, F1a1, F1a1a, and R9c), which will direct further efforts to improve the phylogeny of East Asian mtDNAs.     

Gender and population history: sex bias revealed by studying genetic admixture of Ngazidja population (Comoro Archipelago)

The peopling of Comoro Archipelago is defined by successive waves of migration from three main areas: the East African Coast (Bantu-speaking populations), the Persia and Arabian Peninsula, and Southeast Asia (especially Indonesia). It follows an apparent classic trihybrid admixture model. To better understand the Comorian population admixture dynamics, we analyzed the contributions of these three historical parental components to its genetic pool. To enhance accuracy and reliability, we used both classical and molecular markers. Samples consist of published data: blood group frequencies, 14 KIR genes, 19 mitochondrial DNA SNPs (to highlight female migrations), 14 Y chromosome SNPs (male migrations). We revealed distinct admixture patterns for autosomal and uniparental markers. KIR gene frequencies had never been used to estimate admixture rates, this being a first assessment of their informative power in admixture studies. To avoid major methodological and statistical bias, we determined admixture coefficients through nine well-tried estimators and their associated software programs (ADMIX95, ADMIX, admix 2.0, LEA, LEADMIX, and Mistura). Results from mtDNA and Y chromosome markers point to an important sex-bias in the admixture event. The original Bantu gene pool received a predominant male-mediated contribution from the Arabian Peninsula and Persia, and a female-mediated contribution from Southeast Asia. Admixture rates estimated from autosomal KIR gene markers point also to an unexpected elevated Austronesian contribution.     

Genetic studies of human diversity in East Asia

East Asia is one of the most important regions for studying evolution and genetic diversity of human populations. Recognizing the relevance of characterizing the genetic diversity and structure of East Asian populations for understanding their genetic history and designing and interpreting genetic studies of human diseases, in recent years researchers in China have made substantial efforts to collect samples and generate data especially for markers on Y chromosomes and mtDNA. The hallmark of these efforts is the discovery and confirmation of consistent distinction between northern and southern East Asian populations at genetic markers across the genome. With the confirmation of an African origin for East Asian populations and the observation of a dominating impact of the gene flow entering East Asia from the south in early human settlement, interpretation of the north-south division in this context poses the challenge to the field. Other areas of interest that have been studied include the gene flow between East Asia and its neighbouring regions (i.e. Central Asia, the Sub-continent, America and the Pacific Islands), the origin of Sino-Tibetan populations and expansion of the Chinese.     

KIR genotypic diversity can track ancestries in heterogeneous populations: a potential confounder for disease association studies

Killer cell immunoglobulin-like receptors (KIR) are encoded by highly polymorphic genes that regulate the activation of natural killer (NK) cells and other lymphocyte subsets and likely play key roles in innate and adaptive immunity. Association studies increasingly implicate KIR in disease predisposition and outcome but could be confounded by unknown KIR genetic structure in heterogeneous populations. To examine this, we characterized the diversity of 16 KIR genes in 712 Northern Californians (NC) stratified by self-assigned ethnicities and compared the profiles of KIR polymorphism with other US and global populations using a reference database. Sixty-eight distinct KIR genotypes were characterized: 58 in 457 Caucasians (NCC), 17 in 47 African Americans (NCAA), 21 in 80 Asians (NCA), 20 in 74 Hispanics (NCH), and 18 in 54 "other" ethnicities (NCO). KIR genotype patterns and frequencies in the 4 defined ethnicities were compared with each other and with 34 global populations by phylogenetic analysis. Although there were no population-specific genotypes, the KIR genotype frequency patterns faithfully traced the ancestry of NCC, NCAA, and NCA but not of NCH whose ancestries are known to be more heterogeneous. KIR genotype frequencies can therefore track ethnic ancestries in modern urban populations. Our data emphasize the importance of selecting ethnically matched controls in KIR-based studies to avert spurious associations.     

Inferences on population history of a seed chalcid wasp: invasion success despite a severe founder effect from an unexpected source population

Most invasive species established in Europe originate from either Asia or North America, but little is currently known about the potential of the Anatolian Peninsula (Asia Minor) and/or the Near East to constitute invasion sources. Mediterranean forests are generally fragile ecosystems that can be threatened by invasive organisms coming from different regions of the Mediterranean Basin, but for which historical data are difficult to gather and the phylogeographic patterns are still poorly understood for most terrestrial organisms. In this study, we characterized the genetic structure of Megastigmus schimitscheki, an invasive seed‐feeding insect species originating from the Near East, and elucidated its invasion route in South‐eastern France in the mid 1990s. To disentangle the evolutionary history of this introduction, we gathered samples from the main native regions (Taurus Mountains in Turkey, Lebanon and Cyprus) and from the invaded region that we genotyped using five microsatellite markers and for which we sequenced the mitochondrial Cytochrome Oxidase I gene. We applied a set of population genetic statistics and methods, including approximate Bayesian computation. We proposed a detailed phylogeographic pattern for the Near East populations, and we unambiguously showed that the French invasive populations originated from Cyprus, although the available historical data strongly suggested that Turkey could be the most plausible source area. Interestingly, we could show that the introduced populations were founded from an extremely restricted number of individuals that realized a host switch from Cedrus brevifolia to C. atlantica. Evolutionary hypotheses are discussed to account for this unlikely scenario.     

Linkage disequilibrium organization of the human KIR superlocus: implications for KIR data analyses

An extensive family-based study of linkage disequilibrium (LD) in the killer cell immunoglobulin-like receptors (KIR) cluster was performed. We aimed to describe the LD structure in the KIR gene cluster using a sample of 418 founder haplotypes identified by segregation in a group of 106 families from Northern Ireland. The LD was studied at two levels of polymorphism: the structural level (presence or absence of KIR genes) and the allelic level (between alleles of KIR genes). LD was further assessed using the predictive value of one KIR polymorphism for another one in order to provide an interpretative framework for the LD effect in association studies. In line with previous research, distinct patterns of KIR genetic diversity within the genomic region centromeric to KIR2DL4 (excluding KIR2DL4) and within the telomeric region including KIR2DL4 were found. In a comprehensive PPV/NPV-based LD analysis within the KIR cluster, robust tag markers were found that can be used to identify which genes are concomitantly present or absent and to further identify groups of associated KIR alleles. Several extended KIR haplotypes in the study population were identified (KIR2DS2*POS-KIR2DL2*001-KIR2DL5B*002-KIR2DS3*00103-KIR2DL1*00401; KIR2DL4*011-KIR3DL1/S1*005-KIR2DS4*003-KIR3DL2*003; KIR2DL4*00802-KIR3DL1/S1*004-KIR2DS4*006-KIR3DL2*005; KIR2DL4*00801-KIR3DL1/S1*00101-KIR2DS4*003-KIR3DL2*001; KIR2DL4*00103-KIR3DL1/S1*008-KIR2DS4*003-KIR3DL2*009; KIR2DL4*00102-KIR3DL1/S1*01502/*002-KIR2DS4*00101-KIR3DL2*002; KIR2DL4*00501-KIR3DL1/S1*013-KIR2DL5A*001-KIR2DS5*002-KIR2DS1*002-KIR3DL2*007). The present study provides a rationale for analyzing associations of KIR polymorphisms by taking into account the complex LD structure of the KIR region.     

Predominance of group A KIR haplotypes in Japanese associated with diverse NK cell repertoires of KIR expression

Genomic DNA from a panel of 41 healthy unrelated Japanese individuals was typed for the presence or absence of 16 KIR genes and pseudogenes. Only eight different KIR genotypes were found. Group A haplotypes outnumbered group B haplotypes in frequency by approximately 3:1, with individuals having two group A haplotypes accounting for 56% of the panel. The frequency of A haplotypes in the Japanese is higher than that observed in other populations. Flow cytometric comparison of KIR expression in 19 panel members showed considerable diversity in NK cell repertoire, which was also seen within the group of individuals having two A haplotypes. This diversity is likely due to allelic polymorphism in expressed genes of the A haplotype. In comparison to other populations, the Japanese appear less heterogeneous in KIR genotype as assessed by gene content.     

Recombination networks as genetic markers in a human variation study of the Old World

We have analyzed human genetic diversity in 33 Old World populations including 23 populations obtained through Genographic Project studies. A set of 1,536 SNPs in five X chromosome regions were genotyped in 1,288 individuals (mostly males). We use a novel analysis employing subARG network construction with recombining chromosomal segments. Here, a subARG is constructed independently for each of five gene-free regions across the X chromosome, and the results are aggregated across them. For PCA, MDS and ancestry inference with STRUCTURE, the subARG is processed to obtain feature vectors of samples and pairwise distances between samples. The observed population structure, estimated from the five short X chromosomal segments, supports genome-wide frequency-based analyses: African populations show higher genetic diversity, and the general trend of shared variation is seen across the globe from Africa through Middle East, Europe, Central Asia, Southeast Asia, and East Asia in broad patterns. The recombinational analysis was also compared with established methods based on SNPs and haplotypes. For haplotypes, we also employed a fixed-length approach based on information-content optimization. Our recombinational analysis suggested a southern migration route out of Africa, and it also supports a single, rapid human expansion from Africa to East Asia through South Asia.