Changes in mouse epidermal DNA methylation during development of squamous cell carcinoma in response to UVR.

Squamous cell carcinoma (SCC) is a common skin cancer, often caused by exposure to ultraviolet radiation (UVR). Recent studies have shown that changes in DNA methylation play a crucial role in the development of cancers. However, methylation patterns of SCC are not well characterised. Identifying biomarkers for the risk of developing SCC could be helpful for early detection and diagnosis and can potentially improve treatment and prevention strategies. This study aimed to investigate methylation changes in the epidermis of mice exposed to UVR for 24 weeks. We examined the DNA methylation levels of 260 199 CpGs using the Illumina Infinium Mouse Methylation BeadChip and studied the epidermis of UVR-exposed and unexposed mice every 4 weeks for 24 weeks (n = 39). We identified CpGs with large differences in methylation levels (ß-values) between UVR-exposed and unexposed mice. We also observed differences in the epigenetic age of these mice. We identified CpGs in Rev, Ipmk, Rad51b, Fgfr2, Fgfr3 and Ctnnb1 that may serve as potential biomarkers for SCC risk and could be helpful for the early detection and prevention of SCC. Further investigations are necessary to determine the biological functions and clinical significance of these CpGs.

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.

Genetic investigations of 100 inherited cardiac disease-related genes in deceased individuals with schizophrenia.

Cardiac diseases and sudden cardiac death (SCD) are more prevalent in individuals diagnosed with schizophrenia compared to the general population, with especially coronary artery disease (CAD) as the major cardiovascular cause of death. Antipsychotic medications, genetics, and lifestyle factors may contribute to the increased SCD in individuals with schizophrenia. The role of antipsychotic medications and lifestyle factors have been widely investigated, while the genetic predisposition to inherited cardiac diseases in schizophrenia is poorly understood. In this study, we examined 100 genes associated with inherited cardiomyopathies and cardiac channelopathies in 97 deceased individuals diagnosed with schizophrenia for the prevalence of genetic variants associated with SCD. The deceased individuals had various causes of death and were included in the SURVIVE project, a prospective, autopsy-based study of mentally ill individuals in Denmark. This is the first study of multiple inherited cardiac disease-related genes in deceased individuals with diagnosed schizophrenia to shed light on the genetic predisposition to SCD in individuals with schizophrenia. We found no evidence for an overrepresentation of rare variants with high penetrance in inherited cardiac diseases, following the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG) consensus guidelines. However, we found that the deceased individuals had a statistically significantly increased polygenic burden caused by variants in the investigated heart genes compared to the general population. This indicates that common variants with smaller effects in heart genes may play a role in schizophrenia.

Differential Methylation in the GSTT1 Regulatory Region in Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy.

Sudden cardiac death (SCD) is a diagnostic challenge in forensic medicine. In a relatively large proportion of the SCDs, the deaths remain unexplained after autopsy. This challenge is likely caused by unknown disease mechanisms. Changes in DNA methylation have been associated with several heart diseases, but the role of DNA methylation in SCD is unknown. In this study, we investigated DNA methylation in two SCD subtypes, sudden unexplained death (SUD) and sudden unexpected death in epilepsy (SUDEP). We assessed DNA methylation of more than 850,000 positions in cardiac tissue from nine SUD and 14 SUDEP cases using the Illumina Infinium MethylationEPIC BeadChip. In total, six differently methylated regions (DMRs) between the SUD and SUDEP cases were identified. The DMRs were located in proximity to or overlapping genes encoding proteins that are a part of the glutathione S-transferase (GST) superfamily. Whole genome sequencing (WGS) showed that the DNA methylation alterations were not caused by genetic changes, while whole transcriptome sequencing (WTS) showed that DNA methylation was associated with expression levels of the GSTT1 gene. In conclusion, our results indicate that cardiac DNA methylation is similar in SUD and SUDEP, but with regional differential methylation in proximity to GST genes.

Targeted exon sequencing in deceased schizophrenia patients in Denmark.

Schizophrenia patients have higher mortality rates and lower life expectancy than the general population. However, forensic investigations of their deaths often fail to determine the cause of death, hindering prevention. As schizophrenia is a highly heritable condition and given recent advances in our understanding of the genetics of schizophrenia, it is now possible to investigate how genetic factors may contribute to mortality. We made use of findings from genome-wide association studies (GWAS) to design a targeted panel (PsychPlex) for sequencing of exons of 451 genes near index single nucleotide polymorphisms (SNPs) identified with GWAS. We sequenced the DNA of 95 deceased schizophrenia patients included in SURVIVE, a prospective, autopsy-based study of mentally ill persons in Denmark. We compared the allele frequencies of 1039 SNPs in these cases with the frequencies of 2000 Danes without psychiatric diseases and calculated their deleteriousness (CADD) scores. For 81 SNPs highly associated with schizophrenia and CADD scores above 15, expression profiles in the Genotype-Tissue Expression (GTEx) Project indicated that these variants were in exons, whose expressions are increased in several types of brain tissues, particularly in the cerebellum. Molecular pathway analysis indicated the involvement of 163 different pathways. As for rare SNP variants, most variants were scored as either benign or likely benign with an average of 17 variants of unknown significance per individual and no pathogenic variant. Our results highlight the potential of DNA sequencing of an exon panel to discover genetic factors that may be involved in the development of schizophrenia.

Pigment genes not skin pigmentation affect UVB-induced vitamin D.

Skin pigmentation is believed to contribute to the generally low serum 25-hydroxyvitamin D (25(OH)D) concentrations observed in darker-skinned persons. The influence of measured skin pigmentation on UVB-induced 25(OH)D increase was investigated together with 9 demographic and 13 genetic parameters (pigment SNPs). Forty participants representing a wide range in measured skin pigmentation were exposed to identical UVB doses on identical body areas over nine weeks with weekly measurements of serum 25(OH)D. This study took place in Denmark during winter, a period with negligible ambient UVB, so variation in 25(OH)D synthesis was not influenced by latitude, season, sun and clothing habits. The increase in 25(OH)D concentration displayed considerable variation (range: 2.9 to 139 nmol L-1). Constitutive and facultative skin pigmentation exerted separate influence on the variation of the UVB-induced linear 25(OH)D increase. However, this influence was statistically non-significant in the presence of separate significant pigment SNPs. The variation in the 25(OH)D increase in the combined linear model was not explained by measured skin pigmentation but by sex, height, age and seven SNPs located in the ASIP, MTAP, MIR196A29 and Solute Carrier Family genes. This linear model including individual intercepts and the 10 parameters influencing the slope explained 77.4% of the variation. This study confirmed the influence of sex, age and height on 25(OH)D increase and found that pigment genes provided a better relation to UVB-induced 25(OH)D increase compared to the actual measured skin pigmentation. Therefore, only investigating skin pigmentation obscures other causal parameters for low 25(OH)D.

Serum 25(OH)D levels after oral vitamin D3 supplementation and UVB exposure correlate.

BACKGROUND: The inter-individual variation in 25(OH)D3 increase (Δ25(OH)D3 ) after vitamin D3 supplementation was determined and compared with the UVB irradiation response. METHODS: Nineteen Danish participants received 85 µg vitamin D3 (cholecalciferol) daily for nine weeks with regular serum 25(OH)D3 measurements. These participants had three years earlier taken part in a 9-week controlled UVB study. The Δ25(OH)D3 was not confounded by ambient UVB, BMI or ethnicity. RESULTS: Δ25(OH)D3 was 53 nmol L-1 and almost identical to Δ25(OH)D3 (52 nmol L-1 ) after UVB. Δ25(OH)D3 ranged from 17 to 91 nmol L-1 (span 74 nmol L-1 ) and was about half of that observed after UVB irradiation (span 136 nmol L-1 ). The interquartile ranges for vitamin D3 supplementation (38.8-71.4 nmol L-1 , span: 32.6 nmol L-1 ) and UVB irradiation (35.7-65.4 nmol L-1 , span: 29.7 nmol L-1 ) were similar indicating a comparable response of the two interventions. As the 25(OH)D3 start levels (R2  = 0.398, P = 3.8 × 10-3 ), 25(OH)D3 end levels (R2  = 0.457, P = 1.5 × 10-3 ) and Δ25(OH)D3 (R2  = 0.253, P = 0.028) between both interventions were correlated, this suggested a possible common individual background for the variation. Four pigment SNPs influenced the variation in the vitamin D3 -induced and UVB-induced Δ25(OH)D3 . A combined model including the influence of these four SNPs and the 25(OH)D3 start level explained 86.8% (P = 1.6 × 10-35 ) of the individual variation after vitamin D3 supplementation. CONCLUSION: The inter-individual variation in the two interventions was comparable and had no common demographic but a partly common genetic background.

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.

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.

Shotgun DNA sequencing for human identification: Dynamic SNP selection and likelihood ratio calculations accounting for errors.

Shotgun sequencing is a DNA analysis method that potentially determines the nucleotide sequence of every DNA fragment in a sample, unlike PCR-based genotyping methods that is widely used in forensic genetics and targets predefined short tandem repeats (STRs) or predefined single nucleotide polymorphisms (SNPs). Shotgun DNA sequencing is particularly useful for highly degraded low-quality DNA samples, such as ancient samples or those from crime scenes. Here, we developed a statistical model for human identification using shotgun sequencing data and developed formulas for calculating the evidential weight as a likelihood ratio (LR). The model uses a dynamic set of binary SNP loci and takes the error rate from shotgun sequencing into consideration in a probabilistic manner. To our knowledge, the method is the first to make this possible. Results from replicated shotgun sequencing of buccal swabs (high-quality samples) and hair samples (low-quality samples) were arranged in a genotype-call confusion matrix to estimate the calling error probability by maximum likelihood and Bayesian inference. Different genotype quality filters may be applied to account for genotyping errors. An error probability of zero resulted in the commonly used LR formula for the weight of evidence. Error probabilities above zero reduced the LR contribution of matching genotypes and increased the LR in the case of a mismatch between the genotypes of the trace and the person of interest. In the latter scenario, the LR increased from zero (occurring when the error probability was zero) to low positive values, which allow for the possibility that the mismatch may be due to genotyping errors. We developed an open-source R package, wgsLR, which implements the method, including estimation of the calling error probability and calculation of LR values. The R package includes all formulas used in this paper and the functionalities to generate the formulas.

Evaluation of DNAmAge in paired fresh, frozen, and formalin-fixed paraffin-embedded heart tissues.

The continued development in methylome analysis has enabled a more precise assessment of DNA methylation, but treatment of target tissue prior to analysis may affect DNA analysis. Prediction of age based on methylation levels in the genome (DNAmAge) has gained much interest in disease predisposition (biological age estimation), but also in chronological donor age estimation in crime case samples. Various epigenetic clocks were designed to predict the age. However, it remains unknown how the storage of the tissues affects the DNAmAge estimation. In this study, we investigated the storage method impact of DNAmAge by the comparing the DNAmAge of the two commonly used storage methods, freezing and formalin-fixation and paraffin-embedding (FFPE) to DNAmAge of fresh tissue. This was carried out by comparing paired heart tissue samples of fresh tissue, samples stored by freezing and FFPE to chronological age and whole blood samples from the same individuals. Illumina EPIC beadchip array was used for methylation analysis and the DNAmAge was evaluated with the following epigenetic clocks: Horvath, Hannum, Levine, Horvath skin+blood clock (Horvath2), PedBE, Wu, BLUP, EN, and TL. We observed differences in DNAmAge among the storage conditions. FFPE samples showed a lower DNAmAge compared to that of frozen and fresh samples. Additionally, the DNAmAge of the heart tissue was lower than that of the whole blood and the chronological age. This highlights caution when evaluating DNAmAge for FFPE samples as the results were underestimated compared with fresh and frozen tissue samples. Furthermore, the study also emphasizes the need for a DNAmAge model based on heart tissue samples for an accurate age estimation.

Comparison of whole transcriptome sequencing of fresh, frozen, and formalin-fixed, paraffin-embedded cardiac tissue.

The use of fresh tissue for molecular studies is preferred but often impossible. Instead, frozen or formalin-fixed, paraffin-embedded (FFPE) tissues are widely used and constitute valuable resources for retrospective studies. We assessed the utility of cardiac tissue stored in different ways for gene expression analyses by whole transcriptome sequencing of paired fresh, frozen, and FFPE tissues. RNA extracted from FFPE was highly degraded. Sequencing of RNA from FFPE tissues yielded higher proportions of intronic and intergenic reads compared to RNA from fresh and frozen tissues. The global gene expression profiles varied with the storage conditions, particularly mitochondrial and long non-coding RNAs. However, we observed high correlations among protein-coding transcripts (ρ > 0.94) with the various storage conditions. We did not observe any significant storage effect on the allele-specific gene expression. However, FFPE had statistically significantly (p < 0.05) more discordant variant calls compared to fresh and frozen tissue. In conclusion, we found that frozen and FFPE tissues can be used for reliable gene expression analyses, provided that proper quality control is performed and caution regarding the technical variability is withheld.

Fresh and frozen cardiac tissue are comparable in DNA methylation array ß-values, but formalin-fixed, paraffin-embedded tissue may overestimate DNA methylation levels.

Untreated fresh cardiac tissue is the optimal tissue material for investigating DNA methylation patterns of cardiac biology and diseases. However, fresh tissue is difficult to obtain. Therefore, tissue stored as frozen or formalin-fixed, paraffin-embedded (FFPE) is widely used for DNA methylation studies. It is unknown whether storage conditions alter the DNA methylation in cardiac tissue. In this study, we compared the DNA methylation patterns of fresh, frozen, and FFPE cardiac tissue to investigate if the storage method affected the DNA methylation results. We used the Infinium MethylationEPIC assay to obtain genome-wide methylation levels in fresh, frozen, and FFPE tissues from nine individuals. We found that the DNA methylation levels of 21.4% of the examined CpG sites were overestimated in the FFPE samples compared to that of fresh and frozen tissue, whereas 5.7% were underestimated. Duplicate analyses of the DNA methylation patterns showed high reproducibility (precision) for frozen and FFPE tissues. In conclusion, we found that frozen and FFPE tissues gave reproducible DNA methylation results and that frozen and fresh tissues gave similar results.

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.

Altered gene expression levels of genes related to muscle function in adults with cerebral palsy.

Cerebral palsy (CP) is the most common cause of movement disorders in children. Next generation sequencing (NGS) studies have previously shown that expression levels are fundamentally different in children with CP compared to typically developing (TD). However, given that children are in full development, we might expect gene expression levels to change once maturity is reached. Therefore, the main purpose of this study was to investigate gene expression levels of 93 target genes in adults with CP using NGS on muscle biopsies of the gastrocnemius, taken from 22 participants (n = 12 adults with CP; n = 10 TD adults). Subsequently, we carried out NGS of the mitochondrial genome to identify mtDNA variants, and additionally we studied the mitochondrial content using transmission electron microscopy images of the gastrocnemius muscle. Finally, we compared systemic ion levels between TD adults and adults with CP. Differential gene expression levels were found in genes involved in muscle contraction (MYH1 and MYBPC2), mitochondrial function kATP5J, CYCS and NDUFB6), calcium handling (CAMK2B and ATP2A), metabolism (LPL), muscle signaling (MYC, CREB1, ACVR2B, LMNA and TRIM54), and ECM (TNC). There was no statistical significant difference between CP and TD for mtDNA variant frequencies and mitochondrial content. The ion levels of Ca2+, Na+ and K+ were statistically significantly reduced while the Cl- levels were significant increased in adults with CP compared to TD adults. These results highlight that most transcriptional differences are related to muscle function in adults with CP and that mitochondrial function might be altered but not mitochondrial content.

SNP allele calling of Illumina Infinium Omni5-4 data using the butterfly method.

We introduce a within-sample SNP calling method, called the "butterfly method", that improves the quality of SNP calling with the Illumina Infinium Omni5-4 SNP Kit. This was done by improving how no-calls are determined from allele signal intensities. High confidence of SNP allele calling is extremely important in forensic genetics and clinical diagnostics. This paper is accompanied by two open-source R packages, omni54manifest and snpbeadchip that make SNP calling easy by helping with bookkeeping and giving easy access to meta-information about the SNPs typed with the Illumina Infinium Omni5-4 Kit (including chromosome, probe type, and SNP bases). We compared the results from our method with those obtained with the Illumina GenomeStudio software (which does not provide sample and SNP specific genotype probabilities or other quality measures), and with whole-genome sequencing (WGS). Given the signal intensities, the SNP calling quality was optimised using a threshold for the a posteriori probability of a SNP belonging to a SNP cluster. By lowering the a posteriori probability threshold for no-calls, we obtained a higher call rate than GenomeStudio. Using a higher a posteriori probability threshold, we achieved a higher concordance with the WGS data than GenomeStudio. Our method had SNP call and concordance rates with WGS data of approximately 99%.

Reproducibility of the Infinium methylationEPIC BeadChip assay using low DNA amounts.

The Infinium MethylationEPIC BeadChip (EPIC) is a reliable method for measuring the DNA methylation of more than 850,000 CpG positions. In clinical and forensic settings, it is critical to be able to work with low DNA amounts without risking reduced reproducibility. We evaluated the EPIC for a range of DNA amounts using two-fold serial dilutions investigated on two different days. While the ß-value distributions were generally unaffected by decreasing DNA amounts, the median squared Pearson's correlation coefficient (R2) of between-days ß-value comparisons decreased from 0.994 (500 ng DNA) to 0.957 (16 ng DNA). The median standard deviation of the ß-values was 0.005 and up to 0.017 (median of medians: 0.014) for ß-values around 0.6-0.7. With decreasing amounts of DNA from 500 ng to 16 ng, the percentage of probes with standard deviations ≤ 0.1 decreased from 99.9% to 99.4%. This study showed that high reproducibility results are obtained with DNA amounts in the range 125-500 ng DNA, while DNA amounts equal to 63 ng or below gave less reproducible results.

Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch.

Skin pigmentation is one of the most prominent and variable phenotypes in humans. We compared the alleles of 163 SNPs and indels from the Human Pigmentation (HuPi) AmpliSeq™ Custom panel, and biogeographic ancestry with the quantitative skin pigmentation levels on the upper arm, lower arm, and forehead of 299 Pakistani individuals from three subpopulations: Baloch, Pashtun, and Punjabi. The biogeographic ancestry of each individual was estimated using the Precision ID Ancestry Panel. All individuals were mainly of mixed South-Central Asian and European ancestry. However, the Baloch individuals also had an average proportion of Sub-Saharan African ancestry of approximately 10%, whereas it was <1% in the Punjabi and Pashtun individuals. The pairwise genetic distances between the Pashtun, Punjabi, and Baloch subpopulations based on the ancestry markers were statistically significantly different. Individuals from the Pashtun subpopulation had statistically significantly lower skin pigmentation than individuals from the Punjabi and Baloch subpopulations (p < 0.05). The proportions of European and Sub-Saharan African ancestry and five SNPs (rs1042602, rs10831496, rs1426654, rs16891982, and rs12913832) were statistically significantly associated with skin pigmentation at either the upper arm, lower arm or forehead in the Pakistani population after correction for multiple testing (p < 10-3). A model based on four of these SNPs (rs1426654, rs1042602, rs16891982, and rs12913832) explained 33% of the upper arm skin pigmentation. The four SNPs and the proportions of European and Sub-Saharan African ancestry explained 37% of the upper arm skin pigmentation. Our results indicate that the four likely causative SNPs, rs1426654, rs1042602, rs16891982, and rs12913832 located in SLC24A5, TYR, SLC45A2, and HERC2, respectively, are essential for skin color variation in the admixed Pakistani subpopulations.

Prediction of Eye Colour in Scandinavians Using the EyeColour 11 (EC11) SNP Set.

Description of a perpetrator's eye colour can be an important investigative lead in a forensic case with no apparent suspects. Herein, we present 11 SNPs (Eye Colour 11-EC11) that are important for eye colour prediction and eye colour prediction models for a two-category reporting system (blue and brown) and a three-category system (blue, intermediate, and brown). The EC11 SNPs were carefully selected from 44 pigmentary variants in seven genes previously found to be associated with eye colours in 757 Europeans (Danes, Swedes, and Italians). Mathematical models using three different reporting systems: a quantitative system (PIE-score), a two-category system (blue and brown), and a three-category system (blue, intermediate, brown) were used to rank the variants. SNPs with a sufficient mean variable importance (above 0.3%) were selected for EC11. Eye colour prediction models using the EC11 SNPs were developed using leave-one-out cross-validation (LOOCV) in an independent data set of 523 Norwegian individuals. Performance of the EC11 models for the two- and three-category system was compared with models based on the IrisPlex SNPs and the most important eye colour locus, rs12913832. We also compared model performances with the IrisPlex online tool (IrisPlex Web). The EC11 eye colour prediction models performed slightly better than the IrisPlex and rs12913832 models in all reporting systems and better than the IrisPlex Web in the three-category system. Three important points to consider prior to the implementation of eye colour prediction in a forensic genetic setting are discussed: (1) the reference population, (2) the SNP set, and (3) the reporting strategy.