Modified entropy-based procedure detects gene-gene-interactions in unconventional genetic models

A Novel Detection Method for High-Order SNP Epistatic Interactions Based on Explicit-Encoding-Based Multitasking Harmony Search

To elucidate the genetic basis of complex diseases, it is crucial to discover the single-nucleotide polymorphisms (SNPs) contributing to disease susceptibility. This is particularly challenging for high-order SNP epistatic interactions (HEIs), which exhibit small individual effects but potentially large joint effects. These interactions are difficult to detect due to the vast search space, encompassing billions of possible combinations, and the computational complexity of evaluating them. This study proposes a novel explicit-encoding-based multitasking harmony search algorithm (MTHS-EE-DHEI) specifically designed to address this challenge. The algorithm operates in three stages. First, a harmony search algorithm is employed, utilizing four lightweight evaluation functions, such as Bayesian network and entropy, to efficiently explore potential SNP combinations related to disease status. Second, a G-test statistical method is applied to filter out insignificant SNP combinations. Finally, two machine learning-based methods, multifactor dimensionality reduction (MDR) as well as random forest (RF), are employed to validate the classification performance of the remaining significant SNP combinations. This research aims to demonstrate the effectiveness of MTHS-EE-DHEI in identifying HEIs compared to existing methods, potentially providing valuable insights into the genetic architecture of complex diseases. The performance of MTHS-EE-DHEI was evaluated on twenty simulated disease datasets and three real-world datasets encompassing age-related macular degeneration (AMD), rheumatoid arthritis (RA), and breast cancer (BC). The results demonstrably indicate that MTHS-EE-DHEI outperforms four state-of-the-art algorithms in terms of both detection power and computational efficiency. The source code is available at https://github.com/shouhengtuo/MTHS-EE-DHEI.git .

Detecting gene-gene interactions from GWAS using diffusion kernel principal components

Genes and gene products do not function in isolation but as components of complex networks of macromolecules through physical or biochemical interactions. Dependencies of gene mutations on genetic background (i.e., epistasis) are believed to play a role in understanding molecular underpinnings of complex diseases such as inflammatory bowel disease (IBD). However, the process of identifying such interactions is complex due to for instance the curse of high dimensionality, dependencies in the data and non-linearity. Here, we propose a novel approach for robust and computationally efficient epistasis detection. We do so by first reducing dimensionality, per gene via diffusion kernel principal components (kpc). Subsequently, kpc gene summaries are used for downstream analysis including the construction of a gene-based epistasis network. We show that our approach is not only able to recover known IBD associated genes but also additional genes of interest linked to this difficult gastrointestinal disease.

A Constrained Generalized Functional Linear Model for Multi-Loci Genetic Mapping

In genome-wide association studies (GWAS), efficient incorporation of linkage disequilibria (LD) among densely typed genetic variants into association analysis is a critical yet challenging problem. Functional linear models (FLM), which impose a smoothing structure on the coefficients of correlated covariates, are advantageous in genetic mapping of multiple variants with high LD. Here we propose a novel constrained generalized FLM (cGFLM) framework to perform simultaneous association tests on a block of linked SNPs with various trait types, including continuous, binary and zero-inflated count phenotypes. The new cGFLM applies a set of inequality constraints on the FLM to ensure model identifiability under different genetic codings. The method is implemented via B-splines, and an augmented Lagrangian algorithm is employed for parameter estimation. For hypotheses testing, a test statistic that accounts for the model constraints was derived, following a mixture of chi-square distributions. Simulation results show that cGFLM is effective in identifying causal loci and gene clusters compared to several competing methods based on single markers and SKAT-C. We applied the proposed method to analyze a candidate gene-based COGEND study and a large-scale GWAS data on dental caries risk.

Roles of interacting stress-related genes in lifespan regulation: insights for translating experimental findings to humans

AIM

Experimental studies provided numerous evidence that caloric/dietary restriction may improve health and increase the lifespan of laboratory animals, and that the interplay among molecules that sense cellular stress signals and those regulating cell survival can play a crucial role in cell response to nutritional stressors. However, it is unclear whether the interplay among corresponding genes also plays a role in human health and lifespan.

METHODS

Literature about roles of cellular stressors have been reviewed, such as amino acid deprivation, and the integrated stress response (ISR) pathway in health and aging. Single nucleotide polymorphisms (SNPs) in two candidate genes (GCN2/EIF2AK4 and CHOP/DDIT3) that are closely involved in the cellular stress response to amino acid starvation, have been selected using information from experimental studies. Associations of these SNPs and their interactions with human survival in the Health and Retirement Study data have been estimated. The impact of collective associations of multiple interacting SNP pairs on survival has been evaluated, using a recently developed composite index: the SNP-specific Interaction Polygenic Risk Score (SIPRS).

RESULTS

Significant interactions have been found between SNPs from GCN2/EIF2AK4 and CHOP/DDI3T genes that were associated with survival 85+ compared to survival between ages 75 and 85 in the total sample (males and females combined) and in females only. This may reflect sex differences in genetic regulation of the human lifespan. Highly statistically significant associations of SIPRS [constructed for the rs16970024 (GCN2/EIF2AK4) and rs697221 (CHOP/DDIT3)] with survival in both sexes also been found in this study.

CONCLUSION

Identifying associations of the genetic interactions with human survival is an important step in translating the knowledge from experimental to human aging research. Significant associations of multiple SNPxSNP interactions in ISR genes with survival to the oldest old age that have been found in this study, can help uncover mechanisms of multifactorial regulation of human lifespan and its heterogeneity.