Investigating the correlation between genotypes and hepatic protein expression of CYP2C9, CYP2C19, CYP2D6, and CYP3A5 using postmortem tissue from a Danish population.

The cytochrome (CYP) P450 family of enzymes plays a central role in the metabolism of many drugs. CYP genes are highly polymorphic, which is known to affect protein levels, but for some low frequent CYP genotypes the correlation between genotype and CYP protein expression is less established. In this study, we determined the CYP2C9, CYP2C19, CYP2D6, and CYP3A5 genotypes of 250 Danish individuals included in a postmortem study. For 116 of the individuals, the hepatic CYP protein levels were investigated by a proteomics approach. Overall, we found the postmortem genetic and proteomic data to be in agreement with those of other studies performed on fresh hepatic tissue, showing the usability of postmortem hepatic tissue for this type of investigation. For less investigated genotypes we could corroborate previously found results: 1) statistically significantly lower levels of hepatic CYP2C9 protein in individuals carrying the CYP2C9*3 variant compared to individuals with two wild type (wt) alleles, 2) comparable levels of CYP2C19 in CYP2C19*2/*17 and CYP2C19*1/*2 individuals, 3) reduced CYP2D6 protein levels in heterozygous individuals with the CYP2D6*3, CYP2D6*4, and CYP2D6*5 gene deletion variants, and 4) significantly lower levels of CYP3A5 protein in CYP3A5*3 homozygous individuals compared to individuals that were heterozygous for the CYP3A5*3 allele or homozygous individuals for the wt alleles. In conclusion, the use of postmortem tissue significantly increases the access to human specimens for research purposes, and postmortem proteomics can be used to investigate the link between CYP genotypes and hepatic protein expression. Significance Statement In tissue from a large postmortem cohort (n=250) we determined the CYP2C9, CYP2C19, CYP2D6, and CYP3A5 genotypes. Hepatic CYP protein levels were investigated in 116 individuals using a proteomics approach. For common genotypes, we found results similar to previous knowledge, pointing towards the usability of postmortem tissue. For the less investigated genotypes, we were able to corroborate genotype / protein expression correlations. It is a novel approach to use a large postmortem cohort to investigate genetic / protein expression correlations.

Cleaning protocols in forensic genetic laboratories.

It is pivotal to avoid cross-sample contamination in forensic genetic laboratories and optimal cleaning protocols for the removal of DNA are essential. A survey was performed, and ten forensic genetic laboratories shared their cleaning protocols in pre-PCR and post-PCR laboratories. The cleaning frequencies on different surface areas were somewhat similar, whereas none of the laboratories used the same cleaning reagents. Therefore, the efficiencies of the cleaning protocol utilised were tested and compared. The results showed that freshly made household bleach and Virkon® removed all amplifiable DNA from the surfaces, whereas DNA AWAY™ and the disinfection reagents ethanol, isopropanol, and ChemGene HLD4L did not.

Development of an Okinawa panel for biogeographic inference of Okinawans.

BACKGROUND: The Precision ID Ancestry Panel with 165 SNP markers was unable to differentiate between mainland Japanese and Okinawa Japanese or to distinguish either of them from other East Asian populations. AIM: An Okinawa panel was developed with the aim of further separating Okinawa Japanese individuals from mainland Japanese and other Asian groups. Seventy-five SNPs were selected using the most informative markers from the literature. Further, 22 SNPs were selected to separate Okinawa Japanese at minimum SNPs. SUBJECTS AND METHODS: Samples were collected from 48 unrelated individuals from mainland Japan and 46 unrelated residents of the Okinawa prefecture. Data were evaluated by STRUCTURE, principal component, and GenoGeographer analyses. RESULTS: The 22 SNP set had similar levels of differentiation in STRUCTURE and PCA analyses as the 75 SNP set. GenoGeographer analysis showed that, out of the 46 Okinawa Japanese individuals, the 75 SNP and 22 SNP sets correctly assigned the Okinawan population as the most likely population of origin for 32 and 31 individuals, respectively. CONCLUSION: Neither SNP set could completely differentiate between Okinawa Japanese and other Asian groups, however, these sets should be useful for crime investigation, when the sample, cost and time are limited.

Gene expressions in cerebral palsy subjects reveal structural and functional changes in the gastrocnemius muscle that are closely associated with passive muscle stiffness.

Cerebral palsy (CP) is a non-progressive motor disorder that affects posture and gait due to contracture development. The purpose of this study is to analyze a possible relation between muscle stiffness and gene expression levels in muscle tissue of children with CP. Next-generation sequencing (NGS) of gene transcripts was carried out in muscle biopsies from gastrocnemius muscle (n = 13 children with CP and n = 13 typical developed (TD) children). Passive stiffness of the ankle plantarflexors was measured. Structural changes of the basement membranes and the sarcomere length were measured. Twelve pre-defined gene target sub-categories of muscle function, structure and metabolism showed significant differences between muscle tissue of CP and TD children. Passive stiffness was significantly correlated to gene expression levels of HSPG2 (p = 0.02; R2 = 0.67), PRELP (p = 0.002; R2 = 0.84), RYR3 (p = 0.04; R2 = 0.66), C COL5A3 (p = 0.0007; R2 = 0.88), ASPH (p = 0.002; R2 = 0.82) and COL4A6 (p = 0.03; R2 = 0.97). Morphological differences in the basement membrane were observed between children with CP and TD children. The sarcomere length was significantly increased in children with CP when compared with TD (p = 0.04). These findings show that gene targets in the categories: calcium handling, basement membrane and collagens, were significantly correlated to passive muscle stiffness. A Reactome pathway analysis showed that pathways involved in DNA repair, ECM proteoglycans and ion homeostasis were amongst the most upregulated pathways in CP, while pathways involved in collagen fibril crosslinking, collagen fibril assembly and collagen turnover were amongst the most downregulated pathways when compared with TD children. These results underline that contracture formation and motor impairment in CP is an interplay between multiple factors.

Forensic application and genetic diversity of 21 autosomal STR loci in five major population groups of Pakistan.

OBJECTIVES: Investigation of genetic diversity of the 21 autosomal STR loci included in the GlobalFilerTM PCR Amplification Kit in 529 Pakistani individuals belonging to the Punjabi, Pashtun, Sindhi, Saraiki, and Baloch ethnic groups. Population genetic parameters and forensic informative metrics for each group were evaluated. RESULTS: SE33 showed the greatest power of discrimination in all populations studied. The combined match probability ranged from 8.06E-27 (Saraiki) to 1.05E-26 (Baloch), and the combined power of exclusion ranged from 0.99999999902 (Punjabi) to 0.99999999964 (Pashtun). D12S391 in the Baloch population and D2S441 in the Saraiki population showed deviation from Hardy-Weinberg equilibrium. CONCLUSION: Significant genetic distances were observed between the Punjabi, Pashtun, and Baloch populations. This study supports the utilization of the GlobalFilerTM STR kit for forensic applications in Pakistan.

Genetic analysis of sixteen autosomal STR loci in three Tunisian populations from Makthar, Nabeul and Sousse.

BACKGROUND: Due to its strategic location, Tunisia witnessed the succession and influence of many civilisations throughout history. However, the majority of studies carried out on Tunisia are focussed on Barbarian ethnicity. AIM: To estimate genetic diversity and genetic structure of three Tunisian populations using autosomal STRs. SUBJECTS AND METHODS: 278 individuals were analysed for sixteen STRs. Allele frequencies and statistical parameters were determined. RESULTS: The studied populations showed genetic affinity with geographically close populations. AMOVA showed no genetic difference between the Tunisian populations. Nevertheless, the variance between the populations of the same group was significant, reflecting their heterogeneity even though they came from the same geographical area, and had the same ethnicity and complex demographic history. CONCLUSION: Our results strongly supported the application of autosomal genetic markers in anthropological and forensic studies. The analyses conducted at the 15-loci level provide the resolution to assess the phylogenetic relationships among the populations examined and other geographically targeted worldwide populations, while the results resulting from the 10-loci studies provide an understanding of the relationships and origins of the North African populations. Furthermore, the current report demonstrates that the battery of autosomal STRs reported are useful, providing the power of discrimination for forensic and paternity analyses.

Skin pigmentation and genetic variants in an admixed Brazilian population of primarily European ancestry.

Although many genes have been shown to be associated with human pigmentary traits and forensic prediction assays exist (e.g. HIrisPlex-S), the genetic knowledge about skin colour remains incomplete. The highly admixed Brazilian population is an interesting study population for investigation of the complex genotype-phenotype architecture of human skin colour because of its large variation. Here, we compared variants in 22 pigmentary genes with quantitative skin pigmentation levels on the buttock, arm, and forehead areas of 266 genetically admixed Brazilian individuals. The genetic ancestry of each individual was estimated by typing 46 AIM-InDels. The mean proportion of genetic ancestry was 68.8% European, 20.8% Sub-Saharan African, and 10.4% Native American. A high correlation (adjusted R2 = 0.65, p < 0.05) was observed between nine SNPs and quantitative skin pigmentation using multiple linear regression analysis. The correlations were notably smaller between skin pigmentation and biogeographic ancestry (adjusted R2 = 0.45, p < 0.05), or markers in the leading forensic skin colour prediction system, the HIrisPlex-S (adjusted R2 = 0.54, p < 0.05). Four of the nine SNPs, OCA2 rs1448484 (rank 2), APBA2 rs4424881 (rank 4), MFSD12 rs10424065 (rank 8), and TYRP1 1408799 (rank 9) were not investigated as part of the HIrisPlex-S selection process, and therefore not included in the HIrisPlex-S model. Our results indicate that these SNPs account for a substantial part of the skin colour variation in individuals of admixed ancestry. Hence, we suggest that these SNPs are considered when developing future skin colour prediction models.

Non-invasive prenatal paternity testing using a standard forensic genetic massively parallel sequencing assay for amplification of human identification SNPs.

Prenatal paternity testing often relies on invasive procedures that cause risk to both the mother and the foetus. Non-invasive, prenatal paternity testing by investigating paternally inherited single nucleotide polymorphisms (SNPs) in cell-free foetal DNA (cffDNA) in maternal plasma was performed at consecutive time points during early gestation. Plasma from 15 pregnant women was investigated at consecutive time points from gestational weeks (GWs) 4-20. The Precision ID Identity Panel and an Ion S5 Sequencer was used to analyse the cffDNA. Paternally inherited foetal SNP alleles were detected from GW7. The median foetal fractions were 0%, 3.9%, 5.1%, 5.2%, and 4.7% at GWs 4, 7, 12, 16, and 20, respectively. The corresponding median numbers of detected paternally inherited foetal autosomal SNP alleles were 0, 3, 9, 10, and 12, respectively. The typical (i.e. geometric mean) paternity indices at GW12 and GW20 were 24 (range 0.0035-8389) and 199 (range 5.1-30,137), respectively. The method is very promising. However, the method can be improved by shortening the lengths of the PCR amplicons and increasing the number of SNPs. To our knowledge, this is the first study to successfully identify paternally inherited foetal SNP alleles at consecutive time points in early gestation independently of the foetal gender.

Sequence variants in muscle tissue-related genes may determine the severity of muscle contractures in cerebral palsy.

Muscle contractures are a common complication to cerebral palsy (CP). The purpose of this study was to evaluate whether individuals with CP carry specific gene variants of important structural genes that might explain the severity of muscle contractures. Next-generation-sequencing (NGS) of 96 candidate genes associated with muscle structure and metabolism were analyzed in 43 individuals with CP (Gross Motor Function classification system [GMFCS] I, n=10; GMFCS II, n=14; GMFCS III, n=19) and four control participants. In silico analysis of the identified variants was performed. The variants were classified into four categories ranging from likely benign (VUS0) to highly likely functional effect (VUS3). All individuals with CP were classified and grouped according to their GMFCS level: Statistical comparisons were made between GMFCS groups. Kruskal-Wallis tests showed significantly more VUS2 variants in the genes COL4 (GMFCS I-III; 1, 1, 5, respectively [p < .04]), COL5 (GMFCS I-III; 1, 1, 5 [p < .04]), COL6 (GMFCS I-III; 0, 4, 7 [p < .003]), and COL9 (GMFCS I-III; 1, 1, 5 [p < .04]), in individuals with CP within GMFCS Level III when compared to the other GMFCS levels. Furthermore, significantly more VUS3 variants in COL6 (GMFCS I-III; 0, 5, 2 [p < .01]) and COL7 (GMFCS I-III; 0, 3, 0 [p < .04]) were identified in the GMFCS II level when compared to the other GMFCS levels. The present results highlight several candidate gene variants in different collagen types with likely functional effects in individuals with CP.

Sequencing of mitochondrial genomes using the Precision ID mtDNA Whole Genome Panel.

Massively parallel sequencing offers a fast and cost-effective method for sequencing of the whole mtDNA genome. The Precision ID mtDNA Whole Genome Panel amplifies the entire mtDNA genome in two multiplex PCRs with 81 primer sets using the Ion AmpliSeq™ technology. In this study, the performance of the panel was evaluated by testing different amplification methods (two-in-one or conservative), the number of PCR cycles (21, 23, and 25), and different reaction volumes (recommended volume or half-volume). Furthermore, a dilution series, controlled mtDNA mixtures, and casework samples were also sequenced. The normalised read depths of the individual fragments were highly consistent irrespectively of the amplification methods or reaction volumes used. The sequencing output matched the mixture ratios of the controlled mtDNA mixtures indicating that the sequencing results were a loyal representation of the input DNA. Complete mtDNA profiles were generated from <10 pg genomic DNA. Seven fragments with relatively low read depths and large variations in read depth were identified. One of these fragments covered part of the control region and the hypervariable region II.

ISO 17025 validation of a next-generation sequencing assay for relationship testing.

The HID-Ion AmpliSeq™ Identity Panel is a next-generation sequencing assay with 90 autosomal and 34 Y-chromosome SNPs that are amplified in one PCR step and subsequently sequenced using the Ion Personal Genome Machine (Ion PGM™) System. This assay was validated for relationship testing in our ISO 17025 accredited laboratory in 2015. Here, the essential parts of the validation report submitted to the Danish Accreditation Fund are presented. A total of 100 unrelated Danes were typed in duplicates and the locus balance, heterozygote balance (Hb) and noise levels were analysed in detail. Two loci were disregarded for casework because genotyping was uncertain. Hb for rs7520386 was skewed and high levels of noise were observed in rs576261. Three general acceptance criteria for analysis of single-source samples were defined: (i) sequencing depth > 200 reads, (ii) noise level < 3% and (iii) Hb > 0.3. A Python script named SNPonPGM was developed to assist the analyst by highlighting loci that do not fulfil the general acceptance criteria. Furthermore, SNPonPGM has functions that reduce the hands-on time of the reporting officer to a few minutes per case. Mixtures with DNA from two individuals in a 1:24 ratio were readily identified using the three criteria and the SNPonPGM script.

A study of the peopling of Greenland using next generation sequencing of complete mitochondrial genomes.

OBJECTIVES: The Greenlandic population history is characterized by a number of migrations of people of various ethnicities. In this work, the analysis of the complete mtDNA genome aimed to contribute to the ongoing debate on the origin of current Greenlanders and, at the same time, to address the migration patterns in the Greenlandic population from a female inheritance demographic perspective. METHODS: We investigated the maternal genetic variation in the Greenlandic population by sequencing the whole mtDNA genome in 127 Greenlandic individuals using the Illumina MiSeq® platform. RESULTS: All Greenlandic individuals belonged to the Inuit mtDNA lineages A2a, A2b1, and D4b1a2a1. No European haplogroup was found. DISCUSSION: The mtDNA lineages seem to support the hypothesis that the Inuit in Greenland are descendants from the Thule migration. The results also reinforce the importance of isolation and genetic drift in shaping the genetic diversity in Greenlanders. Based on the mtDNA sequences, the Greenlandic Inuit are phylogenetically close to Siberian groups and Canadian Inuit.

Performance of a 74-Microhaplotype Assay in Kinship Analyses.

Microhaplotypes (MHs) consisting of multiple SNPs and indels on short stretches of DNA are new and interesting loci for forensic genetic investigations. In this study, we analysed 74 previously defined MHs in two of the populations that our laboratory provides with forensic genetic services, Danes and Greenlanders. In addition to the 229 SNPs that originally made up the 74 MHs, 66 SNPs and 3 indels were identified in the two populations, and 45 of these variants were included in new definitions of the MHs, whereas 24 SNPs were considered rare and of little value for case work. The average effective number of alleles (Ae) was 3.2, 3.0, and 2.6 in Danes, West Greenlanders, and East Greenlanders, respectively. High levels of linkage disequilibrium were observed in East Greenlanders, which reflects the characteristics of this population that has a small size, and signs of admixture and substructure. Pairwise kinship simulations of full siblings, half-siblings, first cousins, and unrelated individuals were performed using allele frequencies from MHs, STRs and SNPs from Danish and Greenlandic populations. The MH panel outperformed the currently used STR and SNP marker sets and was able to differentiate siblings from unrelated individuals with a 0% false positive rate and a 1.1% false negative rate using an LR threshold of 10,000 in the Danish population. However, the panel was not able to differentiate half-siblings or first cousins from unrelated individuals. The results generated in this study will be used to implement MHs as investigative markers for relationship testing in our laboratory.

Identification of individuals from low template blood samples using whole transcriptome shotgun sequencing.

Biological trace samples consisting of very few cells pose a challenge to conventional forensic genetic DNA analysis. RNA may be an alternative to DNA when handling low template samples. Whereas each cell only contains two copies of an autosomal DNA segment, the transcriptome retains much of the genomic variation replicated in abundant RNA fragments. In this study, we describe the development of a prototype RNA-based SNP selection set for forensic human identification from low template samples (50 pg gDNA). Whole blood from a subset of the Danish population (41 individuals) and blood stains subjected to degradation at room temperature for up to two weeks were analysed by whole transcriptome shotgun sequencing. Concordance was determined by DNA genotyping with the Infinium Omni5-4 SNP chip. In the 100 protein-coding genes with the most reads, 5214 bi-allelic SNPs with gnomAD minor allele frequencies > 0.1 in the African/African American, East Asian, and (non-Finnish) European populations were identified. Of these, 24 SNPs in 21 genes passed screening in whole blood and degraded blood stains, with a resulting mean match probability of 4.5 ∙ 10-9. Additionally, ancestry informative SNPs and SNPs in genes useful for body fluid identification were identified in the transcriptome. Consequently, shotgun sequencing of RNA from low template samples may be used for a vast host of forensic genetics purposes, including simultaneous human and body fluid identification, leading to direct donor identification in the identified body fluid.

Collaborative exercise: analysis of age estimation using a QIAGEN protocol and the PyroMark Q48 platform.

Human age estimation from trace samples may give important leads early in a police investigation by contributing to the description of the perpetrator. Several molecular biomarkers are available for the estimation of chronological age, and currently, DNA methylation patterns are the most promising. In this study, a QIAGEN age protocol for age estimation was tested by five forensic genetic laboratories. The assay comprised bisulfite treatment of the extracted DNA, amplification of five CpG loci (in the genes of ELOVL2, C1orf132, TRIM59, KLF14, and FHL2), and sequencing of the amplicons using the PyroMark Q48 platform. Blood samples from 49 individuals with ages ranging from 18 to 64 years as well as negative and methylation controls were analyzed. An existing age estimation model was applied to display a mean absolute deviation of 3.62 years within the reference data set. Key points: Age determination as an intelligence tool during investigations can be a powerful tool in forensic genetics.In this study, five laboratories ran 49 samples and obtained a mean absolute deviation of 3.62 years.Five markers were analyzed on a PyroMark Q48 platform.

SNPtotree-Resolving the Phylogeny of SNPs on Non-Recombining DNA.

Genetic variants on non-recombining DNA and the hierarchical order in which they accumulate are commonly of interest. This variant hierarchy can be established and combined with information on the population and geographic origin of the individuals carrying the variants to find population structures and infer migration patterns. Further, individuals can be assigned to the characterized populations, which is relevant in forensic genetics, genetic genealogy, and epidemiologic studies. However, there is currently no straightforward method to obtain such a variant hierarchy. Here, we introduce the software SNPtotree v1.0, which uniquely determines the hierarchical order of variants on non-recombining DNA without error-prone manual sorting. The algorithm uses pairwise variant comparisons to infer their relationships and integrates the combined information into a phylogenetic tree. Variants that have contradictory pairwise relationships or ambiguous positions in the tree are removed by the software. When benchmarked using two human Y-chromosomal massively parallel sequencing datasets, SNPtotree outperforms traditional methods in the accuracy of phylogenetic trees for sequencing data with high amounts of missing information. The phylogenetic trees of variants created using SNPtotree can be used to establish and maintain publicly available phylogeny databases to further explore genetic epidemiology and genealogy, as well as population and forensic genetics.

Prediction of chronological age and its applications in forensic casework: methods, current practices, and future perspectives.

Estimating an individual's age can be relevant in several areas primarily related to the clinical and forensic fields. In the latter, estimation of an individual's chronological age from biological material left by the perpetrator at a crime scene may provide helpful information for police investigation. Estimation of age is also beneficial in immigration cases, where age can affect the person's protection status under the law, or in disaster victim identification to narrow the list of potential missing persons. In the last decade, research has focused on establishing new approaches for age prediction in the forensic field. From the first forensic age estimations based on morphological inspections of macroscopic changes in bone and teeth, the focus has shifted to molecular methods for age estimation. These methods allow the use of samples from human biological material that does not contain morphological age features and can, in theory, be investigated in traces containing only small amounts of biological material. Molecular methods involving DNA analyses are the primary choice and estimation of DNA methylation levels at specific sites in the genome is the most promising tool. This review aims to provide an overview of the status of forensic age prediction using molecular methods, with particular focus in DNA methylation. The frequent challenges that impact forensic age prediction model development will be addressed, together with the importance of validation efforts within the forensic community.

Comparison of global DNA methylation analysis by whole genome bisulfite sequencing and the Infinium Mouse Methylation BeadChip using fresh and fresh-frozen mouse epidermis.

Until recently, studying the murine methylome was restricted to sequencing-based methods. In this study we compared the global DNA methylation levels of hairless mouse epidermis using the recently released Infinium Mouse Methylation BeadChip from Illumina and whole genome bisulphite sequencing (WGBS). We also studied the effect of sample storage conditions by using fresh and fresh-frozen epidermis. The DNA methylation levels of 123,851 CpG sites covered by both the BeadChip and WGBS were compared. DNA methylation levels obtained with WGBS and the BeadChip were strongly correlated (Pearson correlation r = 0.984). We applied a threshold of 15 reads for the WGBS methylation analysis. Even at a threshold of 10 reads, we observed no substantial difference in DNA methylation levels compared with that obtained with the BeadChip. The DNA methylation levels from the fresh and the fresh-frozen samples were strongly correlated when analysed with both the BeadChip (r = 0.999) and WGBS (r = 0.994). We conclude that the two methods of analysis generally work equally well for studies of DNA methylation of mouse epidermis and find that fresh and fresh-frozen epidermis can generally be used equally well. The choice of method will depend on the specific study's aims and the available resources in the laboratory.

Pitfalls and challenges with population assignments of individuals from admixed populations: Applying Genogeographer on Brazilian individuals.

The assignment of individuals to a population can be of importance for the identification of mass disaster victims or criminal offenders in the field of forensic genetics. This assignment is based on biostatistical methods that process data of ancestry informative markers (AIMs), which are selected based on large allele frequency differences between the populations of interest. However, population assignments of individuals with an admixed genetic background are challenging. Admixed individuals are genetic mosaics of chromosomal segments from the parental populations, which may lead to ambiguous or no population assignment. This is problematic since admixture events are a substantial part of human history. In this study, we present challenges of interpreting the evidential weight of population assignments. We used Genogeographer for likelihood ratio (LR) calculations and Brazilians as examples of admixed individuals. Brazilians are a very heterogenous population representing a three-way admixture between Native Americans, Europeans, and Africans. Ancestry informative markers were typed in a total of 589 individuals from Brazil using the Precision ID Ancestry Panel. The Brazilians were assigned to six metapopulations (East Asia, Europe, Middle East, North Africa, South-Central Asia, Sub-Saharan Africa) defined in the Genogeographer software and LRs were calculated if the AIM profile was not an outlier in all metapopulations and simulated two-way (1:1) admixtures of the six metapopulations. Population assignments failed for 55% of the samples. These samples had significantly higher genetic contributions from East Asia, South-Central Asia and Sub-Saharan Africa, and significantly lower genetic contributions from Europe. Most of the individuals with population assignments were assigned to the metapopulations of Middle East (58%) or North Africa (36%), followed by Europe (4%), South-Central Asia (1%), and Sub-Saharan Africa (1%). For 8% of the samples, population assignments were only possible when assignments to simulated two-way (1:1) admixtures of the six metapopulations were considered. Most of these individuals were assigned to two-way admixtures of North Africa, South-Central Asia, or Sub-Saharan Africa. Relatively low median likelihood ratios (LRs<1000) were observed when comparing population likelihoods for Europe, Middle East, North Africa, South-Central Asia, or simulated 1:1 admixtures of these metapopulations. Comparisons including East Asian or Sub-Saharan African populations resulted in larger median LRs (LR>1010). The results suggested that the Precision ID Ancestry Panel provided too little information and that additional markers specifically selected for sub-continental differentiation may be required for accurate population assignment of admixed individuals. Furthermore, a Genogeographer database with additional populations including admixed populations would be advantageous for interpretation of admixed AIM profiles. It would likely increase the number of population assignments and illustrate alternatives to the most likely population, which would be valuable information for the case officer when writing the case report.

Association between Variants in the OCA2-HERC2 Region and Blue Eye Colour in HERC2 rs12913832 AA and AG Individuals.

The OCA2-HERC2 region is strongly associated with human pigmentation, especially eye colour. The HERC2 SNP rs12913832 is currently the best-known predictor for blue and brown eye colour. However, in a previous study we found that 43 of 166 Norwegians with the brown eye colour genotype rs12913832:AA or AG, did not have the expected brown eye colour. In this study, we carried out massively parallel sequencing of a ~500 kbp HERC2-OCA2 region in 94 rs12913832:AA and AG Norwegians (43 blue-eyed and 51 brown-eyed) to search for novel blue eye colour variants. The new candidate variants were subsequently typed in a Norwegian biobank population (total n = 519) for population specific association analysis. We identified five new variants, rs74409036:A, rs78544415:T, rs72714116:T, rs191109490:C and rs551217952:C, to be the most promising candidates for explaining blue eye colour in individuals with the rs12913832:AA and AG genotype. Additionally, we confirmed the association of the missense variants rs74653330:T and rs121918166:T with blue eye colour, and observed lighter skin colour in rs74653330:T individuals. In total, 37 (86%) of the 43 blue-eyed rs12913832:AA and AG Norwegians could potentially be explained by these seven variants, and we suggest including them in future prediction models.