North Asian population relationships in a global context.

Population genetic studies of North Asian ethnic groups have focused on genetic variation of sex chromosomes and mitochondria. Studies of the extensive variation available from autosomal variation have appeared infrequently. We focus on relationships among population samples using new North Asia microhaplotype data. We combined genotypes from our laboratory on 58 microhaplotypes, distributed across 18 autosomes, on 3945 individuals from 75 populations with corresponding data extracted for 26 populations from the Thousand Genomes consortium and for 22 populations from the GenomeAsia 100 K project. A total of 7107 individuals in 122 total populations are analyzed using STRUCTURE, Principal Component Analysis, and phylogenetic tree analyses. North Asia populations sampled in Mongolia include: Buryats, Mongolians, Altai Kazakhs, and Tsaatans. Available Siberians include samples of Yakut, Khanty, and Komi Zyriane. Analyses of all 122 populations confirm many known relationships and show that most populations from North Asia form a cluster distinct from all other groups. Refinement of analyses on smaller subsets of populations reinforces the distinctiveness of North Asia and shows that the North Asia cluster identifies a region that is ancestral to Native Americans.

The population genetics characteristics of a 90 locus panel of microhaplotypes.

Single-nucleotide polymorphisms (SNPs) and small genomic regions with multiple SNPs (microhaplotypes, MHs) are rapidly emerging as novel forensic investigative tools to assist in individual identification, kinship analyses, ancestry inference, and deconvolution of DNA mixtures. Here, we analyzed information for 90 microhaplotype loci in 4009 individuals from 79 world populations in 6 major biogeographic regions. The study included multiplex microhaplotype sequencing (mMHseq) data analyzed for 524 individuals from 16 populations and genotype data for 3485 individuals from 63 populations curated from public repositories. Analyses of the 79 populations revealed excellent characteristics for this 90-plex MH panel for various forensic applications achieving an overall average effective number of allele values (Ae) of 4.55 (range 1.04-19.27) for individualization and mixture deconvolution. Population-specific random match probabilities ranged from a low of 10-115 to a maximum of 10-66. Mean informativeness (In) for ancestry inference was 0.355 (range 0.117-0.883). 65 novel SNPs were detected in 39 of the MHs using mMHseq. Of the 3018 different microhaplotype alleles identified, 1337 occurred at frequencies > 5% in at least one of the populations studied. The 90-plex MH panel enables effective differentiation of population groupings for major biogeographic regions as well as delineation of distinct subgroupings within regions. Open-source, web-based software is available to support validation of this technology for forensic case work analysis and to tailor MH analysis for specific geographical regions.

Genetic relationships of Southwest Asian and Mediterranean populations.

The Southwest Asian, circum-Mediterranean, and Southern European populations (collectively, SWAMSE) together with Northern European populations form one of five "continental" groups of global populations in many analyses of population relationships. This region is of great anthropologic and forensic interest but relationships of large numbers of populations within the region have not been able to be cleanly resolved with autosomal genetic markers. To examine the genetic boundaries to the SWAMSE region and whether internal structure can be detected we have assembled data for a total of 151 separate autosomal genetic markers on populations in this region and other parts of the world for a global set of 95 populations. The markers include 83 ancestry informative SNPs as singletons and 68 microhaplotype loci defined by 204 SNPs. The 151 loci are ancestry informative on a global scale, identifying at least five biogeographic clusters. One of those clusters is a clear grouping of 37 populations containing the SWAMSE plus northern European populations to the exclusion of populations in South Central Asia and populations from farther East. A refined analysis of the 37 populations shows the northern European populations clustering separately from the SWAMSE populations. Within Southwest Asia the Samaritans and Shabaks are distinct outliers. The Yemenite Jews, Saudi, Kuwaiti, Palestinian Arabs, and Southern Tunisians cluster together loosely while the remaining populations from Northern Iraq, Mediterranean Europe, the Caucasus region, and Iran cluster in a more complex graded fashion. The majority of the SWAMSE populations from the mainland of Southwest Asia form a cluster with little internal structure reflecting a very complex history of endogamy and migrations. The set of 151 DNA polymorphisms not only distinguishes major geographical regions globally but can distinguish ancestry to a small degree within geographical regions such as SWAMSE. We discuss forensic characteristics of the polymorphisms and also identify those that rank highest by Rosenberg's In measure for the SWAMSE region populations and for the global set of populations analyzed. DATA AVAILABILITY: Genotypes on all 151 markers for all 3790 individuals typed in the Kidd Lab on the 72 Kidd lab populations have been deposited in the Zenodo archive and can be freely accessed at https://doi.org/10.5281/zenodo.4658892. Some of the data has been made public previously as supplemental files appended to publications. Data for the additional individuals included in the analyses was taken from already public datasets as indicated in the text.

The distinctive geographic patterns of common pigmentation variants at the OCA2 gene.

Oculocutaneous Albinism type 2 (OCA2) is a gene of great interest because of genetic variation affecting normal pigmentation variation in humans. The diverse geographic patterns for variant frequencies at OCA2 have been evident but have not been systematically investigated, especially outside of Europe. Here we examine population genetic variation in and near the OCA2 gene from a worldwide perspective. The very different patterns of genetic variation found across world regions suggest strong selection effects may have been at work over time. For example, analyses involving the variants that affect pigmentation of the iris argue that the derived allele of the rs1800407 single nucleotide polymorphism, which produces a hypomorphic protein, may have contributed to the previously demonstrated positive selection in Europe for the enhancer variant responsible for light eye color. More study is needed on the relationships of the genetic variation at OCA2 to variation in pigmentation in areas beyond Europe.

Validation of novel forensic DNA markers using multiplex microhaplotype sequencing.

Microhaplotypes (MH) are comprised of multiple single nucleotide polymorphisms (SNPs) that are located within 300 bases of genomic sequence. Improved tools are needed to facilitate broader application of microhaplotypes in a diverse range of populations and forensic settings. We designed an assay for multiplex sequencing of 90 microhaplotypes (mMHseq) that include 46 MH loci with high Effective Number of Alleles (Ae) from previous studies [1], and 44 high Ae MH loci containing between four to fourteen SNPs that were identified from the 1000 Genomes (1KG) Project. The unique design of mMHseq integrates a novel method for multiplex amplification from small DNA amounts, and multiplex sequencing of 48 samples in a single MiSeq run to detect all relevant MH variation. Assay performance was evaluated in a cohort of 156 individuals from seven different world populations from Africa, Asia, and Europe. Three of those populations from East Africa (Chagga, Sandawe, and Zaramo) and one from Eastern Europe (Adygei) had sufficient individuals sequenced by the assay to be included in statistical analyses with the 26 1KG populations. For those 30 populations the mean global average Ae was 5.08 (range: 2.7-11.54) and mean informativeness for biogeographic variation (In) was 0.30 (range: 0.08-0.70). Eighty-five novel SNPs were detected in 58 of the 90 microhaplotypes. Open-source, web-based software was developed to visualize haplotype phase data for each microhaplotype and individual. Our approach for multiplex microhaplotype sequencing can be customized and expanded as novel loci are being discovered.

Population relationships based on 170 ancestry SNPs from the combined Kidd and Seldin panels.

The benefits of ancestry informative SNP (AISNP) panels can best accrue and be properly evaluated only as sufficient reference population data become readily accessible. Ideally the set of reference populations should approximate the genetic diversity of human populations worldwide. The Kidd and Seldin AISNP sets are two panels that have separately accumulated thus far the largest and most diverse collections of data on human reference populations from the major continental regions. A recent tally in the ALFRED allele frequency database finds 164 reference populations available for all the 55 Kidd AISNPs and 132 reference populations for all the 128 Seldin AISNPs. Although much more of the genetic diversity in human populations around the world still needs to be documented, 81 populations have genotype data available for all 170 AISNPs in the union of the Kidd and Seldin panels. In this report we examine admixture and principal component analyses on these 81 worldwide populations and some regional subsets of these reference populations to determine how well the combined panel illuminates population relationships. Analyses of this dataset that focused on Native American populations revealed very strong cluster patterns associated with many of the individual populations studied.

Recent Selection on a Class I ADH Locus Distinguishes Southwest Asian Populations Including Ashkenazi Jews.

The derived human alcohol dehydrogenase (ADH)1B*48His allele of the ADH1B Arg48His polymorphism (rs1229984) has been identified as one component of an East Asian specific core haplotype that underwent recent positive selection. Our study has been extended to Southwest Asia and additional markers in East Asia. Fst values (Sewall Wright's fixation index) and long-range haplotype analyses identify a strong signature of selection not only in East Asian but also in Southwest Asian populations. However, except for the ADH2B*48His allele, different core haplotypes occur in Southwest Asia compared to East Asia and the extended haplotypes also differ. Thus, the ADH1B*48His allele, as part of a core haplotype of 10 kb, has undergone recent rapid increases in frequency independently in the two regions after divergence of the respective populations. Emergence of agriculture may be the common factor underlying the evident selection.

Selecting microhaplotypes optimized for different purposes.

Massively parallel sequencing is transforming forensic work by allowing various useful forensic markers, such as STRPs and SNPs, to be multiplexed providing information on ancestry, individual and familial identification, phenotypes for eye/hair/skin pigmentation, and the deconvolution of mixtures. Microhaplotypes also become feasible with massively parallel sequencing, these are DNA segments (smaller than 300 nucleotides) that are selected to contain multiple SNPs unambiguously defining three or more haplotype alleles occurring at common frequencies. The physical extent of a microhaplotype can thus be covered by a single sequence read making these loci phase-known codominant genetic systems. Such microhaplotypes supply significantly more information than a single SNP can. Our efforts to develop useful sets of microhaplotypes have already identified 182 such loci that we have studied on a large number of human populations from around the world. We present various analyses on 83 populations in our ongoing study for a subset of the best microhaplotypes currently available illustrating their characteristics and potential utility for ancestry, identification, and mixture deconvolution.

Improving ancestry distinctions among Southwest Asian populations.

The Kidd Lab panel of 55 AISNPs can provide up to 10 statistically relevant biogeographic groupings of a global set of populations. A second-tier panel would be useful for increasing the accuracy for further differentiation of populations within a specific global grouping. Because recent advances in massively parallel sequencing (MPS) methods allow the genotyping of many more SNPs, we are now identifying additional SNPs to provide refined discrimination among regional subsets of populations; Southwest Asia and the nearby Mediterranean region (SWA) is our current target for such a "second tier" panel. We selected the potentially best SNPs from various sources: our own laboratory database (>4600 SNPs), AISNP panels (Kidd 55 and Seldin 128 SNP panels), and published papers reporting European and SW Asian populations. Rosenberg's Informativeness, Fst, and allele frequency heatmap matrices are used to determine the best SNPs for the region. A total of 2568 individuals, from 39 different populations ranging from North-East Africa through the SW Asia and Europe to the Ural Mountains, were included in the refinement processes and analyses. Heatmap, PCA, Structure (K = 4), and ancestry inference for selected individuals with an in-lab version of FROG-kb analyses indicate that these 86 AISNPs provide the basis for building an improved, optimized panel of AISNPs that collectively provide additional information on differences among populations in that part of the world. Testing this panel with additional populations from the area and with new SNPs and/or microhaplotypes is expected to improve the panel.

Ancestry inference of 96 population samples using microhaplotypes.

Microhaplotypes have become a new type of forensic marker with a great ability to identify and deconvolute mixtures because massively parallel sequencing (MPS) allows the alleles (haplotypes) of the multi-SNP loci to be determined directly for an individual. As originally defined, a microhaplotype locus is a short segment of DNA with two or more SNPs defining three or more haplotypes. The length is short enough, less than about 300 bp, that the read length of current MPS technology can produce a phase-known sequence of each chromosome of an individual. As part of the discovery phase of our studies, data on 130 microhaplotype loci with estimates of haplotype frequency data on 83 populations have been published. To provide a better picture of global allele frequency variation, we have now tested 13 more populations for 65 of the microhaplotype loci from among those with higher levels of inter-population gene frequency variation, including 8 loci not previously published. These loci provide clear distinctions among 6 biogeographic regions and provide some information distinguishing up to 10 clusters of populations.

Evaluating 130 microhaplotypes across a global set of 83 populations.

Today the primary DNA markers used in forensics are short tandem repeat (STR) polymorphisms (STRPs), initially selected because they are highly polymorphic. However, the increasingly common need to deal with samples with a mixture of DNA from two or more individuals sometimes is complicated by the inherent stutter involved with PCR amplification, especially in strongly unbalanced mixtures when the minor component coincides with the stutter range of the major component. Also, the STRPs in use provide little evidence of ancestry of a single source sample beyond broad "continental" resolution. Methodologies for analyzing DNA have become much more powerful in recent years. Massively parallel sequencing (MPS) is a new method being considered for routine use in forensics. Primarily to aid in mixture deconvolution and avoid the issue of stutter, we have begun to investigate a new type of forensic marker, microhaplotype loci, that will provide useful information on mixtures of DNA and on ancestry when typed using massively parallel sequencing (MPS). We have identified 130 loci and estimated their haplotype (allele) frequencies in 83 different population samples. Many of these loci are shown to be highly informative for individual identification and for mixture identification and deconvolution.

Genetic variation in Tunisia in the context of human diversity worldwide.

OBJECTIVES: North Africa has a complex demographic history of migrations from within Africa, Europe, and the Middle East. However, population genetic studies, especially for autosomal genetic markers, are few relative to other world regions. We examined autosomal markers for eight Tunisian and Libyan populations in order to place them in a global context. MATERIALS AND METHODS: Data were collected by TaqMan on 399 autosomal single nucleotide polymorphisms on 331 individuals from Tunisia and Libya. These data were combined with data on the same SNPs previously typed on 2585 individuals from 57 populations from around the world. Where meaningful, close by SNPs were combined into multiallelic haplotypes. Data were evaluated by clustering, principal components, and population tree analyses. For a subset of 102 SNPs, data from the literature on seven additional North African populations were included in analyses. RESULTS: Average heterozygosity of the North African populations is high relative to our global samples, consistent with a complex demographic history. The Tunisian and Libyan samples form a discrete cluster in the global and regional views and can be separated from sub-Sahara, Middle East, and Europe. Within Tunisia the Nebeur and Smar are outlier groups. Across North Africa, pervasive East-West geographical patterns were not found. DISCUSSION: Known historical migrations and invasions did not displace or homogenize the genetic variation in the region but rather enriched it. Even a small region like Tunisia contains considerable genetic diversity. Future studies across North Africa have the potential to increase our understanding of the historical demographic factors influencing the region. Am J Phys Anthropol 161:62-71, 2016. © 2016 The Authors American Journal of Physical Anthropology Published by Wiley Periodicals, Inc.

Criteria for selecting microhaplotypes: mixture detection and deconvolution.

BACKGROUND: DNA sequencing is likely to become a standard typing method in forensics in the near future. We define a microhaplotype to be a locus with two or more single nucleotide polymorphisms (SNPs) that occur within a short segment of DNA (e.g., 200 bp) that can be covered by a single sequence run and collectively define a multiallelic locus. Microhaplotypes can be highly informative for many forensic questions, including detection of mixtures of two or more sources in a DNA sample, a common problem in forensic practice. RESULTS: When all alleles are equally frequent, the probability of detecting three or more alleles in a mixture is at maximum. The classical population genetics concept of effective number of alleles at a locus, termed Ae, converts the unequal allele frequencies at a locus into a value that is equivalent to some number of equally frequent alleles, allowing microhaplotype loci to be ranked. The expectations for the ability to qualitatively detect mixtures are given for different integer values of Ae, and the cumulative probabilities of detecting mixtures based on testing multiple microhaps are shown to exceed 95% with as few as five loci with average Ae values of even slightly greater than 3.0. CONCLUSIONS: Microhaplotypes with Ae values of >3 will be exceedingly useful in ordinary forensic practice. Based on our studies, 3-SNP microhaplotypes will sometimes meet this criterion, but 4-SNP microhaplotypes can even exceed this criterion and have values >4.

Current sequencing technology makes microhaplotypes a powerful new type of genetic marker for forensics.

SNPs that are molecularly very close (<10kb) will generally have extremely low recombination rates, much less than 10(-4). Multiple haplotypes will often exist because of the history of the origins of the variants at the different sites, rare recombinants, and the vagaries of random genetic drift and/or selection. Such multiallelic haplotype loci are potentially important in forensic work for individual identification, for defining ancestry, and for identifying familial relationships. The new DNA sequencing capabilities currently available make possible continuous runs of a few hundred base pairs so that we can now determine the allelic combination of multiple SNPs on each chromosome of an individual, i.e., the phase, for multiple SNPs within a small segment of DNA. Therefore, we have begun to identify regions, encompassing two to four SNPs with an extent of <200bp that define multiallelic haplotype loci. We have identified candidate regions and have collected pilot data on many candidate microhaplotype loci. Here we present 31 microhaplotype loci that have at least three alleles, have high heterozygosity, are globally informative, and are statistically independent at the population level. This study of microhaplotype loci (microhaps) provides proof of principle that such markers exist and validates their usefulness for ancestry inference, lineage-clan-family inference, and individual identification. The true value of microhaplotypes will come with sequencing methods that can establish alleles unambiguously, including disentangling of mixtures, because a single sequencing run on a single strand of DNA will encompass all of the SNPs.

Extensive sequence variation in the 3' untranslated region of the KRAS gene in lung and ovarian cancer cases.

While cancer is a serious health issue, there are very few genetic biomarkers that predict predisposition, prognosis, diagnosis, and treatment response. Recently, sequence variations that disrupt microRNA (miRNA)-mediated regulation of genes have been shown to be associated with many human diseases, including cancer. In an early example, a variant at one particular single nucleotide polymorphism (SNP) in a let-7 miRNA complementary site in the 3' untranslated region (3' UTR) of the KRAS gene was associated with risk and outcome of various cancers. The KRAS oncogene is an important regulator of cellular proliferation, and is frequently mutated in cancers. To discover additional sequence variants in the 3' UTR of KRAS with the potential as genetic biomarkers, we resequenced the complete region of the 3' UTR of KRAS in multiple non-small cell lung cancer and epithelial ovarian cancer cases either by Sanger sequencing or capture enrichment followed by high-throughput sequencing. Here we report a comprehensive list of sequence variations identified in cases, with some potentially dysregulating expression of KRAS by altering putative miRNA complementary sites. Notably, rs712, rs9266, and one novel variant may have a functional role in regulation of KRAS by disrupting complementary sites of various miRNAs, including let-7 and miR-181.

Progress toward an efficient panel of SNPs for ancestry inference.

Many panels of ancestry informative single nucleotide polymorphisms have been proposed in recent years for various purposes including detecting stratification in biomedical studies and determining an individual's ancestry in a forensic context. All of the panels have limitations in their generality and efficiency for routine forensic work. Some panels have used only a few populations to validate them. Some panels are based on very large numbers of SNPs thereby limiting the ability of others to test different populations. We have been working toward an efficient and globally useful panel of ancestry informative markers that is comprised of a small number of highly informative SNPs. We have developed a panel of 55 SNPs analyzed on 73 populations from around the world. We present the details of the panel and discuss its strengths and limitations.

Haplotype structure and positive selection at TLR1.

Toll-like receptor 1, when dimerized with Toll-like receptor 2, is a cell surface receptor that, upon recognition of bacterial lipoproteins, activates the innate immune system. Variants in TLR1 associate with the risk of a variety of medical conditions and diseases, including sepsis, leprosy, tuberculosis, and others. The foremost of these is rs5743618 c.2079T>G(p.(Ile602Ser)), the derived allele of which is associated with reduced risk of sepsis, leprosy, and other diseases. Interestingly, 602Ser, which shows signatures of selection, inhibits TLR1 surface trafficking and subsequent activation of NFκB upon recognition of a ligand. This suggests that reduced TLR1 activity may be beneficial for human health. To better understand TLR1 variation and its link to human health, we have typed all 7 high-frequency missense variants (>5% in at least one population) along with 17 other variants in and around TLR1 in 2548 individuals from 56 populations from around the globe. We have also found additional signatures of selection on missense variants not associated with rs5743618, suggesting that there may be multiple functional alleles under positive selection in this gene.

Worldwide population variation and haplotype analysis at the serotonin transporter gene SLC6A4 and implications for association studies.

BACKGROUND: Variation at the serotonin transporter gene, SLC6A4, has been associated with a variety of neuropsychiatric disorders and could be involved in other health-related phenotypes. METHODS: To determine the extent of variation at SLC6A4, we genotyped 23 markers on approximately 2500 individuals from 47 global populations, including the promoter variable number tandem repeat (VNTR) and 2 single nucleotide polymorphisms (SNPs) immediately flanking its variable region (rs25531 and rs25532), the intron 2 VNTR, and 19 additional SNPs. RESULTS: We observed several rare alleles at the promoter VNTR (some novel) and population-specific distributions of the reported functional SNPs rs25531, rs25532, and rs6355, as well as two alleles at the intron 2 VNTR. Alleles of interest at the VNTRs occurred on specific haplotype backgrounds. The repeat-number variants at the promoter VNTR and the intron 2 VNTR, as well as the putative functional SNPs, showed ethnic variation in frequencies. The more common alleles at the VNTR polymorphisms show wide geographic distributions, whereas rare alleles at both show more restricted distributions. The derived alleles at the two functional SNPs in the promoter VNTR show restricted distributions and occur primarily on different repeat number alleles. CONCLUSIONS: Our findings illustrate significant variation worldwide at SLC6A4 and that the functionally implicated alleles at the SNPs rs25531, rs25532, and rs6355 occur on limited haplotypes and vary significantly in global distribution. Association studies at SLC6A4 cannot a priori extrapolate across populations and should account for the multiple polymorphisms with possible functional variation across this locus, rather than focusing solely on one or two polymorphisms as commonly seen.

Exploring genomic structure differences and similarities between the Greek and European HapMap populations: implications for association studies.

Studies of the genomic structure of the Greek population and Southeastern Europe are limited, despite the central position of the area as a gateway for human migrations into Europe. HapMap has provided a unique tool for the analysis of human genetic variation. Europe is represented by the CEU (Northwestern Europe) and the TSI populations (Tuscan Italians from Southern Europe), which serve as reference for the design of genetic association studies. Furthermore, genetic association findings are often transferred to unstudied populations. Although initial studies support the fact that the CEU can, in general, be used as reference for the selection of tagging SNPs in European populations, this has not been extensively studied across Europe. We set out to explore the genomic structure of the Greek population (56 individuals) and compare it to the HapMap TSI and CEU populations. We studied 1112 SNPs (27 regions, 13 chromosomes). Although the HapMap European populations are, in general, a good reference for the Greek population, regions of population differentiation do exist and results should not be light-heartedly generalized. We conclude that, perhaps due to the individual evolutionary history of each genomic region, geographic proximity is not always a perfect guide for selecting a reference population for an unstudied population.

Expanding data and resources for forensic use of SNPs in individual identification.

The potential value of SNPs for individual identification has been recognized by many researchers and different panels have been proposed. Here we present a new interface in the ALFRED database to access compendia of allele frequencies for several published panels of markers for forensic uses. One of those is our panel of individual identification SNPs (IISNPs) based on samples of 44 populations originating from many parts of the world. Here we also present additional data and additional statistical analyses that continue to support the value of our panel of IISNPs as a universal panel. We also describe initial developments of multiplex methods and various robustness analyses for our 45 marker IISNP panel.