A generalized combinatorial approach for detecting gene-by-gene and gene-by-environment interactions with application to nicotine dependence

Involvement of Single Nucleotide Variants in the Klotho Gene Among Obesity Individuals with and without Type 2 Diabetes Mellitus in the Saudi Population

Purpose

Aging is characterized by the gradual physiological changes and alterations that accumulate over time in the human body. The combination of obesity and ageing can lead to an increased risk of serious health issues or death. Single nucleotide variants (SNVs) in the Klotho gene were commonly studied, including that in type 2 diabetes mellitus (T2DM).

Aim

The aim of this study is to examine the possible effect of SNVs in Klotho on the obese population in Saudi Arabia using middle-aged participants with and without T2DM.

Methods

This study consists of 100 controls and 100 obesity patients, in which 50 had T2DM and the remaining 50 were obese without T2DM. Genotyping was performed with PCR, and Sanger sequencing analysis was used to validate the molecular association.

Results

In this study, rs1207568 (p = 0.001-0.003) and rs9527025 (p = 0.001-0.00004) SNVs were associated with obesity cases. However, none of the genotypes or allele frequencies showed a positive association with the rs564481 SNV (p = 0.344-0.881). The multiple linear regression model showed that waist and hip were associated (p = 0.01-0.02). ANOVA analysis showed age (p = 0.04), hip (p = 0.002), SBP, and TC (p = 0.02) were associated. Finally, SNV (rs1207568 and rs95270250) and obesity (p < 0.001) associations were confirmed through gene multifactor dimensionality reduction analysis with gene-gene interaction, dendrogram, and graphical depletion method.

Conclusion

This study concludes that rs1207568 and rs9527025 SNVs are associated with obesity in the Saudi population. Additional genetical statistics showed significant association between dependent and independent variables. SNVs in Klotho play a role in the Saudi population's susceptibility to obesity.

Gene‒Environment Interaction Between CAST Gene and Eye-Rubbing in the Chinese Keratoconus Cohort Study: A Case-Only Study

Purpose

Keratoconus (KC), characterized by progressive corneal protrusion and thinning, is a complex disease influenced by the combination of genetic and environmental factors. The purpose of this study was to explore potential gene‒environment interaction between the calpastatin (CAST) gene and eye-rubbing in KC.

Methods

A case-only study including 930 patients (676 patients with eye-rubbing and 254 patients without eye-rubbing) from the Chinese Keratoconus (CKC) cohort study was performed in the present study. Genotyping of single nucleotide polymorphism (SNP) was conducted using the Illumina Infinium Human Asian Screening Array (ASA) Beadchip. The gene‒environment interactions between CAST gene and eye-rubbing were analyzed using PLINK version 1.90. The interactions between CAST genotypes and eye-rubbing were analyzed by logistic regression models. The SNP-SNP-environment interactions were analyzed using generalized multifactor dimensionality reduction (GMDR).

Results

Three SNPs in CAST gene, namely, rs26515, rs27991, and rs9314177, reached the significance threshold for interactions (defined as P < 2.272 × 10-3). Notably, the minor alleles of these three SNPs exhibited negative interactions with eye-rubbing in KC. The results of logistic regression models revealed that the minor allele homozygotes and heterozygotes of rs26515, rs27991, and rs9314177 also exhibited negative interactions with eye-rubbing. Furthermore, GMDR analysis revealed the significant SNP-SNP-environment interactions among rs26515, rs27991, rs9314177, and eye-rubbing in KC.

Conclusions

This study identified rs26515, rs27991, and rs9314177 in CAST gene existed gene-environment interactions with eye-rubbing in KC, which is highly important for understanding the underlying biological mechanisms of KC and guiding precision prevention and proper management.

Interaction between HTR2A rs3125 and negative life events in suicide attempts among patients with major depressive disorder: a cross-sectional study

BACKGROUND

Both genetic and environmental factors play crucial roles in the development of major depressive disorder (MDD) and suicide attempts (SA). However, the interaction between both items remains unknown. This study aims to explore the interactions between the genetic variants of the serotonin 2 A receptor (HTR2A) and the nitric oxide synthase 1 (NOS1) and environmental factors in patients who experience MDD and SA.

METHODS

A total of 334 patients with MDD and a history of SA (MDD-SA) were recruited alongside 518 patients with MDD with no history of SA (MDD-NSA), and 716 healthy controls (HC). The demographic data and clinical characteristics were collected. Sequenom mass spectrometry was used to detect eight tag-single nucleotide polymorphisms (tagSNPs) in HTR2A (rs1328683, rs17068986, and rs3125) and NOS1 (rs1123425, rs2682826, rs3741476, rs527590, and rs7959232). Generalized multifactor dimensionality reduction (GMDR) was used to analyze the gene-environment interactions.

RESULTS

Four tagSNPs (rs17068986, rs3125, rs527590, and rs7959232) exhibited significant differences between the three groups. However, these differences were not significant between the MDD-SA and MDD-NSA groups after Bonferroni correction. A logistic regression analysis revealed that negative life events (OR = 1.495, 95%CI: 1.071-2.087, P = 0.018), self-guilt (OR = 2.263, 95%CI: 1.515-3.379, P < 0.001), and negative cognition (OR = 2.252, 95%CI: 1.264-4.013, P = 0.006) were all independently associated with SA in patients with MDD. Furthermore, GMDR analysis indicated a significant interaction between HTR2A rs3125 and negative life events. Negative life events in conjunction with the HTR2A rs3125 CG + GG genotype were associated with a higher SA risk in patients with MDD when compared to the absence of negative life events in conjunction with the CC genotype (OR = 2.547, 95% CI: 1.264-5.131, P = 0.009).

CONCLUSION

Several risk factors and a potential interaction between HTR2A rs3125 and negative life events were identified in patients with SA and MDD. The observed interaction likely modulates the risk of MDD and SA, shedding light on the pathogenesis of SA in patients with MDD.

Searching for gene-gene interactions through variance quantitative trait loci of 29 continuous Taiwan Biobank phenotypes

Introduction: After the era of genome-wide association studies (GWAS), thousands of genetic variants have been identified to exhibit main effects on human phenotypes. The next critical issue would be to explore the interplay between genes, the so-called "gene-gene interactions" (GxG) or epistasis. An exhaustive search for all single-nucleotide polymorphism (SNP) pairs is not recommended because this will induce a harsh penalty of multiple testing. Limiting the search of epistasis on SNPs reported by previous GWAS may miss essential interactions between SNPs without significant marginal effects. Moreover, most methods are computationally intensive and can be challenging to implement genome-wide. Methods: I here searched for GxG through variance quantitative trait loci (vQTLs) of 29 continuous Taiwan Biobank (TWB) phenotypes. A discovery cohort of 86,536 and a replication cohort of 25,460 TWB individuals were analyzed, respectively. Results: A total of 18 nearly independent vQTLs with linkage disequilibrium measure r 2 < 0.01 were identified and replicated from nine phenotypes. 15 significant GxG were found with p-values <1.1E-5 (in the discovery cohort) and false discovery rates <2% (in the replication cohort). Among these 15 GxG, 11 were detected for blood traits including red blood cells, hemoglobin, and hematocrit; 2 for total bilirubin; 1 for fasting glucose; and 1 for total cholesterol (TCHO). All GxG were observed for gene pairs on the same chromosome, except for the APOA5 (chromosome 11)-TOMM40 (chromosome 19) interaction for TCHO. Discussion: This study provided a computationally feasible way to search for GxG genome-wide and applied this approach to 29 phenotypes.

Correlations of SDF-1ɑ and XRCC1 gene polymorphisms with the risk of renal cancer development and bioinformatics studies of SDF-1α and XRCC1 and the prognosis of renal cancer

To study the relationships between stromal cell-derived factor-1 (SDF-1ɑ) and renal cell carcinoma (RCC) susceptibility and the presence of single nucleotide polymorphisms in the human X-ray cross-complementary repair gene (XRCC1). Compare SDF-1 based on RCC related data in the TCGA database α, The expression difference of XRCC1 between RCC tissue and normal tissue; Collect 166 newly diagnosed RCC cases and 166 healthy individuals who underwent physical examinations during the same period, and detect genotype using iMLDR method. The results The rs1801157 locus (C:T) of the SDF-1α gene was not significantly associated with the pathohistological type, the rs1799782 locus (G:A) of the XRCC1 gene was associated with the pathohistological type of RCC, and there were interactions between rs1799782 and smoking, alcohol consumption, pesticide exposure, hair dye, and urine holding. The rs1799782 locus of the XRCC1 gene may be a key factor in the pathogenesis and pathological development of RCC. High SDF-1ɑ expression is a protective factor for the overall survival of patients with RCC, and SDF-1ɑ and XRCC1 may be important for the treatment of RCC.

Genetic Association Between ICAM-1 Gene Variants and Susceptibility to Ischemic Cardiomyopathy

ABSTRACT

The current work was aimed at exploring the association between single nucleotide polymorphisms (SNPs) in the ICAM-1 gene, along with the identification of additional haplotypes and their potential role in the susceptibility to ischemic cardiomyopathy (ICM). The control group underwent a Hardy-Weinberg equilibrium test. The associations of genotypes and alleles with susceptibility to ICM were then analyzed using logistic regression analysis. Subsequently odds ratios (ORs) along with 95% confidence intervals (95% CI) were calculated. Interaction analysis was conducted between these SNPs. Furthermore, linkage disequilibrium analysis and haplotype analysis were performed on SNPs that showed interactions with each other. The incidence of ICM was significantly higher among individuals carrying the T allele of rs3093032 (OR = 2.032, 95% CI, 1.275-3.241, P = 0.003) relative to those with the C allele. In addition, CT genotype carriers had a higher susceptibility to ICM than CC genotype carriers (OR = 2.490, 95% CI, 1.445-4.29, P = 0.001). Furthermore, 3 SNPs (rs3093032, rs923366, rs3093030) exhibited a strong interaction with each other, whereas rs281437 showed no interaction with the other 3 SNPs. Individuals carrying the C rs3093032 -T rs923366 -C rs3093030 haplotype had an elevated risk of ICM compared with those carrying the C rs3093032 -C rs923366 -C rs3093030 haplotype (OR = 2.280, 95% CI, 1.568-3.315, P < 0.001). Moreover, individuals carrying the T rs3093032 -C rs923366 -C rs3093030 haplotype were more susceptible to ICM than those carrying the C rs3093032 -C rs923366 -C rs3093030 haplotype (OR = 2.388, 95% CI, 1.469-3.880, P < 0.001). Regarding rs3093032, the minor alleles and haplotypes are associated with an increased ICM risk: 3 SNPs (rs3093032, rs923366, rs3093030) in ICAM-1 have strong interaction with each other.

Interaction between CYP1A1 gene polymorphism and environment factors on risk of endometrial cancer

BACKGROUND

The purpose of this study was to investigate the impact of single nucleotide polymorphisms (SNPs) of the CYP1A1 gene and the gene-environment interaction on the susceptibility to endometrial cancer in Chinese women.

METHOD

Logistic regression was performed to investigate the association between the four SNPs of the CYP1A1 gene and the risk of endometrial cancer. Generalized multifactor dimensionality reduction (GMDR) was employed to analyze the gene-environmental interaction.

RESULTS

A total of 934 women with a mean age of 61.7 ± 10.5 years were selected, including 310 endometrial cancer patients and 624 normal controls. The frequency of rs4646421- T allele was higher in endometrial cancer patients than normal controls, the T allele of rs4646421 was 28.1% in endometrial cancer patients and 21.0% in normal controls (p < 0.001). Logistic regression analysis showed that the rs4646421 - T allele was associated with increased risk of endometrial cancer, OR (95% CI) were 1.52 (1.11-1.97) and 1.91 (1.35-2.52), respectively. GMDR analysis found a significant two-locus model (p = 0.0107) involving rs4646421 and abdominal obesity (defined by waist circumference), indicating a potential gene-environment interaction between rs4646421 and abdominal obesity. Abdominal obese subjects with rs4646421- CT or TT genotype have the highest risk of endometrial cancer, compared to non-abdominal obese subjects with the rs4646421- CC genotype, the OR (95%CI) was 2.23 (1.62-2.91).

CONCLUSIONS

Both the rs4646421- T allele and the interaction between rs4646421 and abdominal obesity were associated with increased risk of endometrial cancer.

Molecular role of non-exonic variants in CALPAIN 10 gene in polycystic ovarian syndrome in Saudi women

Introduction

Non-diabetic women with polycystic ovarian syndrome (PCOS) often have abnormal insulin regulation. Calpain 10 (CALP10) is a biomarker of type 2 diabetes mellitus, with some of its single-nucleotide polymorphisms (SNPs) influencing PCOS development.

Methods

In this case-control study on 90 women each with and without PCOS, we explored the molecular role of five CALP10 SNPs using biochemical parameters and Sanger sequencing analyses.

Results

Different genetic models, genotypes, and allele frequencies were significantly associated with UCSNP-19 (rs3842570; p=0.01), UCSNP-44 (rs2975760; p=0.009), UCSNP-56 (rs2975762; p<0.0001), and UCSNP-63 (rs5030952; p=0.0003) in women with PCOS. The multiple logistic regression model showed a strong association of CALP10 SNPs with fasting blood glucose (p<0.001). ANOVA showed significant associations with various biochemical parameters such as FSH (p=0.0001) in UCSNP-19 (rs3842570), FI (p=0.002), TG (p=0.01) in UCSNP-56 (rs2975762) and FBG (p=0.001), FI (p=0.004), FSH (p=0.02) & LDLc (p=0.04) in UCSNP-63 (rs5030952) SNPs. Haplotype analysis also revealed significant associations between different combinations of alleles in the studied 5 SNPs in women with PCOS (p<0.05). Generalized multifactor dimensionality reduction analysis showed the best gene-gene interactions among the five SNPs in CALP10I (p<0.05). However, dendrogram and graphical depletion models found no strong association in women with PCOS.

Conclusion

In conclusion, this study confirms rs3842570, rs2975760, rs2975767, and rs5030952 SNPs in CALP10 gene is associated in diagnosed PCOS women in the Saudi Arabia.

Genetic Susceptibility to Tardive Dyskinesia and Cognitive Impairments in Chinese Han Schizophrenia: Role of Oxidative Stress-Related and Adenosine Receptor Genes

Objective

Tardive dyskinesia (TD) is a severe rhythmic movement disorder caused by long-term antipsychotic medication in chronic patients with schizophrenia (SCZ). We aimed to investigate the association between polymorphisms in oxidative stress-related genes (GSTM1, SOD2, NOS1, and NOS3) and adenosine receptor gene (ADORA2A), as well as their interactions, with the occurrence and severity of TD, and cognitive impairments in a Chinese Han population of SCZ patients.

Methods

Two hundred and sixteen SCZ patients were recruited and divided into TD group (n=157) and non-TD group (n=59). DNA extraction was performed by a high-salt method, followed by SNP genotyping using matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS). The severity of TD, psychopathology and cognitive functioning were assessed using the Abnormal Involuntary Movement Scale (AIMS), the Positive and Negative Syndrome Scale (PANSS) and the Repeated Battery for Assessment of Neuropsychological Status (RBANS), respectively.

Results

The combination of GSTM1-rs738491, NOS1-rs738409 and ADORA2A-rs229883 was identified as the best three-point model to predict TD occurrence (p=0.01). Additionally, GSTM-rs738491 CC or NOS3-rs1800779 AG genotypes may be protective factors for psychiatric symptoms in TD patients. TD patients carrying the NOS1-rs738409 AG or ADORA2A-rs229883 TT genotypes exhibited poorer cognitive performance.

Conclusion

Our findings suggest that the interaction of oxidative stress-related genes and adenosine receptor gene may play a role in the susceptibility and severity of TD in Chinese Han SCZ patient. Furthermore, these genes may also affect the psychiatric symptoms and cognitive function of TD patients.

Gene polymorphisms associated with sudden decreases in heart rate during extensive peritoneal lavage with distilled water after gastrectomy

BACKGROUND

Our previous study found that the telomerase-associated protein 1 (TEP1, rs938886 and rs1713449) and homo sapiens RecQ like helicase 5 (RECQL5, rs820196) single nucleotide polymorphisms (SNPs) were associated with changes in heart rate (HR) ≥ 30% during peritoneal lavage with distilled water after gastrectomy. This study established a single tube method for detecting these three SNPs using two-dimensional (2D) polymerase chain reaction (PCR), and investigated whether SNP-SNP and SNP-environment interactions increase the risk of high HR variability (HRV).

AIM

To investigate whether genotypes, genetic patterns, SNP-SNP and SNP-environment interactions were associated with HRV.

METHODS

2D PCR was used to establish a single-tube method to detect TEP1 rs938886 and rs1713449 and RECQL5 rs820196, and the results were compared with those of sanger sequencing. After adjusting for confounders such as age, sex, smoking, hypertension, and thyroid dysfunction, a nonconditional logistic regression model was used to assess the associations between the genotypes and the genetic patterns (codominant, dominant, overdominant, recessive, and additive) of the three SNPs and a risk ≥ 15% or ≥ 30% of a sudden drop in HR during postoperative peritoneal lavage in patients with gastric cancer. Gene-gene and gene-environment interactions were analyzed using generalized multifactor dimensionality reduction.

RESULTS

The coincidence rate between the 2D PCR and sequencing was 100%. When the HRV cutoff value was 15%, the patients with the RECQL5 (rs820196) TC genotype had a higher risk of high HRV than those who had the TT genotype (odds ratio = 1.97; 95%CI: 1.05-3.70; P = 0.045). Under the codominant and overdominant models, the TC genotype of RECQL5 (rs820196) was associated with a higher risk of HR decrease relative to the TT and TT + CC genotypes (P = 0.031 and 0.016, respectively). When the HRV cutoff value was 30%, patients carrying the GC-TC genotypes of rs938886 and rs820196 showed a higher HRV risk when compared with the GG-TT genotype carriers (P = 0.01). In the three-factor model of rs938886, rs820196, and rs1713449, patients carrying the GC-TC-CT genotype had a higher risk of HRV compared with the wild-type GG-TT-CC carriers (P = 0.01). For rs820196, nonsmokers with the TC genotype had a higher HRV risk compared with nonsmokers carrying the TT genotype (P = 0.04). When the HRV cutoff value was 15%, patients carrying the TT-TT and the TC-CT genotypes of rs820196 and rs1713449 showed a higher HRV risk when compared with TT-CC genotype carriers (P = 0.04 and 0.01, respectively). Patients carrying the GC-CT-TC genotypes of rs938886, rs1713449, and rs820196 showed a higher HRV risk compared with GG-CC-TT genotype carriers (P = 0.02). When the HRV cutoff value was 15%, the best-fitting models for the interactions between the SNPs and the environment were the rs820196-smoking (P = 0.022) and rs820196-hypertension (P = 0.043) models. Consistent with the results of the previous grouping, for rs820196, the TC genotype nonsmokers had a higher HRV risk compared with nonsmokers carrying the TT genotype (P = 0.01).

CONCLUSION

The polymorphism of the RECQL5 and TEP1 genes were associated with HRV during peritoneal lavage with distilled water after gastrectomy.

Interactive effects between CDHR3 genotype and rhinovirus species for diagnosis and severity of respiratory tract infections in hospitalized children

Rhinovirus (RV) is the leading pathogen causing childhood wheezing, with rhinovirus C (RV-C) species reported to cause asthma exacerbation. Allele A of single-nucleotide polymorphism (SNP) CDHR3_rs6967330 upregulates epithelial expression of RV-C receptors which results in more severe asthma exacerbations in children. Nevertheless, there are limited data on interactions between CDHR3 variants and their impact on severity of RV-related pediatric respiratory tract infections (RTIs). Medical records of RV-related RTIs in children aged below 18 years who were hospitalized in two public hospitals in 2015-2016 were independently reviewed by two paediatricians. Archived nasopharyngeal aspirates were retrieved for RV detection and sequencing as well as CDHR3 genotyping. HaploView v.5.0 and generalized multifactor dimensionality reduction (GMDR) analysis were employed for haplotypic assignment and gene-environment interaction analyses. Among 1019 studied cases, our results confirmed the relationship between RV-C species and more severe RTIs. Besides the top risk variant rs6967330-A, we identified rs140154310-T to be associated with RV-C susceptibility under the additive model [odds ratio (OR) 2.53, 95% CI 1.15-5.56; P = 0.021]. Rs140154310 was associated with wheezing illness (OR 2.38, 95% CI 1.12-5.04; P = 0.024), with such association being stronger in subjects who wheezed due to RV-C infections (OR 2.71, 95% CI 1.32-5.58; P = 0.007). Haplotype GAG constructed from rs4730125, rs6967330, and rs73195665 was associated with increased risk of RV-C infection (OR 1.71, 95% CI 1.11-2.65; P = 0.016) and oxygen supplementation (OR 1.93, 95% CI 1.13-3.30; P = 0.016). GMDR analyses revealed epistatic interaction between rs140154310 and rs6967330 of CDHR3 for RV-C infection (P = 0.001), RV-C-associated lower RTI (P = 0.004), and RV-C-associated wheeze (P = 0.007). There was synergistic gene-environmental interaction between rs3887998 and RV-C for more severe clinical outcomes (P < 0.001). To conclude, rs140154310-T is another risk variant for RV-C susceptibility and more severe RTIs. Synergistic epistatic interaction is found between CDHR3 SNPs and RV-C for RTI severity, which is likely mediated by susceptibility to RV-C. Haplotypic analysis and GMDR should be included in identifying prediction models of CDHR3 for childhood asthma and RTIs. IMPORTANCE This case-control study investigated the interaction between CDHR3 genotypes and rhinovirus (RV) species on disease severity in Hong Kong children hospitalized for respiratory tract infection (RTI). There were synergistic effects between RV-C and CDHR3 SNPs for RTI severity, which was mainly driven by RV-C. Specifically, rs6967330 and rs140154310 alone and their epistatic interaction were associated with RV-C-related and severe RTIs in our subjects. Therefore, genotyping of CDHR3 SNPs may help physicians formulate prediction models for severity of RV-associated RTIs.

Molecular Effect of Variants in Toll-like Receptor 4 Gene in Saudi Patients with Type 2 Diabetes Mellitus

Single-nucleotide polymorphisms (SNPs) in the Toll-like receptor 4 (TLR4) gene have been documented in type 2 diabetes mellitus (T2DM) and other diseases in the Saudi population. We investigated the relationship between rs11536889, rs4986790, and rs4986791 SNPs in the TLR4 gene and T2DM in the Saudi population; 105 patients with T2DM and 105 healthy controls were analyzed. The TLR4 gene was amplified through PCR, followed by restriction fragment length polymorphism analysis for rs4986791 and Sanger sequencing for rs11536889 and rs4986790 SNPs. The clinical and biochemical characteristics were associated with T2DM (p < 0.05). The rs11536889, rs4986790, and rs4986791 SNPs in control subjects followed the Hardy-Weinberg equilibrium (p > 0.05). Alleles were associated with rs11536889, rs4986791, heterozygous codominant, and dominant models (p < 0.05). However, the rs4986790 SNP was not associated with T2DM (p > 0.05). Logistic regression analysis showed that high-density lipoprotein cholesterol (HDLc) levels were associated with T2DM (p < 0.001). Analysis of variance showed that waist (p = 0.0005) and hip circumferences (p = 0.002) in rs4986790 and rs4986791 SNPs, in SBP (p = 0.001), DBP (p = 0.002), and HDLc levels (p = 0.003), were associated with T2DM subjects. T2DM was also associated with the haplotype (p < 0.001) but not with linkage disequilibrium. The gene-gene interaction was associated with the three SNPs studied in patients with T2DM according to the generalized multifactor dimensionality reduction model (p < 0.0001). Dendrogram and graphical depletion analysis revealed a moderate association in patients with T2DM. The results suggest that rs11536889 and rs4986790 SNPs are genotypically and allelically associated with T2DM in Saudi patients. Future functional studies are recommended to validate the genetic roles of these SNPs in the pathogenesis and progression of diseases.

SAMM50-rs2073082, -rs738491 and -rs3761472 Interactions Enhancement of Susceptibility to Non-Alcoholic Fatty Liver Disease

BACKGROUND AND AIM

Several studies have identified that three SAMM50 polymorphisms (rs2073082, rs738491, rs3761472) are associated with an increased risk of non-alcoholic fatty liver disease (NAFLD). However, the clinical significance of the SAMM50 SNP in relation to NAFLD remains largely unknown. Therefore, we conducted a clinical study and SNP-SNP interaction analysis to further elucidate the effect of the SAMM50 SNP on the progression of NAFLD in the elderly.

METHODS

A total of 1053 patients over the age of 65 years were recruited. Liver fat and fibrosis were detected by abdominal ultrasound or FibroScan, respectively. Genomic DNA was extracted and then genotyped by Fluidigm 96.96 Dynamic Array. Multivariable logistic regression was used to evaluate the association between NAFLD and SNP. SNP-SNP interactions were analyzed using generalized multivariate dimensionality reduction (GMDR).

RESULTS

The risk of NAFLD was substantially higher in people who carried SAMM50-rs2073082 G and -rs738491 T alleles (OR, 1.962; 95% CI, 1.448-2.659; p < 0.001; OR, 1.532; 95% CI, 1.246-1.884; p = 0.021, respectively) compared to noncarriers. Carriers of the rs738491 T and rs3761472 G alleles in the cohort showed a significant increase in liver stiffness measurements (LSM). The combination of the three SNPs showed the highest predictive power for NAFLD. The rs2073082 G allele, rs738491 T allele and rs3761472 G carriers had a two-fold higher risk of NAFLD compared to noncarriers.

CONCLUSIONS

Our research has demonstrated a strong correlation between the genetic polymorphism of SAMM50 and NAFLD in the elderly, which will contribute to a better understanding of the impact of age and genetics on this condition. Additionally, this study provides a potential predictive model for the early clinical warning of NAFLD.

Contribution of gene polymorphisms on 3p25 to salivary gland carcinoma, ameloblastoma, and odontogenic keratocyst in the Chinese Han population

OBJECTIVE

This study aimed to investigate the contribution of gene polymorphisms in 3p25 to salivary gland carcinoma (SGC), ameloblastoma (AM), and odontogenic keratocyst (OKC) in the Chinese Han population.

STUDY DESIGN

Sixteen tag-single nucleotide polymorphisms (SNPs) within 5 genes (SYN2, TIMP4, PPARG, RAF1, and IQSEC1) in 3p25 were genotyped in 411 individuals with or without SGC, AM, and OKC. Genotype, clinical phenotype, and bioinformatics analyses were performed to evaluate the function of candidate SNPs.

RESULTS

SYN2-rs3773364, TIMP4-rs3755724, PPARG-rs10865710, and PPARG-rs1175544 were related to decreased SGC susceptibility, whereas IQSEC1-rs2600322 and IQSEC1-rs2686742 decreased and increased AM risk, respectively. Stratification analysis revealed that the significance of the identified SNPs was stronger in females or individuals younger than 46 years in SGC. PPARG-rs10865710 and PPARG-rs1175544 were associated with lower lymph node metastasis. SYN2-rs3773364 and PPARG-rs1175544 were associated with favorable SGC patient survival. Functional assessments linked PPARG-rs1175544 to PPARG expression regulation. Linkage disequilibrium analysis revealed a haplotype (SYN2-rs3773364-A, TIMP4-rs3817004-A, and TIMP4-rs3755724-C) associated with decreased susceptibility to SGC. Generalized multifactor dimensionality reduction analysis indicated the gene-gene interactions among IQSEC1, TIMP4, and PPARG in SGC, AM, and OKC progression.

CONCLUSIONS

These variants play important roles in the progression of SGC, AM, and OKC in the Chinese Han population and may be considered biomarkers for early diagnosis and prognosis prediction.

Variants in genes related to inflammation and endothelial function can increase the risk for carotid atherosclerosis in southwestern China

Aim

To investigate the potential association between polymorphisms in genes involved in endothelial function, inflammation and carotid atherosclerosis.

Methods

This was a three-center, population-based sectional survey conducted in Sichuan province of southwestern China. We randomly selected 8 different communities in Sichuan, and the residents in each community volunteered to participate in the survey by face-to-face questionnaire. A total of 2,377 residents with high stroke risk population in the 8 communities were included. Carotid atherosclerosis was evaluated by carotid ultrasound, and the 19 single nucleotide polymorphisms (SNPs) in 10 endothelial function as well as inflammation relevant genes were measured in the high stroke risk population. Carotid atherosclerosis was defined by the presence of carotid plaque or any carotid stenosis ≥15% or mean intima-media thickness (IMT) > 0.9 mm. Generalized multifactor dimensionality reduction (GMDR) approach was used to analyze gene-gene interactions among the 19 SNPs.

Results

Among the 2,377 subjects with high stroke risk, 1,028 subjects had carotid atherosclerosis (43.2%), of which 852 (35.8%) cases had carotid plaque, 295 (12.4%) cases had ≥15% carotid stenosis, whereas 445 (18.7%) had mean IMT > 0.9 mm. Multivariate logistic regression revealed that IL1A rs1609682 TT and HABP2 rs7923349 TT served as independent risk factors for carotid atherosclerosis (OR, 1.45, 95% CI: 1.034-2.032, p = 0.031, and OR, 1.829, 95% CI: 1.228-2.723, p = 0.003). GMDR analysis indicated that there was a significant gene-gene interaction found among IL1A rs1609682, ITGA2 rs1991013, and HABP2 rs7923349. After adjusting the covariates, the high-risk interactive genotypes in the 3 variants were significantly associated with a significantly higher risk for carotid atherosclerosis (OR, 2.08, 95% CI: 1.257-5.98, p < 0.001).

Conclusion

The prevalence of carotid atherosclerosis was observed to be extremely high in the high-risk stroke population in southwestern China. There were associations observed between the specific variants in inflammation and endothelial function relevant genes and carotid atherosclerosis. The high-risk interactive genotypes among IL1A rs1609682, ITGA2 rs1991013, and HABP2 rs7923349 significantly increased the risk of carotid atherosclerosis. These results are expected to provide novel strategies for the prevention of carotid atherosclerosis. The gene-gene interactive analysis used in this study may be very helpful to elucidate complex genetic risk factors for carotid atherosclerosis.

Inflammation and endothelial function-related gene polymorphisms are associated with carotid atherosclerosis-A study of community population in Southwest China

OBJECTIVES

To investigate the relationships between 18 single nucleotide polymorphisms with carotid atherosclerosis and whether interactions among these genes were associated with an increased risk of carotid atherosclerosis.

METHODS

Face-to-face surveys were conducted with individuals aged 40 or older in eight communities. A total of 2377 individuals were included in the study. Ultrasound was used to detect carotid atherosclerosis in the included population. 18 loci of 10 genes associated with inflammation and endothelial function were detected. Gene-gene interactions were analyzed using generalized multifactor dimensionality reduction (GMDR).

RESULTS

Among the 2377 subjects, 445 (18.7%) subjects had increased intima-media thickness in the common carotid artery (CCA-IMT), and 398 (16.7%) subjects were detected with vulnerable plaque. In addition, NOS2A rs2297518 polymorphism was associated with increased CCA-IMT, IL1A rs1609682, and HABP2 rs7923349 polymorphisms were associated with vulnerable plaque. Besides, GMDR analysis showed significant gene-gene interactions among TNFSF4 rs1234313, IL1A rs1609682, TLR4 rs1927911, ITGA2 rs1991013, NOS2A rs2297518, IL6R rs4845625, ITGA2 rs4865756, HABP2 rs7923349, NOS2A rs8081248, HABP2 rs932650.

CONCLUSION

The prevalences of increased CCA-IMT and vulnerable plaque were high in Southwestern China's high-risk stroke population. Furthermore, inflammation and endothelial function-related gene polymorphisms were associated with carotid atherosclerosis.

Multifactor dimensionality reduction method identifies novel SNP interactions in the WNT protein interaction networks that are associated with recurrence risk in colorectal cancer

Background

Interactions among genetic variants are rarely studied but may explain a part of the variability in patient outcomes.

Objectives

In this study, we aimed to identify 1 to 3 way interactions among SNPs from five Wnt protein interaction networks that predict the 5-year recurrence risk in a cohort of stage I-III colorectal cancer patients.

Methods

423 patients recruited to the Newfoundland Familial Colorectal Cancer Registry were included. Five Wnt family member proteins (Wnt1, Wnt2, Wnt5a, Wnt5b, and Wnt11) were selected. The BioGRID database was used to identify the proteins interacting with each of these proteins. Genotypes of the SNPs located in the interaction network genes were retrieved from a genome-wide SNP genotype data previously obtained in the patient cohort. The GMDR 0.9 program was utilized to examine 1-, 2-, and 3-SNP interactions using a 5-fold cross validation step. Top GMDR 0.9 models were assessed by permutation testing and, if significant, prognostic associations were verified by multivariable logistic regression models.

Results

GMDR 0.9 has identified novel 1, 2, and 3-way SNP interactions associated with 5-year recurrence risk in colorectal cancer. Nine of these interactions were multi loci interactions (2-way or 3-way). Identified interaction models were able to distinguish patients based on their 5-year recurrence-free status in multivariable regression models. The significance of interactions was the highest in the 3-SNP models. Several of the identified SNPs were eQTLs, indicating potential biological roles of the genes they were associated with in colorectal cancer recurrence.

Conclusions

We identified novel interacting genetic variants that associate with 5-year recurrence risk in colorectal cancer. A significant portion of the genes identified were previously linked to colorectal cancer pathogenesis or progression. These variants and genes are of interest for future functional and prognostic studies. Our results provide further evidence for the utility of GMDR models in identifying novel prognostic biomarkers and the biological importance of the Wnt pathways in colorectal cancer.

MTOR gene polymorphism may be associated with microscopic polyangiitis susceptibility in a Guangxi population of China

OBJECTIVE

Microscopic polyangiitis (MPA) onset is affected by genetic predisposition. Autophagy plays a certain role in antineutrophil cytoplasmic antibody-associated vasculitis developing. A key factor in autophagy regulating, the genetic polymorphism of MTOR gene is essential. The objective was to explore the associations between MTOR gene polymorphism and MPA susceptibility in a Guangxi population of China.

METHODS

A sum of 208 MPA cases and 209 healthy volunteers from Guangxi in this case-control study, four important single nucleotide polymorphism (SNP) loci of MTOR gene including rs3806317, rs1064261, rs1883965 and rs2295080 were examined. Multiplex polymerase chain reaction combined with high-throughput sequencing was performed. Subgroup analysis was evaluated by gender and ethnicity. Linkage disequilibrium and haplotype analysis were tested. Multi-SNPs interaction among mTOR signaling pathway was assessed.

RESULTS

For rs2295080, homozygous mutant GG genotype was associated with a decreased susceptibility of MPA in recessive model (OR = 0.38, 95%CI: 0.14-1.00, p = 0.040), particularly in the subgroup of female (OR = 0.16, 95%CI: 0.03-0.74, p = 0.006) and Han population (OR = 0.32, 95%CI: 0.10-1.00, p = 0.034). Individual carrying G allele was linked with decreasing MPA susceptibility in Han population of Guangxi (OR = 0.65, 95%CI: 0.44-0.97, p = 0.036). In haplotype analysis, the haplotype AAT was correlated with increasing susceptibility of MPA (OR = 1.347, 95%CI: 1.004-1.807, p = 0.046). Moreover, in the multi-SNPs interaction analysis, the six-locus model was identified as the best interaction model (p < 0.05).

CONCLUSION

These findings suggest that rs2295080 polymorphism of MTOR gene may be associated with MPA susceptibility in a Guangxi population of China and G allele might be an important protective factor.

Evidence for Epistatic Interaction between HLA-G and LILRB1 in the Pathogenesis of Nonsegmental Vitiligo

Vitiligo is the most frequent cause of depigmentation worldwide. Genetic association studies have discovered about 50 loci associated with disease, many with immunological functions. Among them is HLA-G, which modulates immunity by interacting with specific inhibitory receptors, mainly LILRB1 and LILRB2. Here we investigated the LILRB1 and LILRB2 association with vitiligo risk and evaluated the possible role of interactions between HLA-G and its receptors in this pathogenesis. We tested the association of the polymorphisms of HLA-G, LILRB1, and LILRB2 with vitiligo using logistic regression along with adjustment by ancestry. Further, methods based on the multifactor dimensionality reduction (MDR) approach (MDR v.3.0.2, GMDR v.0.9, and MB-MDR) were used to detect potential epistatic interactions between polymorphisms from the three genes. An interaction involving rs9380142 and rs2114511 polymorphisms was identified by all methods used. The polymorphism rs9380142 is an HLA-G 3'UTR variant (+3187) with a well-established role in mRNA stability. The polymorphism rs2114511 is located in the exonic region of LILRB1. Although no association involving this SNP has been reported, ChIP-Seq experiments have identified this position as an EBF1 binding site. These results highlight the role of an epistatic interaction between HLA-G and LILRB1 in vitiligo pathogenesis.

Inflammation and endothelial function relevant genetic polymorphisms in carotid stenosis in southwestern China

Aim

To evaluate the relationship between carotid stenosis with variants in genes referred to inflammation and endothelial function.

Methods

There was a multi-center, cross sectional survey in southwestern China. The eight communities were selected at random in southwestern China. The residents aged ≥40 years volunteered to participate in face-to-face survey. Subjects with at least three of the aforementioned eight stroke related risk factors or a history of stroke were classified as high-risk population for stroke. A total of 2,377 subjects were the high-risk population for stroke in the eight communities, and degree of carotid stenosis was assessed by carotid ultrasound. Genotypes of 6 variants in 3 genes related to inflammation and endothelial function were examined. Gene-gene interaction was analyzed by generalized multifactor dimensionality reduction (GMDR).

Results

Carotid stenosis were found in 295 (12.41%) subjects, of whom 51 (17.29%) had moderate or severe stenosis. According to multivariate logistic regression analysis, we found that HABP2rs7923349TT was independent risk factor for carotid stenosis (OR, 1.96, 95% CI: 1.22-3.13, P = 0.005) and ITGA2rs1991013AG and HABP2rs7923349TT were independent risk factors for moderate to severe carotid stenosis (OR, 2.28, 95% CI: 1.28-4.07, P = 0.005; OR, 2.90, 95% CI: 1.19-7.08, P = 0.019). GMDR analysis showed that there was a significant gene-gene interaction between ITGA2 rs4865756 and HABP2 rs7923349, and the high-risk interactive genotype in the two variants was independently associated with a higher risk for carotid stenosis after adjusting the covariates (OR,1. 42, 95% CI 1.10-1.84, P = 0.008).

Conclusions

Prevalence of carotid stenosis was very high in the high-risk stroke population in southwestern China. Variants in genes referred in endothelial function were associated with the carotid stenosis. The high-risk interactive genotype in ITGA2 rs4865756 and HABP2 rs7923349 was independently associated with a higher risk for carotid stenosis.

Overview of frequent pattern mining

Various methods of frequent pattern mining have been applied to genetic problems, specifically, to the combined association of two genotypes (a genotype pattern, or diplotype) at different DNA variants with disease. These methods have the ability to come up with a selection of genotype patterns that are more common in affected than unaffected individuals, and the assessment of statistical significance for these selected patterns poses some unique problems, which are briefly outlined here.

Association of ATG7 gene polymorphisms with microscopic polyangiitis in Chinese individuals

Microscopic polyangiitis (MPA) is a type of antineutrophil cytoplasmic antibody (ANCA)-related vasculitis. Autophagy-related gene 7 (ATG7) protects against complicated disorder states in model organisms, but the way ATG7 dysfunction contributes to MPA remains elusive. This investigation assessed the impacts of ATG7 single-nucleotide polymorphisms (SNPs) on microscopic polyangiitis (MPA) in China. A total of 211 controls and 214 MPA patients were recruited and analyzed. Polymerase chain reaction (PCR) and high-throughput sequencing were adopted to detect the ATG7 SNPs (rs75492008, rs2594966, rs6442260 and rs8154), and stratification analysis, different genetic models and differences in allele and genotype frequencies were evaluated. Haplotype evaluation was performed after linkage disequilibrium (LD) analyses, and interactions between alleles were assessed. Generalized multifactor dimensionality reduction (GMDR) was adopted to analyze SNP-SNP interactions among the four ATG7 SNPs and phosphatidylinositol-4, 5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and unc-nc-like autophagy activating kinase 1 (ULK1) SNPs previously studied by our team. Relationships between ATG7 polymorphisms, disease activity biomarkers and therapeutic effects in MPA were analyzed. Sex stratification analysis of the rs2594966 GG genotype with codominant and recessive models showed OR=3.42, 95% CI [1.19-9.80], P=0.041 and OR=3.31, 95% CI [1.23-8.90], P=0.012, respectively. Haplotype G-G-C-T was related to an increased MPA risk (OR=1.5, 95% CI [0.999-2.266], P=0.029). Permutation testing of GMDR suggested that ATG7 rs6442260 and rs8154, PIK3CA rs1607237, and ULK1 rs4964879 might interact with each other in MPA development (P<0.05). Among 214 MPA patients, 79 available complete follow-up clinical datasets were gathered from September 2009 to October 2020, showing that rs75492008 and rs4964879 affect the correlation between C-reactive protein (CRP) and the erythrocyte sedimentation rate (ESR) in MPA activity. Patients with rs8154 TT and rs1607237 CC genotypes had better clinical treatment effects (P<0.05). Gene polymorphisms may be related to MPA in China, exhibiting correlation with MPA activity indicators, treatment and prognosis.

Examining SNP-SNP interactions and risk of clinical outcomes in colorectal cancer using multifactor dimensionality reduction based methods

Background: SNP interactions may explain the variable outcome risk among colorectal cancer patients. Examining SNP interactions is challenging, especially with large datasets. Multifactor Dimensionality Reduction (MDR)-based programs may address this problem.

Objectives: 1) To compare two MDR-based programs for their utility; and 2) to apply these programs to sets of MMP and VEGF-family gene SNPs in order to examine their interactions in relation to colorectal cancer survival outcomes.

Methods: This study applied two data reduction methods, Cox-MDR and GMDR 0.9, to study one to three way SNP interactions. Both programs were run using a 5-fold cross validation step and the top models were verified by permutation testing. Prognostic associations of the SNP interactions were verified using multivariable regression methods. Eight datasets, including SNPs from MMP family genes (n = 201) and seven sets of VEGF-family interaction networks (n = 1,517 SNPs) were examined.

Results: ∼90 million potential interactions were examined. Analyses in the MMP and VEGF gene family datasets found several novel 1- to 3-way SNP interactions. These interactions were able to distinguish between the patients with different outcome risks (regression p-values 0.03–2.2E-09). The strongest association was detected for a 3-way interaction including CHRM3.rs665159_EPN1.rs6509955_PTGER3.rs1327460 variants.

Conclusion: Our work demonstrates the utility of data reduction methods while identifying potential prognostic markers in colorectal cancer.

A Review of Feature Selection Methods for Machine Learning-Based Disease Risk Prediction

Machine learning has shown utility in detecting patterns within large, unstructured, and complex datasets. One of the promising applications of machine learning is in precision medicine, where disease risk is predicted using patient genetic data. However, creating an accurate prediction model based on genotype data remains challenging due to the so-called “curse of dimensionality” (i.e., extensively larger number of features compared to the number of samples). Therefore, the generalizability of machine learning models benefits from feature selection, which aims to extract only the most “informative” features and remove noisy “non-informative,” irrelevant and redundant features. In this article, we provide a general overview of the different feature selection methods, their advantages, disadvantages, and use cases, focusing on the detection of relevant features (i.e., SNPs) for disease risk prediction.

Association between MIR31HG polymorphisms and the risk of Lumbar disc herniation in Chinese Han population

Lumbar disc herniation (LDH) is a common spinal disease that endangers human health. Genetic factors play a vital role in the progression of LDH. This study aimed to explore the relationship of the MIR31HG polymorphism with LDH risk in the Chinese population. Seven candidate SNPs on MIR31HG in 504 patients with LDH and 503 healthy people were genotyped by Agena MassARRAY platform. Logistic regression was used to calculate the relationship between MIR31HG polymorphism and LDH risk under different genetic models. Multi-factor dimensionality reduction (MDR) analysis was performed to evaluate the SNP-SNP interaction. We found that rs10965059 was significantly associated with a decreased risk of LDH under the dominant (OR = 0.46, 95% CI: 0.34-0.62, P < 0.001), log-additive (OR = 0.59, 95%CI: 0.45-0.76, P < 0.001), and codominant (OR = 0.40, 95%CI: 0.29-0.55, P < 0.001) models in the overall analysis. In the subgroup analyses of age, male, and complications, we found that rs10965059 was associated with a reduced risk of LDH. However, there was no significant correlation between MiR-31HG polymorphisms and risk of LDH in females. In addition, the three SNPs (rs72703442-rs2025327-rs55683539) was mapped to a 26kb LD block with D' >0.96, suggesting a significant linkage disequilibrium presence among each pair SNPs. MDR analysis showed that the best single-locus and multi-locus models for the prediction of LDH risk were rs10965059 and seven-locus models, respectively, and both of them increased LDH risk. Our results shown that in the Chinese Han population, the MIR31HG polymorphism rs10965059 was involved in a risk to symptomatic LDH, which provides a scientific basis for early screening, prevention, diagnosis and treatment of local LDH high-risk populations.

Identification of the associations between genes and quantitative traits using entropy-based kernel density estimation

Genetic associations have been quantified using a number of statistical measures. Entropy-based mutual information may be one of the more direct ways of estimating the association, in the sense that it does not depend on the parametrization. For this purpose, both the entropy and conditional entropy of the phenotype distribution should be obtained. Quantitative traits, however, do not usually allow an exact evaluation of entropy. The estimation of entropy needs a probability density function, which can be approximated by kernel density estimation. We have investigated the proper sequence of procedures for combining the kernel density estimation and entropy estimation with a probability density function in order to calculate mutual information. Genotypes and their interactions were constructed to set the conditions for conditional entropy. Extensive simulation data created using three types of generating functions were analyzed using two different kernels as well as two types of multifactor dimensionality reduction and another probability density approximation method called m-spacing. The statistical power in terms of correct detection rates was compared. Using kernels was found to be most useful when the trait distributions were more complex than simple normal or gamma distributions. A full-scale genomic dataset was explored to identify associations using the 2-h oral glucose tolerance test results and γ-glutamyl transpeptidase levels as phenotypes. Clearly distinguishable single-nucleotide polymorphisms (SNPs) and interacting SNP pairs associated with these phenotypes were found and listed with empirical p-values.

Concurrent outcomes from multiple approaches of epistasis analysis for human body mass index associated loci provide insights into obesity biology

Genome wide association studies (GWAS) have focused on elucidating the genetic architecture of complex traits by assessing single variant effects in additive genetic models, albeit explaining a fraction of the trait heritability. Epistasis has recently emerged as one of the intrinsic mechanisms that could explain part of this missing heritability. We conducted epistasis analysis for genome-wide body mass index (BMI) associated SNPs in Alzheimer’s Disease Neuroimaging Initiative (ADNI) and followed up top significant interacting SNPs for replication in the UK Biobank imputed genotype dataset. We report two pairwise epistatic interactions, between rs2177596 (RHBDD1) and rs17759796 (MAPK1), rs1121980 (FTO) and rs6567160 (MC4R), obtained from a consensus of nine different epistatic approaches. Gene interaction maps and tissue expression profiles constructed for these interacting loci highlights co-expression, co-localisation, physical interaction, genetic interaction, and shared pathways emphasising the neuronal influence in obesity and implicating concerted expression of associated genes in liver, pancreas, and adipose tissues insinuating to metabolic abnormalities characterized by obesity. Detecting epistasis could thus be a promising approach to understand the effect of simultaneously interacting multiple genetic loci in disease aetiology, beyond single locus effects.

Interaction of SAMM50-rs738491, PARVB-rs5764455 and PNPLA3-rs738409 Increases Susceptibility to Nonalcoholic Steatohepatitis

BACKGROUND AND AIMS

Previous studies have reported that the single nucleotide polymorphisms (SNPs) of SAMM50-rs738491, PARVB-rs5764455 and PNPLA3-rs738409 are associated with nonalcoholic fatty liver disease (NAFLD). However, no studies have examined the effect of interactions between these three genotypes to affect liver disease severity. We assessed the effect of these three SNPs on nonalcoholic steatohepatitis (NASH) and also examined the gene-gene interactions in a Chinese population with biopsy-confirmed NAFLD.

METHODS

We enrolled 415 consecutive adult individuals with biopsy-proven NAFLD. Multivariable logistic regression analysis was undertaken to test associations between NASH and SNPs in SAMM50-rs738491, PARVB-rs5764455 and PNPLA3-rs738409. Gene-gene interactions were analyzed by performing a generalized multifactor dimensionality reduction (GMDR) analysis.

RESULTS

The mean ± standard deviation age of these 415 patients was 41.3±12.5 years, and 75.9% were men. Patients with SAMM50-rs738491 TT, PARVB-rs5764455 AA or PNPLA3-rs738409 GG genotypes had a higher risk of NASH, even after adjustment for age, sex and body mass index. GMDR analysis showed that the combination of all three SNPs was the best model for predicting NASH. Additionally, the odds ratio of the haplotype T-A-G for predicting the risk of NASH was nearly three times higher than that of the haplotype G-C-C.

CONCLUSIONS

NAFLD patients carrying the SAMM50-rs738491 TT, PARVB-rs5764455 AA or PNPLA3-rs738409 GG genotypes are at greater risk of NASH. These three SNPs may synergistically interact to increase susceptibility to NASH.

Genetic and Clinical Factors Associated with Opioid Response in Chinese Han Patients with Cancer Pain: An Exploratory Cross-Sectional Study

Introduction

Studies have shown that genetic variation and environmental factors are associated with individual differences in therapeutic efficacy and side effects of opioids. However, the focus of these studies has been on a single factor of single-nucleotide polymorphisms (SNPs) or haplotypes, for which results have rarely been validated. For complex traits, such as cancer pain and opioid response, interactions between multiple genetic variation and environmental factors need to be considered to explain the opioid individual differences.

Methods

We conducted an exploratory two-stage cross-sectional study with 1027 Chinese patients who were taking strong opioid medications for their cancer pain, and genotyped 110 SNPs to explore the association of SNPs, haplotypes, gene–gene and gene–environment interactions with opioid dose, pain relief, and opioid-induced constipation.

Results

Due to the failure to meet Benjamini–Hochberg criteria in the discovery stage or to be validated in replication stage, no association was found between SNPs, haplotypes, paired SNP–SNP interactions or multi-dimensional gene–gene interactions and opioid response. However, for gene–environment interactions, optimal models have been constructed in all phenotypes of opioid response.

Conclusions

This study reveals for the first time that construction of multidimensional gene–environment interactions enables better interpretations of the effect of genetic variation and environmental factors on the opioid response in patients with cancer pain.

Trial registration

Chictr.org.cn, identifier, ChiCTR2000033576.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40122-022-00353-5.

MMP9 SNP and MMP SNP-SNP interactions increase the risk for ischemic stroke in the Han Hakka population

OBJECTIVES

To investigate the association of eight variants of four matrix metalloproteinase (MMP) genes with ischemic stroke (IS) and whether interactions among these single nucleotide polymorphisms (SNPs) increases the risk of IS.

METHODS

Among 547 patients with ischemic stroke and 350 controls, matrix-assisted laser desorption/ionization time of flight mass spectrometry was used to examine eight variants arising from four different genes, including MMP-1 (rs1799750), MMP-2 (rs243865, rs2285053, rs2241145), MMP-9 (rs17576), and MMP-12 (rs660599, rs2276109, and rs652438). Gene-gene interactions were employed using generalized multifactor dimensionality reduction (GMDR) methods.

RESULTS

The frequency of rs17576 was significantly higher in IS patients than in controls (p = .033). Logistic regression analysis revealed the AG and GG genotypes of rs17576 to be associated with a higher risk for IS, with the odds ratio and 95% confidence interval being 2.490 (1.251-4.959) and 2.494 (1.274-4.886), respectively. GMDR analysis showed a significant SNP-SNP interaction between rs17576 and rs660599 (the testing balanced accuracy was 53.70% and cross-validation consistency was 8/10, p = .0107). Logistic regression analysis showed the interaction between rs17576 and rs660599 to be an independent risk factor for IS with an odds ratio of 1.568 and a 95% confidence interval of 1.152-2.135.

CONCLUSION

An MMP-9 rs17576 polymorphism is associated with increased IS risk in the Han Hakka population and interaction between MMP-9 rs17576 and MMP-12 rs660599 is associated with increased IS risk as well.

The relationship between polymorphism of insulin-like growth factor I gene and susceptibility to type 2 diabetes in Uygur population, Xinjiang, China

Background

Type 2 diabetes (T2DM) susceptibility varies among different populations and is affected by gene single nucleotide polymorphism (SNP). Insulin-like growth factor (IGF)-1 gene, which has many SNP loci, is involved in T2DM pathogenesis. However, the relationship of IGF-1 gene polymorphism with T2DM in Uyghur population is less studied.

Objective

To investigate the relationship between T2DM susceptibility and polymorphism of IGF-1 gene in Uyghur population of Xinjiang, China.

Methods

This study enrolled 220 cases (122 males (55.46%) and 98 females (44.54%); mean age of 53.40 ± 10.94 years) of T2DM patients (T2DM group) and 229 (124 males (54.15%) and 105 females (45.85%); mean age of 51.64 ± 10.48 years) healthy controls (control group). Biochemical indexes were determined. IGF-1 gene polymorphism was analyzed by SNP genotyping.

Results

The levels of TG, HDL, LDL, BUN, and Cr were statistically significant between the T2DM group and the control group. In terms of IGF-1 polymorphism, T2DM group had higher frequency of AA genotype (OR = 2.40, 95% CI = 1.19–4.84) and allele A (OR = 1.55, 95% CI = 1.17–2.06) of rs35767 loci, suggesting that rs35767 is related to the occurrence of T2DM. A total of 5 gene interaction models was obtained through analyzing the interaction of 5 SNP loci with the GMDR method. Among them, the two-factor model that included rs35767 locus and rs5742694 locus had statistical difference with a large cross-validation consistency (10/10). The combination of GG/CC, GA/AA, AA/AA, and AA/AC genotype was in high-risk group, whereas the combination of GG/AA, GG/AC, GA/AC and GA/CC genotype was in the low-risk group. The risk of T2DM in the high-risk group was 2.165 times than that of the low-risk group (OR = 2.165, 95% CI = 1.478–3.171).

Conclusion

TG, HDL, LDL, BUN, and Cr are influencing factors of T2DM in Uyghur population. The rs35767 locus of IGF-1 gene may be associated with T2DM in Uyghur population. The high-risk group composing of rs35767 locus and rs5742694 locus has a higher risk of T2DM.

Uncovering complementary sets of variants for predicting quantitative phenotypes

MOTIVATION

Genome-wide association studies show that variants in individual genomic loci alone are not sufficient to explain the heritability of complex, quantitative phenotypes. Many computational methods have been developed to address this issue by considering subsets of loci that can collectively predict the phenotype. This problem can be considered a challenging instance of feature selection in which the number of dimensions (loci that are screened) is much larger than the number of samples. While currently available methods can achieve decent phenotype prediction performance, they either do not scale to large datasets or have parameters that require extensive tuning.

RESULTS

We propose a fast and simple algorithm, Macarons, to select a small, complementary subset of variants by avoiding redundant pairs that are likely to be in linkage disequilibrium. Our method features two interpretable parameters that control the time/performance trade-off without requiring parameter tuning. In our computational experiments, we show that Macarons consistently achieves similar or better prediction performance than state-of-the-art selection methods while having a simpler premise and being at least two orders of magnitude faster. Overall, Macarons can seamlessly scale to the human genome with ∼107 variants in a matter of minutes while taking the dependencies between the variants into account.

AVAILABILITYAND IMPLEMENTATION

Macarons is available in Matlab and Python at https://github.com/serhan-yilmaz/macarons.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Analysis of rs1864182 and rs1864183 variants in ATG10 gene and antineutrophil cytoplasmic autoantibody-associated vasculitis in Chinese Guangxi population

Objectives

To investigate the association of autophagy‐associated gene 10 (ATG10) gene polymorphisms (rs1864182 and rs1864183) with antineutrophil cytoplasmic autoantibody (ANCA)‐associated vasculitis (AAV) in Chinese Guangxi population.

Methods

The single nucleotide polymorphisms (SNPs) of ATG10 rs1864182 and rs1864183 in 395 participants (195 AAVs and 200 healthy controls) were genotyped. Generalized multiple dimensionality reduction (GMDR) was used to analyze the SNP‐SNP interactions among two SNPs of ATG10 gene and other SNPs of autophagy gene previously studied by our research team.

Results

In this study, we found that the two ATG10 SNPs were not associated with AAV risk in Chinese Guangxi population. However, there were statistically significant differences in the incidence of hemoptysis, hematuria, and proteinuria among the three genotypes of ATG10 rs1864182 and rs1864183 (p < 0.05). Moreover, permutation test of GMDR suggested that immunity‐related GTPase M(IRGM) rs4958847, autophagy‐associated gene 7 (ATG7) rs6442260, ATG7 rs2594966, ATG10 rs1864183, protein kinase B(AKT2) rs3730051, and AKT2 rs11552192 might interact with each other in the process of developing AAV (p < 0.05).

Conclusions

Our results indicated that there existed no association between ATG10 SNPs and AAV, and SNP‐SNP interactions among IRGM rs4958847, ATG7 rs6442260, ATG7 rs2594966, ATG10 rs1864183, AKT2 rs3730051, and AKT2 rs11552192 may confer AAV risk in the Chinese Guangxi population.

Effects of SNPs in SOD2 and SOD3 interacted with fluoride exposure on the susceptibility of dental fluorosis

A total of 649 children aged 7-13 years of age were recruited in a cross-sectional study in Tongxu County, China (2017) to assess the effects of interaction between single nucleotide polymorphisms (SNPs) in SOD2 and SOD3 gene and fluoride exposure on dental fluorosis (DF) status. Associations between biomarkers and DF status were evaluated. Logistic regression suggested that the risk of DF in children with rs10370 GG genotype and rs5746136 TT genotype was 1.89-fold and 1.72-fold than that in children with TT/CC genotype, respectively. Increased T-SOD activity was associated with a lower risk of DF (OR = 0.99). The rs2855262*rs10370*UF model was regarded as the optimal interaction model in generalized multifactor dimensionality reduction analyses. Our findings suggested that rs4880 and rs10370 might be useful genetic markers for DF, and there might be interactions among rs10370 in SOD2, rs2855262 in SOD3, and fluoride exposure on DF status.

CYP17A1-ATP2B1 SNPs and Gene-Gene and Gene-Environment Interactions on Essential Hypertension

Background: The association between the CYP17A1 and ATP2B1 SNPs and essential hypertension (referred to as hypertension) is far from being consistent. In addition to the heterogeneity of hypertension resulting in inconsistent results, gene–gene and gene–environment interactions may play a major role in the pathogenesis of hypertension rather than a single gene or environmental factor.

Methods: A case–control study consisting of 1,652 individuals (hypertension, 816; control, 836) was conducted in Maonan ethnic minority of China. Genotyping of the four SNPs was performed by the next-generation sequencing technology.

Results: The frequencies of minor alleles and genotypes of four SNPs were different between the two groups (p < 0.001). According to genetic dominance model analysis, three (rs1004467, rs11191548, and rs17249754) SNPs and two haplotypes (CYP17A1 rs1004467G-rs11191548C and ATP2B1 rs1401982G-rs17249754A) were negatively correlated, whereas rs1401982 SNP and the other two haplotypes (CYP17A1 rs1004467A-rs11191548T and ATP2B1 rs1401982A-rs17249754G) were positively associated with hypertension risk (p ≤ 0.002 for all). Two best significant two-locus models were screened out by GMDR software involving SNP–environment (rs11191548 and BMI ≥ 24 kg/m²) and haplotype–environment (CYP17A1 rs1004467G-rs11191548C and BMI ≥ 24 kg/m²) interactions (p ≤ 0.01). The subjects carrying some genotypes increased the hypertension risk.

Conclusions: Our outcomes implied that the rs1004467, rs11191548, and rs17249754 SNPs and CYP17A1 rs1004467G-rs11191548C and ATP2B1 rs1401982G-rs17249754A haplotypes have protective effects, whereas the rs1401982 SNP and CYP17A1 rs1004467A-rs11191548T and ATP2B1 rs1401982A-rs17249754G haplotypes showed adverse effect on the prevalence of hypertension. Several SNP–environment interactions were also detected.

Spatial rank-based multifactor dimensionality reduction to detect gene-gene interactions for multivariate phenotypes

Background

Identifying interaction effects between genes is one of the main tasks of genome-wide association studies aiming to shed light on the biological mechanisms underlying complex diseases. Multifactor dimensionality reduction (MDR) is a popular approach for detecting gene–gene interactions that has been extended in various forms to handle binary and continuous phenotypes. However, only few multivariate MDR methods are available for multiple related phenotypes. Current approaches use Hotelling’s T² statistic to evaluate interaction models, but it is well known that Hotelling’s T² statistic is highly sensitive to heavily skewed distributions and outliers.

Results

We propose a robust approach based on nonparametric statistics such as spatial signs and ranks. The new multivariate rank-based MDR (MR-MDR) is mainly suitable for analyzing multiple continuous phenotypes and is less sensitive to skewed distributions and outliers. MR-MDR utilizes fuzzy k-means clustering and classifies multi-locus genotypes into two groups. Then, MR-MDR calculates a spatial rank-sum statistic as an evaluation measure and selects the best interaction model with the largest statistic. Our novel idea lies in adopting nonparametric statistics as an evaluation measure for robust inference. We adopt tenfold cross-validation to avoid overfitting. Intensive simulation studies were conducted to compare the performance of MR-MDR with current methods. Application of MR-MDR to a real dataset from a Korean genome-wide association study demonstrated that it successfully identified genetic interactions associated with four phenotypes related to kidney function. The R code for conducting MR-MDR is available at https://github.com/statpark/MR-MDR.

Conclusions

Intensive simulation studies comparing MR-MDR with several current methods showed that the performance of MR-MDR was outstanding for skewed distributions. Additionally, for symmetric distributions, MR-MDR showed comparable power. Therefore, we conclude that MR-MDR is a useful multivariate non-parametric approach that can be used regardless of the phenotype distribution, the correlations between phenotypes, and sample size.

The Interaction of Single Nucleotide Polymorphisms on Fibroblast Growth Factor 19 Superfamily Genes Is Associated With Alcohol Dependence-Related Aggression

Alcohol dependence (AD) is characterized by compulsive alcohol consumption, which involves behavioral impairments such as aggression. Members of fibroblast growth factor (FGF) 19 superfamily, including FGF19, FGF21, and FGF23, are major endocrine mediators that play an important role in alcohol metabolism and alcohol related disorders. The objective of the present study is to explore the possible associations among the interaction of single nucleotide polymorphisms (SNPs) of the FGF 19 superfamily, AD occurrence, and aggression in patients with AD. A total of 956 subjects were enrolled in this study, including 482 AD patients and 474 healthy controls (HCs). Michigan alcoholism screening test (MAST) was used to measure the level of AD, a Chinese version of the Buss–Perry Aggression Questionnaire was used to evaluate the aggressive behavior of subjects, and MassARRAY^@ system was used to genotype rs948992 of FGF19, rs11665841 and rs11665896 of FGF21, rs7955866 and rs11063118 of FGF23. The results showed that AD patients presented a significantly higher level of aggression compared to HCs, and MAST scores were significantly positively associated Buss–Perry aggression scores (r = 0.402, p < 0.001) in AD patients. The interaction of FGF19 rs948992 TC × FGF21 rs11665896 GG presented the high-risk genotype combination predicting the high level of AD. In addition, the interaction of FGF19 rs948992 TC × FGF21 rs11665896 TG × FGF23 rs11063118 TT presented the high-risk genotype combination predicting the high level of aggression in AD patients. Our results added evidence linking the combination of rs948992 TC × rs11665896 TG × rs11063118 TT to aggressive behavior in AD patients and pointed out the potential usefulness of the SNPs of FGF19 superfamily as a predictor for the aggression in AD patients.

An association study in the Taiwan Biobank elicits three novel candidates for cognitive aging in old adults: NCAM1, TTC12 and ZBTB20

The dopamine receptor-related loci have been suggested to be associated with cognitive functions and neurodegenerative diseases. It is unknown whether genetic variants such as single nucleotide polymorphisms (SNPs) in the dopamine receptor-related loci could contribute to cognitive aging independently as well as by virtue of complicated interplays in the elder population. To assess whether SNPs in the dopamine receptor-related loci are associated with cognitive aging in the elder population, we evaluated SNPs in the DRD1, NCAM1-TTC12-ANKK1-DRD2, DRD3-LOC107986115-ZNF80-TIGIT-MIR568-ZBTB20, DRD4, and DRD5-SLC2A9 loci from 25,195 older Taiwanese individuals from the Taiwan Biobank. Mini-Mental State Examination (MMSE) was scrutinized for all participants, where MMSE scores were employed to evaluate cognitive functions. From our analysis, we identified three novel genes for cognitive aging that have not previously been reported: ZBTB20 on chromosome 3 and NCAM1 and TTC12 on chromosome 11. NCAM1 and ZBTB20 are strong candidates for having a role in cognitive aging with mutations in ZBTB20 resulting in intellectual disability, and NCAM1 previously found to be associated with associative memory in humans. Additionally, we found the effects of interplays between physical activity and these three novel genes. Our study suggests that genetic variants in the dopamine receptor-related loci may influence cognitive aging individually and by means of gene-physical activity interactions.

Hyppo-X: A Scalable Exploratory Framework for Analyzing Complex Phenomics Data

Phenomics is an emerging branch of modern biology that uses high throughput phenotyping tools to capture multiple environmental and phenotypic traits, often at massive spatial and temporal scales. The resulting high dimensional data represent a treasure trove of information for providing an in-depth understanding of how multiple factors interact and contribute to the overall growth and behavior of different genotypes. However, computational tools that can parse through such complex data and aid in extracting plausible hypotheses are currently lacking. In this article, we present Hyppo-X, a new algorithmic approach to visually explore complex phenomics data and in the process characterize the role of environment on phenotypic traits. We model the problem as one of unsupervised structure discovery, and use emerging principles from algebraic topology and graph theory for discovering higher-order structures of complex phenomics data. We present an open source software which has interactive visualization capabilities to facilitate data navigation and hypothesis formulation. We test and evaluate Hyppo-X on two real-world plant (maize) data sets. Our results demonstrate the ability of our approach to delineate divergent subpopulation-level behavior. Notably, our approach shows how environmental factors could influence phenotypic behavior, and how that effect varies across different genotypes and different time scales. To the best of our knowledge, this effort provides one of the first approaches to systematically formalize the problem of hypothesis extraction for phenomics data. Considering the infancy of the phenomics field, tools that help users explore complex data and extract plausible hypotheses in a data-guided manner will be critical to future advancements in the use of such data.

Association between TFEB gene polymorphism, gene-environment interaction, and fatty liver disease: A case-control study in China

BACKGROUND

Fatty liver disease (FLD) is a serious public health problem that is rapidly increasing. Evidence indicates that the transcription factor EB (TFEB) gene may be involved in the pathophysiology of FLD; however, whether TEFB polymorphism has an association with FLD remains unclear.

OBJECTIVES

To explore the association among TFEB polymorphism, gene-environment interaction, and FLD and provide epidemiological evidence for clarifying the genetic factors of FLD.

METHODS

This study is a case-control study. Sequenom MassARRAY was applied in genotyping. Logical regression was used to analyze the association between TFEB polymorphism and FLD, and the gene-environment interaction in FLD was evaluated by multiplication and additive interaction models.

RESULTS

(1) The alleles and genotypes of each single nucleotide polymorphism and haplotypes of TFEB in the case and control groups were evenly distributed; no statistically substantial difference was observed. (2) Logistic regression analysis indicated that TFEB polymorphism is not remarkably associated with FLD. (3) In the multiplicative interaction model, rs1015149, rs1062966, rs11754668 and rs2273068 had remarkable interaction with the amount of cigarette smoking. Rs1062966 and rs11754668 also had a considerable interaction body mass index and alcohol intake, respectively. However, no remarkable additive interaction was observed.

CONCLUSION

TFEB polymorphism is not directly associated with FLD susceptibility, but the risk can be changed through gene-environment interaction.

Transcription factor 7-like 2 gene- smoking interaction on the risk of diabetic nephropathy in Chinese Han population

Objectives

To evaluate the relationship between transcription factor 7-like 2 (TCF7L2) gene polymorphism and diabetic nephropathy (DN) risk, as well as the effect of gene-environment interactions on DN risk in Chinese Han population.

Methods

The Hardy-Weinberg equilibrium (HWE) and the relationship between TCF7L2 gene single nucleotide polymorphism (SNPs) and DN susceptibility were evaluated by SNPStats. The interaction among four SNPs and environmental factors were tested by generalized multifactor dimensionality reduction (GMDR). The consistency of cross validation, accuracy of test balance and sign test were calculated to evaluate the interaction of each selection. The logistic regression was used to test the interaction between rs7903146 and current smoking by stratified analysis.

Results

Logistic regression analysis indicated that the DN risk of rs7903146-T allele carriers were obviously higher than that in CC genotype carriers (CT + TT versus CC), adjusted OR (95 %CI) = 1.64 (1.24–2.06). However, we also discovered that people with rs12255372, rs11196205 and rs290487 minor allele had non-significant difference risk of DN compared with people with major allele. The GMDR model found a significant two-locus model (p = 0.0100) including rs7903146 and current smoking, suggesting a potential gene–environment interaction between rs7903146 and current smoking. Compared with never smokers with rs7903146- CC genotype, current smokers with rs7903146- CT or TT genotype had the highest DN risk. After covariate adjustment, OR (95 %CI) was 2.15 (1.58–2.78).

Conclusions

We found a significant relationship of rs7903146-T alleles, and the interaction between rs7903146-T and current smoking with increased DN risk.

MMP-2 and MMP-9 Polymorphisms and Preeclampsia Risk in Tunisian Arabs: A Case-Control Study

The abnormal production of matrix metalloproteinases (MMPs), especially MMP-9 and MMP-2, plays a pivotal role in hypertensive disorders of pregnancy, and as such, can influence the development of preeclampsia. These alterations may result from functional genetic polymorphisms in the promoter region of MMP-9 and MMP-2 genes, which modify MMP-9 and MMP-2 expression. We investigated the association of MMP-9 polymorphism rs3918242 (-1562 C>T) and MMP-2 polymorphism rs2285053 (-735 C>T) with the risk of preeclampsia. This case–control study was conducted on 345 women with preeclampsia and 281 age-matched women with normal pregnancies from Tunisian hospitals. Genomic DNA was extracted from whole blood collected at delivery. Genotypes for -1562 C>T and -735 C>T polymorphisms were performed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). An increased frequency of heterozygous MMP-9 -1562 C/T genotype carriers was observed in women with preeclampsia compared to healthy controls (p = 0.03). In contrast, the MMP-2 -735 C>T polymorphism was not significantly different regarding frequency distribution of the allele and genotype between healthy pregnant women and women with preeclampsia. Our study suggests that the MMP-9 -1562 C/T variant, associated with high MMP-9 production, could be a genetic risk factor for preeclampsia in Tunisian women.

Germline Variants and Genetic Interactions of Several EMT Regulatory Genes Increase the Risk of HBV-Related Hepatocellular Carcinoma

Epithelial-mesenchymal transition (EMT) plays an important role in the development of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). We hypothesized that germline variants in the major EMT regulatory genes (SNAIL1, ZEB1, ZEB2, TWIST1) may influence the development of HBV-related HCC. We included 421 cases of HBsAg-positive patients with HCC, 1371 cases of HBsAg-positive subjects without HCC [patients with chronic hepatitis B (CHB) or liver cirrhosis (LC)] and 618 cases of healthy controls in the case-control study. Genotype, allele, and haplotype associations in the major EMT regulatory genes were tested. Environment-gene and gene-gene interactions were analysed using the non-parametric model-free multifactor dimensionality reduction (MDR) method. The SNAIL1rs4647958T>C was associated with a significantly increased risk of both HCC (CT+CC vs. TT: OR=1.559; 95% confidence interval [CI], 1.073-2.264; P=0.020) and CHB+LC (CT+CC vs. TT: OR=1.509; 95% CI, 1.145-1.988; P=0.003). Carriers of the TWIST1rs2285681G>C (genotypes CT+CC) had an increased risk of HCC (CG+CC vs. GG: OR=1.407; 95% CI, 1.065-1.858; P=0.016). The ZEB2rs3806475T>C was associated with significantly increased risk of both HCC (P _recessive =0.001) and CHB+LC (P _recessive<0.001). The CG haplotype of the rs4647958/rs1543442 haplotype block was associated with significant differences between healthy subjects and HCC patients (P=0.0347). Meanwhile, the CT haplotype of the rs2285681/rs2285682 haplotype block was associated with significant differences between CHB+LC and HCC patients (P=0.0123). In MDR analysis, the combination of TWIST1rs2285681, ZEB2rs3806475, SNAIL1rs4647958 exhibited the most significant association with CHB+LC and Health control in the three-locus model. Our results suggest significant single-gene associations and environment-gene/gene-gene interactions of EMT-related genes with HBV-related HCC.

Association of Extracellular Signal-Regulated Kinase Genes With Myopia: A Longitudinal Study of Chinese Children

Objective

The present study was designed to investigate whether the extracellular signal-regulated kinase (ERK) signaling pathway, a downstream component of dopamine signaling, is involved in myopia among Chinese children.

Methods

During a 3.5-year follow-up, 488 primary school students were enrolled in this study. Non-cycloplegic spherical equivalent refraction (SE) and other ocular parameters were assessed. Four variants of four genes in the ERK signaling pathway were selected: RASGRF1 rs6495367, PTPN5 rs1550870, PTPRR rs11178469, and PDGFRA rs6554163. SNPscan was used to genotype single-nucleotide polymorphisms (SNPs). PLINK software was used to assess the associations of the genetic variants with the occurrence or development of myopia, SE, and other ocular parameters. We created a protein-protein interaction (PPI) network and microRNA (miRNA)-gene network using String and Cytoscape and conducted enrichment analyses on the genes in these networks.

Results

In total, 426 children (baseline age: 7.28 ± 0.26 years; 236 (55.4%) boys and 190 girls) wereenrolled. After adjusting for confounding factors with 10,000 permutations, children with the CT or TT genotype of PTPN5 rs1550870 were more susceptible to myopia than those with the CC genotype (adjusted p = 0.011). Additionally, PTPN5 rs1550870 was correlated with significant myopic shift and increasing axial length (AL) and lens thickness (LT) but had a negative effect on central corneal thickness (CCT). RASGRF1 rs6495367 was negatively associated with myopic shift (additive: adjusted p = 0.034; dominant: adjusted p = 0.020), myopic SE and AL. PDGFRA rs6554163 TA or AA was negatively associated with increasing LT (adjusted p = 0.033). Evaluation of the effects of SNP-SNP combinations on incident myopia revealed a statistically significant one-locus model: PTPN5 rs1550870 [cross-validation consistency (CVC) = 10/10, adjusted p = 0.0107]. The genes in the PPI and miRNA-gene interaction networks were subjected to enrichment analyses, which suggested that these genes are involved mainly in eye development and dopaminergic synapse-related processes.

Conclusion

We identified genetic variants of crucial ERK signaling pathway genes that were significantly correlated with myopia and ocular parameter alterations in Chinese children. A combination of gene and miRNA functional analyses with enrichment analyses highlights the regulatory effects associated with ocular development and dopamine biological functions. This study offers novel clues to understand the role of dopamine in the molecular mechanisms of myopia.

Interactions among variants in P53 apoptotic pathway genes are associated with neurologic deterioration and functional outcome after acute ischemic stroke

OBJECTIVE

Neurologic deterioration (ND) and functional outcome after ischemic stroke (IS) are not accurately predicted by clinical pictures on admission. The aim of present study was to investigate the association of variants in P53 apoptotic pathway genes with ND and functional outcome after IS.

METHODS

Genotypes of nine variants in apoptosis-relevant genes were measured in patients with acute IS. Gene-gene interactions were analyzed by generalized multifactor dimensionality reduction (GMDR). The primary outcome was ND. ND was diagnosed in patients who worsened ≥2 points (National Institutes of Health Stroke Scale [NIHSS] score) within the first 10 days of stroke onset. The secondary outcome was functional status at 90 days after IS as measured by modified Rankin Scale (mRS) score.

RESULTS

A total of 705 enrolled patients, ND occurred in 174 (24.7%) patients, and 184 (26.1%) patients were poor functional outcome (mRS score > 2). Although the nine variants were not significantly associated with ND and functional outcome by univariate analysis, there was a gene-gene interaction among P53rs1042522, MDM-2rs2279744, and MMP-9 rs3918242 using GMDR analysis. The high-risk interaction among the three variants was independently associated with higher risk of ND (HR, 2.04, 95% CI: 1.22-5.64, p = .018) and poor functional outcome (OR, 2.68, 95% CI: 1.68-7.86, p = .004) after adjusting for the covariates.

CONCLUSION

The interactions among P53 rs1042522, MDM-2 rs2279744, and MMP-9 rs3918242 may increase the risk of ND and poor functional outcome and may be considered as a genetic marker of predicting ND and poor functional outcome after stroke.

SPADIS: An Algorithm for Selecting Predictive and Diverse SNPs in GWAS

Phenotypic heritability of complex traits and diseases is seldom explained by individual genetic variants identified in genome-wide association studies (GWAS). Many methods have been developed to select a subset of variant loci, which are associated with or predictive of the phenotype. Selecting connected SNPs on SNP-SNP networks have been proven successful in finding biologically interpretable and predictive SNPs. However, we argue that the connectedness constraint favors selecting redundant features that affect similar biological processes and therefore does not necessarily yield better predictive performance. In this paper, we propose a novel method called SPADIS that favors the selection of remotely located SNPs in order to account for their complementary effects in explaining a phenotype. SPADIS selects a diverse set of loci on a SNP-SNP network. This is achieved by maximizing a submodular set function with a greedy algorithm that ensures a constant factor approximation to the optimal solution. We compare SPADIS to the state-of-the-art method SConES, on a dataset of Arabidopsis Thaliana with continuous flowering time phenotypes. SPADIS has better average phenotype prediction performance in 15 out of 17 phenotypes when the same number of SNPs are selected and provides consistent improvements across multiple networks and settings on average. Moreover, it identifies more candidate genes and runs faster.

Endocannabinoid Gene × Gene Interaction Association to Alcohol Use Disorder in Two Adolescent Cohorts

Genetic markers of the endocannabinoid system have been linked to a variety of addiction-related behaviors that extend beyond cannabis use. In the current study we investigate the relationship between endocannabinoid (eCB) genetic markers and alcohol use disorder (AUD) in European adolescents (14-18 years old) followed in the IMAGEN study (n = 2,051) and explore replication in a cohort of North American adolescents from Canadian Saguenay Youth Study (SYS) (n = 772). Case-control status is represented by a score of more than 7 on the Alcohol Use Disorder Identification Test (AUDIT). First a set-based test method was used to examine if a relationship between the eCB system and AUDIT case/control status exists at the gene level. Using only SNPs that are both independent and significantly associated to case-control status, we perform Fisher's exact test to determine SNP level odds ratios in relation to case-control status and then perform logistic regressions as post-hoc analysis, while considering various covariates. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the most robust SNP×SNP interaction of the five eCB genes with positive AUDIT screen. While no gene-sets were significantly associated to AUDIT scores after correction for multiple tests, in the case/control analysis, 7 SNPs were significantly associated with AUDIT scores of > 7 (p < 0.05; OR<1). Two SNPs remain significant after correction by false discovery rate (FDR): rs9343525 in CNR1 (pcorrected =0.042, OR = 0.73) and rs507961 in MGLL (pcorrected = 0.043, OR = 0.78). Logistic regression showed that both rs9353525 (CNR1) and rs507961 (MGLL) remained significantly associated with positive AUDIT screens (p < 0.01; OR < 1) after correction for multiple covariables and interaction of covariable × SNP. This result was not replicated in the SYS cohort. The GMDR model revealed a significant three-SNP interaction (p = 0.006) involving rs484061 (MGLL), rs4963307 (DAGLA), and rs7766029 (CNR1) predicted case-control status, after correcting for multiple covariables in the IMAGEN sample. A binomial logistic regression of the combination of these three SNPs by phenotype in the SYS cohort showed a result in the same direction as seen in the IMAGEN cohort (BETA = 0.501, p = 0.06). While preliminary, the present study suggests that the eCB system may play a role in the development of AUD in adolescents.

Candidate Genes for Age at Menarche Are Associated With Uterine Leiomyoma

Age at menarche (AAM) is an important marker of the pubertal development and function of the hypothalamic-pituitary-ovarian system. It was reported as a possible factor for a risk of uterine leiomyoma (UL). However, while more than 350 loci for AAM have been determined by genome-wide association studies (GWASs) to date, no studies of these loci for their association with UL have been conducted so far. In this study, we analyzed 52 candidate loci for AAM for possible association with UL in a sample of 569 patients and 981 controls. The results of the study suggested that 23 out of the 52 studied polymorphisms had association with UL. Locus rs7759938 LIN28B was individually associated with the disease according to the dominant model. Twenty loci were associated with UL within 11 most significant models of intergenic interactions. Nine loci involved in 16 most significant models of interactions between single-nucleotide polymorphism (SNP), induced abortions, and chronic endometritis were associated with UL. Among the 23 loci associated with UL, 16 manifested association also with either AAM (7 SNPs) or height and/or body mass index (BMI) (13 SNPs). The above 23 SNPs and 514 SNPs linked to them have non-synonymous, regulatory, and expression quantitative trait locus (eQTL) significance for 35 genes, which play roles in the pathways related to development of the female reproductive organs and hormone-mediated signaling [false discovery rate (FDR) ≤ 0.05]. This is the first study reporting associations of candidate genes for AAM with UL.

Pro-inflammatory cytokine polymorphisms and interactions with dietary alcohol and estrogen, risk factors for invasive breast cancer using a post genome-wide analysis for gene-gene and gene-lifestyle interaction

Molecular and genetic immune-related pathways connected to breast cancer and lifestyles in postmenopausal women are not fully characterized. In this study, we explored the role of pro-inflammatory cytokines such as C-reactive protein (CRP) and interleukin-6 (IL-6) in those pathways at the genome-wide level. With single-nucleotide polymorphisms (SNPs) in the biomarkers and lifestyles together, we further constructed risk profiles to improve predictability for breast cancer. Our earlier genome-wide association gene-environment interaction study used large cohort data from the Women's Health Initiative Database for Genotypes and Phenotypes Study and identified 88 SNPs associated with CRP and IL-6. For this study, we added an additional 68 SNPs from previous GWA studies, and together with 48 selected lifestyles, evaluated for the association with breast cancer risk via a 2-stage multimodal random survival forest and generalized multifactor dimensionality reduction methods. Overall and in obesity strata (by body mass index, waist, waist-to-hip ratio, exercise, and dietary fat intake), we identified the most predictive genetic and lifestyle variables. Two SNPs (SALL1 rs10521222 and HLA-DQA1 rs9271608) and lifestyles, including alcohol intake, lifetime cumulative exposure to estrogen, and overall and visceral obesity, are the most common and strongest predictive markers for breast cancer across the analyses. The risk profile that combined those variables presented their synergistic effect on the increased breast cancer risk in a gene-lifestyle dose-dependent manner. Our study may contribute to improved predictability for breast cancer and suggest potential interventions for the women with the risk genotypes and lifestyles to reduce their breast cancer risk.