Analysis of epigenetic clocks links yoga, sleep, education, reduced meat intake, coffee, and a SOCS2 gene variant to slower epigenetic aging

DNA methylation (DNAm) clocks hold promise for measuring biological age, useful for guiding clinical interventions and forensic identification. This study compared the commonly used DNAm clocks, using DNA methylation and SNP data generated from nearly 1000 human blood or buccal swab samples. We evaluated different preprocessing methods for age estimation, investigated the association of epigenetic age acceleration (EAA) with various lifestyle and sociodemographic factors, and undertook a series of novel genome-wide association analyses for different EAA measures to find associated genetic variants. Our results highlighted the Skin&Blood clock with ssNoob normalization as the most accurate predictor of chronological age. We provided novel evidence for an association between the practice of yoga and a reduction in the pace of aging (DunedinPACE). Increased sleep and physical activity were associated with lower mortality risk score (MRS) in our dataset. University degree, vegetable consumption, and coffee intake were associated with reduced levels of epigenetic aging, whereas smoking, higher BMI, meat consumption, and manual occupation correlated well with faster epigenetic aging, with FitAge, GrimAge, and DunedinPACE clocks showing the most robust associations. In addition, we found a novel association signal for SOCS2 rs73218878 (p = 2.87 × 10-8) and accelerated GrimAge. Our study emphasizes the importance of an optimized DNAm analysis workflow for accurate estimation of epigenetic age, which may influence downstream analyses. The results support the influence of genetic background on EAA. The associated SOCS2 is a member of the suppressor of cytokine signaling family known for its role in human longevity. The reported association between various risk factors and EAA has practical implications for the development of health programs to improve quality of life and reduce premature mortality associated with age-related diseases.

Epigenetic clock in the aorta and age-related endothelial dysfunction in mice

While epigenetic age (EA) of mouse blood can be determined using DNA methylation analysis at three CpG sites in the Prima1, Hsf4 and Kcns1 genes it is not known whether this approach is useful for predicting vascular biological age. In this study we validated the 3-CpG estimator for age prediction in mouse blood, developed a new predictive model for EA in mouse aorta, and assessed whether epigenetic age acceleration (EAA) measured with blood and aorta samples correlates with age-dependent endothelial dysfunction. Endothelial function was characterized in vivo by MRI in 8-96-week-old C57BL/6 mice. Arterial stiffness was measured by USG-doppler. EA-related changes within 41 CpG sites in Prima1, Kcns1 and Hsf4 loci, were analyzed in the aorta and blood using bisulfite amplicon high-throughput sequencing. Progressive age-dependent endothelial dysfunction and changes in arterial stiffness were observed in 36-96-week-old C57BL/6 mice. Methylation levels of the investigated loci correlated with chronological age in blood and the aorta. The new model for EA estimation in aorta included three cytosines located in the Kcns1 and Hsf4, explained R2 = 87.8% of the variation in age, and predicted age with an mean absolute error of 9.6 weeks in the independent test set. EAA in the aorta was associated with endothelial dysfunction in the abdominal aorta and femoral artery what was consistent with the EAA direction estimated in blood samples. The rate of vascular biological ageing in mice, reflected by the age-dependent systemic endothelial dysfunction, could be estimated using DNA methylation measurements at three loci in aorta and blood samples.

Development of an epigenetic age predictor for costal cartilage with a simultaneous somatic tissue differentiation system

Age prediction from DNA has been a topic of interest in recent years due to the promising results obtained when using epigenetic markers. Since DNA methylation gradually changes across the individual's lifetime, prediction models have been developed accordingly for age estimation. The tissue-dependence for this biomarker usually necessitates the development of tissue-specific age prediction models, in this way, multiple models for age inference have been constructed for the most commonly encountered forensic tissues (blood, oral mucosa, semen). The analysis of skeletal remains has also been attempted and prediction models for bone have now been reported. Recently, the VISAGE Enhanced Tool was developed for the simultaneous DNA methylation analysis of 8 age-correlated loci using targeted high-throughput sequencing. It has been shown that this method is compatible with epigenetic age estimation models for blood, buccal cells, and bone. Since when dealing with decomposed cadavers or postmortem samples, cartilage samples are also an important biological source, an age prediction model for cartilage has been generated in the present study based on methylation data collected using the VISAGE Enhanced Tool. In this way, we have developed a forensic cartilage age prediction model using a training set composed of 109 samples (19-74 age range) based on DNA methylation levels from three CpGs in FHL2, TRIM59 and KLF14, using multivariate quantile regression which provides a mean absolute error (MAE) of ± 4.41 years. An independent testing set composed of 72 samples (19-75 age range) was also analyzed and provided an MAE of ± 4.26 years. In addition, we demonstrate that the 8 VISAGE markers, comprising EDARADD, TRIM59, ELOVL2, MIR29B2CHG, PDE4C, ASPA, FHL2 and KLF14, can be used as tissue prediction markers which provide reliable blood, buccal cells, bone, and cartilage differentiation using a developed multinomial logistic regression model. A training set composed of 392 samples (n = 87 blood, n = 86 buccal cells, n = 110 bone and n = 109 cartilage) was used for building the model (correct classifications: 98.72%, sensitivity: 0.988, specificity: 0.996) and validation was performed using a testing set composed of 192 samples (n = 38 blood, n = 36 buccal cells, n = 46 bone and n = 72 cartilage) showing similar predictive success to the training set (correct classifications: 97.4%, sensitivity: 0.968, specificity: 0.991). By developing both a new cartilage age model and a tissue differentiation model, our study significantly expands the use of the VISAGE Enhanced Tool while increasing the amount of DNA methylation-based information obtained from a single sample and a single forensic laboratory analysis. Both models have been placed in the open-access Snipper forensic classification website.

miR-378 influences muscle satellite cells and enhances adipogenic potential of fibro-adipogenic progenitors but does not affect muscle regeneration in the glycerol-induced injury model

Skeletal muscle regeneration relies on the reciprocal interaction between many types of cells. Regenerative capacity may be altered in different disorders. In our study, we investigated whether the deletion of miR-378a (miR-378) affects muscle regeneration. We subjected 6-week-old wild-type (WT) and miR-378 knockout (miR-378-/-) animals to the glycerol-induced muscle injury and performed analyses in various time-points. In miR-378-/- animals, an elevated abundance of muscle satellite cells (mSCs) on day 3 was found. Furthermore, fibro-adipogenic progenitors (FAPs) isolated from the muscle of miR-378-/- mice exhibited enhanced adipogenic potential. At the same time, lack of miR-378 did not affect inflammation, fibrosis, adipose tissue deposition, centrally nucleated fiber count, muscle fiber size, FAP abundance, and muscle contractility at any time point analyzed. To conclude, our study revealed that miR-378 deletion influences the abundance of mSCs and the adipogenic potential of FAPs, but does not affect overall regeneration upon acute, glycerol-induced muscle injury.

Introduction of a multiplex amplicon sequencing assay to quantify DNA methylation in target cytosine markers underlying four selected epigenetic clocks

BACKGROUND

DNA methylation analysis has proven to be a powerful tool for age assessment. However, the implementation of epigenetic age prediction in diagnostics or routine forensic casework requires appropriate laboratory methods. In this study, we aimed to compare the performance of large-scale DNA methylation analysis protocols that show promise in terms of accuracy, throughput, multiplexing capacity, and high sensitivity.

RESULTS

The protocols were designed to target a predefined panel of 161 genomic CG/CA sites from four known estimators of epigenetic age-related parameters, optimized and validated using artificially methylated controls or blood samples. We successfully targeted 96% of these loci using two enrichment protocols: Ion AmpliSeq™, an amplicon-based method integrated with Ion Torrent S5, and SureSelectXT Methyl-Seq, a hybridization-based method followed by MiSeq FGx sequencing. Both protocols demonstrated high accuracy and robustness. Although hybridization assays have greater multiplexing capabilities, the best overall performance was observed for the amplicon-based protocol with the lowest variability in DNA methylation at 25 ng of starting DNA, mean observed marker coverage of ~ 6.7 k reads, and accuracy of methylation quantification with a mean absolute difference between observed and expected methylation beta value of 0.054. The Ion AmpliSeq method correlated strongly with genome-scale EPIC microarray data (R = 0.91) and showed superiority in terms of methylation measurement accuracy. Method-to-method bias was accounted for by the use of linear transformation, which provided a highly accurate prediction of calendar age with a mean absolute error of less than 5 years for the VISAGE and Hannum age clocks used. The pace of aging (PoAm) and the mortality risk score (MRS) estimators included in our panel represent next-generation clocks, were found to have low to moderate correlations with the VISAGE and Hannum models (R < 0.75), and thus may capture different aspects of epigenetic aging.

CONCLUSIONS

We propose a laboratory tool that allows the quantification of DNA methylation in cytosines underlying four different clocks, thus providing broad information on epigenetic aging while maintaining a reasonable number of CpG markers, opening the way to a wide range of applications in forensics, medicine, and healthcare.

Development and evaluations of the ancestry informative markers of the VISAGE Enhanced Tool for Appearance and Ancestry

The VISAGE Enhanced Tool for Appearance and Ancestry (ET) has been designed to combine markers for the prediction of bio-geographical ancestry plus a range of externally visible characteristics into a single massively parallel sequencing (MPS) assay. We describe the development of the ancestry panel markers used in ET, and the enhanced analyses they provide compared to previous MPS-based forensic ancestry assays. As well as established autosomal single nucleotide polymorphisms (SNPs) that differentiate sub-Saharan African, European, East Asian, South Asian, Native American, and Oceanian populations, ET includes autosomal SNPs able to efficiently differentiate populations from Middle East regions. The ability of the ET autosomal ancestry SNPs to distinguish Middle East populations from other continentally defined population groups is such that characteristic patterns for this region can be discerned in genetic cluster analysis using STRUCTURE. Joint cluster membership estimates showing individual co-ancestry that signals North African or East African origins were detected, or cluster patterns were seen that indicate origins from central and Eastern regions of the Middle East. In addition to an augmented panel of autosomal SNPs, ET includes panels of 85 Y-SNPs, 16 X-SNPs and 21 autosomal Microhaplotypes. The Y- and X-SNPs provide a distinct method for obtaining extra detail about co-ancestry patterns identified in males with admixed backgrounds. This study used the 1000 Genomes admixed African and admixed American sample sets to fully explore these enhancements to the analysis of individual co-ancestry. Samples from urban and rural Brazil with contrasting distributions of African, European, and Native American co-ancestry were also studied to gauge the efficiency of combining Y- and X-SNP data for this purpose. The small panel of Microhaplotypes incorporated in ET were selected because they showed the highest levels of haplotype diversity amongst the seven population groups we sought to differentiate. Microhaplotype data was not formally combined with single-site SNP genotypes to analyse ancestry. However, the haplotype sequence reads obtained with ET from these loci creates an effective system for de-convoluting two-contributor mixed DNA. We made simple mixture experiments to demonstrate that when the contributors have different ancestries and the mixture ratios are imbalanced (i.e., not 1:1 mixtures) the ET Microhaplotype panel is an informative system to infer ancestry when this differs between the contributors.

DNAmFitAge: biological age indicator incorporating physical fitness

Physical fitness is a well-known correlate of health and the aging process and DNA methylation (DNAm) data can capture aging via epigenetic clocks. However, current epigenetic clocks did not yet use measures of mobility, strength, lung, or endurance fitness in their construction. We develop blood-based DNAm biomarkers for fitness parameters gait speed (walking speed), maximum handgrip strength, forced expiratory volume in one second (FEV1), and maximal oxygen uptake (VO2max) which have modest correlation with fitness parameters in five large-scale validation datasets (average r between 0.16-0.48). We then use these DNAm fitness parameter biomarkers with DNAmGrimAge, a DNAm mortality risk estimate, to construct DNAmFitAge, a new biological age indicator that incorporates physical fitness. DNAmFitAge is associated with low-intermediate physical activity levels across validation datasets (p = 6.4E-13), and younger/fitter DNAmFitAge corresponds to stronger DNAm fitness parameters in both males and females. DNAmFitAge is lower (p = 0.046) and DNAmVO2max is higher (p = 0.023) in male body builders compared to controls. Physically fit people have a younger DNAmFitAge and experience better age-related outcomes: lower mortality risk (p = 7.2E-51), coronary heart disease risk (p = 2.6E-8), and increased disease-free status (p = 1.1E-7). These new DNAm biomarkers provide researchers a new method to incorporate physical fitness into epigenetic clocks.

miR-378 affects metabolic disturbances in the mdx model of Duchenne muscular dystrophy

Although Duchenne muscular dystrophy (DMD) primarily affects muscle tissues, the alterations to systemic metabolism manifested in DMD patients contribute to the severe phenotype of this fatal disorder. We propose that microRNA-378a (miR-378) alters carbohydrate and lipid metabolism in dystrophic mdx mice. In our study, we utilized double knockout animals which lacked both dystrophin and miR-378 (mdx/miR-378^-/-). RNA sequencing of the liver identified 561 and 194 differentially expressed genes that distinguished mdx versus wild-type (WT) and mdx/miR-378^-/- versus mdx counterparts, respectively. Bioinformatics analysis predicted, among others, carbohydrate metabolism disorder in dystrophic mice, as functionally proven by impaired glucose tolerance and insulin sensitivity. The lack of miR-378 in mdx animals mitigated those effects with a faster glucose clearance in a glucose tolerance test (GTT) and normalization of liver glycogen levels. The absence of miR-378 also restored the expression of genes regulating lipid homeostasis, such as Acly, Fasn, Gpam, Pnpla3, and Scd1. In conclusion, we report for the first time that miR-378 loss results in increased systemic metabolism of mdx mice. Together with our previous finding, demonstrating alleviation of the muscle-related symptoms of DMD, we propose that the inhibition of miR-378 may represent a new strategy to attenuate the multifaceted symptoms of DMD.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

Development and inter-laboratory validation of the VISAGE enhanced tool for age estimation from semen using quantitative DNA methylation analysis

The analysis of DNA methylation has become an established method for chronological age estimation. This has triggered interest in the forensic community to develop new methods for age estimation from biological crime scene material. Various assays are available for age estimation from somatic tissues, the majority from blood. Age prediction from semen requires different DNA methylation markers and the only assays currently developed for forensic analysis are based on SNaPshot or pyrosequencing. Here, we describe a new assay using massively parallel sequencing to analyse 13 candidate CpG sites targeted in two multiplex PCRs. The assay has been validated by five consortium laboratories of the VISible Attributes through GEnomics (VISAGE) project within a collaborative exercise and was tested for reproducible quantification of DNA methylation levels and sensitivity with DNA methylation controls. Furthermore, DNA extracts and stains on Whatman FTA cards from two semen samples were used to evaluate concordance and mimic casework samples. Overall, the assay yielded high read depths (> 1000 reads) at all 13 marker positions. The methylation values obtained indicated robust quantification with an average standard deviation of 2.8% at the expected methylation level of 50% across the 13 markers and a good performance with 50 ng DNA input into bisulfite conversion. The absolute difference of quantifications from one participating laboratory to the mean quantifications of concordance and semen stains of remaining laboratories was approximately 1%. These results demonstrated the assay to be robust and suitable for age estimation from semen in forensic investigations. In addition to the 13-marker assay, a more streamlined protocol combining only five age markers in one multiplex PCR was developed. Preliminary results showed no substantial differences in DNA methylation quantification between the two assays, indicating its applicability with the VISAGE age model for semen developed with data from the complete 13-marker tool.

Impact of excessive alcohol abuse on age prediction using the VISAGE enhanced tool for epigenetic age estimation in blood

DNA methylation-based clocks provide the most accurate age estimates with practical implications for clinical and forensic genetics. However, the effects of external factors that may influence the estimates are poorly studied. Here, we evaluated the effect of alcohol consumption on epigenetic age prediction in a cohort of extreme alcohol abusers. Blood samples from deceased alcohol abusers and age- and sex-matched controls were analyzed using the VISAGE enhanced tool for age prediction from somatic tissues that enables examination of 44 CpGs within eight age markers. Significantly altered DNA methylation was recorded for alcohol abusers in MIR29B2CHG. This resulted in a mean predicted age of 1.4 years higher compared to the controls and this trend increased in older individuals. The association of alcohol abuse with epigenetic age acceleration, as determined by the prediction analysis performed based on MIR29B2CHG, was small but significant (β = 0.190; P-value = 0.007). However, the observed alteration in DNA methylation of MIR29B2CHG had a non-significant effect on age estimation with the VISAGE age prediction model. The mean absolute error in the alcohol-abusing cohort was 3.1 years, compared to 3.3 years in the control group. At the same time, upregulation of MIR29B2CHG expression may have a biological function, which merits further studies.

Searching for improvements in predicting human eye colour from DNA

Increasing understanding of human genome variability allows for better use of the predictive potential of DNA. An obvious direct application is the prediction of the physical phenotypes. Significant success has been achieved, especially in predicting pigmentation characteristics, but the inference of some phenotypes is still challenging. In search of further improvements in predicting human eye colour, we conducted whole-exome (enriched in regulome) sequencing of 150 Polish samples to discover new markers. For this, we adopted quantitative characterization of eye colour phenotypes using high-resolution photographic images of the iris in combination with DIAT software analysis. An independent set of 849 samples was used for subsequent predictive modelling. Newly identified candidates and 114 additional literature-based selected SNPs, previously associated with pigmentation, and advanced machine learning algorithms were used. Whole-exome sequencing analysis found 27 previously unreported candidate SNP markers for eye colour. The highest overall prediction accuracies were achieved with LASSO-regularized and BIC-based selected regression models. A new candidate variant, rs2253104, located in the ARFIP2 gene and identified with the HyperLasso method, revealed predictive potential and was included in the best-performing regression models. Advanced machine learning approaches showed a significant increase in sensitivity of intermediate eye colour prediction (up to 39%) compared to 0% obtained for the original IrisPlex model. We identified a new potential predictor of eye colour and evaluated several widely used advanced machine learning algorithms in predictive analysis of this trait. Our results provide useful hints for developing future predictive models for eye colour in forensic and anthropological studies.

DNA methylation-based age clocks: From age prediction to age reversion

Aging as an irretrievable occurrence throughout the entire life is characterized by a progressive decline in physiological functionality and enhanced disease vulnerability. Numerous studies have demonstrated that epigenetic modifications, particularly DNA methylation (DNAm), correlate with aging and age-related diseases. Several investigations have attempted to predict chronological age using the age-related alterations in the DNAm of certain CpG sites. Here we categorize different studies that tracked the aging process in the DNAm landscape to show how epigenetic age clocks evolved from a chronological age estimator to an indicator of lifespan and healthspan. We also describe the health and disease predictive potential of estimated epigenetic age acceleration regarding different clinical conditions and lifestyle factors. Considering the revealed age-related epigenetic changes, the recent age-reprogramming strategies are discussed which are promising methods for resetting the aging clocks.

Development of the VISAGE enhanced tool and statistical models for epigenetic age estimation in blood, buccal cells and bones

DNA methylation is known as a biomarker for age with applications in forensics. Here we describe the VISAGE (VISible Attributes through GEnomics) Consortium's enhanced tool for epigenetic age estimation in somatic tissues. The tool is based on eight DNA methylation markers (44 CpGs), bisulfite multiplex PCR followed by sequencing on the MiSeq FGx platform, and three statistical prediction models for blood, buccal cells and bones. The model for blood is based on six CpGs from ELOVL2, MIR29B2CHG, KLF14, FHL2, TRIM59 and PDE4C, and predicts age with a mean absolute error (MAE) of 3.2 years, while the model for buccal cells includes five CpGs from PDE4C, MIR29B2CHG, ELOVL2, KLF14 and EDARADD and predicts age with MAE of 3.7 years, and the model for bones has six CpGs from ELOVL2, KLF14, PDE4C and ASPA and predicts age with MAE of 3.4 years. The VISAGE enhanced tool for age estimation in somatic tissues enables reliable collection of DNA methylation data from small amounts of DNA using a sensitive multiplex MPS assay that provides accurate estimation of age in blood, buccal swabs, and bones using the statistical model tailored to each tissue.

Deletion of Mcpip1 in Mcpip1fl/flAlbCre mice recapitulates the phenotype of human primary biliary cholangitis

Primary biliary cholangitis (PBC) is an autoimmune disease characterized by progressive destruction of the intrahepatic bile ducts. The immunopathology of PBC involves excessive inflammation; therefore, negative regulators of inflammatory response, such as Monocyte Chemoattractant Protein-1-Induced Protein-1 (MCPIP1) may play important roles in the development of PBC. The aim of this work was to verify whether Mcpip1 expression protects against development of PBC. Genetic deletion of Zc3h12a was used to characterize the role of Mcpip1 in the pathogenesis of PBC in 6–52-week-old mice. We found that Mcpip1 deficiency in the liver (Mcpip1^fl/flAlb^Cre) recapitulates most of the features of human PBC, in contrast to mice with Mcpip1 deficiency in myeloid cells (Mcpip1^fl/flLysM^Cre mice), which present with robust myeloid cell-driven systemic inflammation. In Mcpip1^fl/flAlb^Cre livers, intrahepatic bile ducts displayed proliferative changes with inflammatory infiltration, bile duct destruction, and fibrosis leading to cholestasis. In plasma, increased concentrations of IgG, IgM, and AMA autoantibodies (anti-PDC-E2) were detected. Interestingly, the phenotype of Mcpip1^fl/flAlb^Cre mice was robust in 6-week-old, but milder in 12–24-week-old mice. Hepatic transcriptome analysis of 6-week-old and 24-week-old Mcpip1^fl/flAlb^Cre mice showed 812 and 8 differentially expressed genes, respectively, compared with age-matched control mice, and revealed a distinct set of genes compared to those previously associated with development of PBC. In conclusion, Mcpip1^fl/flAlb^Cre mice display early postnatal phenotype that recapitulates most of the features of human PBC.

Testing the impact of trait prevalence priors in Bayesian-based genetic prediction modeling of human appearance traits

The prediction of appearance traits by use of solely genetic information has become an established approach and a number of statistical prediction models have already been developed for this purpose. However, given limited knowledge on appearance genetics, currently available models are incomplete and do not include all causal genetic variants as predictors. Therefore such prediction models may benefit from the inclusion of additional information that acts as a proxy for this unknown genetic background. Use of priors, possibly informed by trait category prevalence values in biogeographic ancestry groups, in a Bayesian framework may thus improve the prediction accuracy of previously predicted externally visible characteristics, but has not been investigated as of yet. In this study, we assessed the impact of using trait prevalence-informed priors on the prediction performance in Bayesian models for eye, hair and skin color as well as hair structure and freckles in comparison to the respective prior-free models. Those prior-free models were either similarly defined either very close to the already established ones by using a reduced predictive marker set. However, these differences in the number of the predictive markers should not affect significantly our main outcomes. We observed that such priors often had a strong effect on the prediction performance, but to varying degrees between different traits and also different trait categories, with some categories barely showing an effect. While we found potential for improving the prediction accuracy of many of the appearance trait categories tested by using priors, our analyses also showed that misspecification of those prior values often severely diminished the accuracy compared to the respective prior-free approach. This emphasizes the importance of accurate specification of prevalence-informed priors in Bayesian prediction modeling of appearance traits. However, the existing literature knowledge on spatial prevalence is sparse for most appearance traits, including those investigated here. Due to the limitations in appearance trait prevalence knowledge, our results render the use of trait prevalence-informed priors in DNA-based appearance trait prediction currently infeasible.

A Comparison of Forensic Age Prediction Models Using Data From Four DNA Methylation Technologies

Individual age estimation can be applied to criminal, legal, and anthropological investigations. DNA methylation has been established as the biomarker of choice for age prediction, since it was observed that specific CpG positions in the genome show systematic changes during an individual’s lifetime, with progressive increases or decreases in methylation levels. Subsequently, several forensic age prediction models have been reported, providing average age prediction error ranges of ±3–4 years, using a broad spectrum of technologies and underlying statistical analyses. DNA methylation assessment is not categorical but quantitative. Therefore, the detection platform used plays a pivotal role, since quantitative and semi-quantitative technologies could potentially result in differences in detected DNA methylation levels. In the present study, we analyzed as a shared sample pool, 84 blood-based DNA controls ranging from 18 to 99 years old using four different technologies: EpiTYPER^®, pyrosequencing, MiSeq, and SNaPshot^TM. The DNA methylation levels detected for CpG sites from ELOVL2, FHL2, and MIR29B2 with each system were compared. A restricted three CpG-site age prediction model was rebuilt for each system, as well as for a combination of technologies, based on previous training datasets, and age predictions were calculated accordingly for all the samples detected with the previous technologies. While the DNA methylation patterns and subsequent age predictions from EpiTYPER^®, pyrosequencing, and MiSeq systems are largely comparable for the CpG sites studied, SNaPshot^TM gives bigger differences reflected in higher predictive errors. However, these differences can be reduced by applying a z-score data transformation.

Exploring the possibility of predicting human head hair greying from DNA using whole-exome and targeted NGS data

Background

Greying of the hair is an obvious sign of human aging. In addition to age, sex- and ancestry-specific patterns of hair greying are also observed and the progression of greying may be affected by environmental factors. However, little is known about the genetic control of this process. This study aimed to assess the potential of genetic data to predict hair greying in a population of nearly 1000 individuals from Poland.

Results

The study involved whole-exome sequencing followed by targeted analysis of 378 exome-wide and literature-based selected SNPs. For the selection of predictors, the minimum redundancy maximum relevance (mRMRe) method was used, and then two prediction models were developed. The models included age, sex and 13 unique SNPs. Two SNPs of the highest mRMRe score included whole-exome identified KIF1A rs59733750 and previously linked with hair loss FGF5 rs7680591. The model for greying vs. no greying prediction achieved accuracy of cross-validated AUC = 0.873. In the 3-grade classification cross-validated AUC equalled 0.864 for no greying, 0.791 for mild greying and 0.875 for severe greying. Although these values present fairly accurate prediction, most of the prediction information was brought by age alone. Genetic variants explained < 10% of hair greying variation and the impact of particular SNPs on prediction accuracy was found to be small.

Conclusions

The rate of changes in human progressive traits shows inter-individual variation, therefore they are perceived as biomarkers of the biological age of the organism. The knowledge on the mechanisms underlying phenotypic aging can be of special interest to the medicine, cosmetics industry and forensics. Our study improves the knowledge on the genetics underlying hair greying processes, presents prototype models for prediction and proves hair greying being genetically a very complex trait. Finally, we propose a four-step approach based on genetic and epigenetic data analysis allowing for i) sex determination; ii) genetic ancestry inference; iii) greying-associated SNPs assignment and iv) epigenetic age estimation, all needed for a final prediction of greying.

Angiotensin converting enzyme: A review on expression profile and its association with human disorders with special focus on SARS-CoV-2 infection

Angiotensin-converting enzyme (ACE) and its homologue, ACE2, have been mostly associated with hypertensive disorder. However, recent pandemia of SARS-CoV-2 has put these proteins at the center of attention, as this virus has been shown to exploit ACE2 protein to enter cells. Clear difference in the response of affected patients to this virus has urged researchers to find the molecular basis and pathophysiology of the cell response to this virus. Different levels of expression and function of ACE proteins, underlying disorders, consumption of certain medications and the existence of certain genomic variants within ACE genes are possible explanations for the observed difference in the response of individuals to the SARS-CoV-2 infection. In the current review, we discuss the putative mechanisms for this observation.

The challenge of predicting human pigmentation traits in degraded bone samples with the MPS-based HIrisPlex-S system

Identification of human remains is an important part of human DNA analysis studies. STR and mitochondrial DNA markers are well suited for the analysis of degraded biological samples including bone material. However, these DNA markers may be useless when reference material is not available. In these cases, predictive DNA analysis can support the process of human identification by providing investigative leads. Forensic DNA phenotyping has progressed significantly by offering new methods based on massively parallel sequencing technology, but the frequent degradation processes observed in skeletal remains can make analysis of such samples challenging. In this study, we demonstrate the usefulness of a recently established Ion AmpliSeq^TM HIrisPlex-S panel using Ion Torrent technology for analyzing bone samples that show different levels of DNA degradation. In total, 63 bone samples at post-mortem intervals up to almost 80 years were genotyped and eye, hair and skin colour predictions were performed using the HIrisPlex-S models. Following the recommended coverage thresholds, it was possible to establish full DNA profiles comprising of 41 DNA variants for 35 samples (55.6%). For 5 samples (7.9%) no DNA profiles were generated. The remaining 23 samples (36.5%) produced partial profiles and showed a clear underperformance of 3 HIrisPlex-S SNPs - rs1545397 (OCA2), rs1470608 (OCA2) and rs10756819 (BNC2), all used for skin colour prediction only. None of the 23 samples gave complete genotypes needed for skin colour prediction was obtained, and in 7 of them (25.9%) the 3 underperformed SNPs were the cause. At the same time, the prediction of eye and hair colour using complete IrisPlex and HIrisPlex profiles could be made for these 23 samples in 20 (87.0%) and 12 cases (52.2%), respectively. Complete HIrisPlex-S profiles were generated from as little as 49 pg of template DNA. Five samples for which the HIrisPlex-S analysis failed, consistently failed in standard STR analysis. Importantly, the 3 underperforming SNPs produced significantly lower number of reads in good quality samples. Nonetheless, the AUC loss resulting from missing data for these 3 SNPs is not considered large (≤0.004) and the prediction of pigmentation from partial profiles is also available in the current HPS tool. The study shows that DNA degradation and the resulting loss of data are the most serious challenge to DNA phenotyping of skeletal remains. Although the newly developed HIrisPlex-S panel has been successfully validated in the current research, primer redesign for the 3 underperforming SNPs in the MPS design should be considered in the future.

Modified aging of elite athletes revealed by analysis of epigenetic age markers

Recent progress in epigenomics has led to the development of prediction systems that enable accurate age estimation from DNA methylation data. Our objective was to track responses to intense physical exercise of individual age-correlated DNA methylation markers and to infer their potential impact on the aging processes. The study showed accelerated DNA hypermethylation for two CpG sites in TRIM59 and KLF14. Both markers predicted the investigated elite athletes to be several years older than controls and this effect was more substantial in subjects involved in power sports. Accordingly, the complete 5-CpG model revealed age acceleration of elite athletes (P=1.503x10^-7) and the result was more significant amongst power athletes (P=1.051x10^-9). The modified methylation of TRIM59 and KLF14 in top athletes may be accounted for by the biological roles played by these genes. Their known anti-tumour and anti-inflammatory activities suggests that intense physical training has a complex influence on aging and potentially launches signalling networks that contribute to the observed lower risk of elite athletes to develop cardiovascular disease and cancer.

Non-CYP2D6 Variants Selected by a GWAS Improve the Prediction of Impaired Tamoxifen Metabolism in Patients with Breast Cancer

A certain minimum plasma concentration of (Z)-endoxifen is presumably required for breast cancer patients to benefit from tamoxifen therapy. In this study, we searched for DNA variants that could aid in the prediction of risk for insufficient (Z)-endoxifen exposure. A metabolic ratio (MR) corresponding to the (Z)-endoxifen efficacy threshold level was adopted as a cutoff value for a genome-wide association study comprised of 287 breast cancer patients. Multivariate regression was used to preselect variables exhibiting an independent impact on the MR and develop models to predict below-threshold MR values. In total, 15 single-nucleotide polymorphisms (SNPs) were significantly associated with below-threshold MR values. The strongest association was with rs8138080 (WBP2NL). Two alternative models for MR prediction were developed. The predictive accuracy of Model 1, including rs7245, rs6950784, rs1320308, and the CYP2D6 genotype, was considerably higher than that of the CYP2D6 genotype alone (AUC 0.879 vs 0.758). Model 2, which was developed using the same three SNPs as for Model 1 plus rs8138080, appeared as an interesting alternative to the full CYP2D6 genotype testing. In conclusion, the four novel SNPs, tested alone or in combination with the CYP2D6 genotype, improved the prediction of impaired tamoxifen-to-endoxifen metabolism, potentially allowing for treatment optimization.

Towards broadening Forensic DNA Phenotyping beyond pigmentation: Improving the prediction of head hair shape from DNA

Human head hair shape, commonly classified as straight, wavy, curly or frizzy, is an attractive target for Forensic DNA Phenotyping and other applications of human appearance prediction from DNA such as in paleogenetics. The genetic knowledge underlying head hair shape variation was recently improved by the outcome of a series of genome-wide association and replication studies in a total of 26,964 subjects, highlighting 12 loci of which 8 were novel and introducing a prediction model for Europeans based on 14 SNPs. In the present study, we evaluated the capacity of DNA-based head hair shape prediction by investigating an extended set of candidate SNP predictors and by using an independent set of samples for model validation. Prediction model building was carried out in 9674 subjects (6068 from Europe, 2899 from Asia and 707 of admixed European and Asian ancestries), used previously, by considering a novel list of 90 candidate SNPs. For model validation, genotype and phenotype data were newly collected in 2415 independent subjects (2138 Europeans and 277 non-Europeans) by applying two targeted massively parallel sequencing platforms, Ion Torrent PGM and MiSeq, or the MassARRAY platform. A binomial model was developed to predict straight vs. non-straight hair based on 32 SNPs from 26 genetic loci we identified as significantly contributing to the model. This model achieved prediction accuracies, expressed as AUC, of 0.664 in Europeans and 0.789 in non-Europeans; the statistically significant difference was explained mostly by the effect of one EDAR SNP in non-Europeans. Considering sex and age, in addition to the SNPs, slightly and insignificantly increased the prediction accuracies (AUC of 0.680 and 0.800, respectively). Based on the sample size and candidate DNA markers investigated, this study provides the most robust, validated, and accurate statistical prediction models and SNP predictor marker sets currently available for predicting head hair shape from DNA, providing the next step towards broadening Forensic DNA Phenotyping beyond pigmentation traits.

Hot on the Trail of Genes that Shape Our Fingerprints

Fingerprint patterns have been associated with their ability to identify an individual uniquely, but can uniqueness be understood genetically? Ho et al. point out some key variants that may be responsible for some of the concentric patterns that are observed in digital skin. Furthermore, they propose that one highly associated gene, ADAMTS9-AS2, has a role in epigenetic regulation, a role that may be important in early-stage digit development.

Development of a methylation marker set for forensic age estimation using analysis of public methylation data and the Agena Bioscience EpiTYPER system

Individual age estimation has the potential to provide key information that could enhance and extend DNA intelligence tools. Following predictive tests for externally visible characteristics developed in recent years, prediction of age could guide police investigations and improve the assessment of age-related phenotype expression patterns such as hair colour changes and early onset of male pattern baldness. DNA methylation at CpG positions has emerged as the most promising DNA tests to ascertain the individual age of the donor of a biological contact trace. Although different methodologies are available to detect DNA methylation, EpiTYPER technology (Agena Bioscience, formerly Sequenom) provides useful characteristics that can be applied as a discovery tool in localized regions of the genome. In our study, a total of twenty-two candidate genomic regions, selected from the assessment of publically available data from the Illumina HumanMethylation 450 BeadChip, had a total of 177 CpG sites with informative methylation patterns that were subsequently investigated in detail. From the methylation analyses made, a novel age prediction model based on a multivariate quantile regression analysis was built using the seven highest age-correlated loci of ELOVL2, ASPA, PDE4C, FHL2, CCDC102B, C1orf132 and chr16:85395429. The detected methylation levels in these loci provide a median absolute age prediction error of ±3.07years and a percentage of prediction error relative to the age of 6.3%. We report the predictive performance of the developed model using cross validation of a carefully age-graded training set of 725 European individuals and a test set of 52 monozygotic twin pairs. The multivariate quantile regression age predictor, using the CpG sites selected in this study, has been placed in the open-access Snipper forensic classification website.

Further evidence for population specific differences in the effect of DNA markers and gender on eye colour prediction in forensics

The genetics of eye colour has been extensively studied over the past few years, and the identified polymorphisms have been applied with marked success in the field of Forensic DNA Phenotyping. A picture that arises from evaluation of the currently available eye colour prediction markers shows that only the analysis of HERC2-OCA2 complex has similar effectiveness in different populations, while the predictive potential of other loci may vary significantly. Moreover, the role of gender in the explanation of human eye colour variation should not be neglected in some populations. In the present study, we re-investigated the data for 1020 Polish individuals and using neural networks and logistic regression methods explored predictive capacity of IrisPlex SNPs and gender in this population sample. In general, neural networks provided higher prediction accuracy comparing to logistic regression (AUC increase by 0.02–0.06). Four out of six IrisPlex SNPs were associated with eye colour in the studied population. HERC2 rs12913832, OCA2 rs1800407 and SLC24A4 rs12896399 were found to be the most important eye colour predictors (p < 0.007) while the effect of rs16891982 in SLC45A2 was less significant. Gender was found to be significantly associated with eye colour with males having ~1.5 higher odds for blue eye colour comparing to females (p = 0.002) and was ranked as the third most important factor in blue/non-blue eye colour determination. However, the implementation of gender into the developed prediction models had marginal and ambiguous impact on the overall accuracy of prediction confirming that the effect of gender on eye colour in this population is small. Our study indicated the advantage of neural networks in prediction modeling in forensics and provided additional evidence for population specific differences in the predictive importance of the IrisPlex SNPs and gender.

Evaluation of the predictive capacity of DNA variants associated with straight hair in Europeans

DNA-based prediction of hair morphology, defined as straight, curly or wavy hair, could contribute to an improved description of an unknown offender and allow more accurate forensic reconstructions of physical appearance in the field of forensic DNA phenotyping. Differences in scalp hair morphology are significant at the worldwide scale and within Europe. The only genome-wide association study made to date revealed the Trichohyalin gene (TCHH) to be significantly associated with hair morphology in Europeans and reported weaker associations for WNT10A and FRAS1 genes. We conducted a study that centered on six SNPs located in these three genes with a sample of 528 individuals from Poland. The predictive capacity of the candidate DNA variants was evaluated using logistic regression; classification and regression trees; and neural networks, by applying a 10-fold cross validation procedure. Additionally, an independent test set of 142 males from six European populations was used to verify performance of the developed prediction models. Our study confirmed association of rs11803731 (TCHH), rs7349332 (WNT10A) and rs1268789 (FRAS1) SNPs with hair morphology. The combined genotype risk score for straight hair had an odds ratio of 2.7 and these predictors explained ∼ 8.2% of the total variance. The selected three SNPs were found to predict straight hair with a high sensitivity but low specificity when a 10-fold cross validation procedure was applied and the best results were obtained using the neural networks approach (AUC=0.688, sensitivity=91.2%, specificity=23.0%). Application of the neural networks model with 65% probability threshold on an additional test set gave high sensitivity (81.4%) and improved specificity (50.0%) with a total of 78.7% correct calls, but a high non-classification rate (66.9%). The combined TTGGGG SNP genotype for rs11803731, rs7349332, rs1268789 (European frequency=4.5%) of all six straight hair-associated alleles was identified as the best predictor, giving >80% probability of straight hair. Finally, association testing of 44 SNPs previously identified to be associated with male pattern baldness revealed a suggestive association with hair morphology for rs4679955 on 3q25.1. The study results reported provide the starting point for the development of a predictive test for hair morphology in Europeans. More studies are now needed to discover additional determinants of hair morphology to improve the predictive accuracy of this trait in forensic analysis.

Evaluation of DNA variants associated with androgenetic alopecia and their potential to predict male pattern baldness

Androgenetic alopecia, known in men as male pattern baldness (MPB), is a very conspicuous condition that is particularly frequent among European men and thus contributes markedly to variation in physical appearance traits amongst Europeans. Recent studies have revealed multiple genes and polymorphisms to be associated with susceptibility to MPB. In this study, 50 candidate SNPs for androgenetic alopecia were analyzed in order to verify their potential to predict MPB. Significant associations were confirmed for 29 SNPs from chromosomes X, 1, 5, 7, 18 and 20. A simple 5-SNP prediction model and an extended 20-SNP model were developed based on a discovery panel of 305 males from various European populations fitting one of two distinct phenotype categories. The first category consisted of men below 50 years of age with significant baldness and the second; men aged 50 years or older lacking baldness. The simple model comprised the five best predictors: rs5919324 near AR, rs1998076 in the 20p11 region, rs929626 in EBF1, rs12565727 in TARDBP and rs756853 in HDAC9. The extended prediction model added 15 SNPs from five genomic regions that improved overall prevalence-adjusted predictive accuracy measured by area under the receiver characteristic operating curve (AUC). Both models were evaluated for predictive accuracy using a test set of 300 males reflecting the general European population. Applying a 65% probability threshold, high prediction sensitivity of 87.1% but low specificity of 42.4% was obtained in men aged <50 years. In men aged ≥50, prediction sensitivity was slightly lower at 67.7% while specificity reached 90%. Overall, the AUC=0.761 calculated for men at or above 50 years of age indicates these SNPs offer considerable potential for the application of genetic tests to predict MPB patterns, adding a highly informative predictive system to the emerging field of forensic analysis of externally visible characteristics.

Variants of SCARB1 and VDR Involved in Complex Genetic Interactions May Be Implicated in the Genetic Susceptibility to Clear Cell Renal Cell Carcinoma

The current data are still inconclusive in terms of a genetic component involved in the susceptibility to renal cell carcinoma. Our aim was to evaluate 40 selected candidate polymorphisms for potential association with clear cell renal cell carcinoma (ccRCC) based on independent group of 167 patients and 200 healthy controls. The obtained data were searched for independent effects of particular polymorphisms as well as haplotypes and genetic interactions. Association testing implied position rs4765623 in the SCARB1 gene (OR = 1.688, 95% CI: 1.104-2.582, P = 0.016) and a haplotype in VDR comprising positions rs739837, rs731236, rs7975232, and rs1544410 (P = 0.012) to be the risk factors in the studied population. The study detected several epistatic effects contributing to the genetic susceptibility to ccRCC. Variation in GNAS1 was implicated in a strong synergistic interaction with BIRC5. This effect was part of a model suggested by multifactor dimensionality reduction method including also a synergy between GNAS1 and SCARB1 (P = 0.036). Significance of GNAS1-SCARB1 interaction was further confirmed by logistic regression (P = 0.041), which also indicated involvement of SCARB1 in additional interaction with EPAS1 (P = 0.008) as well as revealing interactions between GNAS1 and EPAS1 (P = 0.016), GNAS1 and MC1R (P = 0.031), GNAS1 and VDR (P = 0.032), and MC1R and VDR (P = 0.035).

Increased risk of developing cutaneous malignant melanoma is associated with variation in pigmentation genes and VDR, and may involve epistatic effects

Cutaneous malignant melanoma (CMM) is a malicious human skin cancer that primarily affects individuals with light pigmentation and heavy sun exposure, but also has a known familial association. Multiple genes and polymorphisms have been reported as low-penetrance susceptibility loci for CMM. Here, we examined 33 candidate polymorphisms located in 11 pigmentation genes and the vitamin D receptor gene (VDR) in a population of 130 cutaneous melanoma patients and 707 healthy controls. The genotypes obtained were evaluated for main association effects and potential gene-gene interactions. MC1R, TYR, VDR and SLC45A2 genes were found to be associated with CMM in our population. The results obtained for major function MC1R mutations were the most significant [with odds ratio (OR)=1.787, confidence interval (CI)=1.320-2.419 and P=1.715(-4)], followed by TYR (rs1393350) (with OR=1.569, CI=1.162-2.118, P=0.003), VDR (GCCC haplotype in rs2238136-rs4516035-rs7139166-rs11568820 block) (with OR=5.653, CI=1.794-17.811, P=0.003) and SLC45A2 (rs16891982) (with OR=0.238, CI=0.057-0.987, P=0.048). The study also detected significant intermolecular epistatic effects between MC1R and TYR, SLC45A2 and VDR, HERC2 and VDR, OCA2 and TPCN2, as well as intramolecular interactions between variants within the genes MC1R and VDR. In the final multivariate logistic regression model for CMM development, only the gene-gene interactions discovered remained significant, showing that epistasis may be an important factor in the risk of melanoma.

Prediction of eye color in the Slovenian population using the IrisPlex SNPs

AIM

To evaluate the accuracy of eye color prediction based on six IrisPlex single nucleotide polymorphisms (SNP) in a Slovenian population sample.

METHODS

Six IrisPlex predictor SNPs (HERC2 - rs12913832, OCA2 - rs1800407, SLC45A2 - rs16891982 and TYR - rs1393350, SLC24A4 - rs12896399, and IRF4 - rs12203592) of 105 individuals were analyzed using single base extension approach and SNaPshot chemistry. The IrisPlex multinomial regression prediction model was used to infer eye color probabilities. The accuracy of the IrisPlex was assessed through the calculation of sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and the area under the receiver characteristic operating curves (AUC).

RESULTS

Blue eye color was observed in 44.7%, brown in 29.6%, and intermediate in 25.7% participants. Prediction accuracy expressed by the AUC was 0.966 for blue, 0.913 for brown, and 0.796 for intermediate eye color. Sensitivity was 93.6% for blue, 58.1% for brown, and 0% for intermediate eye color. Specificity was 93.1% for blue, 89.2% for brown, and 100% for intermediate eye color. PPV was 91.7% for blue and 69.2% for brown color. NPV was 94.7% for blue and 83.5% for brown eye color. These values indicate prediction accuracy comparable to that established in other studies.

CONCLUSION

Blue and brown eye color can be reliably predicted from DNA samples using only six polymorphisms, while intermediate eye color defies prediction, indicating that more research is needed to genetically predict the whole variation of eye color in humans.

Bona fide colour: DNA prediction of human eye and hair colour from ancient and contemporary skeletal remains

BACKGROUND

DNA analysis of ancient skeletal remains is invaluable in evolutionary biology for exploring the history of species, including humans. Contemporary human bones and teeth, however, are relevant in forensic DNA analyses that deal with the identification of perpetrators, missing persons, disaster victims or family relationships. They may also provide useful information towards unravelling controversies that surround famous historical individuals. Retrieving information about a deceased person's externally visible characteristics can be informative in both types of DNA analyses. Recently, we demonstrated that human eye and hair colour can be reliably predicted from DNA using the HIrisPlex system. Here we test the feasibility of the novel HIrisPlex system at establishing eye and hair colour of deceased individuals from skeletal remains of various post-mortem time ranges and storage conditions.

METHODS

Twenty-one teeth between 1 and approximately 800 years of age and 5 contemporary bones were subjected to DNA extraction using standard organic protocol followed by analysis using the HIrisPlex system.

RESULTS

Twenty-three out of 26 bone DNA extracts yielded the full 24 SNP HIrisPlex profile, therefore successfully allowing model-based eye and hair colour prediction. HIrisPlex analysis of a tooth from the Polish general Władysław Sikorski (1881 to 1943) revealed blue eye colour and blond hair colour, which was positively verified from reliable documentation. The partial profiles collected in the remaining three cases (two contemporary samples and a 14th century sample) were sufficient for eye colour prediction.

CONCLUSIONS

Overall, we demonstrate that the HIrisPlex system is suitable, sufficiently sensitive and robust to successfully predict eye and hair colour from ancient and contemporary skeletal remains. Our findings, therefore, highlight the HIrisPlex system as a promising tool in future routine forensic casework involving skeletal remains, including ancient DNA studies, for the prediction of eye and hair colour of deceased individuals.

The impact of mitochondrial and nuclear DNA variants on late-onset Alzheimer's disease risk

We investigated the potential contribution of mitochondrial DNA (mtDNA) variants, haplogroups, and polymorphisms in nuclear genes essential for mitochondrial biogenesis and function (PGC-1α TFAM) to late-onset Alzheimer's disease (LOAD) risk. Epistatic interaction analysis was conducted between the studied variables. Our results demonstrate that mtDNA haplogroups and subhaplogroups with putative role in partial uncoupling of oxidative phosphorylation are significantly associated with a decreased LOAD risk (OR <1). Conversely, mtDNA haplogroup H (p = 0.049) and HV cluster (p = 0.018) are significant LOAD risk factors, which was additionally confirmed by meta-analysis (OR = 1.22, OR = 1.25, respectively). Haplogroup K was demonstrated to exert a neutralizing effect on the high risk associated with APOE4+ status (p = 0.014). Further, two synergistic interactions between subhaplogroup H5 and APOE4 status (p = 0.009) and between TFAM rs1937 and APOE4 status (p < 0.001) were detected, influencing LOAD risk. No interaction pointing to a dual mitochondrial-nuclear genome variation effect on LOAD occurrence was identified.

Gene-gene interactions contribute to eye colour variation in humans

Prediction of phenotypes from genetic data is considered to be the first practical application of data gained from association studies, with potential importance for medicine and the forensic sciences. Multiple genes and polymorphisms have been found to be associated with variation in human pigmentation. Their analysis enables prediction of blue and brown eye colour with a reasonably high accuracy. More accurate prediction, especially in the case of intermediate eye colours, may require better understanding of gene-gene interactions affecting this polygenic trait. Using multifactor dimensionality reduction and logistic regression methods, a study of gene-gene interactions was conducted based on variation in 11 known pigmentation genes examined in a cohort of 718 individuals of European descent. The study revealed significant interactions of a redundant character between the HERC2 and OCA2 genes affecting determination of hazel eye colour and between HERC2 and SLC24A4 affecting determination of blue eye colour. Our research indicates interactive effects of a synergistic character between HERC2 and OCA2, and also provides evidence for a novel strong synergistic interaction between HERC2 and TYRP1, both affecting determination of green eye colour.