Ancestry Informative SNPs Research Articles

Forensic ancestry analysis in two Chinese minority populations using massively parallel sequencing of 165 ancestry‐informative SNPs

Ancestry inference based on SNPs with marked allele frequency differences in diverse populations (called ancestry-informative SNP (AISNP)) is rapidly developed with the technology advancements of massively parallel sequencing. Despite the decade of exploration and broad public interest in the peopling of East-Asians, the genetic landscape of Chinese Silk Road populations based on the AISNPs is still little known. In this work, 206 unrelated individuals from Chinese Uyghur and Hui populations were firstly genotyped by 165 AISNPs (The Precision ID Ancestry Panel) using the Ion Torrent Personal Genome Machine system. The ethnic origin of two investigated populations and population structures and genetic relationships were subsequently investigated. The 165 AISNPs panel not only can differentiate Uyghur and Hui populations but also has potential applications in individual identification. Comprehensive population comparisons and admixture estimates demonstrated a predominantly higher European-related ancestry (36.30%) in Uyghurs than Huis (3.66%). Overall, the Precision ID Ancestry Panel can provide good resolution at the intercontinental level, but has limitations on the genetic homogeneous populations, such as the Hui and Han. Additional population-specific AISNPs remain necessary to get better-scale resolution within geographically proximate populations in East Asia.

High sample throughput genotyping for estimating C-lineage introgression in the dark honeybee: an accurate and cost-effective SNP-based tool

The natural distribution of the honeybee (Apis mellifera L.) has been changed by humans in recent decades to such an extent that the formerly widest-spread European subspecies, Apis mellifera mellifera, is threatened by extinction through introgression from highly divergent commercial strains in large tracts of its range. Conservation efforts for A. m. mellifera are underway in multiple European countries requiring reliable and cost-efficient molecular tools to identify purebred colonies. Here, we developed four ancestry-informative SNP assays for high sample throughput genotyping using the iPLEX Mass Array system. Our customized assays were tested on DNA from individual and pooled, haploid and diploid honeybee samples extracted from different tissues using a diverse range of protocols. The assays had a high genotyping success rate and yielded accurate genotypes. Performance assessed against whole-genome data showed that individual assays behaved well, although the most accurate introgression estimates were obtained for the four assays combined (117 SNPs). The best compromise between accuracy and genotyping costs was achieved when combining two assays (62 SNPs). We provide a ready-to-use cost-effective tool for accurate molecular identification and estimation of introgression levels to more effectively monitor and manage A. m. mellifera conservatories.

Massively parallel sequencing of 165 ancestry informative SNPs in two Chinese Tibetan-Burmese minority ethnicities

The Tibeto-Burman language, one subfamily of the Sino-Tibetan languages, is spoken by over 60 million people all over East Asia. Yet the ethnic origin and genetic architecture of Tibeto-Burman speaking populations remain largely unexplored. In the present study, 169 Chinese individuals from Tibeto-Burman speaking populations (two ethnic groups: Tibetan and Yi) in four different geographic regions in western China were analyzed using the Precision ID Ancestry Panel (165 AISNPs) and the Ion PGM System. The performance and corresponding forensic statistical parameters of this AISNPs panel were investigated. Comprehensive population genetic comparisons (143 populations based on Kidd’ SNPs, 92 populations on the basis of Seldin’ SNPs and 31 populations based on the Precision ID Ancestry Panel) and ancestry inference were further performed. Sequencing performance demonstrated that the Precision ID Ancestry Panel is effective and robust. Forensic characteristics suggested that this panel not only can be used for ancestry estimation of Tibeto-Burman populations but also for individual identification. Tibetan and Yi shared a common genetic ancestry origin but experienced the complex history of gene flow, local adaptation, and isolation, and constructed the specific genetic landscape of human genetic diversity of Highlander and Lowlander populations. Tibetan-Burman populations and other East Asian populations showed sufficient genetic difference and could be distinguished into three distinct groups. Furthermore, analysis of population structure revealed that significant genetic difference was existed inter-continent populations and strong genetic affinity was observed within-continent populations. Additional population-specific AISNPs and a relatively more comprehensive database with sufficient reference population data remain necessary to get better-scale resolution within a geographically proximate populations in East Asia.

RETRACTED ARTICLE: Differentiation analysis for estimating individual ancestry from the Tibetan Plateau by an archaic altitude adaptation EPAS1 haplotype among East Asian populations

Tibetans have adapted to the extreme environment of high altitude for hundreds of generations. A highly differentiated 5-SNP (Single Nucleotide Polymorphism) haplotype motif (AGGAA) on a hypoxic pathway gene, EPAS1, is observed in Tibetans and lowlanders. To evaluate the potential usage of the 5-SNP haplotype in ancestry inference for Tibetan or Tibetan-related populations, we analyzed this haplotype in 1053 individuals of 12 Chinese populations residing on the Tibetan Plateau, peripheral regions of Tibet, and plain regions. These data were integrated with the genotypes from the 1000 Genome populations and populations in a previously reported paper for population structure analyses. We found that populations representing highland and lowland groups have different dominant ancestry components. The core Denisovan haplotype (AGGAA) was observed at a frequency of 72.32% in the Tibetan Plateau, with a frequency range from 9.48 to 21.05% in the peripheral regions and < 2.5% in the plains area. From the individual perspective, 87.57% of the individuals from the Tibetan Plateau carried the archaic haplotype, while < 5% of the Chinese Han people carried the haplotype. Our findings indicate that the 5-SNP haplotype has a special distribution pattern in populations of Tibet and peripheral regions and could be integrated into AISNP (Ancestry Informative Single Nucleotide Polymorphism) panels to enhance ancestry resolution.

A Locus on Chromosome 5 Is Associated with Red Cell Allo-Immunization in Adults with Sickle Cell Disease

Background: Transfusion is a standard treatment for acute and chronic complications of sickle cell disease (SCD). This therapy is limited by a high prevalence of red cell antigen alloimmunization, which is not seen in other frequently transfused patient populations. Reasons for this increased prevalence have not been elucidated. One model based study suggested a genetic predisposition in all patient populations. Other studies have provided evidence that antigen mismatches occur commonly between blood donors of European ancestry and SCD patients and that this is a primary cause for alloimmunization. Accordingly, limited phenotypic antigen matching is a standard practice used to minimize alloimmunization risk. We hypothesized that if these population differences are the biologic basis for this phenomenon that admixture mapping could identify significant susceptibility loci.Methods: A cohort study consisting of adults evaluated at NIH and Howard University between 2001-2015 (n=721) had allo-antibody status confirmed with the blood bank at the hospital where the most recent transfusion was performed. Genomic DNA was genotyped for a panel of 1622 ancestry informative SNPs in 359 patients. Genotypes were used to infer ancestry proportions using Structure software assuming a 3 population model. Ancestry was analyzed by regression along with clinical variables. Genome wide admixture mapping was performed to identify regions of the genome where ancestry differences between cases and controls suggested an association. We used AncestryMap for a log-additive case control study assuming various risk models.Results: Five hundred fifty four subjects (77% of cohort) had a transfusion history with 153 (27.6 %) with allo-antibodies and 65 (11.7 %) with no prior transfusions. Antibodies directed against RH, Kell, MNS and Duffy antigens were the most prevalent. Clinical factors associated with alloimmunization in regression models were only significant for lifetime transfusion history by both logistic (antibody status, β=0.56, P=6.31 x 10-6) and linear (number of antibodies, β=0.35, P=0.0004) regressions. When considering genotype data, Amerindian ancestry was associated with antibody status (logistic, β=4.60, P=0.03) and number of antibodies (linear, β=4.92, P=0.006), while African and European ancestry were not. The association between ancestry and alloimmunization risk remained unchanged when the number of prior transfusions was also included in both multivariable models for alloimmunization status and number of antibodies, P=0.09 and P=0.06, respectively. Genome wide mapping by admixture linkage disequilibrium compared 91 allo-Ab cases to 210 controls with prior transfusions under 10 different risk models (range 0.25-9.0). The most significant result was the 0.25 risk model which was highly suggestive for association (genome log factor 10.96, P<0.001; Figure 1A). This scan identified region of chromosome 5 (5q21.3-5q31.1) defined by 14 snps spanning 21.685 Mb and a highly significant peak LOD score of 13.5 (Figure 1B). While this region contains >540 genes, the linkage region is noteworthy for a cytokine gene cluster implicated in hematopoiesis and immune responses.Conclusions: In a large cohort of adults with SCD, the prevalence of red cell alloimmunization is 27.6%. Amerindian ancestry is a risk factor and suggests a genetic component to this phenotype. The number of prior transfusions is also a significant risk factor. Genome wide admixture mapping identified a strong association at chromosome 5q21.3-5q31.1 using a 0.25 risk model. This chromosome 5 association is particularly interesting given the presence a cytokine gene cluster within a portion of the associated region. Analysis of larger, admixed populations will help to fine map allo-immunization susceptibility variants within specific genes that underlie these loci. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Massively parallel sequencing of forensic STRs and SNPs using the Illumina® ForenSeq™ DNA Signature Prep Kit on the MiSeq FGx™ Forensic Genomics System

The ForenSeq™ DNA Signature Prep Kit (ForenSeq Kit) is designed to detect more than 200 forensically relevant markers in a single reaction on the MiSeq FGx™ Forensic Genomics System (MiSeq FGx System), including Amelogenin, 27 autosomal short tandem repeats (A-STRs), 7 X chromosomal STRs (X-STRs), 24 Y chromosomal STRs (Y-STRs) and 94 identity-informative single nucleotide polymorphisms (iSNPs) with the option to contain 22 phenotypic-informative SNPs (pSNPs) and 56 ancestry-informative SNPs (aSNPs). In this study, we evaluated the MiSeq FGx System on three major parts: methodological optimization (DNA extraction, sample quantification, library normalization, diluted libraries concentration, and sample-to-cell arrangement), massively parallel sequencing (MPS) performance (depth of coverage, sequence coverage ratio, and allele coverage ratio), and ForenSeq Kit characteristics (repeatability and concordance, sensitivity, mixture, stability and case-type samples). Results showed that quantitative polymerase chain reaction (qPCR)-based sample quantification and library normalization and the appropriate number of pooled libraries and concentration of diluted libraries provided a greater level of MPS performance and repeatability. Repeatable and concordant genotypes were obtained by the ForenSeq Kit. Full profiles were obtained from ≥100pg input DNA for STRs and ≥200pg for SNPs. A sample with ≥5% minor contributors was considered as a mixture by imbalanced allele coverage ratio distribution, and full profiles from minor contributors were easily detected between 9:1 and 1:9 mixtures with known reference profiles. The ForenSeq Kit tolerated considerable concentrations of inhibitors like ≤200μM hematin and ≤50μg/ml humic acid, and >56% STR profiles and >88% SNP profiles were obtained from ≥200-bp degraded samples. Also, it was adapted to case-type samples. As a whole, the ForenSeq Kit is a well-performed, robust, reliable, reproducible and highly informative assay, and it can fully meet requirements for human identification. Further, sensitive QC indicator and automated sample comparison function in the ForenSeq™ Universal Analysis Software are quite helpful, so that we can concentrate on questionable genotypes and avoid tedious and time-consuming labor to maximum the time spent in data analysis.

Ancestry prediction in Singapore population samples using the Illumina ForenSeq kit

The ability to predict bio-geographic ancestry can be valuable to generate investigative leads towards solving crimes. Ancestry informative marker (AIM) sets include large numbers of SNPs to predict an ancestral population. Massively parallel sequencing has enabled forensic laboratories to genotype a large number of such markers in a single assay. Illumina’s ForenSeq DNA Signature Kit includes the ancestry informative SNPs reported by Kidd et al.In this study, the ancestry prediction capabilities of the ForenSeq kit through sequencing on the MiSeq FGx were evaluated in 1030 unrelated Singapore population samples of Chinese, Malay and Indian origin. A total of 59 ancestry SNPs and phenotypic SNPs with AIM properties were selected. The bio-geographic ancestry of the 1030 samples, as predicted by Illumina’s ForenSeq Universal Analysis Software (UAS), was determined. 712 of the genotyped samples were used as a training sample set for the generation of an ancestry prediction model using STRUCTURE and Snipper. The performance of the prediction model was tested by both methods with the remaining 318 samples.Ancestry prediction in UAS was able to correctly classify the Singapore Chinese as part of the East Asian cluster, while Indians clustered with Ad-mixed Americans and Malays clustered in-between these two reference populations. Principal component analyses showed that the 59 SNPs were only able to account for 26% of the variation between the Singapore sub-populations. Their discriminatory potential was also found to be lower (GST=0.085) than that reported in ALFRED (FST=0.357). The Snipper algorithm was able to correctly predict bio-geographic ancestry in 91% of Chinese and Indian, and 88% of Malay individuals, while the success rates for the STRUCTURE algorithm were 94% in Chinese, 80% in Malay, and 91% in Indian individuals. Both these algorithms were able to provide admixture proportions when present. Ancestry prediction accuracy (in terms of likelihood ratio) was generally high in the absence of admixture. Misclassification occurred in admixed individuals, who were likely offspring of inter-ethnic marriages, and hence whose self-reported bio-geographic ancestries were dependent on that of their fathers, and in individuals of minority sub-populations with inter-ethnic beliefs. The ancestry prediction capabilities of the 59 SNPs on the ForenSeq kit were reasonably effective in differentiating the Singapore Chinese, Malay and Indian sub-populations, and will be of use for investigative purposes. However, there is potential for more accurate prediction through the evaluation of other AIM sets.

Population and performance analyses of four major populations with Illumina’s FGx Forensic Genomics System

The MiSeq FGx Forensic Genomics System (Illumina) enables amplification and massively parallel sequencing of 59 STRs, 94 identity informative SNPs, 54 ancestry informative SNPs, and 24 phenotypic informative SNPs. Allele frequency and population statistics data were generated for the 172 SNP loci included in this panel on four major population groups (Chinese, African Americans, US Caucasians, and Southwest Hispanics). Single-locus and combined random match probability values were generated for the identity informative SNPs. The average combined STR and identity informative SNP random match probabilities (assuming independence) across all four populations were 1.75E-67 and 2.30E-71 with length-based and sequence-based STR alleles, respectively. Ancestry and phenotype predictions were obtained using the ForenSeq™ Universal Analysis System (UAS; Illumina) based on the ancestry informative and phenotype informative SNP profiles generated for each sample. Additionally, performance metrics, including profile completeness, read depth, relative locus performance, and allele coverage ratios, were evaluated and detailed for the 725 samples included in this study. While some genetic markers included in this panel performed notably better than others, performance across populations was generally consistent. The performance and population data included in this study support that accurate and reliable profiles were generated and provide valuable background information for laboratories considering internal validation studies and implementation.

A panel of 130 autosomal single-nucleotide polymorphisms for ancestry assignment in five Asian populations and in Caucasians.

Ancestry informative single-nucleotide polymorphism (AISNP) panels for differentiating between East and Southeast Asian populations are scarce. This study aimed to identify AISNPs for ancestry assignment of five East and Southeast Asian populations, and Caucasians. We analyzed 145 autosomal SNPs of the 627 DNA samples from individuals of six populations (234 Taiwanese Han, 91 Filipinos, 79 Indonesians, 60 Thais, 71 Vietnamese, and 92 Caucasians) using arrays. The multiple logistic regression model and a multi-tier approach were used for ancestry classification. We observed that 130 AISNPs were effective for classifying the ethnic origins with fair accuracy. Among the 130 AISNPs, 122 were useful for stratification between these five Asian populations and 64 were effective for differentiating between Caucasians and these Asian populations. For differentiation between Caucasians and Asians, an accuracy rate of 100% was achieved in these 627 subjects with 50 optimal AISNPs among the 64 effective SNPs. For classification of the five Asian populations, the accuracy rates of ancestry inference using 20 to 57 SNPs for each of the two Asian populations ranged from 74.1% to 100%. Another 14 degraded DNA samples with incomplete profiling were analyzed, and the ancestry of 12 (85.7%) of those subjects was accurately assigned. We developed a 130-AISNP panel for ethnic origin differentiation between the five East and Southeast Asian populations and Caucasians. This AISNP set may be helpful for individual ancestral assignment of these populations in forensic casework.

Flanking region variation of ForenSeq™ DNA Signature Prep Kit STR and SNP loci in Yavapai Native Americans

Massively parallel sequencing (MPS) offers advantages over current capillary electrophoresis-based analysis of short tandem repeat (STR) loci for human identification testing. In particular STR repeat motif sequence information can be obtained, thereby increasing the discrimination power of some loci. While sequence variation within the repeat region is observed relatively frequently in some of the commonly used STRs, there is an additional degree of variation found in the flanking regions adjacent to the repeat motif. Repeat motif and flanking region sequence variation have been described for major population groups, however, not for more isolated populations. Flanking region sequence variation in STR and single nucleotide polymorphism (SNP) loci in the Yavapai population was analyzed using the ForenSeq™ DNA Signature Prep Kit and STRait Razor v2s. Seven and 14 autosomal STRs and identity-informative single nucleotide polymorphisms (iiSNPs), respectively, had some degree of flanking region variation. Three and four of these identity-informative loci, respectively, showed ≥5% increase in expected heterozygosity. The combined length- and sequence-based random match probabilities (RMPs) for 27 autosomal STRs were 6.11×10−26 and 2.79×10−29, respectively. When combined with 94 iiSNPs (a subset of which became microhaplotypes) the combined RMP was 5.49×10−63. Analysis of length-based and sequence-based autosomal STRs in STRUCTURE indicated that the Yavapai are most similar to the Hispanic population. While producing minimal increase in X- and Y-STR discrimination potential, access to flanking region data enabled identification of one novel X-STR and three Y-STR alleles relative to previous reports. Five ancestry-informative SNPs (aiSNPs) and two phenotype-informative SNPs (piSNPs) exhibited notable flanking region variation.

Evaluation of the Precision ID Ancestry Panel for crime case work: A SNP typing assay developed for typing of 165 ancestral informative markers

The application of massive parallel sequencing (MPS) methodologies in forensic genetics is promising and it is gradually being implemented in forensic genetic case work. One of the major advantages of these technologies is that several traditional electrophoresis assays can be combined into one single MPS assay. This reduces both the amount of sample used and the time of the investigations.This study assessed the utility of the Precision ID Ancestry Panel (Thermo Fisher Scientific, Waltham, USA) in forensic genetics. This assay was developed for the Ion Torrent PGM™ System and genotypes 165 ancestry informative SNPs. The performance of the assay and the accompanying software solution for ancestry inference was assessed by typing 142 Danes and 98 Somalis. Locus balance, heterozygote balance, and noise levels were calculated and future analysis criteria for crime case work were estimated. Overall, the Precision ID Ancestry Panel performed well, and only minor changes to the recommended protocol were implemented. Three out of the 165 loci (rs459920, rs7251928, and rs7722456) had consistently poor performance, mainly due to misalignment of homopolymeric stretches. We suggest that these loci should be excluded from the analyses.The different statistical methods for reporting ancestry in forensic genetic case work are discussed.

Increasing the reference populations for the 55 AISNP panel: the need and benefits

Ancestry inference for an individual can only be as good as the reference populations with allele frequency data on the SNPs being used. If the most relevant ancestral population(s) does not have data available for the SNPs studied, then analyses based on DNA evidence may indicate a quite distantly related population, albeit one among the more closely related of the existing reference populations. We have added reference population allele frequencies for 14 additional population samples (with >1100 individuals studied) to the 125 population samples previously published for the Kidd Lab 55 AISNP panel. Allele frequencies are now publicly available for all 55 SNPs in ALFRED and FROG-kb for a total of 139 population samples. This Kidd Lab panel of 55 ancestry informative SNPs has been incorporated in commercial kits by both ThermoFisher Scientific and Illumina for massively parallel sequencing. Researchers employing those kits will find the enhanced set of reference populations useful.

Using Whole-Genome Sequence Information to Foster Conservation Efforts for the European Dark Honey Bee, Apis mellifera mellifera

Pollination is a key ecosystem service for agricultural systems and Western honey bees, Apis mellifera, are the most important managed pollinators. Major losses of managed honey bee colonies reinforced the need to take advantage of locally adapted subspecies and ecotypes to buffer populations against various stressors. However, introductions of non-native honey bees from distant lineages are likely to undermine respective conservation efforts unless reliable and cost effective tools can be used to identify hybridization. The purpose of this study is to characterize current population structure and genetic diversity, and to assess the degree of admixture between native and introduced honey bees. Moreover, we aim to select a reduced number of genetic markers to improve conservation management strategies. We take advantage of recent developments in next-generation sequencing and network-based clustering to investigate conservation efforts for the native European Dark honey bee, A. m. mellifera, which is threatened by introgression in most of its range. We collected whole-genome sequence information from haploid drones of A. m. mellifera, A. m. carnica and Buckfast sampled throughout Switzerland (N=81), as well as from four Swiss A. m. mellifera conservation areas (N=39) and from one conservatory in the French Alps (N = 31). Population structure analyses based upon 3.375 M genome-wide SNPs discerned samples by subspecies and geographic origin (Switzerland or France). Ancestry inference indicated admixed individuals in all of the protected areas, calling for improved management efforts. After testing different subsets of ancestry informative SNPs using three different selection strategies (FST, PCA-based or at random), as few as 50 SNPs are found to be sufficient to differentiate native from introduced honey bees. Therefore our data suggests that a low-density SNP panel can be a precise and cost-effective tool to support conservation management efforts for managed pollinators.

Prediction of biogeographical ancestry from genotype: a comparison of classifiers.

DNA can provide forensic intelligence regarding a donor's biogeographical ancestry (BGA) and other externally visible characteristics (EVCs). A number of algorithms have been proposed to assign individual human genotypes to a BGA using ancestry informative marker (AIM) panels. This study compares the BGA assignment accuracy of the population clustering program STRUCTURE and three generic classification approaches including a Bayesian algorithm, genetic distance, and multinomial logistic regression (MLR). A selection of 142 ancestry informative single nucleotide polymorphisms (SNPs) were chosen from existing marker panels (SNPforID 34-plex, Eurasiaplex, Seldin, and Kidd's AIM panels) to assess BGA classification at the continental level for Africans, Europeans, East Asians, and Amerindians. A training set of 1093 individuals with self-declared BGA from the 1000 Genomes phase 1 database was used by each classifier to predict BGA in a test set of 516 individuals from the HGDP-CEPH (Stanford) cell line panel. Tests were repeated with 0, 10, 50, 70, and 90% of the genotypes missing. Comparison of the area under the receiver operating characteristic curves (AUROCs) showed high accuracy in STRUCTURE and the generic Bayesian approach. The latter algorithm offers a computationally simpler alternative to STRUCTURE with little loss in accuracy and is suitable for phenotype prediction while STRUCTURE is not.

Amerindian/Asian Ancestry and Mortality Are Associated with Allo-Immunization in Adults with Sickle Cell Disease in a Genome Wide Racial Admixture Study

▪Background: Transfusion is a standard therapy for acute and chronic complications of sickle cell disease (SCD). They are complicated by a high prevalence of red cell antigen alloimmunization, which is not seen in other frequently transfused patients. Reasons for this increased prevalence are unclear. Some studies have suggested genetic predisposition to alloimmunization in all transfused patient populations. Others suggested that mismatched red cell antigens occurring commonly between blood donors of European ancestry and patients with SCD who are of African ancestry is a primary cause in SCD. Accordingly, phenotypic antigen matching and transfusion from African American donors are proposed strategies to minimize alloimmunization risk. We hypothesized that if these population differences are the etiologic basis that analysis of ancestry at the genetic level could identify significant genetic susceptibility contributions and localize susceptibility loci.Methods: The Bethesda Sickle Cell Cohort Study consisting of adults evaluated at NIH from 2001-2015 (n=721) was reviewed for transfusion history and the presence of allo-antibodies confirmed with a blood bank. Genomic DNA was genotyped with an admixture panel consisting of 2251 ancestry informative SNPs in 359 patients and in HapMap controls. A final panel of 1622 SNPs was used to infer ancestry proportions using Structure assuming 3 ancestral populations. Ancestry was analyzed in regression models along with clinical variables. Genome wide admixture mapping was performed to identify regions of the genome where the ancestry proportions differed significantly between cases and controls. We used AncestryMap for a log-additive case control study assuming various risk models.Results: Five hundred fifty four subjects (77% of cohort) had a transfusion history with 153 (27.6 %) with allo-antibodies, 336 (60.6 %) with a history of transfusions and no antibodies and 65 (11.7 %) with no prior transfusions. Sixty five percent had a history of more than 1 antibody. The RH, Kell, MNS and Duffy systems were most commonly associated with antibodies (Fig 1A and B). Clinical factors associated with cases in univariate analysis included transfusion burden (P<0.0001), pain hospitalizations (P=0.0007), hemoglobin (P=0.02), ferritin (P=0.002), alkaline phosphatase (P=0.02) and ESR (P=0.05). Only transfusion history remained significant in multivariable models, both when using logistic (antibody status, β=0.56, P=6.31 x 10-6) or linear (number of antibodies, β=0.35, P=0.0004) regressions. Given the univariate associations with pain hospitalizations and ferritin which are known mortality risk factors, it was not surprising that allo-immunization was associated with mortality when compared to the no antibody and non-transfusion groups (Fig 1C; median follow-up 3.8 years, range 0-14.4). The absence of allo-immunization showed decreased risk of death when limited to only those with prior transfusions (Fig 1D; Hazard Ratio 0.53, P=0.006). When considering genotype data, Asian/Amerindian ancestry was associated with antibody status (logistic, β=4.60, P=0.03) and number of antibodies (linear, β=4.92, P=0.006) by regression, while African and European ancestry were not. The relative ancestry contribution remained unchanged when the number of prior transfusions was also included in both models, P=0.09 and P=0.06, respectively (Fig 1E). Genetic data also indicated that 21% (n=77) of genotyped subjects were not admixed (100% African ancestry) which violates AncestryMap’s assumption that all subjects are admixed. These subjects were removed from the admixture scan: an analysis of 81 cases and 159 controls with prior transfusions is underway to identify specific susceptibility loci.Conclusions: In a large cohort of adults with SCD, the prevalence of red cell alloimmunization is 27.6%, with approximately 60% of antibodies attributable to RH and Kell. Asian/Amerindian ancestry is a risk factor, as is the number of prior transfusions. This study demonstrates genetic susceptibility to this complication which is associated with an ancestral population not typically considered in SCD. Identification of additional Amerindian ancestry SNPs will help to fine map susceptibility loci. Long term studies will be necessary to determine the relative contribution of alloimmunization to mortality in adults with SCD. [Display omitted] DisclosuresKato:Mast Therapeutics: Consultancy; Bayer: Research Funding.

Integrative genomic analysis identifies ancestry-related expression quantitative trait loci on DNA polymerase β and supports the association of genetic ancestry with survival disparities in head and neck squamous cell carcinoma.

African Americans with head and neck squamous cell carcinoma (HNSCC) have a lower survival rate than whites. This study investigated the functional importance of ancestry-informative single-nucleotide polymorphisms (SNPs) in HNSCC and also examined the effect of functionally important genetic elements on racial disparities in HNSCC survival. Ancestry-informative SNPs, RNA sequencing, methylation, and copy number variation data for 316 oral cavity and laryngeal cancer patients were analyzed across 178 DNA repair genes. The results of expression quantitative trait locus (eQTL) analyses were also replicated with a Gene Expression Omnibus (GEO) data set. The effects of eQTLs on overall survival (OS) and disease-free survival (DFS) were evaluated. Five ancestry-related SNPs were identified as cis-eQTLs in the DNA polymerase β (POLB) gene (false discovery rate [FDR] < 0.01). The homozygous/heterozygous genotypes containing the African allele showed higher POLB expression than the homozygous white allele genotype (P < .001). A replication study using a GEO data set validated all 5 eQTLs and also showed a statistically significant difference in POLB expression based on genetic ancestry (P = .002). An association was observed between these eQTLs and OS (P < .037; FDR < 0.0363) as well as DFS (P = .018 to .0629; FDR < 0.079) for oral cavity and laryngeal cancer patients treated with platinum-based chemotherapy and/or radiotherapy. Genotypes containing the African allele were associated with poor OS/DFS in comparison with homozygous genotypes harboring the white allele. Analyses show that ancestry-related alleles could act as eQTLs in HNSCC and support the association of ancestry-related genetic factors with survival disparities in patients diagnosed with oral cavity and laryngeal cancer. Cancer 2017;123:849-60. © 2016 American Cancer Society.

Distinguishing contemporary hybridization from past introgression with postgenomic ancestry-informative SNPs in strongly differentiated Ciona species.

Biological introductions bring into contact species that can still hybridize. The evolutionary outcomes of such secondary contacts may be diverse (e.g. adaptive introgression from or into the introduced species) but are not yet well examined in the wild. The recent secondary contact between the non-native sea squirt Ciona robusta (formerly known as C.intestinalis type A) and its native congener C.intestinalis (formerly known as C.intestinalis type B), in the Western English Channel, provides an excellent case study to examine. To examine contemporary hybridization between the two species, we developed a panel of 310 ancestry-informative SNPs from a population transcriptomic study. Hybridization rates were examined on 449 individuals sampled in eight sites from the sympatric range and five sites from allopatric ranges. The results clearly showed an almost complete absence of contemporary hybridization between the two species in syntopic localities, with only one-first-generation hybrid and no other genotype compatible with recent backcrosses. Despite the almost lack of contemporary hybridization, shared polymorphisms were observed in sympatric and allopatric populations of both species. Furthermore, one allopatric population from SE Pacific exhibited a higher rate of shared polymorphisms compared to all other C. robusta populations. Altogether, these results indicate that the observed level of shared polymorphism is more probably the outcome of ancient gene flow spread afterwards at a worldwide scale. They also emphasize efficient reproductive barriers preventing hybridization between introduced and native species, which suggests hybridization should not impede too much the expansion and the establishment of the non-native species in its introduction range.

More comprehensive forensic genetic marker analyses for accurate human remains identification using massively parallel DNA sequencing.

BackgroundAlthough the primary objective of forensic DNA analyses of unidentified human remains is positive identification, cases involving historical or archaeological skeletal remains often lack reference samples for comparison. Massively parallel sequencing (MPS) offers an opportunity to provide biometric data in such cases, and these cases provide valuable data on the feasibility of applying MPS for characterization of modern forensic casework samples. In this study, MPS was used to characterize 140-year-old human skeletal remains discovered at a historical site in Deadwood, South Dakota, United States. The remains were in an unmarked grave and there were no records or other metadata available regarding the identity of the individual. Due to the high throughput of MPS, a variety of biometric markers could be typed using a single sample.ResultsUsing MPS and suitable forensic genetic markers, more relevant information could be obtained from a limited quantity and quality sample. Results were obtained for 25/26 Y-STRs, 34/34 Y SNPs, 166/166 ancestry-informative SNPs, 24/24 phenotype-informative SNPs, 102/102 human identity SNPs, 27/29 autosomal STRs (plus amelogenin), and 4/8 X-STRs (as well as ten regions of mtDNA). The Y-chromosome (Y-STR, Y-SNP) and mtDNA profiles of the unidentified skeletal remains are consistent with the R1b and H1 haplogroups, respectively. Both of these haplogroups are the most common haplogroups in Western Europe. Ancestry-informative SNP analysis also supported European ancestry. The genetic results are consistent with anthropological findings that the remains belong to a male of European ancestry (Caucasian). Phenotype-informative SNP data provided strong support that the individual had light red hair and brown eyes.ConclusionsThis study is among the first to genetically characterize historical human remains with forensic genetic marker kits specifically designed for MPS. The outcome demonstrates that substantially more genetic information can be obtained from the same initial quantities of DNA as that of current CE-based analyses.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3087-2) contains supplementary material, which is available to authorized users.

Evaluating a subset of ancestry informative SNPs for discriminating among Southwest Asian and circum-Mediterranean populations

Many different published sets of single nucleotide polymorphisms (SNPs) and/or insertion-deletion polymorphisms (InDels) can serve as ancestry informative markers (AIMs) to distinguish among continental regions of the world. For a focus on Southwest Asian ancestry we chose to start with the Kidd Lab panel of 55 ancestry-informative SNPs (AISNPs) because it already provided good global reference data (FROG-kb: frog.med.yale.edu) in a set of 73 population samples distinguishing at least 8 biogeographic clusters of populations. This panel serves as a good first tier ancestry panel. We are now interested in identifying region-specific second tier panels for more refined distinction among populations within each of the global regions. We have begun studying the global region centered on Southwest Asia and the region encompassing the Mediterranean Sea. We have incorporated 10 populations from North Africa, Turkey and Iran and included 31 of the original 73 populations and eleven 1000 Genomes Phase3 populations for a total of 3129 individuals from 52 populations, all typed for the 55 AISNPs. We have then identified the subset of the 55 AISNPs that are most informative for this region of the world using Heatmap, Fst, and Informativeness analyses to eliminate those SNPs essentially redundant or providing no information among populations in this region, reducing the number of SNPs to 32. STRUCTURE and PCA analyses show the remaining 32 SNPs identify the North African cluster and appropriately include the Turkish and Iranian samples with the Southwest Asian cluster. These markers provide the basis for building an improved, optimized panel of AISNPs that provides additional information on differences among populations in this part of the world. The data have also allowed an examination of the accuracy of the ancestry inference based on 32 SNPs for the newly studied populations from this region. The likelihood ratio approach to ancestry inference embodied in FROG-kb provides highly significant population assignments within one order of magnitude for each individual in the Turkish, Iranian, and Tunisian populations.

Inter-laboratory evaluation of the EUROFORGEN Global ancestry-informative SNP panel by massively parallel sequencing using the Ion PGM™

The EUROFORGEN Global ancestry-informative SNP (AIM-SNPs) panel is a forensic multiplex of 128 markers designed to differentiate an individual’s ancestry from amongst the five continental population groups of Africa, Europe, East Asia, Native America, and Oceania. A custom multiplex of AmpliSeq™ PCR primers was designed for the Global AIM-SNPs to perform massively parallel sequencing using the Ion PGM™ system. This study assessed individual SNP genotyping precision using the Ion PGM™, the forensic sensitivity of the multiplex using dilution series, degraded DNA plus simple mixtures, and the ancestry differentiation power of the final panel design, which required substitution of three original ancestry-informative SNPs with alternatives. Fourteen populations that had not been previously analyzed were genotyped using the custom multiplex and these studies allowed assessment of genotyping performance by comparison of data across five laboratories. Results indicate a low level of genotyping error can still occur from sequence misalignment caused by homopolymeric tracts close to the target ignment caused by homopolymeric tracts close to the target SNP, despite careful scrutiny of candidate SNPs at the design stage. Such sequence misalignment required the exclusion of component SNP rs2080161 from the Global AIM-SNPs panel. However, the overall genotyping precision and sensitivity of this custom multiplex indicates the Ion PGM™ assay for the Global AIM-SNPs is highly suitable for forensic ancestry analysis with massively parallel sequencing.