Significant Gene-gene Interaction Research Articles

The impact of genetic variations in the serotonergic system on symptom severity and clinical outcome in functional neurological disorders

ObjectiveWe studied gene-environment, as well as gene-gene interaction to elucidate their effects on symptom severity and predict clinical outcomes in functional neurological disorders (FND). MethodsEighty-five patients with mixed FND were genotyped for ten single-nucleotide polymorphisms (SNP) from seven different stress-related genes. We tested cross-sectionally the association between genotype and the symptomatology of FND (symptom severity assessed with the examiner-based clinical global impression score [CGI] and age of onset). Clinical outcome was assessed in 52 patients who participated in a follow-up clinical visit after eight months (following their individual therapies as usual). We tested longitudinally the association between genotype and clinical outcome in FND. We examined the contribution of each SNP and their interaction between them to FND symptomatology and outcome. ResultsWe identified a nominal association between tryptophan hydroxylase 1 (TPH1) rs1800532 and symptom severity (CGI1) in FND under a codominant model (T/T: ßT/T = 2.31, seT/T = 0.57; G/T: ßG/T = -0.18, seG/T = 0.29, P = 0.035), with minor allele (T) carriers presenting more severe symptoms. An association was identified between TPH1 and clinical outcome, suggesting that major allele (G) carriers were more likely to have an improved outcome under a codominant model (G/T: ORG/T = 0.18, CIG/T = [0.02–1.34]; T/T: ORT/T = 2.08, CIT/T = [0.30–14.53], P = 0.041). Our analyses suggested a significant gene-gene interaction for TPH2 (rs4570625) and OXTR (rs2254298) on symptom severity, and a significant gene-gene interaction for TPH1, TPH2 and BDNF (rs1491850) on clinical outcome. ConclusionFND might arise from a complex interplay between individual predisposing risk genes involved in the serotonergic pathway and their gene-gene interactions.

GWAS-Identified Loci are Associated with Obesity and Type 2 Diabetes Mellitus in Patients with Severe COVID-19.

Comorbidities such as obesity and type 2 diabetes mellitus (T2DM) have emerged as critical risk factors exacerbating the severity and mortality of COVID-19. Meanwhile, numerous genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with increased susceptibility to severe COVID-19. This study investigated whether SNPs previously identified by GWAS as risk factors for severe COVID-19 also correlate with common comorbidities-obesity and T2DM-in hospitalized patients with severe COVID-19. DNA samples from 199 hospitalized COVID-19 patients were genotyped using probe-based PCR for 10 GWAS SNPs previously implicated in severe COVID-19 outcomes (rs143334143 CCHCR1, rs111837807 CCHCR1, rs17078346 SLC6A20-LZTFL1, rs17713054 SLC6A20-LZTFL1, rs7949972 ELF5, rs61882275 ELF5, rs12585036 ATP11A, rs67579710 THBS3, THBS3-AS1, rs12610495 DPP9, rs9636867 IFNAR2). The analysis revealed significant associations between certain SNPs and the increased risk of obesity and T2DM in severe COVID-19 patients. Specifically, rs17713054 SLC6A20-LZTFL1 (risk allele A; odds ratio (OR) = 2.34, 95% confidence interval (CI) = 1.24-4.4, p = 0.007) and rs7949972 ELF5 SNP (risk allele T; OR = 1.79, 95% CI = 1.11-2.91, p = 0.015) were associated with increased risk of obesity. SNP rs9636867 IFNAR2 was associated with a higher risk of T2DM (risk allele G, OR = 8.28, 95% CI = 1.69-40.64, p = 0.027). Using the model-based multifactor dimensionality reduction (MB-MDR) approach, the six most significant gene-gene interaction patterns associated with obesity in severe COVID-19 patients were identified and included five polymorphic loci: rs7949972, rs17713054, rs61882275, rs12585036, and rs143334143, participating in two or more of the most significant G-G interactions (pperm < 0.05). In total, the best models of G-G interactions associated with T2DM in patients with severe COVID-19 included eight polymorphic loci, six of which, rs7949972, rs61882275, rs12585036, rs143334143, rs67579710, and rs12610495, were involved in two or more of the most significant G-G interactions. Our study provides novel insights into the genetic associations between GWAS-identified SNPs and the risk of obesity and T2DM in patients with severe COVID-19.

Haplotype analysis on association between C-reactive protein gene and susceptibility to type 2 diabetes mellitus in Chinese Han population : CRP gene and type 2 diabetes mellitus.

We aimed to evaluate the impact of C-reactive protein (CRP) gene polymorphism, additional gene-gene interaction, and haplotypes on susceptibility to type 2 diabetes mellitus (T2DM). SNPstats online software ( https://www.snpstats.net/start.htm ) was employed to evaluate the association between CRP gene and T2DM risk. High-order interactions among SNPs was tested using generalized multifactor dimensionality reduction, and the testing balanced accuracy, training balanced accuracy and cross-validation consistency were calculated. The SHEsisPlus ( http://shesisplus.bio-x.cn/SHEsis.html ) online software was used for haplotype analysis. A total of 730 T2DM patients and 765 controls were enrolled. The T allele of rs1205 is associated with increased susceptibility to T2DM, OR (95% CI) were 1.51 (1.13-2.01), 1.44 (1.10-1.89) and 1.25 (1.01-1.54) for codominant, dominant and over-dominant models, respectively. We also found that minor allele of rs2794521 is associated with decreased susceptibility to T2DM under codominant and recessive models, OR (95%CI) were 0.38 (0.18-0.79) and 0.37 (0.16-0.80) for codominant and recessive models, respectively. No significant gene-gene interaction existed among CRP gene SNPs, allinteraction p- values were more than 0.05. Haplotype analyses suggested the CGCA haplotype containing rs1205-C, rs1130864-G, rs2794521- C and rs3093059- A allele was associated with decreased risk of T2DM, OR (95% CI) = 0.83 (0.68-0.98), P = 0.047. Minor allele of rs1205 was associated with increased T2DM risk. Minor allele of rs2794521 and the CGCA haplotype were associated with decreased T2DM risk.

Inflammation and endothelial function relevant genetic polymorphisms, carotid atherosclerosis, and vascular events in high-risk stroke population.

To identify the associations of 19 single nucleotide polymorphisms (SNPs) in genes involved in inflammation and endothelial function and carotid atherosclerosis with subsequent ischemic stroke and other vascular events in the high-risk stroke population. This was a multicenter community-based sectional survey and prospective cohort study in Sichuan, southwestern China. Eight communities were randomly selected, and the residents in each community were surveyed using a structured face-to-face questionnaire. Carotid ultrasonography and DNA information were obtained from 2,377 out of 2,893 individuals belonging to a high-risk stroke population. Genotypes of the 19 SNPs in genes involved in inflammation and endothelial function were measured. All the 2,377 subjects were followed up for 4.7 years after the face-to-face survey. The primary outcome was ischemic stroke, and the secondary outcome was a composite of vascular events. Among the 2,377 subjects, 2,205 (92.8%) completed a 4.7-year follow-up, 947 (42.9%) had carotid atherosclerosis [372 (16.9%) carotid vulnerable plaque, 405 (18.4%) mean IMT > 0.9 mm, 285 (12.0%) carotid stenosis ≥15%]. Outcomes occurred in 158 (7.2%) subjects [92 (4.2%) ischemic stroke, 17 (0.8%) hemorrhagic stroke, 48 (2.2%) myocardial infarction, and 26 (1.2%) death] during follow-up. There was a significant gene-gene interaction among ITGA2 rs1991013, IL1A rs1609682, and HABP2 rs7923349 in the 19 SNPs. The multivariate logistic regression model revealed that carotid atherosclerosis and the high-risk interactive genotypes among the three SNPs were independent with a higher risk for ischemic stroke (OR = 2.67, 95% CI: 1.52-6.78, p = 0.004; and OR = 3.11, 95% CI: 2.12-9.27, p < 0.001, respectively) and composite vascular events (OR = 3.04, 95% CI: 1.46-6.35, p < 0.001; and OR = 3.23, 95% CI: 1.97-8.52, p < 0.001, respectively). The prevalence of carotid atherosclerosis was shown to be very high in the high-risk stroke population. Specific SNPs, interactions among them, and carotid atherosclerosis were independently associated with a higher risk of ischemic stroke and other vascular events.

Association of maternal phthalates exposure and metabolic gene polymorphisms with congenital heart diseases: a multicenter case-control study

BackgroundThe majority of congenital heart diseases (CHDs) are thought to result from the interactions of genetics and the environment factors. This study aimed to assess the association of maternal non-occupational phthalates exposure, metabolic gene polymorphisms and their interactions with risk of CHDs in offspring.MethodsA multicenter case-control study of 245 mothers with CHDs infants and 268 control mothers of health infant was conducted from six hospitals. Maternal urinary concentrations of eight phthalate metabolites were measured by ultra-high performance liquid chromatography coupled with tandem mass spectrometry (UHPLC-MS/MS). Twenty single nucleotide polymorphisms (SNPs) in cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and 19 (CYP2C19), uridine diphosphate (UDP) glucuronosyl transferase family 1 member A7 (UGT1A7), family 2 member B7 (UGT2B7) and B15(UGT2B15) genes were genotyped. The multivariate logistic regressions were used to estimate the association between maternal phthalates exposure or gene polymorphisms and risk of CHDs. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the gene–gene and gene–phthalates exposure interactions.ResultsThere was no significant difference in phthalate metabolites concentrations between the cases and controls. No significant positive associations were observed between maternal exposure to phthalates and CHDs. The SNPs of UGT1A7 gene at rs4124874 (under three models, log-additive: aOR = 1.74, 95% CI:1.28–2.37; dominant: aOR = 1.86, 95% CI:1.25–2.78; recessive: aOR = 2.50, 95% CI: 1.26–4.94) and rs887829 (under the recessive model: aOR = 13.66, 95% CI: 1.54–121) were significantly associated with an increased risk of CHDs. Furthermore, the associations between rs4124874 (under log-additive and dominant models) of UGT1A7 were statistically significant after the false discovery rate correction. No significant gene-gene or gene-phthalate metabolites interactions were observed.ConclusionsThe polymorphisms of maternal UGT1A7 gene at rs4124874 and rs887829 were significantly associated with an increased risk of CHDs. More large-scale studies or prospective study designs are needed to confirm or refute our findings in the future.

Evolutionary conservation of putative suicidality-related risk genes that produce diminished motivation corrected by clozapine, lithium and antidepressants.

Genome wide association studies (GWAS) and candidate gene analyses have identified genetic variants and genes that may increase the risk for suicidal thoughts and behaviors (STBs). Important unresolved issues surround these tentative risk variants such as the characteristics of the associated genes and how they might elicit STBs. Putative suicidality-related risk genes (PSRGs) were identified by comprehensive literature search and were characterized with respect to evolutionary conservation, participation in gene interaction networks and associated phenotypes. Evolutionary conservation was established with database searches and BLASTP queries, whereas gene-gene interactions were ascertained with GeneMANIA. We then examined whether mutations in risk-gene counterparts in C. elegans produced a diminished motivation phenotype previously connected to suicide risk factors. From the analysis, 105 risk-gene candidates were identified and found to be: 1) highly conserved during evolution, 2) enriched for essential genes, 3) involved in significant gene-gene interactions, and 4) associated with psychiatric disorders, metabolic disturbances and asthma/allergy. Evaluation of 17 mutant strains with loss-of-function/deletion mutations in PSRG orthologs revealed that 11 mutants showed significant evidence of diminished motivation that manifested as immobility in a foraging assay. Immobility was corrected in some or all of the mutants with clozapine, lithium and tricyclic antidepressant drugs. In addition, 5-HT2 receptor and muscarinic receptor antagonists restored goal-directed behavior in most or all of the mutants. These studies increase confidence in the validity of the PSRGs and provide initial clues about possible mechanisms that mediate STBs.

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

To investigate the potential association between polymorphisms in genes involved in endothelial function, inflammation and carotid atherosclerosis. 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. 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). 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.

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. 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). 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. 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.

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

To evaluate the relationship between carotid stenosis with variants in genes referred to inflammation and endothelial function. 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). 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). 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.

Polymorphisms in drug metabolism genes predict the risk of refractory myasthenia gravis.

Nearly 10% to 20% of myasthenia gravis (MG) patients are refractory to conventional treatment for unclear reasons. The study aimed to explore the relationship between drug metabolism gene polymorphisms and refractory MG. One hundred and thirty-one MG patients (33 in the refractory group; 98 in the non-refractory group) admitted to Tongji Hospital were included in this retrospective study. Improved multiplex ligation detection reaction (iMLDR) was used to genotype 13 polymorphisms (NR3C1 rs17209237, rs9324921; FKBP5 rs1360780, rs4713904, rs9296158; HSP90AA1 rs10873531, rs2298877, rs7160651; MDR1 rs1045642, rs1128503, rs2032582; CYP3A4 rs2242480; and CYP3A5 rs776746). We applied multivariable logistic regression to investigate the association between refractory MG and nucleotide polymorphisms. Generalized multifactor dimensionality reduction (GMDR) was used to examine gene-gene interactions. CC genotype of HSP90AA1 rs7160651 was associated with the increased risk of refractory MG than CT genotype [odds ratio (OR) =0.26; P=0.041] and CT + TT genotype (dominant model, OR =0.24; P=0.022). For CYP3A5 rs776746, AA genotype was associated with refractory MG compared with AG genotype (OR =0.11; P=0.017), GG genotype (OR =0.18; P=0.033), and AG + GG genotype (dominant model, OR =0.16; P=0.020). The frequency of CAT haplotype of HSP90AA1 rs10873531, rs2298877, rs7160651 was less common in refractory patients (OR =0.33; P=0.044). No significant gene-gene interactions were observed. HSP90AA1 rs7160651 and CYP3A5 rs776746 were significantly associated with refractory MG. Further studies are warranted to confirm the results and investigate the use of polymorphisms for treatment individualization.

Association of APEX1 and XRCC1 Gene Polymorphisms With HIV-1 Infection Susceptibility and AIDS Progression in a Northern Chinese MSM Population.

Background: Some studies have shown that the base excision repair (BER) pathway has an effect on HIV-1 replication. APEX1 and XRCC1 as key BER genes may affect DNA repair capacity. However, the roles of single nucleotide polymorphisms (SNPs) in APEX1 and XRCC1 and their impact on HIV-1 infection and AIDS progression remain unclear. Methods: A custom-designed 48-Plex SNPscan Kit was used for detection of single nucleotide polymorphisms. 601 HIV-1-infected men who have sex with men (MSM) and 624 age-matched healthy individuals were recruited in northern China. Four SNPs (rs1130409, rs1760944, rs2307486 and rs3136817) in APEX1 gene and three SNPs (rs1001581, rs25487 and rs25489) in XRCC1 gene were genotyped. The generalized multifactor dimension reduction (GMDR) method was used to identify the SNP-SNP interactions. Results: In this study, rs1130409 G allele, rs1001581 C allele and rs25487 C allele were associated with a higher risk of HIV-1 infection susceptibility (p = 0.020, p = 0.007 and p = 0.032, respectively). The frequencies of APEX1 haplotype TT and XRCC1 haplotype CT showed significant differences between cases and controls (p = 0.0372 and p = 0.0189, respectively). Interestingly, stratified analysis showed that the frequency of rs1001581 C allele was significantly higher in AIDS patients with the CD4+ T-lymphocyte count <200 cells/μl than those with >200 cells/μl (p = 0.022). Moreover, significant gene-gene interactions among rs1130409, rs1001581 and rs25487 were identified by GMDR (p = 0.0107). Specially, individuals with five to six risk alleles have a higher susceptibility to HIV-1 infection than those with zero to two risk alleles (p < 0.001). Conclusion: APEX1 and XRCC1 gene polymorphisms were associated with the susceptibility to HIV-1 infection and AIDS progression in MSM populations in northern China.

Association and epistatic analysis of white matter related genes across the continuum schizophrenia and autism spectrum disorders: The joint effect of NRG1-ErbB genes

Background Schizophrenia-spectrum disorders (SSD) and Autism spectrum disorders (ASD) are neurodevelopmental disorders that share clinical, cognitive, and genetic characteristics, as well as particular white matter (WM) abnormalities. In this study, we aimed to investigate the role of a set of oligodendrocyte/myelin-related (OMR) genes and their epistatic effect on the risk for SSD and ASD. Methods We examined 108 SNPs in a set of 22 OMR genes in 1749 subjects divided into three independent samples (187 SSD trios, 915 SSD cases/control, and 91 ASD trios). Genetic association and gene-gene interaction analyses were conducted with PLINK and MB-MDR, and permutation procedures were implemented in both. Results Some OMR genes showed an association trend with SSD, while after correction, the ones that remained significantly associated were MBP, ERBB3, and AKT1. Significant gene-gene interactions were found between (i) NRG1*MBP (perm p-value = 0.002) in the SSD trios sample, (ii) ERBB3*AKT1 (perm p-value = 0.001) in the SSD case-control sample, and (iii) ERBB3*QKI (perm p-value = 0.0006) in the ASD trios sample. Discussion Our results suggest the implication of OMR genes in the risk for both SSD and ASD and highlight the role of NRG1 and ERBB genes. These findings are in line with the previous evidence and may suggest pathophysiological mechanisms related to NRG1/ERBBs signalling in these disorders.

Evaluation of the Genetic Association and Methylation of Immune Response Pathway Genes with the Risk of Chronic Periodontitis in the Uighur Population.

Aim: The aim of this study was to explore the possible associations between single nucleotide polymorphisms (SNPs) and DNA methylation levels of seven genes in the inflammatory response pathway with susceptibility to chronic periodontitis (CP) among the Uighur population of the Xinjiang Autonomous Region of China. Methods: A total of 444 eligible subjects (279 CP patients and 165 healthy controls) were enrolled in the study. Genomic DNA was obtained from gingival tissue for genotyping eight SNPs and performing methylation measurements of seven genes. Results: SNP rs2070745 in the formyl peptide receptor 1 (FPR1) gene achieved statistical significance in a standard allelic association analysis for CP (p = 0.02). The frequency of the rs2070745 minor allele G was higher in the cases than in controls (0.367 vs. 0.291). Additionally, rs2070745 was significantly associated with CP under the dominant genetic model (p = 0.03). Using logistic regression analysis, rs2070745 was found to be consistently associated with CP under the additive dominant model, and this association remained significant after covariates were taken into account [odds ratio (OR) = 1.49 (1.09-2.05), p = 0.014; OR = 1.58 (1.04-2.40), p = 0.031, respectively]. No significant gene-gene interactions were identified. Although we did not find a polymorphism in interleukin 6 (IL6) associated with CP in our study, the methylation level of a CpG island region located within the promoter region of IL6 was significantly less in CP patients compared with controls (p < 0.05). Conclusions: The genetic polymorphism rs2070745 in FPR1 and the methylation level of the promoter region of IL6 might be associated with CP in the Uighur population of China.

A comparative study on the unified model based multifactor dimensionality reduction methods for identifying gene-gene interactions associated with the survival phenotype

BackgroundFor gene-gene interaction analysis, the multifactor dimensionality reduction (MDR) method has been widely employed to reduce multi-levels of gene-gene interactions into high- or low-risk groups using a binary attribute. For the survival phenotype, the Cox-MDR method has been proposed using a martingale residual of a Cox model since Surv-MDR was first proposed using a log-rank test statistic. Recently, the KM-MDR method was proposed using the Kaplan-Meier median survival time as a classifier. All three methods used the cross-validation procedure to identify single nucleotide polymorphism (SNP) using SNP interactions among all possible SNP pairs. Furthermore, these methods require the permutation test to verify the significance of the selected SNP pairs. However, the unified model-based multifactor dimensionality reduction method (UM-MDR) overcomes this shortcoming of MDR by unifying the significance testing with the MDR algorithm within the framework of the regression model. Neither cross-validation nor permutation testing is required to identify SNP by SNP interactions in the UM-MDR method. The UM-MDR method comprises two steps: in the first step, multi-level genotypes are classified into high- or low-risk groups, and an indicator variable for the high-risk group is defined. In the second step, the significance of the indicator variable of the high-risk group is tested in the regression model included with other adjusting covariates. The Cox-UMMDR method was recently proposed by combining Cox-MDR with UM-MDR to identify gene-gene interactions associated with the survival phenotype. In this study, we propose two simple methods either by combining KM-MDR with UM-MDR, called KM-UMMDR or by modifying Cox-UMMDR by adjusting for the covariate effect in step 1, rather than in step 2, a process called Cox2-UMMDR. The KM-UMMDR method allows the covariate effect to be adjusted for in the regression model of step 2, although KM-MDR cannot adjust for the covariate effect in the classification procedure of step 1. In contrast, Cox2-UMMDR differs from Cox-UMMDR in the sense that the martingale residuals are obtained from a Cox model by adjusting for the covariate effect in step 1 of Cox2-UMMDR whereas Cox-UMMDR adjusts for the covariate effect in the regression model in step 2. We performed simulation studies to compare the power of several methods such as KM-UMMDR, Cox-UMMDR, Cox2-UMMDR, Cox-MDR, and KM-MDR by considering the effect of covariates and the marginal effect of SNPs. We also analyzed a real example of Korean leukemia patient data for illustration and a short discussion is provided.ResultsIn the simulation study, two different scenarios are considered: the first scenario compares the power of the cases with and without the covariate effect. The second scenario is to compare the power of cases with the main effect of SNPs versus without the main effect of SNPs. From the simulation results, Cox-UMMDR performs the best across all scenarios among KM-UMMDR, Cox2-UMMDR, Cox-MDR and KM-MDR. As expected, both Cox-UMMDR and Cox-MDR perform better than KM-UMMDR and KM-MDR when a covariate effect exists because the former adjusts for the covariate effect but the latter cannot. However, Cox2-UMMDR behaves similarly to KM-UMMDR and KM-MDR even though there is a covariate effect. This implies that the covariate effect would be more efficiently adjusted for in the regression model of the second step rather than under the classification procedure of the first step. When there is a main effect of any SNP, Cox-UMMDR, Cox2-UMMDR and KM-UMMDR perform better than Cox-MDR and KM-MDR if the main effects of SNPs are properly adjusted for in the regression model. From the simulation results of two different scenarios, Cox-UMMDR seems to be the most robust when there is either any covariate effect adjusting for or any SNP that has a main effect on the survival phenotype. In addition, the power of all methods decreased as the censoring fraction increased from 0.1 to 0.3, as heritability increased. The power of all methods seems to be greater under MAF = 0.2 than under MAF = 0.4. For illustration, both KM-UMMDR and Cox2-UMMDR were applied to identify SNP by SNP interactions with the survival phenotype to a real dataset of Korean leukemia patients.ConclusionBoth KM-UMMDR and Cox2-UMMDR were easily implemented by combining KM-MDR and Cox-MDR with UM-MDR, respectively, to detect significant gene-gene interactions associated with survival time without cross-validation and permutation testing. The simulation results demonstrate the utility of KM-UMMDR, Cox2-UMMDR and Cox-UMMDR, which outperforms Cox-MDR and KM-MDR when some SNPs with only marginal effects might mask the detection of causal epistasis. In addition, Cox-UMMDR, Cox2-UMMDR and Cox-MDR performed better than KM-UMMDR and KM-MDR when there were potentially confounding covariate effects.

Association of vitamin D receptor gene polymorphisms with gestational diabetes mellitus-a case control study in Wuhan, China

BackgroundGestational diabetes mellitus (GDM) increased risk of perinatal complications for both the women and the fetuses. The association between the vitamin D receptor (VDR) gene polymorphism and GDM has not been thoroughly investigated in Chinese pregnant women. Therefore, we aimed to determine whether VDR gene single nucleotide polymorphisms (SNPs) rs154410, rs7975232, rs731236, rs2228570 and rs739837 contribute to GDM risk in Wuhan, China. Moreover, we aimed to explore their combined effects on the risk of GDM.MethodsPregnant women who had prenatal examinations at 24 to 28 weeks’ gestation in our hospital from January 15, 2018 to March 31, 2019 were included in this case-control study. After exclusion, a total of 1684 pregnant women (826 GDM patients and 858 non-diabetic controls) were recruited. The clinical information and blood samples were collected by trained interviewers and nurses. Genotyping of candidate SNPs was conducted on the Sequenom MassARRAY platform. Statistical analyses including t-test, ANOVA, chi-square test and logistic regression were performed to the data with SPSS Software to evaluate differences in genotype distribution and associations with GDM risk. Multifactor dimensionality reduction method was used to explore the gene-gene interactions on the risk of GDM.ResultsDifferences in age, pre-pregnancy BMI, family history of diabetes and previous history of GDM between the case and control groups were statistically significant (P < 0.05), whereas no significant differences were found in height, gravidity, parity, and age of menarche (P > 0.05). There were no significant differences at genotype distributions of the examined VDR gene SNPs (P > 0.05). After adjusting by age, pre-pregnancy BMI, family history of diabetes, the results of logistic regression analysis showed no associations of the five SNPs with GDM in all the four genotype models(P > 0.05). Furthermore, there were no gene-gene interactions on the GDM risk among the five examined VDR gene SNPs.ConclusionsThe VDR gene SNPs rs154410, rs7975232, rs731236, rs2228570 and rs739837 showed neither significant associations nor gene-gene interactions with GDM in Wuhan, China.

Interactions among IGF-1, AKT2, FOXO1, and FOXO3 variations and between genes and physical activities on physical performance in community-dwelling elders.

This study assessed the interactions among IGF-1, AKT2, FOXO1, and FOXO3 variations and the interactions of gene and physical activity on handgrip strength, arm muscle mass-adjusted handgrip (armGrip), gait speed (GS), timed up and go (TUG), and leg press strength (LPS). Nine single nucleotide polymorphisms (SNPs) containing three IGF-1 SNPs (rs6214, rs5742692, and rs35767), two AKT2 SNPs (rs892119 and rs35817154), two FOXO1 SNPs (rs17446593 and rs10507486), and two FOXO3 SNPs (rs9480865 and rs2153960) were genotyped in 472 unrelated elders with a mean age of 73.8 years. We observed significant interactions of IGF-1 SNP rs6214 and rs35767 with regular physical activity on TUG and GS; and AKT2 SNP rs892119 and FOXO3 SNP rs9480865 with regular physical activity on armGrip. Genotype GG of IGF-1 rs6214 and rs35767 in individuals without regular physical activity had poor performance in TUG and GS, as well as GG of AKT2 rs892119 decreased armGrip in individuals without regular physical activity. After FDR adjustment, no significant gene–gene interactions were found. A sedentary lifestyle may increase the risk of impairing physical performance and regular physical activity is a remedy for sarcopenia, even a little regular physical activity can overcome carrying some risk alleles in this pathway.

Large-Scale Analyses Provide No Evidence for Gene-Gene Interactions Influencing Type 2 Diabetes Risk.

A growing number of genetic loci have been shown to influence individual predisposition to type 2 diabetes (T2D). Despite longstanding interest in understanding whether nonlinear interactions between these risk variants additionally influence T2D risk, the ability to detect significant gene-gene interaction (GGI) effects has been limited to date. To increase power to detect GGI effects, we combined recent advances in the fine-mapping of causal T2D risk variants with the increased sample size available within UK Biobank (375,736 unrelated European participants, including 16,430 with T2D). In addition to conventional single variant–based analysis, we used a complementary polygenic score–based approach, which included partitioned T2D risk scores that capture biological processes relevant to T2D pathophysiology. Nevertheless, we found no evidence in support of GGI effects influencing T2D risk. The current study was powered to detect interactions between common variants with odds ratios >1.2, so these findings place limits on the contribution of GGIs to the overall heritability of T2D.

Interaction of two functional genetic variants LOXL1 rs1048661 and VEGFA rs3025039 on the risk of age-related macular degeneration in Chinese women.

BackgroundCumulative evidence indicates that LOXL1 and VEGF-a play important roles in extracellular matrix formation and angiogenesis, respectively. The disorder of extracellular matrix and angiogenesis are the key factors of pathogenesis of age-related macular degeneration (AMD). We hypothesized that rs1048661 (T>G) in the LOXL1 gene and rs3025039 (C>T) in the VEGFA gene might be associated with risk of AMD.MethodsA total of 533 unrelated Chinese subjects, 286 cases (247 with early AMD and 39 with late neovascular AMD) and 247 controls, were included in the study. The gene sequences of LOXL1 rs1048661 and VEGFA rs3025039 were amplified by polymerase chain reaction and genotyped. Interaction between rs1048661 and rs3025039 on AMD risk was also assessed.ResultsLOXL1 rs1048661 but not VEGFA rs3025039 was associated with a significantly increased risk of AMD. The adjusted odds ratio was 1.6 (95% CI, 1.1–2.5) for rs1048661 TT + GT genotype compared with GG homozygotes in the dominant model analysis. Moreover, there was a significant gene-gene interaction between these two polymorphic loci. In VEGFA rs3025039 CC + CT genotype which indicated sufficient expression of VEGF-a, LOXL1 rs1048661 had odds ratios of 1.7 (95% CI, 1.1–2.7) for early AMD and 3.6 (95% CI, 1.1–12.3) for late neovascular AMD in the dominant model analysis. However, LOXL1 rs1048661 did not confer the risk of AMD in subjects harboring VEGFA rs3025039 TT genotype which indicated decreased expression of VEGF-a.ConclusionsOur findings suggest that LOXL1 rs1048661 (T>G) may be involved in the risk of AMD. In addition, LOXL1 rs1048661 and VEGFA rs3025039 interacted to confer the development of AMD, especially for late-stage neovascular AMD. Our data need to be further validated.

Association of GSTP1, GSTT1 and GSTM1 Gene Variants with Coronary Artery Disease in Iranian Population: A Case-Control Study.

BackgroundCoronary artery disease (CAD) is a multifactorial disease that may be caused by the interaction between environmental and genetic risk factors. Glutathione S-transferases (GSTs) are known to participate in detoxification and metabolism of a wide range of xenobiotic compounds and oxidative stress products. Considering the interaction between environmental and genetic factors in CAD, we investigated the genetic polymorphisms of GSTM1, GSTT1, and GSTP1 in the Iranian population.Patients and MethodsTwo hundred and forty-four CAD cases and 281 healthy controls were studied. The genotype of GSTM1, GSTT1, and GSTP1 genes was determined by multiplex polymerase chain reaction (PCR) and PCR-restriction fragment length polymorphism (PCR-RFLP) techniques. Multivariable logistic regression analysis was used to calculate the odds ratios (ORs) and 95% confidence intervals (CI). Multifactor dimensionality reduction (MDR) analysis was also carried out to analyze the gene–gene and gene–environment interaction.ResultsThe genotype and allele distribution of the three variations were not significantly different between CAD patients and controls (p > 0.05). The subgroup analysis revealed no significant gene–gene interactions or gene–gene combination effects linked to CAD susceptibility. However, MDR analysis selected the GSTM, GSTT pairwise and three genes combination models associated with the susceptibility to CAD. In addition, its result revealed that smoking in combination with GSTM1 (two-way) and GSTT, GSTP (three-way) genes might increase the risk of CAD. Furthermore, a significant interaction between GSTT1-null polymorphism and dyslipidemia was found in multivariable logistic regression analyses in the gene–environmental interactions on CAD risk.ConclusionOur results suggest that the GSTM1, GSTT1 and GSTP1 genetic variations are not directly associated with the susceptibility to CAD in Iranian patients. Due to MDR results, there might be a non-linear association between interactions of two or three genes and smoking with CAD. There is also an association between CAD risk factors and GST variations, which requires supplementary confirmation with larger sample sizes.

Implication of genetic variants in overweight and obesity susceptibility among the young Arab population of the United Arab Emirates.

Overweight and obesity are major risk factors for Type 2 Diabetes Mellitus (T2DM), cardiovascular disease (CVD) and cancer. Genetic predisposition has been shown to play a key role in obesity, and genome-wide association studies (GWAS) have identified multiple loci linked with obesity in various ethnic groups. The aim of this study was to validate the reported genetic variants associated with obesity and overweight in a young UAE Arab population. Twenty-two associated single nucleotide polymorphisms (SNPs) at 11 loci (FTO, MC4R, TMEM18, KCTD15, MTCH2, SH2B1, TFAP2B, GNPDA2, NEGR1, PCSK1 and BDNF) were studied in 392 controls and 318 overweight/obese young Emiratis (aged 18-35years). After adjusting for age and smoking, rs3751812 of the FTO gene was associated with overweight/obesity in male participants (p-value<0.016), while SNPs rs17782313, rs571312 of the MC4R gene and rs12463617 of the TMEM18 gene were significantly associated with overweight/obesity in female participants (p-value=0.001, 0.028, 0.044, respectively). Follow-up association tests and logistic regression revealed the contribution of the FTO rs3751812 and MC4R rs571213 SNPs to the risk of overweight/obesity after adjusting for age, sex and smoking (p-value=0.044, 0.049, respectively). In addition, the FTO rs3751812 was associated with the risk of overweight/obesity after adjusting for the effect of other markers (rs17782313, rs571312, rs2867125, rs6548238 and rs12463617) (p-value=0.035). A significant gene-gene interaction was seen between FTO, MCR4 and TMEM18 (p-value=0.013). Our data demonstrates that rs3751812 of the FTO gene is the key SNP associated with risk of overweight/obesity among the young UAE Arab population, in alignment with previous findings. Our results also indicate that the identified genes stratify with sex and risk of overweight/obesity. In addition to their direct association with overweight/obesity, rs17782313 and rs571312, as well as rs2867125 and rs6548238, may have a modifying effect on the risk of overweight/obesity caused by the rs3751812. Population-specific, sex-specific genetic profiling is important in understanding the heritability of obesity.